Skip to content

Commit

Permalink
Merge pull request #4820 from rappoccio/rappoccio_B2G_DQM_EventSel
Browse files Browse the repository at this point in the history 
Adding event selections to B2G offline DQM workflow
  • Loading branch information
@davidlange6
davidlange6 committed Aug 1, 2014
2 parents 403bad1 + 727f53b commit ba924c8
Show file tree
Hide file tree
Showing 3 changed files with 398 additions and 17 deletions.
26 changes: 25 additions & 1 deletion DQM/Physics/python/B2GDQM_cfi.py
View file
Original file line number Diff line number Diff line change
@@ -1,5 +1,7 @@
import FWCore.ParameterSet.Config as cms

from math import pi

B2GDQM = cms.EDAnalyzer(
"B2GDQM",

Expand All @@ -22,7 +24,29 @@
50.,
100.
),
pfMETCollection = cms.InputTag("pfMet")
pfMETCollection = cms.InputTag("pfMet"),

cmsTagLabel = cms.InputTag("cmsTopTagPFJetsCHS"),
muonSrc = cms.InputTag("muons"),
electronSrc = cms.InputTag("gedGsfElectrons"),

allHadPtCut = cms.double(400.0),
allHadRapidityCut = cms.double(1.0),
allHadDeltaPhiCut = cms.double( pi / 2.0),

muonSelect = cms.string("pt > 45.0 & abs(eta)<2.1 & isGlobalMuon & abs(globalTrack.d0)<1 & abs(globalTrack.dz)<20"),
semiMu_HadJetPtCut = cms.double(400.0),
semiMu_LepJetPtCut = cms.double(30.0),
semiMu_dphiHadCut = cms.double( pi / 2.0),
semiMu_dRMin = cms.double( 0.5 ),
semiMu_ptRel = cms.double( 25.0 ),

elecSelect = cms.string("pt > 45.0 & abs(eta)<2.5 & abs(gsfTrack.d0)<1 & abs(gsfTrack.dz)<20"),
semiE_HadJetPtCut = cms.double(400.0),
semiE_LepJetPtCut = cms.double(30.0),
semiE_dphiHadCut = cms.double( pi / 2.0),
semiE_dRMin = cms.double( 0.5 ),
semiE_ptRel = cms.double( 25.0 )


)
300 changes: 294 additions & 6 deletions DQM/Physics/src/B2GDQM.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -98,7 +98,41 @@ B2GDQM::B2GDQM(const edm::ParameterSet& ps){
jetlabelEnd = jetLabels_.end(); jetlabel != jetlabelEnd; ++jetlabel ) {
jetTokens_.push_back( consumes<edm::View<reco::Jet> >( *jetlabel ) );
}
jetPtMins_ = ps.getParameter<std::vector<double> > ("jetPtMins");
cmsTagLabel_ = ps.getParameter<edm::InputTag >("cmsTagLabel");
cmsTagToken_ = consumes<edm::View<reco::BasicJet> > ( cmsTagLabel_ );


muonToken_ = consumes<edm::View<reco::Muon> >( ps.getParameter<edm::InputTag>("muonSrc") );
electronToken_ = consumes<edm::View<reco::GsfElectron> >( ps.getParameter<edm::InputTag>("electronSrc") );

jetPtMins_ = ps.getParameter<std::vector<double> > ("jetPtMins");
allHadPtCut_ = ps.getParameter<double> ("allHadPtCut");
allHadRapidityCut_ = ps.getParameter<double> ("allHadRapidityCut");
allHadDeltaPhiCut_ = ps.getParameter<double> ("allHadDeltaPhiCut");

semiMu_HadJetPtCut_ = ps.getParameter<double> ("semiMu_HadJetPtCut");
semiMu_LepJetPtCut_ = ps.getParameter<double> ("semiMu_LepJetPtCut");
semiMu_dphiHadCut_ = ps.getParameter<double> ("semiMu_dphiHadCut");
semiMu_dRMin_ = ps.getParameter<double> ("semiMu_dRMin");
semiMu_ptRel_ = ps.getParameter<double> ("semiMu_ptRel");
muonSelect_ = std::shared_ptr<StringCutObjectSelector<reco::Muon> > (
new StringCutObjectSelector<reco::Muon>(
ps.getParameter<std::string>( "muonSelect")
)
);

semiE_HadJetPtCut_ = ps.getParameter<double> ("semiE_HadJetPtCut");
semiE_dphiHadCut_ = ps.getParameter<double> ("semiE_dphiHadCut");
semiE_dRMin_ = ps.getParameter<double> ("semiE_dRMin");
semiE_ptRel_ = ps.getParameter<double> ("semiE_ptRel");
elecSelect_ = std::shared_ptr<StringCutObjectSelector<reco::GsfElectron> > (
new StringCutObjectSelector<reco::GsfElectron>(
ps.getParameter<std::string>( "elecSelect")
)
);



PFJetCorService_ = ps.getParameter<std::string>("PFJetCorService");

// MET
Expand Down Expand Up @@ -164,7 +198,7 @@ void B2GDQM::bookHistos(DQMStore* bei){

bei->cd();

//--- Event Interpretation
//--- Jets

for ( unsigned int icoll = 0; icoll < jetLabels_.size(); ++icoll ) {
std::stringstream ss;
Expand Down Expand Up @@ -194,7 +228,53 @@ void B2GDQM::bookHistos(DQMStore* bei){
pfMet_phi = bei->book1D("pfMet_phi", "Pf Missing p_{T} #phi;#phi (radians)", 35, -3.5, 3.5 );



//--- Mu+Jets
bei->setCurrentFolder("Physics/B2G/SemiMu");
semiMu_muPt = bei->book1D("semiMu_muPt", "Pt of Muon in #mu+Jets Channel (GeV)", 50, 0.0, 1000) ;
semiMu_muEta = bei->book1D("semiMu_muEta", "#eta of Muon in #mu+Jets Channel", 60, -6.0, 6.0) ;
semiMu_muPhi = bei->book1D("semiMu_muPhi", "#phi of Muon in #mu+Jets Channel (radians)", 60, -3.14159, 3.14159) ;
semiMu_muDRMin = bei->book1D("semiMu_muDRMin", "#Delta R(E,nearest jet) in #mu+Jets Channel", 50, 0, 10.0);
semiMu_muPtRel = bei->book1D("semiMu_muPtRel", "p_{T}^{REL} in #mu+Jets Channel", 60, 0, 300.);
semiMu_hadJetDR = bei->book1D("semiMu_hadJetDR", "#Delta R(E,had jet) in #mu+Jets Channel", 50, 0, 10.0);
semiMu_hadJetPt = bei->book1D("semiMu_hadJetPt", "Pt of Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 1000) ;
semiMu_hadJetY = bei->book1D("semiMu_hadJetY", "Rapidity of Leading Hadronic Jet in #mu+Jets Channel", 60, -6.0, 6.0) ;
semiMu_hadJetPhi = bei->book1D("semiMu_hadJetPhi", "#phi of Leading Hadronic Jet in #mu+Jets Channel (radians)",60, -3.14159, 3.14159) ;
semiMu_hadJetMass = bei->book1D("semiMu_hadJetMass", "Mass of Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 500) ;
semiMu_hadJetMinMass = bei->book1D("semiMu_hadJetminMass","Minimum Mass Pairing for Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0 , 250.0);
semiMu_mttbar = bei->book1D("semiMu_mttbar", "Mass of #mu+Jets ttbar Candidate", 100, 0., 5000.);


//--- E+Jets
bei->setCurrentFolder("Physics/B2G/SemiE");
semiE_ePt = bei->book1D("semiE_ePt", "Pt of Electron in e+Jets Channel (GeV)", 50, 0.0, 1000) ;
semiE_eEta = bei->book1D("semiE_eEta", "#eta of Electron in e+Jets Channel", 60, -6.0, 6.0) ;
semiE_ePhi = bei->book1D("semiE_ePhi", "#phi of Electron in e+Jets Channel (radians)", 60, -3.14159, 3.14159) ;
semiE_eDRMin = bei->book1D("semiE_eDRMin", "#Delta R(E,nearest jet) in e+Jets Channel", 50, 0, 10.0);
semiE_ePtRel = bei->book1D("semiE_ePtRel", "p_{T}^{REL} in e+Jets Channel", 60, 0, 300.);
semiE_hadJetDR = bei->book1D("semiE_hadJetDR", "#Delta R(E,had jet) in e+Jets Channel", 50, 0, 10.0);
semiE_hadJetPt = bei->book1D("semiE_hadJetPt", "Pt of Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 1000) ;
semiE_hadJetY = bei->book1D("semiE_hadJetY", "Rapidity of Leading Hadronic Jet in e+Jets Channel", 60, -6.0, 6.0) ;
semiE_hadJetPhi = bei->book1D("semiE_hadJetPhi", "#phi of Leading Hadronic Jet in e+Jets Channel (radians)",60, -3.14159, 3.14159) ;
semiE_hadJetMass = bei->book1D("semiE_hadJetMass", "Mass of Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 500) ;
semiE_hadJetMinMass = bei->book1D("semiE_hadJetminMass","Minimum Mass Pairing for Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0 , 250.0);
semiE_mttbar = bei->book1D("semiE_mttbar", "Mass of e+Jets ttbar Candidate", 100, 0., 5000.);


//--- All-hadronic
bei->setCurrentFolder("Physics/B2G/AllHad");
allHad_pt0 = bei->book1D("allHad_pt0", "Pt of Leading All-Hadronic PFJet (GeV)", 50, 0.0, 1000) ;
allHad_y0 = bei->book1D("allHad_y0", "Rapidity of Leading All-Hadronic PFJet", 60, -6.0, 6.0) ;
allHad_phi0 = bei->book1D("allHad_phi0", "#phi of Leading All-Hadronic PFJet (radians)",60, -3.14159, 3.14159) ;
allHad_mass0 = bei->book1D("allHad_mass0", "Mass of Leading All-Hadronic PFJet (GeV)", 50, 0.0, 500) ;
allHad_minMass0 = bei->book1D("allHad_minMass0","Minimum Mass Pairing for Leading All-Hadronic PFJet (GeV)", 50, 0.0 , 250.0);
allHad_pt1 = bei->book1D("allHad_pt1", "Pt of Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 1000) ;
allHad_y1 = bei->book1D("allHad_y1", "Rapidity of Subleading All-Hadronic PFJet", 60, -6.0, 6.0) ;
allHad_phi1 = bei->book1D("allHad_phi1", "#phi of Subleading All-Hadronic PFJet (radians)",60, -3.14159, 3.14159) ;
allHad_mass1 = bei->book1D("allHad_mass1", "Mass of Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 500) ;
allHad_minMass1 = bei->book1D("allHad_minMass1","Minimum Mass Pairing for Subleading All-Hadronic PFJet (GeV)", 50, 0.0 , 250.0);
allHad_mttbar = bei->book1D("allHad_mttbar", "Mass of All-Hadronic ttbar Candidate", 100, 0., 5000.);


bei->cd();
}

Expand All @@ -204,12 +284,15 @@ void B2GDQM::bookHistos(DQMStore* bei){
//
void B2GDQM::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup){

analyzeEventInterpretation(iEvent, iSetup);
analyzeJets (iEvent, iSetup);
analyzeSemiMu (iEvent, iSetup );
analyzeSemiE (iEvent, iSetup );
analyzeAllHad (iEvent, iSetup );
}


void B2GDQM::analyzeEventInterpretation(const Event & iEvent, const edm::EventSetup& iSetup){

void B2GDQM::analyzeJets(const Event & iEvent, const edm::EventSetup& iSetup){

// Loop over the different types of jets,
// Loop over the jets in that collection,
Expand Down Expand Up @@ -267,7 +350,7 @@ void B2GDQM::analyzeEventInterpretation(const Event & iEvent, const edm::EventSe

// For top-tagging, check the minimum mass pairing
if ( jetLabels_[icoll].label() == "cmsTopTagPFJetsCHS") {
CATopJetHelper helper(171.2, 80.4);
CATopJetHelper helper(173., 80.4);
reco::CATopJetProperties properties = helper(*basicjet);
if ( jet->numberOfDaughters() > 2 ) {
boostedJet_minMass[icoll]->Fill ( properties.minMass );
Expand Down Expand Up @@ -308,6 +391,211 @@ void B2GDQM::analyzeEventInterpretation(const Event & iEvent, const edm::EventSe

}

void B2GDQM::analyzeAllHad(const Event & iEvent, const edm::EventSetup& iSetup){

edm::Handle<edm::View<reco::BasicJet> > jetCollection;
bool validJets = iEvent.getByToken(cmsTagToken_, jetCollection );
if(!validJets) return;

// Require two back-to-back jets at high pt with |delta y| < 1.0
if ( jetCollection->size() < 2 ) return;
edm::Ptr<reco::BasicJet> jet0 = jetCollection->ptrAt(0);
edm::Ptr<reco::BasicJet> jet1 = jetCollection->ptrAt(1);
if (jet0.isAvailable() == false || jet1.isAvailable() == false) return;
if ( jet0->pt() < allHadPtCut_ || jet1->pt() < allHadPtCut_ ) return;
if ( std::abs( jet0->rapidity() - jet1->rapidity() ) > allHadRapidityCut_ ) return;
if ( std::abs( reco::deltaPhi( jet0->phi(), jet1->phi() ) ) < M_PI*0.5 ) return;


CATopJetHelper helper(173., 80.4);

allHad_pt0->Fill( jet0->pt() );
allHad_y0->Fill( jet0->rapidity() );
allHad_phi0->Fill( jet0->phi() );
allHad_mass0->Fill( jet0->mass() );
reco::CATopJetProperties properties0 = helper(*jet0);
if ( jet0->numberOfDaughters() > 2 ) {
allHad_minMass0->Fill ( properties0.minMass );
} else {
allHad_minMass0->Fill ( -1.0 );
}

allHad_pt1->Fill( jet1->pt() );
allHad_y1->Fill( jet1->rapidity() );
allHad_phi1->Fill( jet1->phi() );
allHad_mass1->Fill( jet1->mass() );
reco::CATopJetProperties properties1 = helper(*jet1);
if ( jet1->numberOfDaughters() > 2 ) {
allHad_minMass1->Fill ( properties1.minMass );
} else {
allHad_minMass1->Fill ( -1.0 );
}

auto p4cand = (jet0->p4() + jet1->p4());
allHad_mttbar->Fill( p4cand.mass() );
}


void B2GDQM::analyzeSemiMu(const Event & iEvent, const edm::EventSetup& iSetup){

edm::Handle<edm::View<reco::Muon> > muonCollection;
bool validMuons = iEvent.getByToken( muonToken_, muonCollection );

if( !validMuons) return;
if( muonCollection->size() < 1) return;
reco::Muon const & muon = muonCollection->at(0);
if ( ! (*muonSelect_)(muon) ) return;

edm::Handle<edm::View<reco::BasicJet> > jetCollection;
bool validJets = iEvent.getByToken(cmsTagToken_, jetCollection );
if(!validJets) return;
if ( jetCollection->size() < 2 ) return;


double pt0 = -1.0;
double dRMin = 999.0;
edm::Ptr<reco::BasicJet> hadJet; // highest pt jet with dphi(lep,jet) > pi/2
edm::Ptr<reco::BasicJet> lepJet; // closest jet to lepton with pt > ptMin

for ( auto ijet = jetCollection->begin(), ijetBegin = ijet, ijetEnd = jetCollection->end();
ijet != ijetEnd; ++ijet ) {
// Hadronic jets
if ( std::abs( reco::deltaPhi(muon, *ijet) ) > M_PI*0.5 ) {
if ( ijet->pt() > pt0 && ijet->p() > semiMu_HadJetPtCut_ ) {
hadJet = jetCollection->ptrAt( ijet - ijetBegin);
pt0 = hadJet->pt();
}
}
// Leptonic jets
else if ( ijet->pt() > semiMu_LepJetPtCut_ ) {
auto idRMin = reco::deltaR( muon, *ijet );
if ( idRMin < dRMin ) {
lepJet = jetCollection->ptrAt( ijet - ijetBegin);
dRMin = idRMin;
}
}
}
if ( hadJet.isAvailable() == false || lepJet.isAvailable() == false ) return;

auto lepJetP4 = lepJet->p4();
auto muonP4 = muon.p4();

double tot = lepJetP4.mag2();
double ss = muonP4.Dot(lepJet->p4() );
double per = muonP4.mag2();
if (tot > 0.0) per -= ss*ss/tot;
if (per < 0) per = 0;
double ptRel = per;
bool pass2D = dRMin > semiMu_dRMin_ || ptRel > semiMu_ptRel_ ;

if ( !pass2D) return;

CATopJetHelper helper(173., 80.4);

semiMu_muPt->Fill( muon.pt() );
semiMu_muEta->Fill( muon.eta() );
semiMu_muPhi->Fill( muon.phi() );
semiMu_muDRMin->Fill( dRMin );
semiMu_muPtRel->Fill( ptRel );

semiMu_hadJetDR->Fill( reco::deltaR(muon, *hadJet) );
semiMu_mttbar->Fill( 0.0 );

semiMu_hadJetPt->Fill( hadJet->pt() );
semiMu_hadJetY->Fill( hadJet->rapidity() );
semiMu_hadJetPhi->Fill( hadJet->phi() );
semiMu_hadJetMass->Fill( hadJet->mass() );
reco::CATopJetProperties properties0 = helper(*hadJet);
if ( hadJet->numberOfDaughters() > 2 ) {
semiMu_hadJetMinMass->Fill ( properties0.minMass );
} else {
semiMu_hadJetMinMass->Fill ( -1.0 );
}

}



void B2GDQM::analyzeSemiE(const Event & iEvent, const edm::EventSetup& iSetup){

edm::Handle<edm::View<reco::GsfElectron> > electronCollection;
bool validElectrons = iEvent.getByToken( electronToken_, electronCollection );

if( !validElectrons) return;
if( electronCollection->size() < 1) return;
reco::GsfElectron const & electron = electronCollection->at(0);
if ( ! (*elecSelect_)(electron) ) return;

edm::Handle<edm::View<reco::BasicJet> > jetCollection;
bool validJets = iEvent.getByToken(cmsTagToken_, jetCollection );
if(!validJets) return;
if ( jetCollection->size() < 2 ) return;

double pt0 = -1.0;
double dRMin = 999.0;
edm::Ptr<reco::BasicJet> hadJet; // highest pt jet with dphi(lep,jet) > pi/2
edm::Ptr<reco::BasicJet> lepJet; // closest jet to lepton with pt > ptMin

for ( auto ijet = jetCollection->begin(), ijetBegin = ijet, ijetEnd = jetCollection->end();
ijet != ijetEnd; ++ijet ) {
// Hadronic jets
if ( std::abs( reco::deltaPhi(electron, *ijet) ) > M_PI*0.5 ) {
if ( ijet->pt() > pt0 && ijet->p() > semiE_HadJetPtCut_ ) {
hadJet = jetCollection->ptrAt( ijet - ijetBegin);
pt0 = hadJet->pt();
}
}
// Leptonic jets
else if ( ijet->pt() > semiE_LepJetPtCut_ ) {
auto idRMin = reco::deltaR( electron, *ijet );
if ( idRMin < dRMin ) {
lepJet = jetCollection->ptrAt( ijet - ijetBegin);
dRMin = idRMin;
}
}
}
if ( hadJet.isAvailable() == false || lepJet.isAvailable() == false ) return;


auto lepJetP4 = lepJet->p4();
auto electronP4 = electron.p4();

double tot = lepJetP4.mag2();
double ss = electronP4.Dot(lepJet->p4() );
double per = electronP4.mag2();
if (tot > 0.0) per -= ss*ss/tot;
if (per < 0) per = 0;
double ptRel = per;
bool pass2D = dRMin > semiE_dRMin_ || ptRel > semiE_ptRel_ ;

if ( !pass2D) return;

CATopJetHelper helper(173., 80.4);

semiE_ePt->Fill( electron.pt() );
semiE_eEta->Fill( electron.eta() );
semiE_ePhi->Fill( electron.phi() );
semiE_eDRMin->Fill( dRMin );
semiE_ePtRel->Fill( ptRel );

semiE_hadJetDR->Fill( reco::deltaR(electron, *hadJet) );
semiE_mttbar->Fill( 0.0 );

semiE_hadJetPt->Fill( hadJet->pt() );
semiE_hadJetY->Fill( hadJet->rapidity() );
semiE_hadJetPhi->Fill( hadJet->phi() );
semiE_hadJetMass->Fill( hadJet->mass() );
reco::CATopJetProperties properties0 = helper(*hadJet);
if ( hadJet->numberOfDaughters() > 2 ) {
semiE_hadJetMinMass->Fill ( properties0.minMass );
} else {
semiE_hadJetMinMass->Fill ( -1.0 );
}

}



// -- End Luminosity Block
//
void B2GDQM::endLuminosityBlock(edm::LuminosityBlock const& lumiSeg, edm::EventSetup const& eSetup) {
Expand Down
Loading

0 comments on commit ba924c8

Please sign in to comment.