2 void AddTask_GammaConvV1_2(TString trainConfig = "pp", Bool_t isMC = kFALSE){
4 gSystem->Load("libCore.so");
5 gSystem->Load("libTree.so");
6 gSystem->Load("libGeom.so");
7 gSystem->Load("libVMC.so");
8 gSystem->Load("libPhysics.so");
9 gSystem->Load("libMinuit");
10 gSystem->Load("libSTEERBase");
11 gSystem->Load("libESD");
12 gSystem->Load("libAOD");
13 gSystem->Load("libANALYSIS");
14 gSystem->Load("libANALYSISalice");
15 gSystem->Load("libPWGGAGammaConv.so");
16 gSystem->Load("libCDB.so");
17 gSystem->Load("libSTEER.so");
18 gSystem->Load("libSTEERBase.so");
19 gSystem->Load("libTENDER.so");
20 gSystem->Load("libTENDERSupplies.so");
23 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
25 Error("AddTask_GammaConvV1", "No analysis manager found.");
30 if (trainConfig.Contains("PbPb")) IsHeavyIon=1;
31 else if (trainConfig.Contains("pPb")) IsHeavyIon=2;
33 Bool_t doEtaShift = kFALSE;
34 Int_t forceEtaShift = 0; // Carefull !!! Should be zero otherwise will shift eta cut for all periods
35 // Use doEtaShift flag for pPb or Pbp instead (1: shift +0.465, 2: shift -0.465)
37 TString cutnumber = "";
39 cutnumber = "1000000002084001001500000";
40 } else if (IsHeavyIon==2){
41 cutnumber = "8000000062084001001500000";
42 doEtaShift = kFALSE; // Only needed if Eta Range is small, but default is now +- 5.0
44 cutnumber = "0000000002084000002200000";
47 //========= Add PID Reponse to ANALYSIS manager ====
48 if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
49 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
53 //========= Add V0 Reader to ANALYSIS manager =====
54 AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
56 fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
57 fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
58 fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
61 Error("AddTask_V0ReaderV1", "No analysis manager found.");
64 AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
66 if(inputHandler->IsA()==AliESDInputHandler::Class()){
70 if(inputHandler->IsA()==AliAODInputHandler::Class()){
72 // task->SetUseSatelliteAODs(kTRUE);
75 // Set AnalysisCut Number
76 AliConversionCuts *fCuts=NULL;
78 fCuts= new AliConversionCuts(cutnumber.Data(),cutnumber.Data());
79 if(fCuts->InitializeCutsFromCutString(cutnumber.Data())){
81 fCuts->SelectCollisionCandidates(AliVEvent::kINT7);
82 fCuts->DoEtaShift(doEtaShift);
84 fCuts->ForceEtaShift(forceEtaShift);
85 fV0ReaderV1->SetConversionCuts(fCuts);
86 fCuts->SetFillCutHistograms("",kTRUE);
93 AliLog::SetGlobalLogLevel(AliLog::kInfo);
95 //================================================
97 //================================================
98 // find input container
99 //below the trunk version
100 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
101 //connect input V0Reader
102 mgr->AddTask(fV0ReaderV1);
103 mgr->ConnectInput(fV0ReaderV1,0,cinput);
105 //================================================
106 //========= Add task to the ANALYSIS manager =====
107 //================================================
109 //================================================
110 // find input container
111 AliAnalysisTaskGammaConvV1 *task=NULL;
112 task= new AliAnalysisTaskGammaConvV1("GammaConvV1");
113 task->SetIsHeavyIon(IsHeavyIon);
114 // Cut Numbers to use in Analysis
115 Int_t numberOfCuts = 2;
116 if(trainConfig.Contains("PbPb")) numberOfCuts = 3;
117 else if(trainConfig.Contains("pPb")) numberOfCuts = 2;
118 else numberOfCuts = 3;
120 TString *cutarray = new TString[numberOfCuts];
121 TString *mesonCutArray = new TString[numberOfCuts];
123 if(trainConfig.Contains("PbPb")){
124 cutarray[ 0] = "1000003042092970723220000"; mesonCutArray[ 0] = "01022045000"; // all centralities
125 cutarray[ 1] = "3010003042092970723220000"; mesonCutArray[ 1] = "01022045000"; // most central
126 cutarray[ 2] = "1680003042092970723220000"; mesonCutArray[ 2] = "01022065000"; // peripheral
127 } else if(trainConfig.Contains("pPb")){ //pA needs thighter rapidity cut y < 0.5
128 cutarray[ 0] = "8000000082093172023290000"; mesonCutArray[0] = "01629045000"; //standard cut Pi0 Pb 00-100 shifted Eta 0.4
129 cutarray[ 1] = "8000000072093172023290000"; mesonCutArray[1] = "01627045000"; //standard cut Pi0 Pb 00-100 wo shifted Eta 0.3
131 cutarray[ 0] = "0000011002093663003800000"; mesonCutArray[0] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
132 cutarray[ 1] = "0001011002093663003800000"; mesonCutArray[1] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
133 cutarray[ 2] = "0000011002093260003800000"; mesonCutArray[2] = "01631031009"; //standard cut Gamma pp 2-76TeV
137 TList *ConvCutList = new TList();
138 TList *MesonCutList = new TList();
140 TList *HeaderList = new TList();
141 // TObjString *Header1 = new TObjString("PARAM");
142 // HeaderList->Add(Header1);
143 TObjString *Header2 = new TObjString("BOX");
144 HeaderList->Add(Header2);
145 // TObjString *Header3 = new TObjString("Pythia");
146 // HeaderList->Add(Header3);
147 // TObjString *Header4 = new TObjString("Hijing");
148 // HeaderList->Add(Header4);
151 ConvCutList->SetOwner(kTRUE);
152 AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts];
153 MesonCutList->SetOwner(kTRUE);
154 AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
157 for(Int_t i = 0; i<numberOfCuts; i++){
158 analysisCuts[i] = new AliConversionCuts();
159 analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
160 if (trainConfig.Contains("pPb")){
161 analysisCuts[i]->SelectCollisionCandidates(AliVEvent::kINT7);
163 analysisCuts[i]->DoEtaShift(kTRUE);
164 analysisCuts[i]->SetEtaShift("Pbp");
167 ConvCutList->Add(analysisCuts[i]);
169 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
170 analysisMesonCuts[i] = new AliConversionMesonCuts();
171 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
172 MesonCutList->Add(analysisMesonCuts[i]);
173 analysisMesonCuts[i]->SetFillCutHistograms("");
174 analysisCuts[i]->SetAcceptedHeader(HeaderList);
178 task->SetConversionCutList(numberOfCuts,ConvCutList);
179 task->SetMesonCutList(numberOfCuts,MesonCutList);
180 task->SetMoveParticleAccordingToVertex(kTRUE);
181 task->SetDoMesonAnalysis(kTRUE);
182 if (trainConfig.Contains("pPb") || trainConfig.Contains("pp") )task->SetDoMesonQA(kTRUE); //Attention new switch for Pi0 QA
183 if (trainConfig.Contains("pPb") || trainConfig.Contains("pp") )task->SetDoPhotonQA(kTRUE); //Attention new switch small for Photon QA
186 AliAnalysisDataContainer *coutput =
187 mgr->CreateContainer("GammaConvV1", TList::Class(),
188 AliAnalysisManager::kOutputContainer,"GammaConvV1.root");
191 mgr->ConnectInput(task,0,cinput);
192 mgr->ConnectOutput(task,1,coutput);