Bool_t isEMCalTrain = kFALSE
)
{
+ cout << " ############ MACRO EXECUTION STARTED ############\n";
// #### Detect the demanded trigger with its readable name
TString triggerName(Form("Trigger_%i", trigger));
if (trigger == AliVEvent::kAnyINT)
gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
AliEmcalJetTask* jetFinderTask = AddTaskEmcalJet(usedTracks,"",1,jetRadius,1,minJetTrackPt,0.300); // anti-kt
AliEmcalJetTask* jetFinderTaskKT = AddTaskEmcalJet(usedTracks,"",0,ktJetRadius,1,minJetTrackPt,0.300); // kt
+ cout << " Jet finder tasks added: " << jetFinderTask << " + " << jetFinderTaskKT << endl;
// #### Define external rho task
- AliEmcalJetTask* jetFinderRho = AddTaskEmcalJet(usedTracks,"",1,0.4,1,0.150,0.300); // anti-kt
- AliEmcalJetTask* jetFinderRhoKT = AddTaskEmcalJet(usedTracks,"",0,0.4,1,0.150,0.300); // kt
+ AliEmcalJetTask* jetFinderRho = AddTaskEmcalJet(usedTracks,"",1,0.4,1,minJetTrackPt,0.300); // anti-kt
+ AliEmcalJetTask* jetFinderRhoKT = AddTaskEmcalJet(usedTracks,"",0,0.4,1,minJetTrackPt,0.300); // kt
+ cout << " Jet finder tasks (used for bgrd) added: " << jetFinderRho << " + " << jetFinderRhoKT << endl;
gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskRhoSparse.C");
AliAnalysisTaskRhoSparse* rhotask = AddTaskRhoSparse(jetFinderRhoKT->GetName(), NULL, usedTracks, "", bgrdName.Data(), 0.4,"TPC", 0., 5., 0, 0,2,kFALSE,bgrdName.Data(),kTRUE);
+ cout << " Background task added: " << rhotask << endl;
}
// #### Define analysis task
else
task = new AliAnalysisTaskChargedJetsPA(Form("AnalysisPA_%s", "", triggerName.Data()), usedTracks, "","", doJetProfileAnalysis, doTrackcutAnalysis);
+ cout << " Main task created: " << task << endl;
+
// #### Task preferences
task->SetIsPA(isPA);
task->SetCentralityType(centralityType);
task->SetNumberOfCentralityBins(numberOfCentralityBins);
task->SetDoJetAnalysis(doJetAnalysis);
+ cout << " Settings set." << endl;
// #### Add analysis task
manager->AddTask(task);
+ cout << " Task added to manager" << endl;
manager->ConnectInput(task, 0, manager->GetCommonInputContainer());
+ cout << " Input connected, common input container: " << manager->GetCommonInputContainer() << endl;
manager->ConnectOutput(task, 1, contHistos);
+ cout << " Output connected, contHistos: " << contHistos << endl;
+
if(doJetProfileAnalysis)
{
manager->ConnectOutput(task, 2, contJetProfile);
if(isEMCalTrain)
RequestMemory(task,200*1024);
+ cout << " ############ MACRO EXECUTION SUCCESSFUL, will return " << task << " ############\n";
return task;
}