]>
Commit | Line | Data |
---|---|---|
3feee083 | 1 | //__________________________________________________//\r |
667ca262 | 2 | AliBalance *GetBalanceFunctionObject(const char* analysisLevel = "AOD", \r |
6a3f819e | 3 | const char* centralityName = 0x0,\r |
74b04dd4 | 4 | Double_t centrMin = 0.,\r |
5 | Double_t centrMax = 100.,\r | |
49a3750f | 6 | Bool_t bShuffle = kFALSE,\r |
845fdeca | 7 | Bool_t bHBTcut = kFALSE,\r |
e1d6ee8e | 8 | Bool_t bConversionCut = kFALSE,\r |
9 | Bool_t kUsePID = kTRUE) {\r | |
10a747d2 | 10 | //Function to setup the AliBalance object and return it\r |
3feee083 | 11 | AliBalance *gBalance = new AliBalance();\r |
12 | gBalance->SetAnalysisLevel(analysisLevel);\r | |
6432ac6a | 13 | gBalance->SetShuffle(bShuffle);\r |
49a3750f | 14 | gBalance->SetHBTcut(bHBTcut);\r |
845fdeca | 15 | gBalance->SetConversionCut(bConversionCut);\r |
6a3f819e | 16 | if(centralityName) gBalance->SetCentralityIdentifier(centralityName);\r |
74b04dd4 | 17 | gBalance->SetCentralityInterval(centrMin,centrMax);\r |
3feee083 | 18 | \r |
3feee083 | 19 | //Set all analyses separately\r |
20 | //Rapidity\r | |
33a86282 | 21 | gBalance->SetInterval(AliBalance::kRapidity,-0.8,0.8,64,0.0,1.6); \r |
3feee083 | 22 | //Eta\r |
6432ac6a | 23 | gBalance->SetInterval(AliBalance::kEta,-0.8,0.8,64,0.0,1.6);\r |
3feee083 | 24 | //Qlong\r |
6432ac6a | 25 | gBalance->SetInterval(AliBalance::kQlong,-1,1,100,0.0,2.0);\r |
3feee083 | 26 | //Qout\r |
6432ac6a | 27 | gBalance->SetInterval(AliBalance::kQout,-1,1,100,0.0,2.0);\r |
3feee083 | 28 | //Qside\r |
6432ac6a | 29 | gBalance->SetInterval(AliBalance::kQside,-1,1,100,0.0,2.0);\r |
3feee083 | 30 | //Qinv\r |
e1d6ee8e | 31 | gBalance->SetInterval(AliBalance::kQinv,-1,1,100,0.0,5.0);\r |
3feee083 | 32 | //Phi\r |
33a86282 | 33 | gBalance->SetInterval(AliBalance::kPhi,0.,360.,90,0.,180.0);\r |
6432ac6a | 34 | \r |
35 | //Init the histograms\r | |
36 | gBalance->InitHistograms();\r | |
3feee083 | 37 | \r |
38 | return gBalance;\r | |
39 | }\r | |
40 | \r | |
41 | //__________________________________________________//\r | |
a839b0a3 | 42 | AliESDtrackCuts *GetTrackCutsObject(Double_t ptMin, Double_t ptMax, Double_t etaMin, Double_t etaMax, Double_t maxTPCchi2, Double_t maxDCAz, Double_t maxDCAxy, Int_t minNClustersTPC) {\r |
43 | \r | |
3feee083 | 44 | // only used for ESDs\r |
45 | // Function to setup the AliESDtrackCuts object and return it\r | |
46 | AliESDtrackCuts *cuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r | |
a839b0a3 | 47 | cuts->SetRequireTPCStandAlone(kTRUE); // TPC only cuts! \r |
48 | \r | |
49 | // extra TPC cuts (Syst studies)\r | |
50 | if(minNClustersTPC != -1) cuts->SetMinNClustersTPC(minNClustersTPC);\r | |
51 | else cuts->SetMinNClustersTPC(70); // standard for filter bit 128\r | |
3feee083 | 52 | \r |
a839b0a3 | 53 | if(maxTPCchi2 != -1) cuts->SetMaxChi2PerClusterTPC(maxTPCchi2);\r |
54 | \r | |
55 | // extra DCA cuts (Syst studies) \r | |
56 | if(maxDCAz!=-1 && maxDCAxy != -1){\r | |
57 | cuts->SetMaxDCAToVertexZ(maxDCAz);\r | |
58 | cuts->SetMaxDCAToVertexXY(maxDCAxy);\r | |
59 | }\r | |
60 | \r | |
61 | cuts->SetPtRange(ptMin,ptMax);\r | |
62 | cuts->SetEtaRange(etaMin,etaMax);\r | |
3feee083 | 63 | cuts->DefineHistograms(1);\r |
64 | //cuts->SaveHistograms("trackCuts");\r | |
65 | \r | |
66 | return cuts;\r | |
67 | }\r | |
68 | \r |