//________________________________________________________________
AliAnalysisCuts* SetupTrackCuts(Int_t cutInstance)
{
std::cout << "SetupTrackCuts()" <<std::endl;
//AliAnalysisCuts* trackCuts=0x0;
AliESDtrackCuts *fesdTrackCuts = new AliESDtrackCuts();
//global
fesdTrackCuts->SetPtRange( 0.2 , 100. );
fesdTrackCuts->SetEtaRange( -0.8 , 0.8 );
fesdTrackCuts->SetAcceptKinkDaughters(kFALSE);
fesdTrackCuts->SetRequireSigmaToVertex(kFALSE);
fesdTrackCuts->SetDCAToVertex2D(kFALSE);
fesdTrackCuts->SetMaxDCAToVertexZ(3.);
fesdTrackCuts->SetMaxDCAToVertexXY(1.);
fesdTrackCuts->SetRequireTPCRefit(kTRUE);
fesdTrackCuts->SetRequireITSRefit(kTRUE);
if(cutInstance == 0){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 1){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(130);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 2){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 3){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(130);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 4){
fesdTrackCuts->SetMinNClustersITS(6);
fesdTrackCuts->SetMaxChi2PerClusterITS(2.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 5){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 6){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 7){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 8){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 9){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 10){
fesdTrackCuts->SetMinNClustersITS(5);
// fesdTrackCuts->SetMaxChi2PerClusterITS(6);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(120);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 11){
fesdTrackCuts->SetMinNClustersITS(5);
// fesdTrackCuts->SetMaxChi2PerClusterITS(100);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(120);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 12){
fesdTrackCuts->SetMinNClustersITS(6);
// fesdTrackCuts->SetMaxChi2PerClusterITS(100);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 13){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 14){
fesdTrackCuts->SetMinNClustersITS(6);
fesdTrackCuts->SetMaxChi2PerClusterITS(2.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 15){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 16){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
// fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 17){
fesdTrackCuts->SetMinNClustersITS(3);
// fesdTrackCuts->SetMaxChi2PerClusterITS(4);
// fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 18){
fesdTrackCuts->SetMinNClustersITS(4);
// fesdTrackCuts->SetMaxChi2PerClusterITS(6);
// fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 19){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
return fesdTrackCuts;
}
AliESDtrackCuts *CreateTrackCutsPWGJE(Int_t cutMode) {
//
// Macro to create track cuts for PWG Jet analysis
// User can select a specific set by indicating cutMode
// cutMode has 8 digits: first 4 digits additional cuts, last 4 digits standard cuts
// additional cuts are variations of standard cuts (used for hybrid track selection and QA)
// Numbering starts from 1000 For standard and additional cut numbers
AliESDtrackCuts *trackCuts = new AliESDtrackCuts("AliESDtrackCuts");
TString tag;
Int_t mod = 10000;
Bool_t bStdCutsDefined = kFALSE;
//_____________________________________________________________________
// STANDARD CUTS
//Get standard cuts: last 4 digits of cutMode
Int_t stdCutMode = cutMode%mod;
if(stdCutMode == 1000) {
bStdCutsDefined = kTRUE;
trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE,1);
trackCuts->SetMinNCrossedRowsTPC(120);
trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);
trackCuts->SetMaxChi2PerClusterITS(36);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 1e10);
tag = "Global track RAA analysis QM2011 + Chi2ITS<36";
}
if(stdCutMode == 1001) {
bStdCutsDefined = kTRUE;
// TPC
trackCuts->SetMinNClustersTPC(90);
trackCuts->SetMaxChi2PerClusterTPC(4);
trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
trackCuts->SetAcceptKinkDaughters(kFALSE);
trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
// ITS
trackCuts->SetRequireITSRefit(kTRUE);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
//reject fakes
trackCuts->SetMaxChi2PerClusterITS(36);
trackCuts->SetRequireSigmaToVertex(kFALSE);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 100.);
tag = "Global tracks jet analysis with ITSrefit and NclsIter1=90, noSPD requirement";
}
if(stdCutMode == 1002) {
bStdCutsDefined = kTRUE;
// TPC
trackCuts->SetMinNClustersTPC(80);
trackCuts->SetMaxChi2PerClusterTPC(4);
trackCuts->SetAcceptKinkDaughters(kFALSE);
trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
// ITS
trackCuts->SetRequireITSRefit(kTRUE);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
//reject fakes
trackCuts->SetMaxChi2PerClusterITS(36);
trackCuts->SetRequireSigmaToVertex(kFALSE);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 100.);
tag = "Global tracks jet analysis with ITSrefit and Ncls=80, noSPD requirement";
}
if(stdCutMode == 1003) {
bStdCutsDefined = kTRUE;
// tight global tracks
trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kFALSE,1);
trackCuts->SetMinNClustersTPC(0);
trackCuts->SetMinNCrossedRowsTPC(120);
trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.1);// essentially switches it off
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
trackCuts->SetMaxChi2PerClusterITS(36);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
tag = "Global tracks ITSTPC2010 + NCrossedRows + loose ITS";
}
if(stdCutMode == 1004) {
bStdCutsDefined = kTRUE;
// TPC
trackCuts->SetMinNClustersTPC(70);
trackCuts->SetMaxChi2PerClusterTPC(4);
trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
trackCuts->SetAcceptKinkDaughters(kFALSE);
trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
// ITS
trackCuts->SetRequireITSRefit(kTRUE);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
//reject fakes
trackCuts->SetMaxChi2PerClusterITS(36);
trackCuts->SetRequireSigmaToVertex(kFALSE);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 100.);
tag = "Global tracks jet analysis with ITSrefit and NclsIter1=70, noSPD requirement";
}
if(stdCutMode == 1005) {
bStdCutsDefined = kTRUE;
// TPC
trackCuts->SetMinNClustersTPC(70);
trackCuts->SetMaxChi2PerClusterTPC(4);
trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
trackCuts->SetAcceptKinkDaughters(kFALSE);
trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
// ITS
trackCuts->SetRequireITSRefit(kTRUE);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
//reject fakes
trackCuts->SetMaxChi2PerClusterITS(36);
trackCuts->SetRequireSigmaToVertex(kFALSE);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 1E+15);
tag = "Global tracks jet analysis with ITSrefit and NclsIter1=70, noSPD requirement, no upper pt cut";
}
if(stdCutMode == 1006) {
bStdCutsDefined = kTRUE;
// TPC
TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
trackCuts->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep,20.);
trackCuts->SetMinNClustersTPC(70);
trackCuts->SetMaxChi2PerClusterTPC(4);
trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
trackCuts->SetAcceptKinkDaughters(kFALSE);
trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
// ITS
trackCuts->SetRequireITSRefit(kTRUE);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
//reject fakes
trackCuts->SetMaxChi2PerClusterITS(36);
trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
trackCuts->SetRequireSigmaToVertex(kFALSE);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 1E+15);
tag = "Global tracks jet analysis with ITSrefit and NclsIter1=PtDep, noSPD requirement, no upper pt cut, golden chi2";
}
if(stdCutMode == 1007) {
bStdCutsDefined = kTRUE;
trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE,1);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
//
trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 1E+15);
tag = "Global tracks with AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE)";
}
if(stdCutMode == 1008) {
bStdCutsDefined = kTRUE;
trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE,1);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
//
trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
tag = "Global tracks 2011 with NCrossedRows cut";
}
if(stdCutMode == 2000) {
bStdCutsDefined = kTRUE;
// TPC
trackCuts->SetMinNClustersTPC(90);
trackCuts->SetMaxChi2PerClusterTPC(4);
trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
trackCuts->SetAcceptKinkDaughters(kFALSE);
trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMaxFractionSharedTPCClusters(0.4);
//accept secondaries
trackCuts->SetMaxDCAToVertexXY(2.4);
trackCuts->SetMaxDCAToVertexZ(3.2);
trackCuts->SetDCAToVertex2D(kTRUE);
trackCuts->SetRequireSigmaToVertex(kFALSE);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 100.);
tag = "Global tracks jet analysis, loose cuts, NClsIter1=90, no ITS requirements";
}
if(stdCutMode == 2001) {
bStdCutsDefined = kTRUE;
trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
// trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMinNClustersTPC(70);
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 100.);
tag = "TPConly track cuts, loose cuts, NCls=70, no ITS requirements";
}
if(stdCutMode == 2002) {
bStdCutsDefined = kTRUE;
trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
// trackCuts->SetRequireTPCRefit(kTRUE);
trackCuts->SetMinNClustersTPC(0);
trackCuts->SetMinNCrossedRowsTPC(120);
trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.1);// essentially switches it off
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 100.);
tag = "TPConly track cuts, loose cuts, NCls=70, no ITS requirements";
}
if(!bStdCutsDefined) {
printf("last 4 digits do not represent a predefined set of standard cuts. Returning 0\n");
return 0;
}
//_____________________________________________________________________
// ADDITIONAL CUTS
//Get additional cut mode: first 4 digits of cutMode
Int_t addCutMode = (int)((float)cutMode/(float)mod);
if(addCutMode == 1000) {
trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
tag += " + additonal: SPD any requirement";
}
if(addCutMode == 1001) {
trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);
tag += " + additional: w/o hits in SPD";
}
if(addCutMode == 1002) {
trackCuts->SetMaxChi2PerClusterITS(1E10);
tag += " + additional: maxITSChi2=1e10";
}
if(addCutMode == 1003) {
trackCuts->SetMinNClustersTPC(0);
trackCuts->SetMinNCrossedRowsTPC(0);
trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.);
tag += " + additional: minClusters=0 minCrossedRows=0 minCrossedRowsOverFindable=0";
}
if(addCutMode == 1004) {
trackCuts->SetRequireITSRefit(kFALSE);
tag += " + additional: ITSrefit=kFALSE";
}
if(addCutMode == 1005) {
trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kOff);
tag += " + additional: no SPD requirement (kOff)";
}
Printf("Created track cuts for: %s", tag.Data());
return trackCuts;
}
AliESDtrackCuts* MyTrackCuts(Int_t clusterCut = 1, Bool_t selPrimaries = kTRUE, TString optSyt="DefaultITSTPC2011")
{
Double_t dcaxymax = 2.0;
TString PtDcaFormula = "0.0105+0.0350/pt^1.1";//7sigma
Double_t dcazmax = 2.0;
Double_t minNcls = 50;
Double_t maxX2TPCcls = 4.0;
Double_t maxX2ITScls = 36.0;
Double_t minCrossedRows = 70.0;
Double_t minRatioClsCrRowsOverFCls = 0.8;
if(optSyt.Contains("PtDCAXY5s")) {PtDcaFormula = "0.0075+0.025/pt^1.1";}
if(optSyt.Contains("PtDCAXY6s")) {PtDcaFormula = "0.0090+0.030/pt^1.1";}
if(optSyt.Contains("PtDCAXY7s")) {PtDcaFormula = "0.0105+0.035/pt^1.1";}//Defult
if(optSyt.Contains("PtDCAXY8s")) {PtDcaFormula = "0.0120+0.040/pt^1.1";}
if(optSyt.Contains("PtDCAXY9s")) {PtDcaFormula = "0.0135+0.045/pt^1.1";}
if(optSyt.Contains("FixDCAZ1")) {dcazmax = 1.0;}
if(optSyt.Contains("FixDCAZ2")) {dcazmax = 2.0;}//Defult
if(optSyt.Contains("FixDCAZ3")) {dcazmax = 3.0;}
if(optSyt.Contains("NCrRows60")){minCrossedRows = 60;}
if(optSyt.Contains("NCrRows70")){minCrossedRows = 70;}//Defult
if(optSyt.Contains("NCrRows80")){minCrossedRows = 80;}
if(optSyt.Contains("NCrRows90")){minCrossedRows = 90;}
if(optSyt.Contains("RClsCrRowsOvFCls0.7")){minRatioClsCrRowsOverFCls = 0.7;}
if(optSyt.Contains("RClsCrRowsOvFCls0.8")){minRatioClsCrRowsOverFCls = 0.8;}//Defult
if(optSyt.Contains("RClsCrRowsOvFCls0.9")){minRatioClsCrRowsOverFCls = 0.9;}
if(optSyt.Contains("ChiSqrPerTPCCls3")) {maxX2TPCcls = 3.0;}
if(optSyt.Contains("ChiSqrPerTPCCls4")) {maxX2TPCcls = 4.0;}//Defult
if(optSyt.Contains("ChiSqrPerTPCCls5")) {maxX2TPCcls = 5.0;}
if(optSyt.Contains("ChiSqrPerITSCls30")) {maxX2ITScls = 30.0;}
if(optSyt.Contains("ChiSqrPerITSCls36")) {maxX2ITScls = 36.0;}//Defult
if(optSyt.Contains("ChiSqrPerITSCls45")) {maxX2ITScls = 45.0;}
AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts;
if(optSyt.Contains("DefaultITSTPC2011")){
::Info("Config KSTAR ", Form("Default 2011 ESD track cuts used : %s\n",optSyt.Data()));
if(clusterCut == 0) esdTrackCuts->SetMinNClustersTPC(50);
else if (clusterCut == 1) {
esdTrackCuts->SetMinNCrossedRowsTPC(70);
esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);
}
esdTrackCuts->SetMaxChi2PerClusterTPC(4);
esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
esdTrackCuts->SetRequireTPCRefit(kTRUE);
// ITS
esdTrackCuts->SetRequireITSRefit(kTRUE);
esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
AliESDtrackCuts::kAny);
if(selPrimaries) {
// 7*(0.0015+0.0050/pt^1.1)
esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0105+0.0350/pt^1.1");
esdTrackCuts->SetMaxChi2TPCConstrainedGlobal(36);
}
esdTrackCuts->SetMaxDCAToVertexZ(2);
esdTrackCuts->SetDCAToVertex2D(kFALSE);
esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
esdTrackCuts->SetMaxChi2PerClusterITS(36);
}
else{
::Info("Config KSTAR ", Form("User Defined ESD track cuts used for Sys : %s\n",optSyt.Data()));
if(clusterCut == 0) esdTrackCuts->SetMinNClustersTPC(minNcls);
else if (clusterCut == 1) {
esdTrackCuts->SetMinNCrossedRowsTPC(minCrossedRows);
esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(minRatioClsCrRowsOverFCls);
}
esdTrackCuts->SetMaxChi2PerClusterTPC(maxX2TPCcls);
esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
esdTrackCuts->SetRequireTPCRefit(kTRUE);
// ITS
esdTrackCuts->SetRequireITSRefit(kTRUE);
esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
AliESDtrackCuts::kAny);
if(selPrimaries) {
// 7*(0.0015+0.0050/pt^1.1)
esdTrackCuts->SetMaxDCAToVertexXYPtDep(PtDcaFormula.Data());
esdTrackCuts->SetMaxChi2TPCConstrainedGlobal(36);
}
esdTrackCuts->SetMaxDCAToVertexZ(dcazmax);
esdTrackCuts->SetDCAToVertex2D(kFALSE);
esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
esdTrackCuts->SetMaxChi2PerClusterITS(maxX2ITScls);
}
return esdTrackCuts;
}