2 void AddTask_GammaConvV1(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;
35 TString cutnumber = "";
37 cutnumber = "1000000002084001001500000";
38 } else if (IsHeavyIon==2){
39 cutnumber = "8000000062084001001500000";
40 doEtaShift = kFALSE; // Only needed if Eta Range is small, but default is now +- 5.0
42 cutnumber = "0000000002084000002200000";
45 //========= Add PID Reponse to ANALYSIS manager ====
46 if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
47 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
51 //========= Add V0 Reader to ANALYSIS manager =====
52 AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
54 fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
55 fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
56 fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
59 Error("AddTask_V0ReaderV1", "No analysis manager found.");
62 AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
64 if(inputHandler->IsA()==AliESDInputHandler::Class()){
68 if(inputHandler->IsA()==AliAODInputHandler::Class()){
70 // task->SetUseSatelliteAODs(kTRUE);
73 // Set AnalysisCut Number
74 AliConversionCuts *fCuts=NULL;
76 fCuts= new AliConversionCuts(cutnumber.Data(),cutnumber.Data());
77 if(fCuts->InitializeCutsFromCutString(cutnumber.Data())){
79 fCuts->SelectCollisionCandidates(AliVEvent::kINT7);
80 fCuts->DoEtaShift(doEtaShift);
82 fV0ReaderV1->SetConversionCuts(fCuts);
83 fCuts->SetFillCutHistograms("",kTRUE);
90 AliLog::SetGlobalLogLevel(AliLog::kInfo);
92 //================================================
94 //================================================
95 // find input container
96 //below the trunk version
97 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
98 //connect input V0Reader
99 mgr->AddTask(fV0ReaderV1);
100 mgr->ConnectInput(fV0ReaderV1,0,cinput);
102 //================================================
103 //========= Add task to the ANALYSIS manager =====
104 //================================================
106 //================================================
107 // find input container
108 AliAnalysisTaskGammaConvV1 *task=NULL;
109 task= new AliAnalysisTaskGammaConvV1("GammaConvV1");
110 task->SetIsHeavyIon(IsHeavyIon);
111 // Cut Numbers to use in Analysis
112 Int_t numberOfCuts = 2;
113 if(trainConfig.Contains("PbPb")) numberOfCuts = 3;
114 else if(trainConfig.Contains("pPb")) numberOfCuts = 4;
115 else numberOfCuts = 3;
117 TString *cutarray = new TString[numberOfCuts];
118 TString *mesonCutArray = new TString[numberOfCuts];
120 if(trainConfig.Contains("PbPb")){
121 cutarray[ 0] = "1010001042092970023220000"; mesonCutArray[ 0] = "01522045000"; // Standard cut 0-10
122 cutarray[ 1] = "3120001042092970023220000"; mesonCutArray[ 1] = "01522045000"; // Standard cut 5-10
123 cutarray[ 2] = "1120001042092970023220000"; mesonCutArray[ 2] = "01522045000"; // Standard cut 10-20
124 } else if(trainConfig.Contains("pPb")){ //pA needs thighter rapidity cut y < 0.5
125 cutarray[ 0] = "8020000082093172023290000"; mesonCutArray[0] = "01629045000"; //standard cut Pi0 Pb 00-20 shifted Eta 0.4
126 cutarray[ 1] = "8240000082093172023290000"; mesonCutArray[1] = "01629045000"; //standard cut Pi0 Pb 20-40 shifted Eta 0.4
127 cutarray[ 2] = "8460000082093172023290000"; mesonCutArray[2] = "01629045000"; //standard cut Pi0 Pb 40-60 shifted Eta 0.4
128 cutarray[ 3] = "8600000082093172023290000"; mesonCutArray[3] = "01629045000"; //standard cut Pi0 Pb 60-100 shifted Eta 0.4
130 cutarray[ 0] = "0000012002093663003800000"; mesonCutArray[0] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD , only boxes
131 cutarray[ 1] = "0001012002093663003800000"; mesonCutArray[1] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND , only boxes
132 cutarray[ 2] = "0000012002093260003800000"; mesonCutArray[2] = "01631031009"; //standard cut Gamma pp 2-76TeV , only boxes
136 TList *ConvCutList = new TList();
137 TList *MesonCutList = new TList();
139 TList *HeaderList = new TList();
140 // TObjString *Header1 = new TObjString("PARAM");
141 // HeaderList->Add(Header1);
142 TObjString *Header2 = new TObjString("BOX");
143 HeaderList->Add(Header2);
144 // TObjString *Header3 = new TObjString("Pythia");
145 // HeaderList->Add(Header3);
146 // TObjString *Header4 = new TObjString("Hijing");
147 // HeaderList->Add(Header4);
150 ConvCutList->SetOwner(kTRUE);
151 AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts];
152 MesonCutList->SetOwner(kTRUE);
153 AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
156 for(Int_t i = 0; i<numberOfCuts; i++){
157 analysisCuts[i] = new AliConversionCuts();
158 analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
159 if (trainConfig.Contains("pPb")){
160 analysisCuts[i]->SelectCollisionCandidates(AliVEvent::kINT7);
162 analysisCuts[i]->DoEtaShift(kTRUE);
163 analysisCuts[i]->SetEtaShift("Pbp");
166 ConvCutList->Add(analysisCuts[i]);
168 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
169 analysisMesonCuts[i] = new AliConversionMesonCuts();
170 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
171 MesonCutList->Add(analysisMesonCuts[i]);
172 analysisMesonCuts[i]->SetFillCutHistograms("");
173 analysisCuts[i]->SetAcceptedHeader(HeaderList);
176 task->SetConversionCutList(numberOfCuts,ConvCutList);
177 task->SetMesonCutList(numberOfCuts,MesonCutList);
178 task->SetMoveParticleAccordingToVertex(kTRUE);
179 task->SetDoMesonAnalysis(kTRUE);
180 // if (trainConfig.Contains("pPb") || trainConfig.Contains("pp") )task->SetDoMesonQA(kTRUE); //Attention new switch for Pi0 QA
181 // if (trainConfig.Contains("pPb") || trainConfig.Contains("pp") )task->SetDoPhotonQA(kTRUE); //Attention new switch small for Photon QA
184 AliAnalysisDataContainer *coutput =
185 mgr->CreateContainer("GammaConvV1", TList::Class(),
186 AliAnalysisManager::kOutputContainer,"GammaConvV1.root");
189 mgr->ConnectInput(task,0,cinput);
190 mgr->ConnectOutput(task,1,coutput);