Float_t mass = -99;
Double_t pt = esdtrack->Pt();
+ if(fptc==3)
+ pt = 2*pt;
+
if(TMath::Abs(ptot) < pmax && TMath::Abs(pullTPC) <= 3 && TMath::Abs(pt) < ptmax){
fhBBDeu->Fill(ptot*esdtrack->GetSign(),TPCSignal);
mass = ptot/TMath::Sqrt(gamma*gamma - 1); // using inner TPC mom. as approx.
fhTOF->Fill(ptot*esdtrack->GetSign(),beta);
- if(TMath::Abs(mass) > 2.7)continue;
- if(TMath::Abs(mass) < 1. )continue;
+ if(fptc==1){
+ if(TMath::Abs(mass) > 2.7)continue;
+ if(TMath::Abs(mass) < 1. )continue;
+ }
+ if(fptc==2){
+ if(TMath::Abs(mass) > 5.0)continue;
+ if(TMath::Abs(mass) < 1.8 )continue;
+ }
+ if(fptc==3){
+ if(TMath::Abs(mass) > 5.0)continue;
+ if(TMath::Abs(mass) < 1.8 )continue;
+ }
fhMassTOF->Fill(mass);
}