|
|
|
#include <EwkElecDQM.h>
DQM offline for EWK Electrons
Definition at line 19 of file EwkElecDQM.h.
| EwkElecDQM::EwkElecDQM | ( | const edm::ParameterSet & | cfg | ) |
Definition at line 39 of file EwkElecDQM.cc.
References init_histograms(), isValidHltConfig_, cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), and theDbe.
:
// Input collections
trigTag_(cfg.getUntrackedParameter<edm::InputTag> ("TrigTag", edm::InputTag("TriggerResults::HLT"))),
// muonTag_(cfg.getUntrackedParameter<edm::InputTag> ("MuonTag", edm::InputTag("muons"))),
elecTag_(cfg.getUntrackedParameter<edm::InputTag> ("ElecTag", edm::InputTag("gsfElectrons"))),
metTag_(cfg.getUntrackedParameter<edm::InputTag> ("METTag", edm::InputTag("met"))),
// metIncludesMuons_(cfg.getUntrackedParameter<bool> ("METIncludesMuons", false)),
jetTag_(cfg.getUntrackedParameter<edm::InputTag> ("JetTag", edm::InputTag("sisCone5CaloJets"))),
vertexTag_ (cfg.getUntrackedParameter<edm::InputTag> ("VertexTag", edm::InputTag("offlinePrimaryVertices"))),
// Main cuts
// muonTrig_(cfg.getUntrackedParameter<std::string> ("MuonTrig", "HLT_Mu9")),
//elecTrig_(cfg.getUntrackedParameter<std::vector< std::string > >("ElecTrig", "HLT_Ele10_SW_L1R")),
elecTrig_(cfg.getUntrackedParameter<std::vector< std::string > >("ElecTrig")),
// ptCut_(cfg.getUntrackedParameter<double>("PtCut", 25.)),
ptCut_(cfg.getUntrackedParameter<double>("PtCut", 10.)),
// etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.1)),
etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.4)),
sieieCutBarrel_(cfg.getUntrackedParameter<double>("SieieBarrel", 0.01)),
sieieCutEndcap_(cfg.getUntrackedParameter<double>("SieieEndcap", 0.028)),
detainCutBarrel_(cfg.getUntrackedParameter<double>("DetainBarrel", 0.0071)),
detainCutEndcap_(cfg.getUntrackedParameter<double>("DetainEndcap", 0.0066)),
// isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso", true)),
// isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso", false)),
// isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03", 0.1)),
ecalIsoCutBarrel_(cfg.getUntrackedParameter<double>("EcalIsoCutBarrel", 5.7)),
ecalIsoCutEndcap_(cfg.getUntrackedParameter<double>("EcalIsoCutEndcap", 5.0)),
hcalIsoCutBarrel_(cfg.getUntrackedParameter<double>("HcalIsoCutBarrel", 8.1)),
hcalIsoCutEndcap_(cfg.getUntrackedParameter<double>("HcalIsoCutEndcap", 3.4)),
trkIsoCutBarrel_(cfg.getUntrackedParameter<double>("TrkIsoCutBarrel", 7.2)),
trkIsoCutEndcap_(cfg.getUntrackedParameter<double>("TrkIsoCutEndcap", 5.1)),
mtMin_(cfg.getUntrackedParameter<double>("MtMin", -999999)),
mtMax_(cfg.getUntrackedParameter<double>("MtMax", 999999.)),
metMin_(cfg.getUntrackedParameter<double>("MetMin", -999999.)),
metMax_(cfg.getUntrackedParameter<double>("MetMax", 999999.)),
// acopCut_(cfg.getUntrackedParameter<double>("AcopCut", 2.)),
// Muon quality cuts
// dxyCut_(cfg.getUntrackedParameter<double>("DxyCut", 0.2)),
// normalizedChi2Cut_(cfg.getUntrackedParameter<double>("NormalizedChi2Cut", 10.)),
// trackerHitsCut_(cfg.getUntrackedParameter<int>("TrackerHitsCut", 11)),
// isAlsoTrackerMuon_(cfg.getUntrackedParameter<bool>("IsAlsoTrackerMuon", true)),
// Z rejection
// ptThrForZ1_(cfg.getUntrackedParameter<double>("PtThrForZ1", 20.)),
// ptThrForZ2_(cfg.getUntrackedParameter<double>("PtThrForZ2", 10.)),
// Top rejection
eJetMin_(cfg.getUntrackedParameter<double>("EJetMin", 999999.)),
nJetMax_(cfg.getUntrackedParameter<int>("NJetMax", 999999)),
PUMax_(cfg.getUntrackedParameter<unsigned int>("PUMax", 60)),
PUBinCount_(cfg.getUntrackedParameter<unsigned int>("PUBinCount", 12))
// caloJetCollection_(cfg.getUntrackedParameter<edm:InputTag>("CaloJetCollection","sisCone5CaloJets"))
{
isValidHltConfig_ = false;
// access to dbe
theDbe = Service<DQMStore>().operator->();
theDbe->setCurrentFolder("Physics/EwkElecDQM");
init_histograms();
}
| void EwkElecDQM::analyze | ( | const edm::Event & | ev, |
| const edm::EventSetup & | |||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 361 of file EwkElecDQM.cc.
References calcDeltaPhi(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), detainbarrel_after_, detainbarrel_before_, detainCutBarrel_, detainCutEndcap_, detainendcap_after_, detainendcap_before_, reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::dr03HcalTowerSumEt(), reco::GsfElectron::dr04TkSumPt(), ecalisobarrel_after_, ecalisobarrel_before_, ecalIsoCutBarrel_, ecalIsoCutEndcap_, ecalisoendcap_after_, ecalisoendcap_before_, eJetMin_, elecTag_, elecTrig_, metsig::electron, reco::LeafCandidate::energy(), reco::LeafCandidate::et(), eta(), reco::LeafCandidate::eta(), eta_after_, eta_before_, etaCut_, MonitorElement::Fill(), newFWLiteAna::found, edm::Event::getByLabel(), hcalisobarrel_after_, hcalisobarrel_before_, hcalIsoCutBarrel_, hcalIsoCutEndcap_, hcalisoendcap_after_, hcalisoendcap_before_, hltConfigProvider_, i, getHLTprescales::index, invmass_after_, invmass_before_, invmassPU_afterZ_, invmassPU_before_, reco::Vertex::isValid(), j, metsig::jet, jet_et_after_, jet_et_before_, jet_eta_after_, jet_eta_before_, jetTag_, gen::k, LogTrace, reco::LeafCandidate::massSqr(), CaloMET_cfi::met, met_after_, met_before_, metMax_, metMin_, metTag_, mt_after_, mt_before_, mtMax_, mtMin_, nall, neid, nelectrons_after_, nelectrons_before_, nGoodElectrons, niso, nJetMax_, njets_after_, njets_before_, npvs_afterZ_, npvs_before_, nrec, nsel, reco::LeafCandidate::phi(), HLTConfigProvider::prescaleSize(), HLTConfigProvider::prescaleValue(), reco::LeafCandidate::pt(), pt_after_, pt_before_, ptCut_, reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), sieiebarrel_after_, sieiebarrel_before_, sieieCutBarrel_, sieieCutEndcap_, sieieendcap_after_, sieieendcap_before_, reco::GsfElectron::sigmaIetaIeta(), mathSSE::sqrt(), trig_after_, trig_before_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), patRefSel_triggerSelection_cff::triggerResults, trigNames, trigTag_, trkisobarrel_after_, trkisobarrel_before_, trkIsoCutBarrel_, trkIsoCutEndcap_, trkisoendcap_after_, trkisoendcap_before_, GoodVertex_cfg::vertexCollection, and vertexTag_.
{
// Reset global event selection flags
bool rec_sel = false;
bool eid_sel = false;
bool iso_sel = false;
// bool hlt_sel = false;
bool all_sel = false;
// // Muon collection
// Handle<View<Muon> > muonCollection;
// if (!ev.getByLabel(muonTag_, muonCollection)) {
// LogWarning("") << ">>> Muon collection does not exist !!!";
// return;
// }
// unsigned int muonCollectionSize = muonCollection->size();
// Electron collection
Handle<View<GsfElectron> > electronCollection;
if (!ev.getByLabel(elecTag_, electronCollection)) {
//LogWarning("") << ">>> Electron collection does not exist !!!";
return;
}
unsigned int electronCollectionSize = electronCollection->size();
// Beam spot
Handle<reco::BeamSpot> beamSpotHandle;
if (!ev.getByLabel(InputTag("offlineBeamSpot"), beamSpotHandle)) {
//LogWarning("") << ">>> No beam spot found !!!";
return;
}
// MET
double met_px = 0.;
double met_py = 0.;
Handle<View<MET> > metCollection;
if (!ev.getByLabel(metTag_, metCollection)) {
//LogWarning("") << ">>> MET collection does not exist !!!";
return;
}
const MET& met = metCollection->at(0);
met_px = met.px();
met_py = met.py();
// if (!metIncludesMuons_) {
// for (unsigned int i=0; i<muonCollectionSize; i++) {
// const Muon& mu = muonCollection->at(i);
// if (!mu.isGlobalMuon()) continue;
// met_px -= mu.px();
// met_py -= mu.py();
// }
// }
double met_et = sqrt(met_px*met_px+met_py*met_py);
LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met_px << ", " << met_py << " [GeV]";
met_before_->Fill(met_et);
// Vertices in the event
int npvCount = 0;
Handle<View<reco::Vertex> > vertexCollection;
if (!ev.getByLabel(vertexTag_, vertexCollection)) {
LogError("") << ">>> Vertex collection does not exist !!!";
return;
}
for (unsigned int i=0; i<vertexCollection->size(); i++) {
const Vertex& vertex = vertexCollection->at(i);
if (vertex.isValid()) npvCount++;
}
npvs_before_->Fill(npvCount);
// Trigger
Handle<TriggerResults> triggerResults;
if (!ev.getByLabel(trigTag_, triggerResults)) {
//LogWarning("") << ">>> TRIGGER collection does not exist !!!";
return;
}
const edm::TriggerNames & trigNames = ev.triggerNames(*triggerResults);
bool trigger_fired = false;
/*
for (unsigned int i=0; i<triggerResults->size(); i++) {
if (triggerResults->accept(i)) {
LogTrace("") << "Accept by: " << i << ", Trigger: " << trigNames.triggerName(i);
}
}
*/
// the following gives error on CRAFT08 data where itrig1=19 (vector index out of range)
/*
int itrig1 = trigNames.triggerIndex(muonTrig_);
if (triggerResults->accept(itrig1)) trigger_fired = true;
*/
//suggested replacement: lm250909
for (unsigned int i=0; i<triggerResults->size(); i++)
{
std::string trigName = trigNames.triggerName(i);
bool found=false;
// for (unsigned int j = 0; j < elecTrig_.size(); j++)
// {
// if ( trigName == elecTrig_.at(j) && triggerResults->accept(i))
// {
// trigger_fired = true;
// }
// }
for(unsigned int index=0; index<elecTrig_.size() && found==false; index++) {
size_t trigPath = trigName.find(elecTrig_.at(index)); // 0 if found, pos if not
if (trigPath==0) found=true;
}
if(!found) continue;
bool prescaled=false;
for (unsigned int ps= 0; ps< hltConfigProvider_.prescaleSize(); ps++){
const unsigned int prescaleValue = hltConfigProvider_.prescaleValue(ps, trigName) ;
if (prescaleValue != 1) prescaled =true;
}
if(triggerResults->accept(i) && !prescaled) trigger_fired=true;
}
LogTrace("") << ">>> Trigger bit: " << trigger_fired << " for one of ( " ;
for (unsigned int k = 0; k < elecTrig_.size(); k++)
{
LogTrace("") << elecTrig_.at(k) << " ";
}
LogTrace("") << ")";
trig_before_->Fill(trigger_fired);
// // Loop to reject/control Z->mumu is done separately
// unsigned int nmuonsForZ1 = 0;
// unsigned int nmuonsForZ2 = 0;
// for (unsigned int i=0; i<muonCollectionSize; i++) {
// const Muon& mu = muonCollection->at(i);
// if (!mu.isGlobalMuon()) continue;
// double pt = mu.pt();
// if (pt>ptThrForZ1_) nmuonsForZ1++;
// if (pt>ptThrForZ2_) nmuonsForZ2++;
// }
// LogTrace("") << "> Z rejection: muons above " << ptThrForZ1_ << " [GeV]: " << nmuonsForZ1;
// LogTrace("") << "> Z rejection: muons above " << ptThrForZ2_ << " [GeV]: " << nmuonsForZ2;
// nz1_before_->Fill(nmuonsForZ1);
// nz2_before_->Fill(nmuonsForZ2);
// Jet collection
Handle<View<Jet> > jetCollection;
if (!ev.getByLabel(jetTag_, jetCollection)) {
//LogError("") << ">>> JET collection does not exist !!!";
return;
}
float electron_et = -8.0;
float electron_eta = -8.0;
float electron_phi = -8.0;
float electron2_et = -9.0;
float electron2_eta = -9.0;
float electron2_phi = -9.0;
//need to get some electron info so jets can be cleaned of them
for (unsigned int i=0; i<electronCollectionSize; i++)
{
const GsfElectron& elec = electronCollection->at(i);
if (i < 1)
{
electron_et = elec.pt();
electron_eta = elec.eta();
electron_phi = elec.phi();
}
if (i == 2)
{
electron2_et = elec.pt();
electron2_eta = elec.eta();
electron2_phi = elec.phi();
}
}
float jet_et = -8.0;
float jet_eta = -8.0;
int jet_count = 0;
float jet2_et = -9.0;
unsigned int jetCollectionSize = jetCollection->size();
int njets = 0;
for (unsigned int i=0; i<jetCollectionSize; i++) {
const Jet& jet = jetCollection->at(i);
float jet_current_et = jet.et();
// cout << "jet_current_et " << jet_current_et << endl;
// if it overlaps with electron, it is not a jet
if ( electron_et>0.0 && fabs(jet.eta()-electron_eta ) < 0.2
&& calcDeltaPhi(jet.phi(), electron_phi ) < 0.2)
continue;
if ( electron2_et>0.0&& fabs(jet.eta()-electron2_eta) < 0.2
&& calcDeltaPhi(jet.phi(), electron2_phi) < 0.2)
continue;
// if it has too low Et, throw away
// if (jet_current_et < eJetMin_) continue; //Keep if only want to plot above jet cut
if (jet.et()>eJetMin_)
{
njets++;
jet_count++;
}
if (jet_current_et > jet_et)
{
jet2_et = jet_et; // 2nd highest jet get's et from current highest
jet_et = jet.et(); // current highest jet gets et from the new highest
jet_eta = jet.eta();
}
else if (jet_current_et > jet2_et)
{
jet2_et = jet.et();
}
}
//Fill After all electron cuts (or both before and after)
if (jet_et>10) //don't want low energy "jets"
{
jet_et_before_ ->Fill(jet_et);
// jet2_et_before_ ->Fill(jet2_et);
jet_eta_before_ ->Fill(jet_eta);
}
LogTrace("") << ">>> Total number of jets: " << jetCollectionSize;
LogTrace("") << ">>> Number of jets above " << eJetMin_ << " [GeV]: " << njets;
njets_before_->Fill(njets);
// Start counting
nall++;
// Histograms per event should be done only once, so keep track of them
bool hlt_hist_done = false;
//bool minv_hist_done = false;
bool met_hist_done = false;
// bool nz1_hist_done = false;
// bool nz2_hist_done = false;
bool njets_hist_done = false;
// Central selection criteria
//const int NFLAGS = 13; // number of individual selection criteria
const int NFLAGS = 11; // number of individual selection criteria
// 0: pt cut | rec
// 1: eta cut | rec
// 2: sieie | eid
// 3: detain | eid
// 4: ecal iso | iso
// 5: hcal iso | iso
// 6: trk iso | iso
// 7: trigger fired | hlt/all
bool electron_sel[NFLAGS];
// for invariant mass calculation
// keep track of highest-pt electrons for initial (RECO) electrons
// and "good" electrons (passing all cuts)
// first array dimension is for first or second good electron
// second array dimension is for relevant quantities of good electron
// [0]: 1 for electron found or 0 for not found (if 0, no other quantities filled)
// [1]: mSqr
// [2]: E
// [3]: px
// [4]: py
// [5]: pz
// inv mass = sqrt(m_1^2 + m_2^2 + 2*(E_1*E_2 - (px1*px2 + py1*py2 + pz1+pz2) ) )
double electron[2][6];
double goodElectron[2][6];
nGoodElectrons = 0;
for (unsigned int i = 0; i < 2; i++)
{
for (unsigned int j = 0; j < 6; j++)
{
electron[i][j] = 0.;
goodElectron[i][j] = 0.;
}
}
for (unsigned int i=0; i<electronCollectionSize; i++)
{
for (int j=0; j<NFLAGS; ++j)
{
electron_sel[j] = false;
}
const GsfElectron& elec = electronCollection->at(i);
//if (!mu.isGlobalMuon()) continue;
//if (mu.globalTrack().isNull()) continue;
//if (mu.innerTrack().isNull()) continue;
LogTrace("") << "> elec: processing electron number " << i << "...";
//reco::TrackRef gm = mu.globalTrack();
//reco::TrackRef tk = mu.innerTrack();
// should have stuff for electron track?
if (i < 2)
{
electron[i][0] = 1.;
electron[i][1] = elec.massSqr();
electron[i][2] = elec.energy();
electron[i][3] = elec.px();
electron[i][4] = elec.py();
electron[i][5] = elec.pz();
}
// Pt,eta cuts
double pt = elec.pt();
double px = elec.px();
double py = elec.py();
double eta = elec.eta();
LogTrace("") << "\t... pt, eta: " << pt << " [GeV], " << eta;;
if (pt>ptCut_) electron_sel[0] = true;
if (fabs(eta)<etaCut_) electron_sel[1] = true;
bool isBarrel = false;
bool isEndcap = false;
if (eta < 1.4442 && eta > -1.4442)
{
isBarrel = true;
}
else if ((eta > 1.56 && eta < 2.4) || (eta < -1.56 && eta > -2.4))
{
isEndcap = true;
}
// // d0, chi2, nhits quality cuts
// double dxy = tk->dxy(beamSpotHandle->position());
// double normalizedChi2 = gm->normalizedChi2();
// double trackerHits = tk->numberOfValidHits();
// LogTrace("") << "\t... dxy, normalizedChi2, trackerHits, isTrackerMuon?: " << dxy << " [cm], " << normalizedChi2 << ", " << trackerHits << ", " << mu.isTrackerMuon();
// if (fabs(dxy)<dxyCut_) muon_sel[2] = true;
// if (normalizedChi2<normalizedChi2Cut_) muon_sel[3] = true;
// if (trackerHits>=trackerHitsCut_) muon_sel[4] = true;
// if (mu.isTrackerMuon()) muon_sel[5] = true;
pt_before_->Fill(pt);
eta_before_->Fill(eta);
// dxy_before_->Fill(dxy);
// chi2_before_->Fill(normalizedChi2);
// nhits_before_->Fill(trackerHits);
// tkmu_before_->Fill(mu.isTrackerMuon());
// Electron ID cuts
double sieie = (double) elec.sigmaIetaIeta();
double detain = (double) elec.deltaEtaSuperClusterTrackAtVtx(); // think this is detain
if (sieie < sieieCutBarrel_ && isBarrel) electron_sel[2] = true;
if (sieie < sieieCutEndcap_ && isEndcap) electron_sel[2] = true;
if (detain < detainCutBarrel_ && isBarrel) electron_sel[3] = true;
if (detain < detainCutEndcap_ && isEndcap) electron_sel[3] = true;
if (isBarrel)
{
LogTrace("") << "\t... sieie value " << sieie << " (barrel), pass? " << electron_sel[2];
LogTrace("") << "\t... detain value " << detain << " (barrel), pass? " << electron_sel[3];
}
else if (isEndcap)
{
LogTrace("") << "\t... sieie value " << sieie << " (endcap), pass? " << electron_sel[2];
LogTrace("") << "\t... detain value " << detain << " (endcap), pass? " << electron_sel[2];
}
if (isBarrel)
{
sieiebarrel_before_->Fill(sieie);
detainbarrel_before_->Fill(detain);
}
else if (isEndcap)
{
sieieendcap_before_->Fill(sieie);
detainendcap_before_->Fill(detain);
}
// Isolation cuts
//double isovar = mu.isolationR03().sumPt;
double ecalisovar = elec.dr03EcalRecHitSumEt(); // picked one set!
double hcalisovar = elec.dr03HcalTowerSumEt(); // try others if
double trkisovar = elec.dr04TkSumPt(); // doesn't work
//if (isCombinedIso_) {
//isovar += mu.isolationR03().emEt;
//isovar += mu.isolationR03().hadEt;
//}
//if (isRelativeIso_) isovar /= pt;
if (ecalisovar<ecalIsoCutBarrel_ && isBarrel) electron_sel[4] = true;
if (ecalisovar<ecalIsoCutEndcap_ && isEndcap) electron_sel[4] = true;
if (hcalisovar<hcalIsoCutBarrel_ && isBarrel) electron_sel[5] = true;
if (hcalisovar<hcalIsoCutEndcap_ && isEndcap) electron_sel[5] = true;
if (trkisovar<trkIsoCutBarrel_ && isBarrel) electron_sel[6] = true;
if (trkisovar<trkIsoCutEndcap_ && isEndcap) electron_sel[6] = true;
if (isBarrel)
{
LogTrace("") << "\t... ecal isolation value " << ecalisovar << " (barrel), pass? " << electron_sel[4];
LogTrace("") << "\t... hcal isolation value " << hcalisovar << " (barrel), pass? " << electron_sel[5];
LogTrace("") << "\t... trk isolation value " << trkisovar << " (barrel), pass? " << electron_sel[6];
}
else if (isEndcap)
{
LogTrace("") << "\t... ecal isolation value " << ecalisovar << " (endcap), pass? " << electron_sel[4];
LogTrace("") << "\t... hcal isolation value " << hcalisovar << " (endcap), pass? " << electron_sel[5];
LogTrace("") << "\t... trk isolation value " << trkisovar << " (endcap), pass? " << electron_sel[6];
}
//iso_before_->Fill(isovar);
if (isBarrel)
{
ecalisobarrel_before_->Fill(ecalisovar);
hcalisobarrel_before_->Fill(hcalisovar);
trkisobarrel_before_->Fill(trkisovar);
}
else if (isEndcap)
{
ecalisoendcap_before_->Fill(ecalisovar);
hcalisoendcap_before_->Fill(hcalisovar);
trkisoendcap_before_->Fill(trkisovar);
}
// HLT
if (trigger_fired) electron_sel[7] = true;
// // MET/MT cuts
double w_et = met_et+pt;
double w_px = met_px+px;
double w_py = met_py+py;
double massT = w_et*w_et - w_px*w_px - w_py*w_py;
massT = (massT>0) ? sqrt(massT) : 0;
LogTrace("") << "\t... W mass, W_et, W_px, W_py: " << massT << ", " << w_et << ", " << w_px << ", " << w_py << " [GeV]";
if (massT>mtMin_ && massT<mtMax_) electron_sel[8] = true;
mt_before_->Fill(massT);
if (met_et>metMin_ && met_et<metMax_) electron_sel[9] = true;
// // Acoplanarity cuts
// Geom::Phi<double> deltaphi(mu.phi()-atan2(met_py,met_px));
// double acop = deltaphi.value();
// if (acop<0) acop = - acop;
// acop = M_PI - acop;
// LogTrace("") << "\t... acoplanarity: " << acop;
// if (acop<acopCut_) muon_sel[10] = true;
// acop_before_->Fill(acop);
// // Remaining flags (from global event information)
// if (nmuonsForZ1<1 || nmuonsForZ2<2) muon_sel[11] = true;
if (njets<=nJetMax_) electron_sel[10] = true;
// Collect necessary flags "per electron"
int flags_passed = 0;
bool rec_sel_this = true;
bool eid_sel_this = true;
bool iso_sel_this = true;
bool all_sel_this = true;
for (int j=0; j<NFLAGS; ++j)
{
if (electron_sel[j]) flags_passed += 1;
if (j<2 && !electron_sel[j]) rec_sel_this = false;
if (j<4 && !electron_sel[j]) eid_sel_this = false;
if (j<7 && !electron_sel[j]) iso_sel_this = false;
if (!electron_sel[j]) all_sel_this = false;
}
if (all_sel_this)
{
if (nGoodElectrons < 2)
{
goodElectron[nGoodElectrons][0] = 1.;
goodElectron[nGoodElectrons][1] = elec.massSqr();
goodElectron[nGoodElectrons][2] = elec.energy();
goodElectron[nGoodElectrons][3] = elec.px();
goodElectron[nGoodElectrons][4] = elec.py();
goodElectron[nGoodElectrons][5] = elec.pz();
}
nGoodElectrons++;
}
// // "rec" => pt,eta and quality cuts are satisfied
// if (rec_sel_this) rec_sel = true;
// // "iso" => "rec" AND "muon is isolated"
// if (iso_sel_this) iso_sel = true;
// // "hlt" => "iso" AND "event is triggered"
// if (hlt_sel_this) hlt_sel = true;
// // "all" => "met" AND "Z/top rejection cuts"
// if (all_sel_this) all_sel = true;
// "rec" => pt,eta cuts are satisfied
if (rec_sel_this) rec_sel = true;
// "eid" => "rec" AND "electron passes ID"
if (eid_sel_this) iso_sel = true;
// "iso" => "eid" AND "electron is isolated"
if (iso_sel_this) iso_sel = true;
// "met" => "iso" AND "MET/MT"
// "all" => "met" AND "event is triggered"
if (all_sel_this) all_sel = true;
// Do N-1 histograms now (and only once for global event quantities)
if (flags_passed >= (NFLAGS-1))
{
if (!electron_sel[0] || flags_passed==NFLAGS)
{
pt_after_->Fill(pt);
}
if (!electron_sel[1] || flags_passed==NFLAGS)
{
eta_after_->Fill(eta);
}
if (!electron_sel[2] || flags_passed==NFLAGS)
{
if (isBarrel)
{
sieiebarrel_after_->Fill(sieie);
}
else if (isEndcap)
{
sieieendcap_after_->Fill(sieie);
}
}
if (!electron_sel[3] || flags_passed==NFLAGS)
{
if (isBarrel)
{
detainbarrel_after_->Fill(detain);
}
else if (isEndcap)
{
detainendcap_after_->Fill(detain);
}
}
if (!electron_sel[4] || flags_passed==NFLAGS)
{
if (isBarrel)
{
ecalisobarrel_after_->Fill(ecalisovar);
}
else if (isEndcap)
{
ecalisoendcap_after_->Fill(ecalisovar);
}
}
if (!electron_sel[5] || flags_passed==NFLAGS)
{
if (isBarrel)
{
hcalisobarrel_after_->Fill(hcalisovar);
}
else if (isEndcap)
{
hcalisoendcap_after_->Fill(hcalisovar);
}
}
if (!electron_sel[6] || flags_passed==NFLAGS)
{
if (isBarrel)
{
trkisobarrel_after_->Fill(trkisovar);
}
else if (isEndcap)
{
trkisoendcap_after_->Fill(trkisovar);
}
}
// if (!electron_sel[3] || flags_passed==NFLAGS)
// {
// detain_after_->Fill(detain);
// }
// if (!electron_sel[4] || flags_passed==NFLAGS)
// {
// ecaliso_after_->Fill(trackerHits);
// }
// if (!electron_sel[5] || flags_passed==NFLAGS)
// {
// tkelectr_after_->Fill(electr.isTrackerElectron());
// }
// if (!electron_sel[6] || flags_passed==NFLAGS)
// {
// iso_after_->Fill(isovar);
// }
if (!electron_sel[7] || flags_passed==NFLAGS)
{
if (!hlt_hist_done)
{
trig_after_->Fill(trigger_fired);
}
}
hlt_hist_done = true;
if (!electron_sel[8] || flags_passed==NFLAGS)
{
mt_after_->Fill(massT);
}
if (!electron_sel[9] || flags_passed==NFLAGS)
{
if (!met_hist_done) {
met_after_->Fill(met_et);
}
}
met_hist_done = true;
// if (!muon_sel[10] || flags_passed==NFLAGS)
// acop_after_->Fill(acop);
// if (!muon_sel[11] || flags_passed==NFLAGS)
// if (!nz1_hist_done) nz1_after_->Fill(nmuonsForZ1);
// nz1_hist_done = true;
// if (!muon_sel[11] || flags_passed==NFLAGS)
// if (!nz2_hist_done) nz2_after_->Fill(nmuonsForZ2);
// nz2_hist_done = true;
if (!electron_sel[10] || flags_passed==NFLAGS) {
if (!njets_hist_done)
{
njets_after_->Fill(njets);
if (jet_et>10) //don't want low energy "jets"
{
jet_et_after_ ->Fill(jet_et);
jet_eta_after_ ->Fill(jet_eta);
}
}
}
njets_hist_done = true;
} // end N-1 histos block
} // end loop through electrons
// inv mass = sqrt(m_1^2 + m_2^2 + 2*(E_1*E_2 - (px1*px2 + py1*py2 + pz1+pz2) ) )
double invMass=0;
nelectrons_before_->Fill(electronCollectionSize);
if (electronCollectionSize > 1)
{
invMass = sqrt(electron[0][1] + electron[1][1] + 2*(electron[0][2]*electron[1][2] - (electron[0][3]*electron[1][3] + electron[0][4]*electron[1][4] + electron[0][5]*electron[1][5]) ) );
invmass_before_->Fill(invMass);
invmassPU_before_->Fill(invMass,npvCount);
}
nelectrons_after_->Fill(nGoodElectrons);
if (nGoodElectrons > 1)
{
invMass = sqrt(goodElectron[0][1] + goodElectron[1][1] + 2*(goodElectron[0][2]*goodElectron[1][2] - (goodElectron[0][3]*goodElectron[1][3] + goodElectron[0][4]*goodElectron[1][4] + goodElectron[0][5]*goodElectron[1][5]) ) );
invmass_after_->Fill(invMass);
invmassPU_afterZ_->Fill(invMass,npvCount);
npvs_afterZ_->Fill(npvCount);
}
// Collect final flags
if (rec_sel) nrec++;
if (eid_sel) neid++;
if (iso_sel) niso++;
// if (hlt_sel) nhlt++;
// if (met_sel) nmet++;
if (all_sel) {
nsel++;
LogTrace("") << ">>>> Event ACCEPTED";
} else {
LogTrace("") << ">>>> Event REJECTED";
}
return;
}
| void EwkElecDQM::beginJob | ( | void | ) | [virtual] |
| void EwkElecDQM::beginRun | ( | const edm::Run & | iRun, |
| const edm::EventSetup & | iSet | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 106 of file EwkElecDQM.cc.
References hltConfigProvider_, HLTConfigProvider::init(), isValidHltConfig_, LogTrace, nall, neid, niso, nrec, and nsel.
{
nall = 0;
nsel = 0;
nrec = 0;
neid = 0;
niso = 0;
// nhlt = 0;
// nmet = 0;
// passed as parameter to HLTConfigProvider::init(), not yet used
bool isConfigChanged = false;
// isValidHltConfig_ could be used to short-circuit analyze() in case of problems
isValidHltConfig_ = hltConfigProvider_.init( iRun, iSet, "HLT", isConfigChanged );
LogTrace("") << "isValidHltConfig_=" << isValidHltConfig_ << "\n";
}
| double EwkElecDQM::calcDeltaPhi | ( | double | phi1, |
| double | phi2 | ||
| ) |
Definition at line 1023 of file EwkElecDQM.cc.
References SiPixelRawToDigiRegional_cfi::deltaPhi.
Referenced by analyze().
| void EwkElecDQM::endJob | ( | void | ) | [virtual] |
| void EwkElecDQM::endRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | |||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 276 of file EwkElecDQM.cc.
References cond::ecalcond::all, interpolateCardsSimple::eff, nall, neid, niso, nrec, nsel, and mathSSE::sqrt().
{
// overall
double all = nall;
double esel = nsel/all;
LogVerbatim("") << "\n>>>>>> SELECTION SUMMARY BEGIN >>>>>>>>>>>>>>>";
LogVerbatim("") << "Total number of events analyzed: " << nall << " [events]";
LogVerbatim("") << "Total number of events selected: " << nsel << " [events]";
LogVerbatim("") << "Overall efficiency: " << "(" << setprecision(4) << esel*100. <<" +/- "<< setprecision(2) << sqrt(esel*(1-esel)/all)*100. << ")%";
double erec = nrec/all;
double eeid = neid/all;
double eiso = niso/all;
// double ehlt = nhlt/all;
// double emet = nmet/all;
// general reconstruction step??
double num = nrec;
double eff = erec;
double err = sqrt(eff*(1-eff)/all);
LogVerbatim("") << "Passing Pt/Eta/Quality cuts: " << num << " [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) << err*100. << ")%";
// electron ID step
num = neid;
eff = eeid;
err = sqrt(eff*(1-eff)/all);
double effstep = 0.;
double errstep = 0.;
if (nrec>0) effstep = eeid/erec;
if (nrec>0) errstep = sqrt(effstep*(1-effstep)/nrec);
LogVerbatim("") << "Passing eID cuts: " << num << " [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) << err*100. << ")%, to previous step: (" << setprecision(4) << effstep*100. << " +/- "<< setprecision(2) << errstep*100. <<")%";
// isolation step
num = niso;
eff = eiso;
err = sqrt(eff*(1-eff)/all);
effstep = 0.;
errstep = 0.;
if (neid>0) effstep = eiso/eeid;
if (neid>0) errstep = sqrt(effstep*(1-effstep)/neid);
LogVerbatim("") << "Passing isolation cuts: " << num << " [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) << err*100. << ")%, to previous step: (" << setprecision(4) << effstep*100. << " +/- "<< setprecision(2) << errstep*100. <<")%";
// // trigger step
// num = nhlt;
// eff = ehlt;
// err = sqrt(eff*(1-eff)/all);
// effstep = 0.;
// errstep = 0.;
// if (niso>0) effstep = ehlt/eiso;
// if (niso>0) errstep = sqrt(effstep*(1-effstep)/niso);
// LogVerbatim("") << "Passing HLT criteria: " << num << " [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) << err*100. << ")%, to previous step: (" << setprecision(4) << effstep*100. << " +/- "<< setprecision(2) << errstep*100. <<")%";
// trigger step
num = nsel;
eff = esel;
err = sqrt(eff*(1-eff)/all);
effstep = 0.;
errstep = 0.;
if (niso>0) effstep = esel/eiso;
if (niso>0) errstep = sqrt(effstep*(1-effstep)/niso);
LogVerbatim("") << "Passing HLT criteria: " << num << " [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) << err*100. << ")%, to previous step: (" << setprecision(4) << effstep*100. << " +/- "<< setprecision(2) << errstep*100. <<")%";
// // met/acoplanarity cuts
// num = nmet;
// eff = emet;
// err = sqrt(eff*(1-eff)/all);
// effstep = 0.;
// errstep = 0.;
// if (nhlt>0) effstep = emet/ehlt;
// if (nhlt>0) errstep = sqrt(effstep*(1-effstep)/nhlt);
// LogVerbatim("") << "Passing MET/acoplanarity cuts: " << num << " [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) << err*100. << ")%, to previous step: (" << setprecision(4) << effstep*100. << " +/- "<< setprecision(2) << errstep*100. <<")%";
// // Z/top selection cuts ALSO LAST STEP so "sel" for "selection"
// num = nsel;
// eff = esel;
// err = sqrt(eff*(1-eff)/all);
// effstep = 0.;
// errstep = 0.;
// if (nmet>0) effstep = esel/emet;
// if (nmet>0) errstep = sqrt(effstep*(1-effstep)/nmet);
// LogVerbatim("") << "Passing Z/top rejection cuts: " << num << " [events], (" << setprecision(4) << eff*100. <<" +/- "<< setprecision(2) << err*100. << ")%, to previous step: (" << setprecision(4) << effstep*100. << " +/- "<< setprecision(2) << errstep*100. <<")%";
LogVerbatim("") << ">>>>>> SELECTION SUMMARY END >>>>>>>>>>>>>>>\n";
}
| void EwkElecDQM::init_histograms | ( | ) |
Definition at line 129 of file EwkElecDQM.cc.
References DQMStore::book1D(), DQMStore::book2D(), detainbarrel_after_, detainbarrel_before_, detainendcap_after_, detainendcap_before_, ecalisobarrel_after_, ecalisobarrel_before_, ecalisoendcap_after_, ecalisoendcap_before_, eJetMin_, eta_after_, eta_before_, hcalisobarrel_after_, hcalisobarrel_before_, hcalisoendcap_after_, hcalisoendcap_before_, i, invmass_after_, invmass_before_, invmassPU_afterZ_, invmassPU_before_, jet_et_after_, jet_et_before_, jet_eta_after_, jet_eta_before_, jetTag_, edm::InputTag::label(), met_after_, met_before_, metTag_, mt_after_, mt_before_, nelectrons_after_, nelectrons_before_, njets_after_, njets_before_, npvs_afterZ_, npvs_before_, pt_after_, pt_before_, PUBinCount_, PUMax_, sieiebarrel_after_, sieiebarrel_before_, sieieendcap_after_, sieieendcap_before_, theDbe, trig_after_, trig_before_, trkisobarrel_after_, trkisobarrel_before_, trkisoendcap_after_, and trkisoendcap_before_.
Referenced by EwkElecDQM().
{
char chtitle[256] = "";
for (int i=0; i<2; ++i) {
// pt_before_ = theDbe->book1D("PT_BEFORECUTS","Muon transverse momentum (global muon) [GeV],100,0.,100.);
// pt_after_ = theDbe->book1D("PT_LASTCUT","Muon transverse momentum (global muon) [GeV],100,0.,100.);
pt_before_ = theDbe->book1D("PT_BEFORECUTS","Electron transverse momentum [GeV]",100,0.,100.);
pt_after_ = theDbe->book1D("PT_LASTCUT","Electron transverse momentum [GeV]",100,0.,100.);
eta_before_ = theDbe->book1D("ETA_BEFORECUTS","Electron pseudo-rapidity",50,-2.5,2.5);
eta_after_ = theDbe->book1D("ETA_LASTCUT","Electron pseudo-rapidity",50,-2.5,2.5);
sieiebarrel_before_ = theDbe->book1D("SIEIEBARREL_BEFORECUTS","Electron #sigma_{i#etai#eta} (barrel)",70,0.,0.07);
sieiebarrel_after_ = theDbe->book1D("SIEIEBARREL_LASTCUT","Electron #sigma_{i#etai#eta} (barrel)",70,0.,0.07);
sieieendcap_before_ = theDbe->book1D("SIEIEENDCAP_BEFORECUTS","Electron #sigma_{i#etai#eta} (endcap)",70,0.,0.07);
sieieendcap_after_ = theDbe->book1D("SIEIEENDCAP_LASTCUT","Electron #sigma_{i#etai#eta} (endcap)",70,0.,0.07);
detainbarrel_before_ = theDbe->book1D("DETAINBARREL_BEFORECUTS","Electron #Delta#eta_{in} (barrel)",40,-0.02,0.02);
detainbarrel_after_ = theDbe->book1D("DETAINBARREL_LASTCUT","Electron #Delta#eta_{in} (barrel)",40,-0.02,0.02);
detainendcap_before_ = theDbe->book1D("DETAINENDCAP_BEFORECUTS","Electron #Delta#eta_{in} (endcap)",40,-0.02,0.02);
detainendcap_after_ = theDbe->book1D("DETAINENDCAP_LASTCUT","Electron #Delta#eta_{in} (endcap)",40,-0.02,0.02);
// dxy_before_ = theDbe->book1D("DXY_BEFORECUTS","Muon transverse distance to beam spot [cm]",100,-0.5,0.5);
// dxy_after_ = theDbe->book1D("DXY_LASTCUT","Muon transverse distance to beam spot [cm]",100,-0.5,0.5);
// chi2_before_ = theDbe->book1D("CHI2_BEFORECUTS","Normalized Chi2, inner track fit",100,0.,100.);
// chi2_after_ = theDbe->book1D("CHI2_LASTCUT","Normalized Chi2, inner track fit",100,0.,100.);
// nhits_before_ = theDbe->book1D("NHITS_BEFORECUTS","Number of hits, inner track",40,-0.5,39.5);
// nhits_after_ = theDbe->book1D("NHITS_LASTCUT","Number of hits, inner track",40,-0.5,39.5);
// tkmu_before_ = theDbe->book1D("TKMU_BEFORECUTS","Tracker-muon flag (for global muons)",2,-0.5,1.5);
// tkmu_after_ = theDbe->book1D("TKMU_LASTCUT","Tracker-muon flag (for global muons)",2,-0.5,1.5);
// if (isRelativeIso_) {
// if (isCombinedIso_) {
// snprintf(chtitle, 255, "Relative (combined) isolation variable");
// iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Relative (combined) isolation variable",100, 0., 1.);
// iso_after_ = theDbe->book1D("ISO_LASTCUT","Relative (combined) isolation variable",100, 0., 1.);
// } else {
// snprintf(chtitle, 255, "Relative (tracker) isolation variable");
// iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Relative (tracker) isolation variable",100, 0., 1.);
// iso_after_ = theDbe->book1D("ISO_LASTCUT","Relative (tracker) isolation variable",100, 0., 1.);
// }
// } else {
// if (isCombinedIso_) {
// iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Absolute (combined) isolation variable [GeV]",100, 0., 20.);
// iso_after_ = theDbe->book1D("ISO_LASTCUT","Absolute (combined) isolation variable [GeV]",100, 0., 20.);
// } else {
// iso_before_ = theDbe->book1D("ISO_BEFORECUTS","Absolute (tracker) isolation variable [GeV]",100, 0., 20.);
// iso_after_ = theDbe->book1D("ISO_LASTCUT","Absolute (tracker) isolation variable [GeV]",100, 0., 20.);
// }
// }
ecalisobarrel_before_ = theDbe->book1D("ECALISOBARREL_BEFORECUTS","Absolute electron ECAL isolation variable (barrel) [GeV]",50,0.,50.);
ecalisobarrel_after_ = theDbe->book1D("ECALISOBARREL_LASTCUT","Absolute electron ECAL isolation variable (barrel) [GeV]",50,0.,50.);
ecalisoendcap_before_ = theDbe->book1D("ECALISOENDCAP_BEFORECUTS","Absolute electron ECAL isolation variable (endcap) [GeV]",50,0.,50.);
ecalisoendcap_after_ = theDbe->book1D("ECALISOENDCAP_LASTCUT","Absolute electron ECAL isolation variable (endcap) [GeV]",50,0.,50.);
hcalisobarrel_before_ = theDbe->book1D("HCALISOBARREL_BEFORECUTS","Absolute electron HCAL isolation variable (barrel) [GeV]",50,0.,50.);
hcalisobarrel_after_ = theDbe->book1D("HCALISOBARREL_LASTCUT","Absolute electron HCAL isolation variable (barrel) [GeV]",50,0.,50.);
hcalisoendcap_before_ = theDbe->book1D("HCALISOENDCAP_BEFORECUTS","Absolute electron HCAL isolation variable (endcap) [GeV]",50,0.,50.);
hcalisoendcap_after_ = theDbe->book1D("HCALISOENDCAP_LASTCUT","Absolute electron HCAL isolation variable (endcap) [GeV]",50,0.,50.);
trkisobarrel_before_ = theDbe->book1D("TRKISOBARREL_BEFORECUTS","Absolute electron track isolation variable (barrel) [GeV]",50,0.,50.);
trkisobarrel_after_ = theDbe->book1D("TRKISOBARREL_LASTCUT","Absolute electron track isolation variable (barrel) [GeV]",50,0.,50.);
trkisoendcap_before_ = theDbe->book1D("TRKISOENDCAP_BEFORECUTS","Absolute electron track isolation variable (endcap) [GeV]",50,0.,50.);
trkisoendcap_after_ = theDbe->book1D("TRKISOENDCAP_LASTCUT","Absolute electron track isolation variable (endcap) [GeV]",50,0.,50.);
// snprintf(chtitle, 255, "Trigger response (bit %s)", muonTrig_.data());
// trig_before_ = theDbe->book1D("TRIG_BEFORECUTS",chtitle,2,-0.5,1.5);
// trig_after_ = theDbe->book1D("TRIG_LASTCUT",chtitle,2,-0.5,1.5);
//snprintf(chtitle, 255, "Trigger response (bit %s)", elecTrig_.data());
trig_before_ = theDbe->book1D("TRIG_BEFORECUTS","Trigger response",2,-0.5,1.5); // elecTrig_ is now a vector of strings!
trig_after_ = theDbe->book1D("TRIG_LASTCUT","Trigger response",2,-0.5,1.5);
invmass_before_ = theDbe->book1D("INVMASS_BEFORECUTS","Di-electron invariant mass [GeV]",100,0.,200.);
invmass_after_ = theDbe->book1D("INVMASS_AFTERCUTS","Di-electron invariant mass [GeV]",100,0.,200.);
invmassPU_before_ = theDbe->book2D("INVMASS_PU_BEFORECUTS","Di-electron invariant mass [GeV] vs PU; mass [GeV]; PU count",100,0.,200., PUBinCount_, -0.5, PUMax_+0.5);
invmassPU_afterZ_ = theDbe->book2D("INVMASS_PU_AFTERZCUTS","Di-electron invariant mass [GeV] vs PU; mass [GeV]; PU count",100,0.,200., PUBinCount_, -0.5, PUMax_+0.5);
npvs_before_ = theDbe->book1D("NPVs_BEFORECUTS","Number of Valid Primary Vertices; nGoodPVs",PUMax_+1,-0.5,PUMax_+0.5);
//npvs_afterW_ = theDbe->book1D("NPVs_AFTERWCUTS","Number of Valid Primary Vertices",PUMax_+1,-0.5,PUMax_+0.5);
npvs_afterZ_ = theDbe->book1D("NPVs_AFTERZCUTS","Number of Valid Primary Vertices; nGoodPVs",PUMax_+1,-0.5,PUMax_+0.5);
nelectrons_before_ = theDbe->book1D("NELECTRONS_BEFORECUTS","Number of electrons in event",10,-0.5,9.5);
nelectrons_after_ = theDbe->book1D("NELECTRONS_AFTERCUTS","Number of electrons in event",10,-0.5,9.5);
snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
mt_before_ = theDbe->book1D("MT_BEFORECUTS",chtitle,150,0.,300.);
mt_after_ = theDbe->book1D("MT_LASTCUT",chtitle,150,0.,300.);
// snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
// mt_before_ = theDbe->book1D("MT_BEFORECUTS",chtitle,150,0.,300.);
// mt_after_ = theDbe->book1D("MT_LASTCUT",chtitle,150,0.,300.);
snprintf(chtitle, 255, "Missing transverse energy (%s) [GeV]", metTag_.label().data());
met_before_ = theDbe->book1D("MET_BEFORECUTS",chtitle,100,0.,200.);
met_after_ = theDbe->book1D("MET_LASTCUT",chtitle,100,0.,200.);
// snprintf(chtitle, 255, "MU-MET (%s) acoplanarity", metTag_.label().data());
// acop_before_ = theDbe->book1D("ACOP_BEFORECUTS",chtitle,50,0.,M_PI);
// acop_after_ = theDbe->book1D("ACOP_LASTCUT",chtitle,50,0.,M_PI);
// snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ1_);
// nz1_before_ = theDbe->book1D("NZ1_BEFORECUTS",chtitle,10,-0.5,9.5);
// nz1_after_ = theDbe->book1D("NZ1_LASTCUT",chtitle,10,-0.5,9.5);
// snprintf(chtitle, 255, "Z rejection: number of muons above %.2f GeV", ptThrForZ2_);
// nz2_before_ = theDbe->book1D("NZ2_BEFORECUTS",chtitle,10,-0.5,9.5);
// nz2_after_ = theDbe->book1D("NZ2_LASTCUT",chtitle,10,-0.5,9.5);
snprintf(chtitle, 255, "Number of jets (%s) above %.2f GeV", jetTag_.label().data(), eJetMin_);
njets_before_ = theDbe->book1D("NJETS_BEFORECUTS",chtitle,10,-0.5,9.5);
njets_after_ = theDbe->book1D("NJETS_LASTCUT",chtitle,10,-0.5,9.5);
snprintf(chtitle, 255, "Jet with highest E_{T} (%s)", jetTag_.label().data());
jet_et_before_ = theDbe->book1D("JETET1_BEFORECUTS",chtitle, 20, 0., 200.0);
jet_et_after_ = theDbe->book1D("JETET1_AFTERCUTS",chtitle, 20, 0., 200.0);
snprintf(chtitle, 255, "Eta of Jet with highest E_{T} (%s)", jetTag_.label().data());
jet_eta_before_ = theDbe->book1D("JETETA1_BEFORECUTS",chtitle, 20, -5, 5);
jet_eta_after_ = theDbe->book1D("JETETA1_AFTERCUTS",chtitle, 20, -5, 5);
// snprintf(chtitle, 255, "Jet with 2nd highest E_{T} (%s)", jetTag_.label().data());
// jet2_et_before = theDbe->book1D("JETET2_BEFORECUTS",chtitle, 20, 0., 200.0);
// jet2_et_after = theDbe->book1D("JETET2_AFTERCUTS",chtitle, 20, 0., 200.0);
}
}
MonitorElement* EwkElecDQM::detainbarrel_after_ [private] |
Definition at line 110 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::detainbarrel_before_ [private] |
Definition at line 109 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::detainCutBarrel_ [private] |
Definition at line 48 of file EwkElecDQM.h.
Referenced by analyze().
double EwkElecDQM::detainCutEndcap_ [private] |
Definition at line 49 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::detainendcap_after_ [private] |
Definition at line 113 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::detainendcap_before_ [private] |
Definition at line 112 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::ecalisobarrel_after_ [private] |
Definition at line 128 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::ecalisobarrel_before_ [private] |
Definition at line 127 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::ecalIsoCutBarrel_ [private] |
Definition at line 55 of file EwkElecDQM.h.
Referenced by analyze().
double EwkElecDQM::ecalIsoCutEndcap_ [private] |
Definition at line 56 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::ecalisoendcap_after_ [private] |
Definition at line 131 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::ecalisoendcap_before_ [private] |
Definition at line 130 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::eJetMin_ [private] |
Definition at line 75 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
edm::InputTag EwkElecDQM::elecTag_ [private] |
Definition at line 35 of file EwkElecDQM.h.
Referenced by analyze().
const std::vector<std::string> EwkElecDQM::elecTrig_ [private] |
Definition at line 42 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::eta_after_ [private] |
Definition at line 101 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::eta_before_ [private] |
Definition at line 100 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::etaCut_ [private] |
Definition at line 44 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::hcalisobarrel_after_ [private] |
Definition at line 134 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::hcalisobarrel_before_ [private] |
Definition at line 133 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::hcalIsoCutBarrel_ [private] |
Definition at line 57 of file EwkElecDQM.h.
Referenced by analyze().
double EwkElecDQM::hcalIsoCutEndcap_ [private] |
Definition at line 58 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::hcalisoendcap_after_ [private] |
Definition at line 137 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::hcalisoendcap_before_ [private] |
Definition at line 136 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
Definition at line 82 of file EwkElecDQM.h.
Referenced by analyze(), and beginRun().
MonitorElement* EwkElecDQM::invmass_after_ [private] |
Definition at line 149 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::invmass_before_ [private] |
Definition at line 148 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::invmassPU_afterZ_ [private] |
Definition at line 151 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::invmassPU_before_ [private] |
Definition at line 150 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
bool EwkElecDQM::isValidHltConfig_ [private] |
Definition at line 81 of file EwkElecDQM.h.
Referenced by beginRun(), and EwkElecDQM().
MonitorElement* EwkElecDQM::jet_et_after_ [private] |
Definition at line 178 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::jet_et_before_ [private] |
Definition at line 177 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::jet_eta_after_ [private] |
Definition at line 180 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::jet_eta_before_ [private] |
Definition at line 179 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
edm::InputTag EwkElecDQM::jetTag_ [private] |
Definition at line 38 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::met_after_ [private] |
Definition at line 164 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::met_before_ [private] |
Definition at line 163 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
bool EwkElecDQM::metIncludesMuons_ [private] |
Definition at line 37 of file EwkElecDQM.h.
double EwkElecDQM::metMax_ [private] |
Definition at line 64 of file EwkElecDQM.h.
Referenced by analyze().
double EwkElecDQM::metMin_ [private] |
Definition at line 63 of file EwkElecDQM.h.
Referenced by analyze().
edm::InputTag EwkElecDQM::metTag_ [private] |
Definition at line 36 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::mt_after_ [private] |
Definition at line 161 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::mt_before_ [private] |
Definition at line 160 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::mtMax_ [private] |
Definition at line 62 of file EwkElecDQM.h.
Referenced by analyze().
double EwkElecDQM::mtMin_ [private] |
Definition at line 61 of file EwkElecDQM.h.
Referenced by analyze().
unsigned int EwkElecDQM::nall [private] |
Definition at line 84 of file EwkElecDQM.h.
Referenced by analyze(), beginRun(), and endRun().
unsigned int EwkElecDQM::neid [private] |
Definition at line 86 of file EwkElecDQM.h.
Referenced by analyze(), beginRun(), and endRun().
MonitorElement* EwkElecDQM::nelectrons_after_ [private] |
Definition at line 158 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::nelectrons_before_ [private] |
Definition at line 157 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
unsigned int EwkElecDQM::nGoodElectrons [private] |
Definition at line 93 of file EwkElecDQM.h.
Referenced by analyze().
unsigned int EwkElecDQM::niso [private] |
Definition at line 87 of file EwkElecDQM.h.
Referenced by analyze(), beginRun(), and endRun().
int EwkElecDQM::nJetMax_ [private] |
Definition at line 76 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::njets_after_ [private] |
Definition at line 176 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::njets_before_ [private] |
Definition at line 175 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::npvs_afterZ_ [private] |
Definition at line 155 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::npvs_before_ [private] |
Definition at line 153 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
unsigned int EwkElecDQM::nrec [private] |
Definition at line 85 of file EwkElecDQM.h.
Referenced by analyze(), beginRun(), and endRun().
unsigned int EwkElecDQM::nsel [private] |
Definition at line 90 of file EwkElecDQM.h.
Referenced by analyze(), beginRun(), and endRun().
MonitorElement* EwkElecDQM::pt_after_ [private] |
Definition at line 98 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::pt_before_ [private] |
Definition at line 97 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::ptCut_ [private] |
Definition at line 43 of file EwkElecDQM.h.
Referenced by analyze().
unsigned int EwkElecDQM::PUBinCount_ [private] |
Definition at line 79 of file EwkElecDQM.h.
Referenced by init_histograms().
unsigned int EwkElecDQM::PUMax_ [private] |
Definition at line 79 of file EwkElecDQM.h.
Referenced by init_histograms().
MonitorElement* EwkElecDQM::sieiebarrel_after_ [private] |
Definition at line 104 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::sieiebarrel_before_ [private] |
Definition at line 103 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::sieieCutBarrel_ [private] |
Definition at line 46 of file EwkElecDQM.h.
Referenced by analyze().
double EwkElecDQM::sieieCutEndcap_ [private] |
Definition at line 47 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::sieieendcap_after_ [private] |
Definition at line 107 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::sieieendcap_before_ [private] |
Definition at line 106 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
DQMStore* EwkElecDQM::theDbe [private] |
Definition at line 95 of file EwkElecDQM.h.
Referenced by EwkElecDQM(), and init_histograms().
MonitorElement* EwkElecDQM::trig_after_ [private] |
Definition at line 146 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::trig_before_ [private] |
Definition at line 145 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
edm::InputTag EwkElecDQM::trigTag_ [private] |
Definition at line 33 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::trkisobarrel_after_ [private] |
Definition at line 140 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::trkisobarrel_before_ [private] |
Definition at line 139 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
double EwkElecDQM::trkIsoCutBarrel_ [private] |
Definition at line 59 of file EwkElecDQM.h.
Referenced by analyze().
double EwkElecDQM::trkIsoCutEndcap_ [private] |
Definition at line 60 of file EwkElecDQM.h.
Referenced by analyze().
MonitorElement* EwkElecDQM::trkisoendcap_after_ [private] |
Definition at line 143 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
MonitorElement* EwkElecDQM::trkisoendcap_before_ [private] |
Definition at line 142 of file EwkElecDQM.h.
Referenced by analyze(), and init_histograms().
edm::InputTag EwkElecDQM::vertexTag_ [private] |
Definition at line 39 of file EwkElecDQM.h.
Referenced by analyze().
1.7.3