if (dedx < fMinimalTPCdedx) pass=kFALSE;
Double_t time[9];
track->GetIntegratedTimes(time);
+ if (fCutPID && (fParticleID!=AliPID::kUnknown)) //if kUnknown don't cut on PID
+ {
+ if (!PassesAODpidCut(track)) pass=kFALSE;
+ }
+
if (fQA) {
+ // changed 04062014 used to be filled before possible PID cut
Double_t momTPC = track->GetTPCmomentum();
QAbefore( 1)->Fill(momTPC,dedx);
QAbefore( 5)->Fill(track->Pt(),track->DCA());
if (pass) QAafter( 12)->Fill(track->P(),(track->GetTOFsignal()-time[AliPID::kProton]));
}
- if (fCutPID && (fParticleID!=AliPID::kUnknown)) //if kUnknown don't cut on PID
- {
- if (!PassesAODpidCut(track)) pass=kFALSE;
- }
-
+
return pass;
}
cuts->SetPhiMin( 0 );
cuts->SetPhiMax( TMath::TwoPi() );
cuts->SetApplyRecentering(kTRUE);
- return cuts;
+ cuts->SetVZEROgainEqualizationPerRing(kFALSE);
+ return cuts;
}
//-----------------------------------------------------------------------
AliFlowTrackCuts* AliFlowTrackCuts::GetStandardGlobalTrackCuts2010()
if (!PassesNucleiSelection(track)) pass=kFALSE;
break;
//end part added by Natasha
-
+
+ case kTPCTOFNsigma:
+ if (!PassesTPCTOFNsigmaCut(track)) pass = kFALSE;
+ break;
default:
printf("AliFlowTrackCuts::PassesCuts() this should never be called!\n");
pass=kFALSE;
return select;
}
// end part added by Natasha
+//-----------------------------------------------------------------------
+Bool_t AliFlowTrackCuts::PassesTPCTOFNsigmaCut(const AliAODTrack* track)
+{
+ // do a simple combined cut on the n sigma from tpc and tof
+ // with information of the pid response object (needs pid response task)
+ // stub, not implemented yet
+ if(!track) return kFALSE;
+ return kFALSE;
+
+}
+//-----------------------------------------------------------------------------
+Bool_t AliFlowTrackCuts::PassesTPCTOFNsigmaCut(const AliESDtrack* track)
+{
+ // do a simple combined cut on the n sigma from tpc and tof
+ // with information of the pid response object (needs pid response task)
+ // stub, not implemented yet
+ if(!track) return kFALSE;
+ return kFALSE;
+
+}
+
//-----------------------------------------------------------------------
void AliFlowTrackCuts::SetPriors(Float_t centrCur){
//set priors for the bayesian pid selection
return "TOFbetaSimple";
case kTPCNuclei:
return "TPCnuclei";
+ case kTPCTOFNsigma:
+ return "TPCTOFNsigma";
default:
return "NOPID";
}
kTPCdedx, // asymmetric cuts of TPC dedx signal
kTOFbetaSimple, //simple TOF only cut
kTPCbayesian, //bayesian cutTPC
- kTPCNuclei // added by Natasha for Nuclei
+ kTPCNuclei, // added by Natasha for Nuclei
+ kTPCTOFNsigma // simple cut on combined tpc tof nsigma
};
//setters (interface to AliESDtrackCuts)
// part added by F. Noferini
Bool_t PassesTOFbayesianCut(const AliESDtrack* track);
Bool_t PassesNucleiSelection(const AliESDtrack* track); // added by Natasha
+ Bool_t PassesTPCTOFNsigmaCut(const AliAODTrack* track);
+ Bool_t PassesTPCTOFNsigmaCut(const AliESDtrack* track);
Bool_t TPCTOFagree(const AliVTrack *track);
// end part added by F. Noferini