AddTaskMultipleTrackCutIsoConeAnalysis.C

Go to the documentation of this file.

// Set includes for compilation

#if !defined(__CINT__) || defined(__MAKECINT__)

// Root
#include <TString.h>
#include <TSystem.h>
#include <TROOT.h>

// Analysis
#include "AliIsolationCut.h"
#include "AliAnalysisTaskCaloTrackCorrelation.h"

// Macros
#include "AddTaskCaloTrackCorrBase.C"
#include "ConfigureCaloTrackCorrAnalysis.C"
#include "GetAlienGlobalProductionVariables.C"

#endif

AliAnalysisTaskCaloTrackCorrelation * AddTaskMultipleTrackCutIsoConeAnalysis
(
 TString  calorimeter   = "EMCAL", // "DCAL", "PHOS"
 Bool_t   simulation    = kFALSE,
 Int_t    year          = -1,
 TString  col           = "",
 TString  period        = "",
 Int_t    rejectEMCTrig = 0,
 TString  clustersArray = "",
 TString  gloCutsString = "",
 Bool_t   calibrate     = kFALSE,
 Bool_t   nonLinOn      = kFALSE,
 TString  analysisString= "Photon_MergedPi0_DecayPi0_Isolation_Correlation_QA_Charged",
 Float_t  shshMax       = 0.27,
 Float_t  isoCone       = 0.4,
 Float_t  isoPtTh       = 2,
 Int_t    isoMethod     = AliIsolationCut::kSumPtIC,
 Int_t    isoContent    = AliIsolationCut::kNeutralAndCharged,
 Int_t    leading       = 0,
 Int_t    minCen        = -1,
 Int_t    maxCen        = -1,
 Bool_t   mixOn         = kTRUE,
 TString  outputfile    = "",
 Bool_t   printSettings = kFALSE,
 Int_t    debug         = 0,  
 const char *trigSuffix = "EMC7"
)
{
  printf("AddTaskMultipleTrackCutIsoConeAnalysis::Start configuration\n");
  
#if defined(__CINT__)
  
  printf("AddTaskMultipleTrackCutIsoConeAnalysis::Load macros\n");
  // Load macros
  //
  // Load macros
  //
  gROOT->LoadMacro("$ALICE_PHYSICS/PWGGA/CaloTrackCorrelations/macros/AddTaskCaloTrackCorrBase.C");
  gROOT->LoadMacro("$ALICE_PHYSICS/PWGGA/CaloTrackCorrelations/macros/ConfigureCaloTrackCorrAnalysis.C");
  gROOT->LoadMacro("$ALICE_PHYSICS/PWGGA/CaloTrackCorrelations/macros/GetAlienGlobalProductionVariables.C");

#endif

  // First check the ALIEN environment settings
  //
  GetAlienGlobalProductionVariables(simulation,col,period,year,kTRUE);
  
  // Init base task
  //
  AliAnalysisTaskCaloTrackCorrelation * task = AddTaskCaloTrackCorrBase
  (calorimeter, simulation, year, col, period, rejectEMCTrig, clustersArray, gloCutsString,
   calibrate, nonLinOn, minCen, maxCen, mixOn, outputfile, printSettings, debug, trigSuffix);
  
  if ( !task ) return NULL;
  
  // No need to continue configuration if event is not processed
  if ( !task->GetAnalysisMaker()->IsEventProcessed() ) return task ;
  
  TList * anaList = task->GetAnalysisMaker()->GetListOfAnalysisContainers();
  printf("TList name: %s\n",anaList->GetName());
  
  // Configure the different analysis
  //
  Float_t rMin[] = {-1,0.05,0.1};
  
  // 3 exclusion of isolation cone options, no cut and 2 cuts
  for(Int_t irmin = 0; irmin < 3; irmin++)
  {
    // Add this string to deacticate the correlation without isolation
    // and other histograms in the Photon and Pi0 selection task
    if(irmin == 1 ) analysisString+="_MultiIso";
    
    // Test 3 track matching options (no track mathing)
    for(Int_t itm = 0; itm < 3; itm++)
    {
      
      TString histoString = Form("TM%d",itm);
      if ( irmin > 0 ) histoString+=Form("_Rmin%1.2f",rMin[irmin]);
      
      
      ConfigureCaloTrackCorrAnalysis
      ( anaList, calorimeter, simulation, year, col, analysisString, histoString, 
       shshMax, isoCone, rMin[irmin], isoPtTh, isoMethod, isoContent,
       leading, itm, mixOn, printSettings, debug);
    }
  }
  
  // Execute some control task only
  ConfigureCaloTrackCorrAnalysis
  ( anaList, calorimeter, simulation, year, col, "QA_Charged", "", 
   -1, -1, -1, -1, -1, -1,-1,-1,0, printSettings, debug);
  
  printf("AddTaskMultipleTrackCutIsoConeAnalysis::End configuration\n");
  
  return task;
}