fh2MassVsNConstJet1Tagged[i] = 0;
fh2MassVsNConstJet1TaggedMatch[i] = 0;
- fh3PtJet1VsRatVsLeadPtAllSel[i] = 0;
- fh3PtJet1VsRatVsLeadPtTagged[i] = 0;
- fh3PtJet1VsRatVsLeadPtTaggedMatch[i] = 0;
- fpPtVsRatJet1All[i] = 0;
- fpPtVsRatJet1Tagged[i] = 0;
- fpPtVsRatJet1TaggedMatch[i] = 0;
- fh2RatVsAreaJet1All[i] = 0;
- fh2RatVsAreaJet1Tagged[i] = 0;
- fh2RatVsAreaJet1TaggedMatch[i] = 0;
- fh2RatVsNConstJet1All[i] = 0;
- fh2RatVsNConstJet1Tagged[i] = 0;
- fh2RatVsNConstJet1TaggedMatch[i] = 0;
+ fh3PtJet1VsRatVsLeadPtAllSel[i] = 0;
+ fh3PtJet1VsRatVsLeadPtTagged[i] = 0;
+ fh3PtJet1VsRatVsLeadPtTaggedMatch[i] = 0;
+ fpPtVsRatJet1All[i] = 0;
+ fpPtVsRatJet1Tagged[i] = 0;
+ fpPtVsRatJet1TaggedMatch[i] = 0;
+ fh2RatVsAreaJet1All[i] = 0;
+ fh2RatVsAreaJet1Tagged[i] = 0;
+ fh2RatVsAreaJet1TaggedMatch[i] = 0;
+ fh2RatVsNConstJet1All[i] = 0;
+ fh2RatVsNConstJet1Tagged[i] = 0;
+ fh2RatVsNConstJet1TaggedMatch[i] = 0;
fh2EtMassOverEtRSq[i] = 0;
}
if(jetL) drToLJ = jet1->DeltaR(jetL);
fh3MSubPtRawDRMatch[fCentBin]->Fill(var,ptjet1,drToLJ);
if(jet2) {
- Double_t var2 = -999.;
- if(jet2->Pt()>0. || jet2->Pt()<0.) var2 = jet2->M()/jet2->Pt();
+ Double_t var2 = jet2->M();
+ if(fJetMassVarType==kRatMPt) {
+ if(jet2->Pt()>0. || jet2->Pt()<0.) var2 = jet2->M()/jet2->Pt();
+ }
fh3MSubPtTrueDR[fCentBin]->Fill(var,jet2->Pt(),drToLJ);
fh3MTruePtTrueDR[fCentBin]->Fill(var2,jet2->Pt(),drToLJ);
fh3PtTrueDeltaMDR[fCentBin]->Fill(jet2->Pt(),var-var2,drToLJ);
if(jetL) drToLJ = jet1->DeltaR(jetL);
fh3MSubPtRawDRMatch[fCentBin]->Fill(var,ptjet1,drToLJ);
if(jet2) {
- Double_t var2 = -999.;
- if(jet2->Pt()>0. || jet2->Pt()<0.) var2 = jet2->M()/jet2->Pt();
+ Double_t var2 = jet2->M();
+ if(fJetMassVarType==kRatMPt) {
+ if(jet2->Pt()>0. || jet2->Pt()<0.) var2 = jet2->M()/jet2->Pt();
+ }
fh3MSubPtTrueDR[fCentBin]->Fill(var,jet2->Pt(),drToLJ);
fh3MTruePtTrueDR[fCentBin]->Fill(var2,jet2->Pt(),drToLJ);
fh3PtTrueDeltaMDR[fCentBin]->Fill(jet2->Pt(),var-var2,drToLJ);