Fixing coverity 17819
[u/mrichter/AliRoot.git] / MUON / AliMUONRecoCheck.cxx
CommitLineData
b8dc484b 1/**************************************************************************
48217459 2* Copyright(c) 1998-1999, 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**************************************************************************/
b8dc484b 15
d19b6003 16/* $Id$ */
17
3d1463c8 18//-----------------------------------------------------------------------------
5398f946 19/// \class AliMUONRecoCheck
78649106 20/// Utility class to check reconstruction
5398f946 21/// Reconstructed tracks are compared to reference tracks.
22/// The reference tracks are built from AliTrackReference for the
23/// hit in chamber (0..9) and from kinematics for the vertex parameters.
3d1463c8 24//-----------------------------------------------------------------------------
b8dc484b 25
9bdda5f6 26// 13 Nov 2007:
27// Added a method to create a list of reconstructed AliMUONTrack objects from
28// ESD data. This is necessary since the track objects that are actually created
29// during offline reconstruction are no longer stored to disk.
30// - Artur Szostak <artursz@iafrica.com>
103e6575 31// 25 Jan 2008:
32// Use the new ESDInterface to create MUON objects from ESD data
33// - Philippe Pillot
7c29acf1 34//
9bdda5f6 35
cc056de1 36#include "AliMUONRecoCheck.h"
37#include "AliMUONTrack.h"
103e6575 38#include "AliMUONVTrackStore.h"
39#include "AliMUONVCluster.h"
40#include "AliMUONVClusterStore.h"
cc056de1 41#include "AliMUONConstants.h"
103e6575 42#include "AliMUONESDInterface.h"
933daf4c 43#include "AliMUONTrackParam.h"
9f164762 44#include "AliMUONTriggerTrack.h"
45#include "AliMUONVTriggerTrackStore.h"
b7474abb 46#include "AliMUONTriggerStoreV1.h"
47#include "AliMUONLocalTrigger.h"
96ebe67e 48#include "AliMCEventHandler.h"
49#include "AliMCEvent.h"
b8dc484b 50#include "AliStack.h"
cc056de1 51#include "AliTrackReference.h"
9bdda5f6 52#include "AliLog.h"
53#include "AliESDEvent.h"
54#include "AliESDMuonTrack.h"
b7474abb 55#include "AliMUONDigitStoreV2S.h"
56#include "AliMUONDigit.h"
57#include "AliMpVSegmentation.h"
58#include "AliMpSegmentation.h"
59#include "AliMpPad.h"
96ebe67e 60
9f164762 61#include "AliGeomManager.h"
b7474abb 62#include "AliCDBManager.h"
9f164762 63#include "AliMpCDB.h"
64#include "AliMpDDLStore.h"
65#include "AliMUONCDB.h"
66#include "AliMUONGeometryTransformer.h"
67#include "AliMUONTriggerCircuit.h"
68#include "AliMUONVTrackReconstructor.h"
69#include "AliMUONVTriggerStore.h"
b7474abb 70#include "AliMUONCalibrationData.h"
71#include "AliMUONTriggerElectronics.h"
72
9f164762 73
74#include "TGeoManager.h"
75
61fed964 76#include <TFile.h>
77#include <TTree.h>
cc056de1 78#include <TParticle.h>
5b1afea0 79#include <TParticlePDG.h>
96ebe67e 80#include <Riostream.h>
81
9f164762 82#include "AliMUONRecoCheck.h"
83
5398f946 84/// \cond CLASSIMP
b8dc484b 85ClassImp(AliMUONRecoCheck)
5398f946 86/// \endcond
b8dc484b 87
88//_____________________________________________________________________________
1f2d22de 89AliMUONRecoCheck::AliMUONRecoCheck(const Char_t *esdFileName, const Char_t *pathSim)
48217459 90: TObject(),
96ebe67e 91fMCEventHandler(new AliMCEventHandler()),
61fed964 92fESDEvent(new AliESDEvent()),
93fESDTree (0x0),
94fESDFile (0x0),
96ebe67e 95fCurrentEvent(0),
96fTrackRefStore(0x0),
97fRecoTrackRefStore(0x0),
9f164762 98fRecoTriggerRefStore(0x0),
7c29acf1 99fRecoTrackStore(0x0),
9f164762 100fRecoTriggerTrackStore(0x0),
101fGeometryTransformer(0x0),
102fTriggerCircuit(0x0),
b7474abb 103fCalibrationData(0x0),
104fTriggerElectronics(0x0),
7c29acf1 105fESDEventOwner(kTRUE)
068e4c36 106{
48217459 107 /// Normal ctor
61fed964 108
109 // TrackRefs and Particules
96ebe67e 110 fMCEventHandler->SetInputPath(pathSim);
111 fMCEventHandler->InitIO("");
61fed964 112
113 // ESD MUON Tracks
114 fESDFile = TFile::Open(esdFileName); // open the file
115 if (!fESDFile || !fESDFile->IsOpen()) {
116 AliError(Form("opening ESD file %s failed", esdFileName));
117 fESDFile = 0x0;
118 return;
119 }
120 fESDTree = (TTree*) fESDFile->Get("esdTree"); // get the tree
121 if (!fESDTree) {
122 AliError("no ESD tree found");
123 fESDFile->Close();
124 fESDFile = 0x0;
125 return;
126 }
127 fESDEvent->ReadFromTree(fESDTree); // link fESDEvent to the tree
068e4c36 128}
129
130//_____________________________________________________________________________
7c29acf1 131AliMUONRecoCheck::AliMUONRecoCheck(AliESDEvent *esdEvent, AliMCEventHandler *mcEventHandler)
132: TObject(),
133fMCEventHandler(0),
134fESDEvent(0),
135fESDTree (0x0),
136fESDFile (0x0),
137fCurrentEvent(0),
138fTrackRefStore(0x0),
139fRecoTrackRefStore(0x0),
9f164762 140fRecoTriggerRefStore(0x0),
7c29acf1 141fRecoTrackStore(0x0),
9f164762 142fRecoTriggerTrackStore(0x0),
143fGeometryTransformer(0x0),
144fTriggerCircuit(0x0),
b7474abb 145fCalibrationData(0x0),
146fTriggerElectronics(0x0),
7c29acf1 147fESDEventOwner(kFALSE)
148{
149 /// Normal ctor
150
151 // TrackRefs and Particules
152 fMCEventHandler = mcEventHandler;
153
154 // ESD MUON Tracks
155 fESDEvent = esdEvent;
7c29acf1 156}
157
158//_____________________________________________________________________________
48217459 159AliMUONRecoCheck::~AliMUONRecoCheck()
b8dc484b 160{
48217459 161 /// Destructor
7c29acf1 162 if (fESDEventOwner) {
163 delete fMCEventHandler;
164 delete fESDEvent;
165 if (fESDFile) fESDFile->Close();
166 }
96ebe67e 167 ResetStores();
9f164762 168 delete fGeometryTransformer;
169 delete fTriggerCircuit;
b7474abb 170 delete fTriggerElectronics;
171 delete fCalibrationData;
96ebe67e 172}
173
174//_____________________________________________________________________________
175void AliMUONRecoCheck::ResetStores()
176{
177 /// Deletes all the store objects that have been created and resets the pointers to 0x0
178 delete fTrackRefStore; fTrackRefStore = 0x0;
179 delete fRecoTrackRefStore; fRecoTrackRefStore = 0x0;
9f164762 180 delete fRecoTriggerRefStore; fRecoTriggerRefStore = 0x0;
96ebe67e 181 delete fRecoTrackStore; fRecoTrackStore = 0x0;
9f164762 182 delete fRecoTriggerTrackStore; fRecoTriggerTrackStore = 0x0;
96ebe67e 183}
184
185//_____________________________________________________________________________
9f164762 186Bool_t AliMUONRecoCheck::InitCircuit()
187{
188
189 if ( fTriggerCircuit ) return kTRUE;
b7474abb 190
191 if ( ! InitGeometryTransformer() ) return kFALSE;
9f164762 192
193 fTriggerCircuit = new AliMUONTriggerCircuit(fGeometryTransformer);
194
195 // reset tracker for local trigger to trigger track conversion
196 if ( ! AliMUONESDInterface::GetTracker() )
197 AliMUONESDInterface::ResetTracker();
198
199 return kTRUE;
200}
201
202
203//_____________________________________________________________________________
a99c3449 204Int_t AliMUONRecoCheck::GetRunNumber()
205{
206 /// Return the run number of the current ESD event
207
9f164762 208 if (fESDEventOwner && fRecoTrackStore == 0x0 && fRecoTriggerTrackStore == 0x0) {
a99c3449 209 if (!fESDTree || fESDTree->GetEvent(fCurrentEvent) <= 0) {
210 AliError(Form("fails to read ESD object for event %d: cannot get the run number",fCurrentEvent));
211 return -1;
212 }
213 }
214
215 return fESDEvent->GetRunNumber();
216}
217
218//_____________________________________________________________________________
96ebe67e 219Int_t AliMUONRecoCheck::NumberOfEvents() const
220{
221 /// Return the number of events
7c29acf1 222 if (fESDEventOwner && fESDTree) return fESDTree->GetEntries();
96ebe67e 223 return 0;
b8dc484b 224}
225
48217459 226//_____________________________________________________________________________
a99c3449 227AliMUONVTrackStore* AliMUONRecoCheck::ReconstructedTracks(Int_t event, Bool_t refit)
48217459 228{
61fed964 229 /// Return a track store containing the reconstructed tracks (converted into
a99c3449 230 /// MUONTrack objects) for a given event.
231 /// Track parameters at each clusters are computed or not depending on the flag "refit".
232 /// If not, only the track parameters at first cluster are valid.
7c29acf1 233
234 if (!fESDEventOwner) {
a99c3449 235 if (fRecoTrackStore == 0x0) MakeReconstructedTracks(refit);
7c29acf1 236 return fRecoTrackStore;
237 }
238
96ebe67e 239 if (event != fCurrentEvent) {
240 ResetStores();
241 fCurrentEvent = event;
242 }
243
244 if (fRecoTrackStore != 0x0) return fRecoTrackStore;
61fed964 245 else {
246 if (!fESDTree) return 0x0;
247 if (fESDTree->GetEvent(event) <= 0) {
248 AliError(Form("fails to read ESD object for event %d", event));
249 return 0x0;
9bdda5f6 250 }
a99c3449 251 MakeReconstructedTracks(refit);
61fed964 252 return fRecoTrackStore;
9bdda5f6 253 }
9bdda5f6 254}
255
9f164762 256
257//_____________________________________________________________________________
258AliMUONVTriggerTrackStore* AliMUONRecoCheck::TriggeredTracks(Int_t event)
259{
260 /// Return a track store containing the reconstructed trigger tracks (converted into
261 /// MUONTriggerTrack objects) for a given event.
262
263 if (!fESDEventOwner) {
264 if (fRecoTriggerTrackStore == 0x0) MakeTriggeredTracks();
265 return fRecoTriggerTrackStore;
266 }
267
268 if (event != fCurrentEvent) {
269 ResetStores();
270 fCurrentEvent = event;
271 }
272
273 if (fRecoTriggerTrackStore != 0x0) return fRecoTriggerTrackStore;
274 else {
275 if (!fESDTree) return 0x0;
276 if (fESDTree->GetEvent(event) <= 0) {
277 AliError(Form("fails to read ESD object for event %d", event));
278 return 0x0;
279 }
280 MakeTriggeredTracks();
281 return fRecoTriggerTrackStore;
282 }
283}
284
285
9bdda5f6 286//_____________________________________________________________________________
96ebe67e 287AliMUONVTrackStore* AliMUONRecoCheck::TrackRefs(Int_t event)
b8dc484b 288{
48217459 289 /// Return a track store containing the track references (converted into
290 /// MUONTrack objects) for a given event
7c29acf1 291
292 if (!fESDEventOwner) {
a99c3449 293 if (fTrackRefStore == 0x0) MakeTrackRefs();
7c29acf1 294 return fTrackRefStore;
295 }
296
96ebe67e 297 if (event != fCurrentEvent) {
298 ResetStores();
299 fCurrentEvent = event;
300 }
301
302 if (fTrackRefStore != 0x0) return fTrackRefStore;
303 else {
e473e8b0 304 if (!fMCEventHandler->LoadEvent(event)) {
61fed964 305 AliError(Form("fails to read MC objects for event %d", event));
306 return 0x0;
307 }
96ebe67e 308 MakeTrackRefs();
309 return fTrackRefStore;
310 }
b8dc484b 311}
312
313//_____________________________________________________________________________
9f164762 314AliMUONVTriggerTrackStore* AliMUONRecoCheck::TriggerableTracks(Int_t event)
315{
316 /// Return a trigger track store containing the triggerable track references (converted into
317 /// AliMUONTriggerTrack objects) for a given event
318
319 if (!fESDEventOwner) {
320 if (fRecoTriggerRefStore == 0x0) MakeTriggerableTracks();
321 return fRecoTriggerRefStore;
322 }
323
324 if (event != fCurrentEvent) {
325 ResetStores();
326 fCurrentEvent = event;
327 }
328
329 if (fRecoTriggerRefStore != 0x0) return fRecoTriggerRefStore;
330 else {
e473e8b0 331 if (!fMCEventHandler->LoadEvent(event)) {
9f164762 332 AliError(Form("fails to read MC objects for event %d", event));
333 return 0x0;
334 }
335 MakeTriggerableTracks();
336 return fRecoTriggerRefStore;
337 }
338}
339
340
341//_____________________________________________________________________________
f202486b 342AliMUONVTrackStore* AliMUONRecoCheck::ReconstructibleTracks(Int_t event, UInt_t requestedStationMask, Bool_t request2ChInSameSt45)
b8dc484b 343{
f202486b 344 /// Return a track store containing the reconstructible tracks for a given event,
345 /// according to the mask of requested stations and the minimum number of chambers hit in stations 4 & 5.
7c29acf1 346
347 if (!fESDEventOwner) {
a99c3449 348 if (fRecoTrackRefStore == 0x0) {
349 if (TrackRefs(event) == 0x0) return 0x0;
f202486b 350 MakeReconstructibleTracks(requestedStationMask, request2ChInSameSt45);
a99c3449 351 }
7c29acf1 352 return fRecoTrackRefStore;
353 }
354
96ebe67e 355 if (event != fCurrentEvent) {
356 ResetStores();
357 fCurrentEvent = event;
358 }
359
360 if (fRecoTrackRefStore != 0x0) return fRecoTrackRefStore;
361 else {
362 if (TrackRefs(event) == 0x0) return 0x0;
f202486b 363 MakeReconstructibleTracks(requestedStationMask, request2ChInSameSt45);
96ebe67e 364 return fRecoTrackRefStore;
365 }
48217459 366}
b8dc484b 367
9f164762 368
48217459 369//_____________________________________________________________________________
a99c3449 370void AliMUONRecoCheck::MakeReconstructedTracks(Bool_t refit)
61fed964 371{
372 /// Make reconstructed tracks
103e6575 373 if (!(fRecoTrackStore = AliMUONESDInterface::NewTrackStore())) return;
61fed964 374
b1fea02e 375 // loop over all reconstructed tracks and add them to the store (skip ghosts)
103e6575 376 Int_t nTracks = (Int_t) fESDEvent->GetNumberOfMuonTracks();
b1fea02e 377 for (Int_t iTrack = 0; iTrack < nTracks; iTrack++) {
378 AliESDMuonTrack* esdTrack = fESDEvent->GetMuonTrack(iTrack);
a99c3449 379 if (esdTrack->ContainTrackerData()) AliMUONESDInterface::Add(*esdTrack, *fRecoTrackStore, refit);
b1fea02e 380 }
61fed964 381
382}
383
9f164762 384
385//_____________________________________________________________________________
386void AliMUONRecoCheck::MakeTriggeredTracks()
387{
388 /// Make reconstructed trigger tracks
389 if (!(fRecoTriggerTrackStore = AliMUONESDInterface::NewTriggerTrackStore())) return;
390
391 AliMUONVTriggerStore* tmpTriggerStore = AliMUONESDInterface::NewTriggerStore();
392 if ( ! tmpTriggerStore ) return;
393
394 // loop over all reconstructed tracks and add them to the store (include ghosts)
395 Int_t nTracks = (Int_t) fESDEvent->GetNumberOfMuonTracks();
396 for (Int_t iTrack = 0; iTrack < nTracks; iTrack++) {
397 AliESDMuonTrack* esdTrack = fESDEvent->GetMuonTrack(iTrack);
398 if (esdTrack->ContainTriggerData()) AliMUONESDInterface::Add(*esdTrack, *tmpTriggerStore);
399 }
400
b7474abb 401 if ( InitCircuit() ) {
9f164762 402
b7474abb 403 AliMUONVTrackReconstructor* tracker = AliMUONESDInterface::GetTracker();
404 tracker->EventReconstructTrigger(*fTriggerCircuit, *tmpTriggerStore, *fRecoTriggerTrackStore);
405 }
9f164762 406
407 delete tmpTriggerStore;
408}
409
410//_____________________________________________________________________________
411void AliMUONRecoCheck::TriggerToTrack(const AliMUONLocalTrigger& locTrg, AliMUONTriggerTrack& triggerTrack)
412{
413 /// Make trigger track from local trigger info
414 if ( ! InitCircuit() ) return;
415 AliMUONVTrackReconstructor* tracker = AliMUONESDInterface::GetTracker();
416 tracker->TriggerToTrack(*fTriggerCircuit, locTrg, triggerTrack);
417}
418
419
61fed964 420//_____________________________________________________________________________
96ebe67e 421void AliMUONRecoCheck::MakeTrackRefs()
48217459 422{
423 /// Make reconstructible tracks
103e6575 424 AliMUONVTrackStore *tmpTrackRefStore = AliMUONESDInterface::NewTrackStore();
425 if (!tmpTrackRefStore) return;
48217459 426
96ebe67e 427 Double_t x, y, z, pX, pY, pZ, bendingSlope, nonBendingSlope, inverseBendingMomentum;
428 TParticle* particle;
429 TClonesArray* trackRefs;
430 Int_t nTrackRef = fMCEventHandler->MCEvent()->GetNumberOfTracks();
103e6575 431 AliMUONVClusterStore* cStore = AliMUONESDInterface::NewClusterStore();
432 if (!cStore) return;
433 AliMUONVCluster* hit = cStore->CreateCluster(0,0,0);
b8dc484b 434
96ebe67e 435 // loop over simulated tracks
436 for (Int_t iTrackRef = 0; iTrackRef < nTrackRef; ++iTrackRef) {
437 Int_t nHits = fMCEventHandler->GetParticleAndTR(iTrackRef, particle, trackRefs);
48217459 438
96ebe67e 439 // skip empty trackRefs
440 if (nHits < 1) continue;
441
61fed964 442 // get the particle charge for further calculation
443 TParticlePDG* ppdg = particle->GetPDG();
46df197d 444 Int_t charge = ppdg != NULL ? (Int_t)(ppdg->Charge()/3.0) : 0;
96ebe67e 445
61fed964 446 AliMUONTrack track;
48217459 447
96ebe67e 448 // loop over simulated track hits
449 for (Int_t iHit = 0; iHit < nHits; ++iHit) {
450 AliTrackReference* trackReference = static_cast<AliTrackReference*>(trackRefs->UncheckedAt(iHit));
068e4c36 451
61fed964 452 // skip trackRefs not in MUON
453 if (trackReference->DetectorId() != AliTrackReference::kMUON) continue;
454
96ebe67e 455 // Get track parameters of current hit
456 x = trackReference->X();
457 y = trackReference->Y();
458 z = trackReference->Z();
459 pX = trackReference->Px();
460 pY = trackReference->Py();
461 pZ = trackReference->Pz();
b8dc484b 462
96ebe67e 463 // check chamberId of current trackReference
c59f70b9 464 Int_t detElemId = trackReference->UserId();
465 Int_t chamberId = detElemId / 100 - 1;
96ebe67e 466 if (chamberId < 0 || chamberId >= AliMUONConstants::NTrackingCh()) continue;
48217459 467
96ebe67e 468 // set hit parameters
c59f70b9 469 hit->SetUniqueID(AliMUONVCluster::BuildUniqueID(chamberId, detElemId, iHit));
103e6575 470 hit->SetXYZ(x,y,z);
471 hit->SetErrXY(0.,0.);
b8dc484b 472
96ebe67e 473 // compute track parameters at hit
7cdc832b 474 bendingSlope = 0;
475 nonBendingSlope = 0;
476 inverseBendingMomentum = 0;
96ebe67e 477 if (TMath::Abs(pZ) > 0) {
478 bendingSlope = pY/pZ;
479 nonBendingSlope = pX/pZ;
b8dc484b 480 }
96ebe67e 481 Double_t pYZ = TMath::Sqrt(pY*pY+pZ*pZ);
482 if (pYZ >0) inverseBendingMomentum = 1/pYZ;
483 inverseBendingMomentum *= charge;
484
485 // set track parameters at hit
486 AliMUONTrackParam trackParam;
487 trackParam.SetNonBendingCoor(x);
488 trackParam.SetBendingCoor(y);
489 trackParam.SetZ(z);
490 trackParam.SetBendingSlope(bendingSlope);
491 trackParam.SetNonBendingSlope(nonBendingSlope);
492 trackParam.SetInverseBendingMomentum(inverseBendingMomentum);
493
494 // add track parameters at current hit to the track
103e6575 495 track.AddTrackParamAtCluster(trackParam, *hit, kTRUE);
96ebe67e 496 }
497
61fed964 498 // if none of the track hits was in MUON, goto the next track
499 if (track.GetNClusters() < 1) continue;
500
501 // get track parameters at particle's vertex
502 x = particle->Vx();
503 y = particle->Vy();
504 z = particle->Vz();
505 pX = particle->Px();
506 pY = particle->Py();
507 pZ = particle->Pz();
508
509 // compute rest of track parameters at particle's vertex
510 bendingSlope = 0;
511 nonBendingSlope = 0;
512 inverseBendingMomentum = 0;
513 if (TMath::Abs(pZ) > 0) {
514 bendingSlope = pY/pZ;
515 nonBendingSlope = pX/pZ;
516 }
517 Double_t pYZ = TMath::Sqrt(pY*pY+pZ*pZ);
518 if (pYZ >0) inverseBendingMomentum = 1/pYZ;
519 inverseBendingMomentum *= charge;
520
521 // set track parameters at particle's vertex
522 AliMUONTrackParam trackParamAtVertex;
523 trackParamAtVertex.SetNonBendingCoor(x);
524 trackParamAtVertex.SetBendingCoor(y);
525 trackParamAtVertex.SetZ(z);
526 trackParamAtVertex.SetBendingSlope(bendingSlope);
527 trackParamAtVertex.SetNonBendingSlope(nonBendingSlope);
528 trackParamAtVertex.SetInverseBendingMomentum(inverseBendingMomentum);
529
530 // add track parameters at vertex
531 track.SetTrackParamAtVertex(&trackParamAtVertex);
532
103e6575 533 // store the track
2e2d0c44 534 track.SetUniqueID(iTrackRef);
96ebe67e 535 tmpTrackRefStore->Add(track);
b8dc484b 536 }
b8dc484b 537
96ebe67e 538 CleanMuonTrackRef(tmpTrackRefStore);
48217459 539
103e6575 540 delete hit;
541 delete cStore;
96ebe67e 542 delete tmpTrackRefStore;
b8dc484b 543}
544
545//_____________________________________________________________________________
9f164762 546void AliMUONRecoCheck::MakeTriggerableTracks()
547{
548 /// Make triggerable tracks
549 if (!(fRecoTriggerRefStore = AliMUONESDInterface::NewTriggerTrackStore()))
550 return;
551
b7474abb 552 Double_t x, y, z, slopeX, slopeY, pZ, xLoc, yLoc, zLoc;
9f164762 553 TParticle* particle;
554 TClonesArray* trackRefs;
555 Int_t nTrackRef = fMCEventHandler->MCEvent()->GetNumberOfTracks();
556
557 // loop over simulated tracks
558 for (Int_t iTrackRef = 0; iTrackRef < nTrackRef; ++iTrackRef) {
559 Int_t nHits = fMCEventHandler->GetParticleAndTR(iTrackRef, particle, trackRefs);
560
561 // skip empty trackRefs
562 if (nHits < 1) continue;
563
564 AliMUONTriggerTrack track;
b7474abb 565 AliMUONDigitStoreV2S digitStore;
9f164762 566 Int_t hitsOnTrigger = 0;
567 Int_t currCh = -1;
568
569 // loop over simulated track hits
570 for (Int_t iHit = 0; iHit < nHits; ++iHit) {
571 AliTrackReference* trackReference = static_cast<AliTrackReference*>(trackRefs->UncheckedAt(iHit));
572
573 // skip trackRefs not in MUON
574 if (trackReference->DetectorId() != AliTrackReference::kMUON) continue;
575
576 // check chamberId of current trackReference
577 Int_t detElemId = trackReference->UserId();
578 Int_t chamberId = detElemId / 100 - 1;
579 if (chamberId < AliMUONConstants::NTrackingCh() || chamberId >= AliMUONConstants::NCh() ) continue;
b7474abb 580
581 // Get track parameters of current hit
582 x = trackReference->X();
583 y = trackReference->Y();
584 z = trackReference->Z();
585
586 if ( InitTriggerResponse() ) {
587 fGeometryTransformer->Global2Local(detElemId, x, y, z, xLoc, yLoc, zLoc);
588
589 Int_t nboard = 0;
590 for ( Int_t cath = AliMp::kCath0; cath <= AliMp::kCath1; ++cath )
591 {
592 const AliMpVSegmentation* seg
593 = AliMpSegmentation::Instance()
594 ->GetMpSegmentation(detElemId,AliMp::GetCathodType(cath));
595
596 AliMpPad pad = seg->PadByPosition(xLoc,yLoc,kFALSE);
597 Int_t ix = pad.GetIx();
598 Int_t iy = pad.GetIy();
599
600 if ( !pad.IsValid() ) continue;
601
602 if ( cath == AliMp::kCath0 ) nboard = pad.GetLocalBoardId(0);
603
604 AliMUONDigit* digit = new AliMUONDigit(detElemId,nboard,
605 pad.GetLocalBoardChannel(0),cath);
606 digit->SetPadXY(ix,iy);
607 digit->SetCharge(1.);
608 digitStore.Add(*digit,AliMUONVDigitStore::kDeny);
609 }
610 }
611
9f164762 612 if ( hitsOnTrigger == 0 ) {
9f164762 613 pZ = trackReference->Pz();
614 slopeX = ( pZ == 0. ) ? 99999. : trackReference->Px() / pZ;
615 slopeY = ( pZ == 0. ) ? 99999. : trackReference->Py() / pZ;
616
617 track.SetX11(x);
618 track.SetY11(y);
619 track.SetZ11(z);
620 track.SetSlopeX(slopeX);
621 track.SetSlopeY(slopeY);
622 }
623
624 if ( currCh != chamberId ) {
625 hitsOnTrigger++;
626 currCh = chamberId;
627 }
628
629 } // loop on hits
b7474abb 630
631 if ( hitsOnTrigger < 3 ) continue;
632
633 // Check if the track passes the trigger algorithm
634 if ( InitTriggerResponse() ) {
635 AliMUONTriggerStoreV1 triggerStore;
636 fTriggerElectronics->Digits2Trigger(digitStore,triggerStore);
637
638 TIter next(triggerStore.CreateIterator());
639 AliMUONLocalTrigger* locTrg(0x0);
640
641 Int_t ptCutLevel = 0;
642
643 while ( ( locTrg = static_cast<AliMUONLocalTrigger*>(next()) ) )
644 {
645 if ( locTrg->IsTrigX() && locTrg->IsTrigY() )
646 {
647 ptCutLevel = TMath::Max(ptCutLevel, 1);
648 if ( locTrg->LoHpt() ) ptCutLevel = TMath::Max(ptCutLevel, 3);
649 else if ( locTrg->LoLpt() ) ptCutLevel = TMath::Max(ptCutLevel, 2);
650 } // board is fired
651 } // end of loop on Local Trigger
652 track.SetPtCutLevel(ptCutLevel);
9f164762 653 }
b7474abb 654 track.SetUniqueID(iTrackRef);
655 fRecoTriggerRefStore->Add(track);
9f164762 656 }
657}
658
659
660//_____________________________________________________________________________
96ebe67e 661void AliMUONRecoCheck::CleanMuonTrackRef(const AliMUONVTrackStore *tmpTrackRefStore)
b8dc484b 662{
48217459 663 /// Re-calculate hits parameters because two AliTrackReferences are recorded for
664 /// each chamber (one when particle is entering + one when particle is leaving
665 /// the sensitive volume)
103e6575 666 if (!(fTrackRefStore = AliMUONESDInterface::NewTrackStore())) return;
48217459 667
96ebe67e 668 Double_t maxGasGap = 1.; // cm
669 Double_t x, y, z, pX, pY, pZ, x1, y1, z1, pX1, pY1, pZ1, z2;
670 Double_t bendingSlope,nonBendingSlope,inverseBendingMomentum;
103e6575 671 AliMUONVClusterStore* cStore = AliMUONESDInterface::NewClusterStore();
672 if (!cStore) return;
673 AliMUONVCluster* hit = cStore->CreateCluster(0,0,0);
48217459 674
96ebe67e 675 // create iterator
676 TIter next(tmpTrackRefStore->CreateIterator());
48217459 677
96ebe67e 678 // loop over tmpTrackRef
48217459 679 AliMUONTrack* track;
96ebe67e 680 while ( ( track = static_cast<AliMUONTrack*>(next()) ) ) {
681
682 AliMUONTrack newTrack;
683
684 // loop over tmpTrackRef's hits
685 Int_t iHit1 = 0;
686 Int_t nTrackHits = track->GetNClusters();
687 while (iHit1 < nTrackHits) {
688 AliMUONTrackParam *trackParam1 = (AliMUONTrackParam*) track->GetTrackParamAtCluster()->UncheckedAt(iHit1);
689
690 // get track parameters at hit1
691 x1 = trackParam1->GetNonBendingCoor();
692 y1 = trackParam1->GetBendingCoor();
693 z1 = trackParam1->GetZ();
694 pX1 = trackParam1->Px();
695 pY1 = trackParam1->Py();
696 pZ1 = trackParam1->Pz();
697
698 // prepare new track parameters
699 x = x1;
700 y = y1;
701 z = z1;
702 pX = pX1;
703 pY = pY1;
704 pZ = pZ1;
705
706 // loop over next tmpTrackRef's hits
707 Int_t nCombinedHits = 1;
708 for (Int_t iHit2 = iHit1+1; iHit2 < nTrackHits; iHit2++) {
709 AliMUONTrackParam *trackParam2 = (AliMUONTrackParam*) track->GetTrackParamAtCluster()->UncheckedAt(iHit2);
48217459 710
96ebe67e 711 // get z position of hit2
712 z2 = trackParam2->GetZ();
713
714 // complete new track parameters if hit2 is on the same detection element
715 if ( TMath::Abs(z2-z1) < maxGasGap ) {
716 x += trackParam2->GetNonBendingCoor();
717 y += trackParam2->GetBendingCoor();
718 z += z2;
719 pX += trackParam2->Px();
720 pY += trackParam2->Py();
721 pZ += trackParam2->Pz();
722 nCombinedHits++;
48217459 723 iHit1 = iHit2;
724 }
725
96ebe67e 726 }
727
728 // finalize new track parameters
729 x /= (Double_t)nCombinedHits;
730 y /= (Double_t)nCombinedHits;
731 z /= (Double_t)nCombinedHits;
732 pX /= (Double_t)nCombinedHits;
733 pY /= (Double_t)nCombinedHits;
734 pZ /= (Double_t)nCombinedHits;
735 bendingSlope = 0;
736 nonBendingSlope = 0;
737 inverseBendingMomentum = 0;
738 if (TMath::Abs(pZ) > 0) {
739 bendingSlope = pY/pZ;
740 nonBendingSlope = pX/pZ;
741 }
742 Double_t pYZ = TMath::Sqrt(pY*pY+pZ*pZ);
743 if (pYZ >0) inverseBendingMomentum = 1/pYZ;
744 inverseBendingMomentum *= trackParam1->GetCharge();
b8dc484b 745
96ebe67e 746 // set hit parameters
c59f70b9 747 hit->SetUniqueID(trackParam1->GetClusterPtr()->GetUniqueID());
103e6575 748 hit->SetXYZ(x,y,z);
749 hit->SetErrXY(0.,0.);
48217459 750
96ebe67e 751 // set new track parameters at new hit
752 AliMUONTrackParam trackParam;
753 trackParam.SetNonBendingCoor(x);
754 trackParam.SetBendingCoor(y);
755 trackParam.SetZ(z);
756 trackParam.SetBendingSlope(bendingSlope);
757 trackParam.SetNonBendingSlope(nonBendingSlope);
758 trackParam.SetInverseBendingMomentum(inverseBendingMomentum);
759
760 // add track parameters at current hit to the track
103e6575 761 newTrack.AddTrackParamAtCluster(trackParam, *hit, kTRUE);
b8dc484b 762
763 iHit1++;
96ebe67e 764 }
48217459 765
2e2d0c44 766 newTrack.SetUniqueID(track->GetUniqueID());
96ebe67e 767 newTrack.SetTrackParamAtVertex(track->GetTrackParamAtVertex());
768 fTrackRefStore->Add(newTrack);
b8dc484b 769
96ebe67e 770 }
b8dc484b 771
103e6575 772 delete hit;
773 delete cStore;
b8dc484b 774}
775
776//_____________________________________________________________________________
f202486b 777void AliMUONRecoCheck::MakeReconstructibleTracks(UInt_t requestedStationMask, Bool_t request2ChInSameSt45)
b8dc484b 778{
61fed964 779 /// Isolate the reconstructible tracks
103e6575 780 if (!(fRecoTrackRefStore = AliMUONESDInterface::NewTrackStore())) return;
48217459 781
96ebe67e 782 // create iterator on trackRef
783 TIter next(fTrackRefStore->CreateIterator());
b8dc484b 784
f202486b 785 // loop over trackRef and add reconstructible tracks to fRecoTrackRefStore
48217459 786 AliMUONTrack* track;
96ebe67e 787 while ( ( track = static_cast<AliMUONTrack*>(next()) ) ) {
f202486b 788 if (track->IsValid(requestedStationMask, request2ChInSameSt45)) fRecoTrackRefStore->Add(*track);
789 }
790
791}
792
793//_____________________________________________________________________________
794AliMUONTrack* AliMUONRecoCheck::FindCompatibleTrack(AliMUONTrack &track, AliMUONVTrackStore &trackStore,
795 Int_t &nMatchClusters, Bool_t useLabel, Double_t sigmaCut)
796{
797 /// Return the track from the store matched with the given track (or 0x0) and the number of matched clusters.
798 /// Matching is done by using the MC label of by comparing cluster/TrackRef positions according to the flag "useLabel".
799 /// WARNING: Who match who matters since the matching algorithm uses the *fraction* of matched clusters of the given track
800
801 AliMUONTrack *matchedTrack = 0x0;
802 nMatchClusters = 0;
803
804 if (useLabel) { // by using the MC label
96ebe67e 805
f202486b 806 // get the corresponding simulated track if any
807 Int_t label = track.GetMCLabel();
808 matchedTrack = (AliMUONTrack*) trackStore.FindObject(label);
48217459 809
f202486b 810 // get the fraction of matched clusters
811 if (matchedTrack) {
812 Int_t nClusters = track.GetNClusters();
813 for (Int_t iCl = 0; iCl < nClusters; iCl++)
814 if (((AliMUONTrackParam*) track.GetTrackParamAtCluster()->UncheckedAt(iCl))->GetClusterPtr()->GetMCLabel() == label)
815 nMatchClusters++;
816 }
96ebe67e 817
f202486b 818 } else { // by comparing cluster/TrackRef positions
48217459 819
f202486b 820 // look for the corresponding simulated track if any
821 TIter next(trackStore.CreateIterator());
822 AliMUONTrack* track2;
823 while ( ( track2 = static_cast<AliMUONTrack*>(next()) ) ) {
824
825 // check compatibility
826 Int_t n = 0;
827 if (track.Match(*track2, sigmaCut, n)) {
828 matchedTrack = track2;
829 nMatchClusters = n;
830 break;
831 }
832
833 }
48217459 834
b8dc484b 835 }
f202486b 836
837 return matchedTrack;
838
b8dc484b 839}
840
9f164762 841
842//_____________________________________________________________________________
57e2ad1a 843AliMUONTriggerTrack* AliMUONRecoCheck::FindCompatibleTrack(AliMUONTriggerTrack &track, const AliMUONVTriggerTrackStore &triggerTrackStore,
9f164762 844 Double_t sigmaCut)
845{
846 /// Return the trigger track from the store matched with the given track (or 0x0).
847 /// Matching is done by comparing cluster/TrackRef positions.
848
849 AliMUONTriggerTrack *matchedTrack = 0x0;
850
851 // look for the corresponding simulated track if any
852 TIter next(triggerTrackStore.CreateIterator());
853 AliMUONTriggerTrack* track2;
854 while ( ( track2 = static_cast<AliMUONTriggerTrack*>(next()) ) ) {
9f164762 855 // check compatibility
856 if (track.Match(*track2, sigmaCut)) {
857 matchedTrack = track2;
858 break;
859 }
860 }
861
862 return matchedTrack;
863
864}
865
b7474abb 866
867//____________________________________________________________________________
868Bool_t AliMUONRecoCheck::InitTriggerResponse()
869{
870 /// Initialize trigger electronics
871 /// for building of triggerable tracks from MC
872
873 if ( fTriggerElectronics ) return kTRUE;
874
875 if ( ! InitGeometryTransformer() ) return kFALSE;
876
877 if ( ! InitCalibrationData() ) return kFALSE;
878
879 fTriggerElectronics = new AliMUONTriggerElectronics(fCalibrationData);
880
881 return kTRUE;
882}
883
884
885//____________________________________________________________________________
886Bool_t AliMUONRecoCheck::InitCalibrationData()
887{
888 /// Initialize calibration data
889 if ( ! fCalibrationData ) {
890 if ( !AliMUONCDB::CheckOCDB() ) return kFALSE;
891 fCalibrationData = new AliMUONCalibrationData(AliCDBManager::Instance()->GetRun());
892 }
893 return kTRUE;
894}
895
896
897//____________________________________________________________________________
898Bool_t AliMUONRecoCheck::InitGeometryTransformer()
899{
900 /// Return calibration data
901 /// (create it if necessary)
902 if ( ! fGeometryTransformer ) {
903
904 if ( !AliMUONCDB::CheckOCDB() ) return kFALSE;
905
906 if ( !AliGeomManager::GetGeometry() )
907 AliGeomManager::LoadGeometry();
908
909 if ( !AliMpDDLStore::Instance(false) )
910 AliMpCDB::LoadDDLStore();
911
912 fGeometryTransformer = new AliMUONGeometryTransformer();
913 fGeometryTransformer->LoadGeometryData();
914 }
915
916 return kTRUE;
917}
918