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;
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 = 2;
115 else numberOfCuts = 3;
117 TString *cutarray = new TString[numberOfCuts];
118 TString *mesonCutArray = new TString[numberOfCuts];
120 if(trainConfig.Contains("PbPb")){
121 cutarray[ 0] = "1000003042092970723220000"; mesonCutArray[ 0] = "01022045000"; // all centralities
122 cutarray[ 1] = "3010003042092970723220000"; mesonCutArray[ 1] = "01022045000"; // most central
123 cutarray[ 2] = "1680003042092970723220000"; mesonCutArray[ 2] = "01022065000"; // peripheral
124 } else if(trainConfig.Contains("pPb")){ //pA needs thighter rapidity cut y < 0.5
125 cutarray[ 0] = "8000000082093172023290000"; mesonCutArray[0] = "01629045000"; //standard cut Pi0 Pb 00-100 shifted Eta 0.4
126 cutarray[ 1] = "8000000072093172023290000"; mesonCutArray[1] = "01627045000"; //standard cut Pi0 Pb 00-100 wo shifted Eta 0.3
128 cutarray[ 0] = "0000011002093663003800000"; mesonCutArray[0] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
129 cutarray[ 1] = "0001011002093663003800000"; mesonCutArray[1] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
130 cutarray[ 2] = "0000011002093260003800000"; mesonCutArray[2] = "01631031009"; //standard cut Gamma pp 2-76TeV
134 TList *ConvCutList = new TList();
135 TList *MesonCutList = new TList();
137 TList *HeaderList = new TList();
138 // TObjString *Header1 = new TObjString("PARAM");
139 // HeaderList->Add(Header1);
140 TObjString *Header2 = new TObjString("BOX");
141 HeaderList->Add(Header2);
142 // TObjString *Header3 = new TObjString("Pythia");
143 // HeaderList->Add(Header3);
144 // TObjString *Header4 = new TObjString("Hijing");
145 // HeaderList->Add(Header4);
148 ConvCutList->SetOwner(kTRUE);
149 AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts];
150 MesonCutList->SetOwner(kTRUE);
151 AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
154 for(Int_t i = 0; i<numberOfCuts; i++){
155 analysisCuts[i] = new AliConversionCuts();
156 analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
157 if (trainConfig.Contains("pPb")){
158 analysisCuts[i]->SelectCollisionCandidates(AliVEvent::kINT7);
160 analysisCuts[i]->DoEtaShift(kTRUE);
161 analysisCuts[i]->SetEtaShift("Pbp");
164 ConvCutList->Add(analysisCuts[i]);
166 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
167 analysisMesonCuts[i] = new AliConversionMesonCuts();
168 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
169 MesonCutList->Add(analysisMesonCuts[i]);
170 analysisMesonCuts[i]->SetFillCutHistograms("");
171 analysisCuts[i]->SetAcceptedHeader(HeaderList);
175 task->SetConversionCutList(numberOfCuts,ConvCutList);
176 task->SetMesonCutList(numberOfCuts,MesonCutList);
177 task->SetMoveParticleAccordingToVertex(kTRUE);
178 task->SetDoMesonAnalysis(kTRUE);
179 if (trainConfig.Contains("pPb") || trainConfig.Contains("pp") )task->SetDoMesonQA(kTRUE); //Attention new switch for Pi0 QA
180 if (trainConfig.Contains("pPb") || trainConfig.Contains("pp") )task->SetDoPhotonQA(kTRUE); //Attention new switch small for Photon QA
183 AliAnalysisDataContainer *coutput =
184 mgr->CreateContainer("GammaConvV1", TList::Class(),
185 AliAnalysisManager::kOutputContainer,"GammaConvV1.root");
188 mgr->ConnectInput(task,0,cinput);
189 mgr->ConnectOutput(task,1,coutput);