ConfigureCaloTrackCorrAnalysis.C

Go to the documentation of this file.

// Set includes for compilation

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

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

#include "AliAnaPhoton.h"
#include "AliAnaPi0.h"
#include "AliAnaPi0EbE.h"
#include "AliHistogramRanges.h"
#include "AliAnaParticleIsolation.h"
#include "AliAnaParticleHadronCorrelation.h"
#include "AliAnaChargedParticles.h"
#include "AliAnaCalorimeterQA.h"
#include "AliAnaClusterShapeCorrelStudies.h"
#include "AliAnaGeneratorKine.h"
#include "AliAnalysisTaskCaloTrackCorrelation.h"
#include "AliAnaCaloTrackCorrMaker.h"
#endif

TString kAnaCaloTrackCorr = "";

// Some examples of strings: "Photon_MergedPi0_DecayPi0_Isolation_Correlation_Bkg_QA_Charged_HighMult_MultiIso_PerSM_PerTCard",
TString kAnaCutsString = ""; 

void SetAnalysisCommonParameters(AliAnaCaloTrackCorrBaseClass* ana, TString histoString,
                                 TString calorimeter,   Int_t  year,
                                 TString col,           Bool_t simulation,
                                 Bool_t  printSettings, Int_t  debug)
{
  //
  // Histograms ranges
  //
  AliHistogramRanges* histoRanges = ana->GetHistogramRanges();
  
  if ( histoString != "" ) 
    ana->AddToHistogramsName(Form("%s_%s",  histoString.Data(), (ana->GetAddedHistogramsStringToName()).Data()) );

  histoRanges->SetHistoPtRangeAndNBins(0, 100, 200) ; // Energy and pt histograms
  
  if(calorimeter=="EMCAL")
  {
    ana->SetFirstSMCoveredByTRD(-1);
    
    if ( year == 2010 )
    {
      histoRanges->SetHistoPhiRangeAndNBins(78*TMath::DegToRad(), 122*TMath::DegToRad(), 78) ;
      histoRanges->SetHistoXRangeAndNBins(-230,90,120); // QA
      histoRanges->SetHistoYRangeAndNBins(370,450,40);  // QA
    }
    else if ( year < 2014 )
    {
      histoRanges->SetHistoPhiRangeAndNBins(78*TMath::DegToRad(), 182*TMath::DegToRad(), 104) ;
      histoRanges->SetHistoXRangeAndNBins(-460,90,200); // QA
      histoRanges->SetHistoYRangeAndNBins(100,450,100); // QA  
      
      if     (year == 2011) ana->SetFirstSMCoveredByTRD( 6);
      else if(year == 2012 ||
              year == 2013) ana->SetFirstSMCoveredByTRD( 4);
    }
    else // Run2
    {
      histoRanges->SetHistoPhiRangeAndNBins(78*TMath::DegToRad(), 189*TMath::DegToRad(), 111) ;
      histoRanges->SetHistoXRangeAndNBins(-460,460,230); // QA, revise
      histoRanges->SetHistoYRangeAndNBins(-450,450,225); // QA, revise
    }
    
    histoRanges->SetHistoEtaRangeAndNBins(-0.72, 0.72, 144) ;
  }
  else if(calorimeter=="DCAL")
  {
    histoRanges->SetHistoPhiRangeAndNBins(260*TMath::DegToRad(), 327*TMath::DegToRad(), 67) ;
    histoRanges->SetHistoEtaRangeAndNBins(-0.72, 0.72, 144) ;
    histoRanges->SetHistoXRangeAndNBins(-460,460,230); // QA, revise
    histoRanges->SetHistoYRangeAndNBins(-450,450,225); // QA, revise
  } 
  else if(calorimeter=="PHOS")
  {
    histoRanges->SetHistoPhiRangeAndNBins(250*TMath::DegToRad(), 320*TMath::DegToRad(), 70) ;
    histoRanges->SetHistoEtaRangeAndNBins(-0.13, 0.13, 130) ;
  }
  else if(calorimeter=="CTS")
  {
    ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 200) ;
    ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-1.5, 1.5, 300) ;
  }
  
  histoRanges->SetHistoShowerShapeRangeAndNBins(-0.1, 2.9, 300);
  
  // Invariant mass histo
  histoRanges->SetHistoMassRangeAndNBins(0., 1., 200) ;
  histoRanges->SetHistoAsymmetryRangeAndNBins(0., 1. , 100) ;
  
  // check if time calibration is on
  //histoRanges->SetHistoTimeRangeAndNBins(-1000.,1000,1000);
  histoRanges->SetHistoTimeRangeAndNBins(-400.,400,400);
  histoRanges->SetHistoDiffTimeRangeAndNBins(-200, 200, 800);
  
  // track-cluster residuals
  histoRanges->SetHistoTrackResidualEtaRangeAndNBins(-0.06,0.06,120);
  histoRanges->SetHistoTrackResidualPhiRangeAndNBins(-0.06,0.06,120);
  histoRanges->SetHistodRRangeAndNBins(0.,0.06,60);//QA
  
  // QA, electron, charged
  histoRanges->SetHistoPOverERangeAndNBins(0,1.5,150);
  histoRanges->SetHistodEdxRangeAndNBins(0.,200.,200);
  
  // QA
  histoRanges->SetHistoFinePtRangeAndNBins(0, 10, 200) ; // bining for fhAmpId
  histoRanges->SetHistoVertexDistRangeAndNBins(0.,500.,500);
  histoRanges->SetHistoZRangeAndNBins(-350,350,175);
  histoRanges->SetHistoRRangeAndNBins(430,460,30);
  histoRanges->SetHistoV0SignalRangeAndNBins(0,5000,500);
  histoRanges->SetHistoV0MultiplicityRangeAndNBins(0,5000,500);
  
  // QA, correlation
  if(col=="PbPb")
  {
    histoRanges->SetHistoNClusterCellRangeAndNBins(0,100,100);
    histoRanges->SetHistoNClustersRangeAndNBins(0,500,50);
    histoRanges->SetHistoTrackMultiplicityRangeAndNBins(0,2000,200);
  }
  else
  {
    histoRanges->SetHistoNClusterCellRangeAndNBins(0,50,50);
    histoRanges->SetHistoNClustersRangeAndNBins(0,50,50);
    histoRanges->SetHistoTrackMultiplicityRangeAndNBins(0,200,200);
  }
  
  // xE, zT
  histoRanges->SetHistoRatioRangeAndNBins(0.,2.,200);
  histoRanges->SetHistoHBPRangeAndNBins  (0.,10.,200);
  
  // Isolation
  histoRanges->SetHistoPtInConeRangeAndNBins(0, 50 , 250);
  histoRanges->SetHistoPtSumRangeAndNBins   (0, 100, 250);
  
  //
  // MC histograms?
  //
  if(simulation) ana->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.
  else           ana->SwitchOffDataMC() ;
  
 
  //
  // Specialized histograms on multiplicity
  //
  if ( kAnaCutsString.Contains("HighMult") ) 
    ana->SwitchOnFillHighMultiplicityHistograms();
  else
    ana->SwitchOnFillHighMultiplicityHistograms();
  
  //
  // Debug
  //
  if(printSettings) ana->Print("");
  
  ana->SetDebug(debug); // 10 for lots of messages
}

AliAnaPhoton* ConfigurePhotonAnalysis(TString col,           Bool_t simulation,
                                      TString calorimeter,   Int_t year, Int_t tm,
                                      Bool_t  printSettings, Int_t   debug, 
                                      TString histoString )
{
  AliAnaPhoton *ana = new AliAnaPhoton();
  
  // cluster selection cuts
  
  ana->SwitchOffRealCaloAcceptance();
  
  ana->SwitchOnFiducialCut(); 
  if      ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.7,  80, 187) ; // EMC 
  else if ( calorimeter == "DCAL"  ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.7, 260, 327) ; // DMC  
  
  if ( calorimeter.Contains("CAL") ) ana->GetFiducialCut()->DoEMCALFiducialCut(kTRUE);  
  
  ana->SetCalorimeter(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetCalorimeter(calo);
  }
  
  ana->SwitchOnFillShowerShapeHistograms();  // Filled before photon shower shape selection
  if( kAnaCutsString.Contains("MultiIso") ) ana->SwitchOffFillShowerShapeHistograms();
    
  if ( kAnaCutsString.Contains("PerSM") ) 
    ana->SwitchOnFillShowerShapeHistogramsPerSM(); 
  
  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
  
  if(tm) ana->SwitchOnTrackMatchRejection() ;
  else   ana->SwitchOffTrackMatchRejection() ;
  
  // Fill track matching histograms if track matching activated
  // and only once in case of multiple analysis
  if ( tm && !kAnaCutsString.Contains("MultiIso") ) 
  {
    ana->SwitchOnTMHistoFill() ;
    if(tm > 1) ana->SwitchOnTMTrackPtHistoFill() ;
  }
  
  if(calorimeter == "PHOS")
  {
    ana->SetNCellCut(2);// At least 3 cells
    ana->SetMinPt(0.3);
    ana->SetMinDistanceToBadChannel(2, 4, 5);
    ana->SetTimeCut(-1e10,1e10); // open cut
  }
  else
  {//EMCAL
    ana->SetNCellCut(1);// At least 2 cells
    ana->SetMinEnergy(0.3); // avoid mip peak at E = 260 MeV
    ana->SetMaxEnergy(100);
    ana->SetTimeCut(-1e10,1e10); // open cut, usual time window of [425-825] ns if time recalibration is off
                                 // restrict to less than 100 ns when time calibration is on
    ana->SetMinDistanceToBadChannel(2, 4, 6);
    
    // NLM cut, used in all, exclude clusters with more than 2 maxima
    // Not needed if M02 cut is already strong or clusterizer V2
    ana->SetNLMCut(1, 2) ;
  }
  
  //PID cuts (shower shape)
  ana->SwitchOnCaloPID(); // do PID selection, unless specified in GetCaloPID, selection not based on bayesian
  AliCaloPID* caloPID = ana->GetCaloPID();
  //Not used in bayesian
  
  // EMCAL
  
  //caloPID->SetEMCALLambda0CutMax(0.27);
  caloPID->SetEMCALLambda0CutMax(10); // open, full shower shape needed for isolation studies
  caloPID->SetEMCALLambda0CutMin(0.10);
  
  // Track matching
  // tm = 1, fixed cuts
  caloPID->SetEMCALDEtaCut(0.025);
  caloPID->SetEMCALDPhiCut(0.030);
  
  // pT track dependent cuts
  if(tm > 1) caloPID->SwitchOnEMCTrackPtDepResMatching();
  
  // PHOS
  caloPID->SetPHOSDispersionCut(2.5);
  caloPID->SetPHOSRCut(2.);
  //if(kInputData=="AOD") caloPID->SetPHOSRCut(2000.); // Open cut since dX, dZ not stored
  
  // Branch AOD settings
  ana->SetOutputAODName(Form("PhotonTrigger_%s",kAnaCaloTrackCorr.Data()));
  ana->SetOutputAODClassName("AliCaloTrackParticleCorrelation");
  
  //Set Histograms name tag, bins and ranges
  //ana->AddToHistogramsName(Form("AnaPhoton_TM%d_",tm));
  ana->AddToHistogramsName("AnaPhoton_");
  
  SetAnalysisCommonParameters(ana,histoString,calorimeter,year,col,simulation,printSettings,debug) ; // see method below
  
  if(ana->GetFirstSMCoveredByTRD() > 0)
    printf("ConfigurePhotonAnalysis() >>> Set first SM covered by TRD, SM=%d <<< year %d \n", ana->GetFirstSMCoveredByTRD(),year);
  
  // Number of particle type MC histograms
  ana->FillNOriginHistograms (17); // 18 max
  ana->FillNPrimaryHistograms(6); // 6 max

  return ana;
}

AliAnaPi0EbE* ConfigurePi0EbEAnalysis(TString particle,      Int_t  analysis,
                                      Bool_t useSSIso,       Bool_t useAsy,
                                      TString col,           Bool_t simulation,
                                      TString calorimeter,   Int_t  year,  Int_t tm,
                                      Bool_t  printSettings, Int_t debug, 
                                      TString histoString )
{
  // Configuration of pi0 event by event selection
  
  AliAnaPi0EbE *ana = new AliAnaPi0EbE();
  
  ana->SetDebug(debug);
  
  ana->SetAnalysisType((AliAnaPi0EbE::anaTypes)analysis);
  TString opt = "";
  if(analysis==AliAnaPi0EbE::kIMCaloTracks) opt = "Conv";
  if(analysis==AliAnaPi0EbE::kSSCalo)       opt = "SS";
  
  ana->SwitchOffAllNLMHistoFill();
  ana->SwitchOffSelectedClusterHistoFill();
  
  ana->SwitchOffFillWeightHistograms();
  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
  //if(!kTime && !simulation) ana->SwitchOnFillEMCALBCHistograms();
  
  if(tm) ana->SwitchOnTrackMatchRejection() ;
  else   ana->SwitchOffTrackMatchRejection() ;
  
  if ( tm && !kAnaCutsString.Contains("Photon") &&  !kAnaCutsString.Contains("MultiIso") )
    ana->SwitchOnTMHistoFill() ;
  else
    ana->SwitchOffTMHistoFill() ;
    
  ana->SetCalorimeter(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetCalorimeter(calo);
  }
  
  // Branch AOD settings
  ana->SetOutputAODName(Form("%s%sTrigger_%s",particle.Data(), opt.Data(), kAnaCaloTrackCorr.Data()));
  printf("ConfigurePi0EbEAnalysis() *** Out branch %s***\n",ana->GetOutputAODName().Data());
  ana->SetOutputAODClassName("AliCaloTrackParticleCorrelation");
  
  if(analysis == AliAnaPi0EbE::kIMCaloTracks) ana->SetInputAODGammaConvName("PhotonsCTS");
  
  //Set Histograms name tag, bins and ranges
  
  //ana->AddToHistogramsName(Form("Ana%s%sEbE_TM%d_",particle.Data(),opt.Data(),tm));
  ana->AddToHistogramsName(Form("Ana%s%sEbE_",particle.Data(),opt.Data()));
  
  SetAnalysisCommonParameters(ana,histoString,calorimeter,year,col,simulation,printSettings,debug); // see method below
  
  if(analysis!=AliAnaPi0EbE::kSSCalo)
  {
    ana->SetInputAODName(Form("PhotonTrigger_%s",kAnaCaloTrackCorr.Data()));
    
    ana->SetM02CutForInvMass(0.1,0.35); // Loose SS cut
    
    ana->SwitchOnSelectPairInIsolationCone();
    ana->SetR(0.4);
    ana->SetIsolationCandidateMinPt(5);
    
    if(useSSIso)
    {
      ana->SwitchOnSelectIsolatedDecay();
      //ana->AddToHistogramsName(Form("Ana%s%sEbEIsoDecay_TM%d_",particle.Data(),opt.Data(),tm));
      ana->AddToHistogramsName(Form("Ana%s%sEbEIsoDecay_",particle.Data(),opt.Data()));
      ana->SetOutputAODName(Form("%s%sIsoDecayTrigger_%s",particle.Data(), opt.Data(), kAnaCaloTrackCorr.Data()));
    }
    
    if ( calorimeter.Contains("CAL") && !simulation ) ana->SetPairTimeCut(100);
    
    AliNeutralMesonSelection *nms = ana->GetNeutralMesonSelection();
    nms->SetParticle(particle);
    
    //****
    nms->SetInvMassCutMaxParameters(0,0,0); // Overrule the setting in SetParticle for Pi0 option
                                            //****
    
    // Tighten a bit mass cut with respect to default window
    if(particle=="Pi0") nms->SetInvMassCutRange(0.110,0.160);
    if(particle=="Eta") nms->SetInvMassCutRange(0.520,0.580);
    
    //if(!particle.Contains("SideBand")) nms->SwitchOnAngleSelection();
    //else nms->SwitchOnAngleSelection();
    
    nms->SwitchOffAngleSelection();
    
    if(particle.Contains("Pi0SideBand")) // For pi0, do not consider left band
      nms->SetSideBandCutRanges(-1,0,0.190,0.240);
    
    if(particle.Contains("EtaSideBand")) // For pi0, do not consider left band
      nms->SetSideBandCutRanges(0.410,0.470,0.620,0.680);
    
    nms->KeepNeutralMesonSelectionHistos(kTRUE);
    //nms->SetAngleMaxParam(2,0.2);
    nms->SetHistoERangeAndNBins(0, 20, 80) ;
    //nms->SetHistoIMRangeAndNBins(0, 1, 400);
  }
  else  
  {
    // cluster splitting settings
    ana->SetMinEnergy(6);
    ana->SetMaxEnergy(100.);
    
    ana->SetNLMMinEnergy(0, 10);
    ana->SetNLMMinEnergy(1, 6);
    ana->SetNLMMinEnergy(2, 6);
    
    // NLM cut, used in all, exclude clusters with more than 2 maxima
    ana->SetNLMCut(1, 2) ;
    
    //
    ana->SetMinDistanceToBadChannel(2, 4, 6);
    ana->SwitchOnSplitClusterDistToBad();
    ana->SetTimeCut(-1e10,1e10); // Open time cut
    
    AliCaloPID* caloPID = ana->GetCaloPID();
    
    // Track matching
    // tm = 1, fixed cuts
    caloPID->SetEMCALDEtaCut(0.025);
    caloPID->SetEMCALDPhiCut(0.030);
    
    // pT track dependent cuts
    if(tm > 1) caloPID->SwitchOnEMCTrackPtDepResMatching();
    
    caloPID->SetSplitWidthSigma(3); // cut at 3 sigma of the mean pi0 peak.
    
    if(!useSSIso)
    {
      printf("ConfigurePi0EbEAnalysis() << Do not apply SS cut on merged pi0 analysis >> \n");
      caloPID->SwitchOffSplitShowerShapeCut() ;
      //ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenSS_TM%d_",particle.Data(),opt.Data(),tm));
      ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenSS_",particle.Data(),opt.Data()));
      ana->SetOutputAODName(Form("%s%sTrigger_%s_OpenSS",particle.Data(), opt.Data(), kAnaCaloTrackCorr.Data()));
      caloPID->SetClusterSplittingM02Cut(0.1,10);
    }
    else
    {
      caloPID->SetClusterSplittingM02Cut(0.3,4); // Do the selection in the analysis class and not in the PID method to fill SS histograms
      caloPID->SwitchOnSplitShowerShapeCut() ;
    }
    
    if(useAsy)
    {
      caloPID->SwitchOnSplitAsymmetryCut() ;
      ana->GetCaloPID()->SetSubClusterEnergyMinimum(0,2);
      ana->GetCaloPID()->SetSubClusterEnergyMinimum(1,0.5);
      ana->GetCaloPID()->SetSubClusterEnergyMinimum(2,0.5);
    }
    else
    {
      caloPID->SwitchOffSplitAsymmetryCut() ;
      if(!useSSIso)
      {
        //ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenSS_OpenAsy_TM%d_",particle.Data(),opt.Data(),tm));
        ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenSS_OpenAsy_",particle.Data(),opt.Data()));
        ana->SetOutputAODName(Form("%s%sTrigger_%s_OpenSS_OpenAsy",particle.Data(), opt.Data(), kAnaCaloTrackCorr.Data()));
      }
      else
      {
        //ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenAsy_TM%d_",particle.Data(),opt.Data(),tm));
        ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenAsy_",particle.Data(),opt.Data()));
        ana->SetOutputAODName(Form("%s%sTrigger_%s_OpenAsy",particle.Data(), opt.Data(), kAnaCaloTrackCorr.Data()));
      }
    }
    
    // For Pi0 only if  SwitchOnSimpleSplitMassCut()
    caloPID->SetPi0MassRange(0.10, 0.18);
    caloPID->SetEtaMassRange(0.50, 0.60);
    caloPID->SetPhotonMassRange(0.00, 0.08);
    
    caloPID->SetClusterSplittingMinNCells(6);
    
    //caloPID->SetSplitEnergyFractionMinimum(0, 0.95);
    //caloPID->SetSplitEnergyFractionMinimum(1, 0.95);
    //caloPID->SetSplitEnergyFractionMinimum(2, 0.8);
    
    if(col=="PbPb" || kAnaCaloTrackCorr.Contains("150"))
    {
      caloPID->SetClusterSplittingMinNCells(4);
      //caloPID->SetPi0MassShiftHighECell(0.005);
    }
  }
  
  return  ana;
}


AliAnaPi0* ConfigureInvariantMassAnalysis
(TString col,           Bool_t simulation,
 TString calorimeter,   Bool_t bothCalo,   
 Int_t   year,          Bool_t mixOn,
 Bool_t  printSettings, Int_t  debug, 
 TString histoString )
{
  AliAnaPi0 *ana = new AliAnaPi0();
  
  ana->SetDebug(debug);
  
  // Input delta AOD settings
  ana->SetInputAODName(Form("PhotonTrigger_%s",kAnaCaloTrackCorr.Data()));

  if(bothCalo)
  {
    ana->SwitchOnPairWithOtherDetector();
    
    TString otherDetector = Form("PhotonTrigger_%s",kAnaCaloTrackCorr.Data());
    
    if ( calorimeter == "EMCAL" ) 
      otherDetector.ReplaceAll("PhotonTrigger_EMCAL","PhotonTrigger_PHOS");
    else    
      otherDetector.ReplaceAll("PhotonTrigger_PHOS","PhotonTrigger_EMCAL");
    
    ana->SetOtherDetectorInputName(otherDetector);
  }
  
  // Calorimeter settings
  ana->SetCalorimeter(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetCalorimeter(calo);
  }
  
  // Acceptance plots
  ana->SwitchOnFiducialCut(); 
  if      ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.7,  80, 187) ; // EMC 
  else if ( calorimeter == "DCAL"  ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.7, 260, 327) ; // DMC  
  
  if ( calorimeter.Contains("CAL") ) ana->GetFiducialCut()->DoEMCALFiducialCut(kTRUE);  

  ana->SwitchOffRealCaloAcceptance();
  
  // settings for pp collision mixing
  if(mixOn) ana->SwitchOnOwnMix();
  else      ana->SwitchOffOwnMix();
  
  // Cuts
  if (calorimeter == "EMCAL" ) 
  {
    if(year < 2014) ana->SetPairTimeCut(50);
    else            ana->SetPairTimeCut(200); // REMEMBER to remove this when time calib is on
  }
  
  ana->SetNPIDBits(1);
  ana->SetNAsymCuts(1); // no asymmetry cut, previous studies showed small effect.
                        // In EMCAL assymetry cut prevents combination of assymetric decays which is the main source of pi0 at high E.

  if     (col == "pp"  )
  {
    printf("ConfigureInvariantMassAnalysis() - Set pp configuration\n");
    ana->SetNCentrBin(1);
    ana->SwitchOffTrackMultBins();
    ana->SetNZvertBin(10);
    ana->SetNRPBin(1);
    ana->SetNMaxEvMix(100);
    ana->SetMinPt(0.7);
  }
  else if(col == "PbPb")
  {
    printf("ConfigureInvariantMassAnalysis() - Set PbPb configuration\n");
    ana->SetNCentrBin(10);
    ana->SetNZvertBin(10);
    ana->SetNRPBin(4);
    ana->SetNMaxEvMix(10);
    ana->SetMinPt(1.5);
  }
  else if(col =="pPb")
  {
    printf("ConfigureInvariantMassAnalysis() - Set pPb configuration\n");
    ana->SetNCentrBin(1);
    ana->SetNZvertBin(10);
    ana->SetNRPBin(4);
    ana->SetNMaxEvMix(100);
    ana->SetMinPt(0.7);
  }
  
  // Angle cut, avoid pairs with too large angle
  ana->SwitchOnAngleSelection(); 
  ana->SetAngleMaxCut(TMath::DegToRad()*80.); // EMCal: 4 SM in phi, 2 full SMs in eta
  ana->SetAngleCut(0.017); // Minimum angle open, cell size
  
  if ( kAnaCutsString.Contains("PerSM") ) //&& !bothCalo ) 
    ana->SwitchOnSMCombinations();
  
  ana->SwitchOffMultipleCutAnalysis();
  ana->SwitchOnFillAngleHisto();
  ana->SwitchOnFillOriginHisto(); // MC
  
  // Set Histograms name tag, bins and ranges
  if(!bothCalo) ana->AddToHistogramsName(Form("AnaPi0_"));
  else          ana->AddToHistogramsName(Form("AnaPi0_2Det_"));
  
  SetAnalysisCommonParameters(ana,histoString,calorimeter,year,col,simulation,printSettings,debug); // see method below
  if(printSettings) ana->Print("");
  
  return ana;
}


AliAnaParticleIsolation* ConfigureIsolationAnalysis(TString particle,      Int_t   leading,
                                                    Int_t   partInCone,    Int_t   thresType,
                                                    Float_t cone,          Float_t coneMin,
                                                    Float_t pth,           Bool_t  multi,
                                                    TString col,           Bool_t  simulation,
                                                    TString calorimeter,   Int_t   year,       Int_t tm,
                                                    Bool_t  printSettings, Int_t   debug, 
                                                    TString histoString                 )
{
  AliAnaParticleIsolation *ana = new AliAnaParticleIsolation();
  
  ana->SetDebug(debug);
  
  ana->SetMinPt(5);
  ana->SetCalorimeter(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetCalorimeter(calo);
  }
  
  ana->SwitchOnUEBandSubtractionHistoFill();
  ana->SwitchOffCellHistoFill() ;
  
  ana->SwitchOffLeadingOnly();
  ana->SwitchOffCheckNeutralClustersForLeading();
  if( leading > 0 )   ana->SwitchOnLeadingOnly();
  if( leading == 2 ||
     leading == 4)   ana->SwitchOnCheckNeutralClustersForLeading();
  
  // MC
  ana->SwitchOnPrimariesInConeSelection();
  ana->SwitchOnPrimariesPi0DecayStudy() ;
  
  ana->SwitchOffSSHistoFill();

  if(particle.Contains("Photon"))
  {
    ana->SwitchOnSSHistoFill();

    if(kAnaCutsString.Contains("Decay"))
    {
      ana->SwitchOffDecayTaggedHistoFill() ;
      ana->SetNDecayBits(5);
    }
  }
 
  if ( kAnaCutsString.Contains("PerSM") ) 
  {
    ana->SwitchOnFillHistogramsPerSM(); 
    ana->SwitchOnStudyNCellsCut(); 
  }
  
  if ( kAnaCutsString.Contains("PerTCard") ) 
    ana->SwitchOnFillHistogramsPerTCardIndex(); 
  
  if(kAnaCutsString.Contains("Bkg"))
  {
    printf("ConfigureIsolationAnalysis() *** Activate analysis on PtTrig and Bkg bins >> \n");
    ana->SwitchOnPtTrigBinHistoFill();
    ana->SetNPtTrigBins(6);
    //ana->SetPtTrigLimits(0,8); ana->SetPtTrigLimits(1,12); ana->SetPtTrigLimits(2,16); ana->SetPtTrigLimits(3,25);
    
    ana->SwitchOnBackgroundBinHistoFill();
    ana->SetNBackgroundBins(11);
    //ana->SetBackgroundLimits(0,0); ana->SetBackgroundLimits(1,0.2); ana->SetBackgroundLimits(2,3); ana->SetBackgroundLimits(3,0.4);
  }
  
  if(!tm)  ana->SwitchOnTMHistoFill();
  else     ana->SwitchOffTMHistoFill();
  
  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
  
  // Avoid borders of calorimeter
  ana->SwitchOffRealCaloAcceptance();
  ana->SwitchOnFiducialCut();
  if      ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60,  86, 174) ;
  else if ( calorimeter == "DCAL"  ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60, 266, 314) ; 
  else if ( calorimeter == "PHOS"  ) ana->GetFiducialCut()->SetSimplePHOSFiducialCut (0.10, 266, 314) ; 
  
  // Same Eta as EMCal, cut in phi if EMCAL was triggering
  if(particle=="Hadron"  || particle.Contains("CTS"))
  {
    //if(trigger.Contains("EMC"))
    //  ana->GetFiducialCut()->SetSimpleCTSFiducialCut  (0.6, 260, 360) ;
    //else
    ana->GetFiducialCut()->SetSimpleCTSFiducialCut  (0.6, 0, 360) ;
  }
  
  // Branch AOD settings
  
  ana->SetInputAODName(Form("%sTrigger_%s",particle.Data(),kAnaCaloTrackCorr.Data()));
  ana->SetAODObjArrayName(Form("IC%sTrigger_%s_R%1.1f_ThMin%1.1f",particle.Data(),kAnaCaloTrackCorr.Data(),cone,pth));
  
  //
  // Do settings for main isolation cut class
  //
  AliIsolationCut * ic =  ana->GetIsolationCut();
  ic->SetDebug(debug);
  ic->SetParticleTypeInCone(partInCone);
  ic->SetICMethod(thresType);
  ic->SetPtFraction(0.1);
  ic->SetPtThreshold(0.5); // default, change in next lines
  ic->SetSumPtThreshold(1.0); // default, change in next lines
  
  if(cone > 0 && pth > 0)
  {
    ic->SetConeSize(cone);
    ic->SetMinDistToTrigger(coneMin);    
    ic->SetPtThresholdMax(10000);
    
    if(thresType == AliIsolationCut::kPtThresIC)
    {
      printf("ConfigureIsolationAnalysis() *** Iso *** PtThresMin = %1.1f GeV/c *** R = %1.2f *** R min %1.2f\n",pth,cone,coneMin);
      ic->SetPtThreshold(pth);
    }
    
    if(thresType == AliIsolationCut::kSumPtIC)
    {
      printf("ConfigureIsolationAnalysis() *** Iso *** SumPtMin = %1.1f GeV/c *** R = %1.1f *** R min %1.2f\n",pth,cone,coneMin);
      ic->SetSumPtThreshold(pth);
    }
  }
  else
  {
    if(col=="pp")
    {
      ic->SetPtThreshold(0.5);
      ic->SetSumPtThreshold(1.0) ;
      ic->SetConeSize(0.4);
      ic->SetMinDistToTrigger(-1);    
    }
    if(col=="PbPb")
    {
      ic->SetPtThreshold(3.);
      ic->SetSumPtThreshold(3.0) ;
      ic->SetConeSize(0.3);
      ic->SetMinDistToTrigger(-1);    
    }
  }
  
  
  // Do or not do isolation with previously produced AODs.
  // No effect if use of SwitchOnSeveralIsolation()
  ana->SwitchOffReIsolation();
  
  // Multiple IC
  if(multi)
  {
    ic->SetConeSize(1.);    // Take all for first iteration
    ic->SetPtThreshold(100);// Take all for first iteration
    ic->SetMinDistToTrigger(-1);    

    ana->SwitchOnSeveralIsolation() ;
    //ana->SetAODObjArrayName(Form("MultiIC%sTM%d",particle.Data(),tm));
    ana->SetAODObjArrayName(Form("MultiIC%s",particle.Data()));
    
    ana->SetNCones(3);
    ana->SetNPtThresFrac(2);
    ana->SetConeSizes(0,0.3);       ana->SetConeSizes(1,0.4);       ana->SetConeSizes(2,0.5);
    ana->SetPtThresholds(0, 0.5);   ana->SetPtThresholds(1, 1);     ana->SetPtThresholds(2, 1.5);  ana->SetPtThresholds(3, 2);
    ana->SetPtFractions (0, 0.05) ; ana->SetPtFractions (1, 0.1);   ana->SetPtFractions (2, 0.2) ;  ana->SetPtFractions (3, 0.3) ;
    ana->SetSumPtThresholds(0, 0.5) ; ana->SetSumPtThresholds(1, 1) ; ana->SetSumPtThresholds(2, 1.5);  ana->SetSumPtThresholds(3, 2)  ;
    //ana->SetPtThresholds(0, 0.5);
    
    ana->SwitchOffTMHistoFill();
    ana->SwitchOffSSHistoFill();
  }
  else
    ana->SwitchOffSeveralIsolation() ;
  
  AliCaloPID* caloPID = ana->GetCaloPID();
  
  // Track matching 
  // tm = 1, fixed cuts
  caloPID->SetEMCALDEtaCut(0.025);
  caloPID->SetEMCALDPhiCut(0.030);
  
  // pT track dependent cuts
  if(tm > 1) caloPID->SwitchOnEMCTrackPtDepResMatching();
  
  //Set Histograms name tag, bins and ranges
  
//  if(!multi) ana->AddToHistogramsName(Form("AnaIsol%s_TM%d_"     ,particle.Data(),tm));
//  else       ana->AddToHistogramsName(Form("AnaMultiIsol%s_TM%d_",particle.Data(),tm));  
  if(!multi) ana->AddToHistogramsName(Form("AnaIsol%s_"     ,particle.Data()));
  else       ana->AddToHistogramsName(Form("AnaMultiIsol%s_",particle.Data()));
  
  SetAnalysisCommonParameters(ana,histoString, calorimeter,year,col,simulation,printSettings,debug); // see method below
  
  if(particle=="Hadron"  || particle.Contains("CTS"))
  {
    ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 200) ;
    ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-1.5, 1.5, 300) ;
  }
  
  if(printSettings) ic ->Print("");
  
  return ana;
}

AliAnaParticleHadronCorrelation* ConfigureHadronCorrelationAnalysis(TString particle,      Int_t   leading,
                                                                    Bool_t  bIsolated,     Float_t shshMax,
                                                                    Int_t   partInCone,    Int_t   thresType,
                                                                    Float_t cone,          Float_t coneMin,
                                                                    Float_t pth,           Bool_t  mixOn,
                                                                    TString col,           Bool_t  simulation,
                                                                    TString calorimeter,   Int_t   year,     Int_t tm,
                                                                    Bool_t  printSettings, Int_t   debug, 
                                                                    TString histoString                       )
{
  AliAnaParticleHadronCorrelation *ana = new AliAnaParticleHadronCorrelation();
  
  ana->SetTriggerPtRange(5,100);
  ana->SetAssociatedPtRange(0.2,100);
  ana->SetDeltaPhiCutRange  (TMath::DegToRad()*120.,TMath::DegToRad()*240.);
  
  // Underlying event
  ana->SetUeDeltaPhiCutRange(TMath::DegToRad()*60. ,TMath::DegToRad()*120.);
  ana->SwitchOnSeveralUECalculation();
  
  ana->SwitchOffAbsoluteLeading();  // Select trigger leading particle of all the selected tracks
  ana->SwitchOffNearSideLeading();  // Select trigger leading particle of all the particles at +-90 degrees, default
  ana->SwitchOffCheckNeutralClustersForLeading();
  
  if(leading >  0 && leading <  3 ) ana->SwitchOnAbsoluteLeading();
  if(leading >  2 )                 ana->SwitchOnNearSideLeading();
  if(leading == 2 || leading == 4 ) ana->SwitchOnCheckNeutralClustersForLeading();
  
  ana->SwitchOffFillPtImbalancePerPtABinHistograms();
  ana->SwitchOffCorrelationVzBin() ;
  ana->SwitchOffFillEtaGapHistograms();
    
  ana->SwitchOffPi0TriggerDecayCorr();
  
  if(particle.Contains("Photon"))
  {
    if(kAnaCutsString.Contains("Decay"))
    {
      printf("ConfigureHadronCorrelationAnalysis() *** Activate analysis on Tagged clusters as Decay product >> \n");
      ana->SwitchOnDecayTriggerDecayCorr();
      ana->SetNDecayBits(5);
      ana->SwitchOnInvariantMassHistograms();
      if ( kAnaCutsString.Contains("Bkg") && !bIsolated )
        ana->SwitchOnBackgroundBinsTaggedDecayPtInConeHistograms();
    }
    
    //printf("ConfigureHadronCorrelationAnalysis() *** SET M02 limits in correlation task *** \n");
    ana->SetM02Cut(0.10,shshMax);
  }
  
  if ( kAnaCutsString.Contains("PerSM") ) 
    ana->SwitchOnFillHistogramsPerSM(); 
 
  if ( kAnaCutsString.Contains("PerTCard") ) 
    ana->SwitchOnFillHistogramsPerTCardIndex(); 
  
  ana->SetMCGenType(0,7); // Change to higher to include Tagged decays
  
  ana->SwitchOffLeadHadronSelection(); // Open cuts, just fill histograms
  ana->SwitchOnFillLeadHadronHistograms();
  
  if ( kAnaCutsString.Contains("Bkg") && !bIsolated )
  {
    printf("ConfigureHadronCorrelationAnalysis() *** Activate analysis on PtTrig and Bkg bins,\n"
           " \t make sure isolation runs first even if not used to select trigger >> \n");
    ana->SwitchOnBackgroundBinsPtInConeHistograms();
  }
  
  ana->SetLeadHadronPhiCut(TMath::DegToRad()*130, TMath::DegToRad()*230.);
  ana->SetLeadHadronPtCut(0.5, 1000);
  
  // if triggering on PHOS and EMCAL is on
  ana->SwitchOffNeutralCorr(); // Do only correlation with TPC
                               //ana->SetPi0AODBranchName("Pi0EMCAL_TrigEMC7_Cl_TM1");
  
  ana->SwitchOffHMPIDCorrelation();
  
  ana->SwitchOffFillBradHistograms();
  
  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
  
  ana->SwitchOnFillDeltaEtaPhiPtTrigHistograms();
  
  ana->SetNAssocPtBins(16); // set last bin [20,30] GeV/c
  // See AliAnaParticleCorrelation::InitParameters();
  // Default bins{0.2,0.5,1,2,3,4,5,6,7,8,9,10,12,14,16,20,30,40,50,100} GeV/c
  // If you want to change it:
  //  ana->SetAssocPtBinLimit(0, 1) ;
  //  ana->SetAssocPtBinLimit(1, 2) ;
  //  ana->SetAssocPtBinLimit(2, 3) ;
  //  ana->SetAssocPtBinLimit(3, 4) ;
  //  ana->SetAssocPtBinLimit(4, 5) ;
  //  ana->SetAssocPtBinLimit(5, 8) ;
  //  ana->SetAssocPtBinLimit(6, 10) ;
  //  ana->SetAssocPtBinLimit(7, 100);
  
  ana->SetNTriggerPtBins(5); // set last bin [25,30] GeV/c
  // See AliAnaParticleCorrelation::InitParameters();
  // Default bins{10,12,16,20,25,30,40,50,75,100} GeV/c
  // If you want to change it:
  //  ana->SetTriggerPtBinLimit(0,  5) ;
  //  ana->SetTriggerPtBinLimit(1, 10) ;
  
  ana->SelectIsolated(bIsolated); // do correlation with isolated photons
  
  // Mixing with own pool
  if(mixOn)
  {
    ana->SwitchOnOwnMix();
    ana->SwitchOnFillNeutralInMixedEvent();
    
    if(bIsolated)
    {
      //Do settings for main isolation cut class
      AliIsolationCut * ic =  ana->GetIsolationCut();
      ic->SetDebug(debug);
      
      if(cone >0 && pth > 0)
      {
        //printf("ConfigureHadronCorrelationAnalysis() *** PtThres = %1.1f GeV/c *** R = %1.1f *** R min = %1.2f\n",pth,cone,coneMin);
        ic->SetPtThreshold(pth);
        ic->SetConeSize(cone);
        ic->SetMinDistToTrigger(coneMin);    
      }
      else
      {
        if(col=="pp")
        {
          ic->SetPtThreshold(0.5);
          ic->SetConeSize(0.4);
          ic->SetMinDistToTrigger(-1);    
        }
        if(col=="PbPb")
        {
          ic->SetPtThreshold(3.);
          //ic->SetPtThreshold(1.);
          ic->SetConeSize(0.3);
          ic->SetMinDistToTrigger(-1);    
        }
      }
      
      ic->SetPtFraction(0.1);
      ic->SetSumPtThreshold(1.0) ;
      ic->SetParticleTypeInCone(partInCone);
      ic->SetICMethod(thresType);
    }
  }
  else
    ana->SwitchOffOwnMix();
  
  ana->SetNZvertBin(20);
  
  if(col=="pp")
  {
    ana->SetNMaxEvMix(100);
    ana->SwitchOnTrackMultBins();
    ana->SetNTrackMultBin(10);
    ana->SetNRPBin(1);
  }
  else
  {
    ana->SetNMaxEvMix(10);
    ana->SwitchOffTrackMultBins(); // centrality bins
    ana->SetNCentrBin(12);
    ana->SetNRPBin(3);
    if(kAnaCaloTrackCorr.Contains("60_90"))
    {
      //printf("ConfigureHadronCorrelationAnalysis() *** Set mixing for peripheral ***\n");
      ana->SetNMaxEvMix(50);
      ana->SetNCentrBin(2);
    }
  }
    
  // Avoid borders of calorimeter, same as for isolation
  //
  ana->SwitchOnFiducialCut();
  if      ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60,  86, 174) ;
  else if ( calorimeter == "DCAL"  ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60, 266, 314) ; 
  else if ( calorimeter == "PHOS"  ) ana->GetFiducialCut()->SetSimplePHOSFiducialCut (0.10, 266, 314) ; 
  
  // Track matching, in case of mixing and specially in isolation case
  //
  AliCaloPID* caloPID = ana->GetCaloPID();
  // tm = 1, fixed cuts
  caloPID->SetEMCALDEtaCut(0.025);
  caloPID->SetEMCALDPhiCut(0.030);
  
  // pT track dependent cuts
  if(tm > 1) caloPID->SwitchOnEMCTrackPtDepResMatching();
  
  // Input / output delta AOD settings
  //
  ana->SetInputAODName(Form("%sTrigger_%s",particle.Data(),kAnaCaloTrackCorr.Data()));
  ana->SetAODObjArrayName(Form("%sHadronCorrIso%dTrigger_%s",particle.Data(),bIsolated,kAnaCaloTrackCorr.Data()));
  //ana->SetAODNamepTInConeHisto(Form("IC%s_%s_R%1.1f_ThMin%1.1f"           ,particle.Data(),kAnaCaloTrackCorr.Data(),cone,pth));
  
  //Set Histograms name tag, bins and ranges
  //
  //ana->AddToHistogramsName(Form("Ana%sHadronCorr_Iso%d_TM%d_",particle.Data(),bIsolated,tm));
  ana->AddToHistogramsName(Form("Ana%sHadronCorr_Iso%d_",particle.Data(),bIsolated));
  
  SetAnalysisCommonParameters(ana,histoString,calorimeter,year,col,simulation,printSettings,debug); // see method below
  
  if(particle=="Hadron"  || particle.Contains("CTS"))
  {
    ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 200) ;
    ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-1.5, 1.5, 300) ;
  }
  
  return ana;
}

AliAnaChargedParticles* ConfigureChargedAnalysis
( Bool_t simulation, Bool_t printSettings, Int_t   debug, TString histoString  )
{
  AliAnaChargedParticles *ana = new AliAnaChargedParticles();
  
  ana->SetDebug(debug);
  
  // selection cuts
  
  ana->SetMinPt(0.2);
  ana->SwitchOnFiducialCut();
  ana->GetFiducialCut()->SetSimpleCTSFiducialCut(0.8, 0, 360) ; //more restrictive cut in reader and after in isolation
  
  ana->SwitchOffFillVertexBC0Histograms();
  
  ana->SwitchOffFillPileUpHistograms();
  ana->SwitchOffFillTrackBCHistograms();
  
  // Branch AOD settings
  
  ana->SetOutputAODName(Form("HadronTrigger_%s",kAnaCaloTrackCorr.Data()));
  ana->SetOutputAODClassName("AliCaloTrackParticle"); // use if no correlation done
  
  //Set Histograms name tag, bins and ranges
  
  ana->AddToHistogramsName("AnaHadrons_");
  
  SetAnalysisCommonParameters(ana,histoString,"CTS",2012,"pp",simulation,printSettings,debug); // see method below
  
  return ana;
}

AliAnaClusterShapeCorrelStudies* ConfigureClusterShape
(Int_t tm, TString col , Bool_t  simulation, 
 TString calorimeter, Int_t year, 
 Bool_t  printSettings, Int_t   debug, TString histoString )
{
  AliAnaClusterShapeCorrelStudies *ana = new AliAnaClusterShapeCorrelStudies();
  
  ana->SetM02Min(-1);
  ana->SetNCellsPerClusterMin(-1);
  
  ana->SetCalorimeter(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetCalorimeter(calo);
  }
  
  if(simulation) ana->SetConstantTimeShift(615);
  
  ana->SwitchOffFiducialCut();
  
  ana->SwitchOnStudyClusterShape();
  
  ana->SwitchOnStudyEMCalModuleCells();
  
  ana->SwitchOnStudyColRowFromCellMax() ;

  ana->SwitchOffStudyClusterShapeParam();
  
  ana->SwitchOffStudyMatchedPID() ;
  
  ana->SwitchOffStudyWeight();
  
  ana->SetNCellBinLimits(3); // set to -1 for no analysis on predefined bins in nCell
  ana->SetDistToBadMin(2);
  
  ana->SwitchOffStudyTCardCorrelation() ;
  ana->SwitchOffStudyExotic();
  ana->SwitchOffStudyInvariantMass();
  ana->SwitchOffStudyCellTime() ;
  
  // PID cuts (Track-matching)
  ana->SwitchOnCaloPID(); // do PID selection, unless specified in GetCaloPID, selection not based on bayesian
  AliCaloPID* caloPID = ana->GetCaloPID();
  
  if(tm > 1)
  {
    // track pT dependent cut
    caloPID->SwitchOnEMCTrackPtDepResMatching();
    
    // Begining of histograms name
    ana->AddToHistogramsName("Shape_TMDep_");
  }
  else if ( tm )
  {
    // Fix
    caloPID->SwitchOffEMCTrackPtDepResMatching();
    caloPID->SetEMCALDEtaCut(0.025);
    caloPID->SetEMCALDPhiCut(0.030);
    
    // Begining of histograms name
    ana->AddToHistogramsName("Shape_TMFix_"); 
  }
  
  SetAnalysisCommonParameters(ana,histoString,calorimeter,year,col,simulation,printSettings,debug); // see method below
  
  AliHistogramRanges* histoRanges = ana->GetHistogramRanges();
  histoRanges->SetHistoPtRangeAndNBins(0, 50, 100) ; // Reduce a bit size
  
  return ana;
}

AliAnaCalorimeterQA* ConfigureQAAnalysis(TString col,           Bool_t  simulation,
                                         TString calorimeter,   Int_t   year,
                                         Bool_t  printSettings, Int_t   debug, 
                                         TString histoString  )
{
  AliAnaCalorimeterQA *ana = new AliAnaCalorimeterQA();
  
  ana->SetCalorimeter(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetCalorimeter(calo);
  }
  
  ana->SetTimeCut(-1e10,1e10); // Open time cut
  
  ana->SwitchOffCorrelation(); // make sure you switch in the reader PHOS and EMCAL cells and clusters if option is ON
  
  ana->SwitchOffRealCaloAcceptance();
  
  ana->SwitchOffFiducialCut();
  ana->SwitchOffFillAllTH3Histogram();
  ana->SwitchOffFillAllPositionHistogram();
  ana->SwitchOffFillAllPositionHistogram2();
  ana->SwitchOffStudyBadClusters() ;
  ana->SwitchOffFillAllCellTimeHisto() ;
  
  ana->SwitchOnFillAllTrackMatchingHistogram();
  
  ana->AddToHistogramsName("QA_"); // Begining of histograms name
  
  SetAnalysisCommonParameters(ana,histoString,calorimeter,year,col,simulation,printSettings,debug); // see method below
  
  return ana;
}

AliAnaGeneratorKine* ConfigureGenKineAnalysis(Int_t   thresType,     Float_t pth,    
                                              Float_t cone,          Float_t coneMin,
                                              TString col,           Bool_t  simulation,
                                              TString calorimeter,   Int_t   year,
                                              Bool_t  printSettings, Int_t   debug, 
                                              TString histoString        )
{
  AliAnaGeneratorKine *ana = new AliAnaGeneratorKine();
  
  // Trigger detector, acceptance and pT cut
  ana->SetTriggerDetector(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetTriggerDetector(calo);
  }
  
  ana->SetMinPt(2); // Trigger photon, pi0 minimum pT
  if      ( calorimeter == "EMCAL" ) ana->GetFiducialCutForTrigger()->SetSimpleEMCALFiducialCut(0.60,  86, 174) ;
  else if ( calorimeter == "DCAL"  ) ana->GetFiducialCutForTrigger()->SetSimpleEMCALFiducialCut(0.60, 266, 314) ; 
  else if ( calorimeter == "PHOS"  ) ana->GetFiducialCutForTrigger()->SetSimplePHOSFiducialCut (0.10, 266, 314) ; 
  
  // Particles associated to trigger or isolation cone acceptance and pT cut
  ana->SetCalorimeter(calorimeter);
  if(calorimeter == "DCAL") 
  {
    TString calo = "EMCAL";
    ana->SetCalorimeter(calo);
  }
  
  ana->SetMinChargedPt(0.2);
  ana->SetMinNeutralPt(0.3);
  
  if      ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.70,  80, 174) ;
  else if ( calorimeter == "DCAL"  ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.70, 260, 327) ; 
  else if ( calorimeter == "PHOS"  ) ana->GetFiducialCut()->SetSimplePHOSFiducialCut (0.12, 250, 320) ;   
  ana->GetFiducialCut()->SetSimpleCTSFiducialCut(0.9, 0, 360);
  
  // Isolation paramters
  AliIsolationCut * ic = ana->GetIsolationCut();
  ic->SetDebug(debug);
  ic->SetPtThreshold(pth);
  ic->SetConeSize(cone);
  ic->SetMinDistToTrigger(coneMin);    
  ic->SetSumPtThreshold(1.0) ;
  ic->SetICMethod(thresType);
  
  ana->AddToHistogramsName("AnaGenKine_");
  
  SetAnalysisCommonParameters(ana,histoString,calorimeter,year,col,simulation,printSettings,debug); // see method below
  
  return ana;
}

void ConfigureCaloTrackCorrAnalysis
(
 TList*   anaList       = 0x0, 
 TString  calorimeter   = "EMCAL", // "DCAL", "PHOS"
 Bool_t   simulation    = kFALSE,
 Int_t    year          = 2011,
 TString  col           = "pp",
 TString  analysisString= "Photon_MergedPi0_DecayPi0_Isolation_Correlation_QA_Charged",
 TString  histoString   = "",
 Float_t  shshMax       = 0.27,
 Float_t  isoCone       = 0.4,
 Float_t  isoConeMin    = -1,
 Float_t  isoPtTh       = 2,
 Int_t    isoMethod     = AliIsolationCut::kSumPtIC,
 Int_t    isoContent    = AliIsolationCut::kNeutralAndCharged,
 Int_t    leading       = 0,
 Int_t    tm            = 2,
 Bool_t   mixOn         = kTRUE,
 Bool_t   printSettings = kFALSE,
 Int_t    debug         = 0
)
{
  if ( !anaList ) 
  {
    printf("ConfigureCaloTrackCorrAnalysis() << No analysis list passed. Stop >>\n");
    return;
  }
  
  // Check the passed variables
  //
  printf("ConfigureCaloTrackCorrAnalysis() << Settings: Base string <%s>, Analysis string <%s>, "
         "\n calorimeter <%s>, simulation <%d>, year <%d>, col <%s>, "
         "\n shshMax <%2.2f>, R <%1.2f>, Rmin <%1.2f>, isoPtTh <%2.2f>, isoMethod <%d>,isoContent <%d>,"
         "\n leading <%d>, tm <%d>, mixOn <%d>, printSettings <%d>, debug <%d>\n",
         anaList->GetName(), analysisString.Data(), 
         calorimeter.Data(), simulation, year, col.Data(),
         shshMax,isoCone,isoConeMin,isoPtTh,isoMethod,isoContent,
         leading,tm,mixOn,printSettings,debug);
  
  kAnaCutsString    = analysisString;
  
  kAnaCaloTrackCorr = Form("%s",anaList->GetName());

  if ( analysisString.Contains("Isolation") )
    kAnaCaloTrackCorr+= Form("_Iso_Meth%d_Part%d_Pt%1.2f_R%1.2f",isoMethod,isoContent,isoPtTh,isoCone);
  if ( analysisString.Contains("Corr") && analysisString.Contains("Photon") && shshMax > 0)
    kAnaCaloTrackCorr+= Form("_CorrM02_%1.2f",shshMax);
  
  if ( isoConeMin > 0 ) kAnaCaloTrackCorr+=Form("_Rmin%1.2f",isoConeMin);
  if ( leading    > 0 ) kAnaCaloTrackCorr+=Form("_Lead%d"   ,leading);
  if ( tm         > 0 ) kAnaCaloTrackCorr+=Form("_TM%d"     ,tm);
 
  printf("ConfigureCaloTrackCorrAnalysis() <<<< TMP branch internal NAME: %s >>>>>\n",kAnaCaloTrackCorr.Data());
  
  // #### Configure analysis ####
  

  Int_t n = anaList->GetEntries();//Analysis number, order is important
  
  //
  // Photon analysis
  //
  if ( analysisString.Contains("Photon") )
  {
    anaList->AddAt(ConfigurePhotonAnalysis
                   (col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Photon cluster selection
    
    if ( analysisString.Contains("InvMass") )
      anaList->AddAt(ConfigureInvariantMassAnalysis
                     (col,simulation,calorimeter,/*bothCalo*/kFALSE,
                      year,/*mix*/kTRUE,printSettings,debug,histoString), n++);
    
    if ( analysisString.Contains("DecayPi0") )
    {
      anaList->AddAt(ConfigurePi0EbEAnalysis
                     ("Pi0"        , AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
                      col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Pi0 event by event selection, invariant mass and photon tagging from decay
      
      anaList->AddAt(ConfigurePi0EbEAnalysis
                     ("Eta"        , AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
                      col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Eta event by event selection, invariant mass and photon tagging from decay
      
      anaList->AddAt(ConfigurePi0EbEAnalysis
                     ("Pi0SideBand", AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
                      col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Pi0 out of peak event by event selection, and photon tagging from decay
      
      anaList->AddAt(ConfigurePi0EbEAnalysis
                     ("EtaSideBand", AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
                      col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Eta out of peak event by event selection, and photon tagging from decay
    }
    
    if ( analysisString.Contains("Isolation") )
    {
      anaList->AddAt(ConfigureIsolationAnalysis
                     ("Photon", leading, isoContent,isoMethod,isoCone,isoConeMin,isoPtTh, kFALSE,
                      col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Photon isolation
      
    }
    
    if( analysisString.Contains("Correlation") )
    {
      if ( !analysisString.Contains("MultiIso") )
      {
        anaList->AddAt(ConfigureHadronCorrelationAnalysis
                       ("Photon", leading, kFALSE, shshMax, isoContent,isoMethod,isoCone,isoConeMin,isoPtTh, mixOn,
                        col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Gamma-hadron correlation
      }
      
      if ( analysisString.Contains("Isolation")  )
      {
        anaList->AddAt(ConfigureHadronCorrelationAnalysis
                       ("Photon", leading, kTRUE,  shshMax, isoContent,isoMethod,isoCone,isoConeMin,isoPtTh, mixOn,
                        col,simulation,calorimeter,year,tm,printSettings,debug,histoString) , n++); // Isolated gamma hadron correlation
      }
    } // correlation
  }
  
  //
  // Merged pi0 analysis
  //
  if ( analysisString.Contains("MergedPi0") )
  {
    anaList->AddAt(ConfigurePi0EbEAnalysis
                   ("Pi0", AliAnaPi0EbE::kSSCalo,kTRUE,kTRUE,
                    col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++); // Pi0 event by event selection, cluster splitting
    if ( analysisString.Contains("Isolation") )
    {
      anaList->AddAt(ConfigureIsolationAnalysis
                     ("Pi0SS", leading, isoContent,isoMethod,isoCone,isoConeMin,isoPtTh, kFALSE,
                      col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++);          // Pi0 isolation, cluster splits
      
    }
    
    if( analysisString.Contains("Correlation") )
    {
      if ( !analysisString.Contains("MultiIso") )
      {
        anaList->AddAt(ConfigureHadronCorrelationAnalysis
                       ("Pi0SS", leading, kFALSE, shshMax, isoContent,isoMethod,isoCone,isoConeMin,isoPtTh, mixOn,
                        col,simulation,calorimeter,year,tm,printSettings,debug,histoString), n++);  // Pi0-hadron correlation
      }
      
      if ( analysisString.Contains("Isolation") )
      {
        anaList->AddAt(ConfigureHadronCorrelationAnalysis
                       ("Pi0SS", leading, kTRUE,  shshMax, isoContent,isoMethod,isoCone,isoConeMin,isoPtTh, mixOn,
                        col,simulation,calorimeter,year,tm,printSettings,debug,histoString) , n++); // Isolated pi0-hadron correlation
      }
    } // correlation
  }
  
  // Check the generated kinematics
  if(simulation && analysisString.Contains("Generator")) 
  {
    anaList->AddAt(ConfigureGenKineAnalysis
                   (isoMethod,isoPtTh,isoCone,isoConeMin,
                    col,simulation,calorimeter,year,printSettings,debug,histoString), n++);
  }
  
  // Charged analysis
  if ( analysisString.Contains("Charged") ) 
  { 
    anaList->AddAt(ConfigureChargedAnalysis(simulation,printSettings,debug,histoString), n++); // track selection checks
  }
  
  // Cluster Shape studies
  if ( analysisString.Contains("ClusterShape") ) 
  { 
    anaList->AddAt(ConfigureClusterShape(tm,col,simulation,calorimeter,year,printSettings,debug,histoString) , n++);
  }
  
  // Calo QA
  if ( analysisString.Contains("QA") ) 
  { 
    anaList->AddAt(ConfigureQAAnalysis(col,simulation,calorimeter,year,printSettings,debug,histoString) , n++);
  }
  
  printf("ConfigureCaloTrackCorrAnalysis() << End configuration for %s with total analysis %d>>\n",
         kAnaCaloTrackCorr.Data(), n);
}