PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsPA.C

   1 AliAnalysisTaskChargedJetsPA* AddTaskChargedJetsPA(
   2   Double_t            jetRadius               = 0.4,
   3   Int_t               trigger                 = AliVEvent::kINT7,
   4   Bool_t              isMC                    = kFALSE,
   5   Int_t               ptHardBin               = -1,
   6   Double_t            randomConeR             = 0.4,
   7   Double_t            trackBgrdConeR          = 0.6,
   8   const char*         usedTracks              = "PicoTracks",
   9   const char*         centralityType          = "V0A",
  10   Double_t            trackEtaWindow          = 0.9,
  11   Double_t            vertexWindow            = 10.0,
  12   Double_t            vertexMaxR              = 1.0,
  13   Double_t            minJetPt                = 5.0, // signal jet min pt
  14   Double_t            minBackgroundJetPt      = 0.0, // background jet min pt
  15   Double_t            dijetLeadingMinPt       = 10.0,
  16   Double_t            dijetMaxAngleDev        = 10.0,
  17   Int_t               numberOfPtHardBins      = 0,
  18   const char*         fileEtaCorrectionFactors= "alien:///alice/cern.ch/user/r/rhaake/pA/EtaCorrectionFactors.root",
  19   const char*         externalMacro           = NULL,
  20   Bool_t              isEMCalTrain            = kTRUE
  21 )
  22 {
  23   // #### Detect the demanded trigger with its readable name
  24   TString triggerName(Form("Trigger_%i", trigger));
  25   if (trigger == AliVEvent::kAnyINT)
  26     triggerName = "kAnyINT";
  27   else if (trigger == AliVEvent::kAny)
  28     triggerName = "kAny";
  29   else if(trigger == AliVEvent::kINT7)
  30     triggerName = "kINT7";
  31   else if(trigger == AliVEvent::kMB)
  32     triggerName = "kMB";
  33   else if(trigger == AliVEvent::kEMC7)
  34     triggerName = "kEMC7";
  35   else if(trigger == AliVEvent::kEMCEJE)
  36     triggerName = "kEMCEJE";
  37   else if(trigger == AliVEvent::kEMCEGA)
  38     triggerName = "kEMCEGA";
  39
  40   // #### Define manager and data container names
  41   AliAnalysisManager *manager = AliAnalysisManager::GetAnalysisManager();
  42   if (!manager) {
  43     ::Error("AddTaskChargedJetsPA", "No analysis manager to connect to.");
  44     return NULL;
  45   }
  46
  47   TString stringPtHard("");
  48   if (ptHardBin!=-1)
  49     stringPtHard = Form("_PtHard_%d",ptHardBin);
  50   TString myContName("");
  51   if(isMC)
  52     myContName = Form("AnalysisR0%2.0f_%s_MC%s",jetRadius*100,triggerName.Data(), stringPtHard.Data());
  53   else
  54     myContName = Form("AnalysisR0%2.0f_%s%s",jetRadius*100,triggerName.Data(), stringPtHard.Data());
  55
  56   // #### Add necessary jet finder tasks
  57   gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
  58   AliEmcalJetTask* jetFinderTask = AddTaskEmcalJet(usedTracks,"",1,jetRadius,1,0.150,0.300); // anti-kt
  59   AliEmcalJetTask* jetFinderTaskKT = AddTaskEmcalJet(usedTracks,"",0,jetRadius,1,0.150,0.300); // kt
  60
  61   jetFinderTaskKT->SetMinJetPt(minBackgroundJetPt);
  62   // #### Define analysis task
  63   AliAnalysisTaskChargedJetsPA *task = NULL;
  64   contHistos = manager->CreateContainer(myContName.Data(), TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:ChargedJetsPA", AliAnalysisManager::GetCommonFileName()));
  65   task = new AliAnalysisTaskChargedJetsPA(Form("AnalysisPA_%s_%s", jetFinderTask->GetName(), triggerName.Data()), usedTracks, jetFinderTask->GetName(),jetFinderTaskKT->GetName());
  66
  67   // #### Task preferences
  68   task->SetAcceptanceWindows(trackEtaWindow, vertexWindow, vertexMaxR, jetRadius, jetRadius);
  69   task->SetAnalyzeQA(kTRUE);
  70   task->SetAnalyzeBackground(kTRUE);
  71   task->SetAnalyzeDeprecatedBackgrounds(kTRUE);
  72
  73   task->SetSignalJetMinPt(minJetPt);
  74   task->SetSignalJetMinArea(0.6*jetRadius*jetRadius*TMath::Pi());
  75   task->SetDijetLeadingMinPt(dijetLeadingMinPt);
  76   task->SetDijetMaxAngleDeviation(dijetMaxAngleDev);
  77   task->SetRandConeRadius(randomConeR);
  78   task->SetBackgroundJetMinPt(minBackgroundJetPt);
  79   task->SetTRBackgroundConeRadius(trackBgrdConeR);
  80   task->SelectCollisionCandidates(trigger);
  81   task->SetCentralityType(centralityType);
  82   task->SetUsePtHardBin(ptHardBin);
  83
  84   if(numberOfPtHardBins)
  85     task->SetNumberOfPtHardBins(numberOfPtHardBins);
  86
  87   // #### Add analysis task
  88   manager->AddTask(task);
  89   manager->ConnectInput(task, 0, manager->GetCommonInputContainer());
  90   manager->ConnectOutput(task, 1, contHistos);
  91
  92   if(isEMCalTrain)
  93     RequestMemory(task,200*1024);
  94
  95   // #### Do some nasty piggybacking on demand
  96   if (externalMacro)
  97     gROOT->Macro(externalMacro);
  98
  99   return task;
 100 }