]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGPP/ITS/AliAnalysisTaskITSTrackingCheck.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Covering PWGPP to native cmake
[u/mrichter/AliRoot.git] / PWGPP / ITS / AliAnalysisTaskITSTrackingCheck.cxx
CommitLineData
8d63376d 1/**************************************************************************
2 * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/////////////////////////////////////////////////////////////
17//
18// AliAnalysisTask to extract from ESD tracks the information
19// on ITS tracking efficiency and resolutions.
20//
21// Author: A.Dainese, andrea.dainese@pd.infn.it
22/////////////////////////////////////////////////////////////
23
c2893e82 24#include <TStyle.h>
8d63376d 25#include <TChain.h>
26#include <TTree.h>
27#include <TNtuple.h>
28#include <TBranch.h>
29#include <TClonesArray.h>
30#include <TObjArray.h>
31#include <TH1F.h>
32#include <TH2F.h>
33#include <TCanvas.h>
34#include <TParticle.h>
35
36#include "AliAnalysisTask.h"
37#include "AliAnalysisManager.h"
38
39#include "AliMultiplicity.h"
31a96e36 40#include "AliTriggerConfiguration.h"
41#include "AliCDBManager.h"
42#include "AliCDBEntry.h"
8d63376d 43#include "AliVertexerTracks.h"
44#include "AliESDtrack.h"
45#include "AliExternalTrackParam.h"
46#include "AliESDVertex.h"
47#include "AliESDEvent.h"
48#include "AliESDfriend.h"
49#include "AliESDInputHandler.h"
50#include "AliESDInputHandlerRP.h"
f8777463 51#include "AliESDtrackCuts.h"
8d63376d 52#include "AliTrackPointArray.h"
ea1e04c4 53#include "AliITSRecPoint.h"
8d63376d 54
f8777463 55#include "AliInputEventHandler.h"
8d63376d 56#include "AliMCEventHandler.h"
57#include "AliMCEvent.h"
58#include "AliStack.h"
59#include "AliLog.h"
60
61#include "AliGenEventHeader.h"
62#include "AliAnalysisTaskITSTrackingCheck.h"
63
64
65ClassImp(AliAnalysisTaskITSTrackingCheck)
d98a315a 66AliAnalysisTaskITSTrackingCheck::AliAnalysisTaskITSTrackingCheck() :
67AliAnalysisTaskSE(),
b7fb742d 68fNITSLayers(6),
d98a315a 69fReadMC(kFALSE),
70fReadRPLabels(kFALSE),
71fFillNtuples(kFALSE),
72fUseITSSAforNtuples(kFALSE),
0e808834 73fUsePhysSel(kFALSE),
b7fb742d 74fRequireSPDvtx(kTRUE),
9f024d7a 75fRejPileupSPD(kFALSE),
31a96e36 76fCheckSDDIsIn(0),
d98a315a 77fESD(0),
5ac4a99f 78fMinMult(0),
79fMaxMult(1000000),
31a96e36 80fTriggerClass(""),
81fTrigConfig(0),
b7fb742d 82fUsePtBinsForBox(kFALSE),
31a96e36 83fOCDBPath(""),
d98a315a 84fOutput(0),
0e808834 85fHistNEvents(0),
86fHistNEventsFrac(0),
d98a315a 87fHistNtracks(0),
88fHistNclsITSMI(0),
89fHistNclsITSSA(0),
90fHistNclsITSSAInAcc(0),
91fHistClusterMapITSMI(0),
92fHistClusterMapITSMIok(0),
93fHistClusterMapITSMIbad(0),
94fHistClusterMapITSMIskipped(0),
95fHistClusterMapITSMIoutinz(0),
96fHistClusterMapITSMInorefit(0),
97fHistClusterMapITSMInocls(0),
98fHistClusterMapITSMIokoutinzbad(0),
99fHistClusterMapITSSA(0),
100fHistClusterMapITSSAok(0),
101fHistClusterMapITSSAbad(0),
102fHistClusterMapITSSAskipped(0),
103fHistClusterMapITSSAoutinz(0),
104fHistClusterMapITSSAnorefit(0),
105fHistClusterMapITSSAnocls(0),
106fHistClusterMapITSSAokoutinzbad(0),
107fHistClusterMapITSSAInAcc(0),
108fHistClusterMapITSSAokInAcc(0),
109fHistClusterMapITSSAbadInAcc(0),
110fHistClusterMapITSSAskippedInAcc(0),
111fHistClusterMapITSSAoutinzInAcc(0),
112fHistClusterMapITSSAnorefitInAcc(0),
113fHistClusterMapITSSAnoclsInAcc(0),
114fHistClusterMapITSSAokoutinzbadInAcc(0),
115fHistClusterMapModuleITSSAokInAcc(0),
116fHistClusterMapModuleITSSAbadInAcc(0),
117fHistClusterMapModuleITSSAnoclsInAcc(0),
0e808834 118fHistClusterMapModuleITSMIokInAcc(0),
119fHistClusterMapModuleITSMIbadInAcc(0),
120fHistClusterMapModuleITSMInoclsInAcc(0),
0b2e8130 121fHistNClustersMapModule(0),
dbf573fb 122fHistZatSPDouter0ok(0),
123fHistZatSPDouter1ok(0),
124fHistZatSPDouter2ok(0),
125fHistZatSPDouter3ok(0),
126fHistZatSPDouter0notok(0),
127fHistZatSPDouter1notok(0),
128fHistZatSPDouter2notok(0),
129fHistZatSPDouter3notok(0),
0e808834 130fHistxlocSDDok(0),
131fHistzlocSDDok(0),
132fHistxlocVSmodSDDok(0),
133fHistxlocSDDall(0),
134fHistzlocSDDall(0),
d98a315a 135fHistPhiTPCInAcc(0),
547640ff 136fHistEtaTPCInAcc(0),
d98a315a 137fHistPtTPC(0),
138fHistPtTPCInAcc(0),
c94b80d3 139fHistPtTPCInAccTOFbc0(0),
31a96e36 140fHistPtTPCInAccwSDD(0),
141fHistPtTPCInAccTOFbc0wSDD(0),
142fHistPtTPCInAccwoSDD(0),
143fHistPtTPCInAccTOFbc0woSDD(0),
dbf573fb 144fHistPtTPCInAccMCtwoSPD(0),
145fHistPtTPCInAccMConeSPD(0),
0e808834 146fHistdEdxVSPtTPCInAcc(0),
147fHistdEdxVSPtITSTPCsel(0),
0b2e8130 148fHistTPCclsVSPtTPCInAcc(0),
149fHistTPCclsVSPtITSMISPDInAcc(0),
0e808834 150fHistPtVSphiTPCInAcc(0),
151fHistPtTPCInAccNoTRDout(0),
152fHistPtTPCInAccNoTOFout(0),
153fHistPtTPCInAccWithPtTPCAtInnerWall(0),
154fHistPtTPCInAccWithPtTPCAtVtx(0),
155fHistDeltaPtTPC(0),
156fHistPtTPCInAccP(0),
157fHistPtTPCInAccS(0),
158fHistPtTPCInAccPfromStrange(0),
159fHistPtTPCInAccSfromStrange(0),
160fHistPtTPCInAccSfromMat(0),
d98a315a 161fHistPtITSMI2(0),
162fHistPtITSMI3(0),
163fHistPtITSMI4(0),
164fHistPtITSMI5(0),
165fHistPtITSMI6(0),
b7fb742d 166fHistPtITSMI7(0),
d98a315a 167fHistPtITSMISPD(0),
168fHistPtITSMIoneSPD(0),
b7fb742d 169fHistPtITSMItwoSPD(0),
d98a315a 170fHistPtITSMI2InAcc(0),
171fHistPtITSMI3InAcc(0),
172fHistPtITSMI4InAcc(0),
173fHistPtITSMI5InAcc(0),
174fHistPtITSMI6InAcc(0),
b7fb742d 175fHistPtITSMI7InAcc(0),
d98a315a 176fHistPtITSMISPDInAcc(0),
177fHistPtITSMIoneSPDInAcc(0),
b7fb742d 178fHistPtITSMItwoSPDInAcc(0),
c94b80d3 179fHistPtITSMI2InAccTOFbc0(0),
180fHistPtITSMI3InAccTOFbc0(0),
181fHistPtITSMI4InAccTOFbc0(0),
182fHistPtITSMI5InAccTOFbc0(0),
183fHistPtITSMI6InAccTOFbc0(0),
184fHistPtITSMISPDInAccTOFbc0(0),
185fHistPtITSMIoneSPDInAccTOFbc0(0),
31a96e36 186fHistPtITSMI2InAccwSDD(0),
187fHistPtITSMI3InAccwSDD(0),
188fHistPtITSMI4InAccwSDD(0),
189fHistPtITSMI5InAccwSDD(0),
190fHistPtITSMI6InAccwSDD(0),
191fHistPtITSMISPDInAccwSDD(0),
192fHistPtITSMIoneSPDInAccwSDD(0),
193fHistPtITSMI2InAccTOFbc0wSDD(0),
194fHistPtITSMI3InAccTOFbc0wSDD(0),
195fHistPtITSMI4InAccTOFbc0wSDD(0),
196fHistPtITSMI5InAccTOFbc0wSDD(0),
197fHistPtITSMI6InAccTOFbc0wSDD(0),
198fHistPtITSMISPDInAccTOFbc0wSDD(0),
199fHistPtITSMIoneSPDInAccTOFbc0wSDD(0),
200fHistPtITSMI2InAccwoSDD(0),
201fHistPtITSMI3InAccwoSDD(0),
202fHistPtITSMI4InAccwoSDD(0),
203fHistPtITSMI5InAccwoSDD(0),
204fHistPtITSMI6InAccwoSDD(0),
205fHistPtITSMISPDInAccwoSDD(0),
206fHistPtITSMIoneSPDInAccwoSDD(0),
207fHistPtITSMI2InAccTOFbc0woSDD(0),
208fHistPtITSMI3InAccTOFbc0woSDD(0),
209fHistPtITSMI4InAccTOFbc0woSDD(0),
210fHistPtITSMI5InAccTOFbc0woSDD(0),
211fHistPtITSMI6InAccTOFbc0woSDD(0),
212fHistPtITSMISPDInAccTOFbc0woSDD(0),
213fHistPtITSMIoneSPDInAccTOFbc0woSDD(0),
0b2e8130 214fHistPtITSMIoneSPDInAccShared(0),
215fHistPtITSMIoneSPDInAccSharedSPD(0),
216fHistPtITSMISPD1InAccShared(0),
217fHistPtITSMISPD2InAccShared(0),
218fHistPtITSMIoneSPDInAccSharedFake(0),
219fHistPtITSMIoneSPDInAccSharedSPDFake(0),
547640ff 220fHistPhiITSMI2InAcc(0),
221fHistPhiITSMI3InAcc(0),
222fHistPhiITSMI4InAcc(0),
223fHistPhiITSMI5InAcc(0),
224fHistPhiITSMI6InAcc(0),
b7fb742d 225fHistPhiITSMI7InAcc(0),
547640ff 226fHistPhiITSMISPDInAcc(0),
227fHistPhiITSMIoneSPDInAcc(0),
b7fb742d 228fHistPhiITSMItwoSPDInAcc(0),
547640ff 229fHistEtaITSMI2InAcc(0),
230fHistEtaITSMI3InAcc(0),
231fHistEtaITSMI4InAcc(0),
232fHistEtaITSMI5InAcc(0),
233fHistEtaITSMI6InAcc(0),
b7fb742d 234fHistEtaITSMI7InAcc(0),
547640ff 235fHistEtaITSMISPDInAcc(0),
236fHistEtaITSMIoneSPDInAcc(0),
b7fb742d 237fHistEtaITSMItwoSPDInAcc(0),
0faeb7ba 238fHistPtITSMI2InAccFake(0),
239fHistPtITSMI3InAccFake(0),
240fHistPtITSMI4InAccFake(0),
241fHistPtITSMI5InAccFake(0),
242fHistPtITSMI6InAccFake(0),
b7fb742d 243fHistPtITSMI7InAccFake(0),
0faeb7ba 244fHistPtITSMISPDInAccFake(0),
245fHistPtITSMIoneSPDInAccFake(0),
b7fb742d 246fHistPtITSMItwoSPDInAccFake(0),
6c6fab7b 247fHistPtITSMIoneSPDthreeSDDSSDInAcc(0),
0e808834 248fHistPtITSTPCsel(0),
c94b80d3 249fHistPtITSTPCselTOFbc0(0),
31a96e36 250fHistPtITSTPCselwSDD(0),
251fHistPtITSTPCselTOFbc0wSDD(0),
252fHistPtITSTPCselwoSDD(0),
253fHistPtITSTPCselTOFbc0woSDD(0),
0e808834 254fHistPtITSTPCselP(0),
255fHistPtITSTPCselS(0),
af874c53 256fHistPtITSTPCselFake(0),
0e808834 257fHistPtITSTPCselPfromStrange(0),
258fHistPtITSTPCselSfromStrange(0),
259fHistPtITSTPCselSfromMat(0),
260fHistPtITSMI2InAccP(0),
261fHistPtITSMI3InAccP(0),
262fHistPtITSMI4InAccP(0),
263fHistPtITSMI5InAccP(0),
264fHistPtITSMI6InAccP(0),
b7fb742d 265fHistPtITSMI7InAccP(0),
0e808834 266fHistPtITSMISPDInAccP(0),
267fHistPtITSMIoneSPDInAccP(0),
b7fb742d 268fHistPtITSMItwoSPDInAccP(0),
0e808834 269fHistPtITSMI2InAccS(0),
270fHistPtITSMI3InAccS(0),
271fHistPtITSMI4InAccS(0),
272fHistPtITSMI5InAccS(0),
273fHistPtITSMI6InAccS(0),
b7fb742d 274fHistPtITSMI7InAccS(0),
0e808834 275fHistPtITSMISPDInAccS(0),
276fHistPtITSMIoneSPDInAccS(0),
b7fb742d 277fHistPtITSMItwoSPDInAccS(0),
d98a315a 278fHistPtITSMIokbadoutinz6(0),
279fHistPtITSMIokbadoutinz4InAcc(0),
280fHistPtITSMIokbadoutinz5InAcc(0),
281fHistPtITSMIokbadoutinz6InAcc(0),
282fHistPhiITSMIokbadoutinz6InAcc(0),
0e808834 283fHistRProdVtxInAccP(0),
284fHistRProdVtxInAccS(0),
285fHistd0rphiTPCInAccP150200(0),
286fHistd0rphiTPCInAccP500700(0),
287fHistd0rphiTPCInAccP10001500(0),
288fHistd0rphiTPCInAccS150200(0),
289fHistd0rphiTPCInAccS500700(0),
290fHistd0rphiTPCInAccS10001500(0),
291fHistd0rphiITSMISPDInAccP150200(0),
292fHistd0rphiITSMISPDInAccP500700(0),
293fHistd0rphiITSMISPDInAccP10001500(0),
294fHistd0rphiITSMISPDInAccS150200(0),
295fHistd0rphiITSMISPDInAccS500700(0),
296fHistd0rphiITSMISPDInAccS10001500(0),
297fHistd0rphiITSMIoneSPDInAccP150200(0),
dbf573fb 298fHistd0rphiITSMIoneSPDInAccP350450(0),
0e808834 299fHistd0rphiITSMIoneSPDInAccP500700(0),
300fHistd0rphiITSMIoneSPDInAccP10001500(0),
dbf573fb 301fHistd0rphiITSMIoneSPDInAccP25004000(0),
302fHistd0rphiITSMIoneSPDInAccP40008000(0),
0e808834 303fHistd0zITSMIoneSPDInAccP150200(0),
304fHistd0zITSMIoneSPDInAccP500700(0),
305fHistd0zITSMIoneSPDInAccP10001500(0),
9a3cd95a 306fHistd0zVSetaTPCInAccP10001500(0),
0e808834 307fHistd0rphiVSphiITSMIoneSPDInAccP10001500(0),
308fHistd0rphiVSetaITSMIoneSPDInAccP10001500(0),
309fHistd0rphiITSMIoneSPDInAccS150200(0),
dbf573fb 310fHistd0rphiITSMIoneSPDInAccS350450(0),
0e808834 311fHistd0rphiITSMIoneSPDInAccS500700(0),
312fHistd0rphiITSMIoneSPDInAccS500700from22(0),
313fHistd0rphiITSMIoneSPDInAccS500700from211(0),
314fHistd0rphiITSMIoneSPDInAccS500700from310(0),
315fHistd0rphiITSMIoneSPDInAccS500700from321(0),
dbf573fb 316fHistd0rphiITSMIoneSPDInAccS500700from3122(0),
0e808834 317fHistd0rphiITSMIoneSPDInAccS10001500(0),
dbf573fb 318fHistd0rphiITSMIoneSPDInAccS25004000(0),
319fHistd0rphiITSMIoneSPDInAccS40008000(0),
320fHistd0rphiITSMIoneSPDInAccS150200fromStrange(0),
321fHistd0rphiITSMIoneSPDInAccS150200fromMat(0),
322fHistd0rphiITSMIoneSPDInAccS350450fromStrange(0),
323fHistd0rphiITSMIoneSPDInAccS350450fromMat(0),
324fHistd0rphiITSMIoneSPDInAccS500700fromStrange(0),
325fHistd0rphiITSMIoneSPDInAccS500700fromMat(0),
326fHistd0rphiITSMIoneSPDInAccS10001500fromStrange(0),
327fHistd0rphiITSMIoneSPDInAccS10001500fromMat(0),
328fHistd0rphiITSMIoneSPDInAccS25004000fromStrange(0),
329fHistd0rphiITSMIoneSPDInAccS25004000fromMat(0),
330fHistd0rphiITSMIoneSPDInAccS40008000fromStrange(0),
331fHistd0rphiITSMIoneSPDInAccS40008000fromMat(0),
0e808834 332fHistd0zITSMIoneSPDInAccS150200(0),
333fHistd0zITSMIoneSPDInAccS500700(0),
334fHistd0zITSMIoneSPDInAccS10001500(0),
335fHistPDGMoth(0),
336fHistPDGMoth150200(0),
337fHistPDGMoth500700(0),
338fHistPDGMoth10001500(0),
339fHistPDGTrk(0),
af874c53 340fHistITSRedChi2NonFakePt02(0),
341fHistITSRedChi2FakePt02(0),
342fHistITSRedChi2NonFakePt05(0),
343fHistITSRedChi2FakePt05(0),
344fHistITSRedChi2NonFakePt1(0),
345fHistITSRedChi2FakePt1(0),
d98a315a 346fNtupleESDTracks(0),
347fNtupleITSAlignExtra(0),
0e808834 348fNtupleITSAlignSPDTracklets(0),
349fESDtrackCutsTPC(0),
350fESDtrackCutsITSTPC(0)
d98a315a 351{
352 // Constructor
779c3917 353 for(Int_t i=0; i<11; i++) fCountsPerPtBin[i]=0;
0b2e8130 354
d98a315a 355}
8d63376d 356
357//________________________________________________________________________
358AliAnalysisTaskITSTrackingCheck::AliAnalysisTaskITSTrackingCheck(const char *name) :
d98a315a 359AliAnalysisTaskSE(name),
b7fb742d 360fNITSLayers(6),
2308c814 361fReadMC(kFALSE),
8d63376d 362fReadRPLabels(kFALSE),
2308c814 363fFillNtuples(kFALSE),
3f18d0a6 364fUseITSSAforNtuples(kFALSE),
f8777463 365fUsePhysSel(kFALSE),
b7fb742d 366fRequireSPDvtx(kTRUE),
9f024d7a 367fRejPileupSPD(kFALSE),
31a96e36 368fCheckSDDIsIn(0),
8d63376d 369fESD(0),
5ac4a99f 370fMinMult(0),
371fMaxMult(1000000),
31a96e36 372fTriggerClass(""),
373fTrigConfig(0),
b7fb742d 374fUsePtBinsForBox(kFALSE),
31a96e36 375fOCDBPath(""),
8d63376d 376fOutput(0),
f8777463 377fHistNEvents(0),
378fHistNEventsFrac(0),
8d63376d 379fHistNtracks(0),
380fHistNclsITSMI(0),
381fHistNclsITSSA(0),
2308c814 382fHistNclsITSSAInAcc(0),
8d63376d 383fHistClusterMapITSMI(0),
384fHistClusterMapITSMIok(0),
385fHistClusterMapITSMIbad(0),
386fHistClusterMapITSMIskipped(0),
387fHistClusterMapITSMIoutinz(0),
388fHistClusterMapITSMInorefit(0),
389fHistClusterMapITSMInocls(0),
2308c814 390fHistClusterMapITSMIokoutinzbad(0),
8d63376d 391fHistClusterMapITSSA(0),
392fHistClusterMapITSSAok(0),
393fHistClusterMapITSSAbad(0),
394fHistClusterMapITSSAskipped(0),
395fHistClusterMapITSSAoutinz(0),
396fHistClusterMapITSSAnorefit(0),
397fHistClusterMapITSSAnocls(0),
2308c814 398fHistClusterMapITSSAokoutinzbad(0),
399fHistClusterMapITSSAInAcc(0),
400fHistClusterMapITSSAokInAcc(0),
401fHistClusterMapITSSAbadInAcc(0),
402fHistClusterMapITSSAskippedInAcc(0),
403fHistClusterMapITSSAoutinzInAcc(0),
404fHistClusterMapITSSAnorefitInAcc(0),
405fHistClusterMapITSSAnoclsInAcc(0),
406fHistClusterMapITSSAokoutinzbadInAcc(0),
407fHistClusterMapModuleITSSAokInAcc(0),
408fHistClusterMapModuleITSSAbadInAcc(0),
409fHistClusterMapModuleITSSAnoclsInAcc(0),
f8777463 410fHistClusterMapModuleITSMIokInAcc(0),
411fHistClusterMapModuleITSMIbadInAcc(0),
412fHistClusterMapModuleITSMInoclsInAcc(0),
0b2e8130 413fHistNClustersMapModule(0),
dbf573fb 414fHistZatSPDouter0ok(0),
415fHistZatSPDouter1ok(0),
416fHistZatSPDouter2ok(0),
417fHistZatSPDouter3ok(0),
418fHistZatSPDouter0notok(0),
419fHistZatSPDouter1notok(0),
420fHistZatSPDouter2notok(0),
421fHistZatSPDouter3notok(0),
f8777463 422fHistxlocSDDok(0),
423fHistzlocSDDok(0),
424fHistxlocVSmodSDDok(0),
425fHistxlocSDDall(0),
426fHistzlocSDDall(0),
3f18d0a6 427fHistPhiTPCInAcc(0),
547640ff 428fHistEtaTPCInAcc(0),
8d63376d 429fHistPtTPC(0),
3f18d0a6 430fHistPtTPCInAcc(0),
c94b80d3 431fHistPtTPCInAccTOFbc0(0),
31a96e36 432fHistPtTPCInAccwSDD(0),
433fHistPtTPCInAccTOFbc0wSDD(0),
434fHistPtTPCInAccwoSDD(0),
435fHistPtTPCInAccTOFbc0woSDD(0),
dbf573fb 436fHistPtTPCInAccMCtwoSPD(0),
437fHistPtTPCInAccMConeSPD(0),
f8777463 438fHistdEdxVSPtTPCInAcc(0),
439fHistdEdxVSPtITSTPCsel(0),
0b2e8130 440fHistTPCclsVSPtTPCInAcc(0),
441fHistTPCclsVSPtITSMISPDInAcc(0),
f8777463 442fHistPtVSphiTPCInAcc(0),
443fHistPtTPCInAccNoTRDout(0),
444fHistPtTPCInAccNoTOFout(0),
445fHistPtTPCInAccWithPtTPCAtInnerWall(0),
446fHistPtTPCInAccWithPtTPCAtVtx(0),
447fHistDeltaPtTPC(0),
448fHistPtTPCInAccP(0),
449fHistPtTPCInAccS(0),
450fHistPtTPCInAccPfromStrange(0),
451fHistPtTPCInAccSfromStrange(0),
452fHistPtTPCInAccSfromMat(0),
8d63376d 453fHistPtITSMI2(0),
454fHistPtITSMI3(0),
455fHistPtITSMI4(0),
456fHistPtITSMI5(0),
457fHistPtITSMI6(0),
b7fb742d 458fHistPtITSMI7(0),
8d63376d 459fHistPtITSMISPD(0),
7b4d74b2 460fHistPtITSMIoneSPD(0),
b7fb742d 461fHistPtITSMItwoSPD(0),
3f18d0a6 462fHistPtITSMI2InAcc(0),
463fHistPtITSMI3InAcc(0),
464fHistPtITSMI4InAcc(0),
465fHistPtITSMI5InAcc(0),
466fHistPtITSMI6InAcc(0),
b7fb742d 467fHistPtITSMI7InAcc(0),
3f18d0a6 468fHistPtITSMISPDInAcc(0),
7b4d74b2 469fHistPtITSMIoneSPDInAcc(0),
b7fb742d 470fHistPtITSMItwoSPDInAcc(0),
c94b80d3 471fHistPtITSMI2InAccTOFbc0(0),
472fHistPtITSMI3InAccTOFbc0(0),
473fHistPtITSMI4InAccTOFbc0(0),
474fHistPtITSMI5InAccTOFbc0(0),
475fHistPtITSMI6InAccTOFbc0(0),
476fHistPtITSMISPDInAccTOFbc0(0),
477fHistPtITSMIoneSPDInAccTOFbc0(0),
31a96e36 478fHistPtITSMI2InAccwSDD(0),
479fHistPtITSMI3InAccwSDD(0),
480fHistPtITSMI4InAccwSDD(0),
481fHistPtITSMI5InAccwSDD(0),
482fHistPtITSMI6InAccwSDD(0),
483fHistPtITSMISPDInAccwSDD(0),
484fHistPtITSMIoneSPDInAccwSDD(0),
485fHistPtITSMI2InAccTOFbc0wSDD(0),
486fHistPtITSMI3InAccTOFbc0wSDD(0),
487fHistPtITSMI4InAccTOFbc0wSDD(0),
488fHistPtITSMI5InAccTOFbc0wSDD(0),
489fHistPtITSMI6InAccTOFbc0wSDD(0),
490fHistPtITSMISPDInAccTOFbc0wSDD(0),
491fHistPtITSMIoneSPDInAccTOFbc0wSDD(0),
492fHistPtITSMI2InAccwoSDD(0),
493fHistPtITSMI3InAccwoSDD(0),
494fHistPtITSMI4InAccwoSDD(0),
495fHistPtITSMI5InAccwoSDD(0),
496fHistPtITSMI6InAccwoSDD(0),
497fHistPtITSMISPDInAccwoSDD(0),
498fHistPtITSMIoneSPDInAccwoSDD(0),
499fHistPtITSMI2InAccTOFbc0woSDD(0),
500fHistPtITSMI3InAccTOFbc0woSDD(0),
501fHistPtITSMI4InAccTOFbc0woSDD(0),
502fHistPtITSMI5InAccTOFbc0woSDD(0),
503fHistPtITSMI6InAccTOFbc0woSDD(0),
504fHistPtITSMISPDInAccTOFbc0woSDD(0),
505fHistPtITSMIoneSPDInAccTOFbc0woSDD(0),
0b2e8130 506fHistPtITSMIoneSPDInAccShared(0),
507fHistPtITSMIoneSPDInAccSharedSPD(0),
508fHistPtITSMISPD1InAccShared(0),
509fHistPtITSMISPD2InAccShared(0),
510fHistPtITSMIoneSPDInAccSharedFake(0),
511fHistPtITSMIoneSPDInAccSharedSPDFake(0),
547640ff 512fHistPhiITSMI2InAcc(0),
513fHistPhiITSMI3InAcc(0),
514fHistPhiITSMI4InAcc(0),
515fHistPhiITSMI5InAcc(0),
516fHistPhiITSMI6InAcc(0),
b7fb742d 517fHistPhiITSMI7InAcc(0),
547640ff 518fHistPhiITSMISPDInAcc(0),
519fHistPhiITSMIoneSPDInAcc(0),
b7fb742d 520fHistPhiITSMItwoSPDInAcc(0),
547640ff 521fHistEtaITSMI2InAcc(0),
522fHistEtaITSMI3InAcc(0),
523fHistEtaITSMI4InAcc(0),
524fHistEtaITSMI5InAcc(0),
525fHistEtaITSMI6InAcc(0),
b7fb742d 526fHistEtaITSMI7InAcc(0),
547640ff 527fHistEtaITSMISPDInAcc(0),
528fHistEtaITSMIoneSPDInAcc(0),
b7fb742d 529fHistEtaITSMItwoSPDInAcc(0),
0faeb7ba 530fHistPtITSMI2InAccFake(0),
531fHistPtITSMI3InAccFake(0),
532fHistPtITSMI4InAccFake(0),
533fHistPtITSMI5InAccFake(0),
534fHistPtITSMI6InAccFake(0),
b7fb742d 535fHistPtITSMI7InAccFake(0),
0faeb7ba 536fHistPtITSMISPDInAccFake(0),
537fHistPtITSMIoneSPDInAccFake(0),
b7fb742d 538fHistPtITSMItwoSPDInAccFake(0),
6c6fab7b 539fHistPtITSMIoneSPDthreeSDDSSDInAcc(0),
f8777463 540fHistPtITSTPCsel(0),
c94b80d3 541fHistPtITSTPCselTOFbc0(0),
31a96e36 542fHistPtITSTPCselwSDD(0),
543fHistPtITSTPCselTOFbc0wSDD(0),
544fHistPtITSTPCselwoSDD(0),
545fHistPtITSTPCselTOFbc0woSDD(0),
f8777463 546fHistPtITSTPCselP(0),
547fHistPtITSTPCselS(0),
af874c53 548fHistPtITSTPCselFake(0),
f8777463 549fHistPtITSTPCselPfromStrange(0),
550fHistPtITSTPCselSfromStrange(0),
551fHistPtITSTPCselSfromMat(0),
552fHistPtITSMI2InAccP(0),
553fHistPtITSMI3InAccP(0),
554fHistPtITSMI4InAccP(0),
555fHistPtITSMI5InAccP(0),
556fHistPtITSMI6InAccP(0),
b7fb742d 557fHistPtITSMI7InAccP(0),
f8777463 558fHistPtITSMISPDInAccP(0),
559fHistPtITSMIoneSPDInAccP(0),
b7fb742d 560fHistPtITSMItwoSPDInAccP(0),
f8777463 561fHistPtITSMI2InAccS(0),
562fHistPtITSMI3InAccS(0),
563fHistPtITSMI4InAccS(0),
564fHistPtITSMI5InAccS(0),
565fHistPtITSMI6InAccS(0),
b7fb742d 566fHistPtITSMI7InAccS(0),
f8777463 567fHistPtITSMISPDInAccS(0),
568fHistPtITSMIoneSPDInAccS(0),
b7fb742d 569fHistPtITSMItwoSPDInAccS(0),
3f18d0a6 570fHistPtITSMIokbadoutinz6(0),
c2893e82 571fHistPtITSMIokbadoutinz4InAcc(0),
572fHistPtITSMIokbadoutinz5InAcc(0),
3f18d0a6 573fHistPtITSMIokbadoutinz6InAcc(0),
574fHistPhiITSMIokbadoutinz6InAcc(0),
f8777463 575fHistRProdVtxInAccP(0),
576fHistRProdVtxInAccS(0),
577fHistd0rphiTPCInAccP150200(0),
578fHistd0rphiTPCInAccP500700(0),
579fHistd0rphiTPCInAccP10001500(0),
580fHistd0rphiTPCInAccS150200(0),
581fHistd0rphiTPCInAccS500700(0),
582fHistd0rphiTPCInAccS10001500(0),
583fHistd0rphiITSMISPDInAccP150200(0),
584fHistd0rphiITSMISPDInAccP500700(0),
585fHistd0rphiITSMISPDInAccP10001500(0),
586fHistd0rphiITSMISPDInAccS150200(0),
587fHistd0rphiITSMISPDInAccS500700(0),
588fHistd0rphiITSMISPDInAccS10001500(0),
589fHistd0rphiITSMIoneSPDInAccP150200(0),
dbf573fb 590fHistd0rphiITSMIoneSPDInAccP350450(0),
f8777463 591fHistd0rphiITSMIoneSPDInAccP500700(0),
592fHistd0rphiITSMIoneSPDInAccP10001500(0),
dbf573fb 593fHistd0rphiITSMIoneSPDInAccP25004000(0),
594fHistd0rphiITSMIoneSPDInAccP40008000(0),
f8777463 595fHistd0zITSMIoneSPDInAccP150200(0),
596fHistd0zITSMIoneSPDInAccP500700(0),
597fHistd0zITSMIoneSPDInAccP10001500(0),
9a3cd95a 598fHistd0zVSetaTPCInAccP10001500(0),
f8777463 599fHistd0rphiVSphiITSMIoneSPDInAccP10001500(0),
600fHistd0rphiVSetaITSMIoneSPDInAccP10001500(0),
601fHistd0rphiITSMIoneSPDInAccS150200(0),
dbf573fb 602fHistd0rphiITSMIoneSPDInAccS350450(0),
f8777463 603fHistd0rphiITSMIoneSPDInAccS500700(0),
604fHistd0rphiITSMIoneSPDInAccS500700from22(0),
605fHistd0rphiITSMIoneSPDInAccS500700from211(0),
606fHistd0rphiITSMIoneSPDInAccS500700from310(0),
607fHistd0rphiITSMIoneSPDInAccS500700from321(0),
dbf573fb 608fHistd0rphiITSMIoneSPDInAccS500700from3122(0),
f8777463 609fHistd0rphiITSMIoneSPDInAccS10001500(0),
dbf573fb 610fHistd0rphiITSMIoneSPDInAccS25004000(0),
611fHistd0rphiITSMIoneSPDInAccS40008000(0),
612fHistd0rphiITSMIoneSPDInAccS150200fromStrange(0),
613fHistd0rphiITSMIoneSPDInAccS150200fromMat(0),
614fHistd0rphiITSMIoneSPDInAccS350450fromStrange(0),
615fHistd0rphiITSMIoneSPDInAccS350450fromMat(0),
616fHistd0rphiITSMIoneSPDInAccS500700fromStrange(0),
617fHistd0rphiITSMIoneSPDInAccS500700fromMat(0),
618fHistd0rphiITSMIoneSPDInAccS10001500fromStrange(0),
619fHistd0rphiITSMIoneSPDInAccS10001500fromMat(0),
620fHistd0rphiITSMIoneSPDInAccS25004000fromStrange(0),
621fHistd0rphiITSMIoneSPDInAccS25004000fromMat(0),
622fHistd0rphiITSMIoneSPDInAccS40008000fromStrange(0),
623fHistd0rphiITSMIoneSPDInAccS40008000fromMat(0),
f8777463 624fHistd0zITSMIoneSPDInAccS150200(0),
625fHistd0zITSMIoneSPDInAccS500700(0),
626fHistd0zITSMIoneSPDInAccS10001500(0),
627fHistPDGMoth(0),
628fHistPDGMoth150200(0),
629fHistPDGMoth500700(0),
630fHistPDGMoth10001500(0),
631fHistPDGTrk(0),
af874c53 632fHistITSRedChi2NonFakePt02(0),
633fHistITSRedChi2FakePt02(0),
634fHistITSRedChi2NonFakePt05(0),
635fHistITSRedChi2FakePt05(0),
636fHistITSRedChi2NonFakePt1(0),
637fHistITSRedChi2FakePt1(0),
8d63376d 638fNtupleESDTracks(0),
639fNtupleITSAlignExtra(0),
f8777463 640fNtupleITSAlignSPDTracklets(0),
641fESDtrackCutsTPC(0),
642fESDtrackCutsITSTPC(0)
8d63376d 643{
644 // Constructor
645
f8777463 646 for(Int_t i=0; i<11; i++) fCountsPerPtBin[i]=0;
8d63376d 647
648 // Define input and output slots here
8d63376d 649 // Output slot #0 writes into a TList container
d98a315a 650 DefineOutput(1, TList::Class()); //My private output
8d63376d 651}
652//________________________________________________________________________
653AliAnalysisTaskITSTrackingCheck::~AliAnalysisTaskITSTrackingCheck()
654{
655 // Destructor
656
657 // histograms are in the output list and deleted when the output
658 // list is deleted by the TSelector dtor
659
9c12af5f 660 if (fOutput && !AliAnalysisManager::GetAnalysisManager()->IsProofMode()) {
8d63376d 661 delete fOutput;
662 fOutput = 0;
663 }
664}
8d63376d 665//________________________________________________________________________
d98a315a 666void AliAnalysisTaskITSTrackingCheck::UserCreateOutputObjects()
8d63376d 667{
668 // Create histograms
669 // Called once
670
c2893e82 671 gStyle->SetHistLineWidth(2);
672
f8777463 673 Int_t nPtBins=34;
b7fb742d 674 Float_t xPtBinsDefault[35]={0,0.025,0.05,0.075,0.1,0.125,0.15,0.175,0.2,0.225,0.25,0.275,0.3,0.325,0.35,0.375,0.4,0.5,0.6,0.7,0.8,1.0,1.5,2.,2.5,3,4,5,6,8,10,15,20,25,30};
675 Float_t xPtBinsForBox[35]={0,0.025,0.05,0.075,0.08,0.125,0.15,0.175,0.2,0.225,0.25,0.275,0.280,0.325,0.35,0.375,0.4,0.45,0.55,0.7,0.75,0.9,1.2,1.7,2.5,3.5,4.5,5,6,8,13,15,23,28,35};
676 Float_t xPtBins[35];
677 for(Int_t ib=0; ib<=nPtBins; ib++) {
678 if(fUsePtBinsForBox) {
679 xPtBins[ib]=xPtBinsForBox[ib];
680 } else {
681 xPtBins[ib]=xPtBinsDefault[ib];
682 }
683 }
684
3f18d0a6 685
f8777463 686 for(Int_t i=0; i<11; i++) fCountsPerPtBin[i]=0;
8d63376d 687
688 // Several histograms are more conveniently managed in a TList
689 fOutput = new TList;
690 fOutput->SetOwner();
691
f8777463 692 fHistPDGMoth = new TH1F("fHistPDGMoth","",4000,0,4000);
693 fOutput->Add(fHistPDGMoth);
694 fHistPDGMoth150200 = new TH1F("fHistPDGMoth150200","",4000,0,4000);
695 fOutput->Add(fHistPDGMoth150200);
696 fHistPDGMoth500700 = new TH1F("fHistPDGMoth500700","",4000,0,4000);
697 fOutput->Add(fHistPDGMoth500700);
698 fHistPDGMoth10001500 = new TH1F("fHistPDGMoth10001500","",4000,0,4000);
699 fOutput->Add(fHistPDGMoth10001500);
700 fHistPDGTrk = new TH1F("fHistPDGTrk","",4000,0,4000);
701 fOutput->Add(fHistPDGTrk);
702
703 fHistNEvents = new TH1F("fHistNEvents", "Events: -1 tot, 0 sel, 1 vSPD3D, 2 vSPDZ, 3 vSPD, 4 vTRK; Type; N Events",12, -1.5, 10.5);
704 fHistNEvents->SetMinimum(0);
705 fOutput->Add(fHistNEvents);
706
707 fHistNEventsFrac = (TH1F*)fHistNEvents->Clone("fHistNEventsFrac");
708 fOutput->Add(fHistNEventsFrac);
709
8d63376d 710 fHistNtracks = new TH1F("fHistNtracks", "N ESD tracks; N tracks; Events",5000, -0.5, 4999.5);
711 fHistNtracks->Sumw2();
712 fHistNtracks->SetMinimum(0);
713 fOutput->Add(fHistNtracks);
714
b7fb742d 715 fHistNclsITSMI = new TH1F("fHistNclsITSMI", "N ITS clusters per track (MI); N clusters; Counts",8, -0.5, 7.5);
8d63376d 716 fHistNclsITSMI->Sumw2();
717 fHistNclsITSMI->SetMinimum(0);
718 fOutput->Add(fHistNclsITSMI);
719
2308c814 720 fHistNclsITSSAInAcc = new TH1F("fHistNclsITSSAInAcc", "N ITS clusters per track (SA); N clusters; Counts",7, -0.5, 6.5);
721 fHistNclsITSSAInAcc->Sumw2();
722 fHistNclsITSSAInAcc->SetMinimum(0);
723 fOutput->Add(fHistNclsITSSAInAcc);
724
b7fb742d 725 fHistNclsITSSA = new TH1F("fHistNclsITSSA", "N ITS clusters per track (SA); N clusters; Counts",8, -0.5, 7.5);
8d63376d 726 fHistNclsITSSA->Sumw2();
727 fHistNclsITSSA->SetMinimum(0);
728 fOutput->Add(fHistNclsITSSA);
729
b7fb742d 730 fHistClusterMapITSMI = new TH1F("fHistClusterMapITSMI", "N tracks with point on Layer (MI); Layer; N tracks",7, -0.5, 6.5);
8d63376d 731 fHistClusterMapITSMI->Sumw2();
732 fHistClusterMapITSMI->SetMinimum(0);
733 fOutput->Add(fHistClusterMapITSMI);
734
b7fb742d 735 fHistClusterMapITSSA = new TH1F("fHistClusterMapITSSA", "N tracks with point on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
8d63376d 736 fHistClusterMapITSSA->Sumw2();
737 fHistClusterMapITSSA->SetMinimum(0);
738 fOutput->Add(fHistClusterMapITSSA);
739
b7fb742d 740 fHistClusterMapITSSAInAcc = new TH1F("fHistClusterMapITSSAInAcc", "N tracks with point on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 741 fHistClusterMapITSSAInAcc->Sumw2();
742 fHistClusterMapITSSAInAcc->SetMinimum(0);
743 fOutput->Add(fHistClusterMapITSSAInAcc);
744
b7fb742d 745 fHistClusterMapITSMIok = new TH1F("fHistClusterMapITSMIok", "N tracks with ok on Layer (MI); Layer; N tracks",7, -0.5, 6.5);
8d63376d 746 fHistClusterMapITSMIok->Sumw2();
747 fHistClusterMapITSMIok->SetMinimum(0);
748 fOutput->Add(fHistClusterMapITSMIok);
749
0b2e8130 750
751
b7fb742d 752 fHistClusterMapITSSAokInAcc = new TH1F("fHistClusterMapITSSAokInAcc", "N tracks with ok on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 753 fHistClusterMapITSSAokInAcc->Sumw2();
754 fHistClusterMapITSSAokInAcc->SetMinimum(0);
755 fOutput->Add(fHistClusterMapITSSAokInAcc);
756
757 fHistClusterMapModuleITSSAokInAcc = new TH1F("fHistClusterMapModuleITSSAokInAcc", "N tracks with ok on Module (SA); Module; N tracks",2198, -0.5, 2197.5);
758 fHistClusterMapModuleITSSAokInAcc->SetMinimum(0);
759 fOutput->Add(fHistClusterMapModuleITSSAokInAcc);
760
f8777463 761 fHistClusterMapModuleITSMIokInAcc = new TH1F("fHistClusterMapModuleITSMIokInAcc", "N tracks with ok on Module (MI); Module; N tracks",2198, -0.5, 2197.5);
762 fHistClusterMapModuleITSMIokInAcc->SetMinimum(0);
763 fOutput->Add(fHistClusterMapModuleITSMIokInAcc);
764
0b2e8130 765 fHistNClustersMapModule = new TH1F("fHistNClustersMapModule", "N clusters; Module; N tracks",2198, -0.5, 2197.5);
766 fHistNClustersMapModule->SetMinimum(0);
767 fOutput->Add(fHistNClustersMapModule);
768
b7fb742d 769 fHistClusterMapITSSAok = new TH1F("fHistClusterMapITSSAok", "N tracks with ok on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
8d63376d 770 fHistClusterMapITSSAok->Sumw2();
771 fHistClusterMapITSSAok->SetMinimum(0);
772 fOutput->Add(fHistClusterMapITSSAok);
773
b7fb742d 774 fHistClusterMapITSMIbad = new TH1F("fHistClusterMapITSMIbad", "N tracks with bad on Layer (MI); Layer; N tracks",7, -0.5, 6.5);
8d63376d 775 fHistClusterMapITSMIbad->Sumw2();
776 fHistClusterMapITSMIbad->SetMinimum(0);
777 fOutput->Add(fHistClusterMapITSMIbad);
778
b7fb742d 779 fHistClusterMapITSSAbadInAcc = new TH1F("fHistClusterMapITSSAbadInAcc", "N tracks with bad on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 780 fHistClusterMapITSSAbadInAcc->Sumw2();
781 fHistClusterMapITSSAbadInAcc->SetMinimum(0);
782 fOutput->Add(fHistClusterMapITSSAbadInAcc);
783
784 fHistClusterMapModuleITSSAbadInAcc = new TH1F("fHistClusterMapModuleITSSAbadInAcc", "N tracks with bad on Module (SA); Module; N tracks",2198, -0.5, 2197.5);
785 fHistClusterMapModuleITSSAbadInAcc->SetMinimum(0);
786 fOutput->Add(fHistClusterMapModuleITSSAbadInAcc);
787
f8777463 788 fHistClusterMapModuleITSMIbadInAcc = new TH1F("fHistClusterMapModuleITSMIbadInAcc", "N tracks with bad on Module (MI); Module; N tracks",2198, -0.5, 2197.5);
789 fHistClusterMapModuleITSMIbadInAcc->SetMinimum(0);
790 fOutput->Add(fHistClusterMapModuleITSMIbadInAcc);
791
b7fb742d 792 fHistClusterMapITSSAbad = new TH1F("fHistClusterMapITSSAbad", "N tracks with bad on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
8d63376d 793 fHistClusterMapITSSAbad->Sumw2();
794 fHistClusterMapITSSAbad->SetMinimum(0);
795 fOutput->Add(fHistClusterMapITSSAbad);
796
b7fb742d 797 fHistClusterMapITSMIskipped = new TH1F("fHistClusterMapITSMIskipped", "N tracks with skip on Layer (MI); Layer; N tracks",7, -0.5, 6.5);
8d63376d 798 fHistClusterMapITSMIskipped->Sumw2();
799 fHistClusterMapITSMIskipped->SetMinimum(0);
800 fOutput->Add(fHistClusterMapITSMIskipped);
801
b7fb742d 802 fHistClusterMapITSSAskippedInAcc = new TH1F("fHistClusterMapITSSAskippedInAcc", "N tracks with skip on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 803 fHistClusterMapITSSAskippedInAcc->Sumw2();
804 fHistClusterMapITSSAskippedInAcc->SetMinimum(0);
805 fOutput->Add(fHistClusterMapITSSAskippedInAcc);
806
b7fb742d 807 fHistClusterMapITSSAskipped = new TH1F("fHistClusterMapITSSAskipped", "N tracks with skip on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
8d63376d 808 fHistClusterMapITSSAskipped->Sumw2();
809 fHistClusterMapITSSAskipped->SetMinimum(0);
810 fOutput->Add(fHistClusterMapITSSAskipped);
811
b7fb742d 812 fHistClusterMapITSMIoutinz = new TH1F("fHistClusterMapITSMIoutinz", "N tracks out in z on Layer (MI); Layer; N tracks",7, -0.5, 6.5);
8d63376d 813 fHistClusterMapITSMIoutinz->Sumw2();
814 fHistClusterMapITSMIoutinz->SetMinimum(0);
815 fOutput->Add(fHistClusterMapITSMIoutinz);
816
b7fb742d 817 fHistClusterMapITSSAoutinzInAcc = new TH1F("fHistClusterMapITSSAoutinzInAcc", "N tracks with out in z on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 818 fHistClusterMapITSSAoutinzInAcc->Sumw2();
819 fHistClusterMapITSSAoutinzInAcc->SetMinimum(0);
820 fOutput->Add(fHistClusterMapITSSAoutinzInAcc);
821
b7fb742d 822 fHistClusterMapITSSAoutinz = new TH1F("fHistClusterMapITSSAoutinz", "N tracks with out in z on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
8d63376d 823 fHistClusterMapITSSAoutinz->Sumw2();
824 fHistClusterMapITSSAoutinz->SetMinimum(0);
825 fOutput->Add(fHistClusterMapITSSAoutinz);
826
b7fb742d 827 fHistClusterMapITSSAokoutinzbad = new TH1F("fHistClusterMapITSSAokoutinzbad", "N tracks with cluster or bad zone or out in z (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 828 fHistClusterMapITSSAokoutinzbad->Sumw2();
829 fHistClusterMapITSSAokoutinzbad->SetMinimum(0);
830 fOutput->Add(fHistClusterMapITSSAokoutinzbad);
831
b7fb742d 832 fHistClusterMapITSMIokoutinzbad = new TH1F("fHistClusterMapITSMIokoutinzbad", "N tracks with cluster or bad zone or out in z (MI); Layer; N tracks",7, -0.5, 6.5);
2308c814 833 fHistClusterMapITSMIokoutinzbad->Sumw2();
834 fHistClusterMapITSMIokoutinzbad->SetMinimum(0);
835 fOutput->Add(fHistClusterMapITSMIokoutinzbad);
836
b7fb742d 837 fHistClusterMapITSSAokoutinzbadInAcc = new TH1F("fHistClusterMapITSSAokoutinzbadInAcc", "N tracks with cluster or bad zone or out in z (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 838 fHistClusterMapITSSAokoutinzbadInAcc->Sumw2();
839 fHistClusterMapITSSAokoutinzbadInAcc->SetMinimum(0);
840 fOutput->Add(fHistClusterMapITSSAokoutinzbadInAcc);
841
b7fb742d 842 fHistClusterMapITSMInorefit = new TH1F("fHistClusterMapITSMInorefit", "N tracks with norefit on Layer (MI); Layer; N tracks",7, -0.5, 6.5);
8d63376d 843 fHistClusterMapITSMInorefit->Sumw2();
844 fHistClusterMapITSMInorefit->SetMinimum(0);
845 fOutput->Add(fHistClusterMapITSMInorefit);
846
b7fb742d 847 fHistClusterMapITSSAnorefitInAcc = new TH1F("fHistClusterMapITSSAnorefitInAcc", "N tracks with norefit on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 848 fHistClusterMapITSSAnorefitInAcc->Sumw2();
849 fHistClusterMapITSSAnorefitInAcc->SetMinimum(0);
850 fOutput->Add(fHistClusterMapITSSAnorefitInAcc);
851
b7fb742d 852 fHistClusterMapITSSAnorefit = new TH1F("fHistClusterMapITSSAnorefit", "N tracks with norefit on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
8d63376d 853 fHistClusterMapITSSAnorefit->Sumw2();
854 fHistClusterMapITSSAnorefit->SetMinimum(0);
855 fOutput->Add(fHistClusterMapITSSAnorefit);
856
b7fb742d 857 fHistClusterMapITSMInocls = new TH1F("fHistClusterMapITSMInocls", "N tracks with nocls on Layer (MI); Layer; N tracks",7, -0.5, 6.5);
8d63376d 858 fHistClusterMapITSMInocls->Sumw2();
859 fHistClusterMapITSMInocls->SetMinimum(0);
860 fOutput->Add(fHistClusterMapITSMInocls);
861
b7fb742d 862 fHistClusterMapITSSAnoclsInAcc = new TH1F("fHistClusterMapITSSAnoclsInAcc", "N tracks with nocls on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
2308c814 863 fHistClusterMapITSSAnoclsInAcc->Sumw2();
864 fHistClusterMapITSSAnoclsInAcc->SetMinimum(0);
865 fOutput->Add(fHistClusterMapITSSAnoclsInAcc);
866
867 fHistClusterMapModuleITSSAnoclsInAcc = new TH1F("fHistClusterMapModuleITSSAnoclsInAcc", "N tracks with nocls on Module (SA); Module; N tracks",2198, -0.5, 2197.5);
868 fHistClusterMapModuleITSSAnoclsInAcc->SetMinimum(0);
869 fOutput->Add(fHistClusterMapModuleITSSAnoclsInAcc);
870
f8777463 871 fHistClusterMapModuleITSMInoclsInAcc = new TH1F("fHistClusterMapModuleITSMInoclsInAcc", "N tracks with nocls on Module (MI); Module; N tracks",2198, -0.5, 2197.5);
872 fHistClusterMapModuleITSMInoclsInAcc->SetMinimum(0);
873 fOutput->Add(fHistClusterMapModuleITSMInoclsInAcc);
874
b7fb742d 875 fHistClusterMapITSSAnocls = new TH1F("fHistClusterMapITSSAnocls", "N tracks with nocls on Layer (SA); Layer; N tracks",7, -0.5, 6.5);
8d63376d 876 fHistClusterMapITSSAnocls->Sumw2();
877 fHistClusterMapITSSAnocls->SetMinimum(0);
878 fOutput->Add(fHistClusterMapITSSAnocls);
879
dbf573fb 880
881 fHistZatSPDouter0ok = new TH1F("fHistZatSPDouter0ok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
882 fHistZatSPDouter0ok->Sumw2();
883 fHistZatSPDouter0ok->SetMinimum(0);
884 fOutput->Add(fHistZatSPDouter0ok);
885 fHistZatSPDouter1ok = new TH1F("fHistZatSPDouter1ok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
886 fHistZatSPDouter1ok->Sumw2();
887 fHistZatSPDouter1ok->SetMinimum(0);
888 fOutput->Add(fHistZatSPDouter1ok);
889 fHistZatSPDouter2ok = new TH1F("fHistZatSPDouter2ok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
890 fHistZatSPDouter2ok->Sumw2();
891 fHistZatSPDouter2ok->SetMinimum(0);
892 fOutput->Add(fHistZatSPDouter2ok);
893 fHistZatSPDouter3ok = new TH1F("fHistZatSPDouter3ok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
894 fHistZatSPDouter3ok->Sumw2();
895 fHistZatSPDouter3ok->SetMinimum(0);
896 fOutput->Add(fHistZatSPDouter3ok);
897 fHistZatSPDouter0notok = new TH1F("fHistZatSPDouter0notok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
898 fHistZatSPDouter0notok->Sumw2();
899 fHistZatSPDouter0notok->SetMinimum(0);
900 fOutput->Add(fHistZatSPDouter0notok);
901 fHistZatSPDouter1notok = new TH1F("fHistZatSPDouter1notok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
902 fHistZatSPDouter1notok->Sumw2();
903 fHistZatSPDouter1notok->SetMinimum(0);
904 fOutput->Add(fHistZatSPDouter1notok);
905 fHistZatSPDouter2notok = new TH1F("fHistZatSPDouter2notok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
906 fHistZatSPDouter2notok->Sumw2();
907 fHistZatSPDouter2notok->SetMinimum(0);
908 fOutput->Add(fHistZatSPDouter2notok);
909 fHistZatSPDouter3notok = new TH1F("fHistZatSPDouter3notok", "Track z at SPDouter radius; z [cm]; N tracks",3200,-16,16);
910 fHistZatSPDouter3notok->Sumw2();
911 fHistZatSPDouter3notok->SetMinimum(0);
912 fOutput->Add(fHistZatSPDouter3notok);
913
f8777463 914 fHistxlocSDDok = new TH1F("fHistxlocSDDok", "SDD points; xloc [cm]; N tracks",75, -3.75, 3.75);
915 fHistxlocSDDok->Sumw2();
916 fHistxlocSDDok->SetMinimum(0);
917 fOutput->Add(fHistxlocSDDok);
918
919 fHistxlocVSmodSDDok = new TH2F("fHistxlocVSmodSDDok", "SDD points; module; xloc [cm]",260,239.5,499.5,25, -3.75, 3.75);
920 fOutput->Add(fHistxlocVSmodSDDok);
921
922 fHistzlocSDDok = new TH1F("fHistzlocSDDok", "SDD points; zloc [cm]; N tracks",77, -3.85, 3.85);
923 fHistzlocSDDok->Sumw2();
924 fHistzlocSDDok->SetMinimum(0);
925 fOutput->Add(fHistzlocSDDok);
926
927 fHistxlocSDDall = new TH1F("fHistxlocSDDall", "SDD points; xloc [cm]; N tracks",75, -3.75, 3.75);
928 fHistxlocSDDall->Sumw2();
929 fHistxlocSDDall->SetMinimum(0);
930 fOutput->Add(fHistxlocSDDall);
931
932 fHistzlocSDDall = new TH1F("fHistzlocSDDall", "SDD points; zloc [cm]; N tracks",77, -3.85, 3.85);
933 fHistzlocSDDall->Sumw2();
934 fHistzlocSDDall->SetMinimum(0);
935 fOutput->Add(fHistzlocSDDall);
0b2e8130 936
f8777463 937
3f18d0a6 938 fHistPhiTPCInAcc = new TH1F("fHistPhiTPCInAcc","Azimuthal distribution of TPC tracks; #phi; N tracks",100, 0, 2.*3.1415);
939 fHistPhiTPCInAcc->Sumw2();
940 fHistPhiTPCInAcc->SetMinimum(0);
941 fOutput->Add(fHistPhiTPCInAcc);
942
547640ff 943 fHistEtaTPCInAcc = new TH1F("fHistEtaTPCInAcc","Eta distribution of TPC tracks; #eta; N tracks",100, -1.5, +1.5);
944 fHistEtaTPCInAcc->Sumw2();
945 fHistEtaTPCInAcc->SetMinimum(0);
946 fOutput->Add(fHistEtaTPCInAcc);
947
0b2e8130 948
3f18d0a6 949 fHistPhiITSMIokbadoutinz6InAcc = new TH1F("fHistPhiITSMIokbadoutinz6InAcc","Azimuthal distribution of ITSMI tracks with 6 layers OK; #phi; N tracks",100,0,2.*3.1415);
950 fHistPhiITSMIokbadoutinz6InAcc->Sumw2();
951 fHistPhiITSMIokbadoutinz6InAcc->SetMinimum(0);
952 fOutput->Add(fHistPhiITSMIokbadoutinz6InAcc);
8d63376d 953
3f18d0a6 954 fHistPtTPC = new TH1F("fHistPtTPC","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
8d63376d 955 fHistPtTPC->Sumw2();
956 fHistPtTPC->SetMinimum(0);
957 fOutput->Add(fHistPtTPC);
958
b7fb742d 959 fHistPtITSMI7 = new TH1F("fHistPtITSMI7","pt distribution of ITSMI7 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
960 fHistPtITSMI7->Sumw2();
961 fHistPtITSMI7->SetMinimum(0);
962 fOutput->Add(fHistPtITSMI7);
963
3f18d0a6 964 fHistPtITSMI6 = new TH1F("fHistPtITSMI6","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
8d63376d 965 fHistPtITSMI6->Sumw2();
966 fHistPtITSMI6->SetMinimum(0);
967 fOutput->Add(fHistPtITSMI6);
968
3f18d0a6 969 fHistPtITSMI5 = new TH1F("fHistPtITSMI5","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
8d63376d 970 fHistPtITSMI5->Sumw2();
971 fHistPtITSMI5->SetMinimum(0);
972 fOutput->Add(fHistPtITSMI5);
973
3f18d0a6 974 fHistPtITSMI4 = new TH1F("fHistPtITSMI4","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
8d63376d 975 fHistPtITSMI4->Sumw2();
976 fHistPtITSMI4->SetMinimum(0);
977 fOutput->Add(fHistPtITSMI4);
978
3f18d0a6 979 fHistPtITSMI3 = new TH1F("fHistPtITSMI3","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
8d63376d 980 fHistPtITSMI3->Sumw2();
981 fHistPtITSMI3->SetMinimum(0);
982 fOutput->Add(fHistPtITSMI3);
983
3f18d0a6 984 fHistPtITSMI2 = new TH1F("fHistPtITSMI2","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
8d63376d 985 fHistPtITSMI2->Sumw2();
986 fHistPtITSMI2->SetMinimum(0);
987 fOutput->Add(fHistPtITSMI2);
988
3f18d0a6 989 fHistPtITSMISPD = new TH1F("fHistPtITSMISPD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
8d63376d 990 fHistPtITSMISPD->Sumw2();
991 fHistPtITSMISPD->SetMinimum(0);
992 fOutput->Add(fHistPtITSMISPD);
3f18d0a6 993
7b4d74b2 994 fHistPtITSMIoneSPD = new TH1F("fHistPtITSMIoneSPD","pt distribution of ITSMIoneSPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
995 fHistPtITSMIoneSPD->Sumw2();
996 fHistPtITSMIoneSPD->SetMinimum(0);
997 fOutput->Add(fHistPtITSMIoneSPD);
998
b7fb742d 999 fHistPtITSMItwoSPD = new TH1F("fHistPtITSMItwoSPD","pt distribution of ITSMItwoSPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1000 fHistPtITSMItwoSPD->Sumw2();
1001 fHistPtITSMItwoSPD->SetMinimum(0);
1002 fOutput->Add(fHistPtITSMItwoSPD);
1003
3f18d0a6 1004 fHistPtTPCInAcc = new TH1F("fHistPtTPCInAcc","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1005 fHistPtTPCInAcc->Sumw2();
1006 fHistPtTPCInAcc->SetMinimum(0);
1007 fOutput->Add(fHistPtTPCInAcc);
c94b80d3 1008
1009 fHistPtTPCInAccTOFbc0 = new TH1F("fHistPtTPCInAccTOFbc0","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1010 fHistPtTPCInAccTOFbc0->Sumw2();
1011 fHistPtTPCInAccTOFbc0->SetMinimum(0);
1012 fOutput->Add(fHistPtTPCInAccTOFbc0);
31a96e36 1013
1014 fHistPtTPCInAccwSDD = new TH1F("fHistPtTPCInAccwSDD","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1015 fHistPtTPCInAccwSDD->Sumw2();
1016 fHistPtTPCInAccwSDD->SetMinimum(0);
1017 fOutput->Add(fHistPtTPCInAccwSDD);
1018
1019 fHistPtTPCInAccTOFbc0wSDD = new TH1F("fHistPtTPCInAccTOFbc0wSDD","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1020 fHistPtTPCInAccTOFbc0wSDD->Sumw2();
1021 fHistPtTPCInAccTOFbc0wSDD->SetMinimum(0);
1022 fOutput->Add(fHistPtTPCInAccTOFbc0wSDD);
1023
1024 fHistPtTPCInAccwoSDD = new TH1F("fHistPtTPCInAccwoSDD","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1025 fHistPtTPCInAccwoSDD->Sumw2();
1026 fHistPtTPCInAccwoSDD->SetMinimum(0);
1027 fOutput->Add(fHistPtTPCInAccwoSDD);
1028
1029 fHistPtTPCInAccTOFbc0woSDD = new TH1F("fHistPtTPCInAccTOFbc0woSDD","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1030 fHistPtTPCInAccTOFbc0woSDD->Sumw2();
1031 fHistPtTPCInAccTOFbc0woSDD->SetMinimum(0);
1032 fOutput->Add(fHistPtTPCInAccTOFbc0woSDD);
3f18d0a6 1033
dbf573fb 1034 fHistPtTPCInAccMCtwoSPD = new TH1F("fHistPtTPCInAccMCtwoSPD","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1035 fHistPtTPCInAccMCtwoSPD->Sumw2();
1036 fHistPtTPCInAccMCtwoSPD->SetMinimum(0);
1037 fOutput->Add(fHistPtTPCInAccMCtwoSPD);
1038
1039 fHistPtTPCInAccMConeSPD = new TH1F("fHistPtTPCInAccMConeSPD","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1040 fHistPtTPCInAccMConeSPD->Sumw2();
1041 fHistPtTPCInAccMConeSPD->SetMinimum(0);
1042 fOutput->Add(fHistPtTPCInAccMConeSPD);
1043
f8777463 1044 fHistdEdxVSPtTPCInAcc = new TH2F("fHistdEdxVSPtTPCInAcc","dE/dx distribution of TPC tracks; p_{t} [GeV/c]; TPC dE/dx",100,0,1,50,0,500);
1045 fOutput->Add(fHistdEdxVSPtTPCInAcc);
1046
1047 fHistdEdxVSPtITSTPCsel = new TH2F("fHistdEdxVSPtITSTPCsel","dE/dx distribution of TPC tracks; p_{t} [GeV/c]; ITS dE/dx",100,0,1,50,0,500);
1048 fOutput->Add(fHistdEdxVSPtITSTPCsel);
1049
0b2e8130 1050 fHistTPCclsVSPtTPCInAcc = new TH2F("fHistTPCclsVSPtTPCInAcc","TPC ncls; p_{t} [GeV/c]; nclsTPC",50,0,10,80,0,160);
1051 fOutput->Add(fHistTPCclsVSPtTPCInAcc);
1052 fHistTPCclsVSPtITSMISPDInAcc = new TH2F("fHistTPCclsVSPtITSMISPDInAcc","TPC ncls; p_{t} [GeV/c]; nclsTPC",50,0,10,80,0,160);
1053 fOutput->Add(fHistTPCclsVSPtITSMISPDInAcc);
1054
f8777463 1055 fHistPtVSphiTPCInAcc = new TH2F("fHistPtVSphiTPCInAcc","pt distribution of TPC tracks; phi; p_{t} [GeV/c]",18,0,6.28,20,0,0.5);
1056 fHistPtVSphiTPCInAcc->SetMinimum(0);
1057 fOutput->Add(fHistPtVSphiTPCInAcc);
1058
1059 fHistPtTPCInAccNoTRDout = new TH1F("fHistPtTPCInAccNoTRDout","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1060 fHistPtTPCInAccNoTRDout->Sumw2();
1061 fHistPtTPCInAccNoTRDout->SetMinimum(0);
1062 fOutput->Add(fHistPtTPCInAccNoTRDout);
1063
1064 fHistPtTPCInAccNoTOFout = new TH1F("fHistPtTPCInAccNoTOFout","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1065 fHistPtTPCInAccNoTOFout->Sumw2();
1066 fHistPtTPCInAccNoTOFout->SetMinimum(0);
1067 fOutput->Add(fHistPtTPCInAccNoTOFout);
1068
1069 fHistPtTPCInAccWithPtTPCAtVtx = new TH1F("fHistPtTPCInAccWithPtTPCAtVtx","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1070 fHistPtTPCInAccWithPtTPCAtVtx->Sumw2();
1071 fHistPtTPCInAccWithPtTPCAtVtx->SetMinimum(0);
1072 fOutput->Add(fHistPtTPCInAccWithPtTPCAtVtx);
1073
1074 fHistPtTPCInAccWithPtTPCAtInnerWall = new TH1F("fHistPtTPCInAccWithPtTPCAtInnerWall","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1075 fHistPtTPCInAccWithPtTPCAtInnerWall->Sumw2();
1076 fHistPtTPCInAccWithPtTPCAtInnerWall->SetMinimum(0);
1077 fOutput->Add(fHistPtTPCInAccWithPtTPCAtInnerWall);
1078
dbf573fb 1079 fHistDeltaPtTPC = new TH2F("fHistDeltaPtTPC","pt distribution of TPC tracks; p_{t} [GeV/c]; p_{t} TPC at vtx - p_{t} at inner wall [GeV/c]",100,0,1,50,-0.2,0.2);
f8777463 1080 fHistDeltaPtTPC->SetMinimum(0);
1081 fOutput->Add(fHistDeltaPtTPC);
1082
b7fb742d 1083 fHistPtITSMI7InAcc = new TH1F("fHistPtITSMI7InAcc","pt distribution of ITSMI7 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1084 fHistPtITSMI7InAcc->Sumw2();
1085 fHistPtITSMI7InAcc->SetMinimum(0);
1086 fOutput->Add(fHistPtITSMI7InAcc);
1087
3f18d0a6 1088 fHistPtITSMI6InAcc = new TH1F("fHistPtITSMI6InAcc","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1089 fHistPtITSMI6InAcc->Sumw2();
1090 fHistPtITSMI6InAcc->SetMinimum(0);
1091 fOutput->Add(fHistPtITSMI6InAcc);
1092
1093 fHistPtITSMI5InAcc = new TH1F("fHistPtITSMI5InAcc","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1094 fHistPtITSMI5InAcc->Sumw2();
1095 fHistPtITSMI5InAcc->SetMinimum(0);
1096 fOutput->Add(fHistPtITSMI5InAcc);
1097
1098 fHistPtITSMI4InAcc = new TH1F("fHistPtITSMI4InAcc","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1099 fHistPtITSMI4InAcc->Sumw2();
1100 fHistPtITSMI4InAcc->SetMinimum(0);
1101 fOutput->Add(fHistPtITSMI4InAcc);
8d63376d 1102
3f18d0a6 1103 fHistPtITSMI3InAcc = new TH1F("fHistPtITSMI3InAcc","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1104 fHistPtITSMI3InAcc->Sumw2();
1105 fHistPtITSMI3InAcc->SetMinimum(0);
1106 fOutput->Add(fHistPtITSMI3InAcc);
1107
1108 fHistPtITSMI2InAcc = new TH1F("fHistPtITSMI2InAcc","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1109 fHistPtITSMI2InAcc->Sumw2();
1110 fHistPtITSMI2InAcc->SetMinimum(0);
1111 fOutput->Add(fHistPtITSMI2InAcc);
1112
1113 fHistPtITSMISPDInAcc = new TH1F("fHistPtITSMISPDInAcc","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1114 fHistPtITSMISPDInAcc->Sumw2();
1115 fHistPtITSMISPDInAcc->SetMinimum(0);
1116 fOutput->Add(fHistPtITSMISPDInAcc);
1117
7b4d74b2 1118 fHistPtITSMIoneSPDInAcc = new TH1F("fHistPtITSMIoneSPDInAcc","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1119 fHistPtITSMIoneSPDInAcc->Sumw2();
1120 fHistPtITSMIoneSPDInAcc->SetMinimum(0);
1121 fOutput->Add(fHistPtITSMIoneSPDInAcc);
f8777463 1122
b7fb742d 1123 fHistPtITSMItwoSPDInAcc = new TH1F("fHistPtITSMItwoSPDInAcc","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1124 fHistPtITSMItwoSPDInAcc->Sumw2();
1125 fHistPtITSMItwoSPDInAcc->SetMinimum(0);
1126 fOutput->Add(fHistPtITSMItwoSPDInAcc);
1127
c94b80d3 1128 fHistPtITSMI6InAccTOFbc0 = new TH1F("fHistPtITSMI6InAccTOFbc0","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1129 fHistPtITSMI6InAccTOFbc0->Sumw2();
1130 fHistPtITSMI6InAccTOFbc0->SetMinimum(0);
1131 fOutput->Add(fHistPtITSMI6InAccTOFbc0);
1132
1133 fHistPtITSMI5InAccTOFbc0 = new TH1F("fHistPtITSMI5InAccTOFbc0","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1134 fHistPtITSMI5InAccTOFbc0->Sumw2();
1135 fHistPtITSMI5InAccTOFbc0->SetMinimum(0);
1136 fOutput->Add(fHistPtITSMI5InAccTOFbc0);
1137
1138 fHistPtITSMI4InAccTOFbc0 = new TH1F("fHistPtITSMI4InAccTOFbc0","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1139 fHistPtITSMI4InAccTOFbc0->Sumw2();
1140 fHistPtITSMI4InAccTOFbc0->SetMinimum(0);
1141 fOutput->Add(fHistPtITSMI4InAccTOFbc0);
1142
1143 fHistPtITSMI3InAccTOFbc0 = new TH1F("fHistPtITSMI3InAccTOFbc0","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1144 fHistPtITSMI3InAccTOFbc0->Sumw2();
1145 fHistPtITSMI3InAccTOFbc0->SetMinimum(0);
1146 fOutput->Add(fHistPtITSMI3InAccTOFbc0);
1147
1148 fHistPtITSMI2InAccTOFbc0 = new TH1F("fHistPtITSMI2InAccTOFbc0","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1149 fHistPtITSMI2InAccTOFbc0->Sumw2();
1150 fHistPtITSMI2InAccTOFbc0->SetMinimum(0);
1151 fOutput->Add(fHistPtITSMI2InAccTOFbc0);
1152
1153 fHistPtITSMISPDInAccTOFbc0 = new TH1F("fHistPtITSMISPDInAccTOFbc0","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1154 fHistPtITSMISPDInAccTOFbc0->Sumw2();
1155 fHistPtITSMISPDInAccTOFbc0->SetMinimum(0);
1156 fOutput->Add(fHistPtITSMISPDInAccTOFbc0);
1157
1158 fHistPtITSMIoneSPDInAccTOFbc0 = new TH1F("fHistPtITSMIoneSPDInAccTOFbc0","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1159 fHistPtITSMIoneSPDInAccTOFbc0->Sumw2();
1160 fHistPtITSMIoneSPDInAccTOFbc0->SetMinimum(0);
1161 fOutput->Add(fHistPtITSMIoneSPDInAccTOFbc0);
1162
31a96e36 1163 fHistPtITSMI6InAccwSDD = new TH1F("fHistPtITSMI6InAccwSDD","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1164 fHistPtITSMI6InAccwSDD->Sumw2();
1165 fHistPtITSMI6InAccwSDD->SetMinimum(0);
1166 fOutput->Add(fHistPtITSMI6InAccwSDD);
1167
1168 fHistPtITSMI5InAccwSDD = new TH1F("fHistPtITSMI5InAccwSDD","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1169 fHistPtITSMI5InAccwSDD->Sumw2();
1170 fHistPtITSMI5InAccwSDD->SetMinimum(0);
1171 fOutput->Add(fHistPtITSMI5InAccwSDD);
1172
1173 fHistPtITSMI4InAccwSDD = new TH1F("fHistPtITSMI4InAccwSDD","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1174 fHistPtITSMI4InAccwSDD->Sumw2();
1175 fHistPtITSMI4InAccwSDD->SetMinimum(0);
1176 fOutput->Add(fHistPtITSMI4InAccwSDD);
1177
1178 fHistPtITSMI3InAccwSDD = new TH1F("fHistPtITSMI3InAccwSDD","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1179 fHistPtITSMI3InAccwSDD->Sumw2();
1180 fHistPtITSMI3InAccwSDD->SetMinimum(0);
1181 fOutput->Add(fHistPtITSMI3InAccwSDD);
1182
1183 fHistPtITSMI2InAccwSDD = new TH1F("fHistPtITSMI2InAccwSDD","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1184 fHistPtITSMI2InAccwSDD->Sumw2();
1185 fHistPtITSMI2InAccwSDD->SetMinimum(0);
1186 fOutput->Add(fHistPtITSMI2InAccwSDD);
1187
1188 fHistPtITSMISPDInAccwSDD = new TH1F("fHistPtITSMISPDInAccwSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1189 fHistPtITSMISPDInAccwSDD->Sumw2();
1190 fHistPtITSMISPDInAccwSDD->SetMinimum(0);
1191 fOutput->Add(fHistPtITSMISPDInAccwSDD);
1192
1193 fHistPtITSMIoneSPDInAccwSDD = new TH1F("fHistPtITSMIoneSPDInAccwSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1194 fHistPtITSMIoneSPDInAccwSDD->Sumw2();
1195 fHistPtITSMIoneSPDInAccwSDD->SetMinimum(0);
1196 fOutput->Add(fHistPtITSMIoneSPDInAccwSDD);
1197
1198 fHistPtITSMI6InAccTOFbc0wSDD = new TH1F("fHistPtITSMI6InAccTOFbc0wSDD","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1199 fHistPtITSMI6InAccTOFbc0wSDD->Sumw2();
1200 fHistPtITSMI6InAccTOFbc0wSDD->SetMinimum(0);
1201 fOutput->Add(fHistPtITSMI6InAccTOFbc0wSDD);
1202
1203 fHistPtITSMI5InAccTOFbc0wSDD = new TH1F("fHistPtITSMI5InAccTOFbc0wSDD","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1204 fHistPtITSMI5InAccTOFbc0wSDD->Sumw2();
1205 fHistPtITSMI5InAccTOFbc0wSDD->SetMinimum(0);
1206 fOutput->Add(fHistPtITSMI5InAccTOFbc0wSDD);
1207
1208 fHistPtITSMI4InAccTOFbc0wSDD = new TH1F("fHistPtITSMI4InAccTOFbc0wSDD","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1209 fHistPtITSMI4InAccTOFbc0wSDD->Sumw2();
1210 fHistPtITSMI4InAccTOFbc0wSDD->SetMinimum(0);
1211 fOutput->Add(fHistPtITSMI4InAccTOFbc0wSDD);
1212
1213 fHistPtITSMI3InAccTOFbc0wSDD = new TH1F("fHistPtITSMI3InAccTOFbc0wSDD","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1214 fHistPtITSMI3InAccTOFbc0wSDD->Sumw2();
1215 fHistPtITSMI3InAccTOFbc0wSDD->SetMinimum(0);
1216 fOutput->Add(fHistPtITSMI3InAccTOFbc0wSDD);
1217
1218 fHistPtITSMI2InAccTOFbc0wSDD = new TH1F("fHistPtITSMI2InAccTOFbc0wSDD","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1219 fHistPtITSMI2InAccTOFbc0wSDD->Sumw2();
1220 fHistPtITSMI2InAccTOFbc0wSDD->SetMinimum(0);
1221 fOutput->Add(fHistPtITSMI2InAccTOFbc0wSDD);
1222
1223 fHistPtITSMISPDInAccTOFbc0wSDD = new TH1F("fHistPtITSMISPDInAccTOFbc0wSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1224 fHistPtITSMISPDInAccTOFbc0wSDD->Sumw2();
1225 fHistPtITSMISPDInAccTOFbc0wSDD->SetMinimum(0);
1226 fOutput->Add(fHistPtITSMISPDInAccTOFbc0wSDD);
1227
1228 fHistPtITSMIoneSPDInAccTOFbc0wSDD = new TH1F("fHistPtITSMIoneSPDInAccTOFbc0wSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1229 fHistPtITSMIoneSPDInAccTOFbc0wSDD->Sumw2();
1230 fHistPtITSMIoneSPDInAccTOFbc0wSDD->SetMinimum(0);
1231 fOutput->Add(fHistPtITSMIoneSPDInAccTOFbc0wSDD);
1232
1233 fHistPtITSMI6InAccwoSDD = new TH1F("fHistPtITSMI6InAccwoSDD","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1234 fHistPtITSMI6InAccwoSDD->Sumw2();
1235 fHistPtITSMI6InAccwoSDD->SetMinimum(0);
1236 fOutput->Add(fHistPtITSMI6InAccwoSDD);
1237
1238 fHistPtITSMI5InAccwoSDD = new TH1F("fHistPtITSMI5InAccwoSDD","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1239 fHistPtITSMI5InAccwoSDD->Sumw2();
1240 fHistPtITSMI5InAccwoSDD->SetMinimum(0);
1241 fOutput->Add(fHistPtITSMI5InAccwoSDD);
1242
1243 fHistPtITSMI4InAccwoSDD = new TH1F("fHistPtITSMI4InAccwoSDD","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1244 fHistPtITSMI4InAccwoSDD->Sumw2();
1245 fHistPtITSMI4InAccwoSDD->SetMinimum(0);
1246 fOutput->Add(fHistPtITSMI4InAccwoSDD);
1247
1248 fHistPtITSMI3InAccwoSDD = new TH1F("fHistPtITSMI3InAccwoSDD","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1249 fHistPtITSMI3InAccwoSDD->Sumw2();
1250 fHistPtITSMI3InAccwoSDD->SetMinimum(0);
1251 fOutput->Add(fHistPtITSMI3InAccwoSDD);
1252
1253 fHistPtITSMI2InAccwoSDD = new TH1F("fHistPtITSMI2InAccwoSDD","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1254 fHistPtITSMI2InAccwoSDD->Sumw2();
1255 fHistPtITSMI2InAccwoSDD->SetMinimum(0);
1256 fOutput->Add(fHistPtITSMI2InAccwoSDD);
1257
1258 fHistPtITSMISPDInAccwoSDD = new TH1F("fHistPtITSMISPDInAccwoSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1259 fHistPtITSMISPDInAccwoSDD->Sumw2();
1260 fHistPtITSMISPDInAccwoSDD->SetMinimum(0);
1261 fOutput->Add(fHistPtITSMISPDInAccwoSDD);
1262
1263 fHistPtITSMIoneSPDInAccwoSDD = new TH1F("fHistPtITSMIoneSPDInAccwoSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1264 fHistPtITSMIoneSPDInAccwoSDD->Sumw2();
1265 fHistPtITSMIoneSPDInAccwoSDD->SetMinimum(0);
1266 fOutput->Add(fHistPtITSMIoneSPDInAccwoSDD);
1267
1268 fHistPtITSMI6InAccTOFbc0woSDD = new TH1F("fHistPtITSMI6InAccTOFbc0woSDD","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1269 fHistPtITSMI6InAccTOFbc0woSDD->Sumw2();
1270 fHistPtITSMI6InAccTOFbc0woSDD->SetMinimum(0);
1271 fOutput->Add(fHistPtITSMI6InAccTOFbc0woSDD);
1272
1273 fHistPtITSMI5InAccTOFbc0woSDD = new TH1F("fHistPtITSMI5InAccTOFbc0woSDD","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1274 fHistPtITSMI5InAccTOFbc0woSDD->Sumw2();
1275 fHistPtITSMI5InAccTOFbc0woSDD->SetMinimum(0);
1276 fOutput->Add(fHistPtITSMI5InAccTOFbc0woSDD);
1277
1278 fHistPtITSMI4InAccTOFbc0woSDD = new TH1F("fHistPtITSMI4InAccTOFbc0woSDD","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1279 fHistPtITSMI4InAccTOFbc0woSDD->Sumw2();
1280 fHistPtITSMI4InAccTOFbc0woSDD->SetMinimum(0);
1281 fOutput->Add(fHistPtITSMI4InAccTOFbc0woSDD);
1282
1283 fHistPtITSMI3InAccTOFbc0woSDD = new TH1F("fHistPtITSMI3InAccTOFbc0woSDD","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1284 fHistPtITSMI3InAccTOFbc0woSDD->Sumw2();
1285 fHistPtITSMI3InAccTOFbc0woSDD->SetMinimum(0);
1286 fOutput->Add(fHistPtITSMI3InAccTOFbc0woSDD);
1287
1288 fHistPtITSMI2InAccTOFbc0woSDD = new TH1F("fHistPtITSMI2InAccTOFbc0woSDD","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1289 fHistPtITSMI2InAccTOFbc0woSDD->Sumw2();
1290 fHistPtITSMI2InAccTOFbc0woSDD->SetMinimum(0);
1291 fOutput->Add(fHistPtITSMI2InAccTOFbc0woSDD);
1292
1293 fHistPtITSMISPDInAccTOFbc0woSDD = new TH1F("fHistPtITSMISPDInAccTOFbc0woSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1294 fHistPtITSMISPDInAccTOFbc0woSDD->Sumw2();
1295 fHistPtITSMISPDInAccTOFbc0woSDD->SetMinimum(0);
1296 fOutput->Add(fHistPtITSMISPDInAccTOFbc0woSDD);
1297
1298 fHistPtITSMIoneSPDInAccTOFbc0woSDD = new TH1F("fHistPtITSMIoneSPDInAccTOFbc0woSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1299 fHistPtITSMIoneSPDInAccTOFbc0woSDD->Sumw2();
1300 fHistPtITSMIoneSPDInAccTOFbc0woSDD->SetMinimum(0);
1301 fOutput->Add(fHistPtITSMIoneSPDInAccTOFbc0woSDD);
1302
0b2e8130 1303 fHistPtITSMIoneSPDInAccShared = new TH1F("fHistPtITSMIoneSPDInAccShared","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1304 fHistPtITSMIoneSPDInAccShared->Sumw2();
1305 fHistPtITSMIoneSPDInAccShared->SetMinimum(0);
1306 fOutput->Add(fHistPtITSMIoneSPDInAccShared);
1307
1308 fHistPtITSMIoneSPDInAccSharedSPD = new TH1F("fHistPtITSMIoneSPDInAccSharedSPD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1309 fHistPtITSMIoneSPDInAccSharedSPD->Sumw2();
1310 fHistPtITSMIoneSPDInAccSharedSPD->SetMinimum(0);
1311 fOutput->Add(fHistPtITSMIoneSPDInAccSharedSPD);
1312
1313 fHistPtITSMISPD1InAccShared = new TH1F("fHistPtITSMISPD1InAccShared","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1314 fHistPtITSMISPD1InAccShared->Sumw2();
1315 fHistPtITSMISPD1InAccShared->SetMinimum(0);
1316 fOutput->Add(fHistPtITSMISPD1InAccShared);
1317
1318 fHistPtITSMISPD2InAccShared = new TH1F("fHistPtITSMISPD2InAccShared","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1319 fHistPtITSMISPD2InAccShared->Sumw2();
1320 fHistPtITSMISPD2InAccShared->SetMinimum(0);
1321 fOutput->Add(fHistPtITSMISPD2InAccShared);
1322
1323 fHistPtITSMIoneSPDInAccSharedFake = new TH1F("fHistPtITSMIoneSPDInAccSharedFake","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1324 fHistPtITSMIoneSPDInAccSharedFake->Sumw2();
1325 fHistPtITSMIoneSPDInAccSharedFake->SetMinimum(0);
1326 fOutput->Add(fHistPtITSMIoneSPDInAccSharedFake);
1327
1328 fHistPtITSMIoneSPDInAccSharedSPDFake = new TH1F("fHistPtITSMIoneSPDInAccSharedSPDFake","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1329 fHistPtITSMIoneSPDInAccSharedSPDFake->Sumw2();
1330 fHistPtITSMIoneSPDInAccSharedSPDFake->SetMinimum(0);
1331 fOutput->Add(fHistPtITSMIoneSPDInAccSharedSPDFake);
1332
b7fb742d 1333 fHistPhiITSMI7InAcc = new TH1F("fHistPhiITSMI7InAcc","phi distribution of ITSMI7 tracks; #phi; N tracks",100,0,2*3.1415);
1334 fHistPhiITSMI7InAcc->Sumw2();
1335 fHistPhiITSMI7InAcc->SetMinimum(0);
1336 fOutput->Add(fHistPhiITSMI7InAcc);
1337
547640ff 1338 fHistPhiITSMI6InAcc = new TH1F("fHistPhiITSMI6InAcc","phi distribution of ITSMI6 tracks; #phi; N tracks",100,0,2*3.1415);
1339 fHistPhiITSMI6InAcc->Sumw2();
1340 fHistPhiITSMI6InAcc->SetMinimum(0);
1341 fOutput->Add(fHistPhiITSMI6InAcc);
1342
1343 fHistPhiITSMI5InAcc = new TH1F("fHistPhiITSMI5InAcc","phi distribution of ITSMI5 tracks; #phi; N tracks",100,0,2*3.1415);
1344 fHistPhiITSMI5InAcc->Sumw2();
1345 fHistPhiITSMI5InAcc->SetMinimum(0);
1346 fOutput->Add(fHistPhiITSMI5InAcc);
1347
1348 fHistPhiITSMI4InAcc = new TH1F("fHistPhiITSMI4InAcc","phi distribution of ITSMI4 tracks; #phi; N tracks",100,0,2*3.1415);
1349 fHistPhiITSMI4InAcc->Sumw2();
1350 fHistPhiITSMI4InAcc->SetMinimum(0);
1351 fOutput->Add(fHistPhiITSMI4InAcc);
1352
1353 fHistPhiITSMI3InAcc = new TH1F("fHistPhiITSMI3InAcc","phi distribution of ITSMI3 tracks; #phi; N tracks",100,0,2*3.1415);
1354 fHistPhiITSMI3InAcc->Sumw2();
1355 fHistPhiITSMI3InAcc->SetMinimum(0);
1356 fOutput->Add(fHistPhiITSMI3InAcc);
1357
1358 fHistPhiITSMI2InAcc = new TH1F("fHistPhiITSMI2InAcc","phi distribution of ITSMI2 tracks; #phi; N tracks",100,0,2*3.1415);
1359 fHistPhiITSMI2InAcc->Sumw2();
1360 fHistPhiITSMI2InAcc->SetMinimum(0);
1361 fOutput->Add(fHistPhiITSMI2InAcc);
1362
1363 fHistPhiITSMISPDInAcc = new TH1F("fHistPhiITSMISPDInAcc","phi distribution of ITSMISPD tracks; #phi; N tracks",100,0,2*3.1415);
1364 fHistPhiITSMISPDInAcc->Sumw2();
1365 fHistPhiITSMISPDInAcc->SetMinimum(0);
1366 fOutput->Add(fHistPhiITSMISPDInAcc);
1367
1368 fHistPhiITSMIoneSPDInAcc = new TH1F("fHistPhiITSMIoneSPDInAcc","phi distribution of ITSMISPD tracks; #phi; N tracks",100,0,2*3.1415);
1369 fHistPhiITSMIoneSPDInAcc->Sumw2();
1370 fHistPhiITSMIoneSPDInAcc->SetMinimum(0);
1371 fOutput->Add(fHistPhiITSMIoneSPDInAcc);
1372
b7fb742d 1373 fHistPhiITSMItwoSPDInAcc = new TH1F("fHistPhiITSMItwoSPDInAcc","phi distribution of ITSMISPD tracks; #phi; N tracks",100,0,2*3.1415);
1374 fHistPhiITSMItwoSPDInAcc->Sumw2();
1375 fHistPhiITSMItwoSPDInAcc->SetMinimum(0);
1376 fOutput->Add(fHistPhiITSMItwoSPDInAcc);
1377
547640ff 1378 fHistEtaITSMI6InAcc = new TH1F("fHistEtaITSMI6InAcc","eta distribution of ITSMI6 tracks; #eta; N tracks",100,-1.5,+1.5);
1379 fHistEtaITSMI6InAcc->Sumw2();
1380 fHistEtaITSMI6InAcc->SetMinimum(0);
1381 fOutput->Add(fHistEtaITSMI6InAcc);
b7fb742d 1382
1383 fHistEtaITSMI7InAcc = new TH1F("fHistEtaITSMI7InAcc","eta distribution of ITSMI7 tracks; #eta; N tracks",100,-1.5,+1.5);
1384 fHistEtaITSMI7InAcc->Sumw2();
1385 fHistEtaITSMI7InAcc->SetMinimum(0);
1386 fOutput->Add(fHistEtaITSMI7InAcc);
547640ff 1387
1388 fHistEtaITSMI5InAcc = new TH1F("fHistEtaITSMI5InAcc","eta distribution of ITSMI5 tracks; #eta; N tracks",100,-1.5,+1.5);
1389 fHistEtaITSMI5InAcc->Sumw2();
1390 fHistEtaITSMI5InAcc->SetMinimum(0);
1391 fOutput->Add(fHistEtaITSMI5InAcc);
1392
1393 fHistEtaITSMI4InAcc = new TH1F("fHistEtaITSMI4InAcc","eta distribution of ITSMI4 tracks; #eta; N tracks",100,-1.5,+1.5);
1394 fHistEtaITSMI4InAcc->Sumw2();
1395 fHistEtaITSMI4InAcc->SetMinimum(0);
1396 fOutput->Add(fHistEtaITSMI4InAcc);
1397
1398 fHistEtaITSMI3InAcc = new TH1F("fHistEtaITSMI3InAcc","eta distribution of ITSMI3 tracks; #eta; N tracks",100,-1.5,+1.5);
1399 fHistEtaITSMI3InAcc->Sumw2();
1400 fHistEtaITSMI3InAcc->SetMinimum(0);
1401 fOutput->Add(fHistEtaITSMI3InAcc);
1402
1403 fHistEtaITSMI2InAcc = new TH1F("fHistEtaITSMI2InAcc","eta distribution of ITSMI2 tracks; #eta; N tracks",100,-1.5,+1.5);
1404 fHistEtaITSMI2InAcc->Sumw2();
1405 fHistEtaITSMI2InAcc->SetMinimum(0);
1406 fOutput->Add(fHistEtaITSMI2InAcc);
1407
1408 fHistEtaITSMISPDInAcc = new TH1F("fHistEtaITSMISPDInAcc","eta distribution of ITSMISPD tracks; #eta; N tracks",100,-1.5,+1.5);
1409 fHistEtaITSMISPDInAcc->Sumw2();
1410 fHistEtaITSMISPDInAcc->SetMinimum(0);
1411 fOutput->Add(fHistEtaITSMISPDInAcc);
1412
1413 fHistEtaITSMIoneSPDInAcc = new TH1F("fHistEtaITSMIoneSPDInAcc","eta distribution of ITSMISPD tracks; #eta; N tracks",100,-1.5,+1.5);
1414 fHistEtaITSMIoneSPDInAcc->Sumw2();
1415 fHistEtaITSMIoneSPDInAcc->SetMinimum(0);
1416 fOutput->Add(fHistEtaITSMIoneSPDInAcc);
1417
b7fb742d 1418 fHistEtaITSMItwoSPDInAcc = new TH1F("fHistEtaITSMItwoSPDInAcc","eta distribution of ITSMISPD tracks; #eta; N tracks",100,-1.5,+1.5);
1419 fHistEtaITSMItwoSPDInAcc->Sumw2();
1420 fHistEtaITSMItwoSPDInAcc->SetMinimum(0);
1421 fOutput->Add(fHistEtaITSMItwoSPDInAcc);
1422
547640ff 1423
b7fb742d 1424 fHistPtITSMI7InAccFake = new TH1F("fHistPtITSMI7InAccFake","pt distribution of ITSMI7 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1425 fHistPtITSMI7InAccFake->Sumw2();
1426 fHistPtITSMI7InAccFake->SetMinimum(0);
1427 fOutput->Add(fHistPtITSMI7InAccFake);
1428
0faeb7ba 1429 fHistPtITSMI6InAccFake = new TH1F("fHistPtITSMI6InAccFake","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1430 fHistPtITSMI6InAccFake->Sumw2();
1431 fHistPtITSMI6InAccFake->SetMinimum(0);
1432 fOutput->Add(fHistPtITSMI6InAccFake);
1433
1434 fHistPtITSMI5InAccFake = new TH1F("fHistPtITSMI5InAccFake","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1435 fHistPtITSMI5InAccFake->Sumw2();
1436 fHistPtITSMI5InAccFake->SetMinimum(0);
1437 fOutput->Add(fHistPtITSMI5InAccFake);
1438
1439 fHistPtITSMI4InAccFake = new TH1F("fHistPtITSMI4InAccFake","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1440 fHistPtITSMI4InAccFake->Sumw2();
1441 fHistPtITSMI4InAccFake->SetMinimum(0);
1442 fOutput->Add(fHistPtITSMI4InAccFake);
1443
1444 fHistPtITSMI3InAccFake = new TH1F("fHistPtITSMI3InAccFake","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1445 fHistPtITSMI3InAccFake->Sumw2();
1446 fHistPtITSMI3InAccFake->SetMinimum(0);
1447 fOutput->Add(fHistPtITSMI3InAccFake);
1448
1449 fHistPtITSMI2InAccFake = new TH1F("fHistPtITSMI2InAccFake","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1450 fHistPtITSMI2InAccFake->Sumw2();
1451 fHistPtITSMI2InAccFake->SetMinimum(0);
1452 fOutput->Add(fHistPtITSMI2InAccFake);
1453
1454 fHistPtITSMISPDInAccFake = new TH1F("fHistPtITSMISPDInAccFake","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1455 fHistPtITSMISPDInAccFake->Sumw2();
1456 fHistPtITSMISPDInAccFake->SetMinimum(0);
1457 fOutput->Add(fHistPtITSMISPDInAccFake);
1458
1459 fHistPtITSMIoneSPDInAccFake = new TH1F("fHistPtITSMIoneSPDInAccFake","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1460 fHistPtITSMIoneSPDInAccFake->Sumw2();
1461 fHistPtITSMIoneSPDInAccFake->SetMinimum(0);
1462 fOutput->Add(fHistPtITSMIoneSPDInAccFake);
1463
b7fb742d 1464 fHistPtITSMItwoSPDInAccFake = new TH1F("fHistPtITSMItwoSPDInAccFake","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1465 fHistPtITSMItwoSPDInAccFake->Sumw2();
1466 fHistPtITSMItwoSPDInAccFake->SetMinimum(0);
1467 fOutput->Add(fHistPtITSMItwoSPDInAccFake);
1468
6c6fab7b 1469 fHistPtITSMIoneSPDthreeSDDSSDInAcc = new TH1F("fHistPtITSMIoneSPDthreeSDDSSDInAcc","pt distribution of ITSMI tracks (>0 in SPD, >2 in SDD+SSD); p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1470 fHistPtITSMIoneSPDthreeSDDSSDInAcc->Sumw2();
1471 fHistPtITSMIoneSPDthreeSDDSSDInAcc->SetMinimum(0);
1472 fOutput->Add(fHistPtITSMIoneSPDthreeSDDSSDInAcc);
1473
f8777463 1474 fHistPtITSTPCsel = new TH1F("fHistPtITSTPCsel","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1475 fHistPtITSTPCsel->Sumw2();
1476 fHistPtITSTPCsel->SetMinimum(0);
1477 fOutput->Add(fHistPtITSTPCsel);
1478
c94b80d3 1479 fHistPtITSTPCselTOFbc0 = new TH1F("fHistPtITSTPCselTOFbc0","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1480 fHistPtITSTPCselTOFbc0->Sumw2();
1481 fHistPtITSTPCselTOFbc0->SetMinimum(0);
1482 fOutput->Add(fHistPtITSTPCselTOFbc0);
1483
31a96e36 1484 fHistPtITSTPCselwSDD = new TH1F("fHistPtITSTPCselwSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1485 fHistPtITSTPCselwSDD->Sumw2();
1486 fHistPtITSTPCselwSDD->SetMinimum(0);
1487 fOutput->Add(fHistPtITSTPCselwSDD);
1488
1489 fHistPtITSTPCselTOFbc0wSDD = new TH1F("fHistPtITSTPCselTOFbc0wSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1490 fHistPtITSTPCselTOFbc0wSDD->Sumw2();
1491 fHistPtITSTPCselTOFbc0wSDD->SetMinimum(0);
1492 fOutput->Add(fHistPtITSTPCselTOFbc0wSDD);
1493
1494 fHistPtITSTPCselwoSDD = new TH1F("fHistPtITSTPCselwoSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1495 fHistPtITSTPCselwoSDD->Sumw2();
1496 fHistPtITSTPCselwoSDD->SetMinimum(0);
1497 fOutput->Add(fHistPtITSTPCselwoSDD);
1498
1499 fHistPtITSTPCselTOFbc0woSDD = new TH1F("fHistPtITSTPCselTOFbc0woSDD","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1500 fHistPtITSTPCselTOFbc0woSDD->Sumw2();
1501 fHistPtITSTPCselTOFbc0woSDD->SetMinimum(0);
1502 fOutput->Add(fHistPtITSTPCselTOFbc0woSDD);
1503
f8777463 1504 fHistPtITSTPCselP = new TH1F("fHistPtITSTPCselP","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1505 fHistPtITSTPCselP->Sumw2();
1506 fHistPtITSTPCselP->SetMinimum(0);
1507 fOutput->Add(fHistPtITSTPCselP);
1508
1509 fHistPtITSTPCselS = new TH1F("fHistPtITSTPCselS","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1510 fHistPtITSTPCselS->Sumw2();
1511 fHistPtITSTPCselS->SetMinimum(0);
1512 fOutput->Add(fHistPtITSTPCselS);
1513
af874c53 1514 fHistPtITSTPCselFake = new TH1F("fHistPtITSTPCselFake","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1515 fHistPtITSTPCselFake->Sumw2();
1516 fHistPtITSTPCselFake->SetMinimum(0);
1517 fOutput->Add(fHistPtITSTPCselFake);
1518
f8777463 1519 fHistPtITSTPCselSfromStrange = new TH1F("fHistPtITSTPCselSfromStrange","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1520 fHistPtITSTPCselSfromStrange->Sumw2();
1521 fHistPtITSTPCselSfromStrange->SetMinimum(0);
1522 fOutput->Add(fHistPtITSTPCselSfromStrange);
1523
1524 fHistPtITSTPCselPfromStrange = new TH1F("fHistPtITSTPCselPfromStrange","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1525 fHistPtITSTPCselPfromStrange->Sumw2();
1526 fHistPtITSTPCselPfromStrange->SetMinimum(0);
1527 fOutput->Add(fHistPtITSTPCselPfromStrange);
1528
1529 fHistPtITSTPCselSfromMat = new TH1F("fHistPtITSTPCselSfromMat","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1530 fHistPtITSTPCselSfromMat->Sumw2();
1531 fHistPtITSTPCselSfromMat->SetMinimum(0);
1532 fOutput->Add(fHistPtITSTPCselSfromMat);
1533
1534 fHistPtTPCInAccP = new TH1F("fHistPtTPCInAccP","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1535 fHistPtTPCInAccP->Sumw2();
1536 fHistPtTPCInAccP->SetMinimum(0);
1537 fOutput->Add(fHistPtTPCInAccP);
1538
1539 fHistPtTPCInAccPfromStrange = new TH1F("fHistPtTPCInAccPfromStrange","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1540 fHistPtTPCInAccPfromStrange->Sumw2();
1541 fHistPtTPCInAccPfromStrange->SetMinimum(0);
1542 fOutput->Add(fHistPtTPCInAccPfromStrange);
1543
b7fb742d 1544 fHistPtITSMI7InAccP = new TH1F("fHistPtITSMI7InAccP","pt distribution of ITSMI7 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1545 fHistPtITSMI7InAccP->Sumw2();
1546 fHistPtITSMI7InAccP->SetMinimum(0);
1547 fOutput->Add(fHistPtITSMI7InAccP);
1548
f8777463 1549 fHistPtITSMI6InAccP = new TH1F("fHistPtITSMI6InAccP","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1550 fHistPtITSMI6InAccP->Sumw2();
1551 fHistPtITSMI6InAccP->SetMinimum(0);
1552 fOutput->Add(fHistPtITSMI6InAccP);
1553
1554 fHistPtITSMI5InAccP = new TH1F("fHistPtITSMI5InAccP","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1555 fHistPtITSMI5InAccP->Sumw2();
1556 fHistPtITSMI5InAccP->SetMinimum(0);
1557 fOutput->Add(fHistPtITSMI5InAccP);
1558
1559 fHistPtITSMI4InAccP = new TH1F("fHistPtITSMI4InAccP","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1560 fHistPtITSMI4InAccP->Sumw2();
1561 fHistPtITSMI4InAccP->SetMinimum(0);
1562 fOutput->Add(fHistPtITSMI4InAccP);
1563
1564 fHistPtITSMI3InAccP = new TH1F("fHistPtITSMI3InAccP","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1565 fHistPtITSMI3InAccP->Sumw2();
1566 fHistPtITSMI3InAccP->SetMinimum(0);
1567 fOutput->Add(fHistPtITSMI3InAccP);
1568
1569 fHistPtITSMI2InAccP = new TH1F("fHistPtITSMI2InAccP","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1570 fHistPtITSMI2InAccP->Sumw2();
1571 fHistPtITSMI2InAccP->SetMinimum(0);
1572 fOutput->Add(fHistPtITSMI2InAccP);
1573
1574 fHistPtITSMISPDInAccP = new TH1F("fHistPtITSMISPDInAccP","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1575 fHistPtITSMISPDInAccP->Sumw2();
1576 fHistPtITSMISPDInAccP->SetMinimum(0);
1577 fOutput->Add(fHistPtITSMISPDInAccP);
1578
1579 fHistPtITSMIoneSPDInAccP = new TH1F("fHistPtITSMIoneSPDInAccP","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1580 fHistPtITSMIoneSPDInAccP->Sumw2();
1581 fHistPtITSMIoneSPDInAccP->SetMinimum(0);
1582 fOutput->Add(fHistPtITSMIoneSPDInAccP);
1583
b7fb742d 1584 fHistPtITSMItwoSPDInAccP = new TH1F("fHistPtITSMItwoSPDInAccP","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1585 fHistPtITSMItwoSPDInAccP->Sumw2();
1586 fHistPtITSMItwoSPDInAccP->SetMinimum(0);
1587 fOutput->Add(fHistPtITSMItwoSPDInAccP);
1588
f8777463 1589 fHistPtTPCInAccS = new TH1F("fHistPtTPCInAccS","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1590 fHistPtTPCInAccS->Sumw2();
1591 fHistPtTPCInAccS->SetMinimum(0);
1592 fOutput->Add(fHistPtTPCInAccS);
1593
1594 fHistPtTPCInAccSfromStrange = new TH1F("fHistPtTPCInAccSfromStrange","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1595 fHistPtTPCInAccSfromStrange->Sumw2();
1596 fHistPtTPCInAccSfromStrange->SetMinimum(0);
1597 fOutput->Add(fHistPtTPCInAccSfromStrange);
1598
1599 fHistPtTPCInAccSfromMat = new TH1F("fHistPtTPCInAccSfromMat","pt distribution of TPC tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1600 fHistPtTPCInAccSfromMat->Sumw2();
1601 fHistPtTPCInAccSfromMat->SetMinimum(0);
1602 fOutput->Add(fHistPtTPCInAccSfromMat);
1603
b7fb742d 1604 fHistPtITSMI7InAccS = new TH1F("fHistPtITSMI7InAccS","pt distribution of ITSMI7 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1605 fHistPtITSMI7InAccS->Sumw2();
1606 fHistPtITSMI7InAccS->SetMinimum(0);
1607 fOutput->Add(fHistPtITSMI7InAccS);
1608
f8777463 1609 fHistPtITSMI6InAccS = new TH1F("fHistPtITSMI6InAccS","pt distribution of ITSMI6 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1610 fHistPtITSMI6InAccS->Sumw2();
1611 fHistPtITSMI6InAccS->SetMinimum(0);
1612 fOutput->Add(fHistPtITSMI6InAccS);
1613
1614 fHistPtITSMI5InAccS = new TH1F("fHistPtITSMI5InAccS","pt distribution of ITSMI5 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1615 fHistPtITSMI5InAccS->Sumw2();
1616 fHistPtITSMI5InAccS->SetMinimum(0);
1617 fOutput->Add(fHistPtITSMI5InAccS);
1618
1619 fHistPtITSMI4InAccS = new TH1F("fHistPtITSMI4InAccS","pt distribution of ITSMI4 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1620 fHistPtITSMI4InAccS->Sumw2();
1621 fHistPtITSMI4InAccS->SetMinimum(0);
1622 fOutput->Add(fHistPtITSMI4InAccS);
1623
1624 fHistPtITSMI3InAccS = new TH1F("fHistPtITSMI3InAccS","pt distribution of ITSMI3 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1625 fHistPtITSMI3InAccS->Sumw2();
1626 fHistPtITSMI3InAccS->SetMinimum(0);
1627 fOutput->Add(fHistPtITSMI3InAccS);
1628
1629 fHistPtITSMI2InAccS = new TH1F("fHistPtITSMI2InAccS","pt distribution of ITSMI2 tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1630 fHistPtITSMI2InAccS->Sumw2();
1631 fHistPtITSMI2InAccS->SetMinimum(0);
1632 fOutput->Add(fHistPtITSMI2InAccS);
1633
1634 fHistPtITSMISPDInAccS = new TH1F("fHistPtITSMISPDInAccS","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1635 fHistPtITSMISPDInAccS->Sumw2();
1636 fHistPtITSMISPDInAccS->SetMinimum(0);
1637 fOutput->Add(fHistPtITSMISPDInAccS);
1638
1639 fHistPtITSMIoneSPDInAccS = new TH1F("fHistPtITSMIoneSPDInAccS","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1640 fHistPtITSMIoneSPDInAccS->Sumw2();
1641 fHistPtITSMIoneSPDInAccS->SetMinimum(0);
1642 fOutput->Add(fHistPtITSMIoneSPDInAccS);
b7fb742d 1643
1644 fHistPtITSMItwoSPDInAccS = new TH1F("fHistPtITSMItwoSPDInAccS","pt distribution of ITSMISPD tracks; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1645 fHistPtITSMItwoSPDInAccS->Sumw2();
1646 fHistPtITSMItwoSPDInAccS->SetMinimum(0);
1647 fOutput->Add(fHistPtITSMItwoSPDInAccS);
7b4d74b2 1648
3f18d0a6 1649 fHistPtITSMIokbadoutinz6 = new TH1F("fHistPtITSMIokbadoutinz6","pt distribution of ITSMI tracks with 6 layers OK; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1650 fHistPtITSMIokbadoutinz6->Sumw2();
1651 fHistPtITSMIokbadoutinz6->SetMinimum(0);
1652 fOutput->Add(fHistPtITSMIokbadoutinz6);
8d63376d 1653
c2893e82 1654 fHistPtITSMIokbadoutinz4InAcc = new TH1F("fHistPtITSMIokbadoutinz4InAcc","pt distribution of ITSMI tracks with 4 layers OK; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1655 fHistPtITSMIokbadoutinz4InAcc->Sumw2();
1656 fHistPtITSMIokbadoutinz4InAcc->SetMinimum(0);
1657 fOutput->Add(fHistPtITSMIokbadoutinz4InAcc);
1658
1659 fHistPtITSMIokbadoutinz5InAcc = new TH1F("fHistPtITSMIokbadoutinz5InAcc","pt distribution of ITSMI tracks with 5 layers OK; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1660 fHistPtITSMIokbadoutinz5InAcc->Sumw2();
1661 fHistPtITSMIokbadoutinz5InAcc->SetMinimum(0);
1662 fOutput->Add(fHistPtITSMIokbadoutinz5InAcc);
1663
3f18d0a6 1664 fHistPtITSMIokbadoutinz6InAcc = new TH1F("fHistPtITSMIokbadoutinz6InAcc","pt distribution of ITSMI tracks with 6 layers OK; p_{t} [GeV/c]; N tracks",nPtBins,xPtBins);
1665 fHistPtITSMIokbadoutinz6InAcc->Sumw2();
1666 fHistPtITSMIokbadoutinz6InAcc->SetMinimum(0);
1667 fOutput->Add(fHistPtITSMIokbadoutinz6InAcc);
8d63376d 1668
f8777463 1669 fHistRProdVtxInAccP = new TH1F("fHistRProdVtxInAccP","Radius of production vertex for primaries; r [cm]; N tracks",100,0,10);
1670 fHistRProdVtxInAccP->Sumw2();
1671 fHistRProdVtxInAccP->SetMinimum(0);
1672 fOutput->Add(fHistRProdVtxInAccP);
1673
1674 fHistRProdVtxInAccS = new TH1F("fHistRProdVtxInAccS","Radius of production vertex for secondaries; r [cm]; N tracks",100,0,10);
1675 fHistRProdVtxInAccS->Sumw2();
1676 fHistRProdVtxInAccS->SetMinimum(0);
1677 fOutput->Add(fHistRProdVtxInAccS);
1678
0b2e8130 1679
f8777463 1680 fHistd0rphiTPCInAccP150200 = new TH1F("fHistd0rphiTPCInAccP150200","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-5,5);
1681 fHistd0rphiTPCInAccP150200->Sumw2();
1682 fHistd0rphiTPCInAccP150200->SetMinimum(0);
1683 fOutput->Add(fHistd0rphiTPCInAccP150200);
1684
1685 fHistd0rphiTPCInAccP500700 = new TH1F("fHistd0rphiTPCInAccP500700","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-5,5);
1686 fHistd0rphiTPCInAccP500700->Sumw2();
1687 fHistd0rphiTPCInAccP500700->SetMinimum(0);
1688 fOutput->Add(fHistd0rphiTPCInAccP500700);
1689
1690 fHistd0rphiTPCInAccP10001500 = new TH1F("fHistd0rphiTPCInAccP10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-5,5);
1691 fHistd0rphiTPCInAccP10001500->Sumw2();
1692 fHistd0rphiTPCInAccP10001500->SetMinimum(0);
1693 fOutput->Add(fHistd0rphiTPCInAccP10001500);
1694
1695 fHistd0rphiTPCInAccS150200 = new TH1F("fHistd0rphiTPCInAccS150200","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-5,5);
1696 fHistd0rphiTPCInAccS150200->Sumw2();
1697 fHistd0rphiTPCInAccS150200->SetMinimum(0);
1698 fOutput->Add(fHistd0rphiTPCInAccS150200);
1699
1700 fHistd0rphiTPCInAccS500700 = new TH1F("fHistd0rphiTPCInAccS500700","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-5,5);
1701 fHistd0rphiTPCInAccS500700->Sumw2();
1702 fHistd0rphiTPCInAccS500700->SetMinimum(0);
1703 fOutput->Add(fHistd0rphiTPCInAccS500700);
1704
1705 fHistd0rphiTPCInAccS10001500 = new TH1F("fHistd0rphiTPCInAccS10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-5,5);
1706 fHistd0rphiTPCInAccS10001500->Sumw2();
1707 fHistd0rphiTPCInAccS10001500->SetMinimum(0);
1708 fOutput->Add(fHistd0rphiTPCInAccS10001500);
1709
1710 fHistd0rphiITSMISPDInAccP150200 = new TH1F("fHistd0rphiITSMISPDInAccP150200","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1711 fHistd0rphiITSMISPDInAccP150200->Sumw2();
1712 fHistd0rphiITSMISPDInAccP150200->SetMinimum(0);
1713 fOutput->Add(fHistd0rphiITSMISPDInAccP150200);
1714
1715 fHistd0rphiITSMISPDInAccP500700 = new TH1F("fHistd0rphiITSMISPDInAccP500700","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1716 fHistd0rphiITSMISPDInAccP500700->Sumw2();
1717 fHistd0rphiITSMISPDInAccP500700->SetMinimum(0);
1718 fOutput->Add(fHistd0rphiITSMISPDInAccP500700);
1719
1720 fHistd0rphiITSMISPDInAccP10001500 = new TH1F("fHistd0rphiITSMISPDInAccP10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1721 fHistd0rphiITSMISPDInAccP10001500->Sumw2();
1722 fHistd0rphiITSMISPDInAccP10001500->SetMinimum(0);
1723 fOutput->Add(fHistd0rphiITSMISPDInAccP10001500);
1724
1725 fHistd0rphiITSMISPDInAccS150200 = new TH1F("fHistd0rphiITSMISPDInAccS150200","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1726 fHistd0rphiITSMISPDInAccS150200->Sumw2();
1727 fHistd0rphiITSMISPDInAccS150200->SetMinimum(0);
1728 fOutput->Add(fHistd0rphiITSMISPDInAccS150200);
1729
1730 fHistd0rphiITSMISPDInAccS500700 = new TH1F("fHistd0rphiITSMISPDInAccS500700","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1731 fHistd0rphiITSMISPDInAccS500700->Sumw2();
1732 fHistd0rphiITSMISPDInAccS500700->SetMinimum(0);
1733 fOutput->Add(fHistd0rphiITSMISPDInAccS500700);
1734
1735 fHistd0rphiITSMISPDInAccS10001500 = new TH1F("fHistd0rphiITSMISPDInAccS10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1736 fHistd0rphiITSMISPDInAccS10001500->Sumw2();
1737 fHistd0rphiITSMISPDInAccS10001500->SetMinimum(0);
1738 fOutput->Add(fHistd0rphiITSMISPDInAccS10001500);
1739
1740 fHistd0rphiITSMIoneSPDInAccP150200 = new TH1F("fHistd0rphiITSMIoneSPDInAccP150200","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1741 fHistd0rphiITSMIoneSPDInAccP150200->Sumw2();
1742 fHistd0rphiITSMIoneSPDInAccP150200->SetMinimum(0);
1743 fOutput->Add(fHistd0rphiITSMIoneSPDInAccP150200);
1744
dbf573fb 1745 fHistd0rphiITSMIoneSPDInAccP350450 = new TH1F("fHistd0rphiITSMIoneSPDInAccP350450","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1746 fHistd0rphiITSMIoneSPDInAccP350450->Sumw2();
1747 fHistd0rphiITSMIoneSPDInAccP350450->SetMinimum(0);
1748 fOutput->Add(fHistd0rphiITSMIoneSPDInAccP350450);
1749
f8777463 1750 fHistd0rphiITSMIoneSPDInAccP500700 = new TH1F("fHistd0rphiITSMIoneSPDInAccP500700","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1751 fHistd0rphiITSMIoneSPDInAccP500700->Sumw2();
1752 fHistd0rphiITSMIoneSPDInAccP500700->SetMinimum(0);
1753 fOutput->Add(fHistd0rphiITSMIoneSPDInAccP500700);
1754
1755 fHistd0rphiITSMIoneSPDInAccP10001500 = new TH1F("fHistd0rphiITSMIoneSPDInAccP10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1756 fHistd0rphiITSMIoneSPDInAccP10001500->Sumw2();
1757 fHistd0rphiITSMIoneSPDInAccP10001500->SetMinimum(0);
1758 fOutput->Add(fHistd0rphiITSMIoneSPDInAccP10001500);
1759
dbf573fb 1760 fHistd0rphiITSMIoneSPDInAccP25004000 = new TH1F("fHistd0rphiITSMIoneSPDInAccP25004000","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1761 fHistd0rphiITSMIoneSPDInAccP25004000->Sumw2();
1762 fHistd0rphiITSMIoneSPDInAccP25004000->SetMinimum(0);
1763 fOutput->Add(fHistd0rphiITSMIoneSPDInAccP25004000);
1764
1765 fHistd0rphiITSMIoneSPDInAccP40008000 = new TH1F("fHistd0rphiITSMIoneSPDInAccP40008000","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1766 fHistd0rphiITSMIoneSPDInAccP40008000->Sumw2();
1767 fHistd0rphiITSMIoneSPDInAccP40008000->SetMinimum(0);
1768 fOutput->Add(fHistd0rphiITSMIoneSPDInAccP40008000);
1769
f8777463 1770 fHistd0zITSMIoneSPDInAccP150200 = new TH1F("fHistd0zITSMIoneSPDInAccP150200","Longitudinal imp. par. to VertexTracks for primaries; d_{0} z [cm]; N tracks",300,-1.5,1.5);
1771 fHistd0zITSMIoneSPDInAccP150200->Sumw2();
1772 fHistd0zITSMIoneSPDInAccP150200->SetMinimum(0);
1773 fOutput->Add(fHistd0zITSMIoneSPDInAccP150200);
1774
1775 fHistd0zITSMIoneSPDInAccP500700 = new TH1F("fHistd0zITSMIoneSPDInAccP500700","Longitudinal imp. par. to VertexTracks for primaries; d_{0} z [cm]; N tracks",300,-1.5,1.5);
1776 fHistd0zITSMIoneSPDInAccP500700->Sumw2();
1777 fHistd0zITSMIoneSPDInAccP500700->SetMinimum(0);
1778 fOutput->Add(fHistd0zITSMIoneSPDInAccP500700);
1779
1780 fHistd0zITSMIoneSPDInAccP10001500 = new TH1F("fHistd0zITSMIoneSPDInAccP10001500","Longitudinal imp. par. to VertexTracks for primaries; d_{0} z [cm]; N tracks",300,-1.5,1.5);
1781 fHistd0zITSMIoneSPDInAccP10001500->Sumw2();
1782 fHistd0zITSMIoneSPDInAccP10001500->SetMinimum(0);
1783 fOutput->Add(fHistd0zITSMIoneSPDInAccP10001500);
1784
9a3cd95a 1785 fHistd0zVSetaTPCInAccP10001500 = new TH2F("fHistd0zVSetaTPCInAccP10001500","Long. imp. par. to VertexSPD for TPC tracks; d_{0} z [cm]; eta",100,-4,4,10,-1,1);
1786 fOutput->Add(fHistd0zVSetaTPCInAccP10001500);
1787
f8777463 1788 fHistd0rphiVSphiITSMIoneSPDInAccP10001500 = new TH2F("fHistd0rphiVSphiITSMIoneSPDInAccP10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; phi",30,-0.3,0.3,40,0,2*3.1415);
1789 fOutput->Add(fHistd0rphiVSphiITSMIoneSPDInAccP10001500);
1790
1791 fHistd0rphiVSetaITSMIoneSPDInAccP10001500 = new TH2F("fHistd0rphiVSetaITSMIoneSPDInAccP10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; eta",30,-0.3,0.3,10,-1,1);
1792 fOutput->Add(fHistd0rphiVSetaITSMIoneSPDInAccP10001500);
1793
1794 fHistd0rphiITSMIoneSPDInAccS150200 = new TH1F("fHistd0rphiITSMIoneSPDInAccS150200","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1795 fHistd0rphiITSMIoneSPDInAccS150200->Sumw2();
1796 fHistd0rphiITSMIoneSPDInAccS150200->SetMinimum(0);
1797 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS150200);
1798
dbf573fb 1799 fHistd0rphiITSMIoneSPDInAccS350450 = new TH1F("fHistd0rphiITSMIoneSPDInAccS350450","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1800 fHistd0rphiITSMIoneSPDInAccS350450->Sumw2();
1801 fHistd0rphiITSMIoneSPDInAccS350450->SetMinimum(0);
1802 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS350450);
1803
f8777463 1804 fHistd0rphiITSMIoneSPDInAccS500700 = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1805 fHistd0rphiITSMIoneSPDInAccS500700->Sumw2();
1806 fHistd0rphiITSMIoneSPDInAccS500700->SetMinimum(0);
1807 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700);
1808
1809 fHistd0rphiITSMIoneSPDInAccS500700from22 = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700from22","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1810 fHistd0rphiITSMIoneSPDInAccS500700from22->Sumw2();
1811 fHistd0rphiITSMIoneSPDInAccS500700from22->SetMinimum(0);
1812 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700from22);
1813
1814 fHistd0rphiITSMIoneSPDInAccS500700from211 = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700from211","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1815 fHistd0rphiITSMIoneSPDInAccS500700from211->Sumw2();
1816 fHistd0rphiITSMIoneSPDInAccS500700from211->SetMinimum(0);
1817 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700from211);
1818
1819 fHistd0rphiITSMIoneSPDInAccS500700from310 = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700from310","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1820 fHistd0rphiITSMIoneSPDInAccS500700from310->Sumw2();
1821 fHistd0rphiITSMIoneSPDInAccS500700from310->SetMinimum(0);
1822 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700from310);
1823
1824 fHistd0rphiITSMIoneSPDInAccS500700from321 = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700from321","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1825 fHistd0rphiITSMIoneSPDInAccS500700from321->Sumw2();
1826 fHistd0rphiITSMIoneSPDInAccS500700from321->SetMinimum(0);
1827 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700from321);
1828
dbf573fb 1829 fHistd0rphiITSMIoneSPDInAccS500700from3122 = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700from3122","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1830 fHistd0rphiITSMIoneSPDInAccS500700from3122->Sumw2();
1831 fHistd0rphiITSMIoneSPDInAccS500700from3122->SetMinimum(0);
1832 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700from3122);
1833
f8777463 1834 fHistd0rphiITSMIoneSPDInAccS10001500 = new TH1F("fHistd0rphiITSMIoneSPDInAccS10001500","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1835 fHistd0rphiITSMIoneSPDInAccS10001500->Sumw2();
1836 fHistd0rphiITSMIoneSPDInAccS10001500->SetMinimum(0);
1837 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS10001500);
dbf573fb 1838
1839 fHistd0rphiITSMIoneSPDInAccS25004000 = new TH1F("fHistd0rphiITSMIoneSPDInAccS25004000","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1840 fHistd0rphiITSMIoneSPDInAccS25004000->Sumw2();
1841 fHistd0rphiITSMIoneSPDInAccS25004000->SetMinimum(0);
1842 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS25004000);
1843
1844 fHistd0rphiITSMIoneSPDInAccS40008000 = new TH1F("fHistd0rphiITSMIoneSPDInAccS40008000","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1845 fHistd0rphiITSMIoneSPDInAccS40008000->Sumw2();
1846 fHistd0rphiITSMIoneSPDInAccS40008000->SetMinimum(0);
1847 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS40008000);
1848
1849 fHistd0rphiITSMIoneSPDInAccS150200fromStrange = new TH1F("fHistd0rphiITSMIoneSPDInAccS150200fromStrange","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1850 fHistd0rphiITSMIoneSPDInAccS150200fromStrange->Sumw2();
1851 fHistd0rphiITSMIoneSPDInAccS150200fromStrange->SetMinimum(0);
1852 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS150200fromStrange);
1853
1854 fHistd0rphiITSMIoneSPDInAccS150200fromMat = new TH1F("fHistd0rphiITSMIoneSPDInAccS150200fromMat","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1855 fHistd0rphiITSMIoneSPDInAccS150200fromMat->Sumw2();
1856 fHistd0rphiITSMIoneSPDInAccS150200fromMat->SetMinimum(0);
1857 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS150200fromMat);
1858
1859 fHistd0rphiITSMIoneSPDInAccS350450fromStrange = new TH1F("fHistd0rphiITSMIoneSPDInAccS350450fromStrange","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1860 fHistd0rphiITSMIoneSPDInAccS350450fromStrange->Sumw2();
1861 fHistd0rphiITSMIoneSPDInAccS350450fromStrange->SetMinimum(0);
1862 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS350450fromStrange);
1863
1864 fHistd0rphiITSMIoneSPDInAccS350450fromMat = new TH1F("fHistd0rphiITSMIoneSPDInAccS350450fromMat","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1865 fHistd0rphiITSMIoneSPDInAccS350450fromMat->Sumw2();
1866 fHistd0rphiITSMIoneSPDInAccS350450fromMat->SetMinimum(0);
1867 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS350450fromMat);
1868
1869 fHistd0rphiITSMIoneSPDInAccS500700fromStrange = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700fromStrange","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1870 fHistd0rphiITSMIoneSPDInAccS500700fromStrange->Sumw2();
1871 fHistd0rphiITSMIoneSPDInAccS500700fromStrange->SetMinimum(0);
1872 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700fromStrange);
1873
1874 fHistd0rphiITSMIoneSPDInAccS500700fromMat = new TH1F("fHistd0rphiITSMIoneSPDInAccS500700fromMat","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1875 fHistd0rphiITSMIoneSPDInAccS500700fromMat->Sumw2();
1876 fHistd0rphiITSMIoneSPDInAccS500700fromMat->SetMinimum(0);
1877 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS500700fromMat);
1878
1879 fHistd0rphiITSMIoneSPDInAccS10001500fromStrange = new TH1F("fHistd0rphiITSMIoneSPDInAccS10001500fromStrange","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1880 fHistd0rphiITSMIoneSPDInAccS10001500fromStrange->Sumw2();
1881 fHistd0rphiITSMIoneSPDInAccS10001500fromStrange->SetMinimum(0);
1882 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS10001500fromStrange);
1883
1884 fHistd0rphiITSMIoneSPDInAccS10001500fromMat = new TH1F("fHistd0rphiITSMIoneSPDInAccS10001500fromMat","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1885 fHistd0rphiITSMIoneSPDInAccS10001500fromMat->Sumw2();
1886 fHistd0rphiITSMIoneSPDInAccS10001500fromMat->SetMinimum(0);
1887 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS10001500fromMat);
1888
1889 fHistd0rphiITSMIoneSPDInAccS25004000fromStrange = new TH1F("fHistd0rphiITSMIoneSPDInAccS25004000fromStrange","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1890 fHistd0rphiITSMIoneSPDInAccS25004000fromStrange->Sumw2();
1891 fHistd0rphiITSMIoneSPDInAccS25004000fromStrange->SetMinimum(0);
1892 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS25004000fromStrange);
1893
1894 fHistd0rphiITSMIoneSPDInAccS25004000fromMat = new TH1F("fHistd0rphiITSMIoneSPDInAccS25004000fromMat","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1895 fHistd0rphiITSMIoneSPDInAccS25004000fromMat->Sumw2();
1896 fHistd0rphiITSMIoneSPDInAccS25004000fromMat->SetMinimum(0);
1897 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS25004000fromMat);
1898
1899 fHistd0rphiITSMIoneSPDInAccS40008000fromStrange = new TH1F("fHistd0rphiITSMIoneSPDInAccS40008000fromStrange","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1900 fHistd0rphiITSMIoneSPDInAccS40008000fromStrange->Sumw2();
1901 fHistd0rphiITSMIoneSPDInAccS40008000fromStrange->SetMinimum(0);
1902 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS40008000fromStrange);
1903
1904 fHistd0rphiITSMIoneSPDInAccS40008000fromMat = new TH1F("fHistd0rphiITSMIoneSPDInAccS40008000fromMat","Transverse imp. par. to VertexTracks for primaries; d_{0} rphi [cm]; N tracks",300,-1.5,1.5);
1905 fHistd0rphiITSMIoneSPDInAccS40008000fromMat->Sumw2();
1906 fHistd0rphiITSMIoneSPDInAccS40008000fromMat->SetMinimum(0);
1907 fOutput->Add(fHistd0rphiITSMIoneSPDInAccS40008000fromMat);
1908
f8777463 1909
1910 fHistd0zITSMIoneSPDInAccS150200 = new TH1F("fHistd0zITSMIoneSPDInAccS150200","Longitudinal imp. par. to VertexTracks for secondaries; d_{0} z [cm]; N tracks",300,-1.5,1.5);
1911 fHistd0zITSMIoneSPDInAccS150200->Sumw2();
1912 fHistd0zITSMIoneSPDInAccS150200->SetMinimum(0);
1913 fOutput->Add(fHistd0zITSMIoneSPDInAccS150200);
1914
1915 fHistd0zITSMIoneSPDInAccS500700 = new TH1F("fHistd0zITSMIoneSPDInAccS500700","Longitudinal imp. par. to VertexTracks for secondaries; d_{0} z [cm]; N tracks",300,-1.5,1.5);
1916 fHistd0zITSMIoneSPDInAccS500700->Sumw2();
1917 fHistd0zITSMIoneSPDInAccS500700->SetMinimum(0);
1918 fOutput->Add(fHistd0zITSMIoneSPDInAccS500700);
1919
1920 fHistd0zITSMIoneSPDInAccS10001500 = new TH1F("fHistd0zITSMIoneSPDInAccS10001500","Longitudinal imp. par. to VertexTracks for secondaries; d_{0} z [cm]; N tracks",300,-1.5,1.5);
1921 fHistd0zITSMIoneSPDInAccS10001500->Sumw2();
1922 fHistd0zITSMIoneSPDInAccS10001500->SetMinimum(0);
1923 fOutput->Add(fHistd0zITSMIoneSPDInAccS10001500);
1924
af874c53 1925 fHistITSRedChi2NonFakePt02 = new TH1F("fHistITSRedChi2NonFakePt02","ITS chi2/cluster; ITS #chi^{2}/nclusters; tracks",500,0,100);
1926 fOutput->Add(fHistITSRedChi2NonFakePt02);
1927 fHistITSRedChi2FakePt02 = new TH1F("fHistITSRedChi2FakePt02","ITS chi2/cluster; ITS #chi^{2}/nclusters; tracks",500,0,100);
1928 fOutput->Add(fHistITSRedChi2FakePt02);
1929 fHistITSRedChi2NonFakePt05 = new TH1F("fHistITSRedChi2NonFakePt05","ITS chi2/cluster; ITS #chi^{2}/nclusters; tracks",500,0,100);
1930 fOutput->Add(fHistITSRedChi2NonFakePt05);
1931 fHistITSRedChi2FakePt05 = new TH1F("fHistITSRedChi2FakePt05","ITS chi2/cluster; ITS #chi^{2}/nclusters; tracks",500,0,100);
1932 fOutput->Add(fHistITSRedChi2FakePt05);
1933 fHistITSRedChi2NonFakePt1 = new TH1F("fHistITSRedChi2NonFakePt1","ITS chi2/cluster; ITS #chi^{2}/nclusters; tracks",500,0,100);
1934 fOutput->Add(fHistITSRedChi2NonFakePt1);
1935 fHistITSRedChi2FakePt1 = new TH1F("fHistITSRedChi2FakePt1","ITS chi2/cluster; ITS #chi^{2}/nclusters; tracks",500,0,100);
1936 fOutput->Add(fHistITSRedChi2FakePt1);
1937
f8777463 1938
8d63376d 1939 // ntuples
1940 //
dbf573fb 1941 fNtupleESDTracks = new TNtuple("fNtupleESDTracks","tracks","pt:eta:phi:d0:z0:sigmad0:sigmaz0:ptMC:pdgMC:pdgMCMoth:d0MC:d0MCv:z0MCv:sigmad0MCv:sigmaz0MCv:ITSflag:isPrimary:isTPCSel");
8d63376d 1942 fOutput->Add(fNtupleESDTracks);
1943
c2893e82 1944 fNtupleITSAlignExtra = new TNtuple("fNtupleITSAlignExtra","ITS alignment checks: extra clusters","layer:x:y:z:dxy:dz:xloc:zloc:npoints:pt");
8d63376d 1945 fOutput->Add(fNtupleITSAlignExtra);
1946
c2893e82 1947 fNtupleITSAlignSPDTracklets = new TNtuple("fNtupleITSAlignSPDTracklets","ITS alignment checks: SPD tracklets wrt SPD vertex","phi:theta:z:dxy:dz:pt");
8d63376d 1948 fOutput->Add(fNtupleITSAlignSPDTracklets);
1949
4f7217d1 1950 PostData(1, fOutput);
8d63376d 1951 return;
1952}
1953
1954//________________________________________________________________________
d98a315a 1955void AliAnalysisTaskITSTrackingCheck::UserExec(Option_t *)
8d63376d 1956{
1957 // Main loop
1958 // Called for each event
c0d39d2a 1959
31a96e36 1960
d98a315a 1961 fESD = dynamic_cast<AliESDEvent*>(InputEvent());
1962
1963 if (!fESD) {
1964 Printf("ERROR: fESD not available");
1965 return;
1966 }
1967
0b2e8130 1968 //if(fESD->GetEventNumberInFile()<7) return;
1969
5ac4a99f 1970 // only events in the requested multiplicity range
1971 if(!IsSelectedCentrality()) return;
f8777463 1972
1973 fHistNEvents->Fill(-1);
1974
c0d39d2a 1975
f8777463 1976 Bool_t isSelected = kTRUE;
1977 if(fUsePhysSel) {
1978 isSelected = ((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected();
1979 }
1980 if(!isSelected) return;
1981
31a96e36 1982 TString firedTriggerClasses=fESD->GetFiredTriggerClasses();
1983 if(!firedTriggerClasses.Contains(fTriggerClass.Data())) return;
d98a315a 1984
31a96e36 1985
1986 Bool_t sddIsIn=kTRUE;
1987 if(fCheckSDDIsIn) {
1988
1989 if(!fTrigConfig) {
1990 AliCDBManager* man = AliCDBManager::Instance();
1991 if(fOCDBPath.Contains("OCDB")) { // when running in the QAtrain this is not called (OCBD is already set)
1992 man->SetDefaultStorage(fOCDBPath.Data());
1993 man->SetRun(fESD->GetRunNumber());
1994 }
1995 if(!man) {
1996 AliFatal("CDB not set but needed by AliAnalysisTaskITSTrackingCheck");
1997 return;
1998 }
1999 AliCDBEntry* eT=(AliCDBEntry*)man->Get("GRP/CTP/Config");
2000 if(eT) {
2001 fTrigConfig=(AliTriggerConfiguration*)eT->GetObject();
2002 }
2003 if(!eT || !fTrigConfig) {
2004 AliError("Cannot retrieve CDB entry for GRP/CTP/Config");
2005 return;
2006 }
2007 }
2008
2009 sddIsIn=fESD->IsDetectorInTriggerCluster("ITSSDD",fTrigConfig);
2010 }
2011
2012
2308c814 2013 //if(fESD->GetEventType()!=7) return;
8d63376d 2014
8d63376d 2015 // *********** MC info ***************
2016 TArrayF mcVertex(3);
2017 mcVertex[0]=9999.; mcVertex[1]=9999.; mcVertex[2]=9999.;
2018 Float_t dNchdy=-999.;
2019
2020 TParticle *part=0;
2021 AliESDVertex *vertexMC=0;
2022 AliStack *stack=0;
2023 if (fReadMC) {
2024 AliMCEventHandler *eventHandler = dynamic_cast<AliMCEventHandler*> (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler());
2025 if (!eventHandler) {
2026 Printf("ERROR: Could not retrieve MC event handler");
2027 return;
2028 }
2029
2030 AliMCEvent* mcEvent = eventHandler->MCEvent();
2031 if (!mcEvent) {
2032 Printf("ERROR: Could not retrieve MC event");
2033 return;
2034 }
2035
2036 stack = mcEvent->Stack();
2037 if (!stack) {
2038 AliDebug(AliLog::kError, "Stack not available");
2039 return;
2040 }
2041
2042 AliHeader* header = mcEvent->Header();
2043 if (!header) {
2044 AliDebug(AliLog::kError, "Header not available");
2045 return;
2046 }
2047 AliGenEventHeader* genHeader = header->GenEventHeader();
2048 genHeader->PrimaryVertex(mcVertex);
d1d5a660 2049 //printf("MC vertex: %f %f %f\n",mcVertex[0],mcVertex[1],mcVertex[2]);
8d63376d 2050
2051 Int_t ngenpart = (Int_t)stack->GetNtrack();
d1d5a660 2052 //printf("# generated particles = %d\n",ngenpart);
8d63376d 2053 dNchdy=0;
2054 for(Int_t ip=0; ip<ngenpart; ip++) {
2055 part = (TParticle*)stack->Particle(ip);
2056 // keep only electrons, muons, pions, kaons and protons
2057 Int_t apdg = TMath::Abs(part->GetPdgCode());
2058 if(apdg!=11 && apdg!=13 && apdg!=211 && apdg!=321 && apdg!=2212) continue;
2059 // reject secondaries
af874c53 2060 //if(TMath::Sqrt((part->Vx()-mcVertex[0])*(part->Vx()-mcVertex[0])+(part->Vy()-mcVertex[1])*(part->Vy()-mcVertex[1]))<.01) printf("%f\n",TMath::Sqrt((part->Vx()-mcVertex[0])*(part->Vx()-mcVertex[0])+(part->Vy()-mcVertex[1])*(part->Vy()-mcVertex[1])));
8d63376d 2061 if(TMath::Sqrt((part->Vx()-mcVertex[0])*(part->Vx()-mcVertex[0])+(part->Vy()-mcVertex[1])*(part->Vy()-mcVertex[1]))>0.0010) continue;
2062 // reject incoming protons
2063 Double_t energy = part->Energy();
2064 if(energy>900.) continue;
2065 Double_t pz = part->Pz();
26f8147c 2066 Double_t y = 100.;
2067 if(energy-pz+1.e-13 > 0.) y = 0.5*TMath::Log((energy+pz+1.e-13)/(energy-pz+1.e-13));
8d63376d 2068 if(TMath::Abs(y)<1.0) dNchdy += 0.5; // count 1/2 of particles in |y|<1
2069 }
d1d5a660 2070 //printf("# primary particles = %7.1f\n",dNchdy);
8d63376d 2071 }
2072 // *********** MC info ***************
2073 Double_t mcVtxPos[3]={mcVertex[0],mcVertex[1],mcVertex[2]},mcVtxSigma[3]={0,0,0};
2074 vertexMC = new AliESDVertex(mcVtxPos,mcVtxSigma);
2075
8d63376d 2076
f8777463 2077
2078 //------- event selection --------
2079 Int_t mincontrSPDvtx=1;
0b2e8130 2080 Double_t maxzSPDvtx=5.;
ebcbf764 2081 Double_t maxrSPDvtx=0.5;
dbf573fb 2082 Bool_t useTRKvtx=kTRUE;
2083 Bool_t useSPDvtxifNotTRK=kTRUE;
2084 Int_t mincontrTRKvtx=1;
2085 //Double_t maxzTRKvtx=20.;
2086 //Double_t maxrTRKvtx=1.;
f8777463 2087 //------- TPC track selection --------
0b2e8130 2088 Int_t minclsTPC=70;
f8777463 2089 Double_t maxchi2perTPCcl=4.;
2090 Double_t minEtaInAcc=-0.8; // -0.8
2091 Double_t maxEtaInAcc=0.8; // 0.8
9a3cd95a 2092 Double_t maxdcaxy=2.4;//2.4;
2093 Double_t maxdcaz=3.2;//3.2;
f8777463 2094 AliESDtrackCuts* esdtrackCutsTPC = new AliESDtrackCuts("esdtrackCutsTPC");
2095 esdtrackCutsTPC->SetMaxDCAToVertexXY(maxdcaxy);
2096 esdtrackCutsTPC->SetMaxDCAToVertexZ(maxdcaz);
2097 esdtrackCutsTPC->SetDCAToVertex2D(kTRUE);
2098 esdtrackCutsTPC->SetRequireSigmaToVertex(kFALSE);
2099 esdtrackCutsTPC->SetRequireTPCRefit(kFALSE);// cannot do it because status not copied in AliESDtrack::FillTPCOnlyTrack
2100 esdtrackCutsTPC->SetAcceptKinkDaughters(kFALSE);
2101 esdtrackCutsTPC->SetMinNClustersTPC(minclsTPC);
2102 esdtrackCutsTPC->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
2103 esdtrackCutsTPC->SetEtaRange(minEtaInAcc,maxEtaInAcc);
2104 SetESDtrackCutsTPC(esdtrackCutsTPC);
2105 //------- ITS+TPC track selection --------
2106 Double_t maxdcaxyITSTPC=0.2;
2107 Double_t maxdcazITSTPC=1.e6;
2108 AliESDtrackCuts* esdtrackCutsITSTPC = new AliESDtrackCuts("esdtrackCutsITSTPC");
2109 esdtrackCutsITSTPC->SetMaxDCAToVertexXY(maxdcaxyITSTPC);
2110 esdtrackCutsITSTPC->SetMaxDCAToVertexZ(maxdcazITSTPC);
2111 esdtrackCutsITSTPC->SetDCAToVertex2D(kFALSE);
2112 esdtrackCutsITSTPC->SetRequireSigmaToVertex(kFALSE);
2113 esdtrackCutsITSTPC->SetRequireITSRefit(kTRUE);
2114 esdtrackCutsITSTPC->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
2115 AliESDtrackCuts::kAny);
2116 esdtrackCutsITSTPC->SetAcceptKinkDaughters(kFALSE);
2117 esdtrackCutsITSTPC->SetMinNClustersTPC(minclsTPC);
2118 esdtrackCutsITSTPC->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
2119 esdtrackCutsITSTPC->SetEtaRange(minEtaInAcc,maxEtaInAcc);
2120 SetESDtrackCutsITSTPC(esdtrackCutsITSTPC);
2121 //---------------------------------------
2122
2123
8d63376d 2124 //
2308c814 2125
2126 /*
8d63376d 2127 // ********** Trigger *****************
2128 ULong64_t triggerMask;
2129 ULong64_t spdFO = (1 << 14);
2130 ULong64_t v0left = (1 << 11);
2131 ULong64_t v0right = (1 << 12);
2132
2133 triggerMask=fESD->GetTriggerMask();
2134 // MB1: SPDFO || V0L || V0R
2135 Bool_t eventTriggered = (triggerMask & spdFO || ((triggerMask & v0left) || (triggerMask & v0right)));
2136 //MB2: GFO && V0R
2137 //triggerMask & spdFO && ((triggerMask&v0left) || (triggerMask&v0right))
2138 // ************ Trigger ******************
2139 if(!eventTriggered) return;
2308c814 2140 */
8d63376d 2141
f8777463 2142
2143 fHistNEvents->Fill(0);
2144
8d63376d 2145 // SPD vertex
2146 const AliESDVertex *spdv=fESD->GetPrimaryVertexSPD();
f8777463 2147 // Select good SPD vertices
2148 TString spdvtitle=spdv->GetTitle();
2149 //if(!spdvtitle.Contains("3D")) return;
b7fb742d 2150 if(fRequireSPDvtx) {
2151 if(spdv->GetNContributors()<mincontrSPDvtx ||
e690d4d0 2152 TMath::Abs(spdv->GetZ())>maxzSPDvtx ||
2153 spdv->GetX()*spdv->GetX()+spdv->GetY()*spdv->GetY()>maxrSPDvtx) {
b7fb742d 2154 delete esdtrackCutsTPC; esdtrackCutsTPC=0;
2155 delete esdtrackCutsITSTPC; esdtrackCutsITSTPC=0;
2156 return;
2157 }
f8777463 2158 }
d1d5a660 2159 //spdv->Print();
dbf573fb 2160
9f024d7a 2161 if(fRejPileupSPD) if(fESD->IsPileupFromSPD()) return;
2162
f8777463 2163 //
2164 // Tracks vertex
2165 const AliESDVertex *vertexESD = fESD->GetPrimaryVertexTracks();
d1d5a660 2166 //vertexESD->Print();
f8777463 2167
2168 if(spdvtitle.Contains("3D")) {
2169 fHistNEvents->Fill(1);
2170 fHistNEvents->Fill(3);
2171 } else {
2172 fHistNEvents->Fill(2);
2173 fHistNEvents->Fill(3);
2174 }
2175 if(vertexESD) {
2176 if(vertexESD->GetStatus()) fHistNEvents->Fill(4);
2177 }
2178
8d63376d 2179 Int_t ntracks = fESD->GetNumberOfTracks();
abe2795e 2180 //printf("Tracks # = %d\n",fESD->GetNumberOfTracks());
8d63376d 2181
2182 fHistNtracks->Fill(ntracks);
2183 // Post the data already here
d98a315a 2184 PostData(1, fOutput);
8d63376d 2185
0b2e8130 2186 FillNClustersModuleMap();
2187
8d63376d 2188 Int_t idet,status; Float_t xloc,zloc;
3f18d0a6 2189 Double_t rSPDouter=7.6,rSDDouter=23.9,rSSDouter=43.1;
2190 Double_t zSPDouter=14.1,zSDDouter=29.7,zSSDouter=48.9;
2191
b7fb742d 2192
2193 // loop on tracks to check if there are 7 ITS layers (upgrade)
2194 for(Int_t itr=0; itr<ntracks; itr++) {
2195 AliESDtrack *track = fESD->GetTrack(itr);
2196 if(track->HasPointOnITSLayer(6)) {
2197 fNITSLayers=7;
2198 break;
2199 }
2200 }
2201
8d63376d 2202 // loop on tracks
2203 for(Int_t itr=0; itr<ntracks; itr++) {
2204 AliESDtrack *track = fESD->GetTrack(itr);
af874c53 2205 //if(track->Charge()>0) continue;
2308c814 2206 // remove kink daughters
8d63376d 2207 if(track->GetKinkIndex(0)>0) continue;
2208
c94b80d3 2209
3f18d0a6 2210 // remove tracks not reco in ITS or TPC
2211 if (!(track->GetStatus() & AliESDtrack::kITSin) &&
2212 !(track->GetStatus() & AliESDtrack::kTPCin)) continue;
8d63376d 2213
c94b80d3 2214 // check if tracks if flagged as from bc0 by TOF
2215 Bool_t isTOFbc0 = kFALSE;
2216 if(track->GetTOFBunchCrossing(fESD->GetMagneticField())/2==0) isTOFbc0=kTRUE;
2217
2218 //
f8777463 2219 Bool_t isPrimary=kTRUE,isFromMat=kFALSE,isFromStrange=kFALSE;
2220 Double_t rProdVtx=0,zProdVtx=0;
2221 Int_t pdgTrk=0,pdgMoth=0;
dbf573fb 2222 Int_t nClsMCSPD=0;
f8777463 2223
2224 Int_t trkLabel = TMath::Abs(track->GetLabel());
0faeb7ba 2225 Bool_t isFake=(track->GetLabel()>=0 ? kFALSE : kTRUE);
0b2e8130 2226 Bool_t hasShared=kFALSE,hasSharedSPD=kFALSE;
b7fb742d 2227 for(Int_t i=0;i<fNITSLayers;i++) {
0b2e8130 2228 if(track->HasSharedPointOnITSLayer(i)) {
2229 hasShared=kTRUE;
b7fb742d 2230 if(i<(fNITSLayers-4)) hasSharedSPD=kTRUE;
0b2e8130 2231 }
2232 }
2233
f8777463 2234 // check if it is primary
2235 if(fReadMC && stack) {
2236 isPrimary = stack->IsPhysicalPrimary(trkLabel);
2237 part = (TParticle*)stack->Particle(trkLabel);
2238 rProdVtx = TMath::Sqrt((part->Vx()-mcVertex[0])*(part->Vx()-mcVertex[0])+(part->Vy()-mcVertex[1])*(part->Vy()-mcVertex[1]));
2239 zProdVtx = TMath::Abs(part->Vz()-mcVertex[2]);
2240 //if(rProdVtx<2.8) isPrimary=kTRUE; // this could be tried
2241 pdgTrk = TMath::Abs(part->GetPdgCode());
2242 if(part->GetFirstMother()>=0) {
2243 TParticle* mm=stack->Particle(part->GetFirstMother());
2244 if(mm) pdgMoth = TMath::Abs(mm->GetPdgCode());
2245 }
dbf573fb 2246 if(pdgMoth==310 || pdgMoth==321 || pdgMoth==3122 || pdgMoth==3312) isFromStrange=kTRUE;
f8777463 2247 if(pdgMoth==211 || pdgMoth==22 || pdgMoth==2112 || pdgMoth==2212) isFromMat=kTRUE;
dbf573fb 2248 nClsMCSPD = NumberOfITSClustersMC(trkLabel,240);
f8777463 2249 }
2250
2251 Bool_t itsrefit=kFALSE,tpcrefit=kFALSE,itsfindable=kFALSE,itsfindableAcc=kFALSE;
8d63376d 2252 if ((track->GetStatus() & AliESDtrack::kITSrefit)) itsrefit=kTRUE;
f8777463 2253 if ((track->GetStatus() & AliESDtrack::kTPCrefit)) tpcrefit=kTRUE;
2254 //if ((track->GetStatus() & AliESDtrack::kTPCin)) tpcrefit=kTRUE;
8d63376d 2255
0b2e8130 2256
2257
f8777463 2258 // remove tracks with kTPCin and not kTPCrefit
2259 if ((track->GetStatus() & AliESDtrack::kTPCin) && !tpcrefit) continue;
8d63376d 2260
f8777463 2261 AliESDtrack *trackTPC = 0;
2262 if(tpcrefit) trackTPC = AliESDtrackCuts::GetTPCOnlyTrack(fESD,itr);
674ba68b 2263 // if(trackTPC) trackTPC->RelateToVertex(spdv,fESD->GetMagneticField(),100.); // relate it to the SPD vertex
f8777463 2264
2308c814 2265
2308c814 2266 Int_t nclsITS = track->GetNcls(0);
6c6fab7b 2267 Int_t nclsSDDSSD = 0;
3f18d0a6 2268 Int_t nclsokbadoutinzITS = 0;
2308c814 2269 Bool_t outInZ=kFALSE;
8d63376d 2270
b7fb742d 2271 for(Int_t layer=0; layer<fNITSLayers; layer++) {
2272 if(layer>=(fNITSLayers-4) && track->HasPointOnITSLayer(layer)) nclsSDDSSD++;
dbf573fb 2273 if(layer==0 && !track->HasPointOnITSLayer(1)) continue;
2274 if(layer==1 && !track->HasPointOnITSLayer(0)) continue;
8d63376d 2275 track->GetITSModuleIndexInfo(layer,idet,status,xloc,zloc);
f8777463 2276 if(status<0) continue;
2308c814 2277 if(layer>=2) idet+=240; // add n SPD modules
f8777463 2278 if(layer>=4) idet+=260; // add n SDD modules
2308c814 2279 if(status==4) outInZ=kTRUE;
dbf573fb 2280
2281 if(layer==1 && tpcrefit && itsrefit) {
2282 Double_t xyz[3];
2283 track->GetXYZAt(rSPDouter,fESD->GetMagneticField(),xyz);
2284 if(status==1) {
2285 if(idet%4==0) fHistZatSPDouter0ok->Fill(xyz[2]);
2286 if(idet%4==1) fHistZatSPDouter1ok->Fill(xyz[2]);
2287 if(idet%4==2) fHistZatSPDouter2ok->Fill(xyz[2]);
2288 if(idet%4==3) fHistZatSPDouter3ok->Fill(xyz[2]);
2289 }
2290 if(status==2 || status==5) {
2291 if(idet%4==0) fHistZatSPDouter0notok->Fill(xyz[2]);
2292 if(idet%4==1) fHistZatSPDouter1notok->Fill(xyz[2]);
2293 if(idet%4==2) fHistZatSPDouter2notok->Fill(xyz[2]);
2294 if(idet%4==3) fHistZatSPDouter3notok->Fill(xyz[2]);
2295 }
2296 }
2297 // kill SDD clusters in MC
2298 if(fReadMC && (layer==2 || layer==3)) {
2299 Double_t qITS[4]; track->GetITSdEdxSamples(qITS);
2300 Double_t conv=3.34/5.43;
2301 if(layer==2 && status==1 && qITS[0]*conv < 30.) status=5;
2302 if(layer==3 && status==1 && qITS[1]*conv < 30.) status=5;
2303 }
f8777463 2304 if(tpcrefit) {
dbf573fb 2305 if(trackTPC) if(!fESDtrackCutsTPC->AcceptTrack(trackTPC)) continue;
0b2e8130 2306 if(track->Pt()>0.7 && track->Pt()<2.0) {
2307 if(status==1) {
2308 fHistClusterMapITSMIok->Fill(layer);
0b2e8130 2309 }
2310 if(status==2) fHistClusterMapITSMIbad->Fill(layer);
2311 if(status==3) fHistClusterMapITSMIskipped->Fill(layer);
2312 if(status==4) fHistClusterMapITSMIoutinz->Fill(layer);
2313 if(status==5) fHistClusterMapITSMInocls->Fill(layer);
2314 if(status==6) fHistClusterMapITSMInorefit->Fill(layer);
2315 if(status==1 && !outInZ) fHistClusterMapModuleITSMIokInAcc->Fill(idet);
2316 if((status==2 || status==7) && !outInZ) fHistClusterMapModuleITSMIbadInAcc->Fill(idet);
2317 if(status==5 && !outInZ) fHistClusterMapModuleITSMInoclsInAcc->Fill(idet);
2318
0b2e8130 2319
2320 }
3f18d0a6 2321 if(status==1 || status==2 || status==4) {
2322 fHistClusterMapITSMIokoutinzbad->Fill(layer);
2323 nclsokbadoutinzITS++;
2324 }
dbf573fb 2325 /*if(idet>=80 && idet<120) {
2326 Double_t xyz[3];
2327 track->GetXYZAt(6.8,fESD->GetMagneticField(),xyz);
2328 printf("det %d status %d phi %f z %f zloc %f xloc %f\n",idet,status,TMath::ATan2(xyz[1],xyz[0]),xyz[2],zloc,xloc);
2329 }*/
2330 if((layer==2 || layer==3) && status!=2 && status!=4) {
2331 //Float_t xlocCls;
2332 //Int_t nClsInMod = NumberOfITSClusters(idet,xlocCls);
f8777463 2333 fHistxlocSDDall->Fill(xloc);
2334 fHistzlocSDDall->Fill(zloc);
dbf573fb 2335 if(status==1 /*nClsInMod>0*/) {
f8777463 2336 fHistxlocSDDok->Fill(xloc);
2337 fHistxlocVSmodSDDok->Fill(idet,xloc);
2338 fHistzlocSDDok->Fill(zloc);
0b2e8130 2339
dbf573fb 2340 }
f8777463 2341 }
2308c814 2342 }
2343 if(TESTBIT(track->GetITSClusterMap(),layer)) {
f8777463 2344 if(tpcrefit) {
2308c814 2345 fHistClusterMapITSMI->Fill(layer);
2346 } else {
2347 fHistClusterMapITSSA->Fill(layer);
2348 if(!outInZ) fHistClusterMapITSSAInAcc->Fill(layer);
2349 }
8d63376d 2350 }
2351 }
2352
3f18d0a6 2353 // TPC track findable in ITS
f8777463 2354 if(tpcrefit && trackTPC) {
2355 if(fESDtrackCutsTPC->AcceptTrack(trackTPC)) {
2356 itsfindable=kTRUE;
2357 Double_t zAtSSDouter=100,zAtSDDouter=100,zAtSPDouter=100;
2358 track->GetZAt(rSSDouter,fESD->GetMagneticField(),zAtSSDouter);
2359 track->GetZAt(rSDDouter,fESD->GetMagneticField(),zAtSDDouter);
2360 track->GetZAt(rSPDouter,fESD->GetMagneticField(),zAtSPDouter);
2361 fHistPtTPC->Fill(track->Pt());
dbf573fb 2362 if(TMath::Abs(zAtSSDouter)<1.*zSSDouter &&
2363 TMath::Abs(zAtSDDouter)<1.*zSDDouter &&
2364 TMath::Abs(zAtSPDouter)<1.*zSPDouter) {
0b2e8130 2365 itsfindableAcc=kTRUE;
f8777463 2366 fHistdEdxVSPtTPCInAcc->Fill(track->Pt(),track->GetTPCsignal());
2367 fHistPtTPCInAcc->Fill(track->Pt());
31a96e36 2368 if(fCheckSDDIsIn && sddIsIn) fHistPtTPCInAccwSDD->Fill(track->Pt());
2369 if(fCheckSDDIsIn && !sddIsIn) fHistPtTPCInAccwoSDD->Fill(track->Pt());
2370 if(isTOFbc0) {
2371 fHistPtTPCInAccTOFbc0->Fill(track->Pt());
2372 if(fCheckSDDIsIn && sddIsIn) fHistPtTPCInAccTOFbc0wSDD->Fill(track->Pt());
2373 if(fCheckSDDIsIn && !sddIsIn) fHistPtTPCInAccTOFbc0woSDD->Fill(track->Pt());
2374 }
0b2e8130 2375 fHistTPCclsVSPtTPCInAcc->Fill(trackTPC->Pt(),trackTPC->GetNcls(1));
f8777463 2376 fHistPtVSphiTPCInAcc->Fill(track->Phi(),track->Pt());
2377 if(!(track->GetStatus()&AliESDtrack::kTRDout)) fHistPtTPCInAccNoTRDout->Fill(track->Pt());
2378 if(!(track->GetStatus()&AliESDtrack::kTOFout)) fHistPtTPCInAccNoTOFout->Fill(track->Pt());
2379 fHistPtTPCInAccWithPtTPCAtVtx->Fill(trackTPC->Pt());
2380 Double_t pTPCinnerwall[3];
2381 track->GetInnerPxPyPz(pTPCinnerwall);
2382 Double_t ptTPCinnerwall=TMath::Sqrt(pTPCinnerwall[0]*pTPCinnerwall[0]+pTPCinnerwall[1]*pTPCinnerwall[1]);
2383 fHistPtTPCInAccWithPtTPCAtInnerWall->Fill(ptTPCinnerwall);
dbf573fb 2384 if(track->GetStatus()&AliESDtrack::kITSrefit && !(track->GetStatus()&AliESDtrack::kTOFout)) fHistDeltaPtTPC->Fill(trackTPC->P(),track->P()-trackTPC->P());
547640ff 2385 if(track->Pt()>1) {
2386 fHistPhiTPCInAcc->Fill(track->Phi());
2387 fHistEtaTPCInAcc->Fill(track->Eta());
2388 }
f8777463 2389 if(isPrimary) {
2390 fHistPtTPCInAccP->Fill(track->Pt());
2391 if(pdgTrk==321) fHistPtTPCInAccPfromStrange->Fill(track->Pt());
2392 } else {
2393 fHistPtTPCInAccS->Fill(track->Pt());
2394 if(isFromMat) fHistPtTPCInAccSfromMat->Fill(track->Pt());
2395 if(isFromStrange) fHistPtTPCInAccSfromStrange->Fill(track->Pt());
2396 }
dbf573fb 2397 if(nClsMCSPD==2) fHistPtTPCInAccMCtwoSPD->Fill(track->Pt());
2398 if(nClsMCSPD>=1) fHistPtTPCInAccMConeSPD->Fill(track->Pt());
0b2e8130 2399
f8777463 2400 //if(isPrimary) {fHistRProdVtxInAccP->Fill(rProdVtx);} else {fHistRProdVtxInAccS->Fill(rProdVtx);}
2401 }
3f18d0a6 2402 }
2308c814 2403 }
8d63376d 2404
dbf573fb 2405
2406 if(useTRKvtx) {
2407 // we need the vertex to compute the impact parameters
c0d39d2a 2408 if(!vertexESD) {delete trackTPC; trackTPC=0; continue;}
dbf573fb 2409 if(!(vertexESD->GetStatus()) || vertexESD->GetNContributors()<mincontrTRKvtx) {
2410 if(useSPDvtxifNotTRK) {
2411 vertexESD = fESD->GetPrimaryVertexSPD();
2412 } else {
c0d39d2a 2413 delete trackTPC; trackTPC=0; continue;
dbf573fb 2414 }
2415 }
2416 }
2417
0b2e8130 2418 //if(hasSharedSPD) continue;
d1d5a660 2419 //if(track->GetITSFakeFlag()) continue;
0b2e8130 2420
f8777463 2421 // track prolonged in ITS with different conditions
2422 if(itsrefit) {
2423 if(itsfindable) {
b7fb742d 2424 if(nclsITS==7) fHistPtITSMI7->Fill(track->Pt());
f8777463 2425 if(nclsITS==6) fHistPtITSMI6->Fill(track->Pt());
2426 if(nclsITS==5) fHistPtITSMI5->Fill(track->Pt());
2427 if(nclsITS==4) fHistPtITSMI4->Fill(track->Pt());
2428 if(nclsITS==3) fHistPtITSMI3->Fill(track->Pt());
2429 if(nclsITS==2) fHistPtITSMI2->Fill(track->Pt());
b7fb742d 2430 if(ConditionSPD(track)) {
f8777463 2431 fHistPtITSMISPD->Fill(track->Pt());
0b2e8130 2432 fHistTPCclsVSPtITSMISPDInAcc->Fill(track->Pt(),track->GetNcls(1));
2433 }
b7fb742d 2434 if(ConditionSPDone(track)) {
f8777463 2435 fHistPtITSMIoneSPD->Fill(track->Pt());
b7fb742d 2436 }
2437 if(ConditionSPDtwo(track)) {
2438 fHistPtITSMItwoSPD->Fill(track->Pt());
2439 }
f8777463 2440 if(nclsokbadoutinzITS==6) fHistPtITSMIokbadoutinz6->Fill(track->Pt());
2441 }
2442 if(itsfindableAcc) {
b7fb742d 2443 if(nclsITS==7) {
2444 fHistPtITSMI7InAcc->Fill(track->Pt());
2445 if(track->Pt()>1) {
2446 fHistPhiITSMI7InAcc->Fill(track->Phi());
2447 fHistEtaITSMI7InAcc->Fill(track->Eta());
2448 }
2449 if(isFake) fHistPtITSMI7InAccFake->Fill(track->Pt());
2450 if(isPrimary) {fHistPtITSMI7InAccP->Fill(track->Pt());} else {fHistPtITSMI7InAccS->Fill(track->Pt());}
2451 }
f8777463 2452 if(nclsITS==6) {
2453 fHistPtITSMI6InAcc->Fill(track->Pt());
31a96e36 2454 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI6InAccwSDD->Fill(track->Pt());
2455 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI6InAccwoSDD->Fill(track->Pt());
2456 if(isTOFbc0) {
2457 fHistPtITSMI6InAccTOFbc0->Fill(track->Pt());
2458 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI6InAccTOFbc0wSDD->Fill(track->Pt());
2459 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI6InAccTOFbc0woSDD->Fill(track->Pt());
2460 }
547640ff 2461 if(track->Pt()>1) {
2462 fHistPhiITSMI6InAcc->Fill(track->Phi());
2463 fHistEtaITSMI6InAcc->Fill(track->Eta());
2464 }
0faeb7ba 2465 if(isFake) fHistPtITSMI6InAccFake->Fill(track->Pt());
f8777463 2466 if(isPrimary) {fHistPtITSMI6InAccP->Fill(track->Pt());} else {fHistPtITSMI6InAccS->Fill(track->Pt());}
2467 }
2468 if(nclsITS==5) {
2469 fHistPtITSMI5InAcc->Fill(track->Pt());
31a96e36 2470 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI5InAccwSDD->Fill(track->Pt());
2471 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI5InAccwoSDD->Fill(track->Pt());
2472 if(isTOFbc0) {
2473 fHistPtITSMI5InAccTOFbc0->Fill(track->Pt());
2474 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI5InAccTOFbc0wSDD->Fill(track->Pt());
2475 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI5InAccTOFbc0woSDD->Fill(track->Pt());
2476 }
547640ff 2477 if(track->Pt()>1) {
2478 fHistPhiITSMI5InAcc->Fill(track->Phi());
2479 fHistEtaITSMI5InAcc->Fill(track->Eta());
2480 }
0faeb7ba 2481 if(isFake) fHistPtITSMI5InAccFake->Fill(track->Pt());
f8777463 2482 if(isPrimary) {fHistPtITSMI5InAccP->Fill(track->Pt());} else {fHistPtITSMI5InAccS->Fill(track->Pt());}
2483 }
2484 if(nclsITS==4) {
2485 fHistPtITSMI4InAcc->Fill(track->Pt());
31a96e36 2486 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI4InAccwSDD->Fill(track->Pt());
2487 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI4InAccwoSDD->Fill(track->Pt());
2488 if(isTOFbc0) {
2489 fHistPtITSMI4InAccTOFbc0->Fill(track->Pt());
2490 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI4InAccTOFbc0wSDD->Fill(track->Pt());
2491 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI4InAccTOFbc0woSDD->Fill(track->Pt());
2492 }
547640ff 2493 if(track->Pt()>1) {
2494 fHistPhiITSMI4InAcc->Fill(track->Phi());
2495 fHistEtaITSMI4InAcc->Fill(track->Eta());
2496 }
0faeb7ba 2497 if(isFake) fHistPtITSMI4InAccFake->Fill(track->Pt());
f8777463 2498 if(isPrimary) {fHistPtITSMI4InAccP->Fill(track->Pt());} else {fHistPtITSMI4InAccS->Fill(track->Pt());}
2499 }
2500 if(nclsITS==3) {
2501 fHistPtITSMI3InAcc->Fill(track->Pt());
31a96e36 2502 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI3InAccwSDD->Fill(track->Pt());
2503 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI3InAccwoSDD->Fill(track->Pt());
2504 if(isTOFbc0) {
2505 fHistPtITSMI3InAccTOFbc0->Fill(track->Pt());
2506 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI3InAccTOFbc0wSDD->Fill(track->Pt());
2507 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI3InAccTOFbc0woSDD->Fill(track->Pt());
2508 }
547640ff 2509 if(track->Pt()>1) {
2510 fHistPhiITSMI3InAcc->Fill(track->Phi());
2511 fHistEtaITSMI3InAcc->Fill(track->Eta());
2512 }
0faeb7ba 2513 if(isFake) fHistPtITSMI3InAccFake->Fill(track->Pt());
f8777463 2514 if(isPrimary) {fHistPtITSMI3InAccP->Fill(track->Pt());} else {fHistPtITSMI3InAccS->Fill(track->Pt());}
2515 }
2516 if(nclsITS==2) {
2517 fHistPtITSMI2InAcc->Fill(track->Pt());
31a96e36 2518 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI2InAccwSDD->Fill(track->Pt());
2519 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI2InAccwoSDD->Fill(track->Pt());
2520 if(isTOFbc0) {
2521 fHistPtITSMI2InAccTOFbc0->Fill(track->Pt());
2522 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMI2InAccTOFbc0wSDD->Fill(track->Pt());
2523 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMI2InAccTOFbc0woSDD->Fill(track->Pt());
2524 }
547640ff 2525 if(track->Pt()>1) {
2526 fHistPhiITSMI2InAcc->Fill(track->Phi());
2527 fHistEtaITSMI2InAcc->Fill(track->Eta());
2528 }
0faeb7ba 2529 if(isFake) fHistPtITSMI2InAccFake->Fill(track->Pt());
f8777463 2530 if(isPrimary) {fHistPtITSMI2InAccP->Fill(track->Pt());} else {fHistPtITSMI2InAccS->Fill(track->Pt());}
2531 }
b7fb742d 2532 if(ConditionSPD(track)) {
f8777463 2533 fHistPtITSMISPDInAcc->Fill(track->Pt());
31a96e36 2534 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMISPDInAccwSDD->Fill(track->Pt());
2535 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMISPDInAccwoSDD->Fill(track->Pt());
2536 if(isTOFbc0) {
2537 fHistPtITSMISPDInAccTOFbc0->Fill(track->Pt());
2538 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMISPDInAccTOFbc0wSDD->Fill(track->Pt());
2539 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMISPDInAccTOFbc0woSDD->Fill(track->Pt());
2540 }
547640ff 2541 if(track->Pt()>1) {
2542 fHistPhiITSMISPDInAcc->Fill(track->Phi());
2543 fHistEtaITSMISPDInAcc->Fill(track->Eta());
2544 }
0faeb7ba 2545 if(isFake) fHistPtITSMISPDInAccFake->Fill(track->Pt());
f8777463 2546 if(isPrimary) {fHistPtITSMISPDInAccP->Fill(track->Pt());} else {fHistPtITSMISPDInAccS->Fill(track->Pt());}
2547 }
b7fb742d 2548 if(ConditionSPDtwo(track)) {
2549 fHistPtITSMItwoSPDInAcc->Fill(track->Pt());
2550 if(track->Pt()>1) {
2551 fHistPhiITSMItwoSPDInAcc->Fill(track->Phi());
2552 fHistEtaITSMItwoSPDInAcc->Fill(track->Eta());
2553 }
2554 if(isFake) fHistPtITSMItwoSPDInAccFake->Fill(track->Pt());
2555 if(isPrimary) {fHistPtITSMItwoSPDInAccP->Fill(track->Pt());} else {fHistPtITSMItwoSPDInAccS->Fill(track->Pt());}
2556 }
2557 if(ConditionSPDone(track)) {
f8777463 2558 fHistPtITSMIoneSPDInAcc->Fill(track->Pt());
31a96e36 2559 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMIoneSPDInAccwSDD->Fill(track->Pt());
2560 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMIoneSPDInAccwoSDD->Fill(track->Pt());
2561 if(isTOFbc0) {
2562 fHistPtITSMIoneSPDInAccTOFbc0->Fill(track->Pt());
2563 if(fCheckSDDIsIn && sddIsIn) fHistPtITSMIoneSPDInAccTOFbc0wSDD->Fill(track->Pt());
2564 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSMIoneSPDInAccTOFbc0woSDD->Fill(track->Pt());
2565 }
547640ff 2566 if(track->Pt()>1) {
2567 fHistPhiITSMIoneSPDInAcc->Fill(track->Phi());
2568 fHistEtaITSMIoneSPDInAcc->Fill(track->Eta());
2569 }
0b2e8130 2570 if(track->HasSharedPointOnITSLayer(0) ||
2571 track->HasSharedPointOnITSLayer(1) ||
2572 track->HasSharedPointOnITSLayer(2) ||
2573 track->HasSharedPointOnITSLayer(3) ||
2574 track->HasSharedPointOnITSLayer(4) ||
2575 track->HasSharedPointOnITSLayer(5)) {
2576 fHistPtITSMIoneSPDInAccShared->Fill(track->Pt());
2577 if(isFake) fHistPtITSMIoneSPDInAccSharedFake->Fill(track->Pt());
2578 }
2579 if(track->HasSharedPointOnITSLayer(0) || track->HasSharedPointOnITSLayer(1)) {
2580 fHistPtITSMIoneSPDInAccSharedSPD->Fill(track->Pt());
2581 if(isFake) fHistPtITSMIoneSPDInAccSharedSPDFake->Fill(track->Pt());
2582 }
2583 if(track->HasSharedPointOnITSLayer(0)) fHistPtITSMISPD1InAccShared->Fill(track->Pt());
2584 if(track->HasSharedPointOnITSLayer(1)) fHistPtITSMISPD2InAccShared->Fill(track->Pt());
2585
3df234e9 2586
d8ec295e 2587 Float_t chi2redITS = 0.;
2588 if(track->GetNcls(0)>0) chi2redITS = track->GetITSchi2()/(Float_t)track->GetNcls(0);
0b2e8130 2589
af874c53 2590 if(track->Pt()>0.18 && track->Pt()<0.25) fHistITSRedChi2NonFakePt02->Fill(chi2redITS);
2591 if(track->Pt()>0.4 && track->Pt()<0.6) fHistITSRedChi2NonFakePt05->Fill(chi2redITS);
2592 if(track->Pt()>0.9 && track->Pt()<1.5) fHistITSRedChi2NonFakePt1->Fill(chi2redITS);
2593 if(isFake) {
2594 fHistPtITSMIoneSPDInAccFake->Fill(track->Pt());
2595 if(track->Pt()>0.18 && track->Pt()<0.25) fHistITSRedChi2FakePt02->Fill(chi2redITS);
2596 if(track->Pt()>0.4 && track->Pt()<0.6) fHistITSRedChi2FakePt05->Fill(chi2redITS);
2597 if(track->Pt()>0.9 && track->Pt()<1.5) fHistITSRedChi2FakePt1->Fill(chi2redITS);
2598 }
6c6fab7b 2599 if(nclsSDDSSD>=3) fHistPtITSMIoneSPDthreeSDDSSDInAcc->Fill(track->Pt());
2600
f8777463 2601 if(isPrimary) {fHistPtITSMIoneSPDInAccP->Fill(track->Pt());} else {fHistPtITSMIoneSPDInAccS->Fill(track->Pt());}
2602 }
2603 if(nclsokbadoutinzITS==6) fHistPtITSMIokbadoutinz6InAcc->Fill(track->Pt());
2604 if(nclsokbadoutinzITS==5) fHistPtITSMIokbadoutinz5InAcc->Fill(track->Pt());
2605 if(nclsokbadoutinzITS==4) fHistPtITSMIokbadoutinz4InAcc->Fill(track->Pt());
547640ff 2606 if(nclsokbadoutinzITS==6 && track->Pt()>1) fHistPhiITSMIokbadoutinz6InAcc->Fill(track->Phi());
f8777463 2607 }
2608 }
8d63376d 2609
f8777463 2610 if(tpcrefit) {
2308c814 2611 fHistNclsITSMI->Fill(nclsITS);
2612 } else {
2613 fHistNclsITSSA->Fill(nclsITS);
2614 if(!outInZ) fHistNclsITSSAInAcc->Fill(nclsITS);
2615 }
2616
f8777463 2617
c0d39d2a 2618 if((tpcrefit && fUseITSSAforNtuples) || // only ITS-SA for ntuples
2619 (!tpcrefit && !fUseITSSAforNtuples)) // only ITS-TPC for ntuples
2620 { delete trackTPC; trackTPC=0; continue; }
8d63376d 2621
f8777463 2622 // impact parameter to VertexTracks
e9204849 2623 Float_t d0z0[2]={0,0},covd0z0[3]={0,0,0};
2624 Double_t d0z0TPC[2]={0,0},covd0z0TPC[3]={0,0,0};
dbf573fb 2625 if(useTRKvtx) {
c0d39d2a 2626 if(!track->RelateToVertex(vertexESD,fESD->GetMagneticField(),kVeryBig)) { delete trackTPC; trackTPC=0; continue; }
dbf573fb 2627 } else {
2628 //if(!track->RelateToVertex(vertexMC,fESD->GetMagneticField(),kVeryBig)) continue;
c0d39d2a 2629 if(!track->RelateToVertex(spdv,fESD->GetMagneticField(),kVeryBig)) { delete trackTPC; trackTPC=0; continue; }
dbf573fb 2630 }
f8777463 2631 track->GetImpactParameters(d0z0,covd0z0);
2632 if(trackTPC) trackTPC->PropagateToDCA(spdv,fESD->GetMagneticField(),kVeryBig,d0z0TPC,covd0z0TPC);
c0d39d2a 2633 if(covd0z0[0]<0. || covd0z0[2]<0. || covd0z0TPC[0]<0. || covd0z0TPC[2]<0.) { delete trackTPC; trackTPC=0; continue; }
2634
f8777463 2635
2636 // track that passes final ITS+TPC cuts
2637 if(itsfindableAcc && fESDtrackCutsITSTPC->AcceptTrack(track)) {
2638 fHistPtITSTPCsel->Fill(track->Pt());
31a96e36 2639 if(fCheckSDDIsIn && sddIsIn) fHistPtITSTPCselwSDD->Fill(track->Pt());
2640 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSTPCselwoSDD->Fill(track->Pt());
2641 if(isTOFbc0) {
2642 fHistPtITSTPCselTOFbc0->Fill(track->Pt());
2643 if(fCheckSDDIsIn && sddIsIn) fHistPtITSTPCselTOFbc0wSDD->Fill(track->Pt());
2644 if(fCheckSDDIsIn && !sddIsIn) fHistPtITSTPCselTOFbc0woSDD->Fill(track->Pt());
2645 }
af874c53 2646 if(isFake) fHistPtITSTPCselFake->Fill(track->Pt());
f8777463 2647 fHistdEdxVSPtITSTPCsel->Fill(track->Pt(),track->GetITSsignal());
2648 if(isPrimary) {
2649 fHistPtITSTPCselP->Fill(track->Pt());
2650 if(pdgTrk==321) fHistPtITSTPCselPfromStrange->Fill(track->Pt());
2651 } else {
2652 fHistPtITSTPCselS->Fill(track->Pt());
2653 if(isFromMat) fHistPtITSTPCselSfromMat->Fill(track->Pt());
2654 if(isFromStrange) fHistPtITSTPCselSfromStrange->Fill(track->Pt());
2655 }
2656 }
2657
2658
2659 // fill d0 histos
dbf573fb 2660 if(((!fUseITSSAforNtuples&&itsfindableAcc) || fUseITSSAforNtuples) /*&& track->Charge()>0*/) {
f8777463 2661 if(track->Pt()>0.150 && track->Pt()<0.200) {
2662 if(isPrimary) {
2663 fHistd0rphiTPCInAccP150200->Fill(d0z0TPC[0]);
b7fb742d 2664 if(ConditionSPD(track) && itsrefit)
f8777463 2665 fHistd0rphiITSMISPDInAccP150200->Fill(d0z0[0]);
b7fb742d 2666 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2667 //if(((track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1)) || (!track->HasPointOnITSLayer(0) && track->HasPointOnITSLayer(1))) && itsrefit) {
f8777463 2668 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) fHistd0rphiITSMIoneSPDInAccP150200->Fill(d0z0[0]);
2669 if(TMath::Abs(d0z0[0])<maxdcaxyITSTPC) fHistd0zITSMIoneSPDInAccP150200->Fill(d0z0[1]);
2670 }
2671 } else {
2672 fHistd0rphiTPCInAccS150200->Fill(d0z0TPC[0]);
b7fb742d 2673 if(ConditionSPD(track) && itsrefit)
f8777463 2674 fHistd0rphiITSMISPDInAccS150200->Fill(d0z0[0]);
b7fb742d 2675 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2676 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) {
2677 fHistd0rphiITSMIoneSPDInAccS150200->Fill(d0z0[0]);
2678 if(isFromStrange) fHistd0rphiITSMIoneSPDInAccS150200fromStrange->Fill(d0z0[0]);
2679 if(isFromMat) fHistd0rphiITSMIoneSPDInAccS150200fromMat->Fill(d0z0[0]);
2680 }
f8777463 2681 if(TMath::Abs(d0z0[0])<maxdcaxyITSTPC) fHistd0zITSMIoneSPDInAccS150200->Fill(d0z0[1]);
2682 if(TMath::Abs(d0z0[0])<0.1) fHistPDGMoth150200->Fill(pdgMoth);
2683 }
2684 }
2685 }
dbf573fb 2686 if(track->Pt()>0.350 && track->Pt()<0.450) {
2687 if(isPrimary) {
b7fb742d 2688 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2689 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) fHistd0rphiITSMIoneSPDInAccP350450->Fill(d0z0[0]);
2690 }
2691 } else {
b7fb742d 2692 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2693 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) {
2694 fHistd0rphiITSMIoneSPDInAccS350450->Fill(d0z0[0]);
2695 if(isFromStrange) fHistd0rphiITSMIoneSPDInAccS350450fromStrange->Fill(d0z0[0]);
2696 if(isFromMat) fHistd0rphiITSMIoneSPDInAccS350450fromMat->Fill(d0z0[0]);
2697 }
2698 }
2699 }
2700 }
f8777463 2701 if(track->Pt()>0.500 && track->Pt()<0.700) {
2702 if(isPrimary) {
2703 fHistd0rphiTPCInAccP500700->Fill(d0z0TPC[0]);
b7fb742d 2704 if(ConditionSPD(track) && itsrefit)
f8777463 2705 fHistd0rphiITSMISPDInAccP500700->Fill(d0z0[0]);
b7fb742d 2706 if((ConditionSPDone(track)) && itsrefit) {
f8777463 2707 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) fHistd0rphiITSMIoneSPDInAccP500700->Fill(d0z0[0]);
2708 if(TMath::Abs(d0z0[0])<maxdcaxyITSTPC) fHistd0zITSMIoneSPDInAccP500700->Fill(d0z0[1]);
2709 }
2710 } else {
2711 fHistd0rphiTPCInAccS500700->Fill(d0z0TPC[0]);
b7fb742d 2712 if(ConditionSPD(track) && itsrefit)
f8777463 2713 fHistd0rphiITSMISPDInAccS500700->Fill(d0z0[0]);
b7fb742d 2714 if((ConditionSPDone(track)) && itsrefit) {
f8777463 2715 fHistPDGTrk->Fill(pdgTrk);
2716 fHistPDGMoth->Fill(pdgMoth);
2717 if(TMath::Abs(d0z0[0])<0.1) fHistPDGMoth500700->Fill(pdgMoth);
dbf573fb 2718 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) {
2719 fHistd0rphiITSMIoneSPDInAccS500700->Fill(d0z0[0]);
2720 if(isFromStrange) fHistd0rphiITSMIoneSPDInAccS500700fromStrange->Fill(d0z0[0]);
2721 if(isFromMat) fHistd0rphiITSMIoneSPDInAccS500700fromMat->Fill(d0z0[0]);
2722 }
f8777463 2723 if(TMath::Abs(d0z0[0])<maxdcaxyITSTPC) fHistd0zITSMIoneSPDInAccS500700->Fill(d0z0[1]);
2724 if(pdgMoth==310) fHistd0rphiITSMIoneSPDInAccS500700from310->Fill(d0z0[0]);
2725 if(pdgMoth==321) fHistd0rphiITSMIoneSPDInAccS500700from321->Fill(d0z0[0]);
dbf573fb 2726 if(pdgMoth==3122) fHistd0rphiITSMIoneSPDInAccS500700from3122->Fill(d0z0[0]);
f8777463 2727 if(pdgMoth==211) fHistd0rphiITSMIoneSPDInAccS500700from211->Fill(d0z0[0]);
2728 if(pdgMoth==22) fHistd0rphiITSMIoneSPDInAccS500700from22->Fill(d0z0[0]);
dbf573fb 2729 if(pdgMoth==211 && TMath::Abs(d0z0[0])>0.5) fHistRProdVtxInAccS->Fill(rProdVtx);
f8777463 2730 }
2731 }
2732 }
2733 if(track->Pt()>1.000 && track->Pt()<1.500) {
2734 if(isPrimary) {
2735 fHistd0rphiTPCInAccP10001500->Fill(d0z0TPC[0]);
9a3cd95a 2736 fHistd0zVSetaTPCInAccP10001500->Fill(d0z0TPC[1],trackTPC->Eta());
2737
b7fb742d 2738 if(ConditionSPD(track) && itsrefit)
f8777463 2739 fHistd0rphiITSMISPDInAccP10001500->Fill(d0z0[0]);
b7fb742d 2740 if((ConditionSPDone(track)) && itsrefit) {
f8777463 2741 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) fHistd0rphiITSMIoneSPDInAccP10001500->Fill(d0z0[0]);
2742 if(TMath::Abs(d0z0[0])<maxdcaxyITSTPC) fHistd0zITSMIoneSPDInAccP10001500->Fill(d0z0[1]);
2743 fHistd0rphiVSphiITSMIoneSPDInAccP10001500->Fill(d0z0[0],track->Phi());
2744 fHistd0rphiVSetaITSMIoneSPDInAccP10001500->Fill(d0z0[0],track->Eta());
2745 }
2746 } else {
2747 fHistd0rphiTPCInAccS10001500->Fill(d0z0TPC[0]);
b7fb742d 2748 if(ConditionSPD(track) && itsrefit)
f8777463 2749 fHistd0rphiITSMISPDInAccS10001500->Fill(d0z0[0]);
b7fb742d 2750 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2751 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) {
2752 fHistd0rphiITSMIoneSPDInAccS10001500->Fill(d0z0[0]);
2753 if(isFromStrange) fHistd0rphiITSMIoneSPDInAccS10001500fromStrange->Fill(d0z0[0]);
2754 if(isFromMat) fHistd0rphiITSMIoneSPDInAccS10001500fromMat->Fill(d0z0[0]);
2755 }
f8777463 2756 if(TMath::Abs(d0z0[0])<maxdcaxyITSTPC) fHistd0zITSMIoneSPDInAccS10001500->Fill(d0z0[1]);
2757 if(TMath::Abs(d0z0[0])<0.1) fHistPDGMoth10001500->Fill(pdgMoth);
2758 }
2759 }
2760 }
dbf573fb 2761 if(track->Pt()>2.500 && track->Pt()<4.000) {
2762 if(isPrimary) {
b7fb742d 2763 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2764 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) fHistd0rphiITSMIoneSPDInAccP25004000->Fill(d0z0[0]);
2765 }
2766 } else {
b7fb742d 2767 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2768 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) {
2769 fHistd0rphiITSMIoneSPDInAccS25004000->Fill(d0z0[0]);
2770 if(isFromStrange) fHistd0rphiITSMIoneSPDInAccS25004000fromStrange->Fill(d0z0[0]);
2771 if(isFromMat) fHistd0rphiITSMIoneSPDInAccS25004000fromMat->Fill(d0z0[0]);
2772 }
2773 }
2774 }
2775 }
2776 if(track->Pt()>4.000 && track->Pt()<8.000) {
2777 if(isPrimary) {
b7fb742d 2778 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2779 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) fHistd0rphiITSMIoneSPDInAccP40008000->Fill(d0z0[0]);
2780 }
2781 } else {
b7fb742d 2782 if((ConditionSPDone(track)) && itsrefit) {
dbf573fb 2783 if(TMath::Abs(d0z0[1])<maxdcazITSTPC) {
2784 fHistd0rphiITSMIoneSPDInAccS40008000->Fill(d0z0[0]);
2785 if(isFromStrange) fHistd0rphiITSMIoneSPDInAccS40008000fromStrange->Fill(d0z0[0]);
2786 if(isFromMat) fHistd0rphiITSMIoneSPDInAccS40008000fromMat->Fill(d0z0[0]);
2787 }
2788 }
2789 }
2790 }
2791
8d63376d 2792 }
2793
c0d39d2a 2794 if(trackTPC) { delete trackTPC; trackTPC=0; }
2308c814 2795
f8777463 2796 // encode ITS cluster map, including MC info
2797 Int_t iITSflag=MakeITSflag(track);
8d63376d 2798
2799 // if MC info is available: get particle properties
2800 Float_t ptMC=-999.,pdgMC=-999.,d0MC=-999.;
2801 Double_t d0z0MCv[2]={-999.,-999.},covd0z0MCv[3]={1.,1.,1.};
2802 if(fReadMC) {
2803 part = (TParticle*)stack->Particle(trkLabel);
2804 ptMC=part->Pt();
2805 pdgMC=part->GetPdgCode();
2806 d0MC=ParticleImpParMC(part,vertexMC,0.1*fESD->GetMagneticField());
2807 track->PropagateToDCA(vertexMC,fESD->GetMagneticField(),100.,d0z0MCv,covd0z0MCv);
2808 if(covd0z0MCv[0]<0. || covd0z0MCv[2]<0.) continue;
8d63376d 2809 }
2810
ab846928 2811 Double_t sigmad0MCv=TMath::Sqrt(covd0z0MCv[0]);
2812 if(!itsrefit) sigmad0MCv *= -1.;
f8777463 2813 Float_t isPrimaryFl = (isPrimary ? 1. : 0.);
2814 Float_t isTPCSelFl = (itsfindableAcc ? 1. : 0.);
ab846928 2815
8d63376d 2816 // fill ntuple with track properties
2308c814 2817 if(fFillNtuples && SelectPt(track->Pt())) {
2942f542 2818 Float_t fillArray[22]={static_cast<Float_t>(track->Pt()),static_cast<Float_t>(track->Eta()),static_cast<Float_t>(track->Phi()),static_cast<Float_t>(d0z0[0]),static_cast<Float_t>(d0z0[1]),static_cast<Float_t>(TMath::Sqrt(covd0z0[0])),static_cast<Float_t>(TMath::Sqrt(covd0z0[2])),ptMC,static_cast<Float_t>(pdgMC),static_cast<Float_t>(pdgMoth),static_cast<Float_t>(d0MC),static_cast<Float_t>(d0z0MCv[0]),static_cast<Float_t>(d0z0MCv[1]),static_cast<Float_t>(sigmad0MCv),static_cast<Float_t>(TMath::Sqrt(covd0z0MCv[2])),static_cast<Float_t>(iITSflag),static_cast<Float_t>(isPrimaryFl),static_cast<Float_t>(isTPCSelFl)};
8d63376d 2819 fNtupleESDTracks->Fill(fillArray);
2820 }
2821
2822 //---------------------------------------------
2823 // AliTrackPoints: alignment checks
2824 //
2308c814 2825 if(!fFillNtuples) continue;
ab846928 2826
8d63376d 2827 const AliTrackPointArray *array = track->GetTrackPointArray();
2308c814 2828 if(!array) continue;
8d63376d 2829 AliTrackPoint point;
2830 Int_t pointOnLayer[6]={0,0,0,0,0,0};
2831 Int_t indexAssociated[6]={-1,-1,-1,-1,-1,-1},indexExtra=-1;
2832 Bool_t extra=kFALSE;
2833 Int_t layerId,layerExtra=-1;
2834 for(Int_t ipt=0; ipt<array->GetNPoints(); ipt++) {
2835 array->GetPoint(point,ipt);
2836 Float_t r = TMath::Sqrt(point.GetX()*point.GetX()+point.GetY()*point.GetY());
2837
2838 if(r>3 && r<6) {
2839 layerId = 0;
2840 } else if(r>6 && r<8) {
2841 layerId = 1;
2842 } else if(r>8 && r<18) {
2843 layerId = 2;
3f18d0a6 2844 } else if(r>18 && r<30) {
8d63376d 2845 layerId = 3;
3f18d0a6 2846 } else if(r>30 && r<40) {
8d63376d 2847 layerId = 4;
3f18d0a6 2848 } else if(r>40 && r<50) {
8d63376d 2849 layerId = 5;
2850 } else {
2851 layerId=100;
2852 }
2853
2854 // only ITS points
2855 if(layerId>5) continue;
2856
2857 if(!point.IsExtra()) {
2858 pointOnLayer[layerId]++;
2859 indexAssociated[layerId]=ipt;
2860 } else {
2861 // this is an extra cluster
2862 extra=kTRUE;
2863 layerExtra=layerId;
2864 indexExtra=ipt;
2865 }
2866 } // end loop on AliTrackPoints
2867
2868 TString vtitle = spdv->GetTitle();
2869 if(!vtitle.Contains("3D")) continue;
2870
2871 // SPD tracklet
2872 if(indexAssociated[0]>=0 && indexAssociated[1]>=0) {
2873 AliTrackPoint pointSPD1,pointSPD2;
2874 array->GetPoint(pointSPD1,indexAssociated[0]);
2875 array->GetPoint(pointSPD2,indexAssociated[1]);
2876 Float_t phi=TMath::ATan2(pointSPD2.GetY()-pointSPD1.GetY(),pointSPD2.GetX()-pointSPD1.GetX());
2877 Float_t lambda=TMath::ATan((pointSPD2.GetZ()-pointSPD1.GetZ())/TMath::Sqrt((pointSPD2.GetX()-pointSPD1.GetX())*(pointSPD2.GetX()-pointSPD1.GetX())+(pointSPD2.GetY()-pointSPD1.GetY())*(pointSPD2.GetY()-pointSPD1.GetY())));
2878 Float_t theta=0.5*TMath::Pi()-lambda;
2879 TParticle particle(211,0,0,0,0,0,TMath::Cos(phi),TMath::Sin(phi),TMath::Tan(lambda),10.,pointSPD1.GetX(),pointSPD1.GetY(),pointSPD1.GetZ(),0);
2880 AliESDtrack tracklet(&particle);
2881 Float_t dz[2];
2882 // distance to primary SPD (only if 3D and high multiplicity)
2883 if(spdv->GetNContributors()>10) {
e690d4d0 2884 tracklet.GetDZ(spdv->GetX(),spdv->GetY(),spdv->GetZ(),0,dz);
2885 //tracklet.GetDZ(-0.07,0.25,spdv->GetZ(),0,dz);
c2893e82 2886 fNtupleITSAlignSPDTracklets->Fill(phi,theta,0.5*(pointSPD1.GetZ()+pointSPD2.GetZ()),dz[0],dz[1],track->Pt());
8d63376d 2887 }
2888 }
2889
2890 // distance to extra
3f18d0a6 2891 if(extra && spdv->GetNContributors()>4 && indexAssociated[layerExtra]>-1) {
8d63376d 2892 AliTrackPoint pointExtra,pointAssociated;
2893 array->GetPoint(pointAssociated,indexAssociated[layerExtra]);
2894 array->GetPoint(pointExtra,indexExtra);
e690d4d0 2895 Float_t phiExtra = TMath::ATan2(pointExtra.GetY()-spdv->GetY(),pointExtra.GetX()-spdv->GetX());
2896 Float_t phiAssociated = TMath::ATan2(pointAssociated.GetY()-spdv->GetY(),pointAssociated.GetX()-spdv->GetX());
2897 Float_t rExtra = TMath::Sqrt((pointExtra.GetX()-spdv->GetX())*(pointExtra.GetX()-spdv->GetX())+(pointExtra.GetY()-spdv->GetY())*(pointExtra.GetY()-spdv->GetY()));
2898 Float_t rAssociated = TMath::Sqrt((pointAssociated.GetX()-spdv->GetX())*(pointAssociated.GetX()-spdv->GetX())+(pointAssociated.GetY()-spdv->GetY())*(pointAssociated.GetY()-spdv->GetY()));
8d63376d 2899 Float_t dzExtra[2];
2900 dzExtra[0] = (phiExtra-phiAssociated)*0.5*(rExtra+rAssociated);
e690d4d0 2901 dzExtra[1] = pointExtra.GetZ()-pointAssociated.GetZ()-(rExtra-rAssociated)*(pointAssociated.GetZ()-spdv->GetZ())/rAssociated;
8d63376d 2902 Float_t xlocExtra=-100.,zlocExtra=-100.;
c2893e82 2903 fNtupleITSAlignExtra->Fill(layerExtra,pointExtra.GetX(),pointExtra.GetY(),pointExtra.GetZ(),dzExtra[0],dzExtra[1],xlocExtra,zlocExtra,nclsITS,track->Pt());
8d63376d 2904 }
2905
f8777463 2906 if(trackTPC) { delete trackTPC; trackTPC=0; }
8d63376d 2907 } // end loop on tracks
2908
2909 if(vertexMC) { delete vertexMC; vertexMC=0; }
2910
f8777463 2911 delete esdtrackCutsTPC; esdtrackCutsTPC=0;
2912 delete esdtrackCutsITSTPC; esdtrackCutsITSTPC=0;
2913
0b2e8130 2914 PostData(1, fOutput);
2915
8d63376d 2916 return;
2917}
2918
2919//________________________________________________________________________
2920void AliAnalysisTaskITSTrackingCheck::Terminate(Option_t *)
2921{
2922 // Draw result to the screen
2923 // Called once at the end of the query
2efff183 2924 fOutput = dynamic_cast<TList*> (GetOutputData(1));
8d63376d 2925 if (!fOutput) {
2926 Printf("ERROR: fOutput not available");
2927 return;
2928 }
f8777463 2929 fHistNEvents = dynamic_cast<TH1F*>(fOutput->FindObject("fHistNEvents"));
2930 fHistNEventsFrac = dynamic_cast<TH1F*>(fOutput->FindObject("fHistNEventsFrac"));
2931 if(fHistNEvents && fHistNEventsFrac) {
2932 for(Int_t ibin=2; ibin<=fHistNEvents->GetNbinsX(); ibin++) {
2933 if(fHistNEvents->GetBinContent(2)) fHistNEventsFrac->SetBinContent(ibin,fHistNEvents->GetBinContent(ibin)/fHistNEvents->GetBinContent(2));
2934 }
2935 }
8d63376d 2936
2937 return;
2938}
2939//---------------------------------------------------------------------------
dbf573fb 2940Int_t AliAnalysisTaskITSTrackingCheck::NumberOfITSClustersMC(Int_t label,Int_t nModules) const
8d63376d 2941{
2942 //
2943 // Return number of ITS clusters produced by MC particle with given label
2944 //
2308c814 2945
8d63376d 2946 AliESDInputHandlerRP *esdHRP = dynamic_cast<AliESDInputHandlerRP*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
2947 if(!esdHRP) return -1;
2948 TTree *cTree = (TTree*)esdHRP->GetTreeR("ITS");
2949 if(!cTree) return -1;
dbf573fb 2950 TClonesArray *clusters=0;
8d63376d 2951 cTree->SetBranchAddress("ITSRecPoints",&clusters);
2952 if(!clusters) return -1;
2953
2954 AliITSRecPoint *c=0;
2955 Int_t i,n,icl,lay,ilab;
b7fb742d 2956 Int_t ncls[7]={0,0,0,0,0,0,0};
8d63376d 2957 Int_t nclstot=0;
2958
dbf573fb 2959 for(i=0; i<nModules; i++) {
8d63376d 2960 cTree->GetEvent(i);
2961 n=clusters->GetEntriesFast();
2962 for (icl=0; icl<n; icl++) {
2963 c=(AliITSRecPoint*)clusters->UncheckedAt(icl);
2964 lay=c->GetLayer();
dbf573fb 2965 //printf("mod %d lay %d xloc %f zloc %f\n",i,lay,c->GetDetLocalX(),c->GetDetLocalZ());
8d63376d 2966 for(ilab=0;ilab<3;ilab++) {
2967 if(c->GetLabel(ilab)==label) ncls[lay]++;
2968 }
2969 }
2970 }
b7fb742d 2971 for(i=0;i<fNITSLayers;i++) {
dbf573fb 2972 //printf("%d\n",ncls[i]);
2973 if(ncls[i]) nclstot++;
2974 }
8d63376d 2975
2976 return nclstot;
2977}
2978//---------------------------------------------------------------------------
f8777463 2979Int_t AliAnalysisTaskITSTrackingCheck::NumberOfITSClusters(Int_t idet,Float_t &xloc) const
2980{
2981 //
2982 // Return number of ITS clusters produced by MC particle with given label
2983 //
2984
2985 AliESDInputHandlerRP *esdHRP = dynamic_cast<AliESDInputHandlerRP*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
2986 if(!esdHRP) return -1;
2987 TTree *cTree = (TTree*)esdHRP->GetTreeR("ITS");
2988 if(!cTree) return -1;
2989 TClonesArray *clusters=0; // new TClonesArray("AliITSRecPoint",10000);
2990 cTree->SetBranchAddress("ITSRecPoints",&clusters);
2991 if(!clusters) return -1;
2992
2993 AliITSRecPoint *c=0;
2994 Int_t n,icl;
2995
2996 cTree->GetEvent(idet);
2997 n=clusters->GetEntriesFast();
2998
2999 if(n==1) {
3000 for (icl=0; icl<n; icl++) {
3001 c=(AliITSRecPoint*)clusters->UncheckedAt(icl);
3002 xloc = c->GetDetLocalX();
3003 }
3004 }
3005 return n;
3006}
3007//---------------------------------------------------------------------------
0b2e8130 3008void AliAnalysisTaskITSTrackingCheck::FillNClustersModuleMap()
3009{
3010 //
3011 // Number of clusters per module
3012 //
3013
3014 AliESDInputHandlerRP *esdHRP = dynamic_cast<AliESDInputHandlerRP*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
3015 if(!esdHRP) return;
3016 TTree *cTree = (TTree*)esdHRP->GetTreeR("ITS");
3017 if(!cTree) return;
3018 TClonesArray *clusters=0; // new TClonesArray("AliITSRecPoint",10000);
3019 cTree->SetBranchAddress("ITSRecPoints",&clusters);
3020
3021 for(Int_t idet=0;idet<2198;idet++) {
3022 cTree->GetEvent(idet);
3023 if(!clusters) continue;
3024 fHistNClustersMapModule->SetBinContent(idet,clusters->GetEntriesFast());
3025 }
3026
3027 return;
3028}
3029//---------------------------------------------------------------------------
8d63376d 3030Double_t AliAnalysisTaskITSTrackingCheck::ParticleImpParMC(TParticle *part,
3031 AliESDVertex *vert,
3032 Double_t bzT) const
3033{
3034 //
3035 // Return the MC value of the impact parameter
3036 //
b7fb742d 3037
3038 bzT = TMath::Abs(bzT);
3039
8d63376d 3040
3041 Double_t vx=part->Vx()-vert->GetX();
3042 Double_t vy=part->Vy()-vert->GetY();
3043
3044 Double_t pt=part->Pt();
3045 Double_t px=part->Px();
3046 Double_t py=part->Py();
3047 Double_t charge = (part->GetPdgCode()>0. ? 1. : -1.);
3048 if(TMath::Abs(part->GetPdgCode())<100) charge*=-1.;
3049
b7fb742d 3050
3051 if(TMath::Abs(px)<0.000001) px=0.000001;
f8777463 3052 Double_t rAnd=((10./2.99792458)*pt/bzT)*100.;
8d63376d 3053 Double_t center[3],d0;
3054 center[0]=vx-(1./charge)*rAnd*(py/pt);
3055 center[1]=vy+(1./charge)*rAnd*(px/pt);
3056 center[2]=TMath::Sqrt(center[0]*center[0]+center[1]*center[1]);
3057 d0 = -center[2]+rAnd;
3058
b7fb742d 3059 //if(pt>29) printf("%f %f %f %d %f %f %f %f %f\n",px,py,pt,part->GetPdgCode(),center[0],center[1],center[2],rAnd,d0);
3060
8d63376d 3061 return d0;
3062}
3063//---------------------------------------------------------------------------
3064Bool_t AliAnalysisTaskITSTrackingCheck::SelectPt(Double_t pt)
3065{
3066 //
3067 // Keep only tracks in given pt bins
3068 //
f8777463 3069 Double_t ptlower[11]={0.15,0.29,0.49,0.75,0.9,1.9,3.5,6.5, 9.,19.,27.};
3070 Double_t ptupper[11]={0.16,0.31,0.51,0.85,1.1,2.1,4.5,7.5,11.,21.,33.};
8d63376d 3071
f8777463 3072 for(Int_t i=0; i<11; i++) {
ab846928 3073 if(pt>ptlower[i] && pt<ptupper[i]) {
8d63376d 3074 fCountsPerPtBin[i]++;
3075 return kTRUE;
3076 }
3077 }
3078 return kFALSE;
2308c814 3079 //return kTRUE;
8d63376d 3080}
f8777463 3081//---------------------------------------------------------------------------
3082Int_t AliAnalysisTaskITSTrackingCheck::MakeITSflag(AliESDtrack *track) const {
3083 //
3084 // ITSflag takes the value 0 if the track has no cluster assigned in the SPDs,
3085 // 1 (2) if one cluster is assigned in SPD1(2), 3 if two clusters are present.
3086 // Then the same adding 10,20 or 30 for SDD and 100,200 or 300 for SSD
3087 //
3088 Int_t iITSflag=0;
b7fb742d 3089 if(track->HasPointOnITSLayer(0)) iITSflag+=1;
3090 if(track->HasPointOnITSLayer(fNITSLayers-5)) iITSflag+=2;
3091 if(track->HasPointOnITSLayer(fNITSLayers-4)) iITSflag+=10;
3092 if(track->HasPointOnITSLayer(fNITSLayers-3)) iITSflag+=20;
3093 if(track->HasPointOnITSLayer(fNITSLayers-2)) iITSflag+=100;
3094 if(track->HasPointOnITSLayer(fNITSLayers-1)) iITSflag+=200;
f8777463 3095
3096 // number of associated ITS clusters
3097 iITSflag += 1000*track->GetNcls(0);
3098
3099 // number of associated TPC clusters
3100 iITSflag += 100000*track->GetNcls(1);
3101
3102 // if MC info and is available
3103 // write the number of ITS clusters produced by this track
3104 Int_t nITSclsMC=0;
3105 if(fReadMC && fReadRPLabels) {
3106 nITSclsMC = NumberOfITSClustersMC(TMath::Abs(track->GetLabel()));
3107 if(nITSclsMC>=0) iITSflag += 10000*nITSclsMC;
3108 // flag fake tracks
3109 if(track->GetLabel()<0) iITSflag *= -1;
3110 }
8d63376d 3111
f8777463 3112 return iITSflag;
3113}
5ac4a99f 3114//---------------------------------------------------------------------------
3115Bool_t AliAnalysisTaskITSTrackingCheck::IsSelectedCentrality() const
3116{
3117 //
3118 // check if events is in the required multiplicity range
3119 //
3120
3121 const AliMultiplicity *alimult = fESD->GetMultiplicity();
3122 Int_t ntrklets=1;
2da5db49 3123 Int_t nclsSPDouter=0;
5ac4a99f 3124 if(alimult) {
3125 ntrklets = alimult->GetNumberOfTracklets();
2da5db49 3126 nclsSPDouter = alimult->GetNumberOfITSClusters(1);
5ac4a99f 3127 }
3128
2da5db49 3129 if(nclsSPDouter<fMinMult || nclsSPDouter>fMaxMult) return kFALSE;
3130
5ac4a99f 3131
3132 return kTRUE;
3133}
b7fb742d 3134//---------------------------------------------------------------------------
3135Bool_t AliAnalysisTaskITSTrackingCheck::ConditionSPD(AliESDtrack *track) const
3136{
3137 //
3138 // check if track has points in all inner barrel layers
3139 //
8d63376d 3140
3141
b7fb742d 3142 if(NPointsInnerBarrel(track)==(fNITSLayers-4)) return kTRUE;
3143 return kFALSE;
3144}
3145//---------------------------------------------------------------------------
3146Bool_t AliAnalysisTaskITSTrackingCheck::ConditionSPDone(AliESDtrack *track) const
3147{
3148 //
3149 // check if track has points in at least one inner barrel layer
3150 //
3151
8d63376d 3152
b7fb742d 3153 if(NPointsInnerBarrel(track)>=1) return kTRUE;
3154 return kFALSE;
3155}
3156//---------------------------------------------------------------------------
3157Bool_t AliAnalysisTaskITSTrackingCheck::ConditionSPDtwo(AliESDtrack *track) const
3158{
3159 //
3160 // check if track has points in at least two inner barrel layers
3161 //
3162
3163
3164 if(NPointsInnerBarrel(track)>=2) return kTRUE;
3165 return kFALSE;
3166}
3167//---------------------------------------------------------------------------
3168Int_t AliAnalysisTaskITSTrackingCheck::NPointsInnerBarrel(AliESDtrack *track) const
3169{
3170 //
3171 // N of points in inner barrel layers
3172 //
3173 Int_t n=0;
3174
3175 for(Int_t i=0; i<(fNITSLayers-4); i++) {
3176 if(track->HasPointOnITSLayer(i)) n++;
3177 }
3178
3179 return n;
3180}