]>
Commit | Line | Data |
---|---|---|
a65a7e70 | 1 | AliAnalysisTaskElecHadronCorrel* ConfigHFEElecHadronCorrelPbPb(Bool_t useMC, |
2 | Bool_t EventTrigSelMB=kTRUE, | |
3 | Bool_t TrigSelCen = kTRUE, | |
4 | Double_t CentMin = 0, | |
5 | Double_t CentMax = 7, | |
6 | Double_t TPCNsigMinE = -2, | |
7 | Double_t TPCNsigMaxE = 2, | |
8 | Double_t TPCNsigMinH = -10, | |
9 | Double_t TPCNsigMaxH = -3.5, | |
10 | Double_t SSM02Min = 0.03, | |
11 | Double_t SSM02Max = 0.5, | |
12 | Double_t SSM20Min = 0.03, | |
13 | Double_t SSM20Max = 0.3, | |
14 | Double_t Disp = 1, | |
15 | Double_t EovPMin = 0.8, | |
16 | Double_t EovPMax = 1.2, | |
17 | Double_t InvM = 0.1, | |
18 | const char* CentralityMet = "VOM", | |
19 | Int_t TPCNClsEle = 100, | |
20 | Int_t TPCNClsHad = 80, | |
21 | Bool_t AssoEleITSref=kTRUE, | |
22 | Int_t AssoElecTPCNCls=80, | |
23 | TString ContNameExt = "Central", | |
24 | TString TaskName="hfeCorrl"){ | |
25 | ||
26 | // | |
27 | // HFE standard task configuration | |
28 | // | |
29 | ||
30 | Bool_t kAnalyseTaggedTracks = kTRUE; | |
31 | ||
32 | AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts"); | |
33 | // hfecuts->CreateStandardCuts(); | |
34 | hfecuts->SetMinNClustersTPC(TPCNClsEle); | |
35 | hfecuts->SetMinNClustersITS(3); | |
36 | hfecuts->SetMinNTrackletsTRD(0); | |
37 | hfecuts->SetMinRatioTPCclusters(0.6); | |
38 | // hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable); | |
39 | hfecuts->SetRequireITSPixel(); | |
40 | hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny); | |
41 | hfecuts->SetMaxChi2perClusterITS(-1); | |
42 | hfecuts->SetMaxChi2perClusterTPC(3.5); | |
43 | hfecuts->SetCheckITSLayerStatus(kFALSE); // shud be put back | |
44 | // hfecuts->UnsetVertexRequirement(); | |
45 | hfecuts->SetVertexRange(10.); | |
46 | hfecuts->SetRequireSigmaToVertex(); | |
47 | //hfecuts->SetSigmaToVertex(10); | |
48 | hfecuts->SetTOFPIDStep(kFALSE); | |
49 | // hfecuts->SetQAOn(); | |
50 | hfecuts->SetPtRange(0, 30); | |
51 | ||
52 | TString taskName = TaskName; | |
53 | ||
54 | AliAnalysisTaskElecHadronCorrel *task = new AliAnalysisTaskElecHadronCorrel(taskName); | |
55 | printf("task ------------------------ %p\n ", task); | |
56 | task->SetHFECuts(hfecuts); | |
57 | // task->SetRemovePileUp(kTRUE); | |
58 | // task->SetInvariantMassCut(0.1); | |
59 | ||
60 | task->SetEventTriggerSelectionMB(EventTrigSelMB); | |
61 | task->SetTriggerSelection(TrigSelCen); | |
62 | task->SetCentralityParameters(CentMin, CentMax, CentralityMet); | |
63 | task->SetInvariantMassCut(InvM); | |
64 | task->SetAssoElecTPCNCls(AssoElecTPCNCls); | |
65 | task->SetTPCnsigmaCutsElecSelection(TPCNsigMinE,TPCNsigMaxE); | |
66 | task->SetTPCnsigmaCutsHadSelection(TPCNsigMinH,TPCNsigMaxH); | |
67 | task->SetShowerShapeCutsM02(SSM02Min,SSM02Max); | |
68 | task->SetShowerShapeCutsM20(SSM20Min,SSM20Max); | |
69 | task->SetShowerShapeCutsDisp(0,Disp); | |
70 | task->SetEovPCuts(EovPMin,EovPMax); | |
71 | task->SetRejectKinkMother(kTRUE); | |
72 | task->SetTPCNClsHad(TPCNClsHad); | |
73 | task->SetAssoElecITSrefit(AssoEleITSref); | |
74 | ||
75 | // Define PID | |
76 | AliHFEpid *pid = task->GetPID(); | |
77 | if(useMC) pid->SetHasMCData(kTRUE); | |
78 | pid->AddDetector("TPC", 0); | |
79 | pid->AddDetector("EMCAL", 1); | |
80 | /* | |
81 | // change E/p cuts | |
82 | AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid); | |
83 | emcpid->SetEoPMax(1.2); | |
84 | emcpid->SetEoPMim(0.8); | |
85 | ||
86 | Double_t params[4]; | |
87 | char *cutmodel; | |
88 | if(useMC){ | |
89 | // Monte-Carlo needs modelling of the falling mean with momentum at low momentum | |
90 | // for high momentum it is consistent with a flat -0.94 | |
91 | cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x"; | |
92 | Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246}; | |
93 | for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar]; | |
94 | } else { | |
95 | // Data is consistent with a flat 0.12 | |
96 | cutmodel = "pol0"; | |
97 | //params[0] = -0.0015; | |
98 | //params[0] = -3.0; | |
99 | params[0] = -1.0; //sigma min | |
100 | } | |
101 | // pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); | |
102 | pid->ConfigureTPCasymmetric(0,30,-1,3.0); | |
103 | */ | |
104 | printf("*************************************\n"); | |
105 | printf("Configuring standard Task:\n"); | |
106 | // task->PrintStatus(); | |
107 | pid->PrintStatus(); | |
108 | printf("*************************************\n"); | |
109 | return task; | |
110 | } |