Int_t AddPWG3MuonTrain(Int_t iESDAnalysis=1,
                       Int_t iAODAnalysis=0,
                       Int_t addMuonDistributions=0,
		       Int_t addSingleMuonAnalysis=0,
		       Int_t addMuonHFAnalysis=0,
		       Int_t addDimuonCFContainer=0,
		       Bool_t usePhysicsSelection=kFALSE) {

  // Analysis wagons for PWG3Muon (Roberta)
  
  TString taskName="",loadMacroPath="$ALICE_ROOT/PWG3/muon/";
  Int_t ntasks=0;
  
  if(addMuonDistributions) {
    taskName="AddTaskMuonDistributions.C"; taskName.Prepend(loadMacroPath.Data());
    gROOT->LoadMacro(taskName.Data());
    Bool_t doInvMassFit = kTRUE;
    if(iESDAnalysis){
      AliAnalysisTaskMuonDistributions *esdmuondistributionstask = AddTaskMuonDistributions("ESD",doInvMassFit);	
      if(usePhysicsSelection)esdmuondistributionstask->SelectCollisionCandidates();
    } else if(iAODAnalysis){
      AliAnalysisTaskMuonDistributions *aodmuondistributionstask = AddTaskMuonDistributions("AOD",doInvMassFit);	
    }
    ntasks++;
  }
  

  if(addSingleMuonAnalysis) {
    taskName="AddTaskSingleMuonAnalysis.C"; taskName.Prepend(loadMacroPath.Data());
    gROOT->LoadMacro(taskName.Data());
    AliAnalysisTaskSingleMu *singlemutask = AddTaskSingleMuonAnalysis();	
    ntasks++;
  }

  if(addMuonHFAnalysis) {
    taskName="AddTaskMuonsHF.C"; taskName.Prepend(loadMacroPath.Data());
    gROOT->LoadMacro(taskName.Data());
    Int_t runMode = 0;
    Bool_t isAOD;
    if(iAODAnalysis==1) isAOD=kTRUE;
    else isAOD=kFALSE;
    Bool_t isTree = kFALSE;
    Bool_t isMC   = kFALSE;
    AliAnalysisTaskSEMuonsHF *muonhftask = AddTaskMuonsHF(runMode, isMC, isTree);	
    ntasks++;
  }

  if(addDimuonCFContainer) {
    taskName="AddTaskDimuonCFContainerBuilder.C"; taskName.Prepend(loadMacroPath.Data());
    gROOT->LoadMacro(taskName.Data());
    Int_t runMode = 0;
    Bool_t isAOD;
    if(iAODAnalysis==1) isAOD=kTRUE;
    else isAOD=kFALSE;
    Bool_t isAcceptance = kFALSE;
    Bool_t readMC   = kFALSE;
    Double_t ebeam = 3500.;
    AliAnalysisTaskDimuonCFContainerBuilder *dimuonCFtask =
    AddTaskDimuonCFContainerBuilder(isAOD,readMC,isAcceptance,ebeam);
    ntasks++;
  }

  return ntasks;
}