1 AliAnalysisTask *AddTask_cbaumann_LMEEPbPb2011AOD(Char_t* outputFileName="LMEEoutput.root", Bool_t runAll=kFALSE,Bool_t setMC=kFALSE,Bool_t getFromAlien=kFALSE, Bool_t PIDbaseline=kFALSE, Bool_t rejOnly=kTRUE) {
3 Bool_t bESDANA=kFALSE; //Autodetect via InputHandler
4 //get the current analysis manager
5 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
7 Error("AddTask_cbaumann_LMEEPbPb2011", "No analysis manager found.");
12 // create task and add it to the manager
13 // gSystem->AddIncludePath("$ALICE_ROOT/PWGDQ/dielectron/macrosLMEE");
16 TString configBasePath("$TRAIN_ROOT/cbaumann_dielectron/");
17 TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
18 if (trainRoot.IsNull()) configBasePath= "$ALICE_ROOT/PWGDQ/dielectron/macrosLMEE/";
22 (!gSystem->Exec("alien_cp alien:///alice/cern.ch/user/c/cbaumann/PWGDQ/dielectron/macrosLMEE/ConfigLMEEPbPb2011AOD.C .")) &&
23 (!gSystem->Exec("alien_cp alien:///alice/cern.ch/user/c/cbaumann/PWGDQ/dielectron/macrosLMEE/LMEECutLibAOD.C ."))
25 configBasePath=Form("%s/",gSystem->pwd());
28 TString configFile("ConfigLMEEPbPb2011AOD.C");
29 TString configLMEECutLib("LMEECutLibAOD.C");
31 TString configFilePath(configBasePath+configFile);
32 TString configLMEECutLibPath(configBasePath+configLMEECutLib);
34 //AOD Usage currently tested with separate task, to be merged
35 if (mgr->GetInputEventHandler()->IsA()==AliAODInputHandler::Class()){
36 ::Info("AddTaskLMEEPbPb2011", "no dedicated AOD configuration");
38 else if (mgr->GetInputEventHandler()->IsA()==AliESDInputHandler::Class()){
39 ::Info("AddTaskLMEEPbPb2011AOD","switching on ESD specific code");
44 //Do we have an MC handler?
46 if (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0)
51 if (!gROOT->GetListOfGlobalFunctions()->FindObject(configLMEECutLib.Data()))
52 gROOT->LoadMacro(configLMEECutLibPath.Data());
53 if (!gROOT->GetListOfGlobalFunctions()->FindObject(configFile.Data()))
54 gROOT->LoadMacro(configFilePath.Data());
57 LMEECutLib* cutlib = new LMEECutLib();
58 AliAnalysisTaskMultiDielectron *task=new AliAnalysisTaskMultiDielectron("MultiDiEData");
59 if (!hasMC) task->UsePhysicsSelection();
60 task->SetTriggerMask(AliVEvent::kMB+AliVEvent::kCentral+AliVEvent::kSemiCentral);
61 // task->SelectCollisionCandidates(AliVEvent::kMB+AliVEvent::kCentral+AliVEvent::kSemiCentral);
62 // task->SetRejectPileup();
63 task->SelectCollisionCandidates(AliVEvent::kAny);
64 task->SetEventFilter(cutlib->GetEventCuts(LMEECutLib::kPbPb2011TPCandTOF)); //
67 //load dielectron configuration file
69 //add dielectron analysis with different cuts to the task
73 AliDielectron *lowmass7=ConfigLMEEPbPb2011AOD(7,hasMC,bESDANA);
74 task->AddDielectron(lowmass7);
75 printf("add: %s\n",lowmass7->GetName());
79 AliDielectron *lowmass4=ConfigLMEEPbPb2011AOD(4,hasMC,bESDANA);
80 task->AddDielectron(lowmass4);
81 printf("add: %s\n",lowmass4->GetName());
83 AliDielectron *lowmass1=ConfigLMEEPbPb2011AOD(1,hasMC,bESDANA);
84 task->AddDielectron(lowmass1);
85 printf("add: %s\n",lowmass1->GetName());
88 AliDielectron *lowmass8=ConfigLMEEPbPb2011AOD(8,hasMC,bESDANA);
89 task->AddDielectron(lowmass8);
90 printf("add: %s\n",lowmass8->GetName());
96 //create output container
97 AliAnalysisDataContainer *coutput1 =
98 mgr->CreateContainer("cbaumann_LMEEPbPb2011_tree",
100 AliAnalysisManager::kExchangeContainer,
103 AliAnalysisDataContainer *cOutputHist1 =
104 mgr->CreateContainer("cbaumann_LMEEPbPb2011_out",
106 AliAnalysisManager::kOutputContainer,
109 AliAnalysisDataContainer *cOutputHist2 =
110 mgr->CreateContainer("cbaumann_LMEEPbPb2011_CF",
112 AliAnalysisManager::kOutputContainer,
115 AliAnalysisDataContainer *cOutputHist3 =
116 mgr->CreateContainer("cbaumann_EventStatPbPb2011",
118 AliAnalysisManager::kOutputContainer,
122 mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
123 mgr->ConnectOutput(task, 0, coutput1 );
124 mgr->ConnectOutput(task, 1, cOutputHist1);
125 mgr->ConnectOutput(task, 2, cOutputHist2);
126 mgr->ConnectOutput(task, 3, cOutputHist3);