///
/// Configure the class handling the calorimeter clusters specific methods
///
AliCalorimeterUtils* ConfigureCaloUtils(TString calorimeter, TString trigger,
Bool_t simulation , Bool_t calibrate,
Int_t year , Int_t debugLevel)
{
AliCalorimeterUtils *cu = new AliCalorimeterUtils;
cu->SetDebug(debugLevel);
// Remove clusters close to borders, at least max energy cell is 1 cell away
cu->SetNumberOfCellsFromEMCALBorder(1);
cu->SetNumberOfCellsFromPHOSBorder(2);
// Search of local maxima in cluster
cu->SetLocalMaximaCutE(0.1);
cu->SetLocalMaximaCutEDiff(0.03);
//cu->SwitchOffClusterPlot();
cu->SwitchOffRecalculateClusterTrackMatching();
cu->SwitchOnBadChannelsRemoval() ;
//EMCAL settings
if(!simulation)
cu->SwitchOnLoadOwnEMCALGeometryMatrices();
cu->SwitchOffRecalibration(); // Check the reader if it is taken into account during filtering
cu->SwitchOffRunDepCorrection();
cu->SwitchOffCorrectClusterLinearity();
Bool_t bExotic = kTRUE;
Bool_t bNonLin = kFALSE;
Bool_t bBadMap = kTRUE;
Bool_t bEnCalib = kFALSE;
Bool_t bTiCalib = kFALSE;
if(calibrate && !simulation)
{
cu->SwitchOnRecalibration(); // Check the reader if it is taken into account during filtering
cu->SwitchOffRunDepCorrection();
cu->SwitchOnRecalculateClusterPosition() ;
bEnCalib = kTRUE;
bTiCalib = kTRUE;
}
// AliEMCALRecoUtils * recou = cu->GetEMCALRecoUtils();
//
// gROOT->LoadMacro("$ALICE_PHYSICS/PWGPP/EMCAL/macros/ConfigureEMCALRecoUtils.C");
// ConfigureEMCALRecoUtils(recou,
// simulation,
// bExotic,
// bNonLin,
// bEnCalib,
// bBadMap,
// bTiCalib,
// debugLevel
// );
cu->ConfigureEMCALRecoUtils(simulation,
bExotic,
bNonLin,
bEnCalib,
bBadMap,
bTiCalib,
debugLevel);
//recou->SetExoticCellDiffTimeCut(50.);
if(calorimeter=="PHOS")
{
if(year < 2014) cu->SetNumberOfSuperModulesUsed(3);
else cu->SetNumberOfSuperModulesUsed(4);
}
else
{
Int_t nSM = 20;
Int_t lastEMC = 11;
if (year == 2010) { nSM = 4; lastEMC = 3; }// EMCAL first year
else if (year < 2014) { nSM = 10; lastEMC = 9; }// EMCAL active 2011-2013
cu->SetNumberOfSuperModulesUsed(nSM);
if (trigger.Contains("EMCAL"))
{
cu->SetFirstSuperModuleUsed( 0);
cu->SetLastSuperModuleUsed (lastEMC);
}
else if (trigger.Contains("DCAL"))
{
cu->SetFirstSuperModuleUsed(12);
cu->SetLastSuperModuleUsed (19);
}
else
{
cu->SetFirstSuperModuleUsed(0);
cu->SetLastSuperModuleUsed (cu->GetNumberOfSuperModulesUsed()-1);
}
printf("AddTaskPi0IMGammaCorrQA::CalorimeterUtils() - nSM %d, first %d, last %d\n",
cu->GetNumberOfSuperModulesUsed(),cu->GetFirstSuperModuleUsed(), cu->GetLastSuperModuleUsed());
}
// PHOS
cu->SwitchOffLoadOwnPHOSGeometryMatrices();
if(debugLevel > 0) cu->Print("");
return cu;
}
///
/// Main method calling all the configuration
/// It creates a CaloTrackCorrelations task for calorimeters performance studies,
/// configures it and adds it to the analysis manager.
///
/// The options that can be passed to the macro are:
/// \param suffix : A string with the type of trigger (default: MB, EMC)
/// \param simulation : A bool identifying the data as simulation
/// \param outputFile : A string to change the name of the histograms output file, default is AnalysisResults.root
/// \param year: The year the data was taken, used to configure some histograms
/// \param printSettings : A bool to enable the print of the settings per task
/// \param calibrate: if OADB was updated with calibration parameters not used in reconstruction, apply them here.
///
AliAnalysisTaskCaloTrackCorrelation *AddTaskCalorimeterQA(const char *suffix="default",
Bool_t simulation = kFALSE,
TString outputFile = "",
Int_t year = 2015,
Bool_t printSettings = kFALSE,
Bool_t calibrate = kTRUE)
{
// Get the pointer to the existing analysis manager via the static access method.
//==============================================================================
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
if (!mgr) {
::Error("AddTaskCalorimeterQA", "No analysis manager to connect to.");
return NULL;
}
// Check the analysis type using the event handlers connected to the analysis manager.
//==============================================================================
if (!mgr->GetInputEventHandler()) {
::Error("AddTaskCalorimeter", "This task requires an input event handler");
return NULL;
}
TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
Bool_t kUseKinematics = (mgr->GetMCtruthEventHandler())?kTRUE:kFALSE;
TString ssuffix = suffix;
if(kUseKinematics || simulation)
{
simulation = kTRUE;
printf("AddTaskCalorimeterQA - CAREFUL : Triggered events not checked in simulation!! \n");
if(!ssuffix.Contains("default")) return NULL;
}
// Configure analysis
//===========================================================================
//Reader
//For this particular analysis few things done by the reader.
//Nothing else needs to be set.
AliCaloTrackReader * reader = 0x0;
if (inputDataType.Contains("AOD")) reader = new AliCaloTrackAODReader();
else if(inputDataType.Contains("ESD")) reader = new AliCaloTrackESDReader();
//reader->SetDebug(10);//10 for lots of messages
reader->SwitchOnEMCALCells();
reader->SwitchOnEMCAL();
reader->SwitchOnPHOSCells(); // For correlation plots
reader->SwitchOnPHOS(); // For correlation plots
reader->SetEMCALPtMin(0.);
reader->SwitchOnCTS();
reader->SetCTSPtMin (0.);
reader->SetZvertexCut(10.);
reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches.
reader->SwitchOffWriteDeltaAOD() ;
if(!ssuffix.Contains("default"))
{
reader->SwitchOnTriggerPatchMatching();
reader->SwitchOffBadTriggerEventsRemoval();
reader->SetTriggerPatchTimeWindow(8,9);
//reader->SetEventTriggerL0Threshold(2.);
}
if(!simulation) reader->AnalyzeOnlyPhysicsEvents(); // in case physics selection was not on
if(calibrate && !simulation) reader->SwitchOnClusterRecalculation();
else reader->SwitchOffClusterRecalculation();
if(printSettings) reader->Print("");
// *** Calorimeters Utils ***
AliCalorimeterUtils *cu = new AliCalorimeterUtils;
// Remove clusters close to borders, at least max energy cell is 1 cell away
cu->SetNumberOfCellsFromEMCALBorder(1);
if (year == 2010) cu->SetNumberOfSuperModulesUsed(4); //EMCAL first year
else if (year < 2014) cu->SetNumberOfSuperModulesUsed(10);
else cu->SetNumberOfSuperModulesUsed(20);
cu->SwitchOffRecalibration(); // Check the reader if it is taken into account during filtering
cu->SwitchOffRunDepCorrection();
cu->SwitchOnCorrectClusterLinearity();
Bool_t bExotic = kTRUE;
Bool_t bNonLin = kTRUE;
Bool_t bBadMap = kTRUE;
Bool_t bEnCalib = kFALSE;
Bool_t bTiCalib = kFALSE;
if(calibrate && !simulation)
{
cu->SwitchOnRecalibration(); // Check the reader if it is taken into account during filtering
cu->SwitchOffRunDepCorrection();
cu->SwitchOnRecalculateClusterPosition() ;
bEnCalib = kTRUE;
bTiCalib = kTRUE;
}
AliEMCALRecoUtils* recou = cu->GetEMCALRecoUtils();
gROOT->LoadMacro("$ALICE_PHYSICS/PWGPP/EMCAL/macros/ConfigureEMCALRecoUtils.C");
TString recouStr = Form("(reinterpret_cast<AliEMCALRecoUtils*>(%p))", recou);
gInterpreter->ProcessLine(Form("ConfigureEMCALRecoUtils(%s, %d, %d, %d, %d, %d, %d)", recouStr.Data(),
(Int_t)simulation,
(Int_t)bExotic,
(Int_t)bNonLin,
(Int_t)bEnCalib,
(Int_t)bBadMap,
(Int_t)bTiCalib));
if(bBadMap)
cu->SwitchOnBadChannelsRemoval();
cu->SetDebug(-1);
if(printSettings) cu->Print("");
// ##### Analysis algorithm settings ####
AliAnaCalorimeterQA *emcalQA = new AliAnaCalorimeterQA();
//emcalQA->SetDebug(10); //10 for lots of messages
TString calorimeter = "EMCAL";
emcalQA->SetCalorimeter(calorimeter);
if(simulation)
{
// Access MC stack and fill more histograms
emcalQA->SwitchOnDataMC() ;
emcalQA->SwitchOffStudyBadClusters();
emcalQA->SwitchOffFillAllCellTimeHisto();
}
else
{
emcalQA->SwitchOffDataMC() ;
emcalQA->SwitchOffStudyBadClusters();
emcalQA->SwitchOnFillAllCellTimeHisto();
}
emcalQA->AddToHistogramsName("EMCAL_"); //Begining of histograms name
emcalQA->SwitchOffFiducialCut();
emcalQA->SwitchOnCorrelation();
emcalQA->SwitchOffFillAllTH3Histogram();
emcalQA->SwitchOffFillAllPositionHistogram();
emcalQA->SwitchOffFillAllPositionHistogram2();
//Set Histrograms bins and ranges
emcalQA->GetHistogramRanges()->SetHistoPtRangeAndNBins(0, 50, 100) ;
emcalQA->GetHistogramRanges()->SetHistoFinePtRangeAndNBins(0, 10, 200) ; // bining for fhAmpId
emcalQA->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-0.70, 0.70, 140) ;
if ( year==2010 )
{
emcalQA->GetHistogramRanges()->SetHistoPhiRangeAndNBins(80*TMath::DegToRad(), 120*TMath::DegToRad(), 48) ;
emcalQA->GetHistogramRanges()->SetHistoXRangeAndNBins(-230,90,120);
emcalQA->GetHistogramRanges()->SetHistoYRangeAndNBins(370,450,40);
}
else if ( year < 2014 )
{
emcalQA->GetHistogramRanges()->SetHistoPhiRangeAndNBins(80*TMath::DegToRad(), 180*TMath::DegToRad(), 120) ;
emcalQA->GetHistogramRanges()->SetHistoXRangeAndNBins(-460,90,200);
emcalQA->GetHistogramRanges()->SetHistoYRangeAndNBins(100,450,100);
}
else // Run2
{
emcalQA->GetHistogramRanges()->SetHistoPhiRangeAndNBins(80*TMath::DegToRad(), 327*TMath::DegToRad(), 250) ;
emcalQA->GetHistogramRanges()->SetHistoXRangeAndNBins(-460,460,230);
emcalQA->GetHistogramRanges()->SetHistoYRangeAndNBins(-450,450,225);
}
emcalQA->GetHistogramRanges()->SetHistoMassRangeAndNBins(0., 0.65, 325) ;
emcalQA->GetHistogramRanges()->SetHistoAsymmetryRangeAndNBins(0., 1. , 10 );
emcalQA->GetHistogramRanges()->SetHistoPOverERangeAndNBins(0,2.,50);
emcalQA->GetHistogramRanges()->SetHistodEdxRangeAndNBins(0.,200.,100);
emcalQA->GetHistogramRanges()->SetHistodRRangeAndNBins(0.,0.10,50);
//emcalQA->GetHistogramRanges()->SetHistoTimeRangeAndNBins( 400,900,250);
//emcalQA->GetHistogramRanges()->SetHistoTimeRangeAndNBins(-275,275,250);
emcalQA->GetHistogramRanges()->SetHistoTimeRangeAndNBins(-275,975,250);
emcalQA->GetHistogramRanges()->SetHistoRatioRangeAndNBins(0.,2.,100);
emcalQA->GetHistogramRanges()->SetHistoVertexDistRangeAndNBins(0.,500.,500);
emcalQA->GetHistogramRanges()->SetHistoNClusterCellRangeAndNBins(0,50,50);
emcalQA->GetHistogramRanges()->SetHistoZRangeAndNBins(-350,350,175);
emcalQA->GetHistogramRanges()->SetHistoRRangeAndNBins(430,460,30);
emcalQA->GetHistogramRanges()->SetHistoV0SignalRangeAndNBins(0,5000,100);
emcalQA->GetHistogramRanges()->SetHistoV0MultiplicityRangeAndNBins(0,5000,100);
emcalQA->GetHistogramRanges()->SetHistoTrackMultiplicityRangeAndNBins(0,2500,100);
emcalQA->GetHistogramRanges()->SetHistoShowerShapeRangeAndNBins(0, 3, 120);
emcalQA->GetHistogramRanges()->SetHistoDiffTimeRangeAndNBins(-300, 300, 120);
emcalQA->GetHistogramRanges()->SetHistoTrackResidualEtaRangeAndNBins(-0.075,0.075,50);
emcalQA->GetHistogramRanges()->SetHistoTrackResidualPhiRangeAndNBins(-0.075,0.075,50);
if(printSettings) emcalQA->Print("");
// #### Configure Maker ####
AliAnaCaloTrackCorrMaker * maker = new AliAnaCaloTrackCorrMaker();
if(ssuffix.Contains("default")) maker->SwitchOffDataControlHistograms();
else maker->SwitchOnDataControlHistograms();
maker->SetReader(reader);//pointer to reader
maker->SetCaloUtils(cu); //pointer to calorimeter utils
maker->AddAnalysis(emcalQA,0);
maker->SetAnaDebug(-1) ; // 0 to at least print the event number
maker->SwitchOnHistogramsMaker() ;
maker->SwitchOffAODsMaker() ;
if(simulation)
{
// Calculate the cross section weights, apply them to all histograms
// and fill xsec and trial histo. Sumw2 must be activated.
//maker->GetReader()->GetWeightUtils()->SwitchOnMCCrossSectionCalculation();
//maker->SwitchOnSumw2Histograms();
// For recent productions where the cross sections and trials are not stored in separate file
//maker->GetReader()->GetWeightUtils()->SwitchOnMCCrossSectionFromEventHeader() ;
// Just fill cross section and trials histograms.
maker->GetReader()->GetWeightUtils()->SwitchOnMCCrossSectionHistoFill();
// Add control histogram with pT hard to control aplication of weights
maker->SwitchOnPtHardHistogram();
}
if(printSettings) maker->Print("");
printf("======================== \n");
printf("AddTaskCalorimeterQA() - End Configuration \n");
printf("======================== \n");
// Create task
//===========================================================================
AliAnalysisTaskCaloTrackCorrelation * task = new AliAnalysisTaskCaloTrackCorrelation (Form("CalorimeterPerformance_%s",suffix));
//task->SetConfigFileName(""); //Don't configure the analysis via configuration file.
//task->SetDebugLevel(-1);
task->SetAnalysisMaker(maker);
task->SetBranches("ESD:AliESDRun.,AliESDHeader"); //just a trick to get Constantin's analysis to work
mgr->AddTask(task);
//Create containers
// AliAnalysisDataContainer *cout_pc = mgr->CreateContainer("Calo.Performance",TList::Class(),
// AliAnalysisManager::kOutputContainer, "Calo.Performance.root");
TString cname;
if(outputFile.Length()==0)outputFile = AliAnalysisManager::GetCommonFileName();
cname = Form("CaloQA_%s", suffix);
AliAnalysisDataContainer *cout_pc = mgr->CreateContainer(cname, TList::Class(),
AliAnalysisManager::kOutputContainer,
Form("%s:%s",outputFile.Data(),cname.Data()));
// cname = Form("CaloQACuts_%s", suffix);
// AliAnalysisDataContainer *cout_cuts = mgr->CreateContainer(cname, TList::Class(),
// AliAnalysisManager::kParamContainer,
// Form("%s:%s",outputFile.Data(),cname.Data()));
//Form("%s:PartCorrCuts",outputfile.Data()));
// Create ONLY the output containers for the data produced by the task.
// Get and connect other common input/output containers via the manager as below
//==============================================================================
mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer());
mgr->ConnectOutput (task, 1, cout_pc);
// mgr->ConnectOutput (task, 2, cout_cuts);
return task;
}
AliAnalysisTaskParticleCorrelation *AddTaskPartCorr(TString inputDataType, TString calorimeter, Bool_t kPrintSettings = kFALSE,Bool_t kSimulation = kFALSE, Bool_t outputAOD=kFALSE, Bool_t oldAOD=kFALSE, Int_t debugLevel = 0)
{
// Creates a PartCorr task, configures it and adds it to the analysis manager.
// Get the pointer to the existing analysis manager via the static access method.
//==============================================================================
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
if (!mgr) {
::Error("AddTaskPartCorr", "No analysis manager to connect to.");
return NULL;
}
Bool_t kUseKinematics = kFALSE;
if(kSimulation) {
kUseKinematics = (mgr->GetMCtruthEventHandler())?kTRUE:kFALSE;
if (!kUseKinematics && inputDataType == "AOD") kUseKinematics = kTRUE; //AOD primary should be available ...
}
cout<<"********* ACCESS KINE? "<<kUseKinematics<<endl;
// Configure analysis
//===========================================================================
// *** Reader ***
AliCaloTrackReader * reader =0x0 ;
if(inputDataType.Contains("AOD")) reader = new AliCaloTrackAODReader();
else if(inputDataType=="ESD") reader = new AliCaloTrackESDReader();
else if(inputDataType=="MC" && inputDataType == "ESD") reader = new AliCaloTrackMCReader();
reader->SetDebug(-1);//10 for lots of messages
reader->SwitchOnCTS();
reader->SwitchOffEMCAL();
reader->SwitchOffPHOS();
reader->SwitchOffEMCALCells();
reader->SwitchOffPHOSCells();
reader->SetDataType(2);
//Min particle pT
reader->SetEMCALPtMin(0.3);
reader->SetPHOSPtMin(0.3);
reader->SetCTSPtMin(0.1);
reader->SwitchOnWriteDeltaAOD() ;
// if(oldAOD) reader->SwitchOnOldAODs();
reader->SwitchOnWriteDeltaAOD();
if(kPrintSettings) reader->Print("");
// *** Calorimeters Utils ***
AliCalorimeterUtils *cu = new AliCalorimeterUtils;
//Recalibration
//cu->SwitchOnRecalibration();
//TFile * f = new TFile("RecalibrationFactors.root","read");
//cu->SetEMCALChannelRecalibrationFactors(0,(TH2F*)f->Get("EMCALRecalFactors_SM0"));
//cu->SetEMCALChannelRecalibrationFactors(1,(TH2F*)f->Get("EMCALRecalFactors_SM1"));
//cu->SetEMCALChannelRecalibrationFactors(2,(TH2F*)f->Get("EMCALRecalFactors_SM2"));
//cu->SetEMCALChannelRecalibrationFactors(3,(TH2F*)f->Get("EMCALRecalFactors_SM3"));
//f->Close();
cu->SetDebug(debugLevel);
if(kPrintSettings) cu->Print("");
// ##### Analysis algorithm settings ####
//==================================
// ### Isolation analysis ###
//=================================
//Photon
AliAnaParticleIsolation *anaisol = new AliAnaParticleIsolation();
anaisol->SetDebug(debugLevel);
anaisol->SetMinPt(0);
anaisol->SetInputAODName("ConversionGamma");
anaisol->SetAODObjArrayName("ConversionGamma");
anaisol->SetCalorimeter("CTS");
if(kUseKinematics) anaisol->SwitchOnDataMC() ;//Access MC stack and fill more histograms
else anaisol->SwitchOffDataMC() ;
anaisol->SwitchOffInvariantMass();
AliIsolationCut * ic = anaisol->GetIsolationCut();
ic->SetConeSize(0.4);
ic->SetPtThreshold(0.7);
ic->SetPtFraction(0.1);
ic->SetSumPtThreshold(1.0) ;
ic->SetParticleTypeInCone(AliIsolationCut::kOnlyCharged);
ic->SetICMethod(AliIsolationCut::kSumPtFracIC);
if(kPrintSettings) ic->Print("");
anaisol->SetHistoPtRangeAndNBins(0, 50, 500) ;
anaisol->AddToHistogramsName("AnaIsolPhoton_");
anaisol->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
anaisol->SwitchOnReIsolation();
anaisol->SwitchOffSeveralIsolation();
// anaisol->SetPtThresholds(0, 1.0);
// anaisol->SetPtThresholds(1, 2.0);
// anaisol->SetPtThresholds(2, 3.0);
// anaisol->SetPtThresholds(3, 4.0);
// anaisol->SetPtThresholds(4, 5.0);
// anaisol->SetPtFractions(0, 0.4);
// anaisol->SetPtFractions(1, 0.3);
// anaisol->SetPtFractions(2, 0.1);
// anaisol->SetPtFractions(3, 0.08);
// anaisol->SetPtFractions(4, 0.05);
// anaisol->SetNCones(5);
// anaisol->SetConeSizes(0, 0.50);
// anaisol->SetConeSizes(1, 0.60);
// anaisol->SetConeSizes(2, 0.70);
// anaisol->SetConeSizes(3, 0.80);
// anaisol->SetConeSizes(4, 0.90);
if(kPrintSettings) anaisol->Print("");
// //===========================
// //Correlation analysis
// //===========================
// // ### Correlation with Jet Finder AOD output
// AliAnaParticleJetFinderCorrelation *anacorrjet = new AliAnaParticleJetFinderCorrelation();
// anacorrjet->SetInputAODName("ConversionGamma");//,calorimeter.Data()));
// anacorrjet->SwitchOffFiducialCut();
// anacorrjet->SetDebug(debugLevel);
// anacorrjet->SetConeSize(1);
// anacorrjet->SelectIsolated(kTRUE); // do correlation with isolated photons
// anacorrjet->SetPtThresholdInCone(0.2);
// anacorrjet->SetDeltaPhiCutRange(0.5,5.5);//Mostly Open Cuts
// anacorrjet->SetRatioCutRange(0.01,3); //Mostly Open Cuts
// anacorrjet->UseJetRefTracks(kFALSE); //Not working now
// //Set Histograms bins and ranges
// anacorrjet->SetHistoPtRangeAndNBins(0, 50, 200) ;
// // ana->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
// // ana->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
// if(kPrintSettings) anacorrjet->Print("");
// #### Configure Maker ####
AliAnaPartCorrMaker * maker = new AliAnaPartCorrMaker();
maker->SetReader(reader);//pointer to reader
maker->SetCaloUtils(cu); //pointer to calorimeter utils
Int_t n = 0;//Analysis number, order is important
maker->AddAnalysis(anaisol,n++);
//maker->AddAnalysis(anacorrjet, n++);
maker->SetAnaDebug(-1);
maker->SwitchOnHistogramsMaker() ;
if(kPrintSettings) maker->Print("");
printf("======================== \n");
printf(" End Configuration of PartCorr analysis with detector %s \n",calorimeter.Data());
printf("======================== \n");
// Create task
//===========================================================================
AliAnalysisTaskParticleCorrelation * task = new AliAnalysisTaskParticleCorrelation (Form("PartCorr%s",calorimeter.Data()));
task->SetConfigFileName(""); //Don't configure the analysis via configuration file.
task->SetDebugLevel(debugLevel);
task->SetAnalysisMaker(maker);
if(inputDataType=="ESD" && !kSimulation) task->SelectCollisionCandidates(); //AliPhysicsSelection has to be attached before.
mgr->AddTask(task);
//Create containers
char name[128];
sprintf(name,"PartCorr_%s",calorimeter.Data());
cout<<"Name of task "<<name<<endl;
TString outputfile = AliAnalysisManager::GetCommonFileName();
AliAnalysisDataContainer *cout_pc = mgr->CreateContainer(calorimeter.Data(), TList::Class(),
AliAnalysisManager::kOutputContainer,
Form("%s:PartCorr",outputfile.Data()));
AliAnalysisDataContainer *cout_cuts = mgr->CreateContainer(Form("%sCuts",calorimeter.Data()), TList::Class(),
AliAnalysisManager::kParamContainer,
Form("%s:PartCorrCuts",outputfile.Data()));
// Create ONLY the output containers for the data produced by the task.
// Get and connect other common input/output containers via the manager as below
//==============================================================================
mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer());
// AOD output slot will be used in a different way in future
if(!inputDataType.Contains("delta") && outputAOD) mgr->ConnectOutput (task, 0, mgr->GetCommonOutputContainer());
mgr->ConnectOutput (task, 1, cout_pc);
mgr->ConnectOutput (task, 2, cout_cuts);
return task;
}
