+ Double_t M02u;
+ if(Et<12)
+ M02u = 0.02486*Et*Et - 0.7289*Et + 6.266;
+ else
+ M02u = 14.32/Et - 0.09863;
+ if(M02u<0.65)
+ M02u = 0.65;
+ Double_t M02l = 12.88/Et - 0.3194;
+ if(M02l<0.4)
+ M02l = 0.4;
+ if(dR>0.028){
+ fM02Et->Fill(Et, c->GetM02());
+ if(ceiso-cecore<1)
+ fM02EtCeIso1->Fill(Et, c->GetM02());
+ if(ceiso-cecore<2)
+ fM02EtCeIso2->Fill(Et, c->GetM02());
+ if(ceiso-cecore<5)
+ fM02EtCeIso5->Fill(Et, c->GetM02());
+ if(triso-trcore<1)
+ fM02EtTrIso1->Fill(Et, c->GetM02());
+ if(triso-trcore<2)
+ fM02EtTrIso2->Fill(Et, c->GetM02());
+ if(triso-trcore<5)
+ fM02EtTrIso5->Fill(Et, c->GetM02());
+ if(alliso-allcore<1)
+ fM02EtAllIso1->Fill(Et, c->GetM02());
+ if(alliso-allcore<2)
+ fM02EtAllIso2->Fill(Et, c->GetM02());
+ if(alliso-allcore<5)
+ fM02EtAllIso5->Fill(Et, c->GetM02());
+ if(c->GetM02()>0.1 && c->GetM02()<0.3){
+ fCeIsoVsEtPho->Fill(Et, ceiso - cecore - ceisoue);
+ fTrIsoVsEtPho->Fill(Et, triso - trcore - trisoue);
+ fAllIsoVsEtPho->Fill(Et, alliso - allcore - allisoue);
+ }
+ if(c->GetM02()>M02l && c->GetM02()<M02u){
+ fCeIsoVsEtPi0->Fill(Et, ceiso - cecore - ceisoue);
+ fTrIsoVsEtPi0->Fill(Et, triso - trcore - trisoue);
+ fAllIsoVsEtPi0->Fill(Et, alliso - allcore - allisoue);
+ }
+ }
+ else
+ {