]>
Commit | Line | Data |
---|---|---|
86eaf9ae | 1 | void AddTask_GammaConvDalitzV1_pPb( Int_t trainConfig = 1, |
2 | Bool_t isMC = kFALSE, //run MC | |
39e8c2ce | 3 | Bool_t enableQAMesonTask = kTRUE, //enable QA in AliAnalysisTaskGammaConvDalitzV1 |
4 | Bool_t enableDoMesonChic = kFALSE, // enable additional Chic analysis | |
f315320f | 5 | TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input |
6 | Bool_t doWeighting = kFALSE, //enable Weighting | |
7 | TString generatorName = "DPMJET", | |
39e8c2ce | 8 | TString cutnumberAODBranch = "0000000060084001001500000" |
9 | ) { | |
10 | ||
11 | ||
12 | ||
209b710e | 13 | cout<<"*********Parameters*******"<<endl; |
14 | cout<<"trainConfig: "<<trainConfig<<endl; | |
15 | cout<<"isMC: "<<isMC<<endl; | |
16 | cout<<"enableQAMesonTask: "<<enableQAMesonTask<<endl; | |
17 | cout<<"enableDoMesonChic: "<<enableDoMesonChic<<endl; | |
18 | cout<<"fileNameInputForWeighting: "<<fileNameInputForWeighting.Data()<<endl; | |
19 | cout<<"doWeighting: "<<doWeighting<<endl; | |
20 | cout<<"generatorName: "<<generatorName.Data()<<endl; | |
21 | cout<<"cutnumberAODBranch: "<<cutnumberAODBranch.Data()<<endl; | |
39e8c2ce | 22 | |
23 | // ================= Load Librariers ================================= | |
24 | gSystem->Load("libCore.so"); | |
25 | gSystem->Load("libTree.so"); | |
26 | gSystem->Load("libGeom.so"); | |
27 | gSystem->Load("libVMC.so"); | |
28 | gSystem->Load("libPhysics.so"); | |
29 | gSystem->Load("libMinuit"); | |
30 | gSystem->Load("libSTEERBase"); | |
31 | gSystem->Load("libESD"); | |
32 | gSystem->Load("libAOD"); | |
33 | gSystem->Load("libANALYSIS"); | |
34 | gSystem->Load("libANALYSISalice"); | |
35 | gSystem->Load("libPWGGAGammaConv.so"); | |
36 | gSystem->Load("libCDB.so"); | |
37 | gSystem->Load("libSTEER.so"); | |
38 | gSystem->Load("libSTEERBase.so"); | |
39 | gSystem->Load("libTENDER.so"); | |
40 | gSystem->Load("libTENDERSupplies.so"); | |
41 | ||
42 | ||
43 | cout<<"Entro 0"<<endl; | |
44 | ||
45 | // ================== GetAnalysisManager =============================== | |
46 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); | |
47 | if (!mgr) { | |
86eaf9ae | 48 | Error(Form("AddTask_GammaConvDalitzV1_pPb_%i",trainConfig), "No analysis manager found."); |
39e8c2ce | 49 | return ; |
50 | } | |
51 | ||
52 | // ================== GetInputEventHandler ============================= | |
53 | AliVEventHandler *inputHandler=mgr->GetInputEventHandler(); | |
54 | ||
55 | //========= Add PID Reponse to ANALYSIS manager ==== | |
56 | if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){ | |
57 | gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C"); | |
58 | AddTaskPIDResponse(isMC); | |
59 | } | |
60 | ||
61 | //========= Set Cutnumber for V0Reader ================================ | |
ae4f2cfb | 62 | TString ConvCutnumber = "8000000060084001001500000000"; //Online V0 finder |
86eaf9ae | 63 | TString ElecCuts = "9000540000000200000"; //Electron Cuts |
5a104c5e | 64 | Bool_t doEtaShift = kFALSE; |
39e8c2ce | 65 | |
66 | ||
67 | ||
68 | AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer(); | |
69 | ||
70 | //========= Add V0 Reader to ANALYSIS manager if not yet existent ===== | |
71 | if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){ | |
72 | AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1"); | |
73 | ||
74 | fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE); | |
75 | fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output | |
76 | fV0ReaderV1->SetUseAODConversionPhoton(kTRUE); | |
77 | ||
78 | if (!mgr) { | |
79 | Error("AddTask_V0ReaderV1", "No analysis manager found."); | |
80 | return; | |
81 | } | |
82 | ||
83 | // Set AnalysisCut Number | |
84 | AliConversionCuts *fCuts=NULL; | |
85 | if( ConvCutnumber !=""){ | |
86 | fCuts= new AliConversionCuts(ConvCutnumber.Data(),ConvCutnumber.Data()); | |
87 | fCuts->SetPreSelectionCutFlag(kTRUE); | |
88 | if(fCuts->InitializeCutsFromCutString(ConvCutnumber.Data())){ | |
5a104c5e | 89 | fCuts->DoEtaShift(doEtaShift); |
90 | fV0ReaderV1->SetConversionCuts(fCuts); | |
39e8c2ce | 91 | fCuts->SetFillCutHistograms("",kTRUE); |
92 | } | |
93 | } | |
94 | if(inputHandler->IsA()==AliAODInputHandler::Class()){ | |
95 | // AOD mode | |
96 | fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data())); | |
97 | } | |
98 | fV0ReaderV1->Init(); | |
99 | ||
100 | AliLog::SetGlobalLogLevel(AliLog::kInfo); | |
101 | ||
102 | //connect input V0Reader | |
103 | mgr->AddTask(fV0ReaderV1); | |
104 | mgr->ConnectInput(fV0ReaderV1,0,cinput); | |
105 | ||
106 | } | |
107 | ||
108 | //================================================ | |
109 | //========= Add Electron Selector ================ | |
110 | ||
111 | ||
112 | if( !(AliDalitzElectronSelector*)mgr->GetTask("ElectronSelector") ){ | |
113 | ||
114 | AliDalitzElectronSelector *fElectronSelector = new AliDalitzElectronSelector("ElectronSelector"); | |
115 | ||
116 | // Set AnalysisCut Number | |
117 | ||
118 | AliDalitzElectronCuts *fElecCuts=0; | |
119 | ||
5a104c5e | 120 | //ElecCuts = "900054000000020000"; |
39e8c2ce | 121 | |
122 | if( ElecCuts!=""){ | |
123 | ||
124 | fElecCuts= new AliDalitzElectronCuts(ElecCuts.Data(),ElecCuts.Data()); | |
125 | ||
126 | if(fElecCuts->InitializeCutsFromCutString(ElecCuts.Data())){ | |
127 | ||
128 | fElectronSelector->SetDalitzElectronCuts(fElecCuts); | |
129 | ||
130 | fElecCuts->SetFillCutHistograms("",kTRUE); | |
131 | ||
132 | } | |
133 | ||
134 | } | |
135 | ||
136 | fElectronSelector->Init(); | |
137 | mgr->AddTask(fElectronSelector); | |
5a104c5e | 138 | |
139 | AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer(); | |
140 | ||
141 | //connect input V0Reader | |
142 | ||
143 | mgr->ConnectInput (fElectronSelector,0,cinput1); | |
39e8c2ce | 144 | |
145 | } | |
146 | ||
147 | ||
148 | ||
149 | cout<<"Entro"<<endl; | |
150 | //================================================ | |
151 | //========= Add task to the ANALYSIS manager ===== | |
152 | //================================================ | |
153 | // find input container | |
154 | ||
155 | ||
156 | ||
157 | AliAnalysisTaskGammaConvDalitzV1 *task=NULL; | |
158 | ||
86eaf9ae | 159 | task= new AliAnalysisTaskGammaConvDalitzV1(Form("GammaConvDalitzV1_%i",trainConfig)); |
39e8c2ce | 160 | |
161 | task->SetIsHeavyIon(2); | |
162 | task->SetIsMC(isMC); | |
163 | ||
164 | ||
165 | ||
166 | // Cut Numbers to use in Analysis | |
899d014b | 167 | Int_t numberOfCuts = 4; |
39e8c2ce | 168 | |
169 | TString *ConvCutarray = new TString[numberOfCuts]; | |
170 | ||
171 | TString *ElecCutarray = new TString[numberOfCuts]; | |
172 | ||
173 | TString *MesonCutarray = new TString[numberOfCuts]; | |
174 | ||
175 | Bool_t doEtaShiftIndCuts = kFALSE; | |
39e8c2ce | 176 | TString stringShift = ""; |
177 | ||
178 | // Shifting in pPb direction | |
179 | ||
899d014b | 180 | doEtaShiftIndCuts = kFALSE; |
39e8c2ce | 181 | stringShift = "pPb"; |
182 | ||
86eaf9ae | 183 | |
2279d237 | 184 | |
2279d237 | 185 | |
899d014b | 186 | if( trainConfig == 1 ) { // No eta shift |Y| < 0.8 |
cdd7cb8b | 187 | |
899d014b | 188 | ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
189 | ConvCutarray[1] = "8000011002093603007900000000"; ElecCutarray[1] = "9047540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 15 | |
190 | ConvCutarray[2] = "8000011002093603007800000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 20 | |
191 | ConvCutarray[3] = "8000011002093603007100000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 50 | |
cdd7cb8b | 192 | |
899d014b | 193 | } else if( trainConfig == 2 ) { // No eta shift |Y| < 0.8 |
194 | ||
195 | ConvCutarray[0] = "8000011002093603002200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.7 | |
196 | ConvCutarray[1] = "8000011002093603003200000000"; ElecCutarray[1] = "9047540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.5 | |
197 | ConvCutarray[2] = "8000011002093653007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 0.3 GeV Low and 3.5 High momentum | |
198 | ConvCutarray[3] = "8000011002093601007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 0.5 GeV Low and 5.0 High momentum | |
cdd7cb8b | 199 | |
899d014b | 200 | } else if( trainConfig == 3 ) { // No eta shift |Y| < 0.8 |
201 | ||
202 | ConvCutarray[0] = "8000011002093803007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 2.0sigmas Low and 1 High momentum | |
203 | ConvCutarray[1] = "8000011002093603007200000000"; ElecCutarray[1] = "9043540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and 0 High momentum | |
204 | ConvCutarray[2] = "8000011002093603007200000000"; ElecCutarray[2] = "9047740023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 0.3 GeV Low and 3.5 High momentum | |
205 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047520023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 0.5 GeV Low and 5.0 High momentum | |
cdd7cb8b | 206 | |
899d014b | 207 | } else if( trainConfig == 4 ) { // No eta shift |Y| < 0.8 |
208 | ||
209 | ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9042540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and -1 High momentum | |
210 | ConvCutarray[1] = "8000011002093603007200000000"; ElecCutarray[1] = "9047540013310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + SPD first layer | |
211 | ConvCutarray[2] = "8000011002093603007200000000"; ElecCutarray[2] = "9047540023330262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + PsiPair cut 0.52 | |
212 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031005009000"; //standard cut Pi0 pPb 00-100 Standard cut + Alpha cut < 0.7 | |
cdd7cb8b | 213 | |
899d014b | 214 | } else if( trainConfig == 5 ) { // No eta shift |Y| < 0.8 |
215 | ||
216 | ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9037540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -5,5 | |
217 | ConvCutarray[1] = "8000011002093603007200000000"; ElecCutarray[1] = "9057540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -3,5 | |
218 | ConvCutarray[2] = "8000011002091603007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -5,5 | |
219 | ConvCutarray[3] = "8000011002092603007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -3,5 | |
cdd7cb8b | 220 | |
899d014b | 221 | } else if ( trainConfig == 6 ) { //No eta shift |Y| < 0.8 |
cdd7cb8b | 222 | |
899d014b | 223 | ConvCutarray[0] = "8000011042093603007200000000"; ElecCutarray[0] = "9047540023510262371"; MesonCutarray[0] = "01032035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Y < 0.70 and prim and sec e |eta| < 0.75 //NOTE revisar |
224 | ConvCutarray[1] = "8000011002093603007200000000"; ElecCutarray[1] = "9047540023310263371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single prim Pt cut > 0.150 | |
225 | ConvCutarray[2] = "8000011002093603007200000000"; ElecCutarray[2] = "9047540023310263171"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + Single prim Pt cut > 0.100 | |
226 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540023310262271"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + DCAxy < 1 cm | |
cdd7cb8b | 227 | |
228 | ||
899d014b | 229 | } else if ( trainConfig == 7 ) { // No eta shift |Y| < 0.8 |
230 | ||
231 | ConvCutarray[0] = "8000011002493603007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.075 | |
232 | ConvCutarray[1] = "8000011002193603007200000000"; ElecCutarray[1] = "9047540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.100 | |
233 | ConvCutarray[2] = "8000011002083603007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls sec > 0.35 | |
234 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540027310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls prim > 0.60 | |
cdd7cb8b | 235 | |
cdd7cb8b | 236 | |
899d014b | 237 | } else if ( trainConfig == 8 ) { //No eta shift |Y| < 0.8 |
cdd7cb8b | 238 | |
899d014b | 239 | ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9047540023310262381"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.015 < InvMass(e+,e-) < 0.050 |
240 | ConvCutarray[1] = "8000011002093603007200000000"; ElecCutarray[1] = "9047540023310262391"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.025 < InvMass(e+,e-) < 0.035 | |
241 | ConvCutarray[2] = "8000011002093603001200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + qT < 0.1 | |
242 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 | |
cdd7cb8b | 243 | |
899d014b | 244 | } else if ( trainConfig == 9 ) { //No eta shift |Y| < 0.8 |
245 | ||
cdd7cb8b | 246 | |
899d014b | 247 | ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9047540023310272371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Stardad cut +100 events background |
248 | ConvCutarray[1] = "8000011002093603007200000000"; ElecCutarray[1] = "9047540023310162371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Background method V0 multiplicity | |
249 | ConvCutarray[2] = "8000011002093603007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035000000"; //standard cut Pi0 PbPb 00-100 + No extra smearing | |
250 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540025310262171"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 + Old Standard | |
cdd7cb8b | 251 | |
899d014b | 252 | } else if( trainConfig == 10 ) { // No eta shift |Y| < 0.8 + AddedSignals |
cdd7cb8b | 253 | |
899d014b | 254 | ConvCutarray[0] = "8000012002093603007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
255 | ConvCutarray[1] = "8000012002093603007900000000"; ElecCutarray[1] = "9047540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 15 | |
256 | ConvCutarray[2] = "8000012002093603007800000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 20 | |
257 | ConvCutarray[3] = "8000012002093603007100000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 50 | |
cdd7cb8b | 258 | |
899d014b | 259 | } else if( trainConfig == 11 ) { // No eta shift |Y| < 0.8 + AddedSignals |
260 | ||
261 | ConvCutarray[0] = "8000012002093603002200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.7 | |
262 | ConvCutarray[1] = "8000012002093603003200000000"; ElecCutarray[1] = "9047540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.5 | |
263 | ConvCutarray[2] = "8000012002093653007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 0.3 GeV Low and 3.5 High momentum | |
264 | ConvCutarray[3] = "8000012002093601007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 0.5 GeV Low and 5.0 High momentum | |
cdd7cb8b | 265 | |
899d014b | 266 | } else if( trainConfig == 12 ) { // No eta shift |Y| < 0.8 + AddedSignals |
267 | ||
268 | ConvCutarray[0] = "8000012002093803007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 2.0sigmas Low and 1 High momentum | |
269 | ConvCutarray[1] = "8000012002093603007200000000"; ElecCutarray[1] = "9043540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and 0 High momentum | |
270 | ConvCutarray[2] = "8000012002093603007200000000"; ElecCutarray[2] = "9047740023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 0.3 GeV Low and 3.5 High momentum | |
271 | ConvCutarray[3] = "8000012002093603007200000000"; ElecCutarray[3] = "9047520023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 0.5 GeV Low and 5.0 High momentum | |
cdd7cb8b | 272 | |
899d014b | 273 | } else if( trainConfig == 13 ) { // No eta shift |Y| < 0.8 + AddedSignals |
274 | ||
275 | ConvCutarray[0] = "8000012002093603007200000000"; ElecCutarray[0] = "9042540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and -1 High momentum | |
276 | ConvCutarray[1] = "8000012002093603007200000000"; ElecCutarray[1] = "9047540013310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + SPD first layer | |
277 | ConvCutarray[2] = "8000012002093603007200000000"; ElecCutarray[2] = "9047540023330262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + PsiPair cut 0.52 | |
278 | ConvCutarray[3] = "8000012002093603007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031005009000"; //standard cut Pi0 pPb 00-100 Standard cut + Alpha cut < 0.7 | |
cdd7cb8b | 279 | |
899d014b | 280 | } else if( trainConfig == 14 ) { // No eta shift |Y| < 0.8 + AddedSignals |
281 | ||
282 | ConvCutarray[0] = "8000012002093603007200000000"; ElecCutarray[0] = "9037540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -5,5 | |
283 | ConvCutarray[1] = "8000012002093603007200000000"; ElecCutarray[1] = "9057540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -3,5 | |
284 | ConvCutarray[2] = "8000012002091603007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -5,5 | |
285 | ConvCutarray[3] = "8000012002092603007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -3,5 | |
cdd7cb8b | 286 | |
899d014b | 287 | } else if ( trainConfig == 15 ) { //No eta shift |Y| < 0.8 + AddedSignals |
cdd7cb8b | 288 | |
899d014b | 289 | ConvCutarray[0] = "8000012042093603007200000000"; ElecCutarray[0] = "9047540023510262371"; MesonCutarray[0] = "01032035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Y < 0.70 and prim and sec e |eta| < 0.75 //NOTE revisar |
290 | ConvCutarray[1] = "8000012002093603007200000000"; ElecCutarray[1] = "9047540023310263371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single prim Pt cut > 0.150 | |
291 | ConvCutarray[2] = "8000012002093603007200000000"; ElecCutarray[2] = "9047540023310263171"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + Single prim Pt cut > 0.100 | |
292 | ConvCutarray[3] = "8000012002093603007200000000"; ElecCutarray[3] = "9047540023310262271"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + DCAxy < 1 cm | |
cdd7cb8b | 293 | |
cdd7cb8b | 294 | |
899d014b | 295 | } else if ( trainConfig == 16 ) { // No eta shift |Y| < 0.8 + AddedSignals |
296 | ||
297 | ConvCutarray[0] = "8000011002493603007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.075 | |
298 | ConvCutarray[1] = "8000011002193603007200000000"; ElecCutarray[1] = "9047540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.100 | |
299 | ConvCutarray[2] = "8000011002083603007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls sec > 0.35 | |
300 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540027310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls prim > 0.60 | |
cdd7cb8b | 301 | |
302 | ||
899d014b | 303 | } else if ( trainConfig == 17 ) { //No eta shift |Y| < 0.8 + AddedSignals |
cdd7cb8b | 304 | |
899d014b | 305 | ConvCutarray[0] = "8000012002093603007200000000"; ElecCutarray[0] = "9047540023310262381"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.015 < InvMass(e+,e-) < 0.050 |
306 | ConvCutarray[1] = "8000012002093603007200000000"; ElecCutarray[1] = "9047540023310262391"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.025 < InvMass(e+,e-) < 0.035 | |
307 | ConvCutarray[2] = "8000012002093603001200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + qT < 0.1 | |
308 | ConvCutarray[3] = "8000012002093603007200000000"; ElecCutarray[3] = "9047540023310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 | |
309 | ||
310 | } else if ( trainConfig == 18 ) { //No eta shift |Y| < 0.8 + AddedSignals | |
cdd7cb8b | 311 | |
899d014b | 312 | ConvCutarray[0] = "8000012002093603007200000000"; ElecCutarray[0] = "9047540023310272371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Stardad cut +100 events background |
313 | ConvCutarray[1] = "8000012002093603007200000000"; ElecCutarray[1] = "9047540023310162371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Background method V0 multiplicity | |
314 | ConvCutarray[2] = "8000012002093603007200000000"; ElecCutarray[2] = "9047540023310262371"; MesonCutarray[2] = "01031035000000"; //standard cut Pi0 PbPb 00-100 + No extra smearing | |
315 | ConvCutarray[3] = "8000012002093603007200000000"; ElecCutarray[3] = "9047540025310262171"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 + Old Standard | |
cdd7cb8b | 316 | |
899d014b | 317 | } else if ( trainConfig == 19 ) { |
318 | ||
319 | ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 | |
320 | ConvCutarray[1] = "8000011032093603007200000000"; ElecCutarray[1] = "9047540023910262371"; MesonCutarray[1] = "01033035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 | |
321 | ConvCutarray[2] = "8000011042093603007200000000"; ElecCutarray[2] = "9047540023510262371"; MesonCutarray[2] = "01032035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75 | |
322 | ConvCutarray[3] = "8000011012093603007200000000"; ElecCutarray[3] = "9047540023610262371"; MesonCutarray[3] = "01034035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60 | |
f315320f | 323 | |
899d014b | 324 | } else if ( trainConfig == 20 ) { |
325 | ||
326 | ConvCutarray[0] = "8000012002093603007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 | |
327 | ConvCutarray[1] = "8000012032093603007200000000"; ElecCutarray[1] = "9047540023910262371"; MesonCutarray[1] = "01033035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 | |
328 | ConvCutarray[2] = "8000012042093603007200000000"; ElecCutarray[2] = "9047540023510262371"; MesonCutarray[2] = "01032035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75 | |
329 | ConvCutarray[3] = "8000012012093603007200000000"; ElecCutarray[3] = "9047540023610262371"; MesonCutarray[3] = "01034035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60 | |
58326f11 | 330 | |
2d49ea4c | 331 | } else if ( trainConfig == 21 ) { |
332 | ||
333 | ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9047540043310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 3Cls ITS | |
334 | ConvCutarray[1] = "8000011002093603007200000000"; ElecCutarray[1] = "9047540053310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4Cls ITS | |
335 | ConvCutarray[2] = "8000011002093603007200000000"; ElecCutarray[2] = "9047540063310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 5Cls ITS | |
336 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540073310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4Cls ITS no Any | |
337 | ||
338 | ||
339 | } else if ( trainConfig == 22 ) { | |
340 | ||
341 | ConvCutarray[0] = "8000012002093603007200000000"; ElecCutarray[0] = "9047540043310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 3 ITScls | |
342 | ConvCutarray[1] = "8000012002093603007200000000"; ElecCutarray[1] = "9047540053310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4 ITScls | |
343 | ConvCutarray[2] = "8000012002093603007200000000"; ElecCutarray[2] = "9047540063310262371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 5 ITScls | |
344 | ConvCutarray[3] = "8000012002093603007200000000"; ElecCutarray[3] = "9047540073310262371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4 ITScls no Any | |
345 | ||
58326f11 | 346 | } else if ( trainConfig == 23 ) { |
347 | ||
348 | ConvCutarray[0] = "8000011002493603007200000000"; ElecCutarray[0] = "9047540023310262371"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt > 0.075 | |
349 | ConvCutarray[1] = "8000011002193603007200000000"; ElecCutarray[1] = "9047540023310262371"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt > 0.100 | |
350 | ConvCutarray[2] = "8000011002093603007200000000"; ElecCutarray[2] = "9047540023310263371"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt{e} > 0.150 | |
351 | ConvCutarray[3] = "8000011002093603007200000000"; ElecCutarray[3] = "9047540023310265371"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt{e} > 0.175 | |
352 | ||
353 | ||
899d014b | 354 | } |
355 | ||
356 | ||
357 | ||
39e8c2ce | 358 | |
2d49ea4c | 359 | |
360 | ||
361 | ||
39e8c2ce | 362 | TList *ConvCutList = new TList(); |
363 | TList *MesonCutList = new TList(); | |
364 | TList *ElecCutList = new TList(); | |
365 | ||
366 | TList *HeaderList = new TList(); | |
367 | TObjString *Header1 = new TObjString("pi0_1"); | |
368 | HeaderList->Add(Header1); | |
5a104c5e | 369 | TObjString *Header3 = new TObjString("eta_2"); |
370 | HeaderList->Add(Header3); | |
39e8c2ce | 371 | |
372 | ConvCutList->SetOwner(kTRUE); | |
373 | AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts]; | |
374 | MesonCutList->SetOwner(kTRUE); | |
375 | AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts]; | |
376 | ElecCutList->SetOwner(kTRUE); | |
377 | AliDalitzElectronCuts **analysisElecCuts = new AliDalitzElectronCuts*[numberOfCuts]; | |
378 | ||
379 | ||
380 | ||
381 | for(Int_t i = 0; i<numberOfCuts; i++){ | |
382 | ||
383 | ||
384 | analysisCuts[i] = new AliConversionCuts(); | |
2d49ea4c | 385 | |
f315320f | 386 | |
58326f11 | 387 | if ( ( trainConfig >= 1 && trainConfig <= 9 ) || trainConfig == 19 || trainConfig == 21 || trainConfig == 23 ){ |
899d014b | 388 | |
389 | if (doWeighting){ | |
390 | if (generatorName.CompareTo("DPMJET")==0){ | |
cdd7cb8b | 391 | analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A"); |
899d014b | 392 | } else if (generatorName.CompareTo("HIJING")==0){ |
cdd7cb8b | 393 | analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_MBV0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A"); |
899d014b | 394 | } |
395 | } | |
396 | } | |
2d49ea4c | 397 | else if ( ( trainConfig >= 10 && trainConfig <= 18 ) || trainConfig == 20 || trainConfig == 22 ){ |
899d014b | 398 | |
399 | if (doWeighting){ | |
cdd7cb8b | 400 | analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A"); |
899d014b | 401 | } |
402 | ||
cdd7cb8b | 403 | } |
cdd7cb8b | 404 | |
2d49ea4c | 405 | if( ! analysisCuts[i]->InitializeCutsFromCutString(ConvCutarray[i].Data()) ) { |
406 | cout<<"ERROR: analysisCuts [" <<i<<"]"<<endl; | |
407 | return 0; | |
408 | } else { | |
39e8c2ce | 409 | |
2d49ea4c | 410 | if (doEtaShiftIndCuts) { |
2279d237 | 411 | |
2d49ea4c | 412 | analysisCuts[i]->DoEtaShift(doEtaShiftIndCuts); |
413 | analysisCuts[i]->SetEtaShift(stringShift); | |
414 | ||
415 | } | |
416 | ConvCutList->Add(analysisCuts[i]); | |
417 | analysisCuts[i]->SetFillCutHistograms("",kFALSE); | |
418 | analysisCuts[i]->SetAcceptedHeader(HeaderList); | |
419 | } | |
39e8c2ce | 420 | |
421 | ||
422 | ||
423 | analysisMesonCuts[i] = new AliConversionMesonCuts(); | |
424 | ||
425 | if( ! analysisMesonCuts[i]->InitializeCutsFromCutString(MesonCutarray[i].Data()) ) { | |
426 | cout<<"ERROR: analysisMesonCuts [ " <<i<<" ] "<<endl; | |
427 | return 0; | |
428 | } | |
429 | else { | |
5a104c5e | 430 | MesonCutList->Add(analysisMesonCuts[i]); |
431 | analysisMesonCuts[i]->SetFillCutHistograms(""); | |
39e8c2ce | 432 | } |
433 | ||
434 | ||
435 | TString cutName( Form("%s_%s_%s",ConvCutarray[i].Data(),ElecCutarray[i].Data(),MesonCutarray[i].Data() ) ); | |
436 | ||
437 | ||
438 | analysisElecCuts[i] = new AliDalitzElectronCuts(); | |
439 | if( !analysisElecCuts[i]->InitializeCutsFromCutString(ElecCutarray[i].Data())) { | |
440 | ||
441 | cout<< "ERROR: analysisElecCuts [ " <<i<<" ] "<<endl; | |
442 | return 0; | |
443 | } | |
444 | else { | |
445 | ElecCutList->Add(analysisElecCuts[i]); | |
446 | analysisElecCuts[i]->SetFillCutHistograms("",kFALSE,cutName); | |
447 | } | |
f315320f | 448 | |
39e8c2ce | 449 | |
450 | } | |
451 | ||
452 | ||
453 | task->SetConversionCutList(numberOfCuts,ConvCutList); | |
454 | task->SetMesonCutList(MesonCutList); | |
455 | task->SetElectronCutList(ElecCutList); | |
456 | ||
457 | task->SetMoveParticleAccordingToVertex(kTRUE); | |
2d49ea4c | 458 | task->SetProductionVertextoVGamma(kTRUE); |
39e8c2ce | 459 | |
460 | ||
461 | if(enableQAMesonTask) task->SetDoMesonQA(kTRUE); | |
5a104c5e | 462 | if(enableDoMesonChic) task->SetDoChicAnalysis(kTRUE); |
39e8c2ce | 463 | |
464 | //connect containers | |
465 | AliAnalysisDataContainer *coutput = | |
86eaf9ae | 466 | mgr->CreateContainer(Form("GammaConvDalitzV1_%i",trainConfig), TList::Class(), |
467 | AliAnalysisManager::kOutputContainer,Form("GammaConvV1Dalitz_%i.root",trainConfig)); | |
39e8c2ce | 468 | |
469 | mgr->AddTask(task); | |
470 | mgr->ConnectInput(task,0,cinput); | |
471 | mgr->ConnectOutput(task,1,coutput); | |
472 | ||
473 | return; | |
474 | ||
475 | } |