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