2 // Implementation file for implementation of data analysis aft 900 GeV
4 // Author: A. Pulvirenti
12 #include "TParticle.h"
14 #include "TLorentzVector.h"
17 #include "AliESDpid.h"
18 #include "AliESDEvent.h"
19 #include "AliESDVertex.h"
20 #include "AliESDtrack.h"
22 #include "AliMCEvent.h"
23 #include "AliTOFT0maker.h"
24 #include "AliTOFcalib.h"
25 #include "AliCDBManager.h"
26 #include "AliITSPIDResponse.h"
27 #include "AliRsnMonitorTrack.h"
28 #include "AliRsnDaughter.h"
29 #include "AliRsnEvent.h"
31 #include "AliRsnAnalysisMonitorTask.h"
33 ClassImp(AliRsnAnalysisMonitorTask)
35 //__________________________________________________________________________________________________
36 AliRsnAnalysisMonitorTask::AliRsnAnalysisMonitorTask(const char *name) :
37 AliAnalysisTaskSE(name),
43 fTOFcalibrateESD(kFALSE),
44 fTOFcorrectTExp(kFALSE),
48 fEventCuts("eventCuts", AliRsnCut::kEvent),
49 fTrackCuts("trackCuts", AliRsnCut::kDaughter)
56 DefineOutput(1, TTree::Class());
59 //__________________________________________________________________________________________________
60 AliRsnAnalysisMonitorTask::AliRsnAnalysisMonitorTask(const AliRsnAnalysisMonitorTask& copy) :
61 AliAnalysisTaskSE(copy),
67 fTOFcalibrateESD(kFALSE),
68 fTOFcorrectTExp(kFALSE),
72 fEventCuts(copy.fEventCuts),
73 fTrackCuts(copy.fTrackCuts)
80 //__________________________________________________________________________________________________
81 AliRsnAnalysisMonitorTask& AliRsnAnalysisMonitorTask::operator=(const AliRsnAnalysisMonitorTask& copy)
84 // Assignment operator
87 fTOFcalibrateESD = copy.fTOFcalibrateESD;
88 fTOFcorrectTExp = copy.fTOFcorrectTExp;
89 fTOFuseT0 = copy.fTOFuseT0;
90 fTOFtuneMC = copy.fTOFtuneMC;
91 fTOFresolution = copy.fTOFresolution;
96 //__________________________________________________________________________________________________
97 AliRsnAnalysisMonitorTask::~AliRsnAnalysisMonitorTask()
103 if (fOut) delete fOut;
104 if (fESDpid) delete fESDpid;
105 if (fTrack) delete fTrack;
108 //__________________________________________________________________________________________________
109 void AliRsnAnalysisMonitorTask::UserCreateOutputObjects()
112 // Create the output data container
115 // setup TPC response
116 fESDpid = new AliESDpid;
117 fESDpid->GetTPCResponse().SetBetheBlochParameters(fTPCpar[0], fTPCpar[1], fTPCpar[2], fTPCpar[3], fTPCpar[4]);
119 // setup TOF maker & calibration
120 fTOFcalib = new AliTOFcalib;
121 fTOFmaker = new AliTOFT0maker(fESDpid, fTOFcalib);
122 fTOFmaker->SetTimeResolution(fTOFresolution);
124 // create output branch object
125 fTrack = new AliRsnMonitorTrack;
127 // create output tree
129 fOut = new TTree("rsnMonitor", "Informations on single tracks for cut checking");
130 fOut->Branch("tracks", "AliRsnMonitorTrack", &fTrack);
133 //__________________________________________________________________________________________________
134 void AliRsnAnalysisMonitorTask::UserExec(Option_t *)
137 // Main execution function.
138 // Fills the fHEvents data member with the following legenda:
139 // 0 -- event OK, prim vertex with tracks
140 // 1 -- event OK, prim vertex with SPD
141 // 2 -- event OK but vz large
145 // retrieve ESD event and related stack (if available)
146 AliESDEvent *esd = dynamic_cast<AliESDEvent*>(fInputEvent);
147 AliStack *stack = 0x0;
151 if (fMCEvent) stack = fMCEvent->Stack();
153 // create interface objects to AliRsnEvent to check event cuts
156 event.SetRefMC(fMCEvent);
157 if (!fEventCuts.IsSelected(&event)) return;
160 Int_t type = EventEval(esd);
162 // if processable, then process it
163 if (type == 0) ProcessESD(esd, esd->GetPrimaryVertexTracks(), stack);
164 else if (type == 1) ProcessESD(esd, esd->GetPrimaryVertexSPD() , stack);
167 // update histogram container
171 //__________________________________________________________________________________________________
172 void AliRsnAnalysisMonitorTask::Terminate(Option_t *)
179 //__________________________________________________________________________________________________
180 Int_t AliRsnAnalysisMonitorTask::EventEval(AliESDEvent *esd)
183 // Checks if the event is good for analysis.
185 // ---> 0 if a good primary vertex with tracks was found,
186 // ---> 1 if a good SPD primary vertex was found
187 // ---> 2 otherwise (event to be rejected)
188 // In any case, adds an entry to the TTree, to keep trace of all events.
191 // get the best primary vertex:
192 // first try that with tracks, then the SPD one
193 const AliESDVertex *vTrk = esd->GetPrimaryVertexTracks();
194 const AliESDVertex *vSPD = esd->GetPrimaryVertexSPD();
195 if (vTrk->GetNContributors() > 0) {
197 } else if (vSPD->GetNContributors() > 0) {
204 //__________________________________________________________________________________________________
205 void AliRsnAnalysisMonitorTask::ProcessESD
206 (AliESDEvent *esd, const AliESDVertex *v, AliStack *stack)
209 // Process the ESD container, to read all tracks and copy their useful values.
210 // All info are stored into an AliRsnMonitorTrack object and saved into the
211 // TClonesArray which is one of the branches of the output TTree.
214 // create interfacr objects
216 AliRsnDaughter daughter;
218 event.SetRefMC(fMCEvent);
221 // create the response function and initialize it to MC or not
222 // depending if the AliStack object is there or not
223 Bool_t isMC = (stack != 0x0);
224 AliITSPIDResponse itsrsp(isMC);
226 // TOF stuff #1: init OCDB
227 Int_t run = esd->GetRunNumber();
228 AliCDBManager *cdb = AliCDBManager::Instance();
229 cdb->SetDefaultStorage("raw://");
231 // TOF stuff #2: init calibration
232 fTOFcalib->SetCorrectTExp(fTOFcorrectTExp);
234 // TOF stuff #3: calibrate
235 if (fTOFcalibrateESD) fTOFcalib->CalibrateESD(esd);
236 if (fTOFtuneMC) fTOFmaker->TuneForMC(esd);
238 fTOFmaker->ComputeT0TOF(esd);
239 fTOFmaker->ApplyT0TOF(esd);
240 fESDpid->MakePID(esd, kFALSE, 0.);
243 // loop on all tracks
244 Int_t i, k, nITS, ntracks = esd->GetNumberOfTracks();;
245 Bool_t isTPC, isITSSA, isTOF;
246 Float_t b[2], bCov[3];
249 for (i = 0; i < ntracks; i++) {
250 AliESDtrack *track = esd->GetTrack(i);
251 event.SetDaughter(daughter, i, AliRsnDaughter::kTrack);
253 // reset the output object
254 // 'usable' flag will need to be set to 'ok'
258 fTrack->CutsPassed() = fTrackCuts.IsSelected(&daughter);
260 // skip NULL pointers, kink daughters and tracks which
261 // cannot be propagated to primary vertex
262 if (!track) continue;
263 if ((Int_t)track->GetKinkIndex(0) > 0) continue;
264 if (!track->RelateToVertex(v, esd->GetMagneticField(), kVeryBig)) continue;
266 // get MC info if possible
267 if (stack) fTrack->AdoptMC(TMath::Abs(track->GetLabel()), stack);
270 fTrack->Status() = (UInt_t)track->GetStatus();
271 fTrack->Length() = (Double_t)track->GetIntegratedLength();
272 fTrack->Charge() = (Int_t)track->Charge();
273 fTrack->PrecX() = (Double_t)track->Px();
274 fTrack->PrecY() = (Double_t)track->Py();
275 fTrack->PrecZ() = (Double_t)track->Pz();
277 // evaluate some flags from the status to decide what to do next in some points
278 isTPC = ((fTrack->Status() & AliESDtrack::kTPCin) != 0);
279 isITSSA = ((fTrack->Status() & AliESDtrack::kTPCin) == 0 && (fTrack->Status() & AliESDtrack::kITSrefit) != 0 && (fTrack->Status() & AliESDtrack::kITSpureSA) == 0 && (fTrack->Status() & AliESDtrack::kITSpid) != 0);
280 isTOF = ((fTrack->Status() & AliESDtrack::kTOFout) != 0 && (fTrack->Status() & AliESDtrack::kTIME) != 0);
282 // accept only tracks which are TPC+ITS or ITS standalone
284 fTrack->ITSsa() = kTRUE;
285 fTrack->TOFok() = kFALSE;
287 fTrack->ITSsa() = kFALSE;
288 fTrack->TOFok() = isTOF;
292 // get DCA to primary vertex
293 track->GetImpactParameters(b, bCov);
294 fTrack->DCAr() = (Double_t)b[0];
295 fTrack->DCAz() = (Double_t)b[1];
298 for (k = 0; k < 6; k++) {
299 fTrack->ITSmap(k) = track->HasPointOnITSLayer(k);
302 fTrack->ITSchi2() = track->GetITSchi2();
303 fTrack->ITSsignal() = track->GetITSsignal();
304 nITS = fTrack->SSDcount() + fTrack->SDDcount();
305 for (k = 0; k < AliPID::kSPECIES; k++) {
306 fTrack->ITSnsigma(k) = itsrsp.GetNumberOfSigmas(fTrack->Prec(), fTrack->ITSsignal(), (AliPID::EParticleType)k, nITS, kTRUE);
312 fTrack->TPCcount() = (Int_t)track->GetTPCclusters(0);
313 fTrack->TPCchi2() = (Double_t)track->GetTPCchi2();
314 fTrack->TPCsignal() = (Double_t)track->GetTPCsignal();
315 fTrack->PtpcX() = fTrack->PtpcY() = fTrack->PtpcZ() = 1E10;
316 if (track->GetInnerParam()) {
317 fTrack->PtpcX() = track->GetInnerParam()->Px();
318 fTrack->PtpcY() = track->GetInnerParam()->Py();
319 fTrack->PtpcZ() = track->GetInnerParam()->Pz();
320 for (k = 0; k < AliPID::kSPECIES; k++) {
321 fTrack->TPCnsigma(k) = fESDpid->NumberOfSigmasTPC(track, (AliPID::EParticleType)k);
328 track->GetIntegratedTimes(time);
329 fTrack->TOFsignal() = (Double_t)track->GetTOFsignal();
330 for (k = 0; k < AliPID::kSPECIES; k++) {
331 fTrack->TOFref(k) = time[k];
332 fTrack->TOFsigma(k) = (Double_t)fTOFmaker->GetExpectedSigma(fTrack->Prec(), time[k], AliPID::ParticleMass(k));
336 // add entry to TTree