PWGGA/EMCALTasks/macros/AddTaskSDMGammaMC.C

   1 void AddTaskAliAnalysisTaskGammaMC()
   2 {
   3   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   4   if (!mgr) {
   5     ::Error("AddTaskAliGammaMC", "No analysis manager to connect to.");
   6     return NULL;
   7   }
   8
   9   //#####################################################
  10   // Private Recalibrator:
  11      Int_t GoodTasks  [4] = {1,0,0,0};
  12      Int_t RecalScheme[4] = {0,3,5,7};
  13   //
  14   // 0: no recalibration!
  15   // 1: constant scale factor.
  16   // 2: Symmetric Decay Method
  17   // 3: J's fit to LHC12f1a MC single photons, 4 Aug 2013
  18   // 4: J's fit to the test beam data, 4 Aug 2013
  19   // 5: Based on kSDM/kTBCv3 (for MC)
  20   // 6: kBeamTestCorrectedv2 - in AliROOT!
  21   // 7: kPi0MCv3 - in AliROOT!
  22   //
  23   //#####################################################
  24
  25   AliAnalysisTaskGammaMC* task[4];
  26   AliAnalysisDataContainer*  coutput[4];
  27   char saythis[500];
  28
  29   for(int i=0; i<4; i++){
  30     if(GoodTasks[i]==0)
  31       continue;
  32
  33     sprintf(saythis,"GammaMCTask_%d",i);
  34     task[i] = new AliAnalysisTaskGammaMC(saythis);
  35     task[i]->SelectCollisionCandidates(AliVEvent::kINT7);//LHC13b/c
  36     //task[i]->SelectCollisionCandidates(AliVEvent::kMB);
  37     task[i]->SetRecalScheme(RecalScheme[i]);
  38
  39     task[i]->SetFidPhiMinMax(1.39626, 3.15);// my defaults: 1.39626, 3.15
  40     //task[i]->SetFidPhiMinMax(1.39626, 2.10); //pPb no TRD
  41     //task[i]->SetFidPhiMinMax(2.10, 3.15); //pPb with TRD
  42     task[i]->SetFidEtaMinMax(-0.65, 0.65);
  43
  44     mgr->AddTask(task[i]);
  45     sprintf(saythis,"cont_AliAnalysisTaskGammaMC_%d",RecalScheme[i]);
  46     coutput[i] =
  47       mgr->CreateContainer(saythis,
  48                            TList::Class(),
  49                            AliAnalysisManager::kOutputContainer,
  50                            "GammaMCTask.root");
  51     mgr->ConnectInput (task[i],0,mgr->GetCommonInputContainer());
  52     mgr->ConnectOutput(task[i],1,coutput[i]);
  53
  54     RequestMemory(task[i],1000*1024); // request 1.0GB memory for task
  55
  56     //AliMCEventHandler* handler = new AliMCEventHandler;
  57     //handler->SetReadTR(kFALSE);
  58     //mgr->SetMCtruthEventHandler(handler);
  59   }
  60
  61
  62 }