]>
Commit | Line | Data |
---|---|---|
1 | AliAnalysisTaskMuonDistributions *AddTaskMuonDistributions(const char *kAnalysisType, Bool_t doFitInvMass){ | |
2 | ||
3 | //**************************************************************************************** | |
4 | // Add task class. | |
5 | // The attached class prepares and draws some kinematical distributions of muons/dimuons | |
6 | // Roberta | |
7 | //**************************************************************************************** | |
8 | ||
9 | printf("Creating Task for Muon/Dimuon Histos\n"); | |
10 | ||
11 | // Get the pointer to the existing analysis manager via the static access method. | |
12 | //============================================================================== | |
13 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); | |
14 | if (!mgr) { | |
15 | ::Error("AddTaskMuonDistributions", "No analysis manager to connect to."); | |
16 | return NULL; | |
17 | } | |
18 | ||
19 | TString outputfile = AliAnalysisManager::GetCommonFileName(); | |
20 | outputfile += ":PWG3Muon_Dimuon"; | |
21 | ||
22 | AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("Dimuon",TList::Class(),AliAnalysisManager::kOutputContainer,outputfile); | |
23 | ||
24 | AliAnalysisTaskMuonDistributions *MuonDistributionsTask = new AliAnalysisTaskMuonDistributions("AliAnalysisTaskMuonDistributions"); | |
25 | MuonDistributionsTask->SetAnalysisType(kAnalysisType); | |
26 | // | |
27 | // define by hand the beam energy | |
28 | // | |
29 | MuonDistributionsTask->SetBeamEnergy(5000.); | |
30 | // | |
31 | // define fits limits | |
32 | // | |
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); | |
37 | // | |
38 | // perform fit to the invariant mass spectrum | |
39 | // | |
40 | MuonDistributionsTask->FitInvariantMassSpectrum(doFitInvMass); | |
41 | ||
42 | mgr->AddTask(MuonDistributionsTask); | |
43 | ||
44 | mgr->ConnectInput(MuonDistributionsTask,0,mgr->GetCommonInputContainer()); | |
45 | mgr->ConnectOutput(MuonDistributionsTask,1,coutput1); | |
46 | ||
47 | return MuonDistributionsTask; | |
48 | } |