3 Bool_t AddTrackCutsLHC10h(AliAnalysisTaskESDfilter* esdFilter);
4 Bool_t AddTrackCutsLHC11h(AliAnalysisTaskESDfilter* esdFilter);
5 Bool_t enableTPCOnlyAODTracks = kTRUE;
8 AliAnalysisTaskESDfilter *AddTaskESDFilterPWG4Train(Bool_t useKineFilter=kTRUE,
9 Bool_t writeMuonAOD=kFALSE,
10 Bool_t writeDimuonAOD=kFALSE,
11 Bool_t usePhysicsSelection=kFALSE,
12 Bool_t useCentralityTask=kFALSE)
14 // Creates a filter task and adds it to the analysis manager.
16 // Get the pointer to the existing analysis manager via the static access method.
17 //==============================================================================
18 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
20 ::Error("AddTaskESDFilter", "No analysis manager to connect to.");
24 // This task requires an ESD input handler and an AOD output handler.
25 // Check this using the analysis manager.
26 //===============================================================================
27 TString type = mgr->GetInputEventHandler()->GetDataType();
28 if (!type.Contains("ESD")) {
29 ::Error("AddTaskESDFilter", "ESD filtering task needs the manager to have an ESD input handler.");
32 // Check if AOD output handler exist.
33 AliAODHandler *aod_h = (AliAODHandler*)mgr->GetOutputEventHandler();
35 ::Error("AddTaskESDFilter", "ESD filtering task needs the manager to have an AOD output handler.");
38 // Check if MC handler is connected in case kine filter requested
39 AliMCEventHandler *mcH = (AliMCEventHandler*)mgr->GetMCtruthEventHandler();
40 if (!mcH && useKineFilter) {
41 ::Error("AddTaskESDFilter", "No MC handler connected while kine filtering requested");
45 // Create the task, add it to the manager and configure it.
46 //===========================================================================
47 // Barrel tracks filter
48 AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
50 // Make the AOD a little bit lighter and filtering faster
52 esdfilter->DisableCascades();
53 // esdfilter->DisableV0s();
54 esdfilter->DisableKinks();
55 // esdfilter->DisableTracks();
56 esdfilter->DisablePmdClusters();
57 // esdfilter->DisableCaloClusters();
58 // esdfilter->DisableCells();
59 esdfilter->DisableTracklets();
60 esdfilter->SetWriteHybridGlobalConstrainedOnly(kTRUE);
62 mgr->AddTask(esdfilter);
64 if(usePhysicsSelection){
65 esdfilter->SelectCollisionCandidates(AliVEvent::kAny);
68 // Filtering of MC particles (decays conversions etc)
69 // this task has to go AFTER all other filter tasks
70 // since it fills the AODMC array with all
71 // selected MC Particles, only this way we have the
72 // AODMCparticle information available for following tasks
73 AliAnalysisTaskMCParticleFilter *kinefilter = 0;
75 kinefilter = new AliAnalysisTaskMCParticleFilter("Particle Kine Filter");
76 if(usePhysicsSelection)kinefilter->SelectCollisionCandidates(AliVEvent::kAny);
77 mgr->AddTask(kinefilter);
81 // Filter with cuts on V0s
82 AliESDv0Cuts* esdV0Cuts = new AliESDv0Cuts("Standard V0 Cuts pp", "ESD V0 Cuts");
83 esdV0Cuts->SetMinRadius(0.2);
84 esdV0Cuts->SetMaxRadius(200);
85 esdV0Cuts->SetMinDcaPosToVertex(0.05);
86 esdV0Cuts->SetMinDcaNegToVertex(0.05);
87 esdV0Cuts->SetMaxDcaV0Daughters(1.5);
88 esdV0Cuts->SetMinCosinePointingAngle(0.99);
89 AliAnalysisFilter* v0Filter = new AliAnalysisFilter("v0Filter");
90 v0Filter->AddCuts(esdV0Cuts);
93 Bool_t bSuccess = false;
94 TString runPeriod = AliAnalysisManager::GetGlobalStr("kJetRunPeriod",bSuccess);
96 if(bSuccess&&runPeriod.Contains("LHC11h")){
97 AddTrackCutsLHC11h(esdfilter);
100 Printf("%s%d: Creating Default track cuts 10h",(char*)__FILE__,__LINE__);
101 AddTrackCutsLHC10h(esdfilter);
104 esdfilter->SetV0Filter(v0Filter);
106 // Enable writing of Muon AODs
107 // esdmuonfilter->SetWriteMuonAOD(writeMuonAOD);
109 // Enable writing of Dimuon AODs
110 // esdmuonfilter->SetWriteDimuonAOD(writeDimuonAOD);
112 // Create ONLY the output containers for the data produced by the task.
113 // Get and connect other common input/output containers via the manager as below
114 //==============================================================================
116 mgr->ConnectInput (esdfilter, 0, mgr->GetCommonInputContainer());
117 mgr->ConnectOutput (esdfilter, 0, mgr->GetCommonOutputContainer());
120 // mgr->ConnectInput (esdmuonfilter, 0, mgr->GetCommonInputContainer());
122 mgr->ConnectInput (kinefilter, 0, mgr->GetCommonInputContainer());
123 mgr->ConnectOutput (kinefilter, 0, mgr->GetCommonOutputContainer());
124 AliAnalysisDataContainer *coutputEx = mgr->CreateContainer("cFilterList", TList::Class(),
125 AliAnalysisManager::kOutputContainer,"pyxsec_hists.root");
126 mgr->ConnectOutput (kinefilter, 1,coutputEx);
128 if (useCentralityTask) {
129 mgr->ConnectInput (ctask, 0, mgr->GetCommonInputContainer());
130 mgr->ConnectOutput(ctask, 0, mgr->GetCommonOutputContainer());
137 Bool_t AddTrackCutsLHC10h(AliAnalysisTaskESDfilter* esdfilter){
139 Printf("%s%d: Creating Track Cuts for LH10h",(char*)__FILE__,__LINE__);
141 // Cuts on primary tracks
142 AliESDtrackCuts* esdTrackCutsL = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
144 // ITS stand-alone tracks
145 AliESDtrackCuts* esdTrackCutsITSsa = new AliESDtrackCuts("ITS stand-alone Track Cuts", "ESD Track Cuts");
146 esdTrackCutsITSsa->SetRequireITSStandAlone(kTRUE);
148 // Pixel OR necessary for the electrons
149 AliESDtrackCuts *itsStrong = new AliESDtrackCuts("ITSorSPD", "pixel requirement for ITS");
150 itsStrong->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
153 // PID for the electrons
154 AliESDpidCuts *electronID = new AliESDpidCuts("Electrons", "Electron PID cuts");
155 electronID->SetTPCnSigmaCut(AliPID::kElectron, 3.);
157 // tighter cuts on primary particles for high pT tracks
158 // take the standard cuts, which include already
159 // ITSrefit and use only primaries...
161 // ITS cuts for new jet analysis
162 gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/CreateTrackCutsPWG4.C");
163 AliESDtrackCuts* esdTrackCutsHG0 = CreateTrackCutsPWG4(10001006);
165 // throw out tracks with too low number of clusters in
166 // the first pass (be consistent with TPC only tracks)
167 // N.B. the number off crossed rows still acts on the tracks after
168 // all iterations if we require tpc standalone, number of clusters
169 // and chi2 TPC cuts act on track after the first iteration
170 // esdTrackCutsH0->SetRequireTPCStandAlone(kTRUE);
171 // esdTrackCutsH0->SetMinNClustersTPC(80); // <--- first pass
174 // the complement to the one with SPD requirement
175 AliESDtrackCuts* esdTrackCutsHG1 = CreateTrackCutsPWG4(10011006);
177 // the tracks that must not be taken pass this cut and
179 AliESDtrackCuts* esdTrackCutsHG2 = CreateTrackCutsPWG4(10021006);
181 // standard cuts also used in R_AA analysis
182 AliESDtrackCuts* esdTrackCutsH2 = CreateTrackCutsPWG4(1000);
184 AliESDtrackCuts* esdTrackCutsGCOnly = CreateTrackCutsPWG4(10041006);
187 AliESDtrackCuts* esdTrackCutsTPCCOnly = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
188 esdTrackCutsTPCCOnly->SetMinNClustersTPC(70);
190 // Compose the filter
191 AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
193 trackFilter->AddCuts(esdTrackCutsL);
195 trackFilter->AddCuts(esdTrackCutsITSsa);
197 trackFilter->AddCuts(itsStrong);
198 itsStrong->SetFilterMask(1); // AND with Standard track cuts
200 trackFilter->AddCuts(electronID);
201 electronID->SetFilterMask(4); // AND with Pixel Cuts
203 trackFilter->AddCuts(esdTrackCutsHG0);
205 trackFilter->AddCuts(esdTrackCutsHG1);
207 trackFilter->AddCuts(esdTrackCutsHG2);
209 AliESDtrackCuts* esdTrackCutsHG0_tmp = new AliESDtrackCuts(*esdTrackCutsHG0); // avoid double delete
210 trackFilter->AddCuts(esdTrackCutsHG0_tmp); // add once more for tpc only tracks
212 trackFilter->AddCuts(esdTrackCutsGCOnly);
214 AliESDtrackCuts* esdTrackCutsHG1_tmp = new AliESDtrackCuts(*esdTrackCutsHG1); // avoid double delete
215 trackFilter->AddCuts(esdTrackCutsHG1_tmp); // add once more for tpc only tracks
217 trackFilter->AddCuts(esdTrackCutsH2); // add r_aa cuts
219 trackFilter->AddCuts(esdTrackCutsTPCCOnly); // add QM TPC only track cuts
222 esdfilter->SetGlobalConstrainedFilterMask(1<<8|1<<9); // these tracks are written out as global constrained tracks
223 esdfilter->SetHybridFilterMaskGlobalConstrainedGlobal((1<<4)); // these normal global tracks will be marked as hybrid
225 // esdfilter->SetTPCConstrainedFilterMask(1<<11); // these tracks are written out as tpc constrained tracks
227 esdfilter->SetTrackFilter(trackFilter);
232 Bool_t AddTrackCutsLHC11h(AliAnalysisTaskESDfilter* esdfilter){
235 Printf("%s%d: Creating Track Cuts LHC11h",(char*)__FILE__,__LINE__);
237 // Cuts on primary tracks
238 AliESDtrackCuts* esdTrackCutsL = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
240 // ITS stand-alone tracks
241 AliESDtrackCuts* esdTrackCutsITSsa = new AliESDtrackCuts("ITS stand-alone Track Cuts", "ESD Track Cuts");
242 esdTrackCutsITSsa->SetRequireITSStandAlone(kTRUE);
244 // Pixel OR necessary for the electrons
245 AliESDtrackCuts *itsStrong = new AliESDtrackCuts("ITSorSPD", "pixel requirement for ITS");
246 itsStrong->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
249 // PID for the electrons
250 AliESDpidCuts *electronID = new AliESDpidCuts("Electrons", "Electron PID cuts");
251 electronID->SetTPCnSigmaCut(AliPID::kElectron, 3.);
253 // standard cuts with very loose DCA
254 AliESDtrackCuts* esdTrackCutsH = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE);
255 esdTrackCutsH->SetMaxDCAToVertexXY(2.4);
256 esdTrackCutsH->SetMaxDCAToVertexZ(3.2);
257 esdTrackCutsH->SetDCAToVertex2D(kTRUE);
259 // standard cuts with tight DCA cut
260 AliESDtrackCuts* esdTrackCutsH2 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011();
262 // standard cuts with tight DCA but with requiring the first SDD cluster instead of an SPD cluster
263 // tracks selected by this cut are exclusive to those selected by the previous cut
264 AliESDtrackCuts* esdTrackCutsH3 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011();
265 esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);
266 esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSDD, AliESDtrackCuts::kFirst);
268 // TPC only tracks: Optionally enable the writing of TPConly information
269 // constrained to SPD vertex in the filter below
270 AliESDtrackCuts* esdTrackCutsTPCOnly = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
271 // The following line is needed for 2010 PbPb reprocessing and pp, but not for 2011 PbPb
272 //esdTrackCutsTPCOnly->SetMinNClustersTPC(70);
274 // Extra cuts for hybrids
275 // first the global tracks we want to take
276 AliESDtrackCuts* esdTrackCutsHTG = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE);
277 esdTrackCutsHTG->SetName("Global Hybrid tracks, loose DCA");
278 esdTrackCutsHTG->SetMaxDCAToVertexXY(2.4);
279 esdTrackCutsHTG->SetMaxDCAToVertexZ(3.2);
280 esdTrackCutsHTG->SetDCAToVertex2D(kTRUE);
281 esdTrackCutsHTG->SetMaxChi2TPCConstrainedGlobal(36);
283 // Than the complementary tracks which will be stored as global
284 // constraint, complement is done in the ESDFilter task
285 AliESDtrackCuts* esdTrackCutsHTGC = new AliESDtrackCuts(*esdTrackCutsHTG);
286 esdTrackCutsHTGC->SetName("Global Constraint Hybrid tracks, loose DCA no it requirement");
287 esdTrackCutsHTGC->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kOff);
288 esdTrackCutsHTGC->SetRequireITSRefit(kFALSE);
290 // standard cuts with tight DCA cut, using cluster cut instead of crossed rows (a la 2010 default)
291 AliESDtrackCuts* esdTrackCutsH2Cluster = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kTRUE, 0);
293 // Compose the filter
294 AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
296 trackFilter->AddCuts(esdTrackCutsL);
298 trackFilter->AddCuts(esdTrackCutsITSsa);
300 trackFilter->AddCuts(itsStrong);
301 itsStrong->SetFilterMask(1); // AND with Standard track cuts
303 trackFilter->AddCuts(electronID);
304 electronID->SetFilterMask(4); // AND with Pixel Cuts
306 trackFilter->AddCuts(esdTrackCutsH);
308 trackFilter->AddCuts(esdTrackCutsH2);
310 trackFilter->AddCuts(esdTrackCutsH3);
312 trackFilter->AddCuts(esdTrackCutsTPCOnly);
313 if(enableTPCOnlyAODTracks)esdfilter->SetTPCOnlyFilterMask(128);
314 // 256, 1 << 8 Global Hybrids
315 trackFilter->AddCuts(esdTrackCutsHTG);
316 esdfilter->SetHybridFilterMaskGlobalConstrainedGlobal((1<<8)); // these normal global tracks will be marked as hybrid
317 // 512, 1<< 9 GlobalConstraint Hybrids
318 trackFilter->AddCuts(esdTrackCutsHTGC);
319 esdfilter->SetGlobalConstrainedFilterMask(1<<9); // these tracks are written out as global constrained tracks
320 esdfilter->SetWriteHybridGlobalConstrainedOnly(kTRUE); // write only the complement
322 trackFilter->AddCuts(esdTrackCutsH2Cluster);
323 esdfilter->SetTrackFilter(trackFilter);