1 AliAnalysisTaskSEDvsMultiplicity *AddTaskDvsMultiplicity(Int_t system=0,
5 TString finDirname="Loose",
7 TString finAnObjname="AnalysisCuts",
8 TString estimatorFilename="",
10 Bool_t subtractDau=kFALSE,
11 Bool_t NchWeight=kFALSE,
12 Int_t recoEstimator = AliAnalysisTaskSEDvsMultiplicity::kNtrk10,
13 Int_t MCEstimator = AliAnalysisTaskSEDvsMultiplicity::kEta10,
14 Bool_t isPPbData=kFALSE)
17 // Test macro for the AliAnalysisTaskSE for D+ candidates
18 //Invariant mass histogram and
19 // association with MC truth (using MC info in AOD)
20 // R. Bala, bala@to.infn.it
21 // Get the pointer to the existing analysis manager via the static access method.
22 //==============================================================================
24 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
26 ::Error("AddTaskDvsMultiplicity", "No analysis manager to connect to.");
29 Bool_t stdcuts=kFALSE;
31 if( filename.EqualTo("") ) {
34 filecuts=TFile::Open(filename.Data());
35 if(!filecuts ||(filecuts&& !filecuts->IsOpen())){
36 AliFatal("Input file not found : check your cut object");
42 AliRDHFCuts *analysiscuts=0x0;
47 analysiscuts = new AliRDHFCutsDplustoKpipi();
48 if (system == 0) analysiscuts->SetStandardCutsPP2010();
49 else analysiscuts->SetStandardCutsPbPb2011();
51 else analysiscuts = (AliRDHFCutsDplustoKpipi*)filecuts->Get(finAnObjname);
53 }else if(pdgMeson==421){
55 analysiscuts = new AliRDHFCutsD0toKpi();
56 if (system == 0) analysiscuts->SetStandardCutsPP2010();
57 else analysiscuts->SetStandardCutsPbPb2011();
59 else analysiscuts = (AliRDHFCutsD0toKpi*)filecuts->Get(finAnObjname);
61 }else if(pdgMeson==413){
63 analysiscuts = new AliRDHFCutsDStartoKpipi();
64 if (system == 0) analysiscuts->SetStandardCutsPP2010();
65 else analysiscuts->SetStandardCutsPbPb2011();
67 else analysiscuts = (AliRDHFCutsDStartoKpipi*)filecuts->Get(finAnObjname);
71 AliAnalysisTaskSEDvsMultiplicity *dMultTask = new AliAnalysisTaskSEDvsMultiplicity("dMultAnalysis",pdgMeson,analysiscuts,isPPbData);
72 dMultTask->SetReadMC(readMC);
73 dMultTask->SetDebugLevel(0);
74 dMultTask->SetUseBit(kTRUE);
75 dMultTask->SetDoImpactParameterHistos(kFALSE);
76 dMultTask->SetSubtractTrackletsFromDaughters(subtractDau);
77 dMultTask->SetMultiplicityEstimator(recoEstimator);
78 dMultTask->SetMCPrimariesEstimator(MCEstimator);
79 dMultTask->SetMCOption(MCOption);
80 if(isPPbData) dMultTask->SetIsPPbData();
83 TH1F *hNchPrimaries = (TH1F*)filecuts->Get("hGenPrimaryParticlesInelGt0");
85 dMultTask->UseMCNchWeight(true);
86 dMultTask->SetHistoNchWeight(hNchPrimaries);
88 AliFatal("Histogram for multiplicity weights not found");
94 dMultTask->SetMassLimits(1.5648,2.1648);
95 dMultTask->SetNMassBins(200);
96 }else if(pdgMeson==411)dMultTask->SetMassLimits(pdgMeson,0.2);
98 if(estimatorFilename.EqualTo("") ) {
99 printf("Estimator file not provided, multiplcity corrected histograms will not be filled\n");
101 if ( isPPbData) { //Only use two profiles if pPb
102 const Char_t* periodNames[2] = {"LHC13b", "LHC13c"};
103 TProfile* multEstimatorAvg[2];
104 TFile* fileEstimator=TFile::Open(estimatorFilename.Data());
106 AliFatal("File with multiplicity estimator not found\n");
110 for(Int_t ip=0; ip<2; ip++) {
111 multEstimatorAvg[ip] = (TProfile*)(fileEstimator->Get(Form("SPDmult10_%s",periodNames[ip]))->Clone(Form("SPDmult10_%s_clone",periodNames[ip])));
113 dMultTask->SetMultiplVsZProfileLHC13b(multEstimatorAvg[0]);
114 dMultTask->SetMultiplVsZProfileLHC13c(multEstimatorAvg[1]);
115 dMultTask->SetReferenceMultiplcity(refMult);
119 const Char_t* periodNames[4] = {"LHC10b", "LHC10c", "LHC10d", "LHC10e"};
120 TProfile* multEstimatorAvg[4];
121 TFile* fileEstimator=TFile::Open(estimatorFilename.Data());
123 AliFatal("File with multiplicity estimator not found\n");
127 for(Int_t ip=0; ip<4; ip++) {
128 multEstimatorAvg[ip] = (TProfile*)(fileEstimator->Get(Form("SPDmult10_%s",periodNames[ip]))->Clone(Form("SPDmult10_%s_clone",periodNames[ip])));
130 dMultTask->SetMultiplVsZProfileLHC10b(multEstimatorAvg[0]);
131 dMultTask->SetMultiplVsZProfileLHC10c(multEstimatorAvg[1]);
132 dMultTask->SetMultiplVsZProfileLHC10d(multEstimatorAvg[2]);
133 dMultTask->SetMultiplVsZProfileLHC10e(multEstimatorAvg[3]);
134 dMultTask->SetReferenceMultiplcity(refMult);
137 mgr->AddTask(dMultTask);
139 // Create containers for input/output
141 TString inname = "cinput";
142 TString outname = "coutput";
143 TString cutsname = "coutputCuts";
144 TString normname = "coutputNorm";
145 TString profname = "coutputProf";
147 inname += Name.Data();
148 outname += Name.Data();
149 cutsname += Name.Data();
150 normname += Name.Data();
151 profname += Name.Data();
152 inname += finDirname.Data();
153 outname += finDirname.Data();
154 cutsname += finDirname.Data();
155 normname += finDirname.Data();
156 profname += finDirname.Data();
158 AliAnalysisDataContainer *cinput = mgr->CreateContainer(inname,TChain::Class(),AliAnalysisManager::kInputContainer);
160 TString outputfile = AliAnalysisManager::GetCommonFileName();
161 outputfile += ":PWG3_D2H_DMult_";
162 outputfile += Name.Data();
163 outputfile += finDirname.Data();
165 AliAnalysisDataContainer *coutputCuts = mgr->CreateContainer(cutsname,TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
166 AliAnalysisDataContainer *coutput = mgr->CreateContainer(outname,TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
167 AliAnalysisDataContainer *coutputNorm = mgr->CreateContainer(normname,TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
168 AliAnalysisDataContainer *coutputProf = mgr->CreateContainer(profname,TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
170 mgr->ConnectInput(dMultTask,0,mgr->GetCommonInputContainer());
172 mgr->ConnectOutput(dMultTask,1,coutput);
174 mgr->ConnectOutput(dMultTask,2,coutputCuts);
176 mgr->ConnectOutput(dMultTask,3,coutputNorm);
178 mgr->ConnectOutput(dMultTask,4,coutputProf);