Removing old segmentation
[u/mrichter/AliRoot.git] / MUON / AliMUON.cxx
CommitLineData
4c039060 1/**************************************************************************
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 **************************************************************************/
4c039060 15
88cb7938 16/* $Id$ */
17
a9e2aefa 18
19///////////////////////////////////////////////
fe4da5cc 20// Manager and hits classes for set:MUON //
21////////////////////////////////////////////////
22
88cb7938 23#include "Riostream.h"
24
25#include <AliPDG.h>
a897a37a 26#include <TBRIK.h>
88cb7938 27#include <TCanvas.h>
28#include <TDirectory.h>
29#include <TFile.h>
94de3818 30#include <TGeometry.h>
88cb7938 31#include <TMinuit.h>
fe4da5cc 32#include <TNode.h>
88cb7938 33#include <TNtuple.h>
fe4da5cc 34#include <TObjArray.h>
88cb7938 35#include <TObject.h>
36#include <TObjectTable.h>
37#include <TPad.h>
a897a37a 38#include <TParticle.h>
39#include <TROOT.h>
88cb7938 40#include <TRandom.h>
41#include <TRotMatrix.h>
42#include <TTUBE.h>
a9e2aefa 43#include <TTUBE.h>
88cb7938 44#include <TTree.h>
45#include <TVector.h>
46#include <TVirtualMC.h>
fe4da5cc 47
d4bb94a1 48//#include "AliHeader.h"
88cb7938 49#include "AliLoader.h"
dbb4e169 50#include "AliRunDigitizer.h"
307d9d04 51#include "AliMC.h"
fac3a970 52#include "AliRun.h"
fe4da5cc 53#include "AliMUON.h"
88cb7938 54#include "AliMUONChamberTrigger.h"
88cb7938 55#include "AliMUONConstants.h"
a9e2aefa 56#include "AliMUONHit.h"
88cb7938 57#include "AliMUONMerger.h"
a9e2aefa 58#include "AliMUONPadHit.h"
a9e2aefa 59#include "AliMUONRawCluster.h"
88cb7938 60#include "AliMUONTransientDigit.h"
ecfa008b 61#include "AliMUONTriggerCircuit.h"
d4bb94a1 62#include "AliMUONGeometryBuilder.h"
e118b27e 63#include "AliMUONCommonGeometryBuilder.h"
d1cd2474 64#include "AliMUONVGeometryBuilder.h"
e118b27e 65#include "AliMUONGeometryDEIndexing.h"
a713db22 66#include "AliMUONGeometrySegmentation.h"
dbb4e169 67#include "AliMUONDigitizerv2.h"
d1775029 68#include "AliMUONSDigitizerv1.h"
39e39255 69#include "AliMUONRawData.h"
002920d1 70#include "AliMUONFactoryV2.h"
8c343c7c 71#include "AliLog.h"
a9e2aefa 72
73// Defaults parameters for Z positions of chambers
74// taken from values for "stations" in AliMUON::AliMUON
75// const Float_t zch[7]={528, 690., 975., 1249., 1449., 1610, 1710.};
76// and from array "dstation" in AliMUONv1::CreateGeometry
77// Float_t dstation[5]={20., 20., 20, 20., 20.};
78// for tracking chambers,
79// according to (Z1 = zch - dstation) and (Z2 = zch + dstation)
80// for the first and second chambers in the station, respectively,
81// and from "DTPLANES" in AliMUONv1::CreateGeometry
82// const Float_t DTPLANES = 15.;
83// for trigger chambers,
84// according to (Z1 = zch) and (Z2 = zch + DTPLANES)
85// for the first and second chambers in the station, respectively
fe4da5cc 86
fe4da5cc 87ClassImp(AliMUON)
30178c30 88
ce3f5e87 89//__________________________________________________________________
fe4da5cc 90AliMUON::AliMUON()
30178c30 91 : AliDetector(),
92 fNCh(0),
93 fNTrackingCh(0),
94 fMUONData(0),
95 fSplitLevel(0),
96 fChambers(0),
30178c30 97 fTriggerCircuits(0),
d4bb94a1 98 fGeometryBuilder(0),
002920d1 99 fSegmentationType(2),// set to 2 default wise new seg
e118b27e 100 fDEIndexing(0),
30178c30 101 fAccCut(kFALSE),
102 fAccMin(0.),
103 fAccMax(0.),
104 fMaxStepGas(0.),
105 fMaxStepAlu(0.),
106 fMaxDestepGas(0.),
107 fMaxDestepAlu(0.),
108 fMaxIterPad(0),
109 fCurIterPad(0),
8109731a 110 fMerger(0),
111 fFactory(0)
fe4da5cc 112{
de05461e 113// Default Constructor
114//
1bd26093 115 fIshunt = 0;
fe4da5cc 116}
30178c30 117
ce3f5e87 118//__________________________________________________________________
fe4da5cc 119AliMUON::AliMUON(const char *name, const char *title)
30178c30 120 : AliDetector(name,title),
121 fNCh(AliMUONConstants::NCh()),
122 fNTrackingCh(AliMUONConstants::NTrackingCh()),
123 fMUONData(0),
124 fSplitLevel(0),
125 fChambers(0),
30178c30 126 fTriggerCircuits(0),
d4bb94a1 127 fGeometryBuilder(0),
002920d1 128 fSegmentationType(2),// set to 2 default wise new seg
e118b27e 129 fDEIndexing(0),
30178c30 130 fAccCut(kFALSE),
131 fAccMin(0.),
132 fAccMax(0.),
133 fMaxStepGas(0.1),
134 fMaxStepAlu(0.1),
135 fMaxDestepGas(-1), // Negatives values are ignored by geant3 CONS200
136 fMaxDestepAlu(-1), // in the calculation of the tracking parameters
137 fMaxIterPad(0),
138 fCurIterPad(0),
8109731a 139 fMerger(0),
140 fFactory(0)
fe4da5cc 141{
f665c1ea 142
30178c30 143 fIshunt = 0;
a897a37a 144
ce3f5e87 145 SetMarkerColor(kRed);//
e118b27e 146
147 // Geometry builder
148 fGeometryBuilder = new AliMUONGeometryBuilder(this);
149
150 // Common geometry definitions
151 fGeometryBuilder
152 ->AddBuilder(new AliMUONCommonGeometryBuilder(this));
153
154 // Define the global transformation:
155 // Transformation from the old ALICE coordinate system to a new one:
156 // x->-x, z->-z
157 TGeoRotation* rotGlobal
158 = new TGeoRotation("rotGlobal", 90., 180., 90., 90., 180., 0.);
159 fGeometryBuilder
160 ->SetGlobalTransformation (TGeoCombiTrans(0., 0., 0., rotGlobal));
161
162 // Detection elements indexing
163 fDEIndexing = new AliMUONGeometryDEIndexing();
164
a9e2aefa 165//
ce3f5e87 166// Creating List of Chambers
a9e2aefa 167 Int_t ch;
f665c1ea 168 fChambers = new TObjArray(AliMUONConstants::NCh());
d1cd2474 169
a9e2aefa 170 // Loop over stations
f665c1ea 171 for (Int_t st = 0; st < AliMUONConstants::NCh() / 2; st++) {
a9e2aefa 172 // Loop over 2 chambers in the station
ce3f5e87 173 for (Int_t stCH = 0; stCH < 2; stCH++) {
174 //
175 //
176 // Default Parameters for Muon Tracking Stations
177 ch = 2 * st + stCH;
178 if (ch < AliMUONConstants::NTrackingCh()) {
179 fChambers->AddAt(new AliMUONChamber(ch),ch);
180 } else {
181 fChambers->AddAt(new AliMUONChamberTrigger(ch),ch);
182 }
183 AliMUONChamber* chamber = (AliMUONChamber*) fChambers->At(ch);
d1cd2474 184 //chamber->SetGid(0);
ce3f5e87 185 // Default values for Z of chambers
186 chamber->SetZ(AliMUONConstants::DefaultChamberZ(ch));
187 //
188 chamber->InitGeo(AliMUONConstants::DefaultChamberZ(ch));
189 // Set chamber inner and outer radius to default
f4a7360f 190 chamber->SetRInner(AliMUONConstants::Dmin(st)/2.);
191 chamber->SetROuter(AliMUONConstants::Dmax(st)/2.);
ce3f5e87 192 //
193 } // Chamber stCH (0, 1) in
a9e2aefa 194 } // Station st (0...)
ce3f5e87 195
ce3f5e87 196 // cp new design of AliMUONTriggerDecision
197 fTriggerCircuits = new TObjArray(AliMUONConstants::NTriggerCircuit());
198 for (Int_t circ=0; circ<AliMUONConstants::NTriggerCircuit(); circ++) {
199 fTriggerCircuits->AddAt(new AliMUONTriggerCircuit(),circ);
200 }
fe4da5cc 201}
30178c30 202
ce3f5e87 203//____________________________________________________________________
30178c30 204AliMUON::AliMUON(const AliMUON& rMUON)
205 : AliDetector(rMUON)
a9e2aefa 206{
30178c30 207// Protected copy constructor
208
8c343c7c 209 AliFatal("Not implemented.");
a9e2aefa 210}
30178c30 211
ce3f5e87 212//____________________________________________________________________
fe4da5cc 213AliMUON::~AliMUON()
214{
de05461e 215// Destructor
8c343c7c 216 AliDebug(1,"Calling AliMUON destructor");
ce3f5e87 217 fIshunt = 0;
218 if (fMerger) delete fMerger;
d1cd2474 219
c6df4ef2 220 if (fChambers){
221 fChambers->Delete();
222 delete fChambers;
223 }
224 if (fTriggerCircuits){
225 fTriggerCircuits->Delete();
226 delete fTriggerCircuits;
227 }
228 delete fMUONData;
d4bb94a1 229 delete fGeometryBuilder;
e118b27e 230 delete fDEIndexing;
8109731a 231 delete fFactory;
d1cd2474 232}
30178c30 233
234//________________________________________________________________________
235AliMUON& AliMUON::operator = (const AliMUON& rhs)
236{
237// Protected assignement operator
238
239 if (this == &rhs) return *this;
240
8c343c7c 241 AliFatal("Not implemented.");
30178c30 242
243 return *this;
244}
245
d1cd2474 246//_____________________________________________________________________________
247void AliMUON::AddGeometryBuilder(AliMUONVGeometryBuilder* geomBuilder)
248{
249// Adds the geometry builder to the list
250// ---
251
d4bb94a1 252 fGeometryBuilder->AddBuilder(geomBuilder);
a897a37a 253}
c9d10ab5 254
ce3f5e87 255//____________________________________________________________________
fe4da5cc 256void AliMUON::BuildGeometry()
257{
de05461e 258// Geometry for event display
c9d10ab5 259
a713db22 260 if (!fSegmentationType) {
261 AliFatal("No Segmentation Type defined.");
262 return;
263 }
264
002920d1 265
ae89d471 266// for (Int_t i = 0; i < AliMUONConstants::NCh(); i++)
267// this->Chamber(i).SegmentationModel2(1)->Draw("eventdisplay");// to be check !
268
002920d1 269
fe4da5cc 270}
c9d10ab5 271
ce3f5e87 272//__________________________________________________________________
273void AliMUON::SetTreeAddress()
a897a37a 274{
ce3f5e87 275 GetMUONData()->SetLoader(fLoader);
d1775029 276 // GetMUONData()->MakeBranch("D,S,RC");
277 // GetMUONData()->SetTreeAddress("H,D,S,RC");
52c9bc11 278 GetMUONData()->SetTreeAddress("H");
307d9d04 279 if (fHits != GetMUONData()->Hits()) {
ae31c826 280 if ( gAlice->GetMCApp() )
6d2f7cbd 281 if ( gAlice->GetMCApp()->GetHitLists() ) {
282 fHits = GetMUONData()->Hits();
283 gAlice->GetMCApp()->AddHitList(fHits); // For purifyKine, only necessary when Hit list is created in AliMUONData
284 }
307d9d04 285 }
6d2f7cbd 286 fHits = GetMUONData()->Hits(); // Added by Ivana to use the methods FisrtHit, NextHit of AliDetector
a9e2aefa 287}
288
fe4da5cc 289//___________________________________________
a9e2aefa 290void AliMUON::SetChambersZ(const Float_t *Z)
291{
292 // Set Z values for all chambers (tracking and trigger)
293 // from the array pointed to by "Z"
f665c1ea 294 for (Int_t ch = 0; ch < AliMUONConstants::NCh(); ch++)
2682e810 295 ((AliMUONChamber*) fChambers->At(ch))->SetZ(Z[ch]);
f665c1ea 296 return;
a9e2aefa 297}
ce3f5e87 298//_________________________________________________________________
a9e2aefa 299void AliMUON::SetChambersZToDefault()
fe4da5cc 300{
a9e2aefa 301 // Set Z values for all chambers (tracking and trigger)
302 // to default values
f665c1ea 303 SetChambersZ(AliMUONConstants::DefaultChamberZ());
a9e2aefa 304 return;
fe4da5cc 305}
ce3f5e87 306//_________________________________________________________________
a897a37a 307void AliMUON::SetChargeSlope(Int_t id, Float_t p1)
fe4da5cc 308{
de05461e 309// Set the inverse charge slope for chamber id
ce3f5e87 310 Int_t i=2*(id-1); //PH ((AliMUONChamber*) (*fChambers)[i])->SetSigmaIntegration(p1);
311 //PH ((AliMUONChamber*) (*fChambers)[i+1])->SetSigmaIntegration(p1);
2682e810 312 ((AliMUONChamber*) fChambers->At(i))->SetChargeSlope(p1);
313 ((AliMUONChamber*) fChambers->At(i+1))->SetChargeSlope(p1);
fe4da5cc 314}
ce3f5e87 315//__________________________________________________________________
a897a37a 316void AliMUON::SetChargeSpread(Int_t id, Float_t p1, Float_t p2)
fe4da5cc 317{
de05461e 318// Set sigma of charge spread for chamber id
fe4da5cc 319 Int_t i=2*(id-1);
2682e810 320 ((AliMUONChamber*) fChambers->At(i))->SetChargeSpread(p1,p2);
321 ((AliMUONChamber*) fChambers->At(i+1))->SetChargeSpread(p1,p2);
fe4da5cc 322}
ce3f5e87 323//___________________________________________________________________
a897a37a 324void AliMUON::SetSigmaIntegration(Int_t id, Float_t p1)
fe4da5cc 325{
de05461e 326// Set integration limits for charge spread
fe4da5cc 327 Int_t i=2*(id-1);
2682e810 328 ((AliMUONChamber*) fChambers->At(i))->SetSigmaIntegration(p1);
329 ((AliMUONChamber*) fChambers->At(i+1))->SetSigmaIntegration(p1);
fe4da5cc 330}
331
ce3f5e87 332//__________________________________________________________________
d09fafb0 333void AliMUON::SetMaxAdc(Int_t id, Int_t p1)
fe4da5cc 334{
de05461e 335// Set maximum number for ADCcounts (saturation)
fe4da5cc 336 Int_t i=2*(id-1);
2682e810 337 ((AliMUONChamber*) fChambers->At(i))->SetMaxAdc(p1);
338 ((AliMUONChamber*) fChambers->At(i+1))->SetMaxAdc(p1);
fe4da5cc 339}
d4bb94a1 340
ce3f5e87 341//__________________________________________________________________
a897a37a 342void AliMUON::SetMaxStepGas(Float_t p1)
fe4da5cc 343{
de05461e 344// Set stepsize in gas
ce3f5e87 345 fMaxStepGas=p1;
fe4da5cc 346}
ce3f5e87 347//__________________________________________________________________
a897a37a 348void AliMUON::SetMaxStepAlu(Float_t p1)
fe4da5cc 349{
de05461e 350// Set step size in Alu
fe4da5cc 351 fMaxStepAlu=p1;
352}
ce3f5e87 353//__________________________________________________________________
a897a37a 354void AliMUON::SetMaxDestepGas(Float_t p1)
fe4da5cc 355{
de05461e 356// Set maximum step size in Gas
fe4da5cc 357 fMaxDestepGas=p1;
358}
ce3f5e87 359//__________________________________________________________________
a897a37a 360void AliMUON::SetMaxDestepAlu(Float_t p1)
fe4da5cc 361{
de05461e 362// Set maximum step size in Alu
ce3f5e87 363 fMaxDestepAlu=p1;
fe4da5cc 364}
d4bb94a1 365
d1cd2474 366//____________________________________________________________________
367Float_t AliMUON::GetMaxStepGas() const
368{
369// Return stepsize in gas
370
371 return fMaxStepGas;
372}
373
374//____________________________________________________________________
375Float_t AliMUON::GetMaxStepAlu() const
376{
377// Return step size in Alu
378
379 return fMaxStepAlu;
380}
381
382//____________________________________________________________________
383Float_t AliMUON::GetMaxDestepGas() const
384{
385// Return maximum step size in Gas
386
387 return fMaxDestepGas;
388}
389
390//____________________________________________________________________
391Float_t AliMUON::GetMaxDestepAlu() const
392{
393// Return maximum step size in Gas
394
395 return fMaxDestepAlu;
396}
d4bb94a1 397
ce3f5e87 398//____________________________________________________________________
883746f3 399 void AliMUON::SetAlign(Bool_t align)
a713db22 400{
883746f3 401 // Sets option for alignement to geometry builder
402
403 fGeometryBuilder->SetAlign(align);
404}
405
406//____________________________________________________________________
a713db22 407void AliMUON::SetSegmentationModel(Int_t id, Int_t isec, AliMUONGeometrySegmentation* segmentation)
408{
409// Set the segmentation for chamber id cathode isec
410 ((AliMUONChamber*) fChambers->At(id))->SetSegmentationModel(isec, segmentation);
fe4da5cc 411
412}
ce3f5e87 413//____________________________________________________________________
a9e2aefa 414void AliMUON::SetResponseModel(Int_t id, AliMUONResponse *response)
fe4da5cc 415{
de05461e 416// Set the response for chamber id
2682e810 417 ((AliMUONChamber*) fChambers->At(id))->SetResponseModel(response);
fe4da5cc 418}
ce3f5e87 419//____________________________________________________________________
fe4da5cc 420void AliMUON::SetNsec(Int_t id, Int_t nsec)
421{
de05461e 422// Set number of segmented cathods for chamber id
2682e810 423 ((AliMUONChamber*) fChambers->At(id))->SetNsec(nsec);
fe4da5cc 424}
85a5290f 425//____________________________________________________________________
c92eb8ad 426AliDigitizer* AliMUON::CreateDigitizer(AliRunDigitizer* manager) const
85a5290f 427{
dbb4e169 428 return new AliMUONDigitizerv2(manager);
85a5290f 429}
ce3f5e87 430//_____________________________________________________________________
2ab0c725 431void AliMUON::SDigits2Digits()
432{
d963c261 433
434// write TreeD here
435
436 if (!fMerger) {
8c343c7c 437 AliDebug(1,"Create default AliMUONMerger ");
438 AliDebug(1," no merging, just digitization of 1 event will be done");
439 fMerger = new AliMUONMerger();
2ab0c725 440 }
d963c261 441 fMerger->Init();
442 fMerger->Digitise();
443 char hname[30];
88cb7938 444 // sprintf(hname,"TreeD%d",fLoader->GetHeader()->GetEvent());
445 fLoader->TreeD()->Write(hname,TObject::kOverwrite);
446 fLoader->TreeD()->Reset();
2ab0c725 447}
a9e2aefa 448
d1775029 449//_____________________________________________________________________
450void AliMUON::Hits2SDigits()
451{
dbb4e169 452 // Adaption of AliMUONSDigitizerv1 to be excuted by the AliSimulation framework
453 AliRunLoader* runLoader = fLoader->GetRunLoader();
454 AliRunDigitizer * manager = new AliRunDigitizer(1,1);
e191bb57 455 manager->SetInputStream(0,runLoader->GetFileName(),AliConfig::GetDefaultEventFolderName());
dbb4e169 456 AliMUONDigitizer * dMUON = new AliMUONSDigitizerv1(manager);
457 fLoader->LoadHits("READ");
458 for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
459 runLoader->GetEvent(iEvent);
460 dMUON->Exec("");
461 }
462 fLoader->UnloadHits();
d1775029 463}
39e39255 464//_____________________________________________________________________
465void AliMUON::Digits2Raw()
466{
467 // convert digits of the current event to raw data
468 AliMUONRawData* rawData;
469
470 rawData = new AliMUONRawData(fLoader);
8c343c7c 471 if (!rawData->WriteRawData()) AliInfo("pb writting raw data");
39e39255 472 delete rawData;
473 return;
474}
ce3f5e87 475//_______________________________________________________________________
3fa6cfdd 476AliLoader* AliMUON::MakeLoader(const char* topfoldername)
477{
478//builds standard getter (AliLoader type)
479//if detector wants to use castomized getter, it must overload this method
480
8c343c7c 481
482 AliDebug(1,Form("Creating standard getter for detector %s. Top folder is %s.",
483 GetName(),topfoldername));
128843d9 484 fLoader = new AliLoader(GetName(),topfoldername);
ce3f5e87 485 fMUONData = new AliMUONData(fLoader,GetName(),GetName());
d652f85c 486 fMUONData->SetSplitLevel(fSplitLevel);
3fa6cfdd 487 return fLoader;
488}
ce3f5e87 489//_______________________________________________________________________
a9e2aefa 490AliMUONPadHit* AliMUON::FirstPad(AliMUONHit* hit, TClonesArray *clusters)
fe4da5cc 491{
5a35ac13 492// to be removed
fe4da5cc 493 // Initialise the pad iterator
494 // Return the address of the first padhit for hit
a897a37a 495 TClonesArray *theClusters = clusters;
fe4da5cc 496 Int_t nclust = theClusters->GetEntriesFast();
1bd28025 497 if (nclust && hit->PHlast() > 0) {
30178c30 498 fMaxIterPad=hit->PHlast();
499 fCurIterPad=hit->PHfirst();
500 return (AliMUONPadHit*) clusters->UncheckedAt(fCurIterPad-1);
fe4da5cc 501 } else {
502 return 0;
503 }
504}
ce3f5e87 505//_______________________________________________________________________
a9e2aefa 506AliMUONPadHit* AliMUON::NextPad(TClonesArray *clusters)
fe4da5cc 507{
5a35ac13 508 // To be removed
de05461e 509// Get next pad (in iterator)
510//
30178c30 511 fCurIterPad++;
512 if (fCurIterPad <= fMaxIterPad) {
513 return (AliMUONPadHit*) clusters->UncheckedAt(fCurIterPad-1);
fe4da5cc 514 } else {
515 return 0;
516 }
517}
ce3f5e87 518//_______________________________________________________________________
a897a37a 519
520AliMUONRawCluster *AliMUON::RawCluster(Int_t ichamber, Int_t icathod, Int_t icluster)
521{
de05461e 522//
523// Return rawcluster (icluster) for chamber ichamber and cathode icathod
524// Obsolete ??
ce3f5e87 525 TClonesArray *muonRawCluster = GetMUONData()->RawClusters(ichamber);
a897a37a 526 ResetRawClusters();
88cb7938 527 TTree *treeR = fLoader->TreeR();
a9e2aefa 528 Int_t nent=(Int_t)treeR->GetEntries();
529 treeR->GetEvent(nent-2+icathod-1);
530 //treeR->GetEvent(icathod);
531 //Int_t nrawcl = (Int_t)muonRawCluster->GetEntriesFast();
a897a37a 532
a9e2aefa 533 AliMUONRawCluster * mRaw = (AliMUONRawCluster*)muonRawCluster->UncheckedAt(icluster);
a897a37a 534 //printf("RawCluster _ nent nrawcl icluster mRaw %d %d %d%p\n",nent,nrawcl,icluster,mRaw);
535
536 return mRaw;
537}
ce3f5e87 538//________________________________________________________________________
1cae9436 539void AliMUON::SetMerger(AliMUONMerger* merger)
540{
541// Set pointer to merger
542 fMerger = merger;
543}
ce3f5e87 544//________________________________________________________________________
1cae9436 545AliMUONMerger* AliMUON::Merger()
546{
547// Return pointer to merger
548 return fMerger;
549}
ce01996f 550/* PH Commented out waiting for correct implementation
ce3f5e87 551//________________________________________________________________________
48b32e42 552void AliMUON::RemapTrackHitIDs(Int_t* map)
553{
554// Remaps the track numbers in the hits arrays, so that they correspond
555// to the entry indices in the Kine tree.
556// The correspondance is not direct. To get the real index into the Kine tree
557// compute the particle index as follows:
558//
559// num_primaries = AliStack::GetNprimary();
560// num_tracks = AliStack::GetNtracks();
561// track = AliMUONHit::Track()
562//
563// if (track < num_primaries)
564// particleindex = track + num_tracks - num_primaries;
565// else
566// particleindex = track - num_primaries;
567
568 // Remap the track numbers based on the specified map.
569 AliMUONData* data = GetMUONData();
570 TClonesArray* hits = data->Hits();
571 for (Int_t i = 0; i < hits->GetEntriesFast(); i++)
572 {
573 AliMUONHit* hit = static_cast<AliMUONHit*>( hits->At(i) );
574 hit->SetTrack( map[hit->Track()] );
575 };
576};
ce01996f 577*/