-AliAnalysisTaskSpectraAllChAOD* AddTaskSpectraAllChAOD(\r
- Bool_t mc=kFALSE,\r
- Double_t CentCutMin=0,\r
- Double_t CentCutMax=100,\r
- Double_t QvecCutMin=0,\r
- Double_t QvecCutMax=100,\r
- Double_t EtaMin=-0.8,\r
- Double_t EtaMax=0.8,\r
- Double_t pt=50.,\r
- Double_t ptTofMatch=.6,\r
- UInt_t trkbit=1,\r
- Double_t DCA=100000,\r
- UInt_t minNclsTPC=70,\r
- Double_t nsigmacut=5.,\r
- Int_t PIDtype=3,\r
- TString opt=""){\r
- \r
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
- if (!mgr) \r
- {\r
- ::Error("AddAliAnalysisTaskSpectraAllChAOD", "No analysis manager to connect to.");\r
- return NULL;\r
- } \r
- \r
- // Check the analysis type using the event handlers connected to the analysis manager.\r
- //==============================================================================\r
- if (!mgr->GetInputEventHandler()) \r
- {\r
- ::Error("AliAnalysisTaskSpectraAllChAOD", "This task requires an input event handler");\r
- return NULL;\r
- } \r
- \r
- TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
- if(type.Contains("ESD"))\r
- {\r
- ::Error("AliAnalysisTaskSpectraAllChAOD", "This task requires to run on AOD");\r
- return NULL;\r
- }\r
- \r
- AliSpectraAODTrackCuts * trcuts = new AliSpectraAODTrackCuts(Form("TrackCuts%s",opt.Data())); \r
- trcuts->SetDCA(DCA);\r
- trcuts->SetTrackBits(trkbit);\r
- trcuts->SetPt(pt);\r
- trcuts->SetPtTOFMatching(ptTofMatch); \r
- trcuts->SetEta(EtaMin,EtaMax);\r
- trcuts->SetMinTPCcls(minNclsTPC);\r
- trcuts->PrintCuts();\r
- \r
- AliSpectraAODEventCuts * evcuts = new AliSpectraAODEventCuts(Form("EventCuts%s",opt.Data()));\r
- evcuts->SetQVectorCut(QvecCutMin,QvecCutMax);\r
- evcuts->SetCentralityCutMax(CentCutMax); \r
- evcuts->SetCentralityCutMin(CentCutMin);\r
- if(mc==1)evcuts->SetIsMC(kTRUE);\r
- evcuts->PrintCuts();\r
- \r
- AliHelperPID *pid=new AliHelperPID();\r
- pid->SetName(Form("HelperPID%s",opt.Data()));\r
- pid->SetNSigmaCut(nsigmacut);\r
- pid->SetPIDType(PIDtype);\r
- if(PIDtype==3){\r
- AliPIDCombined *pidc=new AliPIDCombined();\r
- pidc->SetDefaultTPCPriors();\r
- pid->SetPIDCombined(pidc);\r
- }\r
- \r
- AliAnalysisTaskSpectraAllChAOD *task = new AliAnalysisTaskSpectraAllChAOD(Form("TaskAODSpectraCent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f_TrBit%d%s", \r
- CentCutMin,\r
- CentCutMax,\r
- QvecCutMin,\r
- QvecCutMax,\r
- EtaMin,\r
- EtaMax,\r
- trkbit,\r
- opt.Data()));\r
- task->SetEventCuts(evcuts);\r
- task->SetTrackCuts(trcuts);\r
- task->SetHelperPID(pid);\r
- if(mc==1)task->SetIsMC(kTRUE);\r
- \r
- TString outputFileName = AliAnalysisManager::GetCommonFileName();\r
- \r
- TString typeofdata=mc?"MC":"Data";\r
- \r
- outputFileName += Form(":SpectraESE_%s%s",typeofdata.Data(),opt.Data());\r
- \r
- cout<<"outputFileName: "<<outputFileName<<endl;\r
- AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer(); \r
- AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chist%s",opt.Data()), TList::Class(), AliAnalysisManager::kOutputContainer,outputFileName);\r
- AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcut%s",opt.Data()), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,outputFileName);\r
- AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcut%s",opt.Data()), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer, outputFileName);\r
- AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpid%s",opt.Data()), AliHelperPID::Class(), AliAnalysisManager::kOutputContainer, outputFileName);\r
- mgr->AddTask(task);\r
- \r
- mgr->ConnectInput(task, 0, cinput);\r
- mgr->ConnectOutput(task, 1, coutputpt1);\r
- mgr->ConnectOutput(task, 2, coutputpt2);\r
- mgr->ConnectOutput(task, 3, coutputpt3);\r
- mgr->ConnectOutput(task, 4, coutputpt4);\r
- \r
- return task;\r
-}\r
+AliAnalysisTaskSpectraAllChAOD* AddTaskSpectraAllChAOD(
+ Bool_t mc=kFALSE,
+ Double_t CentCutMin=0,
+ Double_t CentCutMax=100,
+ Double_t QvecCutMin=0,
+ Double_t QvecCutMax=100,
+ Double_t EtaMin=-0.8,
+ Double_t EtaMax=0.8,
+ Double_t pt=50.,
+ Double_t ptTofMatch=.6,
+ UInt_t trkbit=1,
+ Double_t DCA=100000,
+ UInt_t minNclsTPC=70,
+ Double_t nsigmacut=5.,
+ Int_t PIDtype=3,
+ TString opt=""){
+
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ if (!mgr)
+ {
+ ::Error("AddAliAnalysisTaskSpectraAllChAOD", "No analysis manager to connect to.");
+ return NULL;
+ }
+
+ // Check the analysis type using the event handlers connected to the analysis manager.
+ //==============================================================================
+ if (!mgr->GetInputEventHandler())
+ {
+ ::Error("AliAnalysisTaskSpectraAllChAOD", "This task requires an input event handler");
+ return NULL;
+ }
+
+ TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
+ if(type.Contains("ESD"))
+ {
+ ::Error("AliAnalysisTaskSpectraAllChAOD", "This task requires to run on AOD");
+ return NULL;
+ }
+
+ AliSpectraAODTrackCuts * trcuts = new AliSpectraAODTrackCuts(Form("TrackCuts%s",opt.Data()));
+ trcuts->SetDCA(DCA);
+ trcuts->SetTrackBits(trkbit);
+ trcuts->SetPt(pt);
+ trcuts->SetPtTOFMatching(ptTofMatch);
+ trcuts->SetEta(EtaMin,EtaMax);
+ trcuts->SetMinTPCcls(minNclsTPC);
+ trcuts->PrintCuts();
+
+ AliSpectraAODEventCuts * evcuts = new AliSpectraAODEventCuts(Form("EventCuts%s",opt.Data()));
+ evcuts->SetQVectorCut(QvecCutMin,QvecCutMax);
+ evcuts->SetCentralityCutMax(CentCutMax);
+ evcuts->SetCentralityCutMin(CentCutMin);
+ if(mc==1)evcuts->SetIsMC(kTRUE);
+ evcuts->PrintCuts();
+
+ AliHelperPID *pid=new AliHelperPID();
+ pid->SetName(Form("HelperPID%s",opt.Data()));
+ pid->SetNSigmaCut(nsigmacut);
+ pid->SetPIDType(PIDtype);
+ if(PIDtype==3){
+ AliPIDCombined *pidc=new AliPIDCombined();
+ pidc->SetDefaultTPCPriors();
+ pid->SetPIDCombined(pidc);
+ }
+
+ AliAnalysisTaskSpectraAllChAOD *task = new AliAnalysisTaskSpectraAllChAOD(Form("TaskAODSpectraCent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f_TrBit%d%s",
+ CentCutMin,
+ CentCutMax,
+ QvecCutMin,
+ QvecCutMax,
+ EtaMin,
+ EtaMax,
+ trkbit,
+ opt.Data()));
+ task->SetEventCuts(evcuts);
+ task->SetTrackCuts(trcuts);
+ task->SetHelperPID(pid);
+ if(mc==1)task->SetIsMC(kTRUE);
+
+ TString outputFileName = AliAnalysisManager::GetCommonFileName();
+
+ TString typeofdata=mc?"MC":"Data";
+
+ outputFileName += Form(":SpectraESE_%s%s",typeofdata.Data(),opt.Data());
+
+ cout<<"outputFileName: "<<outputFileName<<endl;
+ AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
+ AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chist%s",opt.Data()), TList::Class(), AliAnalysisManager::kOutputContainer,outputFileName);
+ AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcut%s",opt.Data()), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,outputFileName);
+ AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcut%s",opt.Data()), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer, outputFileName);
+ AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpid%s",opt.Data()), AliHelperPID::Class(), AliAnalysisManager::kOutputContainer, outputFileName);
+ mgr->AddTask(task);
+
+ mgr->ConnectInput(task, 0, cinput);
+ mgr->ConnectOutput(task, 1, coutputpt1);
+ mgr->ConnectOutput(task, 2, coutputpt2);
+ mgr->ConnectOutput(task, 3, coutputpt3);
+ mgr->ConnectOutput(task, 4, coutputpt4);
+
+ return task;
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff