| 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 ================================= |
| 230b7aff |
24 | gSystem->Load("libCore"); |
| 25 | gSystem->Load("libTree"); |
| 26 | gSystem->Load("libGeom"); |
| 27 | gSystem->Load("libVMC"); |
| 28 | gSystem->Load("libPhysics"); |
| 39e8c2ce |
29 | gSystem->Load("libMinuit"); |
| 30 | gSystem->Load("libSTEERBase"); |
| 31 | gSystem->Load("libESD"); |
| 32 | gSystem->Load("libAOD"); |
| 33 | gSystem->Load("libANALYSIS"); |
| 34 | gSystem->Load("libANALYSISalice"); |
| 230b7aff |
35 | gSystem->Load("libPWGGAGammaConv"); |
| 36 | gSystem->Load("libCDB"); |
| 37 | gSystem->Load("libSTEER"); |
| 38 | gSystem->Load("libSTEERBase"); |
| 39 | gSystem->Load("libTender"); |
| 40 | gSystem->Load("libTenderSupplies"); |
| 39e8c2ce |
41 | |
| 42 | |
| 43 | cout<<"Entro 0"<<endl; |
| 4fac3e9a |
44 | |
| 45 | Int_t isHeavyIon = 2; |
| 39e8c2ce |
46 | |
| 47 | // ================== GetAnalysisManager =============================== |
| 48 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); |
| 49 | if (!mgr) { |
| 86eaf9ae |
50 | Error(Form("AddTask_GammaConvDalitzV1_pPb_%i",trainConfig), "No analysis manager found."); |
| 39e8c2ce |
51 | return ; |
| 52 | } |
| 53 | |
| 54 | // ================== GetInputEventHandler ============================= |
| 55 | AliVEventHandler *inputHandler=mgr->GetInputEventHandler(); |
| 56 | |
| 57 | //========= Add PID Reponse to ANALYSIS manager ==== |
| 58 | if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){ |
| 59 | gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C"); |
| 60 | AddTaskPIDResponse(isMC); |
| 61 | } |
| 62 | |
| 63 | //========= Set Cutnumber for V0Reader ================================ |
| 4fac3e9a |
64 | TString cutnumberEvent = "8000000"; |
| 65 | |
| 3013611f |
66 | TString cutnumberPhoton="060000084001001500000000"; //Online V0 finder |
| 4fac3e9a |
67 | |
| 30b9a400 |
68 | TString ElecCuts = "90005400000002000000"; //Electron Cuts |
| 4fac3e9a |
69 | |
| 5a104c5e |
70 | Bool_t doEtaShift = kFALSE; |
| 4fac3e9a |
71 | |
| 39e8c2ce |
72 | |
| 73 | |
| 74 | AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer(); |
| 75 | |
| 76 | //========= Add V0 Reader to ANALYSIS manager if not yet existent ===== |
| 77 | if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){ |
| 4fac3e9a |
78 | AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1"); |
| 79 | |
| 80 | fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE); |
| 81 | fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output |
| 82 | fV0ReaderV1->SetUseAODConversionPhoton(kTRUE); |
| 83 | |
| 84 | if (!mgr) { |
| 85 | Error("AddTask_V0ReaderV1", "No analysis manager found."); |
| 86 | return; |
| 87 | } |
| 88 | |
| 89 | AliConvEventCuts *fEventCuts=NULL; |
| 90 | if(cutnumberEvent!=""){ |
| 91 | fEventCuts= new AliConvEventCuts(cutnumberEvent.Data(),cutnumberEvent.Data()); |
| 92 | fEventCuts->SetPreSelectionCutFlag(kTRUE); |
| 93 | if(fEventCuts->InitializeCutsFromCutString(cutnumberEvent.Data())){ |
| 94 | fEventCuts->DoEtaShift(doEtaShift); |
| 95 | fV0ReaderV1->SetEventCuts(fEventCuts); |
| 96 | fEventCuts->SetFillCutHistograms("",kTRUE); |
| 97 | } |
| 98 | } |
| 39e8c2ce |
99 | |
| 4fac3e9a |
100 | // Set AnalysisCut Number |
| 101 | AliConversionPhotonCuts *fCuts=NULL; |
| 102 | if(cutnumberPhoton!=""){ |
| 103 | fCuts= new AliConversionPhotonCuts(cutnumberPhoton.Data(),cutnumberPhoton.Data()); |
| 104 | fCuts->SetPreSelectionCutFlag(kTRUE); |
| 105 | fCuts->SetIsHeavyIon(isHeavyIon); |
| 106 | if(fCuts->InitializeCutsFromCutString(cutnumberPhoton.Data())){ |
| 107 | fV0ReaderV1->SetConversionCuts(fCuts); |
| 108 | fCuts->SetFillCutHistograms("",kTRUE); |
| 109 | } |
| 110 | } |
| 111 | if(inputHandler->IsA()==AliAODInputHandler::Class()){ |
| 112 | // AOD mode |
| 113 | cout << "AOD handler: adding " << cutnumberAODBranch.Data() << " as conversion branch" << endl; |
| 114 | fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data())); |
| 115 | } |
| 116 | fV0ReaderV1->Init(); |
| 117 | |
| 118 | AliLog::SetGlobalLogLevel(AliLog::kInfo); |
| 119 | |
| 120 | //connect input V0Reader |
| 121 | mgr->AddTask(fV0ReaderV1); |
| 122 | mgr->ConnectInput(fV0ReaderV1,0,cinput); |
| 123 | |
| 124 | } |
| 39e8c2ce |
125 | |
| 126 | //================================================ |
| 127 | //========= Add Electron Selector ================ |
| 128 | |
| 129 | |
| 130 | if( !(AliDalitzElectronSelector*)mgr->GetTask("ElectronSelector") ){ |
| 131 | |
| 132 | AliDalitzElectronSelector *fElectronSelector = new AliDalitzElectronSelector("ElectronSelector"); |
| 133 | |
| 134 | // Set AnalysisCut Number |
| 135 | |
| 136 | AliDalitzElectronCuts *fElecCuts=0; |
| 137 | |
| 4fac3e9a |
138 | |
| 39e8c2ce |
139 | |
| 140 | if( ElecCuts!=""){ |
| 141 | |
| 142 | fElecCuts= new AliDalitzElectronCuts(ElecCuts.Data(),ElecCuts.Data()); |
| 143 | |
| 144 | if(fElecCuts->InitializeCutsFromCutString(ElecCuts.Data())){ |
| 145 | |
| 146 | fElectronSelector->SetDalitzElectronCuts(fElecCuts); |
| 147 | |
| 148 | fElecCuts->SetFillCutHistograms("",kTRUE); |
| 149 | |
| 150 | } |
| 151 | |
| 152 | } |
| 153 | |
| 154 | fElectronSelector->Init(); |
| 155 | mgr->AddTask(fElectronSelector); |
| 5a104c5e |
156 | |
| 157 | AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer(); |
| 158 | |
| 159 | //connect input V0Reader |
| 160 | |
| 161 | mgr->ConnectInput (fElectronSelector,0,cinput1); |
| 39e8c2ce |
162 | |
| 163 | } |
| 164 | |
| 165 | |
| 166 | |
| 167 | cout<<"Entro"<<endl; |
| 168 | //================================================ |
| 169 | //========= Add task to the ANALYSIS manager ===== |
| 170 | //================================================ |
| 171 | // find input container |
| 172 | |
| 173 | |
| 174 | |
| 175 | AliAnalysisTaskGammaConvDalitzV1 *task=NULL; |
| 176 | |
| 86eaf9ae |
177 | task= new AliAnalysisTaskGammaConvDalitzV1(Form("GammaConvDalitzV1_%i",trainConfig)); |
| 39e8c2ce |
178 | |
| 4fac3e9a |
179 | task->SetIsHeavyIon(isHeavyIon); |
| 39e8c2ce |
180 | task->SetIsMC(isMC); |
| 181 | |
| 182 | |
| 183 | |
| 184 | // Cut Numbers to use in Analysis |
| 899d014b |
185 | Int_t numberOfCuts = 4; |
| 39e8c2ce |
186 | |
| 4fac3e9a |
187 | |
| 188 | |
| 189 | TString *eventCutArray = new TString[numberOfCuts]; |
| 190 | TString *photonCutArray = new TString[numberOfCuts]; |
| 39e8c2ce |
191 | TString *ElecCutarray = new TString[numberOfCuts]; |
| 39e8c2ce |
192 | TString *MesonCutarray = new TString[numberOfCuts]; |
| 193 | |
| 194 | Bool_t doEtaShiftIndCuts = kFALSE; |
| 39e8c2ce |
195 | TString stringShift = ""; |
| 196 | |
| 197 | // Shifting in pPb direction |
| 198 | |
| 899d014b |
199 | doEtaShiftIndCuts = kFALSE; |
| 39e8c2ce |
200 | stringShift = "pPb"; |
| 201 | |
| 86eaf9ae |
202 | |
| 2279d237 |
203 | |
| 2279d237 |
204 | |
| 899d014b |
205 | if( trainConfig == 1 ) { // No eta shift |Y| < 0.8 |
| cdd7cb8b |
206 | |
| 3013611f |
207 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
| 208 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007900000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 15 |
| 209 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007800000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 20 |
| 210 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007100000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 50 |
| cdd7cb8b |
211 | |
| 899d014b |
212 | } else if( trainConfig == 2 ) { // No eta shift |Y| < 0.8 |
| 4fac3e9a |
213 | |
| 3013611f |
214 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603002200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.7 |
| 215 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603003200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.5 |
| 216 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093653007200000000"; ElecCutarray[2] = "90475400233102623710"; 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 |
| 217 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093601007200000000"; ElecCutarray[3] = "90475400233102623710"; 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 |
218 | |
| 899d014b |
219 | } else if( trainConfig == 3 ) { // No eta shift |Y| < 0.8 |
| 220 | |
| 3013611f |
221 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093803007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 2.0sigmas Low and 1 High momentum |
| 222 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90435400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and 0 High momentum |
| 223 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90477400233102623710"; 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 |
| 224 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475200233102623710"; 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 |
225 | |
| 899d014b |
226 | } else if( trainConfig == 4 ) { // No eta shift |Y| < 0.8 |
| 227 | |
| 3013611f |
228 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90425400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and -1 High momentum |
| 229 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400133102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + SPD first layer |
| 230 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233302623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + PsiPair cut 0.52 |
| 231 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031005009000"; //standard cut Pi0 pPb 00-100 Standard cut + Alpha cut < 0.7 |
| cdd7cb8b |
232 | |
| 899d014b |
233 | } else if( trainConfig == 5 ) { // No eta shift |Y| < 0.8 |
| 234 | |
| 3013611f |
235 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90375400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -5,5 |
| 236 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90575400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -3,5 |
| 237 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000091603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -5,5 |
| 238 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000092603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -3,5 |
| cdd7cb8b |
239 | |
| 899d014b |
240 | } else if ( trainConfig == 6 ) { //No eta shift |Y| < 0.8 |
| cdd7cb8b |
241 | |
| 3013611f |
242 | eventCutArray[0]="8000011"; photonCutArray[0] = "042000093603007200000000"; ElecCutarray[0] = "90475400235102623710"; 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 |
| 243 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102633710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single prim Pt cut > 0.150 |
| 244 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102631710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + Single prim Pt cut > 0.100 |
| 245 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102622710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + DCAxy < 1 cm |
| cdd7cb8b |
246 | |
| 247 | |
| 899d014b |
248 | } else if ( trainConfig == 7 ) { // No eta shift |Y| < 0.8 |
| 249 | |
| 3013611f |
250 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000493603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.075 |
| 251 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000193603007200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.100 |
| 252 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000083603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls sec > 0.35 |
| 253 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400273102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls prim > 0.60 |
| cdd7cb8b |
254 | |
| cdd7cb8b |
255 | |
| 899d014b |
256 | } else if ( trainConfig == 8 ) { //No eta shift |Y| < 0.8 |
| cdd7cb8b |
257 | |
| 3013611f |
258 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623810"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.015 < InvMass(e+,e-) < 0.050 |
| 259 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102623910"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.025 < InvMass(e+,e-) < 0.035 |
| 260 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603001200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + qT < 0.1 |
| 261 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
| cdd7cb8b |
262 | |
| 899d014b |
263 | } else if ( trainConfig == 9 ) { //No eta shift |Y| < 0.8 |
| 264 | |
| cdd7cb8b |
265 | |
| 3013611f |
266 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102723710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Stardad cut +100 events background |
| 267 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233101623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Background method V0 multiplicity |
| 268 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035000000"; //standard cut Pi0 PbPb 00-100 + No extra smearing |
| 269 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400253102621710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 + Old Standard |
| cdd7cb8b |
270 | |
| 899d014b |
271 | } else if( trainConfig == 10 ) { // No eta shift |Y| < 0.8 + AddedSignals |
| cdd7cb8b |
272 | |
| 3013611f |
273 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
| 274 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007900000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 15 |
| 275 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007800000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 20 |
| 276 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007100000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Chi2 < 50 |
| cdd7cb8b |
277 | |
| 899d014b |
278 | } else if( trainConfig == 11 ) { // No eta shift |Y| < 0.8 + AddedSignals |
| 279 | |
| 3013611f |
280 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603002200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.7 |
| 281 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603003200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Qt < 0.5 |
| 282 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093653007200000000"; ElecCutarray[2] = "90475400233102623710"; 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 |
| 283 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093601007200000000"; ElecCutarray[3] = "90475400233102623710"; 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 |
284 | |
| 899d014b |
285 | } else if( trainConfig == 12 ) { // No eta shift |Y| < 0.8 + AddedSignals |
| 286 | |
| 3013611f |
287 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093803007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec sec 2.0sigmas Low and 1 High momentum |
| 288 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90435400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and 0 High momentum |
| 289 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90477400233102623710"; 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 |
| 290 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475200233102623710"; 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 |
291 | |
| 899d014b |
292 | } else if( trainConfig == 13 ) { // No eta shift |Y| < 0.8 + AddedSignals |
| 293 | |
| 3013611f |
294 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90425400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec primary 2.0sigmas Low and -1 High momentum |
| 295 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400133102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + SPD first layer |
| 296 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233302623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + PsiPair cut 0.52 |
| 297 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031005009000"; //standard cut Pi0 pPb 00-100 Standard cut + Alpha cut < 0.7 |
| cdd7cb8b |
298 | |
| 899d014b |
299 | } else if( trainConfig == 14 ) { // No eta shift |Y| < 0.8 + AddedSignals |
| 300 | |
| 3013611f |
301 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90375400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -5,5 |
| 302 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90575400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary electron -3,5 |
| 303 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000091603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -5,5 |
| 304 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000092603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -3,5 |
| cdd7cb8b |
305 | |
| 899d014b |
306 | } else if ( trainConfig == 15 ) { //No eta shift |Y| < 0.8 + AddedSignals |
| cdd7cb8b |
307 | |
| 3013611f |
308 | eventCutArray[0]="8000012"; photonCutArray[0] = "042000093603007200000000"; ElecCutarray[0] = "90475400235102623710"; 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 |
| 309 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102633710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single prim Pt cut > 0.150 |
| 310 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102631710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + Single prim Pt cut > 0.100 |
| 311 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102622710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + DCAxy < 1 cm |
| cdd7cb8b |
312 | |
| cdd7cb8b |
313 | |
| 4fac3e9a |
314 | } else if ( trainConfig == 16 ) { // No eta shift |Y| < 0.8 |
| 899d014b |
315 | |
| 3013611f |
316 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000493603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.075 |
| 317 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000193603007200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Single sec Pt cut > 0.100 |
| 318 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000083603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls sec > 0.35 |
| 319 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400273102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Findable Cls prim > 0.60 |
| cdd7cb8b |
320 | |
| 321 | |
| 899d014b |
322 | } else if ( trainConfig == 17 ) { //No eta shift |Y| < 0.8 + AddedSignals |
| cdd7cb8b |
323 | |
| 3013611f |
324 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623810"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.015 < InvMass(e+,e-) < 0.050 |
| 325 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102623910"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.025 < InvMass(e+,e-) < 0.035 |
| 326 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603001200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + qT < 0.1 |
| 327 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
| 899d014b |
328 | |
| 329 | } else if ( trainConfig == 18 ) { //No eta shift |Y| < 0.8 + AddedSignals |
| cdd7cb8b |
330 | |
| 3013611f |
331 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102723710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Stardad cut +100 events background |
| 332 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233101623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + Background method V0 multiplicity |
| 333 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035000000"; //standard cut Pi0 PbPb 00-100 + No extra smearing |
| 334 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400253102621710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 + Old Standard |
| cdd7cb8b |
335 | |
| 899d014b |
336 | } else if ( trainConfig == 19 ) { |
| 337 | |
| 3013611f |
338 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
| 339 | eventCutArray[1]="8000011"; photonCutArray[1] = "032000093603007200000000"; ElecCutarray[1] = "90475400239102623710"; 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 |
| 340 | eventCutArray[2]="8000011"; photonCutArray[2] = "042000093603007200000000"; ElecCutarray[2] = "90475400235102623710"; 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 |
| 341 | eventCutArray[3]="8000011"; photonCutArray[3] = "012000093603007200000000"; ElecCutarray[3] = "90475400236102623710"; 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 |
342 | |
| 899d014b |
343 | } else if ( trainConfig == 20 ) { |
| 344 | |
| 3013611f |
345 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 |
| 346 | eventCutArray[1]="8000012"; photonCutArray[1] = "032000093603007200000000"; ElecCutarray[1] = "90475400239102623710"; 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 |
| 347 | eventCutArray[2]="8000012"; photonCutArray[2] = "042000093603007200000000"; ElecCutarray[2] = "90475400235102623710"; 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 |
| 348 | eventCutArray[3]="8000012"; photonCutArray[3] = "012000093603007200000000"; ElecCutarray[3] = "90475400236102623710"; 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 |
349 | |
| 2d49ea4c |
350 | } else if ( trainConfig == 21 ) { |
| 351 | |
| 3013611f |
352 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400433102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 3Cls ITS |
| 353 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400533102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4Cls ITS |
| 354 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400633102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 5Cls ITS |
| 355 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400733102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4Cls ITS no Any |
| 2d49ea4c |
356 | |
| 357 | |
| 358 | } else if ( trainConfig == 22 ) { |
| 359 | |
| 3013611f |
360 | eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400433102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 3 ITScls |
| 361 | eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400533102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4 ITScls |
| 362 | eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400633102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 5 ITScls |
| 363 | eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400733102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + 4 ITScls no Any |
| 2d49ea4c |
364 | |
| 58326f11 |
365 | } else if ( trainConfig == 23 ) { |
| 366 | |
| 3013611f |
367 | eventCutArray[0]="8000011"; photonCutArray[0] = "002000493603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt > 0.075 |
| 368 | eventCutArray[1]="8000011"; photonCutArray[1] = "002000193603007200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt > 0.100 |
| 369 | eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102633710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt{e} > 0.150 |
| 370 | eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102653710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + Pt{e} > 0.175 |
| 58326f11 |
371 | |
| 372 | |
| 7a5301ed |
373 | } else if ( trainConfig == 24 ) { |
| 374 | |
| 3013611f |
375 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400833202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + new psiPair cut 0.60, 0.0 0.12 |
| 376 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400133202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kFirst + new psiPair cut + 0.60 0.0 0.12 |
| 377 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 |
| 378 | eventCutArray[3] = "8000011"; photonCutArray[3] = "005000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 10cm 0.60, 0.0 0.12 |
| 7a5301ed |
379 | |
| 380 | } else if ( trainConfig == 25 ) { |
| 381 | |
| 3013611f |
382 | eventCutArray[0] = "8000011"; photonCutArray[0] = "008000093603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 12.5cm 0.60, 0.0 0.12 |
| 383 | eventCutArray[1] = "8000011"; photonCutArray[1] = "006000093603007200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 20 cm 0.60, 0.0 0.12 |
| 384 | eventCutArray[2] = "8000011"; photonCutArray[2] = "007000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 35 cm 0.60, 0.0 0.12 |
| 385 | eventCutArray[3] = "8000011"; photonCutArray[3] = "009000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 7.5 cm 0.60, 0.0 0.12 |
| 7a5301ed |
386 | |
| 387 | } else if ( trainConfig == 26 ) { |
| 388 | |
| 3013611f |
389 | eventCutArray[0] = "8000011"; photonCutArray[0] = "008000093603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 12.5cm 0.60, 0.0 0.12 |
| 390 | eventCutArray[1] = "8000011"; photonCutArray[1] = "006000093603007200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 20 cm 0.60, 0.0 0.12 |
| 391 | eventCutArray[2] = "8000011"; photonCutArray[2] = "007000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 35 cm 0.60, 0.0 0.12 |
| 392 | eventCutArray[3] = "8000011"; photonCutArray[3] = "009000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut + gammaR > 7.5 cm 0.60, 0.0 0.12 |
| 7a5301ed |
393 | |
| 394 | } else if ( trainConfig == 27 ) { |
| 395 | |
| 3013611f |
396 | eventCutArray[0]= "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202613710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Ptprim > 100 MeV |
| 397 | eventCutArray[1]= "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233202633710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Ptprim > 150 MeV |
| 398 | eventCutArray[2]= "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202603710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Ptprim > 75 MeV |
| 399 | eventCutArray[3]= "8000011"; photonCutArray[3] = "002000493603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Ptsec > 75 MeV |
| 7a5301ed |
400 | |
| 401 | } else if ( trainConfig == 28 ) { |
| 402 | |
| 3013611f |
403 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000193603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Pt sec > 100 MeV |
| 404 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233202653710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Pt prim > 175 GeV |
| 405 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603001200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Qt > 0.1 |
| 406 | eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603002200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + Qt > 0.07 |
| 7a5301ed |
407 | |
| 408 | } else if ( trainConfig == 29 ) { |
| 409 | |
| 3013611f |
410 | eventCutArray[0]= "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90375400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx primary electron -5,5 |
| 411 | eventCutArray[1]= "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90575400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx primary electron -3,5 |
| 412 | eventCutArray[2]= "8000011"; photonCutArray[2] = "002000091603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx secondary electron -5,5 |
| 413 | eventCutArray[3]= "8000011"; photonCutArray[3] = "002000092603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx secondary electron -3,5 |
| 7a5301ed |
414 | |
| 415 | } else if ( trainConfig == 30 ) { |
| 416 | |
| 3013611f |
417 | eventCutArray[0]= "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90435400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx pion rejec primary 2.0 sigmas Low and 0 High momentum |
| 418 | eventCutArray[1]= "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90425400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx pion rejec primary 2.0 sigmas Low and -1 High momentum |
| 419 | eventCutArray[2]= "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90477400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx pion rejec primary 0.3 GeV Low and 3.5 High momentum |
| 420 | eventCutArray[3]= "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475200233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + dEdx pion rejec primary 0.5 GeV Low and 5.0 High momentum*/ |
| 7a5301ed |
421 | |
| 422 | } else if ( trainConfig == 31 ) { |
| 423 | |
| 3013611f |
424 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093653007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair cut + dEdx pion rejec sec 0.3 GeV Low and 3.5 High momentum |
| 425 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093601007200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair cut + dEdx pion rejec sec 0.5 GeV Low and 5.0 High momentum |
| 426 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093803007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair cut + dEdx pion rejec sec 2.0 sigmas Low and 1 High momentum |
| 427 | eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007900000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair cut + Chi2 < 15 |
| 7a5301ed |
428 | |
| 429 | } else if ( trainConfig == 32 ) { |
| 430 | |
| 3013611f |
431 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007800000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair cut + Chi2 < 20 |
| 432 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007100000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair cut + Chi2 < 50 |
| 433 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031005009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair cut + Alpha < 0.7 |
| 434 | eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233002623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + no psi pair + weights |
| 7a5301ed |
435 | |
| 436 | } else if ( trainConfig == 33 ) { |
| 437 | |
| 3013611f |
438 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202622710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + DCAxy < 1 cm |
| 439 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233202623810"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + 0.015 < InvMass(e+,e-) < 0.050 |
| 440 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202623910"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + 0.025 < InvMass(e+,e-) < 0.035 |
| 441 | eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233202723710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + 100 events background |
| 7a5301ed |
442 | |
| 443 | } else if ( trainConfig == 34 ) { |
| 444 | |
| 3013611f |
445 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233201623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + Background method V0 multiplicity |
| 446 | eventCutArray[1] = "8000011"; photonCutArray[1] = "032000093603007200000000"; ElecCutarray[1] = "90475400239202623710"; MesonCutarray[1] = "01033035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 |
| 447 | eventCutArray[2] = "8000011"; photonCutArray[2] = "042000093603007200000000"; ElecCutarray[2] = "90475400235202623710"; MesonCutarray[2] = "01032035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75 |
| 448 | eventCutArray[3] = "8000011"; photonCutArray[3] = "012000093603007200000000"; ElecCutarray[3] = "90475400236202623710"; MesonCutarray[3] = "01034035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60 |
| 7a5301ed |
449 | |
| 450 | } else if ( trainConfig == 35 ) { |
| 451 | |
| 3013611f |
452 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 |
| 453 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603003200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + Qt < 0.05 |
| 454 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000083603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + Findable Cls > 0.35 Secondary |
| 455 | eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400273202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + Findable Cls > 0.60 primary |
| 7a5301ed |
456 | |
| 457 | } else if ( trainConfig == 36 ) { |
| 458 | |
| 3013611f |
459 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400433202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + ITScls >= 3 |
| 460 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400533202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + ITScls >= 4 |
| 461 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400633202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + ITScls >= 5 |
| 462 | eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035000000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + No extra smearing |
| 7a5301ed |
463 | |
| 464 | } else if ( trainConfig == 37 ) { |
| 465 | |
| 3013611f |
466 | eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202623700"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 |
| 467 | eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400833202623700"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kBoth + new psiPair cut 0.60, 0.0 0.12 |
| 468 | eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400133202623700"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kFirst + new psiPair Cut 0.60, 0.0 0.12 |
| 469 | eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233002623700"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + no psi pair + no weights |
| 7a5301ed |
470 | |
| 899d014b |
471 | } |
| 7a5301ed |
472 | |
| 899d014b |
473 | |
| 474 | |
| 475 | |
| 39e8c2ce |
476 | |
| 2d49ea4c |
477 | |
| 478 | |
| 4fac3e9a |
479 | TList *EventCutList = new TList(); |
| 39e8c2ce |
480 | TList *ConvCutList = new TList(); |
| 481 | TList *MesonCutList = new TList(); |
| 482 | TList *ElecCutList = new TList(); |
| 483 | |
| 484 | TList *HeaderList = new TList(); |
| 485 | TObjString *Header1 = new TObjString("pi0_1"); |
| 486 | HeaderList->Add(Header1); |
| 5a104c5e |
487 | TObjString *Header3 = new TObjString("eta_2"); |
| 488 | HeaderList->Add(Header3); |
| 39e8c2ce |
489 | |
| 4fac3e9a |
490 | EventCutList->SetOwner(kTRUE); |
| 491 | AliConvEventCuts **analysisEventCuts = new AliConvEventCuts*[numberOfCuts]; |
| 39e8c2ce |
492 | ConvCutList->SetOwner(kTRUE); |
| 4fac3e9a |
493 | AliConversionPhotonCuts **analysisCuts = new AliConversionPhotonCuts*[numberOfCuts]; |
| 39e8c2ce |
494 | MesonCutList->SetOwner(kTRUE); |
| 495 | AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts]; |
| 496 | ElecCutList->SetOwner(kTRUE); |
| 497 | AliDalitzElectronCuts **analysisElecCuts = new AliDalitzElectronCuts*[numberOfCuts]; |
| 498 | |
| 499 | |
| 500 | |
| 501 | for(Int_t i = 0; i<numberOfCuts; i++){ |
| 502 | |
| 4fac3e9a |
503 | analysisEventCuts[i] = new AliConvEventCuts(); |
| f315320f |
504 | |
| 7a5301ed |
505 | if ( ( trainConfig >= 1 && trainConfig <= 9 ) || trainConfig == 19 || trainConfig == 21 || trainConfig == 23 || ( trainConfig >= 24 && trainConfig <=36 ) ){ |
| 899d014b |
506 | |
| 507 | if (doWeighting){ |
| 508 | if (generatorName.CompareTo("DPMJET")==0){ |
| 4fac3e9a |
509 | analysisEventCuts[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 |
510 | } else if (generatorName.CompareTo("HIJING")==0){ |
| 4fac3e9a |
511 | analysisEventCuts[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 |
512 | } |
| 513 | } |
| 514 | } |
| 2d49ea4c |
515 | else if ( ( trainConfig >= 10 && trainConfig <= 18 ) || trainConfig == 20 || trainConfig == 22 ){ |
| 899d014b |
516 | |
| 4fac3e9a |
517 | if (doWeighting){ |
| 518 | analysisEventCuts[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"); |
| 519 | } |
| 899d014b |
520 | |
| 4fac3e9a |
521 | } |
| cdd7cb8b |
522 | |
| 4fac3e9a |
523 | |
| 524 | analysisEventCuts[i]->InitializeCutsFromCutString(eventCutArray[i].Data()); |
| 525 | if (doEtaShiftIndCuts) { |
| 526 | analysisEventCuts[i]->DoEtaShift(doEtaShiftIndCuts); |
| 527 | analysisEventCuts[i]->SetEtaShift(stringShift); |
| 528 | } |
| 529 | EventCutList->Add(analysisEventCuts[i]); |
| 530 | analysisEventCuts[i]->SetFillCutHistograms("",kFALSE); |
| 531 | analysisEventCuts[i]->SetAcceptedHeader(HeaderList); |
| 532 | |
| 533 | analysisCuts[i] = new AliConversionPhotonCuts(); |
| 534 | |
| 535 | if( ! analysisCuts[i]->InitializeCutsFromCutString(photonCutArray[i].Data()) ) { |
| 536 | cout<<"ERROR: analysisCuts [" <<i<<"]"<<endl; |
| 537 | return 0; |
| 538 | } |
| 539 | analysisCuts[i]->SetIsHeavyIon(isHeavyIon); |
| 540 | ConvCutList->Add(analysisCuts[i]); |
| 541 | analysisCuts[i]->SetFillCutHistograms("",kFALSE); |
| 2d49ea4c |
542 | |
| 4fac3e9a |
543 | |
| 544 | |
| 39e8c2ce |
545 | |
| 546 | analysisMesonCuts[i] = new AliConversionMesonCuts(); |
| 547 | |
| 548 | if( ! analysisMesonCuts[i]->InitializeCutsFromCutString(MesonCutarray[i].Data()) ) { |
| 549 | cout<<"ERROR: analysisMesonCuts [ " <<i<<" ] "<<endl; |
| 550 | return 0; |
| 551 | } |
| 4fac3e9a |
552 | MesonCutList->Add(analysisMesonCuts[i]); |
| 553 | analysisMesonCuts[i]->SetFillCutHistograms(""); |
| 554 | |
| 39e8c2ce |
555 | |
| 556 | |
| c524a081 |
557 | //TString cutName( Form("%s_%s_%s",photonCutArray[i].Data(),ElecCutarray[i].Data(),MesonCutarray[i].Data() ) ); |
| 39e8c2ce |
558 | |
| 559 | |
| 560 | analysisElecCuts[i] = new AliDalitzElectronCuts(); |
| 561 | if( !analysisElecCuts[i]->InitializeCutsFromCutString(ElecCutarray[i].Data())) { |
| 562 | |
| 563 | cout<< "ERROR: analysisElecCuts [ " <<i<<" ] "<<endl; |
| 564 | return 0; |
| 565 | } |
| 4fac3e9a |
566 | ElecCutList->Add(analysisElecCuts[i]); |
| c524a081 |
567 | analysisElecCuts[i]->SetFillCutHistograms("",kFALSE,ElecCutarray[i].Data()); |
| 4fac3e9a |
568 | |
| 39e8c2ce |
569 | |
| 570 | } |
| 571 | |
| 4fac3e9a |
572 | task->SetEventCutList(numberOfCuts,EventCutList); |
| 39e8c2ce |
573 | task->SetConversionCutList(numberOfCuts,ConvCutList); |
| 574 | task->SetMesonCutList(MesonCutList); |
| 575 | task->SetElectronCutList(ElecCutList); |
| 576 | |
| 577 | task->SetMoveParticleAccordingToVertex(kTRUE); |
| 2d49ea4c |
578 | task->SetProductionVertextoVGamma(kTRUE); |
| 39e8c2ce |
579 | |
| 580 | |
| 581 | if(enableQAMesonTask) task->SetDoMesonQA(kTRUE); |
| 5a104c5e |
582 | if(enableDoMesonChic) task->SetDoChicAnalysis(kTRUE); |
| 39e8c2ce |
583 | |
| 584 | //connect containers |
| 585 | AliAnalysisDataContainer *coutput = |
| 86eaf9ae |
586 | mgr->CreateContainer(Form("GammaConvDalitzV1_%i",trainConfig), TList::Class(), |
| 587 | AliAnalysisManager::kOutputContainer,Form("GammaConvV1Dalitz_%i.root",trainConfig)); |
| 39e8c2ce |
588 | |
| 589 | mgr->AddTask(task); |
| 590 | mgr->ConnectInput(task,0,cinput); |
| 591 | mgr->ConnectOutput(task,1,coutput); |
| 592 | |
| 593 | return; |
| 594 | |
| 595 | } |
git://git.uio.no / u / mrichter / AliRoot.git / blame