Get the analysis.
1580 Int_t JetLoPass = 0;
1581 Int_t JetHiPass = 0;
1583 if (triggerResults.
isValid()){
1586 const unsigned int nTrig(triggerNames.
size());
1587 for (
unsigned int i=0;
i<nTrig;++
i)
1601 std::cout <<
">>> Trigger Lo = " << JetLoPass
1602 <<
" Hi = " << JetHiPass
1611 if (!vertexHandle.
isValid()) {
1612 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1613 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1616 if ( vertexHandle.
isValid() ){
1618 numPV = vertexCollection.size();
1628 if (!gtReadoutRecord.
isValid()) {
1629 LogInfo(
"JetAnalyzer") <<
"JetAnalyzer: Could not find GT readout record" << std::endl;
1630 if (
verbose_)
std::cout <<
"JetAnalyzer: Could not find GT readout record product" << std::endl;
1633 bool techTriggerResultBxE =
false;
1634 bool techTriggerResultBxF =
false;
1635 bool techTriggerResultBx0 =
false;
1637 if (!gtReadoutRecord.
isValid()) {
1638 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find GT readout record" << std::endl;
1639 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find GT readout record product" << std::endl;
1642 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBxE = gtReadoutRecord->technicalTriggerWord(-2);
1643 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBxF = gtReadoutRecord->technicalTriggerWord(-1);
1644 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBx0 = gtReadoutRecord->technicalTriggerWord();
1649 if(techTriggerResultBx0!=0){
1678 bool pass_Z_selection=
false;
1692 for (
unsigned int i=0;
i<Muons->size();
i++){
1693 bool pass_muon_id=
false;
1694 bool pass_muon_iso=
false;
1695 double dxy=fabs((*Muons)[
i].muonBestTrack()->dxy());
1696 double dz=fabs((*Muons)[
i].muonBestTrack()->dz());
1698 dxy=fabs((*Muons)[
i].muonBestTrack()->dxy((*vertexHandle)[0].
position()));
1699 dz=fabs((*Muons)[
i].muonBestTrack()->dz((*vertexHandle)[0].
position()));
1701 if((*Muons)[
i].
pt()>20 && fabs((*Muons)[
i].
eta())<2.3){
1702 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){
1706 float muonIsoPFdb = ((*Muons)[
i].pfIsolationR04().sumChargedHadronPt
1707 +
std::max(0., (*Muons)[
i].pfIsolationR04().sumNeutralHadronEt + (*Muons)[
i].pfIsolationR04().sumPhotonEt - 0.5*(*Muons)[
i].pfIsolationR04().sumPUPt))/(*Muons)[
i].pt();
1708 if(muonIsoPFdb<0.12){
1712 if(pass_muon_id && pass_muon_iso){
1713 if((*Muons)[
i].
pt()>pt0){
1714 mu_index1=mu_index0;
1717 pt0=(*Muons)[
i].pt();
1718 }
else if ((*Muons)[
i].
pt()>pt1){
1720 pt1=(*Muons)[
i].pt();
1725 if(mu_index0>=0 && mu_index1>=0){
1726 if((*Muons)[mu_index0].
charge()*(*Muons)[mu_index1].
charge()<0){
1727 zCand=(*Muons)[mu_index0].polarP4()+(*Muons)[mu_index1].polarP4();
1728 if(fabs(zCand.M()-91.)<20 && zCand.Pt()>30){
1729 pass_Z_selection=
true;
1746 bool jetCollectionIsValid =
false;
1753 if (
jetCleaningFlag_ && (!jetCollectionIsValid || !bPrimaryVertex || !dcsDecision))
return;
1755 unsigned int collSize=-1;
1765 if (
isPFJet_) collSize=pfJets->size();
1770 std::vector<Jet> corJets;
1781 bool cleaned_first_jet=
false;
1782 bool cleaned_second_jet=
false;
1785 int ind1_mu_vetoed=-1;
1786 double pt1_mu_vetoed=-1;
1787 int ind2_mu_vetoed=-1;
1788 double pt2_mu_vetoed=-1;
1789 bool cleaned_first_jet_mu_vetoed=
false;
1790 bool cleaned_second_jet_mu_vetoed=
false;
1799 bool pass_correction_flag=
false;
1803 pass_correction_flag=
true;
1807 pass_correction_flag=
true;
1810 for (
unsigned int ijet=0; ijet<collSize; ijet++) {
1813 bool pass_uncorrected=
false;
1814 bool pass_corrected=
false;
1816 correctedJet=(*caloJets)[ijet];
1822 correctedJet=(*pfJets)[ijet];
1825 correctedJet=(*patJets)[ijet];
1828 pass_uncorrected=
true;
1831 pass_uncorrected=
true;
1835 scale = jetCorr->correction((*caloJets)[ijet]);
1838 scale = jetCorr->correction((*pfJets)[ijet]);
1843 pass_corrected=
true;
1849 recoJets.push_back(correctedJet);
1852 if(!pass_corrected && !pass_uncorrected){
1855 bool jetpassid=
true;
1856 bool Thiscleaned=
true;
1862 reco::JetID jetID = (*jetID_ValueMap_Handle)[calojetref];
1866 Thiscleaned=jetpassid;
1868 if(Thiscleaned && pass_corrected){
1876 Thiscleaned=jetpassid;
1878 if(Thiscleaned && pass_uncorrected){
1885 if(Thiscleaned && pass_corrected){
1890 if (fabs((*caloJets)[ijet].
eta()) <= 1.3) {
1893 }
else if(fabs((*caloJets)[ijet].
eta()) <3.0){
1914 Thiscleaned = jetpassid;
1915 JetIDWPU = jetpassid;
1917 if(Thiscleaned && pass_uncorrected){
1923 if(Thiscleaned && pass_corrected){
1924 if(fabs(correctedJet.
eta()) <= 1.3) {
1925 if(correctedJet.
pt()<=50.){
1928 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
1931 if(correctedJet.
pt()>140.){
1937 }
else if(fabs(correctedJet.
eta()) <= 3) {
1938 if(correctedJet.
pt()<=50.){
1941 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
1944 if(correctedJet.
pt()>140.){
1950 }
else if(fabs(correctedJet.
eta()) <= 5) {
1951 if(correctedJet.
pt()<=50.){
1954 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
1957 if(correctedJet.
pt()>140.){
1969 int puidmvaflag=-10;
1970 int puidcutflag=-10;
1972 puidmva=(*puJetIdMva)[pfjetref];
1973 puidcut=(*puJetId)[pfjetref];
1974 puidmvaflag=(*puJetIdFlagMva)[pfjetref];
1975 puidcutflag=(*puJetIdFlag)[pfjetref];
1981 Thiscleaned = jetpassid;
1984 if(Thiscleaned && pass_uncorrected){
1990 if(Thiscleaned && pass_corrected){
2033 if (correctedJet.
pt()<= 50) {
2038 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2043 if (correctedJet.
pt()>140.) {
2048 if (fabs(correctedJet.
eta()) <= 1.3) {
2050 if (correctedJet.
pt()<=50.) {
2070 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2090 if (correctedJet.
pt()>140.) {
2113 }
else if(fabs(correctedJet.
eta()) <= 3) {
2115 if (correctedJet.
pt()<=50.) {
2135 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2155 if (correctedJet.
pt()>140.) {
2182 if (correctedJet.
pt()<=50.) {
2196 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2210 if (correctedJet.
pt()>140.) {
2226 if(techTriggerResultBx0 && techTriggerResultBxE && techTriggerResultBxF){
2228 if(fabs(correctedJet.
eta()) <= 1.3) {
2233 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2238 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2243 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2247 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2253 if(techTriggerResultBx0 && techTriggerResultBxF){
2255 if(fabs(correctedJet.
eta()) <= 1.3) {
2260 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2265 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2270 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2274 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2280 if(techTriggerResultBx0 && !techTriggerResultBxE && !techTriggerResultBxF){
2282 if(fabs(correctedJet.
eta()) <= 1.3) {
2287 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2292 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2297 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2301 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2307 if(techTriggerResultBx0 && !techTriggerResultBxF){
2309 if(fabs(correctedJet.
eta()) <= 1.3) {
2314 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2319 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2324 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2328 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2357 if(pass_Z_selection){
2358 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 ){
2359 if(correctedJet.
pt()>pt1_mu_vetoed){
2360 pt2_mu_vetoed=pt1_mu_vetoed;
2361 ind2_mu_vetoed=ind1_mu_vetoed;
2362 cleaned_second_jet_mu_vetoed=cleaned_first_jet_mu_vetoed;
2363 pt1_mu_vetoed=correctedJet.
pt();
2364 ind1_mu_vetoed=ijet;
2365 cleaned_first_jet_mu_vetoed=JetIDWPU;
2366 }
else if(correctedJet.
pt()>pt2_mu_vetoed){
2367 pt2_mu_vetoed=correctedJet.
pt();
2368 ind2_mu_vetoed=ijet;
2369 cleaned_second_jet_mu_vetoed=JetIDWPU;
2374 if(correctedJet.
pt()>pt1){
2380 cleaned_second_jet=cleaned_first_jet;
2381 pt1=correctedJet.
pt();
2383 cleaned_first_jet=JetIDWPU;
2384 }
else if(correctedJet.
pt()>pt2){
2388 pt2=correctedJet.
pt();
2390 cleaned_second_jet=JetIDWPU;
2391 }
else if(correctedJet.
pt()>pt3){
2392 pt3=correctedJet.
pt();
2396 if(!pass_corrected){
2405 if(fabs(correctedJet.
eta())<3.0){
2411 if(fabs(correctedJet.
eta())<3.0){
2417 if(Thiscleaned && pass_corrected){
2468 if (fabs(correctedJet.
eta()) <= 1.3) {
2475 }
else if (fabs(correctedJet.
eta()) <= 3) {
2503 if(cleaned_first_jet){
2508 if(numofjets>1 && cleaned_second_jet) {
2509 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2511 dphi=2*acos(-1.)-dphi;
2524 if(
jetCleaningFlag_ && recoJets.size()>1 && cleaned_first_jet && cleaned_second_jet ){
2531 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2533 dphi=2*acos(-1.)-dphi;
2541 reco::JetID jetID1 = (*jetID_ValueMap_Handle)[calojetref1];
2543 reco::JetID jetID2 = (*jetID_ValueMap_Handle)[calojetref2];
2641 float QGLikelihood=-10;
2646 QGmulti=(*qgMultiplicity)[pfjetref];
2647 QGLikelihood=(*qgLikelihood)[pfjetref];
2648 QGptD=(*qgptD)[pfjetref];
2649 QGaxis2=(*qgaxis2)[pfjetref];
2650 if(fabs(recoJets[0].
eta())<1.3){
2652 if (recoJets[0].
pt()>=20. && recoJets[0].
pt()<=50.) {
2658 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2664 if (recoJets[0].
pt()>140.) {
2670 }
else if(fabs(recoJets[0].
eta())<3.0){
2671 if (recoJets[0].
pt()>20. && recoJets[0].
pt()<=50.) {
2677 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2683 if (recoJets[0].
pt()>140.) {
2690 if (recoJets[0].
pt()>20. && recoJets[0].
pt()<=50.) {
2696 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2702 if (recoJets[0].
pt()>140.) {
2710 QGmulti=(*qgMultiplicity)[pfjetref1];
2711 QGLikelihood=(*qgLikelihood)[pfjetref1];
2712 QGptD=(*qgptD)[pfjetref1];
2713 QGaxis2=(*qgaxis2)[pfjetref1];
2714 if(fabs(recoJets[1].
eta())<1.3){
2716 if (recoJets[1].
pt()>=20. && recoJets[1].
pt()<=50.) {
2722 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
2728 if (recoJets[1].
pt()>140.) {
2734 }
else if(fabs(recoJets[1].
eta())<3.0){
2735 if (recoJets[1].
pt()>20. && recoJets[1].
pt()<=50.) {
2741 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
2747 if (recoJets[1].
pt()>140.) {
2754 if (recoJets[1].
pt()>20. && recoJets[1].
pt()<=50.) {
2760 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
2766 if (recoJets[1].
pt()>140.) {
2858 if (fabs(recoJets[0].
eta() < 1.4)) {
2859 double pt_dijet = (recoJets[0].pt() + recoJets[1].pt())/2;
2865 int randJet =iEvent.
id().
event()%2;
2866 if (fabs(recoJets[1].
eta() < 1.4)) {
2882 bool thirdJetCut =
true;
2895 double dijetAsymmetry =(recoJets[jet1].pt() - recoJets[jet2].pt()) / (recoJets[jet1].
pt() + recoJets[jet2].pt());
2905 pt_barrel = recoJets[jet1].pt();
2906 pt_probe = recoJets[jet2].pt();
2916 bool thirdJetCut =
true;
2919 thirdJetCut =
false;
2923 double dijetBalance = (pt_probe - pt_barrel) / pt_dijet;
2932 if(pass_Z_selection && ind1_mu_vetoed>=0 && pt1_mu_vetoed>12 && cleaned_first_jet_mu_vetoed &&
isPFJet_){
2933 bool pass_second_jet_mu_vetoed=
false;
2934 if(cleaned_second_jet_mu_vetoed){
2935 if(ind2_mu_vetoed>=0 && pt2_mu_vetoed/zCand.Pt()<0.2){
2936 pass_second_jet_mu_vetoed=
true;
2939 if(pass_second_jet_mu_vetoed){
2942 recoJet1=(*pfJets)[ind1_mu_vetoed];
2947 scale = jetCorr->correction((*caloJets)[ind1_mu_vetoed]);
2950 scale = jetCorr->correction((*pfJets)[ind1_mu_vetoed]);
2954 double dphi=fabs(recoJet1.
phi()-zCand.Phi());
2956 dphi=2*acos(-1.)-dphi;
2962 if(fabs(dphi-acos(-1.))<0.34){
2970 if(!calometcoll.
isValid())
return;
2971 met=&(calometcoll->front());
2975 if(!pfmetcoll.
isValid())
return;
2976 met=&(pfmetcoll->front());
2985 if(recoJet1.
pt()>20){
3010 double MPF=1.+(met->
px()*zCand.Px()+met->
py()*zCand.Py())/(zCand.Pt()*zCand.Pt());
3011 if(fabs(recoJet1.
eta())<1.3){
3019 }
else if (zCand.Pt()<140){
3031 }
else if (zCand.Pt()<75){
3033 }
else if (zCand.Pt()<150){
3035 }
else if (zCand.Pt()<290){
3041 }
else if(fabs(recoJet1.
eta())<3.0){
3048 }
else if (zCand.Pt()<140){
3060 }
else if (zCand.Pt()<75){
3062 }
else if (zCand.Pt()<150){
3064 }
else if (zCand.Pt()<290){
3077 }
else if (zCand.Pt()<140){
3089 }
else if (zCand.Pt()<100){
3097 float QGLikelihood=-10;
3102 QGmulti=(*qgMultiplicity)[pfjetref];
3103 QGLikelihood=(*qgLikelihood)[pfjetref];
3104 QGptD=(*qgptD)[pfjetref];
3105 QGaxis2=(*qgaxis2)[pfjetref];
3106 if(fabs(recoJet1.
eta())<1.3){
3108 if (recoJet1.
pt()>=20. && recoJet1.
pt()<=50.) {
3114 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3120 if (recoJet1.
pt()>140.) {
3126 }
else if(fabs(recoJet1.
eta())<3.0){
3127 if (recoJet1.
pt()>20. && recoJet1.
pt()<=50.) {
3133 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3139 if (recoJet1.
pt()>140.) {
3146 if (recoJet1.
pt()>20. && recoJet1.
pt()<=50.) {
3152 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3158 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
bool isPFMuon(HWW &, int index, bool requireSamePt=true, float dpt_max=1.0)
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_