1 AliAnalysisTaskMuonDistributions *AddTaskMuonDistributions(const char *kAnalysisType, Bool_t doFitInvMass){
3 //****************************************************************************************
5 // The attached class prepares and draws some kinematical distributions of muons/dimuons
7 //****************************************************************************************
9 printf("Creating Task for Muon/Dimuon Histos\n");
11 // Get the pointer to the existing analysis manager via the static access method.
12 //==============================================================================
13 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
15 ::Error("AddTaskMuonDistributions", "No analysis manager to connect to.");
19 TString outputfile = AliAnalysisManager::GetCommonFileName();
20 outputfile += ":PWG3Muon_Dimuon";
22 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("Dimuon",TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
24 AliAnalysisTaskMuonDistributions *MuonDistributionsTask = new AliAnalysisTaskMuonDistributions("AliAnalysisTaskMuonDistributions");
25 MuonDistributionsTask->SetAnalysisType(kAnalysisType);
27 // define by hand the beam energy
29 MuonDistributionsTask->SetBeamEnergy(5000.);
33 MuonDistributionsTask->SetInvMassFitLimits(2.,5.5);
34 MuonDistributionsTask->SetPsiFitLimits(2.9,3.3);
35 MuonDistributionsTask->SetPsiPFitLimits(3.3,4.2);
36 MuonDistributionsTask->SetBckFitLimits(2.,2.8);
38 // perform fit to the invariant mass spectrum
40 MuonDistributionsTask->FitInvariantMassSpectrum(doFitInvMass);
42 mgr->AddTask(MuonDistributionsTask);
44 mgr->ConnectInput(MuonDistributionsTask,0,mgr->GetCommonInputContainer());
45 mgr->ConnectOutput(MuonDistributionsTask,1,coutput1);
47 return MuonDistributionsTask;