if((pi0->E()+pi0->Pz())/(pi0->E()-pi0->Pz())<=0){
if(hist)hist->Fill(cutIndex);
cutIndex++;
+// cout << "undefined rapidity" << endl;
return kFALSE;
}
else{
cutIndex++;
if(abs(pi0->Rapidity()-fRapidityShift)>fRapidityCutMeson){
if(hist)hist->Fill(cutIndex);
+// cout << abs(pi0->Rapidity()-fRapidityShift) << ">" << fRapidityCutMeson << endl;
return kFALSE;
}
}
// Opening Angle Cut
//fOpeningAngle=2*TMath::ATan(0.134/pi0->P());// physical minimum opening angle
if( pi0->GetOpeningAngle() < fOpeningAngle){
+// cout << pi0->GetOpeningAngle() << "<" << fOpeningAngle << endl;
if(hist)hist->Fill(cutIndex);
return kFALSE;
}
// Alpha Max Cut
if(pi0->GetAlpha()>fAlphaCutMeson){
+// cout << pi0->GetAlpha() << ">" << fAlphaCutMeson << endl;
if(hist)hist->Fill(cutIndex);
return kFALSE;
}
// Alpha Min Cut
if(pi0->GetAlpha()<fAlphaMinCutMeson){
+// cout << pi0->GetAlpha() << "<" << fAlphaMinCutMeson << endl;
if(hist)hist->Fill(cutIndex);
return kFALSE;
}
if (hDCARMesonPrimVtxBefore)hDCARMesonPrimVtxBefore->Fill(pi0->GetDCARMotherPrimVtx());
if (pi0->GetDCABetweenPhotons() > fDCAGammaGammaCut){
+// cout << pi0->GetDCABetweenPhotons() << ">" << fDCAGammaGammaCut << endl;
if(hist)hist->Fill(cutIndex);
return kFALSE;
}
cutIndex++;
if (pi0->GetDCARMotherPrimVtx() > fDCARMesonPrimVtxCut){
+// cout << pi0->GetDCARMotherPrimVtx() << ">" << fDCARMesonPrimVtxCut << endl;
if(hist)hist->Fill(cutIndex);
return kFALSE;
}
if (hDCAZMesonPrimVtxBefore)hDCAZMesonPrimVtxBefore->Fill(pi0->GetDCAZMotherPrimVtx());
if (abs(pi0->GetDCAZMotherPrimVtx()) > fDCAZMesonPrimVtxCut){
+// cout << pi0->GetDCAZMotherPrimVtx() << ">" << fDCAZMesonPrimVtxCut << endl;
if(hist)hist->Fill(cutIndex);
return kFALSE;
}