//______________________________________________________________________________
void AddAnalysisTasks(Int_t merge){
// Add all analysis task wagons to the train
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
//
// Tender and supplies. Needs to be called for every event.
//
if ( VAR_OCDB_SNAPSHOT )
{
AliCDBManager::Instance()->SetDefaultStorage("local:///donotexist");
// explicitely wrong default storage
// so we take _EVERYTHING_ from the snapshot.
AliCDBManager::Instance()->SetRun(0);
AliCDBManager::Instance()->SetSnapshotMode("OCDB_rec.root");
}
if ( VAR_MAKE_COMPACT_ESD )
{
gROOT->LoadMacro("$ALICE_PHYSICS/PWG/muondep/AddTaskCompactTreeMaker.C");
AddTaskCompactTreeMaker();
}
if (iMUONCDBConnect) {
if (iMUONCDBConnect > 1) gROOT->LoadMacro("AddTaskMuonCDBConnect.C");
else gROOT->LoadMacro("$ALICE_PHYSICS/PWG/muondep/AddTaskMuonCDBConnect.C");
AliAnalysisTaskMuonCDBConnect *cdbConnect = AddTaskMuonCDBConnect();
if (!defaultStorage.IsNull()) cdbConnect->SetDefaultStorage(defaultStorage.Data());
if (!alignStorage.IsNull() || alignVersion >= 0 || alignSubVersion >= 0)
cdbConnect->SetAlignStorage(alignStorage.Data(), alignVersion, alignSubVersion);
if (!recoParamStorage.IsNull()) cdbConnect->SetRecoParamStorage(recoParamStorage.Data());
cdbConnect->LoadMagField();
cdbConnect->LoadGeometry();
cdbConnect->LoadMapping();
}
UInt_t offlineTriggerMask = 0;
if (usePhysicsSelection) {
// Physics selection task
gROOT->LoadMacro("$ALICE_PHYSICS/OADB/macros/AddTaskPhysicsSelection.C");
mgr->RegisterExtraFile("event_stat.root");
AliPhysicsSelectionTask *physSelTask = AddTaskPhysicsSelection(kTRUE);
offlineTriggerMask = AliVEvent::kAny;
}
// Centrality (only Pb-Pb)
if (iCollision && useCentrality) {
gROOT->LoadMacro("$ALICE_PHYSICS/OADB/COMMON/MULTIPLICITY/macros/AddTaskMultSelection.C");
AliMultSelectionTask *mult = AddTaskMultSelection(kFALSE);
}
// track selection
AliMuonTrackCuts *trackCuts = 0x0;
if (iMUONEfficiency) {
trackCuts = new AliMuonTrackCuts("stdCuts", "stdCuts");
trackCuts->SetAllowDefaultParams();
trackCuts->SetFilterMask(AliMuonTrackCuts::kMuMatchLpt | AliMuonTrackCuts::kMuEta | AliMuonTrackCuts::kMuThetaAbs | AliMuonTrackCuts::kMuPdca);
trackCuts->SetIsMC(kTRUE);
}
if (iMUONRefit) {
if (iMUONRefit > 1) gROOT->LoadMacro("AddTaskMuonRefit.C");
else gROOT->LoadMacro("$ALICE_PHYSICS/PWG/muondep/AddTaskMuonRefit.C");
AliAnalysisTaskMuonRefit* refit = AddTaskMuonRefit(-1., -1., kTRUE, -1., -1.);
if (!defaultStorage.IsNull()) refit->SetDefaultStorage(defaultStorage.Data());
if (!alignStorage.IsNull()) refit->SetAlignStorage(alignStorage.Data());
refit->RemoveMonoCathodClusters(kTRUE, kFALSE);
}
if (iMUONRefitVtx) {
if (iMUONRefitVtx > 1) gROOT->LoadMacro("AddTaskMuonRefitVtx.C");
else gROOT->LoadMacro("$ALICE_PHYSICS/PWG/muondep/AddTaskMuonRefitVtx.C");
AliAnalysisTaskMuonRefitVtx* refitVtx = AddTaskMuonRefitVtx(kFALSE, kFALSE, kTRUE);
if (!defaultStorage.IsNull()) refitVtx->SetDefaultStorage(defaultStorage.Data());
}
if(iESDMCLabelAddition) {
if(iESDMCLabelAddition > 1) gROOT->LoadMacro("AddTaskESDMCLabelAddition.C");
else gROOT->LoadMacro("$ALICE_PHYSICS/PWG/muondep/AddTaskESDMCLabelAddition.C");
AliAnalysisTaskESDMCLabelAddition *esdmclabel = AddTaskESDMCLabelAddition();
if (!defaultStorage.IsNull()) esdmclabel->SetDefaultStorage(defaultStorage.Data());
if (!alignStorage.IsNull()) esdmclabel->SetAlignStorage(alignStorage.Data());
if (!recoParamStorage.IsNull()) esdmclabel->SetRecoParamStorage(recoParamStorage.Data());
esdmclabel->DecayAsFake(kTRUE);
}
if (iMUONQA) {
if (iMUONQA > 1) gROOT->LoadMacro("AddTaskMuonQA.C");
else gROOT->LoadMacro("$ALICE_PHYSICS/PWGPP/PilotTrain/AddTaskMuonQA.C");
AliAnalysisTaskMuonQA* muonQA = AddTaskMuonQA(kFALSE);
if (usePhysicsSelection) muonQA->SelectCollisionCandidates(offlineTriggerMask);
muonQA->SetTrackCuts(trackCuts);
}
if (useMC && useTR && iMUONPerformance) {
if (iMUONPerformance > 1) gROOT->LoadMacro("AddTaskMuonPerformance.C");
else gROOT->LoadMacro("$ALICE_PHYSICS/PWGPP/MUON/dep/AddTaskMuonPerformance.C");
AliAnalysisTaskMuonPerformance* muonPerformance = AddTaskMuonPerformance();
if (usePhysicsSelection) muonPerformance->SelectCollisionCandidates(offlineTriggerMask);
if (!defaultStorage.IsNull()) muonPerformance->SetDefaultStorage(defaultStorage.Data());
if (!alignStorage.IsNull()) muonPerformance->SetAlignStorage(alignStorage.Data());
if (!recoParamStorage.IsNull()) muonPerformance->SetRecoParamStorage(recoParamStorage.Data());
muonPerformance->UseMCKinematics(kTRUE);
muonPerformance->SetMCTrigLevelFromMatchTrk(kTRUE);
}
if (iMUONEfficiency) {
if (iMUONEfficiency > 1) gROOT->LoadMacro("AddTaskMUONTrackingEfficiency.C");
else gROOT->LoadMacro("$ALICE_PHYSICS/PWGPP/MUON/dep/AddTaskMUONTrackingEfficiency.C");
AliAnalysisTaskMuonTrackingEff* muonEfficiency = AddTaskMUONTrackingEfficiency(kFALSE,kTRUE,"");
if (usePhysicsSelection) muonEfficiency->SelectCollisionCandidates(offlineTriggerMask);
if (!defaultStorage.IsNull()) muonEfficiency->SetDefaultStorage(defaultStorage.Data());
if (!alignStorage.IsNull()) muonEfficiency->SetAlignStorage(alignStorage.Data());
if (!recoParamStorage.IsNull()) muonEfficiency->SetRecoParamStorage(recoParamStorage.Data());
muonEfficiency->SetMuonTrackCuts(*trackCuts);
muonEfficiency->SetMuonPtCut(VAR_EFFTASK_PTMIN);
muonEfficiency->UseMCLabel(kTRUE);
muonEfficiency->EnableDisplay(kFALSE);
}
TString addExtraTasks = VAR_EXTRATASKS_CONFIGMACRO;
if (!addExtraTasks.IsNull()) gROOT->ProcessLineSync(TString::Format(".x %s",addExtraTasks.Data()));
if (iESDfilter) {
// ESD filter task configuration.
gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/ESDfilter/macros/AddTaskESDFilter.C");
if (iMUONcopyAOD) {
printf("Registering delta AOD file\n");
mgr->RegisterExtraFile("AliAOD.Muons.root");
}
AliAnalysisTaskESDfilter *taskesdfilter = AddTaskESDFilter(useKFILTER, iMUONcopyAOD, kFALSE, kFALSE /*usePhysicsSelection*/,kFALSE,kTRUE,kTRUE,kTRUE,1500,VAR_MUONMCMODE,kFALSE,kFALSE,kTRUE,kFALSE); // others
taskesdfilter->DisablePmdClusters();
taskesdfilter->DisableCaloClusters();
taskesdfilter->DisableCells();
taskesdfilter->DisableCaloTrigger("PHOS");
taskesdfilter->DisableCaloTrigger("EMCAL");
taskesdfilter->SetPropagateTrackToEMCal(kFALSE);
if ( 0 && VAR_USE_ITS_RECO ) /* 0 for the moment to get this macro running also with AliRoot <= .... */
{
AliAnalysisTaskESDMuonFilter* muFilter = mgr->GetTask("ESD Muon Filter");
if ( !muFilter )
{
std::cout << "ERROR : got a NULL muFilter ! so I cannot ask to keep SPD tracklets !" << std::endl;
}
else
{
muFilter->SetWithSPDtracklets(kTRUE);
}
}
}
}
