Get the analysis.
1757 Int_t JetLoPass = 0;
1758 Int_t JetHiPass = 0;
1760 if (triggerResults.
isValid()){
1763 const unsigned int nTrig(triggerNames.
size());
1764 for (
unsigned int i=0;
i<nTrig;++
i)
1778 std::cout <<
">>> Trigger Lo = " << JetLoPass
1779 <<
" Hi = " << JetHiPass
1788 if (!vertexHandle.
isValid()) {
1789 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1790 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1793 if ( vertexHandle.
isValid() ){
1795 numPV = vertexCollection.size();
1805 if (!gtReadoutRecord.
isValid()) {
1806 LogInfo(
"JetAnalyzer") <<
"JetAnalyzer: Could not find GT readout record" << std::endl;
1807 if (
verbose_)
std::cout <<
"JetAnalyzer: Could not find GT readout record product" << std::endl;
1811 bool techTriggerResultBxF =
false;
1812 bool techTriggerResultBx0 =
false;
1814 if (!gtReadoutRecord.
isValid()) {
1815 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find GT readout record" << std::endl;
1816 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find GT readout record product" << std::endl;
1820 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBxF = gtReadoutRecord->technicalTriggerWord(-1);
1821 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBx0 = gtReadoutRecord->technicalTriggerWord();
1826 if(techTriggerResultBx0!=0){
1855 bool pass_Z_selection=
false;
1869 for (
unsigned int i=0;
i<Muons->size();
i++){
1870 bool pass_muon_id=
false;
1871 bool pass_muon_iso=
false;
1872 double dxy=fabs((*Muons)[
i].muonBestTrack()->dxy());
1873 double dz=fabs((*Muons)[
i].muonBestTrack()->dz());
1875 dxy=fabs((*Muons)[
i].muonBestTrack()->dxy((*vertexHandle)[0].
position()));
1876 dz=fabs((*Muons)[
i].muonBestTrack()->dz((*vertexHandle)[0].
position()));
1878 if((*Muons)[
i].
pt()>20 && fabs((*Muons)[
i].
eta())<2.3){
1879 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){
1883 float muonIsoPFdb = ((*Muons)[
i].pfIsolationR04().sumChargedHadronPt
1884 +
std::max(0., (*Muons)[
i].pfIsolationR04().sumNeutralHadronEt + (*Muons)[
i].pfIsolationR04().sumPhotonEt - 0.5*(*Muons)[
i].pfIsolationR04().sumPUPt))/(*Muons)[
i].pt();
1885 if(muonIsoPFdb<0.12){
1889 if(pass_muon_id && pass_muon_iso){
1890 if((*Muons)[
i].
pt()>pt0){
1891 mu_index1=mu_index0;
1894 pt0=(*Muons)[
i].pt();
1895 }
else if ((*Muons)[
i].
pt()>pt1){
1897 pt1=(*Muons)[
i].pt();
1902 if(mu_index0>=0 && mu_index1>=0){
1903 if((*Muons)[mu_index0].
charge()*(*Muons)[mu_index1].
charge()<0){
1904 zCand=(*Muons)[mu_index0].polarP4()+(*Muons)[mu_index1].polarP4();
1905 if(fabs(zCand.M()-91.)<20 && zCand.Pt()>30){
1906 pass_Z_selection=
true;
1923 bool jetCollectionIsValid =
false;
1930 if (
jetCleaningFlag_ && (!jetCollectionIsValid || !bPrimaryVertex || !dcsDecision))
return;
1932 unsigned int collSize=-1;
1942 if (
isPFJet_) collSize=pfJets->size();
1947 std::vector<Jet> corJets;
1958 bool cleaned_first_jet=
false;
1959 bool cleaned_second_jet=
false;
1962 int ind1_mu_vetoed=-1;
1963 double pt1_mu_vetoed=-1;
1964 int ind2_mu_vetoed=-1;
1965 double pt2_mu_vetoed=-1;
1966 bool cleaned_first_jet_mu_vetoed=
false;
1967 bool cleaned_second_jet_mu_vetoed=
false;
1976 bool pass_correction_flag=
false;
1980 pass_correction_flag=
true;
1984 pass_correction_flag=
true;
1987 for (
unsigned int ijet=0; ijet<collSize; ijet++) {
1990 bool pass_uncorrected=
false;
1991 bool pass_corrected=
false;
1993 correctedJet=(*caloJets)[ijet];
1999 correctedJet=(*pfJets)[ijet];
2002 correctedJet=(*patJets)[ijet];
2005 pass_uncorrected=
true;
2008 pass_uncorrected=
true;
2012 scale = jetCorr->correction((*caloJets)[ijet]);
2015 scale = jetCorr->correction((*pfJets)[ijet]);
2020 pass_corrected=
true;
2026 recoJets.push_back(correctedJet);
2028 bool jetpassid=
true;
2029 bool Thiscleaned=
true;
2035 reco::JetID jetID = (*jetID_ValueMap_Handle)[calojetref];
2039 Thiscleaned=jetpassid;
2041 if(Thiscleaned && pass_corrected){
2049 Thiscleaned=jetpassid;
2051 if(Thiscleaned && pass_uncorrected){
2058 if(Thiscleaned && pass_corrected){
2063 if (fabs((*caloJets)[ijet].
eta()) <= 1.3) {
2066 }
else if(fabs((*caloJets)[ijet].
eta()) <3.0){
2087 Thiscleaned = jetpassid;
2088 JetIDWPU = jetpassid;
2090 if(Thiscleaned && pass_uncorrected){
2096 if(Thiscleaned && pass_corrected){
2099 if(fabs(correctedJet.
eta()) <= 1.3) {
2100 if(correctedJet.
pt()<=50.){
2103 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2106 if(correctedJet.
pt()>140.){
2113 }
else if(fabs(correctedJet.
eta()) <= 3) {
2114 if(correctedJet.
pt()<=50.){
2117 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2120 if(correctedJet.
pt()>140.){
2127 }
else if(fabs(correctedJet.
eta()) <= 5) {
2128 if(correctedJet.
pt()<=50.){
2131 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2134 if(correctedJet.
pt()>140.){
2147 int puidmvaflag=-10;
2148 int puidcutflag=-10;
2149 puidmva=(*puJetIdMva)[pfjetref];
2150 puidcut=(*puJetId)[pfjetref];
2151 puidmvaflag=(*puJetIdFlagMva)[pfjetref];
2152 puidcutflag=(*puJetIdFlag)[pfjetref];
2155 Thiscleaned = jetpassid;
2158 if(Thiscleaned && pass_uncorrected){
2164 if(Thiscleaned && pass_corrected){
2207 if (correctedJet.
pt()<= 50) {
2212 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2217 if (correctedJet.
pt()>140.) {
2222 if (fabs(correctedJet.
eta()) <= 1.3) {
2224 if (correctedJet.
pt()<=50.) {
2244 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2264 if (correctedJet.
pt()>140.) {
2287 }
else if(fabs(correctedJet.
eta()) <= 3) {
2289 if (correctedJet.
pt()<=50.) {
2309 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2329 if (correctedJet.
pt()>140.) {
2356 if (correctedJet.
pt()<=50.) {
2370 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2384 if (correctedJet.
pt()>140.) {
2428 if(techTriggerResultBx0 && techTriggerResultBxF){
2430 if(fabs(correctedJet.
eta()) <= 1.3) {
2435 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2440 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2445 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2449 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2482 if(techTriggerResultBx0 && !techTriggerResultBxF){
2484 if(fabs(correctedJet.
eta()) <= 1.3) {
2489 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2494 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2499 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2503 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2532 if(pass_Z_selection){
2533 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 ){
2534 if(correctedJet.
pt()>pt1_mu_vetoed){
2535 pt2_mu_vetoed=pt1_mu_vetoed;
2536 ind2_mu_vetoed=ind1_mu_vetoed;
2537 cleaned_second_jet_mu_vetoed=cleaned_first_jet_mu_vetoed;
2538 pt1_mu_vetoed=correctedJet.
pt();
2539 ind1_mu_vetoed=ijet;
2540 cleaned_first_jet_mu_vetoed=JetIDWPU;
2541 }
else if(correctedJet.
pt()>pt2_mu_vetoed){
2542 pt2_mu_vetoed=correctedJet.
pt();
2543 ind2_mu_vetoed=ijet;
2544 cleaned_second_jet_mu_vetoed=JetIDWPU;
2549 if(correctedJet.
pt()>pt1){
2555 cleaned_second_jet=cleaned_first_jet;
2556 pt1=correctedJet.
pt();
2558 cleaned_first_jet=JetIDWPU;
2559 }
else if(correctedJet.
pt()>pt2){
2563 pt2=correctedJet.
pt();
2565 cleaned_second_jet=JetIDWPU;
2566 }
else if(correctedJet.
pt()>pt3){
2567 pt3=correctedJet.
pt();
2571 if(!pass_corrected){
2580 if(fabs(correctedJet.
eta())<3.0){
2586 if(fabs(correctedJet.
eta())<3.0){
2592 if(Thiscleaned && pass_corrected){
2619 if((*patJets)[ijet].hasTagInfo(
"caTop")){
2621 if ( tagInfo != 0 ) {
2627 if((*patJets)[ijet].hasSubjets(
"CMSTopTag")){
2630 if((*patJets)[ijet].hasSubjets(
"CMSTopTag") && (*patJets)[ijet].subjets(
"CMSTopTag").size()>0){
2635 if((*patJets)[ijet].subjets(
"CMSTopTag").size()>1){
2640 if((*patJets)[ijet].subjets(
"CMSTopTag").size()>2){
2645 if((*patJets)[ijet].subjets(
"CMSTopTag").size()>3){
2654 if((*patJets)[ijet].hasSubjets(
"SoftDrop")){
2657 if((*patJets)[ijet].hasSubjets(
"SoftDrop") && (*patJets)[ijet].subjets(
"SoftDrop").size()>0){
2662 if((*patJets)[ijet].subjets(
"SoftDrop").size()>1){
2677 if((*patJets)[ijet].hasTagInfo(
"caTop")){
2679 if ( tagInfo_boosted != 0 ) {
2685 if((*patJets)[ijet].hasSubjets(
"CMSTopTag")){
2688 if((*patJets)[ijet].hasSubjets(
"CMSTopTag") && (*patJets)[ijet].subjets(
"CMSTopTag").size()>0){
2693 if((*patJets)[ijet].subjets(
"CMSTopTag").size()>1){
2698 if((*patJets)[ijet].subjets(
"CMSTopTag").size()>2){
2703 if((*patJets)[ijet].subjets(
"CMSTopTag").size()>3){
2712 if((*patJets)[ijet].hasSubjets(
"SoftDrop")){
2715 if((*patJets)[ijet].hasSubjets(
"SoftDrop") && (*patJets)[ijet].subjets(
"SoftDrop").size()>0){
2720 if((*patJets)[ijet].subjets(
"SoftDrop").size()>1){
2761 if (fabs(correctedJet.
eta()) <= 1.3) {
2768 }
else if (fabs(correctedJet.
eta()) <= 3) {
2796 if(cleaned_first_jet){
2801 if(numofjets>1 && cleaned_second_jet) {
2802 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2804 dphi=2*acos(-1.)-dphi;
2817 if(
jetCleaningFlag_ && recoJets.size()>1 && cleaned_first_jet && cleaned_second_jet ){
2824 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2826 dphi=2*acos(-1.)-dphi;
2834 reco::JetID jetID1 = (*jetID_ValueMap_Handle)[calojetref1];
2836 reco::JetID jetID2 = (*jetID_ValueMap_Handle)[calojetref2];
2934 float QGLikelihood=-10;
2939 QGmulti=(*qgMultiplicity)[pfjetref];
2940 QGLikelihood=(*qgLikelihood)[pfjetref];
2941 QGptD=(*qgptD)[pfjetref];
2942 QGaxis2=(*qgaxis2)[pfjetref];
2943 if(fabs(recoJets[0].
eta())<1.3){
2945 if (recoJets[0].
pt()>=20. && recoJets[0].
pt()<=50.) {
2951 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2957 if (recoJets[0].
pt()>140.) {
2963 }
else if(fabs(recoJets[0].
eta())<3.0){
2964 if (recoJets[0].
pt()>20. && recoJets[0].
pt()<=50.) {
2970 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2976 if (recoJets[0].
pt()>140.) {
2983 if (recoJets[0].
pt()>20. && recoJets[0].
pt()<=50.) {
2989 if (recoJets[0].
pt()>50. && recoJets[0].
pt()<=140.) {
2995 if (recoJets[0].
pt()>140.) {
3003 QGmulti=(*qgMultiplicity)[pfjetref1];
3004 QGLikelihood=(*qgLikelihood)[pfjetref1];
3005 QGptD=(*qgptD)[pfjetref1];
3006 QGaxis2=(*qgaxis2)[pfjetref1];
3007 if(fabs(recoJets[1].
eta())<1.3){
3009 if (recoJets[1].
pt()>=20. && recoJets[1].
pt()<=50.) {
3015 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
3021 if (recoJets[1].
pt()>140.) {
3027 }
else if(fabs(recoJets[1].
eta())<3.0){
3028 if (recoJets[1].
pt()>20. && recoJets[1].
pt()<=50.) {
3034 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
3040 if (recoJets[1].
pt()>140.) {
3047 if (recoJets[1].
pt()>20. && recoJets[1].
pt()<=50.) {
3053 if (recoJets[1].
pt()>50. && recoJets[1].
pt()<=140.) {
3059 if (recoJets[1].
pt()>140.) {
3151 if (fabs(recoJets[0].
eta() < 1.4)) {
3152 double pt_dijet = (recoJets[0].pt() + recoJets[1].pt())/2;
3158 int randJet =iEvent.
id().
event()%2;
3159 if (fabs(recoJets[1].
eta() < 1.4)) {
3175 bool thirdJetCut =
true;
3188 double dijetAsymmetry =(recoJets[jet1].pt() - recoJets[jet2].pt()) / (recoJets[jet1].
pt() + recoJets[jet2].pt());
3198 pt_barrel = recoJets[jet1].pt();
3199 pt_probe = recoJets[jet2].pt();
3209 bool thirdJetCut =
true;
3212 thirdJetCut =
false;
3216 double dijetBalance = (pt_probe - pt_barrel) / pt_dijet;
3225 if(pass_Z_selection && ind1_mu_vetoed>=0 && pt1_mu_vetoed>12 && cleaned_first_jet_mu_vetoed &&
isPFJet_){
3226 bool pass_second_jet_mu_vetoed=
false;
3227 if(cleaned_second_jet_mu_vetoed){
3228 if(ind2_mu_vetoed>=0 && pt2_mu_vetoed/zCand.Pt()<0.2){
3229 pass_second_jet_mu_vetoed=
true;
3232 if(pass_second_jet_mu_vetoed){
3235 recoJet1=(*pfJets)[ind1_mu_vetoed];
3240 scale = jetCorr->correction((*caloJets)[ind1_mu_vetoed]);
3243 scale = jetCorr->correction((*pfJets)[ind1_mu_vetoed]);
3247 double dphi=fabs(recoJet1.
phi()-zCand.Phi());
3249 dphi=2*acos(-1.)-dphi;
3255 if(fabs(dphi-acos(-1.))<0.34){
3263 if(!calometcoll.
isValid())
return;
3264 met=&(calometcoll->front());
3268 if(!pfmetcoll.
isValid())
return;
3269 met=&(pfmetcoll->front());
3278 if(recoJet1.
pt()>20){
3303 double MPF=1.+(met->
px()*zCand.Px()+met->
py()*zCand.Py())/(zCand.Pt()*zCand.Pt());
3304 if(fabs(recoJet1.
eta())<1.3){
3312 }
else if (zCand.Pt()<140){
3324 }
else if (zCand.Pt()<75){
3326 }
else if (zCand.Pt()<150){
3328 }
else if (zCand.Pt()<290){
3334 }
else if(fabs(recoJet1.
eta())<3.0){
3341 }
else if (zCand.Pt()<140){
3353 }
else if (zCand.Pt()<75){
3355 }
else if (zCand.Pt()<150){
3357 }
else if (zCand.Pt()<290){
3370 }
else if (zCand.Pt()<140){
3382 }
else if (zCand.Pt()<100){
3390 float QGLikelihood=-10;
3395 QGmulti=(*qgMultiplicity)[pfjetref];
3396 QGLikelihood=(*qgLikelihood)[pfjetref];
3397 QGptD=(*qgptD)[pfjetref];
3398 QGaxis2=(*qgaxis2)[pfjetref];
3399 if(fabs(recoJet1.
eta())<1.3){
3401 if (recoJet1.
pt()>=20. && recoJet1.
pt()<=50.) {
3407 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3413 if (recoJet1.
pt()>140.) {
3419 }
else if(fabs(recoJet1.
eta())<3.0){
3420 if (recoJet1.
pt()>20. && recoJet1.
pt()<=50.) {
3426 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3432 if (recoJet1.
pt()>140.) {
3439 if (recoJet1.
pt()>20. && recoJet1.
pt()<=50.) {
3445 if (recoJet1.
pt()>50. && recoJet1.
pt()<=140.) {
3451 if (recoJet1.
pt()>140.) {
MonitorElement * mSubJet1_CMSTopTag_eta_boosted
MonitorElement * mqgLikelihood_lowPt_Barrel
edm::EDGetTokenT< edm::ValueMap< float > > qgaxis2Token_
MonitorElement * mCHEn_mediumPt_Barrel
MonitorElement * mSubJet3_CMSTopTag_mass
MonitorElement * mEFrac_EndCap
MonitorElement * mNHEn_lowPt_Barrel
MonitorElement * mJetZBalance_lowZPt_J_Barrel
MonitorElement * mSubJet2_CMSTopTag_eta_boosted
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 * mSubJet2_CMSTopTag_phi_boosted
MonitorElement * mDeltaPt_Z_j1_over_ZPt_30_55_J_Barrel
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
MonitorElement * mSubJet1_SoftDrop_phi_boosted
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
MonitorElement * mtau2_over_tau1
static bool jetSortingRule(reco::Jet x, reco::Jet y)
const properties_type & properties() const
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 * mtau2_over_tau1_boosted
MonitorElement * mePhFracEndCapPlus_BXm1Filled
MonitorElement * mJetZBalance_highZPt_J_Barrel
MonitorElement * mSubJet2_SoftDrop_phi_boosted
MonitorElement * mEmEnergyInHF
MonitorElement * mtau3_over_tau2
MonitorElement * mPhiFirst
MonitorElement * mCutPUJIDDiscriminant_lowPt_Forward
MonitorElement * mCATopTag_nSubJets
MonitorElement * mMPF_mediumZPt_J_Forward
MonitorElement * mNeutralHadronEnergy
MonitorElement * mHFEFrac_mediumPt_Forward
MonitorElement * mPt_Forward
MonitorElement * mCHFrac_highPt_EndCap
MonitorElement * mNHFrac_highPt_EndCap
MonitorElement * mSubJet4_CMSTopTag_phi
MonitorElement * mSubJet3_CMSTopTag_pt
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_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * mHFrac_Forward
MonitorElement * mEmEnergyInEE
MonitorElement * mMass_Barrel
MonitorElement * mMass_Forward
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 * mElEn_highPt_Barrel
edm::EDGetTokenT< reco::CaloJetCollection > caloJetsToken_
MonitorElement * mHFHFrac
MonitorElement * mPrunedMass_boosted
MonitorElement * mPhEn_lowPt_EndCap
Base class for all types of Jets.
MonitorElement * mSubJet1_SoftDrop_mass
double balanceThirdJetCut_
MonitorElement * mMediumCutPUJIDPassFractionVSeta
MonitorElement * mHFEEn_lowPt_Forward
edm::EDGetTokenT< edm::ValueMap< float > > qgLikelihoodToken_
MonitorElement * mSubJet4_CMSTopTag_mass
MonitorElement * meNHFracEndCapMinus_BXm1Filled
MonitorElement * mPhi_EndCap
JetIDSelectionFunctor jetIDFunctor
MonitorElement * mPrunedMass
MonitorElement * mMPF_J_EndCap
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 * mSubJet1_CMSTopTag_mass
MonitorElement * mN90Hits
MonitorElement * mSubJet1_SoftDrop_eta_boosted
PFJetIDSelectionFunctor pfjetIDFunctor
virtual double phi() const final
momentum azimuthal angle
MonitorElement * mSubJet1_SoftDrop_pt_boosted
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 * mPhEn_highPt_Barrel
MonitorElement * mCHFrac_profile
MonitorElement * mSubJet2_SoftDrop_eta_boosted
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 * mCATopTag_minMass
MonitorElement * mMass_lowPt_Barrel
MonitorElement * mEFrac_profile
MonitorElement * mePhFracBarrel_BXm1Filled
MonitorElement * mDeltaPt_Z_j1_over_ZPt_150_290_J_EndCap
MonitorElement * mSubJet2_SoftDrop_pt_boosted
MonitorElement * mHFEMFrac
MonitorElement * mConstituents_Forward
MonitorElement * mPt_CaloJet
edm::EDGetTokenT< std::vector< reco::Vertex > > vertexToken_
MonitorElement * mMVAPUJIDDiscriminant_lowPt_Forward
MonitorElement * mpTD_lowPt_Forward
MonitorElement * mtau3_over_tau2_boosted
MonitorElement * mnSubJetsSoftDrop
JetMETDQMDCSFilter * DCSFilterForDCSMonitoring_
MonitorElement * mePhFracEndCapPlus_BXm1Empty
MonitorElement * mMPF_highZPt_J_Barrel
MonitorElement * mMultiplicityQG_lowPt_EndCap
MonitorElement * meHFHFracPlus_BXm1Filled
MonitorElement * meNHFracEndCapPlus_BXm1Filled
MonitorElement * mSubJet2_CMSTopTag_mass
MonitorElement * mCHFrac_lowPt_EndCap
MonitorElement * mMuMultiplicity_highPt_Barrel
MonitorElement * mNHFracVSeta_mediumPt
MonitorElement * mChargedMultiplicity_profile
MonitorElement * mAxis2_mediumPt_Barrel
MonitorElement * mHFEFrac_lowPt_Forward
MonitorElement * mMass_EndCap
MonitorElement * meHFEMFracPlus_BXm1Filled
MonitorElement * mNHFracVSeta_lowPt
MonitorElement * mLooseJIDPassFractionVSpt
MonitorElement * mAxis2_highPt_Forward
JetMETDQMDCSFilter * DCSFilterForJetMonitoring_
MonitorElement * mpTD_mediumPt_Barrel
MonitorElement * mPt_uncor
MonitorElement * mDeltaPt_Z_j1_over_ZPt_75_150_J_Barrel
MonitorElement * mEMF_CaloJet
MonitorElement * mChargedMuEnergy
MonitorElement * mLooseMVAPUJIDPassFractionVSpt
MonitorElement * mMuMultiplicity_mediumPt_EndCap
MonitorElement * mNHFrac_profile
MonitorElement * mConstituents_profile
MonitorElement * mpTD_highPt_EndCap
MonitorElement * mEmEnergyInEB
MonitorElement * mDijetBalance
MonitorElement * mDeltaPt_Z_j1_over_ZPt_55_75_J_EndCap
MonitorElement * mSubJet1_CMSTopTag_mass_boosted
MonitorElement * mCATopTag_topMass_boosted
MonitorElement * mJ1Pt_over_ZPt_J_Barrel
MonitorElement * mJ1Pt_over_ZPt_J_EndCap
static bool passJetId(int flag, Id level)
MonitorElement * mTightCutPUJIDPassFractionVSeta
MonitorElement * mMediumMVAPUJIDPassFractionVSeta
MonitorElement * mSubJet3_CMSTopTag_eta_boosted
MonitorElement * mePtEndCapPlus_BXm1Empty
MonitorElement * mSubJet1_SoftDrop_mass_boosted
MonitorElement * mChargedHadronEnergy_profile
MonitorElement * mPhi_Barrel
MonitorElement * mPhFrac_lowPt_Barrel
MonitorElement * mSubJet2_CMSTopTag_pt_boosted
MonitorElement * mSubJet2_SoftDrop_eta
MonitorElement * mAxis2_highPt_Barrel
MonitorElement * mMass_lowPt_EndCap
edm::EDGetTokenT< reco::CaloMETCollection > caloMetToken_
MonitorElement * mConstituents_EndCap
MonitorElement * mSubJet1_CMSTopTag_phi
MonitorElement * meEta_BXm1Empty
MonitorElement * mHFHFrac_mediumPt_Forward
MonitorElement * mPhFrac_mediumPt_Barrel
MonitorElement * mFilteredMass
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > gtToken_
MonitorElement * mJetZBalance_mediumZPt_J_Forward
MonitorElement * mSubJet1_CMSTopTag_pt
MonitorElement * mSoftDropMass
MonitorElement * mePtEndCapMinus_BXm1Empty
MonitorElement * mNeutralMultiplicity_profile
MonitorElement * mMultiplicityQG_highPt_EndCap
MonitorElement * meHFHFracMinus_BXm1Filled
MonitorElement * cleanupME
MonitorElement * mPhEn_mediumPt_Barrel
MonitorElement * mTrimmedMass_boosted
MonitorElement * mSubJet2_SoftDrop_mass_boosted
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_
MonitorElement * mSubJet2_SoftDrop_pt
edm::EDGetTokenT< pat::JetCollection > patJetsToken_
MonitorElement * mConstituents_Barrel
MonitorElement * mCATopTag_minMass_boosted
MonitorElement * mJetZBalance_lowZPt_J_EndCap
MonitorElement * meHFHFracMinus_BXm1Empty
MonitorElement * meCHFracEndCapMinus_BXm1Empty
MonitorElement * mNHEn_lowPt_EndCap
MonitorElement * mSubJet4_CMSTopTag_pt_boosted
MonitorElement * mMuEn_lowPt_EndCap
MonitorElement * meHFEMFracMinus_BXm1Filled
MonitorElement * mNeutMultiplicity_lowPt_Barrel
MonitorElement * mSubJet3_CMSTopTag_phi_boosted
MonitorElement * mpTD_highPt_Barrel
MonitorElement * mPhFrac_mediumPt_EndCap
MonitorElement * mSubJet4_CMSTopTag_phi_boosted
MonitorElement * mJ1Pt_over_ZPt_highZPt_J_Forward
MonitorElement * mMuonMultiplicity_profile
MonitorElement * mSubJet2_CMSTopTag_pt
MonitorElement * mCutPUJIDDiscriminant_mediumPt_Forward
MonitorElement * mCutPUJIDDiscriminant_highPt_EndCap
MonitorElement * mChargedHadronEnergy
MonitorElement * mPhi_uncor
MonitorElement * mSubJet2_CMSTopTag_eta
MonitorElement * mCutPUJIDDiscriminant_highPt_Barrel
MonitorElement * mChargedEmEnergy
MonitorElement * mNeutMultiplicity_mediumPt_Forward
MonitorElement * mSubJet1_CMSTopTag_eta
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
virtual double py() const final
y coordinate of momentum vector
MonitorElement * mPt_profile
edm::EDGetTokenT< edm::ValueMap< float > > qgptDToken_
MonitorElement * mPhFracVSeta_lowPt
MonitorElement * mChMultiplicity_lowPt_Barrel
edm::InputTag jetCorrectorTag_
MonitorElement * mnSubJetsCMSTopTag_boosted
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 * mCATopTag_topMass
MonitorElement * mDeltaPt_Z_j1_over_ZPt_75_150_J_EndCap
MonitorElement * mSubJet3_CMSTopTag_eta
MonitorElement * mChMultiplicity_lowPt_EndCap
MonitorElement * mMuEn_mediumPt_Barrel
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken_
MonitorElement * mSubJet3_CMSTopTag_phi
MonitorElement * mMuMultiplicity_highPt_EndCap
MonitorElement * mAxis2_highPt_EndCap
MonitorElement * mMVAPUJIDDiscriminant_highPt_Barrel
MonitorElement * mSubJet1_SoftDrop_pt
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 * 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 * mCutPUJIDDiscriminant_mediumPt_Barrel
MonitorElement * mPhFrac_highPt_EndCap
MonitorElement * mPhFracVSeta_mediumPt
MonitorElement * mMPF_lowZPt_J_Barrel
MonitorElement * meHFEMFracPlus_BXm1Empty
MonitorElement * mSubJet3_CMSTopTag_pt_boosted
double deltaR(double eta1, double eta2, double phi1, double phi2)
MonitorElement * mnSubJetsSoftDrop_boosted
MonitorElement * meCHFracBarrel_BXm1Empty
MonitorElement * mePtBarrel_BXm1Empty
MonitorElement * mNHFracVSpT_Barrel
MonitorElement * mHadEnergyInHF
edm::EDGetTokenT< reco::MuonCollection > MuonsToken_
MonitorElement * mSubJet1_SoftDrop_eta
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 * mSubJet2_CMSTopTag_phi
MonitorElement * mAxis2_lowPt_EndCap
MonitorElement * mElEn_mediumPt_Barrel
MonitorElement * mMVAPUJIDDiscriminant_highPt_Forward
TObject * getRootObject(void) const
MonitorElement * mSubJet3_CMSTopTag_mass_boosted
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
MonitorElement * mNHFrac_highPt_Barrel
MonitorElement * mqgLikelihood_highPt_Forward
MonitorElement * mMediumCutPUJIDPassFractionVSpt
MonitorElement * mePhFracEndCapMinus_BXm1Empty
MonitorElement * mSubJet4_CMSTopTag_eta_boosted
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 * 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 * mSubJet2_SoftDrop_phi
MonitorElement * mMuMultiplicity_lowPt_Barrel
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 * mSubJet4_CMSTopTag_mass_boosted
MonitorElement * mMultiplicityQG_lowPt_Barrel
MonitorElement * mMVAPUJIDDiscriminant_lowPt_EndCap
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
MonitorElement * mTrimmedMass
static int position[264][3]
MonitorElement * mePhFracEndCapMinus_BXm1Filled
MonitorElement * mHFHEn_highPt_Forward
MonitorElement * mePtEndCapPlus_BXm1Filled
MonitorElement * mNHEn_highPt_EndCap
MonitorElement * mePtForwardPlus_BXm1Filled
MonitorElement * mSubJet1_CMSTopTag_phi_boosted
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
virtual double px() const final
x coordinate of momentum vector
MonitorElement * mElEn_lowPt_Barrel
MonitorElement * mMass_highPt_Barrel
MonitorElement * mHadEnergyInHE
std::vector< std::string > highPtJetExpr_
MonitorElement * mMass_mediumPt_Barrel
MonitorElement * meHFEMFracMinus_BXm1Empty
MonitorElement * mCATopTag_nSubJets_boosted
MonitorElement * mChMultiplicity_highPt_EndCap
MonitorElement * mConstituents
MonitorElement * mNHFrac_lowPt_Barrel
MonitorElement * mSubJet4_CMSTopTag_eta
MonitorElement * meEta_BXm1Filled
MonitorElement * mJetEnergyCorrVSEta
MonitorElement * mPhFrac_highPt_Barrel
virtual double eta() const final
momentum pseudorapidity
MonitorElement * mCHFracVSpT_Barrel
MonitorElement * mPhFracVSpT_EndCap
MonitorElement * mChMultiplicity_mediumPt_Barrel
MonitorElement * mMPF_J_Barrel
MonitorElement * mJetZBalance_highZPt_J_EndCap
MonitorElement * mMultiplicityQG_lowPt_Forward
MonitorElement * mEta_profile
MonitorElement * mPt_EndCap
MonitorElement * mCHFracVSeta_lowPt
MonitorElement * mpTD_mediumPt_Forward
MonitorElement * mMVAPUJIDDiscriminant_mediumPt_EndCap
MonitorElement * mCHEn_mediumPt_EndCap
MonitorElement * mTightCutPUJIDPassFractionVSpt
MonitorElement * mSubJet2_SoftDrop_mass
MonitorElement * mSubJet4_CMSTopTag_pt
MonitorElement * mNeutMultiplicity_highPt_Forward
MonitorElement * mDPhiZJet
MonitorElement * mqgLikelihood_highPt_Barrel
MonitorElement * mqgLikelihood_mediumPt_Forward
MonitorElement * mePtForwardMinus_BXm1Filled
MonitorElement * mLooseCutPUJIDPassFractionVSpt
MonitorElement * mMuonMultiplicity
MonitorElement * mPhFracVSpT_Barrel
MonitorElement * mSubJet1_CMSTopTag_pt_boosted
MonitorElement * mHFEMFrac_profile
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
MonitorElement * mnSubJetsCMSTopTag
MonitorElement * mSubJet2_CMSTopTag_mass_boosted
MonitorElement * mConstituents_uncor
MonitorElement * mNeutMultiplicity_mediumPt_Barrel
MonitorElement * meCHFracEndCapMinus_BXm1Filled
MonitorElement * mCutPUJIDDiscriminant_highPt_Forward
MonitorElement * mNeutMultiplicity_lowPt_Forward
MonitorElement * mHFEEn_mediumPt_Forward
MonitorElement * mSubJet1_SoftDrop_phi
MonitorElement * mNeutMultiplicity_highPt_Barrel
MonitorElement * verticesME
MonitorElement * mPhFracVSeta_highPt
MonitorElement * mAxis2_mediumPt_EndCap
MonitorElement * mFilteredMass_boosted
virtual double pt() const final
transverse momentum
MonitorElement * mpTD_highPt_Forward
MonitorElement * mePtForwardMinus_BXm1Empty
MonitorElement * mHFHFrac_profile
MonitorElement * mJ1Pt_over_ZPt_mediumZPt_J_Barrel
MonitorElement * mSoftDropMass_boosted
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
