|
|
|
#include <JetAnalyzer.h>
DQM monitoring source for Calo Jets
Definition at line 41 of file JetAnalyzer.h.
| JetAnalyzer::JetAnalyzer | ( | const edm::ParameterSet & | pSet | ) |
Constructor.
Definition at line 22 of file JetAnalyzer.cc.
References Parameters::parameters.
{
parameters = pSet;
_leadJetFlag = 0;
_JetLoPass = 0;
_JetHiPass = 0;
_ptThreshold = 20.;
_asymmetryThirdJetCut = 5.;
_balanceThirdJetCut = 0.2;
_n90HitsMin =0;
_fHPDMax=1.;
_resEMFMin=0.;
_n90HitsMinLoose =0;
_fHPDMaxLoose=1.;
_resEMFMinLoose=0.;
_n90HitsMinTight =0;
_fHPDMaxTight=1.;
_resEMFMinTight=0.;
_sigmaEtaMinTight=-999.;
_sigmaPhiMinTight=-999.;
}
| JetAnalyzer::~JetAnalyzer | ( | ) | [virtual] |
| void JetAnalyzer::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const reco::CaloJetCollection & | caloJets, | ||
| const int | numPV | ||
| ) |
Get the analysis.
Definition at line 273 of file JetAnalyzer.cc.
References corr, diffTreeTool::diff, dPhi(), eta, edm::EventID::event(), edm::EventBase::id(), metsig::jet, bTagSequences_cff::jetID, LogTrace, phi, and mathSSE::sqrt().
{
int numofjets=0;
double fstPhi=0.;
double sndPhi=0.;
double diff = 0.;
double corr = 0.;
double dphi = -999. ;
bool thiscleaned=false;
bool Loosecleaned=false;
bool Tightcleaned=false;
bool thisemfclean=true;
bool emfcleanLoose=true;
bool emfcleanTight=true;
srand( iEvent.id().event() % 10000);
if (makedijetselection == 1) {
//Dijet selection - careful: the pT is uncorrected!
//if(makedijetselection==1 && caloJets.size()>=2){
if(caloJets.size()>=2){
double dphiDJ = -999. ;
bool emfcleanLooseFirstJet=true;
bool emfcleanLooseSecondJet=true;
bool emfcleanTightFirstJet=true;
bool emfcleanTightSecondJet=true;
bool LoosecleanedFirstJet = false;
bool LoosecleanedSecondJet = false;
bool TightcleanedFirstJet = false;
bool TightcleanedSecondJet = false;
//both jets pass pt threshold
if ((caloJets.at(0)).pt() > _ptThreshold && (caloJets.at(1)).pt() > _ptThreshold ) {
if(fabs((caloJets.at(0)).eta())<3. && fabs((caloJets.at(1)).eta())<3. ){
//calculate dphi
dphiDJ = fabs((caloJets.at(0)).phi()-(caloJets.at(1)).phi());
if (dphiDJ > 3.14) dphiDJ=fabs(dphiDJ -6.28 );
//fill DPhi histo (before cutting)
if (mDPhi) mDPhi->Fill (dphiDJ);
//dphi cut
if(fabs(dphiDJ)>2.1){
//JetID
emfcleanLooseFirstJet=true;
emfcleanTightFirstJet=true;
emfcleanLooseSecondJet=true;
emfcleanTightSecondJet=true;
//jetID for first jet
jetID->calculate(iEvent, (caloJets.at(0)));
if(jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(0)).eta())<2.6) emfcleanLooseFirstJet=false;
if(jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(0)).eta())<2.6) emfcleanTightFirstJet=false;
if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseFirstJet) LoosecleanedFirstJet=true;
if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt((caloJets.at(0)).etaetaMoment())>_sigmaEtaMinTight && sqrt((caloJets.at(0)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightFirstJet) TightcleanedFirstJet=true;
//fill the JID variables histograms BEFORE you cut on them
if (mN90Hits) mN90Hits->Fill (jetID->n90Hits());
if (mfHPD) mfHPD->Fill (jetID->fHPD());
if (mresEMF) mresEMF->Fill (jetID->restrictedEMF());
if (mfRBX) mfRBX->Fill (jetID->fRBX());
//jetID for second jet
jetID->calculate(iEvent, (caloJets.at(1)));
if(jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(1)).eta())<2.6) emfcleanLooseSecondJet=false;
if(jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(1)).eta())<2.6) emfcleanTightSecondJet=false;
if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseSecondJet) LoosecleanedSecondJet=true;
if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt((caloJets.at(1)).etaetaMoment())>_sigmaEtaMinTight && sqrt((caloJets.at(1)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightSecondJet) TightcleanedSecondJet=true;
//fill the JID variables histograms BEFORE you cut on them
if (mN90Hits) mN90Hits->Fill (jetID->n90Hits());
if (mfHPD) mfHPD->Fill (jetID->fHPD());
if (mresEMF) mresEMF->Fill (jetID->restrictedEMF());
if (mfRBX) mfRBX->Fill (jetID->fRBX());
if(LoosecleanedFirstJet && LoosecleanedSecondJet) { //only if both jets are (loose) cleaned
//fill histos for first jet
if (mPt) mPt->Fill ((caloJets.at(0)).pt());
if (mEta) mEta->Fill ((caloJets.at(0)).eta());
if (mPhi) mPhi->Fill ((caloJets.at(0)).phi());
if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(0)).eta(),(caloJets.at(0)).phi());
if (mConstituents) mConstituents->Fill ((caloJets.at(0)).nConstituents());
if (mHFrac) mHFrac->Fill ((caloJets.at(0)).energyFractionHadronic());
if (mEFrac) mEFrac->Fill ((caloJets.at(0)).emEnergyFraction());
//if (mE) mE->Fill ((caloJets.at(0)).energy());
//if (mP) mP->Fill ((caloJets.at(0)).p());
//if (mMass) mMass->Fill ((caloJets.at(0)).mass());
if (mMaxEInEmTowers) mMaxEInEmTowers->Fill ((caloJets.at(0)).maxEInEmTowers());
if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(0)).maxEInHadTowers());
if (mN90Hits) mN90Hits->Fill (jetID->n90Hits());
if (mfHPD) mfHPD->Fill (jetID->fHPD());
if (mresEMF) mresEMF->Fill (jetID->restrictedEMF());
if (mfRBX) mfRBX->Fill (jetID->fRBX());
//sigmaeta and sigmaphi only used in the tight selection.
//fill the histos for them AFTER the loose selection
// if (msigmaEta) msigmaEta->Fill(sqrt((caloJets.at(0)).etaetaMoment()));
// if (msigmaPhi) msigmaPhi->Fill(sqrt((caloJets.at(0)).phiphiMoment()));
//fill histos for second jet
if (mPt) mPt->Fill ((caloJets.at(1)).pt());
if (mEta) mEta->Fill ((caloJets.at(1)).eta());
if (mPhi) mPhi->Fill ((caloJets.at(1)).phi());
if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(1)).eta(),(caloJets.at(1)).phi());
if (mConstituents) mConstituents->Fill ((caloJets.at(1)).nConstituents());
if (mHFrac) mHFrac->Fill ((caloJets.at(1)).energyFractionHadronic());
if (mEFrac) mEFrac->Fill ((caloJets.at(1)).emEnergyFraction());
//if (mE) mE->Fill ((caloJets.at(1)).energy());
//if (mP) mP->Fill ((caloJets.at(1)).p());
//if (mMass) mMass->Fill ((caloJets.at(1)).mass());
if (mMaxEInEmTowers) mMaxEInEmTowers->Fill ((caloJets.at(1)).maxEInEmTowers());
if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(1)).maxEInHadTowers());
//sigmaeta and sigmaphi only used in the tight selection.
//fill the histos for them AFTER the loose selection
// if (msigmaEta) msigmaEta->Fill(sqrt((caloJets.at(1)).etaetaMoment()));
// if (msigmaPhi) msigmaPhi->Fill(sqrt((caloJets.at(1)).phiphiMoment()));
// Fill NPV profiles
//----------------------------------------------------------------
for (int ijet=0; ijet<2; ijet++) {
if (mPt_profile) mPt_profile ->Fill(numPV, (caloJets.at(ijet)).pt());
if (mEta_profile) mEta_profile ->Fill(numPV, (caloJets.at(ijet)).eta());
if (mPhi_profile) mPhi_profile ->Fill(numPV, (caloJets.at(ijet)).phi());
if (mConstituents_profile) mConstituents_profile->Fill(numPV, (caloJets.at(ijet)).nConstituents());
if (mHFrac_profile) mHFrac_profile ->Fill(numPV, (caloJets.at(ijet)).energyFractionHadronic());
if (mEFrac_profile) mEFrac_profile ->Fill(numPV, (caloJets.at(ijet)).emEnergyFraction());
}
}
//let's see how many of these jets passed the JetID cleaning
if(fillJIDPassFrac==1) {
if(LoosecleanedFirstJet) {
mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),1.);
mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.);
} else {
mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),0.);
mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.);
}
if(LoosecleanedSecondJet) {
mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),1.);
mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.);
} else {
mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),0.);
mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.);
}
//TIGHT JID
if(TightcleanedFirstJet) {
mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),1.);
mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.);
} else {
mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),0.);
mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.);
}
if(TightcleanedSecondJet) {
mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),1.);
mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.);
} else {
mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),0.);
mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.);
}
}//if fillJIDPassFrac
}// FABS DPHI < 2.1
}// fabs eta < 3
}// pt jets > threshold
//now do the dijet balance and asymmetry calculations
if (fabs(caloJets.at(0).eta() < 1.4)) {
double pt_dijet = (caloJets.at(0).pt() + caloJets.at(1).pt())/2;
double dPhi = fabs((caloJets.at(0)).phi()-(caloJets.at(1)).phi());
if (dPhi > 3.14) dPhi=fabs(dPhi -6.28 );
if (dPhi > 2.7) {
double pt_probe;
double pt_barrel;
int jet1, jet2;
int randJet = rand() % 2;
if (fabs(caloJets.at(1).eta() < 1.4)) {
if (randJet) {
jet1 = 0;
jet2 = 1;
}
else {
jet1 = 1;
jet2 = 0;
}
/***Di-Jet Asymmetry****
* leading jets eta < 1.4
* leading jets dphi > 2.7
* pt_third jet < threshold
* A = (pt_1 - pt_2)/(pt_1 + pt_2)
* jets 1 and two are randomly ordered
*/
bool thirdJetCut = true;
for (unsigned int third = 2; third < caloJets.size(); ++third)
if (caloJets.at(third).pt() > _asymmetryThirdJetCut)
thirdJetCut = false;
if (thirdJetCut) {
double dijetAsymmetry = (caloJets.at(jet1).pt() - caloJets.at(jet2).pt()) / (caloJets.at(jet1).pt() + caloJets.at(jet2).pt());
mDijetAsymmetry->Fill(dijetAsymmetry);
}// end restriction on third jet pt in asymmetry calculation
}
else {
jet1 = 0;
jet2 = 1;
}
pt_barrel = caloJets.at(jet1).pt();
pt_probe = caloJets.at(jet2).pt();
//dijet balance cuts
/***Di-Jet Balance****
* pt_dijet = (pt_probe+pt_barrel)/2
* leading jets dphi > 2.7
* reject evnets where pt_third/pt_dijet > 0.2
* pv selection
* B = (pt_probe - pt_barrel)/pt_dijet
* select probe randomly from 2 jets if both leading jets are in the barrel
*/
bool thirdJetCut = true;
for (unsigned int third = 2; third < caloJets.size(); ++third)
if (caloJets.at(third).pt()/pt_dijet > _balanceThirdJetCut)
thirdJetCut = false;
if (thirdJetCut) {
double dijetBalance = (pt_probe - pt_barrel) / pt_dijet;
mDijetBalance->Fill(dijetBalance);
}// end restriction on third jet pt ratio in balance calculation
}// dPhi > 2.7
}// leading jet eta cut for asymmetry and balance calculations
}//jet size >= 2
}// do dijet selection
else {
for (reco::CaloJetCollection::const_iterator jet = caloJets.begin(); jet!=caloJets.end(); ++jet) {
LogTrace(jetname)<<"[JetAnalyzer] Analyze Calo Jet";
Loosecleaned=false;
Tightcleaned=false;
if (jet == caloJets.begin()) {
fstPhi = jet->phi();
_leadJetFlag = 1;
} else {
_leadJetFlag = 0;
}
if (jet == (caloJets.begin()+1)) sndPhi = jet->phi();
//jetID
jetID->calculate(iEvent, *jet);
//minimal (uncorrected!) pT cut
if (jet->pt() > _ptThreshold) {
// if (msigmaEta) msigmaEta->Fill(sqrt(jet->etaetaMoment()));
// if (msigmaPhi) msigmaPhi->Fill(sqrt(jet->phiphiMoment()));
//cleaning to use for filling histograms
thisemfclean=true;
if(jetID->restrictedEMF()<_resEMFMin && fabs(jet->eta())<2.6) thisemfclean=false;
if(jetID->n90Hits()>=_n90HitsMin && jetID->fHPD()<_fHPDMax && thisemfclean) thiscleaned=true;
//loose and tight cleaning, used to fill the JetIDPAssFraction histos
if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLoose) Loosecleaned=true;
if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt(jet->etaetaMoment())>_sigmaEtaMinTight && sqrt(jet->phiphiMoment())>_sigmaPhiMinTight && emfcleanTight) Tightcleaned=true;
if(fillJIDPassFrac==1) {
if(Loosecleaned) {
mLooseJIDPassFractionVSeta->Fill(jet->eta(),1.);
mLooseJIDPassFractionVSpt->Fill(jet->pt(),1.);
} else {
mLooseJIDPassFractionVSeta->Fill(jet->eta(),0.);
mLooseJIDPassFractionVSpt->Fill(jet->pt(),0.);
}
//TIGHT
if(Tightcleaned) {
mTightJIDPassFractionVSeta->Fill(jet->eta(),1.);
mTightJIDPassFractionVSpt->Fill(jet->pt(),1.);
} else {
mTightJIDPassFractionVSeta->Fill(jet->eta(),0.);
mTightJIDPassFractionVSpt->Fill(jet->pt(),0.);
}
}
//eventually we could define the "cleaned" flag differently for e.g. HF
if(thiscleaned) {
numofjets++ ;
jetME->Fill(1);
// Leading jet
// Histograms are filled once per event
if (_leadJetFlag == 1) {
if (mEtaFirst) mEtaFirst->Fill (jet->eta());
if (mPhiFirst) mPhiFirst->Fill (jet->phi());
//if (mEFirst) mEFirst->Fill (jet->energy());
if (mPtFirst) mPtFirst->Fill (jet->pt());
}
// --- Passed the low pt jet trigger
if (_JetLoPass == 1) {
/* if (fabs(jet->eta()) <= 1.3) {
if (mPt_Barrel_Lo) mPt_Barrel_Lo->Fill(jet->pt());
if (mEta_Lo) mEta_Lo->Fill(jet->eta());
if (mPhi_Barrel_Lo) mPhi_Barrel_Lo->Fill(jet->phi());
}
if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) {
if (mPt_EndCap_Lo) mPt_EndCap_Lo->Fill(jet->pt());
if (mEta_Lo) mEta_Lo->Fill(jet->eta());
if (mPhi_EndCap_Lo) mPhi_EndCap_Lo->Fill(jet->phi());
}
if (fabs(jet->eta()) > 3.0) {
if (mPt_Forward_Lo) mPt_Forward_Lo->Fill(jet->pt());
if (mEta_Lo) mEta_Lo->Fill(jet->eta());
if (mPhi_Forward_Lo) mPhi_Forward_Lo->Fill(jet->phi());
} */
//if (mEta_Lo) mEta_Lo->Fill (jet->eta());
if (mPhi_Lo) mPhi_Lo->Fill (jet->phi());
if (mPt_Lo) mPt_Lo->Fill (jet->pt());
}
// --- Passed the high pt jet trigger
if (_JetHiPass == 1) {
if (fabs(jet->eta()) <= 1.3) {
if (mPt_Barrel_Hi && jet->pt()>100.) mPt_Barrel_Hi->Fill(jet->pt());
if (mEta_Hi && jet->pt()>100.) mEta_Hi->Fill(jet->eta());
if (mPhi_Barrel_Hi) mPhi_Barrel_Hi->Fill(jet->phi());
//if (mConstituents_Barrel_Hi) mConstituents_Barrel_Hi->Fill(jet->nConstituents());
//if (mHFrac_Barrel_Hi) mHFrac_Barrel_Hi->Fill(jet->energyFractionHadronic());
}
if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) {
if (mPt_EndCap_Hi && jet->pt()>100.) mPt_EndCap_Hi->Fill(jet->pt());
if (mEta_Hi && jet->pt()>100.) mEta_Hi->Fill(jet->eta());
if (mPhi_EndCap_Hi) mPhi_EndCap_Hi->Fill(jet->phi());
//if (mConstituents_EndCap_Hi) mConstituents_EndCap_Hi->Fill(jet->nConstituents());
//if (mHFrac_EndCap_Hi) mHFrac_EndCap_Hi->Fill(jet->energyFractionHadronic());
}
if (fabs(jet->eta()) > 3.0) {
if (mPt_Forward_Hi && jet->pt()>100.) mPt_Forward_Hi->Fill(jet->pt());
if (mEta_Hi && jet->pt()>100.) mEta_Hi->Fill(jet->eta());
if (mPhi_Forward_Hi) mPhi_Forward_Hi->Fill(jet->phi());
//if (mConstituents_Forward_Hi) mConstituents_Forward_Hi->Fill(jet->nConstituents());
//if (mHFrac_Forward_Hi) mHFrac_Forward_Hi->Fill(jet->energyFractionHadronic());
}
if (mEta_Hi && jet->pt()>100.) mEta_Hi->Fill (jet->eta());
if (mPhi_Hi) mPhi_Hi->Fill (jet->phi());
if (mPt_Hi) mPt_Hi->Fill (jet->pt());
}
if (mPt) mPt->Fill (jet->pt());
if (mPt_1) mPt_1->Fill (jet->pt());
if (mPt_2) mPt_2->Fill (jet->pt());
if (mPt_3) mPt_3->Fill (jet->pt());
if (mEta) mEta->Fill (jet->eta());
if (mPhi) mPhi->Fill (jet->phi());
if (mPhiVSEta) mPhiVSEta->Fill(jet->eta(),jet->phi());
if (mConstituents) mConstituents->Fill (jet->nConstituents());
if (mHFrac) mHFrac->Fill (jet->energyFractionHadronic());
if (mEFrac) mEFrac->Fill (jet->emEnergyFraction());
// Fill NPV profiles
//--------------------------------------------------------------------
if (mPt_profile) mPt_profile ->Fill(numPV, jet->pt());
if (mEta_profile) mEta_profile ->Fill(numPV, jet->eta());
if (mPhi_profile) mPhi_profile ->Fill(numPV, jet->phi());
if (mConstituents_profile) mConstituents_profile->Fill(numPV, jet->nConstituents());
if (mHFrac_profile) mHFrac_profile ->Fill(numPV, jet->energyFractionHadronic());
if (mEFrac_profile) mEFrac_profile ->Fill(numPV, jet->emEnergyFraction());
if (fabs(jet->eta()) <= 1.3) {
if (mPt_Barrel) mPt_Barrel->Fill (jet->pt());
if (mPhi_Barrel) mPhi_Barrel->Fill (jet->phi());
//if (mE_Barrel) mE_Barrel->Fill (jet->energy());
if (mConstituents_Barrel) mConstituents_Barrel->Fill(jet->nConstituents());
if (mHFrac_Barrel) mHFrac_Barrel->Fill(jet->energyFractionHadronic());
if (mEFrac_Barrel) mEFrac_Barrel->Fill(jet->emEnergyFraction());
}
if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) {
if (mPt_EndCap) mPt_EndCap->Fill (jet->pt());
if (mPhi_EndCap) mPhi_EndCap->Fill (jet->phi());
//if (mE_EndCap) mE_EndCap->Fill (jet->energy());
if (mConstituents_EndCap) mConstituents_EndCap->Fill(jet->nConstituents());
if (mHFrac_EndCap) mHFrac_EndCap->Fill(jet->energyFractionHadronic());
if (mEFrac_EndCap) mEFrac_EndCap->Fill(jet->emEnergyFraction());
}
if (fabs(jet->eta()) > 3.0) {
if (mPt_Forward) mPt_Forward->Fill (jet->pt());
if (mPhi_Forward) mPhi_Forward->Fill (jet->phi());
//if (mE_Forward) mE_Forward->Fill (jet->energy());
if (mConstituents_Forward) mConstituents_Forward->Fill(jet->nConstituents());
if (mHFrac_Forward) mHFrac_Forward->Fill(jet->energyFractionHadronic());
if (mEFrac_Forward) mEFrac_Forward->Fill(jet->emEnergyFraction());
}
//if (mE) mE->Fill (jet->energy());
//if (mP) mP->Fill (jet->p());
// if (mMass) mMass->Fill (jet->mass());
if (mMaxEInEmTowers) mMaxEInEmTowers->Fill (jet->maxEInEmTowers());
if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (jet->maxEInHadTowers());
if (mHadEnergyInHO) mHadEnergyInHO->Fill (jet->hadEnergyInHO());
if (mHadEnergyInHB) mHadEnergyInHB->Fill (jet->hadEnergyInHB());
if (mHadEnergyInHF) mHadEnergyInHF->Fill (jet->hadEnergyInHF());
if (mHadEnergyInHE) mHadEnergyInHE->Fill (jet->hadEnergyInHE());
if (mEmEnergyInEB) mEmEnergyInEB->Fill (jet->emEnergyInEB());
if (mEmEnergyInEE) mEmEnergyInEE->Fill (jet->emEnergyInEE());
if (mEmEnergyInHF) mEmEnergyInHF->Fill (jet->emEnergyInHF());
if (mN90Hits) mN90Hits->Fill (jetID->n90Hits());
if (mfHPD) mfHPD->Fill (jetID->fHPD());
if (mresEMF) mresEMF->Fill (jetID->restrictedEMF());
if (mfRBX) mfRBX->Fill (jetID->fRBX());
//calculate correctly the dphi
if(numofjets>1) {
diff = fabs(fstPhi - sndPhi);
corr = 2*acos(-1.) - diff;
if(diff < acos(-1.)) {
dphi = diff;
} else {
dphi = corr;
}
}
}
}//pt cut
}
if (mNJets) mNJets->Fill (numofjets);
if (mDPhi && dphi>-998.) mDPhi->Fill (dphi);
if (mNJets_profile) mNJets_profile->Fill(numPV, numofjets);
}//not dijet
}
| void JetAnalyzer::beginJob | ( | DQMStore * | dbe | ) | [virtual] |
Inizialize parameters for histo binning.
Implements JetAnalyzerBase.
Definition at line 50 of file JetAnalyzer.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), jptDQMConfig_cff::eMax, jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, bTagSequences_cff::jetID, LogTrace, nbinsPV, Parameters::parameters, jptDQMConfig_cff::phiMax, jptDQMConfig_cff::phiMin, jptDQMConfig_cff::pMax, jptDQMConfig_cff::ptMax, PtMinSelector_cfg::ptMin, PVlow, PVup, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().
{
jetname = "jetAnalyzer";
LogTrace(jetname)<<"[JetAnalyzer] Parameters initialization";
dbe->setCurrentFolder("JetMET/Jet/"+_source);
jetME = dbe->book1D("jetReco", "jetReco", 3, 1, 4);
jetME->setBinLabel(1,"CaloJets",1);
//
jetID = new reco::helper::JetIDHelper(parameters.getParameter<ParameterSet>("JetIDParams"));
//
fillJIDPassFrac = parameters.getParameter<int>("fillJIDPassFrac");
makedijetselection = parameters.getParameter<int>("makedijetselection");
// monitoring of eta parameter
etaBin = parameters.getParameter<int>("etaBin");
etaMin = parameters.getParameter<double>("etaMin");
etaMax = parameters.getParameter<double>("etaMax");
// monitoring of phi paramater
phiBin = parameters.getParameter<int>("phiBin");
phiMin = parameters.getParameter<double>("phiMin");
phiMax = parameters.getParameter<double>("phiMax");
// monitoring of the transverse momentum
ptBin = parameters.getParameter<int>("ptBin");
ptMin = parameters.getParameter<double>("ptMin");
ptMax = parameters.getParameter<double>("ptMax");
//
eBin = parameters.getParameter<int>("eBin");
eMin = parameters.getParameter<double>("eMin");
eMax = parameters.getParameter<double>("eMax");
//
pBin = parameters.getParameter<int>("pBin");
pMin = parameters.getParameter<double>("pMin");
pMax = parameters.getParameter<double>("pMax");
//
_ptThreshold = parameters.getParameter<double>("ptThreshold");
_asymmetryThirdJetCut = parameters.getParameter<double>("asymmetryThirdJetCut");
_balanceThirdJetCut = parameters.getParameter<double>("balanceThirdJetCut");
_n90HitsMin = parameters.getParameter<int>("n90HitsMin");
_fHPDMax = parameters.getParameter<double>("fHPDMax");
_resEMFMin = parameters.getParameter<double>("resEMFMin");
_sigmaEtaMinTight = parameters.getParameter<double>("sigmaEtaMinTight");
_sigmaPhiMinTight = parameters.getParameter<double>("sigmaPhiMinTight");
_n90HitsMinLoose = parameters.getParameter<int>("n90HitsMinLoose");
_fHPDMaxLoose = parameters.getParameter<double>("fHPDMaxLoose");
_resEMFMinLoose = parameters.getParameter<double>("resEMFMinLoose");
_n90HitsMinTight = parameters.getParameter<int>("n90HitsMinTight");
_fHPDMaxTight = parameters.getParameter<double>("fHPDMaxTight");
_resEMFMinTight = parameters.getParameter<double>("resEMFMinTight");
// Generic jet parameters
mPt = dbe->book1D("Pt", "pt", ptBin, ptMin, ptMax);
mEta = dbe->book1D("Eta", "eta", etaBin, etaMin, etaMax);
mPhi = dbe->book1D("Phi", "phi", phiBin, phiMin, phiMax);
mConstituents = dbe->book1D("Constituents", "# of constituents", 50, 0, 100);
mHFrac = dbe->book1D("HFrac", "HFrac", 120, -0.1, 1.1);
mEFrac = dbe->book1D("EFrac", "EFrac", 120, -0.1, 1.1);
// Book NPV profiles
//----------------------------------------------------------------------------
mPt_profile = dbe->bookProfile("Pt_profile", "pt", nbinsPV, PVlow, PVup, ptBin, ptMin, ptMax);
mEta_profile = dbe->bookProfile("Eta_profile", "eta", nbinsPV, PVlow, PVup, etaBin, etaMin, etaMax);
mPhi_profile = dbe->bookProfile("Phi_profile", "phi", nbinsPV, PVlow, PVup, phiBin, phiMin, phiMax);
mConstituents_profile = dbe->bookProfile("Constituents_profile", "# of constituents", nbinsPV, PVlow, PVup, 50, 0, 100);
mHFrac_profile = dbe->bookProfile("HFrac_profile", "HFrac", nbinsPV, PVlow, PVup, 120, -0.1, 1.1);
mEFrac_profile = dbe->bookProfile("EFrac_profile", "EFrac", nbinsPV, PVlow, PVup, 120, -0.1, 1.1);
if (makedijetselection != 1)
mNJets_profile = dbe->bookProfile("NJets_profile", "number of jets", nbinsPV, PVlow, PVup, 100, 0, 100);
// Set NPV profiles x-axis title
//----------------------------------------------------------------------------
mPt_profile ->setAxisTitle("nvtx",1);
mEta_profile ->setAxisTitle("nvtx",1);
mPhi_profile ->setAxisTitle("nvtx",1);
mConstituents_profile->setAxisTitle("nvtx",1);
mHFrac_profile ->setAxisTitle("nvtx",1);
mEFrac_profile ->setAxisTitle("nvtx",1);
if (makedijetselection != 1) {
mNJets_profile->setAxisTitle("nvtx",1);
}
//mE = dbe->book1D("E", "E", eBin, eMin, eMax);
//mP = dbe->book1D("P", "P", pBin, pMin, pMax);
// mMass = dbe->book1D("Mass", "Mass", 100, 0, 25);
//
mPhiVSEta = dbe->book2D("PhiVSEta", "PhiVSEta", 50, etaMin, etaMax, 24, phiMin, phiMax);
if(makedijetselection!=1){
mPt_1 = dbe->book1D("Pt1", "Pt1", 20, 0, 100);
mPt_2 = dbe->book1D("Pt2", "Pt2", 60, 0, 300);
mPt_3 = dbe->book1D("Pt3", "Pt3", 100, 0, 5000);
// Low and high pt trigger paths
mPt_Lo = dbe->book1D("Pt_Lo", "Pt (Pass Low Pt Jet Trigger)", 20, 0, 100);
//mEta_Lo = dbe->book1D("Eta_Lo", "Eta (Pass Low Pt Jet Trigger)", etaBin, etaMin, etaMax);
mPhi_Lo = dbe->book1D("Phi_Lo", "Phi (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
mPt_Hi = dbe->book1D("Pt_Hi", "Pt (Pass Hi Pt Jet Trigger)", 60, 0, 300);
mEta_Hi = dbe->book1D("Eta_Hi", "Eta (Pass Hi Pt Jet Trigger)", etaBin, etaMin, etaMax);
mPhi_Hi = dbe->book1D("Phi_Hi", "Phi (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
mNJets = dbe->book1D("NJets", "number of jets", 100, 0, 100);
//mPt_Barrel_Lo = dbe->book1D("Pt_Barrel_Lo", "Pt Barrel (Pass Low Pt Jet Trigger)", 20, 0, 100);
//mPhi_Barrel_Lo = dbe->book1D("Phi_Barrel_Lo", "Phi Barrel (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
mConstituents_Barrel = dbe->book1D("Constituents_Barrel", "Constituents Barrel above", 50, 0, 100);
mHFrac_Barrel = dbe->book1D("HFrac_Barrel", "HFrac Barrel", 100, 0, 1);
mEFrac_Barrel = dbe->book1D("EFrac_Barrel", "EFrac Barrel", 110, -0.05, 1.05);
//mPt_EndCap_Lo = dbe->book1D("Pt_EndCap_Lo", "Pt EndCap (Pass Low Pt Jet Trigger)", 20, 0, 100);
//mPhi_EndCap_Lo = dbe->book1D("Phi_EndCap_Lo", "Phi EndCap (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
mConstituents_EndCap = dbe->book1D("Constituents_EndCap", "Constituents EndCap", 50, 0, 100);
mHFrac_EndCap = dbe->book1D("HFrac_Endcap", "HFrac EndCap", 100, 0, 1);
mEFrac_EndCap = dbe->book1D("EFrac_Endcap", "EFrac EndCap", 110, -0.05, 1.05);
//mPt_Forward_Lo = dbe->book1D("Pt_Forward_Lo", "Pt Forward (Pass Low Pt Jet Trigger)", 20, 0, 100);
//mPhi_Forward_Lo = dbe->book1D("Phi_Forward_Lo", "Phi Forward (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
mConstituents_Forward = dbe->book1D("Constituents_Forward", "Constituents Forward", 50, 0, 100);
mHFrac_Forward = dbe->book1D("HFrac_Forward", "HFrac Forward", 100, 0, 1);
mEFrac_Forward = dbe->book1D("EFrac_Forward", "EFrac Forward", 110, -0.05, 1.05);
mPt_Barrel_Hi = dbe->book1D("Pt_Barrel_Hi", "Pt Barrel (Pass Hi Pt Jet Trigger)", 60, 0, 300);
mPhi_Barrel_Hi = dbe->book1D("Phi_Barrel_Hi", "Phi Barrel (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
//mConstituents_Barrel_Hi = dbe->book1D("Constituents_Barrel_Hi", "Constituents Barrel (Pass Hi Pt Jet Trigger)", 50, 0, 100);
//mHFrac_Barrel_Hi = dbe->book1D("HFrac_Barrel_Hi", "HFrac Barrel (Pass Hi Pt Jet Trigger)", 100, 0, 1);
mPt_EndCap_Hi = dbe->book1D("Pt_EndCap_Hi", "Pt EndCap (Pass Hi Pt Jet Trigger)", 60, 0, 300);
mPhi_EndCap_Hi = dbe->book1D("Phi_EndCap_Hi", "Phi EndCap (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
//mConstituents_EndCap_Hi = dbe->book1D("Constituents_EndCap_Hi", "Constituents EndCap (Pass Hi Pt Jet Trigger)", 50, 0, 100);
//mHFrac_EndCap_Hi = dbe->book1D("HFrac_EndCap_Hi", "HFrac EndCap (Pass Hi Pt Jet Trigger)", 100, 0, 1);
mPt_Forward_Hi = dbe->book1D("Pt_Forward_Hi", "Pt Forward (Pass Hi Pt Jet Trigger)", 60, 0, 300);
mPhi_Forward_Hi = dbe->book1D("Phi_Forward_Hi", "Phi Forward (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
//mConstituents_Forward_Hi = dbe->book1D("Constituents_Forward_Hi", "Constituents Forward (Pass Hi Pt Jet Trigger)", 50, 0, 100);
//mHFrac_Forward_Hi = dbe->book1D("HFrac_Forward_Hi", "HFrac Forward (Pass Hi Pt Jet Trigger)", 100, 0, 1);
mPhi_Barrel = dbe->book1D("Phi_Barrel", "Phi_Barrel", phiBin, phiMin, phiMax);
//mE_Barrel = dbe->book1D("E_Barrel", "E_Barrel", eBin, eMin, eMax);
mPt_Barrel = dbe->book1D("Pt_Barrel", "Pt_Barrel", ptBin, ptMin, ptMax);
mPhi_EndCap = dbe->book1D("Phi_EndCap", "Phi_EndCap", phiBin, phiMin, phiMax);
//mE_EndCap = dbe->book1D("E_EndCap", "E_EndCap", eBin, eMin, 2*eMax);
mPt_EndCap = dbe->book1D("Pt_EndCap", "Pt_EndCap", ptBin, ptMin, ptMax);
mPhi_Forward = dbe->book1D("Phi_Forward", "Phi_Forward", phiBin, phiMin, phiMax);
//mE_Forward = dbe->book1D("E_Forward", "E_Forward", eBin, eMin, 4*eMax);
mPt_Forward = dbe->book1D("Pt_Forward", "Pt_Forward", ptBin, ptMin, ptMax);
// Leading Jet Parameters
mEtaFirst = dbe->book1D("EtaFirst", "EtaFirst", 100, -5, 5);
mPhiFirst = dbe->book1D("PhiFirst", "PhiFirst", 70, -3.5, 3.5);
//mEFirst = dbe->book1D("EFirst", "EFirst", 100, 0, 1000);
mPtFirst = dbe->book1D("PtFirst", "PtFirst", 100, 0, 500);
if(fillJIDPassFrac==1){//fillJIDPassFrac defines a collection of cleaned jets, for which we will want to fill the cleaning passing fraction
mLooseJIDPassFractionVSeta = dbe->bookProfile("LooseJIDPassFractionVSeta","LooseJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2);
mLooseJIDPassFractionVSpt = dbe->bookProfile("LooseJIDPassFractionVSpt","LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
mTightJIDPassFractionVSeta = dbe->bookProfile("TightJIDPassFractionVSeta","TightJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2);
mTightJIDPassFractionVSpt = dbe->bookProfile("TightJIDPassFractionVSpt","TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
}
}
// CaloJet specific
mMaxEInEmTowers = dbe->book1D("MaxEInEmTowers", "MaxEInEmTowers", 100, 0, 100);
mMaxEInHadTowers = dbe->book1D("MaxEInHadTowers", "MaxEInHadTowers", 100, 0, 100);
if(makedijetselection!=1) {
mHadEnergyInHO = dbe->book1D("HadEnergyInHO", "HadEnergyInHO", 100, 0, 10);
mHadEnergyInHB = dbe->book1D("HadEnergyInHB", "HadEnergyInHB", 100, 0, 50);
mHadEnergyInHF = dbe->book1D("HadEnergyInHF", "HadEnergyInHF", 100, 0, 50);
mHadEnergyInHE = dbe->book1D("HadEnergyInHE", "HadEnergyInHE", 100, 0, 100);
mEmEnergyInEB = dbe->book1D("EmEnergyInEB", "EmEnergyInEB", 100, 0, 50);
mEmEnergyInEE = dbe->book1D("EmEnergyInEE", "EmEnergyInEE", 100, 0, 50);
mEmEnergyInHF = dbe->book1D("EmEnergyInHF", "EmEnergyInHF", 120, -20, 100);
}
mDPhi = dbe->book1D("DPhi", "dPhi btw the two leading jets", 100, 0., acos(-1.));
//JetID variables
mresEMF = dbe->book1D("resEMF", "resEMF", 50, 0., 1.);
mN90Hits = dbe->book1D("N90Hits", "N90Hits", 50, 0., 50);
mfHPD = dbe->book1D("fHPD", "fHPD", 50, 0., 1.);
mfRBX = dbe->book1D("fRBX", "fRBX", 50, 0., 1.);
// msigmaEta = dbe->book1D("sigmaEta", "sigmaEta", 50, 0., 1.);
// msigmaPhi = dbe->book1D("sigmaPhi", "sigmaPhi", 50, 0., 0.5);
if(makedijetselection==1) {
mDijetAsymmetry = dbe->book1D("DijetAsymmetry", "DijetAsymmetry", 100, -1., 1.);
mDijetBalance = dbe->book1D("DijetBalance", "DijetBalance", 100, -2., 2.);
if (fillJIDPassFrac==1) {
mLooseJIDPassFractionVSeta = dbe->bookProfile("LooseJIDPassFractionVSeta","LooseJIDPassFractionVSeta",50, -3., 3.,0.,1.2);
mLooseJIDPassFractionVSpt = dbe->bookProfile("LooseJIDPassFractionVSpt","LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
mTightJIDPassFractionVSeta = dbe->bookProfile("TightJIDPassFractionVSeta","TightJIDPassFractionVSeta",50, -3., 3.,0.,1.2);
mTightJIDPassFractionVSpt = dbe->bookProfile("TightJIDPassFractionVSpt","TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
}
}
}
| void JetAnalyzer::endJob | ( | void | ) |
Finish up a job.
Definition at line 267 of file JetAnalyzer.cc.
References bTagSequences_cff::jetID.
{
delete jetID;
}
| int JetAnalyzer::getLeadJetFlag | ( | ) | [inline] |
| void JetAnalyzer::setJetHiPass | ( | int | pass | ) | [inline] |
| void JetAnalyzer::setJetLoPass | ( | int | pass | ) | [inline] |
| void JetAnalyzer::setLeadJetFlag | ( | int | flag | ) | [inline] |
| void JetAnalyzer::setSource | ( | std::string | source | ) | [inline] |
Definition at line 62 of file JetAnalyzer.h.
References _source, and LaserTracksInput_cfi::source.
double JetAnalyzer::_asymmetryThirdJetCut [private] |
Definition at line 96 of file JetAnalyzer.h.
double JetAnalyzer::_balanceThirdJetCut [private] |
Definition at line 97 of file JetAnalyzer.h.
double JetAnalyzer::_fHPDMax [private] |
Definition at line 102 of file JetAnalyzer.h.
double JetAnalyzer::_fHPDMaxLoose [private] |
Definition at line 107 of file JetAnalyzer.h.
double JetAnalyzer::_fHPDMaxTight [private] |
Definition at line 110 of file JetAnalyzer.h.
int JetAnalyzer::_JetHiPass [private] |
Definition at line 91 of file JetAnalyzer.h.
Referenced by setJetHiPass().
int JetAnalyzer::_JetLoPass [private] |
Definition at line 90 of file JetAnalyzer.h.
Referenced by setJetLoPass().
int JetAnalyzer::_leadJetFlag [private] |
Definition at line 92 of file JetAnalyzer.h.
Referenced by getLeadJetFlag(), and setLeadJetFlag().
int JetAnalyzer::_n90HitsMin [private] |
Definition at line 104 of file JetAnalyzer.h.
int JetAnalyzer::_n90HitsMinLoose [private] |
Definition at line 109 of file JetAnalyzer.h.
int JetAnalyzer::_n90HitsMinTight [private] |
Definition at line 112 of file JetAnalyzer.h.
double JetAnalyzer::_ptThreshold [private] |
Definition at line 94 of file JetAnalyzer.h.
double JetAnalyzer::_resEMFMin [private] |
Definition at line 103 of file JetAnalyzer.h.
double JetAnalyzer::_resEMFMinLoose [private] |
Definition at line 108 of file JetAnalyzer.h.
double JetAnalyzer::_resEMFMinTight [private] |
Definition at line 111 of file JetAnalyzer.h.
double JetAnalyzer::_sigmaEtaMinTight [private] |
Definition at line 113 of file JetAnalyzer.h.
double JetAnalyzer::_sigmaPhiMinTight [private] |
Definition at line 114 of file JetAnalyzer.h.
std::string JetAnalyzer::_source [private] |
Definition at line 86 of file JetAnalyzer.h.
Referenced by setSource().
int JetAnalyzer::_theend [private] |
Definition at line 93 of file JetAnalyzer.h.
int JetAnalyzer::eBin [private] |
Definition at line 129 of file JetAnalyzer.h.
double JetAnalyzer::eMax [private] |
Definition at line 131 of file JetAnalyzer.h.
double JetAnalyzer::eMin [private] |
Definition at line 130 of file JetAnalyzer.h.
int JetAnalyzer::etaBin [private] |
Definition at line 117 of file JetAnalyzer.h.
double JetAnalyzer::etaMax [private] |
Definition at line 119 of file JetAnalyzer.h.
double JetAnalyzer::etaMin [private] |
Definition at line 118 of file JetAnalyzer.h.
int JetAnalyzer::fillJIDPassFrac [private] |
Definition at line 106 of file JetAnalyzer.h.
reco::helper::JetIDHelper* JetAnalyzer::jetID [private] |
Definition at line 141 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::jetME [private] |
Definition at line 138 of file JetAnalyzer.h.
std::string JetAnalyzer::jetname [private] |
Definition at line 85 of file JetAnalyzer.h.
int JetAnalyzer::makedijetselection [private] |
Definition at line 99 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mConstituents [private] |
Definition at line 163 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mConstituents_Barrel [private] |
Definition at line 179 of file JetAnalyzer.h.
Definition at line 195 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mConstituents_EndCap [private] |
Definition at line 184 of file JetAnalyzer.h.
Definition at line 199 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mConstituents_Forward [private] |
Definition at line 189 of file JetAnalyzer.h.
Definition at line 203 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mConstituents_profile [private] |
Definition at line 268 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mDijetAsymmetry [private] |
Definition at line 259 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mDijetBalance [private] |
Definition at line 258 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mDPhi [private] |
Definition at line 216 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEFrac [private] |
Definition at line 165 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEFrac_Barrel [private] |
Definition at line 181 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEFrac_EndCap [private] |
Definition at line 186 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEFrac_Forward [private] |
Definition at line 191 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEFrac_profile [private] |
Definition at line 270 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEmEnergyInEB [private] |
Definition at line 232 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEmEnergyInEE [private] |
Definition at line 233 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEmEnergyInHF [private] |
Definition at line 234 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEta [private] |
Definition at line 161 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEta_Hi [private] |
Definition at line 253 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEta_profile [private] |
Definition at line 266 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEtaFirst [private] |
Definition at line 219 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mfHPD [private] |
Definition at line 238 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mfRBX [private] |
Definition at line 239 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHB [private] |
Definition at line 229 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHE [private] |
Definition at line 231 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHF [private] |
Definition at line 230 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHO [private] |
Definition at line 228 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac [private] |
Definition at line 164 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Barrel [private] |
Definition at line 180 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Barrel_Hi [private] |
Definition at line 196 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_EndCap [private] |
Definition at line 185 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_EndCap_Hi [private] |
Definition at line 200 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Forward [private] |
Definition at line 190 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Forward_Hi [private] |
Definition at line 204 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_profile [private] |
Definition at line 269 of file JetAnalyzer.h.
Definition at line 243 of file JetAnalyzer.h.
Definition at line 244 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mMaxEInEmTowers [private] |
Definition at line 226 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mMaxEInHadTowers [private] |
Definition at line 227 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mN90Hits [private] |
Definition at line 237 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mNJets [private] |
Definition at line 215 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mNJets_profile [private] |
Definition at line 264 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi [private] |
Definition at line 162 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Barrel [private] |
Definition at line 169 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Barrel_Hi [private] |
Definition at line 194 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_EndCap [private] |
Definition at line 172 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_EndCap_Hi [private] |
Definition at line 198 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Forward [private] |
Definition at line 175 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Forward_Hi [private] |
Definition at line 202 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Hi [private] |
Definition at line 254 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Lo [private] |
Definition at line 250 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_profile [private] |
Definition at line 267 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhiFirst [private] |
Definition at line 220 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhiVSEta [private] |
Definition at line 166 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt [private] |
Definition at line 157 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_1 [private] |
Definition at line 158 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_2 [private] |
Definition at line 159 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_3 [private] |
Definition at line 160 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Barrel [private] |
Definition at line 168 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Barrel_Hi [private] |
Definition at line 193 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_EndCap [private] |
Definition at line 171 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_EndCap_Hi [private] |
Definition at line 197 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Forward [private] |
Definition at line 174 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Forward_Hi [private] |
Definition at line 201 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Hi [private] |
Definition at line 255 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Lo [private] |
Definition at line 251 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_profile [private] |
Definition at line 265 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPtFirst [private] |
Definition at line 222 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mresEMF [private] |
Definition at line 240 of file JetAnalyzer.h.
Definition at line 245 of file JetAnalyzer.h.
Definition at line 246 of file JetAnalyzer.h.
edm::ParameterSet JetAnalyzer::parameters [private] |
Definition at line 83 of file JetAnalyzer.h.
int JetAnalyzer::pBin [private] |
Definition at line 133 of file JetAnalyzer.h.
int JetAnalyzer::phiBin [private] |
Definition at line 121 of file JetAnalyzer.h.
double JetAnalyzer::phiMax [private] |
Definition at line 123 of file JetAnalyzer.h.
double JetAnalyzer::phiMin [private] |
Definition at line 122 of file JetAnalyzer.h.
double JetAnalyzer::pMax [private] |
Definition at line 135 of file JetAnalyzer.h.
double JetAnalyzer::pMin [private] |
Definition at line 134 of file JetAnalyzer.h.
int JetAnalyzer::ptBin [private] |
Definition at line 125 of file JetAnalyzer.h.
double JetAnalyzer::ptMax [private] |
Definition at line 127 of file JetAnalyzer.h.
double JetAnalyzer::ptMin [private] |
Definition at line 126 of file JetAnalyzer.h.
Definition at line 88 of file JetAnalyzer.h.
1.7.3