Float_t netConeArea = TMath::Pi()*(fIsoConeR*fIsoConeR - 0.04*0.04);
GetCeIso(clsVec, ceiso, cephiband, cecore);
GetTrIso(clsVec, triso, trphiband, trcore);
- ceiso -= Et;
alliso = ceiso + triso;
allphiband = cephiband + trphiband;
allcore = cecore + trcore;
Float_t trisoue = trphiband/phibandArea*netConeArea;
Float_t allisoue = allphiband/phibandArea*netConeArea;
fM02Et->Fill(Et, c->GetM02());
- if(ceiso<1)
+ if(ceiso-cecore<1)
fM02EtCeIso1->Fill(Et, c->GetM02());
- if(ceiso<2)
+ if(ceiso-cecore<2)
fM02EtCeIso2->Fill(Et, c->GetM02());
- if(ceiso<5)
+ if(ceiso-cecore<5)
fM02EtCeIso5->Fill(Et, c->GetM02());
- if(triso<1)
+ if(triso-trcore<1)
fM02EtTrIso1->Fill(Et, c->GetM02());
- if(triso<2)
+ if(triso-trcore<2)
fM02EtTrIso2->Fill(Et, c->GetM02());
- if(triso<5)
+ if(triso-trcore<5)
fM02EtTrIso5->Fill(Et, c->GetM02());
- if(alliso<1)
+ if(alliso-allcore<1)
fM02EtAllIso1->Fill(Et, c->GetM02());
- if(alliso<2)
+ if(alliso-allcore<2)
fM02EtAllIso2->Fill(Et, c->GetM02());
- if(alliso<5)
+ if(alliso-allcore<5)
fM02EtAllIso5->Fill(Et, c->GetM02());
if(c->GetM02()>0.1 && c->GetM02()<0.3){
fCeIsoVsEtPho->Fill(Et, ceiso - cecore - ceisoue);
Double_t dR = TMath::Sqrt(pow(c->GetTrackDx(),2)+pow(c->GetTrackDz(),2));
if(dR<0.014){
fM02EtTM->Fill(Et, c->GetM02());
- if(ceiso<1)
+ if(ceiso-cecore<1)
fM02EtCeIso1TM->Fill(Et, c->GetM02());
- if(ceiso<2)
+ if(ceiso-cecore<2)
fM02EtCeIso2TM->Fill(Et, c->GetM02());
- if(ceiso<5)
+ if(ceiso-cecore<5)
fM02EtCeIso5TM->Fill(Et, c->GetM02());
- if(triso<1)
+ if(triso-trcore<1)
fM02EtTrIso1TM->Fill(Et, c->GetM02());
- if(triso<2)
+ if(triso-trcore<2)
fM02EtTrIso2TM->Fill(Et, c->GetM02());
- if(triso<5)
+ if(triso-trcore<5)
fM02EtTrIso5TM->Fill(Et, c->GetM02());
- if(alliso<1)
+ if(alliso-allcore<1)
fM02EtAllIso1TM->Fill(Et, c->GetM02());
- if(alliso<2)
+ if(alliso-allcore<2)
fM02EtAllIso2TM->Fill(Et, c->GetM02());
- if(alliso<5)
+ if(alliso-allcore<5)
fM02EtAllIso5TM->Fill(Et, c->GetM02());
if(c->GetM02()>0.1 && c->GetM02()<0.3){
fCeIsoVsEtPhoTM->Fill(Et, ceiso);