1 //=========================================================================//
3 // Analysis AddTask for Particle Ratio Fluctuation Study //
4 // Author: Deepika Rathee || Satyajit Jenara //
5 // drathee@cern.ch || sjena@cern.ch
6 // Thu Jun 19 11:44:51 CEST 2014
8 //=========================================================================//
10 void AddAliEbyEPidBayesTTaskMC(Double_t ptl = 0,
11 Int_t AODfilterBit = 768,
13 Int_t requestTofPid = 0,
14 Double_t nSigmaCut = 5.,
15 TString centralityEstimator = "V0M",
16 TString ctaskname = "2010") {
20 TString taskname = "CfHmBayes_";
21 taskname += ctaskname;
23 taskname += Form("%d",isMC);
25 taskname += centralityEstimator;
27 taskname += Form("%d",AODfilterBit);
30 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
32 ::Error("AddTaskFluctuations", "No analysis manager to connect to.");
36 if (!mgr->GetInputEventHandler()) {
37 ::Error("AddTaskFluctuations", "This task requires an input event handler");
40 TString type = mgr->GetInputEventHandler()->GetDataType();
41 TString basefilename = AliAnalysisManager::GetCommonFileName();
43 AliHelperPID* help = new AliHelperPID();
44 help->SetNSigmaCut(nSigmaCut);
45 help->SetPIDType(pidtype);
46 if (requestTofPid) help->SetfRequestTOFPID(requestTofPid);
47 if (ptl != 0 ) help->SetfPtTOFPID(ptl);
48 //help->SetisMC(isMC);
51 AliPIDCombined *pidc=new AliPIDCombined();
52 pidc->SetDefaultTPCPriors();
53 help->SetPIDCombined(pidc);
57 /* AliHelperPID *pid=new AliHelperPID();
58 pid->SetName(Form("HelperPID"));
61 pid->SetfPtTOFPID(0.3);
62 AliPIDCombined *pidc=new AliPIDCombined();
63 pidc->SetDefaultTPCPriors();
64 pid->SetPIDCombined(pidc);
71 AliEbyEPidTTaskMC *taskqa;
72 taskqa = new AliEbyEPidTTaskMC(taskname.Data());
73 taskqa->SetAODtrackCutBit(AODfilterBit);
74 taskqa->SetCentralityEstimator(centralityEstimator.Data());
75 taskqa->SetAnalysisType(isMC);
76 taskqa->SetHelperPID(help);
79 AliAnalysisDataContainer *coutqa = mgr->CreateContainer(Form("%s_QA",taskname.Data()),TList::Class(),AliAnalysisManager::kOutputContainer, Form("%s",basefilename.Data()));
80 AliAnalysisDataContainer *coutt = mgr->CreateContainer("fEventTree",TTree::Class(),AliAnalysisManager::kOutputContainer, Form("%s",basefilename.Data()));
82 mgr->ConnectInput(taskqa, 0, mgr->GetCommonInputContainer());
83 mgr->ConnectOutput(taskqa, 1, coutqa);
84 mgr->ConnectOutput(taskqa, 2, coutt);