PWG4/macros/AddTaskCaloTrackCorrM.C

   1 AliAnalysisTaskParticleCorrelationM *AddTaskCaloTrackCorrM(TString data, TString calorimeter, Bool_t kPrintSettings = kFALSE)
   2 {
   3   // Creates a CaloTrackCorr task, configures it and adds it to the analysis manager.
   4
   5   // Get the pointer to the existing analysis manager via the static access method.
   6   //==============================================================================
   7   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   8   if (!mgr) {
   9     ::Error("AddTaskCaloTrackCorr", "No analysis manager to connect to.");
  10     return NULL;
  11   }
  12   Bool_t kUseKinematics = kFALSE;
  13   kUseKinematics = (mgr->GetMCtruthEventHandler())?kTRUE:kFALSE;
  14
  15   cout<<"********* ACCESS KINE? "<<kUseKinematics<<endl;
  16
  17   // Configure analysis
  18   //===========================================================================
  19
  20   //Reader
  21   AliCaloTrackReader * reader = 0x0;
  22   if(data.Contains("AOD")) reader = new AliCaloTrackAODReader();
  23   else if(data=="ESD") reader = new AliCaloTrackESDReader();
  24   else if(data=="MC" || (kUseKinematics && data == "ESD")) reader = new AliCaloTrackMCReader();
  25   reader->SetDebug(-1);//10 for lots of messages
  26   reader->SwitchOnCTS();
  27
  28   if(calorimeter == "EMCAL") {
  29     reader->SwitchOnEMCALCells();
  30     reader->SwitchOnEMCAL();
  31   }
  32   if(calorimeter == "PHOS") {
  33     reader->SwitchOnPHOSCells();
  34     reader->SwitchOnPHOS();
  35   }
  36
  37    reader->SwitchOnSuspiciousClustersRemoval();  //EMCAL
  38
  39   // for case data="deltaAOD", no need to fill the EMCAL/PHOS cluster lists
  40   if(data.Contains("delta")){
  41     reader->SwitchOffEMCAL();
  42     reader->SwitchOffPHOS();
  43     reader->SwitchOffEMCALCells();
  44     reader->SwitchOffPHOSCells();
  45   }
  46
  47   if(kUseKinematics){
  48     if(data == "ESD"){
  49       reader->SwitchOnStack();
  50       reader->SwitchOffAODMCParticles();
  51     }
  52     else if(data == "AOD"){
  53       reader->SwitchOffStack();
  54       reader->SwitchOnAODMCParticles();
  55     }
  56   }
  57
  58   reader->SetZvertexCut(10.);
  59
  60   //Min particle pT
  61   reader->SetEMCALPtMin(0.5);
  62   reader->SetPHOSPtMin(0.5);
  63   reader->SetCTSPtMin(0.1);
  64   reader->SwitchOffWriteDeltaAOD()  ;
  65   if(kPrintSettings) reader->Print("");
  66
  67   // ##### Analysis algorithm settings ####
  68   AliCaloPID * pid = new AliCaloPID();
  69   pid->SetDispersionCut(1.5);
  70   pid->SetTOFCut(5.e-9);
  71   pid->SetDebug(-1);
  72   if(kPrintSettings) pid->Print("");
  73   //Fiducial cut
  74 //  AliFiducialCut * fidCut1stYear = new AliFiducialCut();
  75 //  fidCut1stYear->DoCTSFiducialCut(kTRUE) ;
  76 //  fidCut1stYear->SetSimpleCTSFiducialCut(0.8,0.,360.);
  77 //  fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
  78 //  fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
  79 //  fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
  80 //  fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,300.);
  81   // -------------------------------------------------
  82   // --- Isolation Cut ---
  83   // -------------------------------------------------
  84   AliIsolationCut * ic = new AliIsolationCut();
  85   ic->SetConeSize(0.4);
  86   //ic->SetPtThreshold(0.7);
  87   ic->SetPtFraction(0.1);
  88   ic->SetPtThreshold(0.7) ;
  89   ic->SetSumPtThreshold(1.0) ;
  90   //choose different method for IC:
  91   //kPtThresIC, kSumPtIC, kPtFracIC, kSumPtFracIC
  92   //  ic->SetICMethod(AliIsolationCut::kPtThresIC);
  93   ic->SetICMethod(AliIsolationCut::kSumPtFracIC);
  94   //particle in cone: kNeutralAndCharged=0, kOnlyNeutral=1, kOnlyCharged=2
  95   ic->SetParticleTypeInCone(AliIsolationCut::kOnlyCharged);
  96   if(kPrintSettings) ic->Print("");
  97
  98   //analysis with calorimeter triggers
  99   if(calorimeter=="PHOS" || calorimeter=="EMCAL") {
 100     AliCalorimeterUtils * cu = new AliCalorimeterUtils();
 101     cu->SwitchOnBadChannelsRemoval();
 102
 103     cu->SetNumberOfCellsFromEMCALBorder(1) ; //nEMCAL);
 104     cu->SetNumberOfCellsFromPHOSBorder(2) ; //nPHOS);
 105
 106     // -----------------------------------
 107     // --- Photon and Pi0  Analysis ---
 108     // -----------------------------------
 109
 110     AliAnaPhoton *anaphoton = new AliAnaPhoton();
 111     anaphoton->SetDebug(-1); //10 for lots of messages
 112     anaphoton->SetMinPt(1.0);
 113     anaphoton->SetCaloPID(pid);
 114     anaphoton->SetCalorimeter(calorimeter);
 115     if(kUseKinematics) anaphoton->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 116     else  anaphoton->SwitchOffDataMC() ;
 117     anaphoton->SwitchOffCaloPID();
 118     anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
 119     anaphoton->SwitchOnFiducialCut();
 120     AliFiducialCut * fidCut1stYear = anaphoton->GetFiducialCut();
 121     fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
 122     fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
 123     fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
 124     fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,320.);
 125     anaphoton->SwitchOnTrackMatchRejection();
 126     //settings for different multiplicity analysis
 127     anaphoton->SwitchOffEventSelection() ;
 128     anaphoton->SetMultiplicity(80, 120);
 129
 130     if(calorimeter == "EMCAL"){
 131       anaphoton->SetNCellCut(1);
 132       anaphoton->SetMinPt(0.3);
 133       //if(!kUseKinematics) anaphoton->SetTimeCut(525, 725);
 134       anaphoton->SetMinDistanceToBadChannel(1, 2, 3);
 135     }
 136     else{
 137       anaphoton->SetMinPt(0.3);
 138       anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
 139       anaphoton->SetNCellCut(2);
 140     }
 141     if(!data.Contains("delta")) {
 142       anaphoton->SetOutputAODName(Form("Triggers%s",calorimeter.Data()));
 143       anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
 144     }
 145     else anaphoton->SetOutputAODName(Form("Triggers%s",calorimeter.Data()));
 146     anaphoton->AddToHistogramsName("AnaPhoton_");
 147     //Set Histograms bins and ranges
 148     anaphoton->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 149     //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 150     //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 151
 152     if(kPrintSettings) anaphoton->Print("");
 153     AliAnaPi0 *anapi0 = new AliAnaPi0();
 154     anapi0->SetDebug(-1);//10 for lots of messages
 155     anapi0->SetInputAODName(Form("Triggers%s",calorimeter.Data()));
 156     anapi0->AddToHistogramsName("AnaPi0_");
 157     anapi0->SetCaloPID(pid);
 158     anapi0->SetCalorimeter(calorimeter);
 159     anapi0->SwitchOnFiducialCut();
 160     AliFiducialCut * fidCut1stYear = anapi0->GetFiducialCut();
 161     fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
 162     fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
 163     fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
 164     fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,320.);
 165     anapi0->SetNPID(1);
 166     //settings for different multiplicity analysis
 167     anapi0->SwitchOnEventSelection() ;
 168     anapi0->SetNCentrBin(1);
 169     anapi0->SetMultiplicity(80, 120);
 170     anapi0->SetMultiBin(1);
 171     anapi0->SwitchOffDataMC() ;//Access MC stack and fill more histograms
 172     if(calorimeter=="PHOS") anapi0->SetNumberOfModules(3); //PHOS first year
 173     else if(calorimeter=="EMCAL") anapi0->SetNumberOfModules(4); //EMCAL first year
 174     anapi0->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 175     //anapi0->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 176     //anapi0->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.8, 0.8, 200) ;
 177     anapi0->GetHistogramRanges()->SetHistoMassRangeAndNBins(0., 1.0, 100) ;
 178     anapi0->GetHistogramRanges()->SetHistoAsymmetryRangeAndNBins(0., 1. , 10) ;
 179     if(kPrintSettings) anapi0->Print("");
 180
 181     // -------------------------------------------------
 182     // --- Pi0 EbE Analysis ---
 183     // -------------------------------------------------
 184
 185     AliNeutralMesonSelection *nms = new AliNeutralMesonSelection();
 186     nms->SetInvMassCutRange(0.05, 0.16)     ;
 187     nms->SwitchOffAngleSelection() ;
 188     nms->KeepNeutralMesonSelectionHistos(kTRUE);
 189     //Set Histrograms bins and ranges
 190     nms->SetHistoERangeAndNBins(0, 50, 100) ;
 191     nms->SetHistoPtRangeAndNBins(0, 50, 100) ;
 192     //      nms->GetHistogramRanges()->SetHistoAngleRangeAndNBins(0, 0.3, 100) ;
 193     //      nsm->GetHistogramRanges()->SetHistoIMRangeAndNBins(0, 0.4, 100) ;
 194
 195     AliAnaPi0EbE *anapi0ebe = new AliAnaPi0EbE();
 196     anapi0ebe->SetDebug(-1);//10 for lots of messages
 197     anapi0ebe->SetAnalysisType(AliAnaPi0EbE::kIMCalo);
 198     anapi0ebe->SetMinPt(0.5);
 199     anapi0ebe->SetCalorimeter(calorimeter);
 200     anapi0ebe->SetInputAODName(Form("Triggers%s",calorimeter.Data()));
 201     if(!data.Contains("delta")) {
 202       anapi0ebe->SetOutputAODName(Form("Pi0s%s",calorimeter.Data()));
 203       anapi0ebe->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
 204     }
 205     else  anapi0ebe->SetOutputAODName(Form("Pi0s%s",calorimeter.Data()));
 206
 207     if(kUseKinematics) anapi0ebe->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 208     else  anapi0ebe->SwitchOffDataMC() ;
 209     anapi0ebe->SwitchOffEventSelection() ;
 210     anapi0ebe->SetMultiplicity(80, 120);
 211     anapi0ebe->SetMultiBin(1);
 212     anapi0ebe->SetNeutralMesonSelection(nms);
 213     //Set Histrograms bins and ranges
 214     anapi0ebe->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 215     //      anapi0->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 216     //      anapi0->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 217     if(kPrintSettings) anapi0ebe->Print("");
 218
 219     // ### Pi0 Correlation with hadrons, not isolated
 220     AliAnaParticleHadronCorrelation *anacorrhadronpi0 = new AliAnaParticleHadronCorrelation();
 221     anacorrhadronpi0->SetInputAODName(Form("Pi0s%s",calorimeter.Data()));
 222     anacorrhadronpi0->AddToHistogramsName("AnaHadronCorrPi0_");
 223     anacorrhadronpi0->SetDebug(-1);
 224     anacorrhadronpi0->SwitchOffCaloPID();
 225     anacorrhadronpi0->SwitchOnFiducialCut();
 226     AliFiducialCut * fidCut1stYear = anacorrhadronpi0->GetFiducialCut();
 227     fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
 228     fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
 229     fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
 230     fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,320.);
 231     fidCut1stYear->DoCTSFiducialCut(kTRUE) ;
 232     fidCut1stYear->SetSimpleCTSFiducialCut(0.8,0.,360.);
 233     anacorrhadronpi0->SwitchOnDecayCorr();
 234     anacorrhadronpi0->SetPtCutRange(0.5,50);
 235     anacorrhadronpi0->SetDeltaPhiCutRange(1.5,4.5);
 236     anacorrhadronpi0->SelectIsolated(kFALSE); // do correlation with non isolated pi0
 237     anacorrhadronpi0->SetMultiplicity(80, 100);
 238     anacorrhadronpi0->SetMultiBin(1);
 239     anacorrhadronpi0->SwitchOffNeutralCorr();
 240     anacorrhadronpi0->SwitchOffEventSelection();
 241     anacorrhadronpi0->SwitchOnSeveralUECalculation();
 242     anacorrhadronpi0->SetUeDeltaPhiCutRange(TMath::Pi()/3, 2*TMath::Pi()/3);
 243     if(kUseKinematics) anacorrhadronpi0->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 244     else  anacorrhadronpi0->SwitchOffDataMC() ;
 245     //if(calorimeter=="PHOS"){
 246     //  //Correlate with particles in EMCAL
 247     //  anacorrhadronpi0->SwitchOnCaloPID();
 248     //  anacorrhadronpi0->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
 249     //}
 250     //Set Histograms bins and ranges
 251     anacorrhadronpi0->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 252     //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 253     //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 254     if(kPrintSettings) anacorrhadronpi0->Print("");
 255
 256     AliAnaParticleIsolation *anaisolpi0 = new AliAnaParticleIsolation();
 257     anaisolpi0->SetDebug(-1);
 258     anaisolpi0->SetMinPt(2.);
 259     anaisolpi0->SetInputAODName(Form("Pi0s%s",calorimeter.Data()));
 260     anaisolpi0->AddToHistogramsName("AnaIsolPi0_");
 261     anaisolpi0->SetCalorimeter(calorimeter);
 262     if(kUseKinematics) anaisolpi0->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 263     else  anaisolpi0->SwitchOffDataMC() ;
 264     //Select clusters with no pair, if both clusters with pi0 mass
 265     anaisolpi0->SwitchOffInvariantMass();
 266     //anaisol->SetNeutralMesonSelection(nms);
 267     //Do isolation cut
 268     anaisolpi0->SetIsolationCut(ic);
 269     //Do or not do isolation with previously produced AODs.
 270     //No effect if use of SwitchOnSeveralIsolation()
 271     anaisolpi0->SwitchOffReIsolation();
 272     //Multiple IC
 273     anaisolpi0->SwitchOffSeveralIsolation() ;
 274     //Set Histograms bins and ranges
 275     anaisolpi0->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 276     //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 277     //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 278     if(kPrintSettings) anaisol->Print("");
 279
 280     // ### Pi0 Correlation with hadrons, isolated
 281     AliAnaParticleHadronCorrelation *anacorrhadronisopi0 = new AliAnaParticleHadronCorrelation();
 282     anacorrhadronisopi0->SetInputAODName(Form("Pi0s%s",calorimeter.Data()));
 283     anacorrhadronisopi0->AddToHistogramsName("AnaHadronCorrIsoPi0_");
 284     anacorrhadronisopi0->SetDebug(-1);
 285     anacorrhadronisopi0->SwitchOffCaloPID();
 286     anacorrhadronisopi0->SwitchOnFiducialCut();
 287     AliFiducialCut * fidCut1stYear = anacorrhadronisopi0->GetFiducialCut();
 288     fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
 289     fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
 290     fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
 291     fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,320.);
 292     fidCut1stYear->DoCTSFiducialCut(kTRUE) ;
 293     fidCut1stYear->SetSimpleCTSFiducialCut(0.8,0.,360.);
 294     anacorrhadronisopi0->SwitchOnDecayCorr();
 295     anacorrhadronisopi0->SetPtCutRange(0.5,50);
 296     anacorrhadronisopi0->SetDeltaPhiCutRange(1.5,4.5);
 297     anacorrhadronisopi0->SelectIsolated(kTRUE); // do correlation with isolated pi0
 298     anacorrhadronisopi0->SetMultiplicity(80, 100);
 299     anacorrhadronisopi0->SetMultiBin(1);
 300     anacorrhadronisopi0->SwitchOffNeutralCorr();
 301     anacorrhadronisopi0->SwitchOffEventSelection();
 302     anacorrhadronisopi0->SwitchOnSeveralUECalculation();
 303     anacorrhadronisopi0->SetUeDeltaPhiCutRange(TMath::Pi()/3, 2*TMath::Pi()/3);
 304     if(kUseKinematics) anacorrhadronisopi0->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 305     else  anacorrhadronisopi0->SwitchOffDataMC() ;
 306     //if(calorimeter=="PHOS"){
 307     //  //Correlate with particles in EMCAL
 308     //  anacorrhadronpi0->SwitchOnCaloPID();
 309     //  anacorrhadronpi0->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
 310     //}
 311     //Set Histograms bins and ranges
 312     anacorrhadronisopi0->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 313     //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 314     //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 315     if(kPrintSettings) anacorrhadronisopi0->Print("");
 316
 317   } //analysis in calorimeter
 318   else {
 319     //---charge particle trigger---------------------
 320     //------------------------------------------------
 321
 322     AliAnaChargedParticles *anacharge = new AliAnaChargedParticles();
 323     anacharge->SetDebug(-1); //10 for lots of messages
 324     anacharge->SetMinPt(1.0);
 325     // anacharge->SetCaloPID(pid);
 326     if(kUseKinematics) anacharge->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 327     else  anacharge->SwitchOffDataMC() ;
 328     anacharge->SwitchOffCaloPID();
 329     anacharge->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
 330     anacharge->SwitchOnFiducialCut();
 331     AliFiducialCut * fidCut1stYear = anacharge->GetFiducialCut();
 332     fidCut1stYear->DoCTSFiducialCut(kTRUE) ;
 333     fidCut1stYear->SetSimpleCTSFiducialCut(0.8,0.,360.);
 334     if(!data.Contains("delta")) {
 335       anacharge->SetOutputAODName(Form("Triggers%s",calorimeter.Data()));
 336       anacharge->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
 337     }
 338     else anacharge->SetOutputAODName(Form("Triggers%s",calorimeter.Data()));
 339     anacharge->AddToHistogramsName("AnaCharge_");
 340     //Set Histograms bins and ranges
 341     anacharge->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 342     //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 343     //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 344
 345     if(kPrintSettings) anacharge->Print("");
 346   }
 347
 348   // -------------------------------------------------
 349   // ---  Correlation Analysis with non-isolated triggers ---
 350   // -------------------------------------------------
 351   // ### Correlation with hadrons
 352   AliAnaParticleHadronCorrelation *anacorrhadron = new AliAnaParticleHadronCorrelation();
 353   anacorrhadron->SetInputAODName(Form("Triggers%s",calorimeter.Data()));
 354   anacorrhadron->AddToHistogramsName("AnaHadronCorrTrig_");
 355   anacorrhadron->SetDebug(-1);
 356   anacorrhadron->SwitchOnCaloPID();
 357   anacorrhadron->SwitchOnFiducialCut();
 358   AliFiducialCut * fidCut1stYear = anacorrhadron->GetFiducialCut();
 359   fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
 360   fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
 361   fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
 362   fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,320.);
 363   fidCut1stYear->DoCTSFiducialCut(kTRUE) ;
 364   fidCut1stYear->SetSimpleCTSFiducialCut(0.8,0.,360.);
 365   anacorrhadron->SetPtCutRange(0.5,50);
 366   anacorrhadron->SetDeltaPhiCutRange(1.5,4.5);
 367   anacorrhadron->SetMultiplicity(80, 120);
 368   anacorrhadron->SetMultiBin(1);
 369   anacorrhadron->SwitchOffNeutralCorr();
 370   anacorrhadron->SwitchOffEventSelection();
 371   anacorrhadron->SwitchOnSeveralUECalculation();
 372   anacorrhadron->SetUeDeltaPhiCutRange(TMath::Pi()/3, 2*TMath::Pi()/3);
 373   anacorrhadron->SelectIsolated(kFALSE); // do correlation with isolated photons
 374   if(kUseKinematics) anacorrhadron->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 375   else  anacorrhadron->SwitchOffDataMC() ;
 376   //if(calorimeter=="PHOS"){
 377   //Correlate with particles in EMCAL
 378   //anacorrhadron->SwitchOnCaloPID();
 379   //anacorrhadron->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
 380   //}
 381   //Set Histograms bins and ranges
 382   anacorrhadron->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 383   //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 384   //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 385   if(kPrintSettings) anacorrhadron->Print("");
 386
 387
 388
 389   // ### Isolation analysis ###
 390   AliAnaParticleIsolation *anaisol = new AliAnaParticleIsolation();
 391   anaisol->SetDebug(-1);
 392   anaisol->SetMinPt(2.0);
 393   anaisol->SetInputAODName(Form("Triggers%s",calorimeter.Data()));
 394   anaisol->AddToHistogramsName("AnaIsolTrig_");
 395   anaisol->SetCalorimeter(calorimeter);
 396   if(kUseKinematics) anaisol->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 397   else  anaisol->SwitchOffDataMC() ;
 398   //Select clusters with no pair, if both clusters with pi0 mass
 399   anaisol->SwitchOffInvariantMass();
 400   //anaisol->SetNeutralMesonSelection(nms);
 401   //Do isolation cut
 402   anaisol->SetIsolationCut(ic);
 403   //Do or not do isolation with previously produced AODs.
 404   //No effect if use of SwitchOnSeveralIsolation()
 405   anaisol->SwitchOffReIsolation();
 406   //Multiple IC
 407   anaisol->SwitchOffSeveralIsolation() ;
 408   //Set Histograms bins and ranges
 409   anaisol->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 410   //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 411   //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 412
 413   if(kPrintSettings) anaisol->Print("");
 414
 415   // -------------------------------------------------
 416   // ---  Correlation Analysis with isolated triggers ---
 417   // -------------------------------------------------
 418   // ### Correlation with hadrons
 419   AliAnaParticleHadronCorrelation *anacorrisohadron = new AliAnaParticleHadronCorrelation();
 420   anacorrisohadron->SetInputAODName(Form("Triggers%s",calorimeter.Data()));
 421   anacorrisohadron->AddToHistogramsName("AnaHadronCorrIsoTrig_");
 422   anacorrisohadron->SetDebug(-1);
 423   anacorrisohadron->SwitchOffCaloPID();
 424   anacorrisohadron->SwitchOnFiducialCut();
 425   AliFiducialCut * fidCut1stYear = anacorrisohadron->GetFiducialCut();
 426   fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
 427   fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
 428   fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
 429   fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,320.);
 430   fidCut1stYear->DoCTSFiducialCut(kTRUE) ;
 431   fidCut1stYear->SetSimpleCTSFiducialCut(0.8,0.,360.);
 432   anacorrisohadron->SetPtCutRange(0.5,50);
 433   anacorrisohadron->SetDeltaPhiCutRange(1.5,4.5);
 434   anacorrisohadron->SetMultiplicity(80, 100);
 435   anacorrisohadron->SetMultiBin(1);
 436   anacorrisohadron->SwitchOffNeutralCorr();
 437   anacorrisohadron->SwitchOffEventSelection();
 438   anacorrisohadron->SwitchOnSeveralUECalculation();
 439   anacorrisohadron->SetUeDeltaPhiCutRange(TMath::Pi()/3, 2*TMath::Pi()/3);
 440   anacorrisohadron->SelectIsolated(kTRUE); // do correlation with isolated photons
 441   if(kUseKinematics) anacorrisohadron->SwitchOnDataMC() ;//Access MC stack and fill more histograms
 442   else  anacorrisohadron->SwitchOffDataMC() ;
 443   //if(calorimeter=="PHOS"){
 444   //Correlate with particles in EMCAL
 445   //anacorrhadron->SwitchOnCaloPID();
 446   //anacorrhadron->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
 447   //}
 448   //Set Histograms bins and ranges
 449   anacorrisohadron->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
 450   //      ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
 451   //      ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
 452   if(kPrintSettings) anacorrisohadron->Print("");
 453
 454   // #### Configure Maker ####
 455   AliAnaCaloTrackCorrMaker * maker = new AliAnaCaloTrackCorrMaker();
 456   maker->SetReader(reader);//pointer to reader
 457   Int_t n = 0;//Analysis number, order is important
 458   // Particle selection analysis
 459   if(calorimeter=="PHOS" || calorimeter=="EMCAL") {
 460     maker->SetCaloUtils(cu);  //pointer to calorimeter utils
 461     maker->AddAnalysis(anaphoton,n++);
 462     maker->AddAnalysis(anapi0,n++);
 463     maker->AddAnalysis(anapi0ebe,n++);
 464     maker->AddAnalysis(anacorrhadronpi0,n++);
 465     maker->AddAnalysis(anaisolpi0,n++);
 466     maker->AddAnalysis(anacorrhadronisopi0,n++);
 467   }
 468    if(calorimeter=="CTS")
 469      maker->AddAnalysis(anacharge,n++);
 470   // Correlation analysis
 471   maker->AddAnalysis(anacorrhadron,n++);
 472   // Isolation analysis
 473   maker->AddAnalysis(anaisol,n++);
 474   // Correlation analysis with isolated triggers
 475   maker->AddAnalysis(anacorrisohadron,n);
 476   maker->SetAnaDebug(-1)  ;
 477   maker->SwitchOnHistogramsMaker()  ;
 478   if(data.Contains("delta")) maker->SwitchOffAODsMaker()  ;
 479   else                       maker->SwitchOnAODsMaker()  ;
 480
 481   if(kPrintSettings) maker->Print("");
 482
 483   printf("======================== \n");
 484   printf(" End Configuration of CaloTrackCorr analysis with detector %s \n",calorimeter.Data());
 485   printf("======================== \n");
 486
 487   // Create task
 488   //===========================================================================
 489  // AliAnalysisTaskParticleCorrelationM * task = new AliAnalysisTaskParticleCorrelationM(Form("CaloTrackCorr%s",calorimeter.Data()));
 490   AliAnalysisTaskParticleCorrelationM * task = new AliAnalysisTaskParticleCorrelationM("CaloTrackCorr");
 491   task->SetConfigFileName(""); //Don't configure the analysis via configuration file.
 492   //task->SetDebugLevel(-1);
 493   task->SetAnalysisMaker(maker);
 494   //if(!kSimulation)task->SelectCollisionCandidates(); //AliPhysicsSelection has to be attached before.
 495   mgr->AddTask(task);
 496
 497 //  char name[128];
 498 //  sprintf(name,"CaloTrackCorr_%s",calorimeter.Data());
 499 //  cout<<"Name of task "<<name<<endl;
 500   //AliAnalysisDataContainer *cout_pc = mgr->CreateContainer(Form(name),TList::Class(),
 501   //                                       AliAnalysisManager::kOutputContainer, Form("CaloTrackCorr_%s.root",calorimeter.Data()));
 502
 503   TString outputfile = AliAnalysisManager::GetCommonFileName();
 504   outputfile.ReplaceAll(".root","") ;
 505   outputfile.Append("M.root") ;
 506   //  AliAnalysisDataContainer *cout_pc = mgr->CreateContainer(Form("CaloTrackCorr_%s",calorimeter.Data()),  TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:CaloTrackCorr_%s",outputfile.Data(),calorimeter.Data()));
 507   AliAnalysisDataContainer *cout_pc   = mgr->CreateContainer(calorimeter.Data(), TList::Class(),
 508                                                              AliAnalysisManager::kOutputContainer,
 509                                                              Form("%s:CaloTrackCorr",outputfile.Data()));
 510
 511   AliAnalysisDataContainer *cout_cuts = mgr->CreateContainer(Form("%sCuts",calorimeter.Data()), TList::Class(),
 512                                                              AliAnalysisManager::kParamContainer,
 513                                                              Form("%s:CaloTrackCorrCuts",outputfile.Data()));
 514   // Create ONLY the output containers for the data produced by the task.
 515   // Get and connect other common input/output containers via the manager as below
 516   //==============================================================================
 517   mgr->ConnectInput  (task, 0, mgr->GetCommonInputContainer());
 518   // AOD output slot will be used in a different way in future
 519  //if(!data.Contains("delta")) mgr->ConnectOutput (task, 0, mgr->GetCommonOutputContainer());
 520  // mgr->ConnectOutput (task, 0, mgr->GetCommonOutputContainer());
 521   mgr->ConnectOutput (task, 1, cout_pc);
 522   mgr->ConnectOutput (task, 2, cout_cuts);
 523
 524   return task;
 525 }
 526
 527