Get the analysis.
1577 Int_t JetLoPass = 0;
1578 Int_t JetHiPass = 0;
1580 if (triggerResults.
isValid()){
1583 const unsigned int nTrig(triggerNames.
size());
1584 for (
unsigned int i=0;
i<nTrig;++
i)
1598 std::cout <<
">>> Trigger Lo = " << JetLoPass
1599 <<
" Hi = " << JetHiPass
1608 if (!vertexHandle.
isValid()) {
1609 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1610 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1613 if ( vertexHandle.
isValid() ){
1615 numPV = vertexCollection.size();
1625 if (!gtReadoutRecord.
isValid()) {
1626 LogInfo(
"JetAnalyzer") <<
"JetAnalyzer: Could not find GT readout record" << std::endl;
1627 if (
verbose_)
std::cout <<
"JetAnalyzer: Could not find GT readout record product" << std::endl;
1630 bool techTriggerResultBxE =
false;
1631 bool techTriggerResultBxF =
false;
1632 bool techTriggerResultBx0 =
false;
1634 if (!gtReadoutRecord.
isValid()) {
1635 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find GT readout record" << std::endl;
1636 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find GT readout record product" << std::endl;
1639 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBxE = gtReadoutRecord->technicalTriggerWord(-2);
1640 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBxF = gtReadoutRecord->technicalTriggerWord(-1);
1641 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBx0 = gtReadoutRecord->technicalTriggerWord();
1646 if(techTriggerResultBx0!=0){
1675 bool pass_Z_selection=
false;
1689 for (
unsigned int i=0;
i<Muons->size();
i++){
1690 bool pass_muon_id=
false;
1691 bool pass_muon_iso=
false;
1692 double dxy=fabs((*Muons)[
i].muonBestTrack()->dxy());
1693 double dz=fabs((*Muons)[
i].muonBestTrack()->dz());
1695 dxy=fabs((*Muons)[
i].muonBestTrack()->dxy((*vertexHandle)[0].
position()));
1696 dz=fabs((*Muons)[
i].muonBestTrack()->dz((*vertexHandle)[0].
position()));
1698 if((*Muons)[
i].
pt()>20 && fabs((*Muons)[
i].
eta())<2.3){
1699 if((*Muons)[
i].isGlobalMuon() && (*Muons)[
i].isPFMuon() && (*Muons)[
i].globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && (*Muons)[
i].numberOfMatchedStations() > 1 && dxy < 0.2 && (*Muons)[
i].numberOfMatchedStations() > 1 && dz<0.5 && (*Muons)[
i].
innerTrack()->hitPattern().numberOfValidPixelHits() > 0 && (*Muons)[
i].
innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5){
1703 float muonIsoPFdb = ((*Muons)[
i].pfIsolationR04().sumChargedHadronPt
1704 +
std::max(0., (*Muons)[
i].pfIsolationR04().sumNeutralHadronEt + (*Muons)[
i].pfIsolationR04().sumPhotonEt - 0.5*(*Muons)[
i].pfIsolationR04().sumPUPt))/(*Muons)[
i].pt();
1705 if(muonIsoPFdb<0.12){
1709 if(pass_muon_id && pass_muon_iso){
1710 if((*Muons)[
i].
pt()>pt0){
1711 mu_index1=mu_index0;
1714 pt0=(*Muons)[
i].pt();
1715 }
else if ((*Muons)[
i].
pt()>pt1){
1717 pt1=(*Muons)[
i].pt();
1722 if(mu_index0>=0 && mu_index1>=0){
1723 if((*Muons)[mu_index0].
charge()*(*Muons)[mu_index1].
charge()<0){
1724 zCand=(*Muons)[mu_index0].polarP4()+(*Muons)[mu_index1].polarP4();
1725 if(fabs(zCand.M()-91.)<20 && zCand.Pt()>30){
1726 pass_Z_selection=
true;
1743 bool jetCollectionIsValid =
false;
1750 if (
jetCleaningFlag_ && (!jetCollectionIsValid || !bPrimaryVertex || !dcsDecision))
return;
1752 unsigned int collSize=-1;
1762 if (
isPFJet_) collSize=pfJets->size();
1767 std::vector<Jet> corJets;
1778 bool cleaned_first_jet=
false;
1779 bool cleaned_second_jet=
false;
1782 int ind1_mu_vetoed=-1;
1783 double pt1_mu_vetoed=-1;
1784 int ind2_mu_vetoed=-1;
1785 double pt2_mu_vetoed=-1;
1786 bool cleaned_first_jet_mu_vetoed=
false;
1787 bool cleaned_second_jet_mu_vetoed=
false;
1796 bool pass_correction_flag=
false;
1800 pass_correction_flag=
true;
1804 pass_correction_flag=
true;
1807 for (
unsigned int ijet=0; ijet<collSize; ijet++) {
1810 bool pass_uncorrected=
false;
1811 bool pass_corrected=
false;
1813 correctedJet=(*caloJets)[ijet];
1819 correctedJet=(*pfJets)[ijet];
1822 correctedJet=(*patJets)[ijet];
1825 pass_uncorrected=
true;
1828 pass_uncorrected=
true;
1832 scale = jetCorr->correction((*caloJets)[ijet]);
1835 scale = jetCorr->correction((*pfJets)[ijet]);
1840 pass_corrected=
true;
1846 recoJets.push_back(correctedJet);
1848 bool jetpassid=
true;
1849 bool Thiscleaned=
true;
1855 reco::JetID jetID = (*jetID_ValueMap_Handle)[calojetref];
1859 Thiscleaned=jetpassid;
1861 if(Thiscleaned && pass_corrected){
1869 Thiscleaned=jetpassid;
1871 if(Thiscleaned && pass_uncorrected){
1878 if(Thiscleaned && pass_corrected){
1883 if (fabs((*caloJets)[ijet].
eta()) <= 1.3) {
1886 }
else if(fabs((*caloJets)[ijet].
eta()) <3.0){
1907 Thiscleaned = jetpassid;
1908 JetIDWPU = jetpassid;
1910 if(Thiscleaned && pass_uncorrected){
1916 if(Thiscleaned && pass_corrected){
1917 if(fabs(correctedJet.
eta()) <= 1.3) {
1918 if(correctedJet.
pt()<=50.){
1921 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
1924 if(correctedJet.
pt()>140.){
1930 }
else if(fabs(correctedJet.
eta()) <= 3) {
1931 if(correctedJet.
pt()<=50.){
1934 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
1937 if(correctedJet.
pt()>140.){
1943 }
else if(fabs(correctedJet.
eta()) <= 5) {
1944 if(correctedJet.
pt()<=50.){
1947 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
1950 if(correctedJet.
pt()>140.){
1962 int puidmvaflag=-10;
1963 int puidcutflag=-10;
1964 puidmva=(*puJetIdMva)[pfjetref];
1965 puidcut=(*puJetId)[pfjetref];
1966 puidmvaflag=(*puJetIdFlagMva)[pfjetref];
1967 puidcutflag=(*puJetIdFlag)[pfjetref];
1983 Thiscleaned = jetpassid;
1986 if(Thiscleaned && pass_uncorrected){
1992 if(Thiscleaned && pass_corrected){
2035 if (correctedJet.
pt()<= 50) {
2040 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2045 if (correctedJet.
pt()>140.) {
2050 if (fabs(correctedJet.
eta()) <= 1.3) {
2052 if (correctedJet.
pt()<=50.) {
2072 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2092 if (correctedJet.
pt()>140.) {
2115 }
else if(fabs(correctedJet.
eta()) <= 3) {
2117 if (correctedJet.
pt()<=50.) {
2137 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2157 if (correctedJet.
pt()>140.) {
2184 if (correctedJet.
pt()<=50.) {
2198 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2212 if (correctedJet.
pt()>140.) {
2228 if(techTriggerResultBx0 && techTriggerResultBxE && techTriggerResultBxF){
2230 if(fabs(correctedJet.
eta()) <= 1.3) {
2235 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2240 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2245 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2249 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2255 if(techTriggerResultBx0 && techTriggerResultBxF){
2257 if(fabs(correctedJet.
eta()) <= 1.3) {
2262 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2267 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2272 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2276 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2282 if(techTriggerResultBx0 && !techTriggerResultBxE && !techTriggerResultBxF){
2284 if(fabs(correctedJet.
eta()) <= 1.3) {
2289 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2294 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2299 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2303 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2309 if(techTriggerResultBx0 && !techTriggerResultBxF){
2311 if(fabs(correctedJet.
eta()) <= 1.3) {
2316 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2321 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2326 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2330 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2359 if(pass_Z_selection){
2360 if(
deltaR((*Muons)[mu_index0].
eta(),(*Muons)[mu_index0].
phi(),correctedJet.
eta(),correctedJet.
phi())>0.2 &&
deltaR((*Muons)[mu_index1].eta(),(*Muons)[mu_index1].phi(),correctedJet.
eta(),correctedJet.
phi())>0.2 ){
2361 if(correctedJet.
pt()>pt1_mu_vetoed){
2362 pt2_mu_vetoed=pt1_mu_vetoed;
2363 ind2_mu_vetoed=ind1_mu_vetoed;
2364 cleaned_second_jet_mu_vetoed=cleaned_first_jet_mu_vetoed;
2365 pt1_mu_vetoed=correctedJet.
pt();
2366 ind1_mu_vetoed=ijet;
2367 cleaned_first_jet_mu_vetoed=JetIDWPU;
2368 }
else if(correctedJet.
pt()>pt2_mu_vetoed){
2369 pt2_mu_vetoed=correctedJet.
pt();
2370 ind2_mu_vetoed=ijet;
2371 cleaned_second_jet_mu_vetoed=JetIDWPU;
2376 if(correctedJet.
pt()>pt1){
2382 cleaned_second_jet=cleaned_first_jet;
2383 pt1=correctedJet.
pt();
2385 cleaned_first_jet=JetIDWPU;
2386 }
else if(correctedJet.
pt()>pt2){
2390 pt2=correctedJet.
pt();
2392 cleaned_second_jet=JetIDWPU;
2393 }
else if(correctedJet.
pt()>pt3){
2394 pt3=correctedJet.
pt();
2398 if(!pass_corrected){
2407 if(fabs(correctedJet.
eta())<3.0){
2413 if(fabs(correctedJet.
eta())<3.0){
2419 if(Thiscleaned && pass_corrected){
2470 if (fabs(correctedJet.
eta()) <= 1.3) {
2477 }
else if (fabs(correctedJet.
eta()) <= 3) {
2505 if(cleaned_first_jet){
2510 if(numofjets>1 && cleaned_second_jet) {
2511 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2513 dphi=2*acos(-1.)-dphi;
2526 if(
jetCleaningFlag_ && recoJets.size()>1 && cleaned_first_jet && cleaned_second_jet ){
2533 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2535 dphi=2*acos(-1.)-dphi;
2543 reco::JetID jetID1 = (*jetID_ValueMap_Handle)[calojetref1];
2545 reco::JetID jetID2 = (*jetID_ValueMap_Handle)[calojetref2];
2643 float QGLikelihood=-10;
2648 QGmulti=(*qgMultiplicity)[pfjetref];
2649 QGLikelihood=(*qgLikelihood)[pfjetref];
2650 QGptD=(*qgptD)[pfjetref];
2651 QGaxis2=(*qgaxis2)[pfjetref];
2652 if(fabs(recoJets[0].
eta())<1.3){
2654 if (recoJets[0].
pt()>=20. && recoJets[0].
pt()<=50.) {
2660 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2666 if (recoJets[0].
pt()>140.) {
2672 }
else if(fabs(recoJets[0].
eta())<3.0){
2673 if (recoJets[0].
pt()>20. && recoJets[0].
pt()<=50.) {
2679 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2685 if (recoJets[0].
pt()>140.) {
2692 if (recoJets[0].
pt()>20. && recoJets[0].
pt()<=50.) {
2698 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2704 if (recoJets[0].
pt()>140.) {
2712 QGmulti=(*qgMultiplicity)[pfjetref1];
2713 QGLikelihood=(*qgLikelihood)[pfjetref1];
2714 QGptD=(*qgptD)[pfjetref1];
2715 QGaxis2=(*qgaxis2)[pfjetref1];
2716 if(fabs(recoJets[1].
eta())<1.3){
2718 if (recoJets[1].
pt()>=20. && recoJets[1].
pt()<=50.) {
2724 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
2730 if (recoJets[1].
pt()>140.) {
2736 }
else if(fabs(recoJets[1].
eta())<3.0){
2737 if (recoJets[1].
pt()>20. && recoJets[1].
pt()<=50.) {
2743 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
2749 if (recoJets[1].
pt()>140.) {
2756 if (recoJets[1].
pt()>20. && recoJets[1].
pt()<=50.) {
2762 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
2768 if (recoJets[1].
pt()>140.) {
2860 if (fabs(recoJets[0].
eta() < 1.4)) {
2861 double pt_dijet = (recoJets[0].pt() + recoJets[1].pt())/2;
2867 int randJet =iEvent.
id().
event()%2;
2868 if (fabs(recoJets[1].
eta() < 1.4)) {
2884 bool thirdJetCut =
true;
2897 double dijetAsymmetry =(recoJets[jet1].pt() - recoJets[jet2].pt()) / (recoJets[jet1].
pt() + recoJets[jet2].pt());
2907 pt_barrel = recoJets[jet1].pt();
2908 pt_probe = recoJets[jet2].pt();
2918 bool thirdJetCut =
true;
2921 thirdJetCut =
false;
2925 double dijetBalance = (pt_probe - pt_barrel) / pt_dijet;
2934 if(pass_Z_selection && ind1_mu_vetoed>=0 && pt1_mu_vetoed>12 && cleaned_first_jet_mu_vetoed &&
isPFJet_){
2935 bool pass_second_jet_mu_vetoed=
false;
2936 if(cleaned_second_jet_mu_vetoed){
2937 if(ind2_mu_vetoed>=0 && pt2_mu_vetoed/zCand.Pt()<0.2){
2938 pass_second_jet_mu_vetoed=
true;
2941 if(pass_second_jet_mu_vetoed){
2944 recoJet1=(*pfJets)[ind1_mu_vetoed];
2949 scale = jetCorr->correction((*caloJets)[ind1_mu_vetoed]);
2952 scale = jetCorr->correction((*pfJets)[ind1_mu_vetoed]);
2956 double dphi=fabs(recoJet1.
phi()-zCand.Phi());
2958 dphi=2*acos(-1.)-dphi;
2964 if(fabs(dphi-acos(-1.))<0.34){
2972 if(!calometcoll.
isValid())
return;
2973 met=&(calometcoll->front());
2977 if(!pfmetcoll.
isValid())
return;
2978 met=&(pfmetcoll->front());
2987 if(recoJet1.
pt()>20){
3012 double MPF=1.+(met->
px()*zCand.Px()+met->
py()*zCand.Py())/(zCand.Pt()*zCand.Pt());
3013 if(fabs(recoJet1.
eta())<1.3){
3021 }
else if (zCand.Pt()<140){
3033 }
else if (zCand.Pt()<75){
3035 }
else if (zCand.Pt()<150){
3037 }
else if (zCand.Pt()<290){
3043 }
else if(fabs(recoJet1.
eta())<3.0){
3050 }
else if (zCand.Pt()<140){
3062 }
else if (zCand.Pt()<75){
3064 }
else if (zCand.Pt()<150){
3066 }
else if (zCand.Pt()<290){
3079 }
else if (zCand.Pt()<140){
3091 }
else if (zCand.Pt()<100){
3099 float QGLikelihood=-10;
3104 QGmulti=(*qgMultiplicity)[pfjetref];
3105 QGLikelihood=(*qgLikelihood)[pfjetref];
3106 QGptD=(*qgptD)[pfjetref];
3107 QGaxis2=(*qgaxis2)[pfjetref];
3108 if(fabs(recoJet1.
eta())<1.3){
3110 if (recoJet1.
pt()>=20. && recoJet1.
pt()<=50.) {
3116 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3122 if (recoJet1.
pt()>140.) {
3128 }
else if(fabs(recoJet1.
eta())<3.0){
3129 if (recoJet1.
pt()>20. && recoJet1.
pt()<=50.) {
3135 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3141 if (recoJet1.
pt()>140.) {
3148 if (recoJet1.
pt()>20. && recoJet1.
pt()<=50.) {
3154 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3160 if (recoJet1.
pt()>140.) {
MonitorElement * mqgLikelihood_lowPt_Barrel
edm::EDGetTokenT< edm::ValueMap< float > > qgaxis2Token_
MonitorElement * mCHEn_mediumPt_Barrel
MonitorElement * mEFrac_EndCap
MonitorElement * mNHEn_lowPt_Barrel
MonitorElement * mJetZBalance_lowZPt_J_Barrel
MonitorElement * mMass_lowPt_Forward
EventNumber_t event() const
MonitorElement * mMass_mediumPt_Forward
MonitorElement * mePtForwardPlus_BXm1Empty
MonitorElement * mMass_highPt_Forward
MonitorElement * mHadEnergyInHB
MonitorElement * mDeltaPt_Z_j1_over_ZPt_290_J_Barrel
MonitorElement * mHFEEn_highPt_Forward
MonitorElement * mDeltaPt_Z_j1_over_ZPt_30_55_J_Barrel
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
MonitorElement * mMultiplicityQG_highPt_Forward
MonitorElement * mCHEn_highPt_Barrel
MonitorElement * mPhi_Forward
MonitorElement * mChMultiplicity_highPt_Barrel
MonitorElement * mLooseCutPUJIDPassFractionVSeta
MonitorElement * mDeltaPt_Z_j1_over_ZPt_30_55_J_EndCap
MonitorElement * mJetZBalance_lowZPt_J_Forward
MonitorElement * mChMultiplicity_mediumPt_EndCap
MonitorElement * mNeutralEmEnergy
static bool jetSortingRule(reco::Jet x, reco::Jet y)
MonitorElement * meCHFracBarrel_BXm2BXm1Filled
MonitorElement * mDeltaPt_Z_j1_over_ZPt_55_100_J_Forward
MonitorElement * mPhEn_mediumPt_EndCap
virtual void scaleEnergy(double fScale)
scale energy of the jet
MonitorElement * mePhFracEndCapPlus_BXm1Filled
MonitorElement * mJetZBalance_highZPt_J_Barrel
MonitorElement * mEmEnergyInHF
MonitorElement * mePtForwardPlus_BXm2BXm1Filled
MonitorElement * mePtEndCapPlus_BXm2BXm1Empty
MonitorElement * mPhiFirst
MonitorElement * mCutPUJIDDiscriminant_lowPt_Forward
MonitorElement * mMPF_mediumZPt_J_Forward
MonitorElement * mNeutralHadronEnergy
MonitorElement * mHFEFrac_mediumPt_Forward
MonitorElement * mPt_Forward
MonitorElement * mCHFrac_highPt_EndCap
MonitorElement * mNHFrac_highPt_EndCap
bool fill_jet_high_level_histo
MonitorElement * mqgLikelihood_highPt_EndCap
MonitorElement * mChargedEmEnergy_profile
MonitorElement * mJ1Pt_over_ZPt_lowZPt_J_Barrel
edm::EDGetTokenT< reco::JetCorrector > jetCorrectorToken_
std::vector< std::string > lowPtJetExpr_
MonitorElement * meHFEMFracMinus_BXm2BXm1Empty
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * mHFrac_Forward
MonitorElement * mEmEnergyInEE
MonitorElement * mHFHFrac_highPt_Forward
MonitorElement * meNHFracBarrel_BXm1Filled
MonitorElement * mHFHEn_lowPt_Forward
MonitorElement * mMultiplicityQG_mediumPt_EndCap
MonitorElement * mNHFrac_lowPt_EndCap
MonitorElement * mCHFrac_lowPt_Barrel
edm::EDGetTokenT< edm::ValueMap< int > > cutBasedPUIDToken_
MonitorElement * meCHFracEndCapMinus_BXm2BXm1Empty
MonitorElement * mElEn_highPt_Barrel
edm::EDGetTokenT< reco::CaloJetCollection > caloJetsToken_
MonitorElement * mHFHFrac
MonitorElement * mPhEn_lowPt_EndCap
Base class for all types of Jets.
double balanceThirdJetCut_
MonitorElement * mMediumCutPUJIDPassFractionVSeta
MonitorElement * mHFEEn_lowPt_Forward
edm::EDGetTokenT< edm::ValueMap< float > > qgLikelihoodToken_
MonitorElement * mePhFracEndCapPlus_BXm2BXm1Empty
MonitorElement * meNHFracEndCapMinus_BXm1Filled
MonitorElement * mPhi_EndCap
JetIDSelectionFunctor jetIDFunctor
MonitorElement * mMPF_J_EndCap
MonitorElement * meCHFracEndCapMinus_BXm2BXm1Filled
MonitorElement * mqgLikelihood_mediumPt_Barrel
edm::EDGetTokenT< reco::PFMETCollection > pfMetToken_
MonitorElement * mPhFrac_lowPt_EndCap
MonitorElement * mMass_highPt_EndCap
MonitorElement * mMVAPUJIDDiscriminant_mediumPt_Barrel
MonitorElement * mNHFracVSeta_highPt
MonitorElement * mDeltaPt_Z_j1_over_ZPt_290_J_EndCap
MonitorElement * mMPF_highZPt_J_Forward
MonitorElement * mCutPUJIDDiscriminant_lowPt_Barrel
MonitorElement * mN90Hits
PFJetIDSelectionFunctor pfjetIDFunctor
MonitorElement * mJetZBalance_mediumZPt_J_Barrel
MonitorElement * mCHEn_lowPt_EndCap
edm::EDGetTokenT< reco::PFJetCollection > pfJetsToken_
MonitorElement * meNHFracEndCapPlus_BXm1Empty
MonitorElement * mCutPUJIDDiscriminant_lowPt_EndCap
std::vector< Vertex > VertexCollection
collection of Vertex objects
MonitorElement * mCHFracVSeta_highPt
MonitorElement * meCHFracEndCapPlus_BXm2BXm1Filled
MonitorElement * mPhEn_highPt_Barrel
MonitorElement * mCHFrac_profile
MonitorElement * mNeutMultiplicity_lowPt_EndCap
MonitorElement * mNHEn_mediumPt_EndCap
edm::EDGetTokenT< edm::ValueMap< int > > mvaPUIDToken_
MonitorElement * mpTD_mediumPt_EndCap
MonitorElement * meCHFracEndCapPlus_BXm1Filled
MonitorElement * mMPF_lowZPt_J_EndCap
MonitorElement * mNeutralMultiplicity
Strings::size_type size() const
edm::EDGetTokenT< edm::ValueMap< float > > mvaFullPUDiscriminantToken_
MonitorElement * mePtForwardMinus_BXm2BXm1Filled
MonitorElement * mMass_lowPt_Barrel
MonitorElement * mEFrac_profile
MonitorElement * meCHFracBarrel_BXm2BXm1Empty
MonitorElement * mePhFracBarrel_BXm1Filled
MonitorElement * mDeltaPt_Z_j1_over_ZPt_150_290_J_EndCap
MonitorElement * mHFEMFrac
MonitorElement * mConstituents_Forward
edm::EDGetTokenT< std::vector< reco::Vertex > > vertexToken_
MonitorElement * mMVAPUJIDDiscriminant_lowPt_Forward
MonitorElement * mpTD_lowPt_Forward
JetMETDQMDCSFilter * DCSFilterForDCSMonitoring_
MonitorElement * mePhFracEndCapPlus_BXm1Empty
MonitorElement * mMPF_highZPt_J_Barrel
MonitorElement * mMultiplicityQG_lowPt_EndCap
MonitorElement * meHFHFracPlus_BXm1Filled
MonitorElement * meNHFracEndCapPlus_BXm1Filled
MonitorElement * mCHFrac_lowPt_EndCap
virtual double eta() const
momentum pseudorapidity
virtual double pt() const
transverse momentum
MonitorElement * mMuMultiplicity_highPt_Barrel
MonitorElement * mNHFracVSeta_mediumPt
MonitorElement * mChargedMultiplicity_profile
MonitorElement * mAxis2_mediumPt_Barrel
MonitorElement * mHFEFrac_lowPt_Forward
MonitorElement * meHFEMFracPlus_BXm1Filled
MonitorElement * mNHFracVSeta_lowPt
MonitorElement * mLooseJIDPassFractionVSpt
MonitorElement * mAxis2_highPt_Forward
JetMETDQMDCSFilter * DCSFilterForJetMonitoring_
MonitorElement * mpTD_mediumPt_Barrel
MonitorElement * mPt_uncor
MonitorElement * meEta_BXm2BXm1Filled
MonitorElement * mDeltaPt_Z_j1_over_ZPt_75_150_J_Barrel
MonitorElement * mChargedMuEnergy
MonitorElement * mLooseMVAPUJIDPassFractionVSpt
MonitorElement * meHFHFracMinus_BXm2BXm1Filled
MonitorElement * mMuMultiplicity_mediumPt_EndCap
MonitorElement * mNHFrac_profile
MonitorElement * mConstituents_profile
MonitorElement * mpTD_highPt_EndCap
MonitorElement * mEmEnergyInEB
MonitorElement * mDijetBalance
MonitorElement * meNHFracEndCapPlus_BXm2BXm1Empty
MonitorElement * mDeltaPt_Z_j1_over_ZPt_55_75_J_EndCap
MonitorElement * mJ1Pt_over_ZPt_J_Barrel
MonitorElement * mJ1Pt_over_ZPt_J_EndCap
static bool passJetId(int flag, Id level)
MonitorElement * meHFEMFracPlus_BXm2BXm1Empty
MonitorElement * mTightCutPUJIDPassFractionVSeta
MonitorElement * mMediumMVAPUJIDPassFractionVSeta
MonitorElement * mePtEndCapPlus_BXm1Empty
MonitorElement * mChargedHadronEnergy_profile
MonitorElement * mPhi_Barrel
MonitorElement * mPhFrac_lowPt_Barrel
MonitorElement * mAxis2_highPt_Barrel
MonitorElement * mMass_lowPt_EndCap
MonitorElement * mePtForwardMinus_BXm2BXm1Empty
edm::EDGetTokenT< reco::CaloMETCollection > caloMetToken_
MonitorElement * mConstituents_EndCap
MonitorElement * meEta_BXm1Empty
MonitorElement * mHFHFrac_mediumPt_Forward
MonitorElement * mPhFrac_mediumPt_Barrel
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > gtToken_
MonitorElement * mJetZBalance_mediumZPt_J_Forward
MonitorElement * mePhFracBarrel_BXm2BXm1Filled
MonitorElement * mePtEndCapMinus_BXm1Empty
MonitorElement * mNeutralMultiplicity_profile
MonitorElement * mMultiplicityQG_highPt_EndCap
MonitorElement * meHFHFracMinus_BXm1Filled
MonitorElement * cleanupME
MonitorElement * mPhEn_mediumPt_Barrel
MonitorElement * mLooseMVAPUJIDPassFractionVSeta
MonitorElement * mNeutMultiplicity_mediumPt_EndCap
MonitorElement * mNeutralFraction
MonitorElement * mNHFracVSpT_EndCap
MonitorElement * mNeutralEmEnergy_profile
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
MonitorElement * mNHFrac_mediumPt_Barrel
MonitorElement * mPhEn_lowPt_Barrel
edm::EDGetTokenT< edm::ValueMap< float > > cutBasedPUDiscriminantToken_
edm::EDGetTokenT< pat::JetCollection > patJetsToken_
MonitorElement * mConstituents_Barrel
MonitorElement * mJetZBalance_lowZPt_J_EndCap
MonitorElement * meHFHFracMinus_BXm1Empty
MonitorElement * meCHFracEndCapMinus_BXm1Empty
MonitorElement * mNHEn_lowPt_EndCap
MonitorElement * mMuEn_lowPt_EndCap
MonitorElement * mePhFracBarrel_BXm2BXm1Empty
MonitorElement * meHFEMFracMinus_BXm1Filled
MonitorElement * mNeutMultiplicity_lowPt_Barrel
MonitorElement * mpTD_highPt_Barrel
MonitorElement * mPhFrac_mediumPt_EndCap
MonitorElement * mJ1Pt_over_ZPt_highZPt_J_Forward
MonitorElement * mMuonMultiplicity_profile
MonitorElement * mCutPUJIDDiscriminant_mediumPt_Forward
MonitorElement * mCutPUJIDDiscriminant_highPt_EndCap
MonitorElement * mChargedHadronEnergy
MonitorElement * mePhFracEndCapPlus_BXm2BXm1Filled
MonitorElement * mPhi_uncor
MonitorElement * mePhFracEndCapMinus_BXm2BXm1Empty
MonitorElement * mCutPUJIDDiscriminant_highPt_Barrel
MonitorElement * meHFEMFracPlus_BXm2BXm1Filled
MonitorElement * meNHFracBarrel_BXm2BXm1Filled
MonitorElement * mChargedEmEnergy
MonitorElement * mNeutMultiplicity_mediumPt_Forward
MonitorElement * mMass_mediumPt_EndCap
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
MonitorElement * mPhFrac_profile
MonitorElement * meNHFracBarrel_BXm1Empty
MonitorElement * mJ1Pt_over_ZPt_highZPt_J_Barrel
MonitorElement * mPt_profile
edm::EDGetTokenT< edm::ValueMap< float > > qgptDToken_
MonitorElement * mPhFracVSeta_lowPt
MonitorElement * mChMultiplicity_lowPt_Barrel
edm::InputTag jetCorrectorTag_
MonitorElement * mDeltaPt_Z_j1_over_ZPt_55_75_J_Barrel
MonitorElement * mEtaFirst
MonitorElement * mElEn_lowPt_EndCap
MonitorElement * mTightMVAPUJIDPassFractionVSeta
MonitorElement * mNHFrac_mediumPt_EndCap
MonitorElement * mEFrac_Forward
static std::string const triggerResults
MonitorElement * mDeltaPt_Z_j1_over_ZPt_75_150_J_EndCap
MonitorElement * mChMultiplicity_lowPt_EndCap
MonitorElement * mMuEn_mediumPt_Barrel
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken_
MonitorElement * mMuMultiplicity_highPt_EndCap
MonitorElement * mePtEndCapMinus_BXm2BXm1Empty
MonitorElement * mAxis2_highPt_EndCap
MonitorElement * mMVAPUJIDDiscriminant_highPt_Barrel
MonitorElement * mNeutralHadronEnergy_profile
MonitorElement * mHFrac_Barrel
bool filter(const edm::Event &evt, const edm::EventSetup &es)
MonitorElement * mNJets_profile
MonitorElement * mMultiplicityQG_highPt_Barrel
MonitorElement * mHFEFracVSpT_Forward
MonitorElement * mMPF_mediumZPt_J_Barrel
MonitorElement * mePhFracEndCapMinus_BXm2BXm1Filled
MonitorElement * mePhFracBarrel_BXm1Empty
MonitorElement * mCHEn_lowPt_Barrel
MonitorElement * mHFEFrac_highPt_Forward
MonitorElement * mJ1Pt_over_ZPt_J_Forward
MonitorElement * mElEn_highPt_EndCap
MonitorElement * mCutPUJIDDiscriminant_mediumPt_EndCap
MonitorElement * mJ1Pt_over_ZPt_mediumZPt_J_Forward
MonitorElement * mJ1Pt_over_ZPt_mediumZPt_J_EndCap
MonitorElement * meNHFracEndCapMinus_BXm2BXm1Empty
MonitorElement * mCutPUJIDDiscriminant_mediumPt_Barrel
MonitorElement * mPhFrac_highPt_EndCap
MonitorElement * mPhFracVSeta_mediumPt
MonitorElement * mMPF_lowZPt_J_Barrel
MonitorElement * meHFEMFracPlus_BXm1Empty
double deltaR(double eta1, double eta2, double phi1, double phi2)
MonitorElement * mePtBarrel_BXm2BXm1Filled
MonitorElement * meCHFracBarrel_BXm1Empty
MonitorElement * mePtBarrel_BXm1Empty
MonitorElement * mNHFracVSpT_Barrel
MonitorElement * mHadEnergyInHF
edm::EDGetTokenT< reco::MuonCollection > MuonsToken_
MonitorElement * mMPF_lowZPt_J_Forward
MonitorElement * mDijetAsymmetry
MonitorElement * mNeutMultiplicity_highPt_EndCap
MonitorElement * mpTD_lowPt_Barrel
MonitorElement * mHFHEn_mediumPt_Forward
MonitorElement * mJetEnergyCorr
MonitorElement * mJetZBalance_highZPt_J_Forward
MonitorElement * mAxis2_lowPt_EndCap
MonitorElement * mElEn_mediumPt_Barrel
MonitorElement * mMVAPUJIDDiscriminant_highPt_Forward
MonitorElement * mePtEndCapMinus_BXm2BXm1Filled
TObject * getRootObject(void) const
edm::EDGetTokenT< edm::ValueMap< reco::JetID > > jetID_ValueMapToken_
MonitorElement * mChargedMultiplicity
MonitorElement * mJ1Pt_over_ZPt_highZPt_J_EndCap
edm::EDGetTokenT< edm::ValueMap< int > > qgMultiplicityToken_
MonitorElement * mNHEn_highPt_Barrel
std::map< std::string, MonitorElement * > map_of_MEs
MonitorElement * mEFrac_Barrel
double asymmetryThirdJetCut_
MonitorElement * mCHFracVSeta_mediumPt
T const * product() const
MonitorElement * mePtBarrel_BXm1Filled
MonitorElement * mCHFrac_mediumPt_EndCap
virtual double px() const
x coordinate of momentum vector
MonitorElement * mNHFrac_highPt_Barrel
MonitorElement * mqgLikelihood_highPt_Forward
MonitorElement * mMediumCutPUJIDPassFractionVSpt
MonitorElement * mePhFracEndCapMinus_BXm1Empty
MonitorElement * mMuMultiplicity_mediumPt_Barrel
std::string const & triggerName(unsigned int index) const
MonitorElement * mHFrac_EndCap
MonitorElement * mJetEnergyCorrVSPt
MonitorElement * mElEn_mediumPt_EndCap
MonitorElement * mAxis2_lowPt_Barrel
MonitorElement * mPhiVSEta
MonitorElement * mMPF_mediumZPt_J_EndCap
MonitorElement * mqgLikelihood_mediumPt_EndCap
MonitorElement * meHFHFracPlus_BXm2BXm1Filled
MonitorElement * mePtEndCapMinus_BXm1Filled
MonitorElement * mChargedMuEnergy_profile
MonitorElement * mLooseJIDPassFractionVSptNoHF
MonitorElement * mEta_uncor
MonitorElement * mMediumMVAPUJIDPassFractionVSpt
MonitorElement * mCHEn_highPt_EndCap
MonitorElement * mMVAPUJIDDiscriminant_highPt_EndCap
MonitorElement * mJ1Pt_over_ZPt_lowZPt_J_EndCap
MonitorElement * mMPF_J_Forward
MonitorElement * mCHFracVSpT_EndCap
MonitorElement * mHFHFracVSpT_Forward
MonitorElement * mMVAPUJIDDiscriminant_lowPt_Barrel
MonitorElement * mpTD_lowPt_EndCap
MonitorElement * mNHEn_mediumPt_Barrel
MonitorElement * mJetZBalance_mediumZPt_J_EndCap
MonitorElement * mMuMultiplicity_lowPt_Barrel
MonitorElement * mePtForwardPlus_BXm2BXm1Empty
pat::strbitset getBitTemplate() const
Get an empty bitset with the proper names.
MonitorElement * mMuEn_highPt_EndCap
MonitorElement * mLooseJIDPassFractionVSeta
MonitorElement * mHFrac_profile
MonitorElement * mMPF_highZPt_J_EndCap
MonitorElement * mPhEn_highPt_EndCap
MonitorElement * mMultiplicityQG_lowPt_Barrel
MonitorElement * mMVAPUJIDDiscriminant_lowPt_EndCap
MonitorElement * meHFHFracMinus_BXm2BXm1Empty
MonitorElement * mMuEn_mediumPt_EndCap
MonitorElement * mqgLikelihood_lowPt_EndCap
MonitorElement * mAxis2_mediumPt_Forward
MonitorElement * mMuEn_lowPt_Barrel
virtual int nConstituents() const
of constituents
MonitorElement * mTightMVAPUJIDPassFractionVSpt
MonitorElement * mCHFrac_highPt_Barrel
MonitorElement * mMuEn_highPt_Barrel
MonitorElement * meCHFracEndCapPlus_BXm1Empty
static int position[264][3]
MonitorElement * mePhFracEndCapMinus_BXm1Filled
MonitorElement * mHFHEn_highPt_Forward
MonitorElement * mePtEndCapPlus_BXm1Filled
MonitorElement * mNHEn_highPt_EndCap
MonitorElement * mePtForwardPlus_BXm1Filled
MonitorElement * mCHFrac_mediumPt_Barrel
MonitorElement * mAxis2_lowPt_Forward
MonitorElement * mMultiplicityQG_mediumPt_Barrel
edm::InputTag theTriggerResultsLabel_
MonitorElement * mPtFirst
MonitorElement * mDeltaPt_Z_j1_over_ZPt_30_55_J_Forward
MonitorElement * meNHFracEndCapMinus_BXm1Empty
MonitorElement * mPhi_profile
MonitorElement * mPt_Barrel
MonitorElement * meCHFracBarrel_BXm1Filled
MonitorElement * mElEn_lowPt_Barrel
MonitorElement * mMass_highPt_Barrel
MonitorElement * meNHFracEndCapMinus_BXm2BXm1Filled
MonitorElement * mHadEnergyInHE
std::vector< std::string > highPtJetExpr_
MonitorElement * mMass_mediumPt_Barrel
MonitorElement * meHFEMFracMinus_BXm1Empty
MonitorElement * mChMultiplicity_highPt_EndCap
MonitorElement * meHFHFracPlus_BXm2BXm1Empty
MonitorElement * mConstituents
MonitorElement * mNHFrac_lowPt_Barrel
MonitorElement * meEta_BXm1Filled
MonitorElement * mJetEnergyCorrVSEta
MonitorElement * mPhFrac_highPt_Barrel
MonitorElement * mCHFracVSpT_Barrel
MonitorElement * mPhFracVSpT_EndCap
MonitorElement * mChMultiplicity_mediumPt_Barrel
MonitorElement * mMPF_J_Barrel
MonitorElement * mJetZBalance_highZPt_J_EndCap
MonitorElement * mMultiplicityQG_lowPt_Forward
MonitorElement * meCHFracEndCapPlus_BXm2BXm1Empty
MonitorElement * mEta_profile
MonitorElement * mPt_EndCap
MonitorElement * mCHFracVSeta_lowPt
MonitorElement * meEta_BXm2BXm1Empty
MonitorElement * mpTD_mediumPt_Forward
MonitorElement * mMVAPUJIDDiscriminant_mediumPt_EndCap
MonitorElement * mCHEn_mediumPt_EndCap
MonitorElement * mTightCutPUJIDPassFractionVSpt
MonitorElement * mNeutMultiplicity_highPt_Forward
MonitorElement * mDPhiZJet
MonitorElement * mqgLikelihood_highPt_Barrel
MonitorElement * mqgLikelihood_mediumPt_Forward
virtual double phi() const
momentum azimuthal angle
MonitorElement * mePtForwardMinus_BXm1Filled
MonitorElement * mLooseCutPUJIDPassFractionVSpt
MonitorElement * mMuonMultiplicity
MonitorElement * mPhFracVSpT_Barrel
MonitorElement * mHFEMFrac_profile
MonitorElement * mePtEndCapPlus_BXm2BXm1Filled
MonitorElement * mMultiplicityQG_mediumPt_Forward
MonitorElement * mqgLikelihood_lowPt_Forward
MonitorElement * mDeltaPt_Z_j1_over_ZPt_100_J_Forward
MonitorElement * mHFHFrac_lowPt_Forward
MonitorElement * mMVAPUJIDDiscriminant_mediumPt_Forward
MonitorElement * mJ1Pt_over_ZPt_lowZPt_J_Forward
MonitorElement * mMuMultiplicity_lowPt_EndCap
MonitorElement * mZJetAsymmetry
MonitorElement * meHFHFracPlus_BXm1Empty
virtual double py() const
y coordinate of momentum vector
MonitorElement * mConstituents_uncor
MonitorElement * mNeutMultiplicity_mediumPt_Barrel
MonitorElement * meCHFracEndCapMinus_BXm1Filled
MonitorElement * mCutPUJIDDiscriminant_highPt_Forward
MonitorElement * mNeutMultiplicity_lowPt_Forward
MonitorElement * mHFEEn_mediumPt_Forward
MonitorElement * meHFEMFracMinus_BXm2BXm1Filled
MonitorElement * mNeutMultiplicity_highPt_Barrel
MonitorElement * verticesME
MonitorElement * mPhFracVSeta_highPt
MonitorElement * mAxis2_mediumPt_EndCap
MonitorElement * mePtBarrel_BXm2BXm1Empty
MonitorElement * mpTD_highPt_Forward
MonitorElement * mePtForwardMinus_BXm1Empty
MonitorElement * meNHFracBarrel_BXm2BXm1Empty
MonitorElement * meNHFracEndCapPlus_BXm2BXm1Filled
MonitorElement * mHFHFrac_profile
MonitorElement * mJ1Pt_over_ZPt_mediumZPt_J_Barrel
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.
MonitorElement * mDeltaPt_Z_j1_over_ZPt_150_290_J_Barrel
MonitorElement * mHadEnergyInHO
edm::InputTag mInputCollection_
Public Member Functions inherited from DQMEDAnalyzer
1.8.5
