Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Adding event selections to B2G offline DQM workflow #4820

Merged
merged 5 commits into from Aug 1, 2014
Merged

Adding event selections to B2G offline DQM workflow #4820

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
11c28cd
Adding first attempt at actual event selections for B2G DQM
rappoccio Jul 29, 2014
47fe0d3
Adding python config for event selections
rappoccio Jul 29, 2014
1065645
Adding simple electron selections
rappoccio Jul 29, 2014
2833bcf
Merged rappoccio_B2G_DQM_EventSel from repository rappoccio
rappoccio Jul 29, 2014
727f53b
Reversing a cut that was previously the wrong comparison
rappoccio Jul 29, 2014
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
26 changes: 25 additions & 1 deletion DQM/Physics/python/B2GDQM_cfi.py
View file
Open in desktop
@@ -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
Open in desktop
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