1 /* created by fbellini@cern.ch on 29/04/2013 */
2 /* last modified by fbellini on 29/04/2013 */
4 // UInt_t kTriggerMuonAll = AliVEvent::kMUL7 | AliVEvent::kMUSH7 | AliVEvent::kMUU7 | AliVEvent::kMUS7;
5 // UInt_t kTriggerMuonBarell = AliVEvent::kMUU7;
6 // UInt_t kTriggerEMC = AliVEvent::kEMC7;
7 // UInt_t kTriggerHM = AliVEvent::kHighMult;
8 UInt_t kTriggerInt = AliVEvent::kAnyINT;
9 UInt_t kTriggerMask = kTriggerInt;
11 AliAnalysisTaskSE * AddTaskTOFqaID(UInt_t triggerMask = kTriggerMask,
12 Bool_t flagEnableAdvancedCheck=kFALSE,
13 Bool_t useStdCuts2011 = kTRUE,
15 Short_t absPdgCode = 0)
17 // Task for checking TOF QA
19 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
21 ::Error("AddTask", "No analysis manager to connect to.");
25 // Check the analysis type using the event handlers connected to the analysis manager.
26 if (!mgr->GetInputEventHandler()) {
27 ::Error("AddTask", "This task requires an input event handler");
30 TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
32 if (isMC && (absPdgCode<0)) {
33 ::Error("AddTask","Invalid selected PDG code for MC. Set absPdgCode=0 for no specie selection.");
37 AliAnalysisTaskTOFqaID *task = new AliAnalysisTaskTOFqaID(Form("taskTOFqaID_%i",absPdgCode));
38 task->EnableAdvancedCheck(flagEnableAdvancedCheck);
39 task->SetSelectMCspecies(isMC, absPdgCode);
40 task->SelectCollisionCandidates(triggerMask);
41 //AliLog::SetClassDebugLevel("AliAnalysisTaskTOFqaID",4);
44 /* cuts used for QA in 2010 p-p */
46 AliESDtrackCuts* esdTrackCutsLoose2010 = new AliESDtrackCuts("AliESDtrackCuts", "esdTrackCutsLoose2010");
47 esdTrackCutsLoose2010->SetMinNClustersTPC(70);
48 esdTrackCutsLoose2010->SetMaxChi2PerClusterTPC(3.5);
49 esdTrackCutsLoose2010->SetMaxCovDiagonalElements(2, 2, 0.5, 0.5, 2);
50 esdTrackCutsLoose2010->SetRequireTPCRefit(kTRUE);
51 esdTrackCutsLoose2010->SetMaxDCAToVertexXY(3.0);
52 esdTrackCutsLoose2010->SetMaxDCAToVertexZ(3.0);
53 esdTrackCutsLoose2010->SetRequireSigmaToVertex(kTRUE);
54 esdTrackCutsLoose2010->SetAcceptKinkDaughters(kFALSE);
55 esdTrackCutsLoose2010->SetMaxNsigmaToVertex(4.0);
58 /* standard cuts ITS-TPC 2010 */
59 AliESDtrackCuts* esdTrackCutsStd2010 = new AliESDtrackCuts("AliESDtrackCuts", "Standard2010");
60 esdTrackCutsStd2010->GetStandardITSTPCTrackCuts2010(kTRUE,0);
62 // esdTrackCutsStd2010->SetMinNClustersTPC(70);
63 // esdTrackCutsStd2010->SetMaxChi2PerClusterTPC(4);
64 // esdTrackCutsStd2010->SetAcceptKinkDaughters(kFALSE);
65 // esdTrackCutsStd2010->SetRequireTPCRefit(kTRUE);
67 // esdTrackCutsStd2010->SetRequireITSRefit(kTRUE);
68 // esdTrackCutsStd2010->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
69 // AliESDtrackCuts::kAny);
70 // esdTrackCutsStd2010->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01");//selects primaries
71 // esdTrackCutsStd2010->SetMaxDCAToVertexZ(2);
72 // esdTrackCutsStd2010->SetDCAToVertex2D(kFALSE);
73 // esdTrackCutsStd2010->SetRequireSigmaToVertex(kFALSE);
76 /* standard cuts ITS-TPC 2011 */
77 AliESDtrackCuts* esdTrackCutsStd2011 = new AliESDtrackCuts("AliESDtrackCuts", "Standard2011");
78 esdTrackCutsStd2011->GetStandardITSTPCTrackCuts2011(kTRUE,0);
80 AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
82 trackFilter->AddCuts(esdTrackCutsStd2011);
84 trackFilter->AddCuts(esdTrackCutsStd2010);
85 task->SetTrackFilter(trackFilter);
87 TString partName(task->GetSpeciesName(absPdgCode));
89 // Create containers for input/output
90 AliAnalysisDataContainer *cInputTOFqa = mgr->CreateContainer("cInputTOFqa",TChain::Class(),AliAnalysisManager::kInputContainer);
91 AliAnalysisDataContainer *cGeneralTOFqa = mgr->CreateContainer(Form("base_%s",partName.Data()),TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:TOF_Performance",mgr->GetCommonFileName(), partName.Data()));
92 AliAnalysisDataContainer *cTimeZeroTOFqa = mgr->CreateContainer(Form("timeZero_%s",partName.Data()),TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:TOF_Performance",mgr->GetCommonFileName(),partName.Data()));
93 AliAnalysisDataContainer *cPIDTOFqa = mgr->CreateContainer(Form("pid_%s",partName.Data()),TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:TOF_Performance",mgr->GetCommonFileName(),partName.Data()));
94 AliAnalysisDataContainer *cTRDcheckTOFqa = mgr->CreateContainer(Form("trd_%s",partName.Data()),TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:TOF_Performance",mgr->GetCommonFileName(),partName.Data()));
95 AliAnalysisDataContainer *cTriggerTOFqa = mgr->CreateContainer(Form("trigger_%s",partName.Data()),TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:TOF_Performance",mgr->GetCommonFileName(),partName.Data()));
98 mgr->ConnectInput(task, 0,mgr->GetCommonInputContainer());
99 mgr->ConnectOutput(task, 1, cGeneralTOFqa);
100 mgr->ConnectOutput(task, 2, cTimeZeroTOFqa);
101 mgr->ConnectOutput(task, 3, cPIDTOFqa);
102 mgr->ConnectOutput(task, 4, cTRDcheckTOFqa);
103 mgr->ConnectOutput(task, 5, cTriggerTOFqa);