29 _asymmetryThirdJetCut = 5.;
30 _balanceThirdJetCut = 0.2;
40 _sigmaEtaMinTight=-999.;
41 _sigmaPhiMinTight=-999.;
52 jetname =
"jetAnalyzer";
54 LogTrace(jetname)<<
"[JetAnalyzer] Parameters initialization";
57 jetME = dbe->
book1D(
"jetReco",
"jetReco", 3, 1, 4);
64 fillJIDPassFrac =
parameters.getParameter<
int>(
"fillJIDPassFrac");
65 makedijetselection =
parameters.getParameter<
int>(
"makedijetselection");
68 etaBin =
parameters.getParameter<
int>(
"etaBin");
73 phiBin =
parameters.getParameter<
int>(
"phiBin");
84 eMin =
parameters.getParameter<
double>(
"eMin");
89 pMin =
parameters.getParameter<
double>(
"pMin");
93 _ptThreshold =
parameters.getParameter<
double>(
"ptThreshold");
94 _asymmetryThirdJetCut =
parameters.getParameter<
double>(
"asymmetryThirdJetCut");
95 _balanceThirdJetCut =
parameters.getParameter<
double>(
"balanceThirdJetCut");
96 _n90HitsMin =
parameters.getParameter<
int>(
"n90HitsMin");
97 _fHPDMax =
parameters.getParameter<
double>(
"fHPDMax");
98 _resEMFMin =
parameters.getParameter<
double>(
"resEMFMin");
99 _sigmaEtaMinTight =
parameters.getParameter<
double>(
"sigmaEtaMinTight");
100 _sigmaPhiMinTight =
parameters.getParameter<
double>(
"sigmaPhiMinTight");
102 _n90HitsMinLoose =
parameters.getParameter<
int>(
"n90HitsMinLoose");
103 _fHPDMaxLoose =
parameters.getParameter<
double>(
"fHPDMaxLoose");
104 _resEMFMinLoose =
parameters.getParameter<
double>(
"resEMFMinLoose");
105 _n90HitsMinTight =
parameters.getParameter<
int>(
"n90HitsMinTight");
106 _fHPDMaxTight =
parameters.getParameter<
double>(
"fHPDMaxTight");
107 _resEMFMinTight =
parameters.getParameter<
double>(
"resEMFMinTight");
111 mPt = dbe->
book1D(
"Pt",
"pt", ptBin, ptMin, ptMax);
112 mEta = dbe->
book1D(
"Eta",
"eta", etaBin, etaMin, etaMax);
113 mPhi = dbe->
book1D(
"Phi",
"phi", phiBin, phiMin, phiMax);
114 mConstituents = dbe->
book1D(
"Constituents",
"# of constituents", 50, 0, 100);
115 mHFrac = dbe->
book1D(
"HFrac",
"HFrac", 120, -0.1, 1.1);
116 mEFrac = dbe->
book1D(
"EFrac",
"EFrac", 120, -0.1, 1.1);
128 if (makedijetselection != 1)
135 mEta_profile ->setAxisTitle(
"nvtx",1);
136 mPhi_profile ->setAxisTitle(
"nvtx",1);
137 mConstituents_profile->setAxisTitle(
"nvtx",1);
138 mHFrac_profile ->setAxisTitle(
"nvtx",1);
139 mEFrac_profile ->setAxisTitle(
"nvtx",1);
141 if (makedijetselection != 1) {
142 mNJets_profile->setAxisTitle(
"nvtx",1);
150 mPhiVSEta = dbe->
book2D(
"PhiVSEta",
"PhiVSEta", 50, etaMin, etaMax, 24, phiMin, phiMax);
151 if(makedijetselection!=1){
152 mPt_1 = dbe->
book1D(
"Pt1",
"Pt1", 20, 0, 100);
153 mPt_2 = dbe->
book1D(
"Pt2",
"Pt2", 60, 0, 300);
154 mPt_3 = dbe->
book1D(
"Pt3",
"Pt3", 100, 0, 5000);
156 mPt_Lo = dbe->
book1D(
"Pt_Lo",
"Pt (Pass Low Pt Jet Trigger)", 20, 0, 100);
158 mPhi_Lo = dbe->
book1D(
"Phi_Lo",
"Phi (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
160 mPt_Hi = dbe->
book1D(
"Pt_Hi",
"Pt (Pass Hi Pt Jet Trigger)", 60, 0, 300);
161 mEta_Hi = dbe->
book1D(
"Eta_Hi",
"Eta (Pass Hi Pt Jet Trigger)", etaBin, etaMin, etaMax);
162 mPhi_Hi = dbe->
book1D(
"Phi_Hi",
"Phi (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
163 mNJets = dbe->
book1D(
"NJets",
"number of jets", 100, 0, 100);
167 mConstituents_Barrel = dbe->
book1D(
"Constituents_Barrel",
"Constituents Barrel above", 50, 0, 100);
168 mHFrac_Barrel = dbe->
book1D(
"HFrac_Barrel",
"HFrac Barrel", 100, 0, 1);
169 mEFrac_Barrel = dbe->
book1D(
"EFrac_Barrel",
"EFrac Barrel", 110, -0.05, 1.05);
173 mConstituents_EndCap = dbe->
book1D(
"Constituents_EndCap",
"Constituents EndCap", 50, 0, 100);
174 mHFrac_EndCap = dbe->
book1D(
"HFrac_Endcap",
"HFrac EndCap", 100, 0, 1);
175 mEFrac_EndCap = dbe->
book1D(
"EFrac_Endcap",
"EFrac EndCap", 110, -0.05, 1.05);
179 mConstituents_Forward = dbe->
book1D(
"Constituents_Forward",
"Constituents Forward", 50, 0, 100);
180 mHFrac_Forward = dbe->
book1D(
"HFrac_Forward",
"HFrac Forward", 100, 0, 1);
181 mEFrac_Forward = dbe->
book1D(
"EFrac_Forward",
"EFrac Forward", 110, -0.05, 1.05);
183 mPt_Barrel_Hi = dbe->
book1D(
"Pt_Barrel_Hi",
"Pt Barrel (Pass Hi Pt Jet Trigger)", 60, 0, 300);
184 mPhi_Barrel_Hi = dbe->
book1D(
"Phi_Barrel_Hi",
"Phi Barrel (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
188 mPt_EndCap_Hi = dbe->
book1D(
"Pt_EndCap_Hi",
"Pt EndCap (Pass Hi Pt Jet Trigger)", 60, 0, 300);
189 mPhi_EndCap_Hi = dbe->
book1D(
"Phi_EndCap_Hi",
"Phi EndCap (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
193 mPt_Forward_Hi = dbe->
book1D(
"Pt_Forward_Hi",
"Pt Forward (Pass Hi Pt Jet Trigger)", 60, 0, 300);
194 mPhi_Forward_Hi = dbe->
book1D(
"Phi_Forward_Hi",
"Phi Forward (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
198 mPhi_Barrel = dbe->
book1D(
"Phi_Barrel",
"Phi_Barrel", phiBin, phiMin, phiMax);
200 mPt_Barrel = dbe->
book1D(
"Pt_Barrel",
"Pt_Barrel", ptBin, ptMin, ptMax);
202 mPhi_EndCap = dbe->
book1D(
"Phi_EndCap",
"Phi_EndCap", phiBin, phiMin, phiMax);
204 mPt_EndCap = dbe->
book1D(
"Pt_EndCap",
"Pt_EndCap", ptBin, ptMin, ptMax);
206 mPhi_Forward = dbe->
book1D(
"Phi_Forward",
"Phi_Forward", phiBin, phiMin, phiMax);
208 mPt_Forward = dbe->
book1D(
"Pt_Forward",
"Pt_Forward", ptBin, ptMin, ptMax);
211 mEtaFirst = dbe->
book1D(
"EtaFirst",
"EtaFirst", 100, -5, 5);
212 mPhiFirst = dbe->
book1D(
"PhiFirst",
"PhiFirst", 70, -3.5, 3.5);
214 mPtFirst = dbe->
book1D(
"PtFirst",
"PtFirst", 100, 0, 500);
215 if(fillJIDPassFrac==1){
216 mLooseJIDPassFractionVSeta = dbe->
bookProfile(
"LooseJIDPassFractionVSeta",
"LooseJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2);
217 mLooseJIDPassFractionVSpt = dbe->
bookProfile(
"LooseJIDPassFractionVSpt",
"LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
218 mTightJIDPassFractionVSeta = dbe->
bookProfile(
"TightJIDPassFractionVSeta",
"TightJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2);
219 mTightJIDPassFractionVSpt = dbe->
bookProfile(
"TightJIDPassFractionVSpt",
"TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
225 mMaxEInEmTowers = dbe->
book1D(
"MaxEInEmTowers",
"MaxEInEmTowers", 100, 0, 100);
226 mMaxEInHadTowers = dbe->
book1D(
"MaxEInHadTowers",
"MaxEInHadTowers", 100, 0, 100);
227 if(makedijetselection!=1) {
228 mHadEnergyInHO = dbe->
book1D(
"HadEnergyInHO",
"HadEnergyInHO", 100, 0, 10);
229 mHadEnergyInHB = dbe->
book1D(
"HadEnergyInHB",
"HadEnergyInHB", 100, 0, 50);
230 mHadEnergyInHF = dbe->
book1D(
"HadEnergyInHF",
"HadEnergyInHF", 100, 0, 50);
231 mHadEnergyInHE = dbe->
book1D(
"HadEnergyInHE",
"HadEnergyInHE", 100, 0, 100);
232 mEmEnergyInEB = dbe->
book1D(
"EmEnergyInEB",
"EmEnergyInEB", 100, 0, 50);
233 mEmEnergyInEE = dbe->
book1D(
"EmEnergyInEE",
"EmEnergyInEE", 100, 0, 50);
234 mEmEnergyInHF = dbe->
book1D(
"EmEnergyInHF",
"EmEnergyInHF", 120, -20, 100);
236 mDPhi = dbe->
book1D(
"DPhi",
"dPhi btw the two leading jets", 100, 0., acos(-1.));
240 mresEMF = dbe->
book1D(
"resEMF",
"resEMF", 50, 0., 1.);
241 mN90Hits = dbe->
book1D(
"N90Hits",
"N90Hits", 50, 0., 50);
242 mfHPD = dbe->
book1D(
"fHPD",
"fHPD", 50, 0., 1.);
243 mfRBX = dbe->
book1D(
"fRBX",
"fRBX", 50, 0., 1.);
248 if(makedijetselection==1) {
249 mDijetAsymmetry = dbe->
book1D(
"DijetAsymmetry",
"DijetAsymmetry", 100, -1., 1.);
250 mDijetBalance = dbe->
book1D(
"DijetBalance",
"DijetBalance", 100, -2., 2.);
251 if (fillJIDPassFrac==1) {
252 mLooseJIDPassFractionVSeta = dbe->
bookProfile(
"LooseJIDPassFractionVSeta",
"LooseJIDPassFractionVSeta",50, -3., 3.,0.,1.2);
253 mLooseJIDPassFractionVSpt = dbe->
bookProfile(
"LooseJIDPassFractionVSpt",
"LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
254 mTightJIDPassFractionVSeta = dbe->
bookProfile(
"TightJIDPassFractionVSeta",
"TightJIDPassFractionVSeta",50, -3., 3.,0.,1.2);
255 mTightJIDPassFractionVSpt = dbe->
bookProfile(
"TightJIDPassFractionVSpt",
"TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
281 double dphi = -999. ;
282 bool thiscleaned=
false;
283 bool Loosecleaned=
false;
284 bool Tightcleaned=
false;
285 bool thisemfclean=
true;
286 bool emfcleanLoose=
true;
287 bool emfcleanTight=
true;
289 srand( iEvent.
id().
event() % 10000);
292 if (makedijetselection == 1) {
295 if(caloJets.size()>=2){
296 double dphiDJ = -999. ;
297 bool emfcleanLooseFirstJet=
true;
298 bool emfcleanLooseSecondJet=
true;
299 bool emfcleanTightFirstJet=
true;
300 bool emfcleanTightSecondJet=
true;
301 bool LoosecleanedFirstJet =
false;
302 bool LoosecleanedSecondJet =
false;
303 bool TightcleanedFirstJet =
false;
304 bool TightcleanedSecondJet =
false;
306 if ((caloJets.at(0)).pt() > _ptThreshold && (caloJets.at(1)).pt() > _ptThreshold ) {
307 if(fabs((caloJets.at(0)).
eta())<3. && fabs((caloJets.at(1)).
eta())<3. ){
309 dphiDJ = fabs((caloJets.at(0)).
phi()-(caloJets.at(1)).
phi());
310 if (dphiDJ > 3.14) dphiDJ=fabs(dphiDJ -6.28 );
312 if (mDPhi) mDPhi->Fill (dphiDJ);
314 if(fabs(dphiDJ)>2.1){
316 emfcleanLooseFirstJet=
true;
317 emfcleanTightFirstJet=
true;
318 emfcleanLooseSecondJet=
true;
319 emfcleanTightSecondJet=
true;
321 jetID->calculate(iEvent, (caloJets.at(0)));
322 if(
jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(0)).
eta())<2.6) emfcleanLooseFirstJet=
false;
323 if(
jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(0)).
eta())<2.6) emfcleanTightFirstJet=
false;
324 if(
jetID->n90Hits()>=_n90HitsMinLoose &&
jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseFirstJet) LoosecleanedFirstJet=
true;
325 if(
jetID->n90Hits()>=_n90HitsMinTight &&
jetID->fHPD()<_fHPDMaxTight &&
sqrt((caloJets.at(0)).etaetaMoment())>_sigmaEtaMinTight &&
sqrt((caloJets.at(0)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightFirstJet) TightcleanedFirstJet=
true;
327 if (mN90Hits) mN90Hits->Fill (
jetID->n90Hits());
328 if (mfHPD) mfHPD->Fill (
jetID->fHPD());
329 if (mresEMF) mresEMF->Fill (
jetID->restrictedEMF());
330 if (mfRBX) mfRBX->Fill (
jetID->fRBX());
333 jetID->calculate(iEvent, (caloJets.at(1)));
334 if(
jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(1)).
eta())<2.6) emfcleanLooseSecondJet=
false;
335 if(
jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(1)).
eta())<2.6) emfcleanTightSecondJet=
false;
336 if(
jetID->n90Hits()>=_n90HitsMinLoose &&
jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseSecondJet) LoosecleanedSecondJet=
true;
337 if(
jetID->n90Hits()>=_n90HitsMinTight &&
jetID->fHPD()<_fHPDMaxTight &&
sqrt((caloJets.at(1)).etaetaMoment())>_sigmaEtaMinTight &&
sqrt((caloJets.at(1)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightSecondJet) TightcleanedSecondJet=
true;
339 if (mN90Hits) mN90Hits->Fill (
jetID->n90Hits());
340 if (mfHPD) mfHPD->Fill (
jetID->fHPD());
341 if (mresEMF) mresEMF->Fill (
jetID->restrictedEMF());
342 if (mfRBX) mfRBX->Fill (
jetID->fRBX());
344 if(LoosecleanedFirstJet && LoosecleanedSecondJet) {
346 if (mPt) mPt->Fill ((caloJets.at(0)).pt());
347 if (mEta) mEta->Fill ((caloJets.at(0)).
eta());
348 if (mPhi) mPhi->Fill ((caloJets.at(0)).
phi());
349 if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(0)).
eta(),(caloJets.at(0)).
phi());
350 if (mConstituents) mConstituents->Fill ((caloJets.at(0)).nConstituents());
351 if (mHFrac) mHFrac->Fill ((caloJets.at(0)).energyFractionHadronic());
352 if (mEFrac) mEFrac->Fill ((caloJets.at(0)).emEnergyFraction());
356 if (mMaxEInEmTowers) mMaxEInEmTowers->Fill ((caloJets.at(0)).maxEInEmTowers());
357 if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(0)).maxEInHadTowers());
358 if (mN90Hits) mN90Hits->Fill (
jetID->n90Hits());
359 if (mfHPD) mfHPD->Fill (
jetID->fHPD());
360 if (mresEMF) mresEMF->Fill (
jetID->restrictedEMF());
361 if (mfRBX) mfRBX->Fill (
jetID->fRBX());
367 if (mPt) mPt->Fill ((caloJets.at(1)).pt());
368 if (mEta) mEta->Fill ((caloJets.at(1)).
eta());
369 if (mPhi) mPhi->Fill ((caloJets.at(1)).
phi());
370 if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(1)).
eta(),(caloJets.at(1)).
phi());
371 if (mConstituents) mConstituents->Fill ((caloJets.at(1)).nConstituents());
372 if (mHFrac) mHFrac->Fill ((caloJets.at(1)).energyFractionHadronic());
373 if (mEFrac) mEFrac->Fill ((caloJets.at(1)).emEnergyFraction());
377 if (mMaxEInEmTowers) mMaxEInEmTowers->Fill ((caloJets.at(1)).maxEInEmTowers());
378 if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(1)).maxEInHadTowers());
387 for (
int ijet=0; ijet<2; ijet++) {
389 if (mPt_profile) mPt_profile ->Fill(numPV, (caloJets.at(ijet)).pt());
390 if (mEta_profile) mEta_profile ->Fill(numPV, (caloJets.at(ijet)).
eta());
391 if (mPhi_profile) mPhi_profile ->Fill(numPV, (caloJets.at(ijet)).
phi());
392 if (mConstituents_profile) mConstituents_profile->Fill(numPV, (caloJets.at(ijet)).nConstituents());
393 if (mHFrac_profile) mHFrac_profile ->Fill(numPV, (caloJets.at(ijet)).energyFractionHadronic());
394 if (mEFrac_profile) mEFrac_profile ->Fill(numPV, (caloJets.at(ijet)).emEnergyFraction());
400 if(fillJIDPassFrac==1) {
401 if(LoosecleanedFirstJet) {
402 mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).
eta(),1.);
403 mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.);
405 mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).
eta(),0.);
406 mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.);
408 if(LoosecleanedSecondJet) {
409 mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).
eta(),1.);
410 mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.);
412 mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).
eta(),0.);
413 mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.);
416 if(TightcleanedFirstJet) {
417 mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).
eta(),1.);
418 mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.);
420 mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).
eta(),0.);
421 mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.);
423 if(TightcleanedSecondJet) {
424 mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).
eta(),1.);
425 mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.);
427 mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).
eta(),0.);
428 mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.);
436 if (fabs(caloJets.at(0).eta() < 1.4)) {
437 double pt_dijet = (caloJets.at(0).pt() + caloJets.at(1).pt())/2;
439 double dPhi = fabs((caloJets.at(0)).
phi()-(caloJets.at(1)).
phi());
440 if (dPhi > 3.14) dPhi=fabs(dPhi -6.28 );
447 int randJet =
rand() % 2;
449 if (fabs(caloJets.at(1).eta() < 1.4)) {
466 bool thirdJetCut =
true;
467 for (
unsigned int third = 2; third < caloJets.size(); ++third)
468 if (caloJets.at(third).pt() > _asymmetryThirdJetCut)
471 double dijetAsymmetry = (caloJets.at(jet1).pt() - caloJets.at(jet2).pt()) / (caloJets.at(jet1).pt() + caloJets.at(jet2).pt());
472 mDijetAsymmetry->Fill(dijetAsymmetry);
481 pt_barrel = caloJets.at(jet1).pt();
482 pt_probe = caloJets.at(jet2).pt();
493 bool thirdJetCut =
true;
494 for (
unsigned int third = 2; third < caloJets.size(); ++third)
495 if (caloJets.at(third).pt()/pt_dijet > _balanceThirdJetCut)
498 double dijetBalance = (pt_probe - pt_barrel) / pt_dijet;
499 mDijetBalance->Fill(dijetBalance);
506 for (reco::CaloJetCollection::const_iterator
jet = caloJets.begin();
jet!=caloJets.end(); ++
jet) {
507 LogTrace(jetname)<<
"[JetAnalyzer] Analyze Calo Jet";
510 if (
jet == caloJets.begin()) {
516 if (
jet == (caloJets.begin()+1)) sndPhi =
jet->phi();
520 if (
jet->pt() > _ptThreshold) {
525 if(
jetID->restrictedEMF()<_resEMFMin && fabs(
jet->eta())<2.6) thisemfclean=
false;
526 if(
jetID->n90Hits()>=_n90HitsMin &&
jetID->fHPD()<_fHPDMax && thisemfclean) thiscleaned=
true;
528 if(
jetID->n90Hits()>=_n90HitsMinLoose &&
jetID->fHPD()<_fHPDMaxLoose && emfcleanLoose) Loosecleaned=
true;
529 if(
jetID->n90Hits()>=_n90HitsMinTight &&
jetID->fHPD()<_fHPDMaxTight &&
sqrt(
jet->etaetaMoment())>_sigmaEtaMinTight &&
sqrt(
jet->phiphiMoment())>_sigmaPhiMinTight && emfcleanTight) Tightcleaned=
true;
531 if(fillJIDPassFrac==1) {
533 mLooseJIDPassFractionVSeta->Fill(
jet->eta(),1.);
534 mLooseJIDPassFractionVSpt->Fill(
jet->pt(),1.);
536 mLooseJIDPassFractionVSeta->Fill(
jet->eta(),0.);
537 mLooseJIDPassFractionVSpt->Fill(
jet->pt(),0.);
541 mTightJIDPassFractionVSeta->Fill(
jet->eta(),1.);
542 mTightJIDPassFractionVSpt->Fill(
jet->pt(),1.);
544 mTightJIDPassFractionVSeta->Fill(
jet->eta(),0.);
545 mTightJIDPassFractionVSpt->Fill(
jet->pt(),0.);
555 if (_leadJetFlag == 1) {
556 if (mEtaFirst) mEtaFirst->Fill (
jet->eta());
557 if (mPhiFirst) mPhiFirst->Fill (
jet->phi());
559 if (mPtFirst) mPtFirst->Fill (
jet->pt());
562 if (_JetLoPass == 1) {
579 if (mPhi_Lo) mPhi_Lo->Fill (
jet->phi());
580 if (mPt_Lo) mPt_Lo->Fill (
jet->pt());
584 if (_JetHiPass == 1) {
585 if (fabs(
jet->eta()) <= 1.3) {
586 if (mPt_Barrel_Hi &&
jet->pt()>100.) mPt_Barrel_Hi->Fill(
jet->pt());
587 if (mEta_Hi &&
jet->pt()>100.) mEta_Hi->Fill(
jet->eta());
588 if (mPhi_Barrel_Hi) mPhi_Barrel_Hi->Fill(
jet->phi());
592 if ( (fabs(
jet->eta()) > 1.3) && (fabs(
jet->eta()) <= 3) ) {
593 if (mPt_EndCap_Hi &&
jet->pt()>100.) mPt_EndCap_Hi->Fill(
jet->pt());
594 if (mEta_Hi &&
jet->pt()>100.) mEta_Hi->Fill(
jet->eta());
595 if (mPhi_EndCap_Hi) mPhi_EndCap_Hi->Fill(
jet->phi());
599 if (fabs(
jet->eta()) > 3.0) {
600 if (mPt_Forward_Hi &&
jet->pt()>100.) mPt_Forward_Hi->Fill(
jet->pt());
601 if (mEta_Hi &&
jet->pt()>100.) mEta_Hi->Fill(
jet->eta());
602 if (mPhi_Forward_Hi) mPhi_Forward_Hi->Fill(
jet->phi());
607 if (mEta_Hi &&
jet->pt()>100.) mEta_Hi->Fill (
jet->eta());
608 if (mPhi_Hi) mPhi_Hi->Fill (
jet->phi());
609 if (mPt_Hi) mPt_Hi->Fill (
jet->pt());
612 if (mPt) mPt->Fill (
jet->pt());
613 if (mPt_1) mPt_1->Fill (
jet->pt());
614 if (mPt_2) mPt_2->Fill (
jet->pt());
615 if (mPt_3) mPt_3->Fill (
jet->pt());
616 if (mEta) mEta->Fill (
jet->eta());
617 if (mPhi) mPhi->Fill (
jet->phi());
619 if (mPhiVSEta) mPhiVSEta->Fill(
jet->eta(),
jet->phi());
621 if (mConstituents) mConstituents->Fill (
jet->nConstituents());
622 if (mHFrac) mHFrac->Fill (
jet->energyFractionHadronic());
623 if (mEFrac) mEFrac->Fill (
jet->emEnergyFraction());
628 if (mPt_profile) mPt_profile ->Fill(numPV,
jet->pt());
629 if (mEta_profile) mEta_profile ->Fill(numPV,
jet->eta());
630 if (mPhi_profile) mPhi_profile ->Fill(numPV,
jet->phi());
631 if (mConstituents_profile) mConstituents_profile->Fill(numPV,
jet->nConstituents());
632 if (mHFrac_profile) mHFrac_profile ->Fill(numPV,
jet->energyFractionHadronic());
633 if (mEFrac_profile) mEFrac_profile ->Fill(numPV,
jet->emEnergyFraction());
636 if (fabs(
jet->eta()) <= 1.3) {
637 if (mPt_Barrel) mPt_Barrel->Fill (
jet->pt());
638 if (mPhi_Barrel) mPhi_Barrel->Fill (
jet->phi());
640 if (mConstituents_Barrel) mConstituents_Barrel->Fill(
jet->nConstituents());
641 if (mHFrac_Barrel) mHFrac_Barrel->Fill(
jet->energyFractionHadronic());
642 if (mEFrac_Barrel) mEFrac_Barrel->Fill(
jet->emEnergyFraction());
644 if ( (fabs(
jet->eta()) > 1.3) && (fabs(
jet->eta()) <= 3) ) {
645 if (mPt_EndCap) mPt_EndCap->Fill (
jet->pt());
646 if (mPhi_EndCap) mPhi_EndCap->Fill (
jet->phi());
648 if (mConstituents_EndCap) mConstituents_EndCap->Fill(
jet->nConstituents());
649 if (mHFrac_EndCap) mHFrac_EndCap->Fill(
jet->energyFractionHadronic());
650 if (mEFrac_EndCap) mEFrac_EndCap->Fill(
jet->emEnergyFraction());
652 if (fabs(
jet->eta()) > 3.0) {
653 if (mPt_Forward) mPt_Forward->Fill (
jet->pt());
654 if (mPhi_Forward) mPhi_Forward->Fill (
jet->phi());
656 if (mConstituents_Forward) mConstituents_Forward->Fill(
jet->nConstituents());
657 if (mHFrac_Forward) mHFrac_Forward->Fill(
jet->energyFractionHadronic());
658 if (mEFrac_Forward) mEFrac_Forward->Fill(
jet->emEnergyFraction());
665 if (mMaxEInEmTowers) mMaxEInEmTowers->Fill (
jet->maxEInEmTowers());
666 if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (
jet->maxEInHadTowers());
668 if (mHadEnergyInHO) mHadEnergyInHO->Fill (
jet->hadEnergyInHO());
669 if (mHadEnergyInHB) mHadEnergyInHB->Fill (
jet->hadEnergyInHB());
670 if (mHadEnergyInHF) mHadEnergyInHF->Fill (
jet->hadEnergyInHF());
671 if (mHadEnergyInHE) mHadEnergyInHE->Fill (
jet->hadEnergyInHE());
672 if (mEmEnergyInEB) mEmEnergyInEB->Fill (
jet->emEnergyInEB());
673 if (mEmEnergyInEE) mEmEnergyInEE->Fill (
jet->emEnergyInEE());
674 if (mEmEnergyInHF) mEmEnergyInHF->Fill (
jet->emEnergyInHF());
676 if (mN90Hits) mN90Hits->Fill (
jetID->n90Hits());
677 if (mfHPD) mfHPD->Fill (
jetID->fHPD());
678 if (mresEMF) mresEMF->Fill (
jetID->restrictedEMF());
679 if (mfRBX) mfRBX->Fill (
jetID->fRBX());
683 diff = fabs(fstPhi - sndPhi);
684 corr = 2*acos(-1.) -
diff;
685 if(diff < acos(-1.)) {
695 if (mNJets) mNJets->Fill (numofjets);
696 if (mDPhi && dphi>-998.) mDPhi->Fill (dphi);
699 if (mNJets_profile) mNJets_profile->Fill(numPV, numofjets);
EventNumber_t event() const
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
JetAnalyzer(const edm::ParameterSet &)
Constructor.
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
virtual ~JetAnalyzer()
Destructor.
double dPhi(double phi1, double phi2)
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
void endJob()
Finish up a job.
void beginJob(DQMStore *dbe)
Inizialize parameters for histo binning.
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
void setCurrentFolder(const std::string &fullpath)
void analyze(const edm::Event &, const edm::EventSetup &, const reco::CaloJetCollection &caloJets, const int numPV)
Get the analysis.
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects