1 AliAnalysisTaskSpectraAllChAOD* AddTaskSpectraAllChAOD(
4 Double_t CentCutMax=100,
6 Double_t QvecCutMax=100,
10 Double_t ptTofMatch=.6,
14 Double_t nsigmacut=5.,
18 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
21 ::Error("AddAliAnalysisTaskSpectraAllChAOD", "No analysis manager to connect to.");
25 // Check the analysis type using the event handlers connected to the analysis manager.
26 //==============================================================================
27 if (!mgr->GetInputEventHandler())
29 ::Error("AliAnalysisTaskSpectraAllChAOD", "This task requires an input event handler");
33 TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
34 if(type.Contains("ESD"))
36 ::Error("AliAnalysisTaskSpectraAllChAOD", "This task requires to run on AOD");
40 AliSpectraAODTrackCuts * trcuts = new AliSpectraAODTrackCuts(Form("TrackCuts%s",opt.Data()));
42 trcuts->SetTrackBits(trkbit);
44 trcuts->SetPtTOFMatching(ptTofMatch);
45 trcuts->SetEta(EtaMin,EtaMax);
46 trcuts->SetMinTPCcls(minNclsTPC);
49 AliSpectraAODEventCuts * evcuts = new AliSpectraAODEventCuts(Form("EventCuts%s",opt.Data()));
50 evcuts->SetQVectorCut(QvecCutMin,QvecCutMax);
51 evcuts->SetCentralityCutMax(CentCutMax);
52 evcuts->SetCentralityCutMin(CentCutMin);
53 if(mc==1)evcuts->SetIsMC(kTRUE);
56 AliHelperPID *pid=new AliHelperPID();
57 pid->SetName(Form("HelperPID%s",opt.Data()));
58 pid->SetNSigmaCut(nsigmacut);
59 pid->SetPIDType(PIDtype);
60 pid->SetfPtTOFPID(ptTofMatch);
62 AliPIDCombined *pidc=new AliPIDCombined(Form("PIDCombined%s",opt.Data()),Form("PIDCombined%s",opt.Data()));
63 pidc->SetDefaultTPCPriors();
64 pid->SetPIDCombined(pidc);//not implemented
67 AliAnalysisTaskSpectraAllChAOD *task = new AliAnalysisTaskSpectraAllChAOD(Form("TaskAODSpectraCent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f_TrBit%d%s",
76 task->SetEventCuts(evcuts);
77 task->SetTrackCuts(trcuts);
78 task->SetHelperPID(pid);
79 if(mc==1)task->SetIsMC(kTRUE);
81 TString outputFileName = AliAnalysisManager::GetCommonFileName();
83 TString typeofdata=mc?"MC":"Data";
85 outputFileName += Form(":SpectraESE_%s%s",typeofdata.Data(),opt.Data());
87 cout<<"outputFileName: "<<outputFileName<<endl;
88 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
89 AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chist%s",opt.Data()), TList::Class(), AliAnalysisManager::kOutputContainer,outputFileName);
90 AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcut%s",opt.Data()), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,outputFileName);
91 AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcut%s",opt.Data()), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer, outputFileName);
92 AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpid%s",opt.Data()), AliHelperPID::Class(), AliAnalysisManager::kOutputContainer, outputFileName);
95 mgr->ConnectInput(task, 0, cinput);
96 mgr->ConnectOutput(task, 1, coutputpt1);
97 mgr->ConnectOutput(task, 2, coutputpt2);
98 mgr->ConnectOutput(task, 3, coutputpt3);
99 mgr->ConnectOutput(task, 4, coutputpt4);