Updates to RAW package. Used during the Computing DC.
[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"
8109731a 70#include "AliMUONFactory.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),
b2312f80 99 fSegmentationType(1),// set to 1 default wise old 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),
b2312f80 128 fSegmentationType(1),// set to 1 default wise old 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
265 if (fSegmentationType == 1) {
266 for (Int_t i=0; i<7; i++) {
267 for (Int_t j=0; j<2; j++) {
268 Int_t id=2*i+j+1;
269 this->Chamber(id-1).SegmentationModel(1)->Draw("eventdisplay");
270 }
271 }
272 } else {
ae89d471 273// for (Int_t i = 0; i < AliMUONConstants::NCh(); i++)
274// this->Chamber(i).SegmentationModel2(1)->Draw("eventdisplay");// to be check !
275
276 }
fe4da5cc 277}
c9d10ab5 278
ce3f5e87 279//__________________________________________________________________
280void AliMUON::SetTreeAddress()
a897a37a 281{
ce3f5e87 282 GetMUONData()->SetLoader(fLoader);
d1775029 283 // GetMUONData()->MakeBranch("D,S,RC");
284 // GetMUONData()->SetTreeAddress("H,D,S,RC");
52c9bc11 285 GetMUONData()->SetTreeAddress("H");
307d9d04 286 if (fHits != GetMUONData()->Hits()) {
ae31c826 287 if ( gAlice->GetMCApp() )
6d2f7cbd 288 if ( gAlice->GetMCApp()->GetHitLists() ) {
289 fHits = GetMUONData()->Hits();
290 gAlice->GetMCApp()->AddHitList(fHits); // For purifyKine, only necessary when Hit list is created in AliMUONData
291 }
307d9d04 292 }
6d2f7cbd 293 fHits = GetMUONData()->Hits(); // Added by Ivana to use the methods FisrtHit, NextHit of AliDetector
a9e2aefa 294}
295
ce3f5e87 296//____________________________________________________________________
a9e2aefa 297void AliMUON::SetPadSize(Int_t id, Int_t isec, Float_t p1, Float_t p2)
298{
de05461e 299// Set the pad size for chamber id and cathode isec
a9e2aefa 300 Int_t i=2*(id-1);
2682e810 301 ((AliMUONChamber*) fChambers->At(i)) ->SetPadSize(isec,p1,p2);
302 ((AliMUONChamber*) fChambers->At(i+1))->SetPadSize(isec,p1,p2);
a897a37a 303}
304
fe4da5cc 305//___________________________________________
a9e2aefa 306void AliMUON::SetChambersZ(const Float_t *Z)
307{
308 // Set Z values for all chambers (tracking and trigger)
309 // from the array pointed to by "Z"
f665c1ea 310 for (Int_t ch = 0; ch < AliMUONConstants::NCh(); ch++)
2682e810 311 ((AliMUONChamber*) fChambers->At(ch))->SetZ(Z[ch]);
f665c1ea 312 return;
a9e2aefa 313}
ce3f5e87 314//_________________________________________________________________
a9e2aefa 315void AliMUON::SetChambersZToDefault()
fe4da5cc 316{
a9e2aefa 317 // Set Z values for all chambers (tracking and trigger)
318 // to default values
f665c1ea 319 SetChambersZ(AliMUONConstants::DefaultChamberZ());
a9e2aefa 320 return;
fe4da5cc 321}
ce3f5e87 322//_________________________________________________________________
a897a37a 323void AliMUON::SetChargeSlope(Int_t id, Float_t p1)
fe4da5cc 324{
de05461e 325// Set the inverse charge slope for chamber id
ce3f5e87 326 Int_t i=2*(id-1); //PH ((AliMUONChamber*) (*fChambers)[i])->SetSigmaIntegration(p1);
327 //PH ((AliMUONChamber*) (*fChambers)[i+1])->SetSigmaIntegration(p1);
2682e810 328 ((AliMUONChamber*) fChambers->At(i))->SetChargeSlope(p1);
329 ((AliMUONChamber*) fChambers->At(i+1))->SetChargeSlope(p1);
fe4da5cc 330}
ce3f5e87 331//__________________________________________________________________
a897a37a 332void AliMUON::SetChargeSpread(Int_t id, Float_t p1, Float_t p2)
fe4da5cc 333{
de05461e 334// Set sigma of charge spread for chamber id
fe4da5cc 335 Int_t i=2*(id-1);
2682e810 336 ((AliMUONChamber*) fChambers->At(i))->SetChargeSpread(p1,p2);
337 ((AliMUONChamber*) fChambers->At(i+1))->SetChargeSpread(p1,p2);
fe4da5cc 338}
ce3f5e87 339//___________________________________________________________________
a897a37a 340void AliMUON::SetSigmaIntegration(Int_t id, Float_t p1)
fe4da5cc 341{
de05461e 342// Set integration limits for charge spread
fe4da5cc 343 Int_t i=2*(id-1);
2682e810 344 ((AliMUONChamber*) fChambers->At(i))->SetSigmaIntegration(p1);
345 ((AliMUONChamber*) fChambers->At(i+1))->SetSigmaIntegration(p1);
fe4da5cc 346}
347
ce3f5e87 348//__________________________________________________________________
d09fafb0 349void AliMUON::SetMaxAdc(Int_t id, Int_t p1)
fe4da5cc 350{
de05461e 351// Set maximum number for ADCcounts (saturation)
fe4da5cc 352 Int_t i=2*(id-1);
2682e810 353 ((AliMUONChamber*) fChambers->At(i))->SetMaxAdc(p1);
354 ((AliMUONChamber*) fChambers->At(i+1))->SetMaxAdc(p1);
fe4da5cc 355}
d4bb94a1 356
ce3f5e87 357//__________________________________________________________________
a897a37a 358void AliMUON::SetMaxStepGas(Float_t p1)
fe4da5cc 359{
de05461e 360// Set stepsize in gas
ce3f5e87 361 fMaxStepGas=p1;
fe4da5cc 362}
ce3f5e87 363//__________________________________________________________________
a897a37a 364void AliMUON::SetMaxStepAlu(Float_t p1)
fe4da5cc 365{
de05461e 366// Set step size in Alu
fe4da5cc 367 fMaxStepAlu=p1;
368}
ce3f5e87 369//__________________________________________________________________
a897a37a 370void AliMUON::SetMaxDestepGas(Float_t p1)
fe4da5cc 371{
de05461e 372// Set maximum step size in Gas
fe4da5cc 373 fMaxDestepGas=p1;
374}
ce3f5e87 375//__________________________________________________________________
a897a37a 376void AliMUON::SetMaxDestepAlu(Float_t p1)
fe4da5cc 377{
de05461e 378// Set maximum step size in Alu
ce3f5e87 379 fMaxDestepAlu=p1;
fe4da5cc 380}
d4bb94a1 381
d1cd2474 382//____________________________________________________________________
383Float_t AliMUON::GetMaxStepGas() const
384{
385// Return stepsize in gas
386
387 return fMaxStepGas;
388}
389
390//____________________________________________________________________
391Float_t AliMUON::GetMaxStepAlu() const
392{
393// Return step size in Alu
394
395 return fMaxStepAlu;
396}
397
398//____________________________________________________________________
399Float_t AliMUON::GetMaxDestepGas() const
400{
401// Return maximum step size in Gas
402
403 return fMaxDestepGas;
404}
405
406//____________________________________________________________________
407Float_t AliMUON::GetMaxDestepAlu() const
408{
409// Return maximum step size in Gas
410
411 return fMaxDestepAlu;
412}
d4bb94a1 413
ce3f5e87 414//____________________________________________________________________
883746f3 415 void AliMUON::SetAlign(Bool_t align)
a713db22 416{
883746f3 417 // Sets option for alignement to geometry builder
418
419 fGeometryBuilder->SetAlign(align);
420}
421
422//____________________________________________________________________
a30a000f 423void AliMUON::SetSegmentationModel(Int_t id, Int_t isec, AliSegmentation *segmentation)
fe4da5cc 424{
de05461e 425// Set the segmentation for chamber id cathode isec
2682e810 426 ((AliMUONChamber*) fChambers->At(id))->SetSegmentationModel(isec, segmentation);
a713db22 427}
428//____________________________________________________________________
429void AliMUON::SetSegmentationModel(Int_t id, Int_t isec, AliMUONGeometrySegmentation* segmentation)
430{
431// Set the segmentation for chamber id cathode isec
432 ((AliMUONChamber*) fChambers->At(id))->SetSegmentationModel(isec, segmentation);
fe4da5cc 433
434}
ce3f5e87 435//____________________________________________________________________
a9e2aefa 436void AliMUON::SetResponseModel(Int_t id, AliMUONResponse *response)
fe4da5cc 437{
de05461e 438// Set the response for chamber id
2682e810 439 ((AliMUONChamber*) fChambers->At(id))->SetResponseModel(response);
fe4da5cc 440}
ce3f5e87 441//____________________________________________________________________
fe4da5cc 442void AliMUON::SetNsec(Int_t id, Int_t nsec)
443{
de05461e 444// Set number of segmented cathods for chamber id
2682e810 445 ((AliMUONChamber*) fChambers->At(id))->SetNsec(nsec);
fe4da5cc 446}
85a5290f 447//____________________________________________________________________
c92eb8ad 448AliDigitizer* AliMUON::CreateDigitizer(AliRunDigitizer* manager) const
85a5290f 449{
dbb4e169 450 return new AliMUONDigitizerv2(manager);
85a5290f 451}
ce3f5e87 452//_____________________________________________________________________
2ab0c725 453void AliMUON::SDigits2Digits()
454{
d963c261 455
456// write TreeD here
457
458 if (!fMerger) {
8c343c7c 459 AliDebug(1,"Create default AliMUONMerger ");
460 AliDebug(1," no merging, just digitization of 1 event will be done");
461 fMerger = new AliMUONMerger();
2ab0c725 462 }
d963c261 463 fMerger->Init();
464 fMerger->Digitise();
465 char hname[30];
88cb7938 466 // sprintf(hname,"TreeD%d",fLoader->GetHeader()->GetEvent());
467 fLoader->TreeD()->Write(hname,TObject::kOverwrite);
468 fLoader->TreeD()->Reset();
2ab0c725 469}
a9e2aefa 470
d1775029 471//_____________________________________________________________________
472void AliMUON::Hits2SDigits()
473{
dbb4e169 474 // Adaption of AliMUONSDigitizerv1 to be excuted by the AliSimulation framework
475 AliRunLoader* runLoader = fLoader->GetRunLoader();
476 AliRunDigitizer * manager = new AliRunDigitizer(1,1);
e191bb57 477 manager->SetInputStream(0,runLoader->GetFileName(),AliConfig::GetDefaultEventFolderName());
dbb4e169 478 AliMUONDigitizer * dMUON = new AliMUONSDigitizerv1(manager);
479 fLoader->LoadHits("READ");
480 for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
481 runLoader->GetEvent(iEvent);
482 dMUON->Exec("");
483 }
484 fLoader->UnloadHits();
d1775029 485}
39e39255 486//_____________________________________________________________________
487void AliMUON::Digits2Raw()
488{
489 // convert digits of the current event to raw data
490 AliMUONRawData* rawData;
491
492 rawData = new AliMUONRawData(fLoader);
8c343c7c 493 if (!rawData->WriteRawData()) AliInfo("pb writting raw data");
39e39255 494 delete rawData;
495 return;
496}
ce3f5e87 497//_______________________________________________________________________
3fa6cfdd 498AliLoader* AliMUON::MakeLoader(const char* topfoldername)
499{
500//builds standard getter (AliLoader type)
501//if detector wants to use castomized getter, it must overload this method
502
8c343c7c 503
504 AliDebug(1,Form("Creating standard getter for detector %s. Top folder is %s.",
505 GetName(),topfoldername));
128843d9 506 fLoader = new AliLoader(GetName(),topfoldername);
ce3f5e87 507 fMUONData = new AliMUONData(fLoader,GetName(),GetName());
d652f85c 508 fMUONData->SetSplitLevel(fSplitLevel);
3fa6cfdd 509 return fLoader;
510}
ce3f5e87 511//_______________________________________________________________________
a9e2aefa 512AliMUONPadHit* AliMUON::FirstPad(AliMUONHit* hit, TClonesArray *clusters)
fe4da5cc 513{
5a35ac13 514// to be removed
fe4da5cc 515 // Initialise the pad iterator
516 // Return the address of the first padhit for hit
a897a37a 517 TClonesArray *theClusters = clusters;
fe4da5cc 518 Int_t nclust = theClusters->GetEntriesFast();
1bd28025 519 if (nclust && hit->PHlast() > 0) {
30178c30 520 fMaxIterPad=hit->PHlast();
521 fCurIterPad=hit->PHfirst();
522 return (AliMUONPadHit*) clusters->UncheckedAt(fCurIterPad-1);
fe4da5cc 523 } else {
524 return 0;
525 }
526}
ce3f5e87 527//_______________________________________________________________________
a9e2aefa 528AliMUONPadHit* AliMUON::NextPad(TClonesArray *clusters)
fe4da5cc 529{
5a35ac13 530 // To be removed
de05461e 531// Get next pad (in iterator)
532//
30178c30 533 fCurIterPad++;
534 if (fCurIterPad <= fMaxIterPad) {
535 return (AliMUONPadHit*) clusters->UncheckedAt(fCurIterPad-1);
fe4da5cc 536 } else {
537 return 0;
538 }
539}
ce3f5e87 540//_______________________________________________________________________
a897a37a 541
542AliMUONRawCluster *AliMUON::RawCluster(Int_t ichamber, Int_t icathod, Int_t icluster)
543{
de05461e 544//
545// Return rawcluster (icluster) for chamber ichamber and cathode icathod
546// Obsolete ??
ce3f5e87 547 TClonesArray *muonRawCluster = GetMUONData()->RawClusters(ichamber);
a897a37a 548 ResetRawClusters();
88cb7938 549 TTree *treeR = fLoader->TreeR();
a9e2aefa 550 Int_t nent=(Int_t)treeR->GetEntries();
551 treeR->GetEvent(nent-2+icathod-1);
552 //treeR->GetEvent(icathod);
553 //Int_t nrawcl = (Int_t)muonRawCluster->GetEntriesFast();
a897a37a 554
a9e2aefa 555 AliMUONRawCluster * mRaw = (AliMUONRawCluster*)muonRawCluster->UncheckedAt(icluster);
a897a37a 556 //printf("RawCluster _ nent nrawcl icluster mRaw %d %d %d%p\n",nent,nrawcl,icluster,mRaw);
557
558 return mRaw;
559}
ce3f5e87 560//________________________________________________________________________
1cae9436 561void AliMUON::SetMerger(AliMUONMerger* merger)
562{
563// Set pointer to merger
564 fMerger = merger;
565}
ce3f5e87 566//________________________________________________________________________
1cae9436 567AliMUONMerger* AliMUON::Merger()
568{
569// Return pointer to merger
570 return fMerger;
571}
ce01996f 572/* PH Commented out waiting for correct implementation
ce3f5e87 573//________________________________________________________________________
48b32e42 574void AliMUON::RemapTrackHitIDs(Int_t* map)
575{
576// Remaps the track numbers in the hits arrays, so that they correspond
577// to the entry indices in the Kine tree.
578// The correspondance is not direct. To get the real index into the Kine tree
579// compute the particle index as follows:
580//
581// num_primaries = AliStack::GetNprimary();
582// num_tracks = AliStack::GetNtracks();
583// track = AliMUONHit::Track()
584//
585// if (track < num_primaries)
586// particleindex = track + num_tracks - num_primaries;
587// else
588// particleindex = track - num_primaries;
589
590 // Remap the track numbers based on the specified map.
591 AliMUONData* data = GetMUONData();
592 TClonesArray* hits = data->Hits();
593 for (Int_t i = 0; i < hits->GetEntriesFast(); i++)
594 {
595 AliMUONHit* hit = static_cast<AliMUONHit*>( hits->At(i) );
596 hit->SetTrack( map[hit->Track()] );
597 };
598};
ce01996f 599*/