MUON/AliMUONv1.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *      SigmaEffect_thetadegrees                                                                  *
   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 purpeateose. It is      *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /* $Id$ */
  17
  18 /////////////////////////////////////////////////////////
  19 //  Manager and hits classes for set:MUON version 1    //
  20 /////////////////////////////////////////////////////////
  21
  22 #include <TRandom.h>
  23 #include <TF1.h>
  24 #include <TClonesArray.h>
  25 #include <TRandom.h>
  26 #include <TVirtualMC.h>
  27 #include <TGeoMatrix.h>
  28
  29 #include "AliMUONv1.h"
  30 #include "AliConst.h"
  31 #include "AliMUONConstants.h"
  32 #include "AliMUONSegFactory.h"
  33 #include "AliMUONResponseFactory.h"
  34 #include "AliMUONSegmentation.h"
  35 #include "AliMUONHit.h"
  36 #include "AliMUONTriggerCircuit.h"
  37 #include "AliMUONGeometryBuilder.h"
  38 #include "AliMUONGeometry.h"
  39 #include "AliMUONGeometryTransformer.h"
  40 #include "AliMUONGeometryModule.h"
  41 #include "AliMUONGeometrySVMap.h"
  42 #include "AliMUONGeometryDetElement.h"
  43 #include "AliMagF.h"
  44 #include "AliRun.h"
  45 #include "AliMC.h"
  46 #include "AliLog.h"
  47
  48 #include <string>
  49
  50 ClassImp(AliMUONv1)
  51
  52 //___________________________________________
  53 AliMUONv1::AliMUONv1()
  54   : AliMUON(),
  55     fStepManagerVersionOld(kTRUE),
  56     fStepManagerVersionDE(kTRUE),
  57     fAngleEffect(kTRUE),
  58     fStepMaxInActiveGas(0.6),
  59     fStepSum(0x0),
  60     fDestepSum(0x0),
  61     fTrackMomentum(),
  62     fTrackPosition(),
  63     fElossRatio(0x0),
  64     fAngleEffect10(0x0),
  65     fAngleEffectNorma(0x0)
  66 {
  67 // Default constructor
  68         AliDebug(1,Form("default (empty) ctor this=%p",this));
  69 }
  70
  71 //___________________________________________
  72 AliMUONv1::AliMUONv1(const char *name, const char *title)
  73   : AliMUON(name,title),
  74     fStepManagerVersionOld(kTRUE),
  75     fStepManagerVersionDE(kTRUE),
  76     fAngleEffect(kTRUE),
  77     fStepMaxInActiveGas(0.6),
  78     fStepSum(0x0),
  79     fDestepSum(0x0),
  80     fTrackMomentum(),
  81     fTrackPosition(),
  82     fElossRatio(0x0),
  83     fAngleEffect10(0x0),
  84     fAngleEffectNorma(0x0)
  85 {
  86 // Standard onstructor
  87
  88         AliDebug(1,Form("ctor this=%p",this));
  89
  90     // By default include all stations
  91
  92     fStepSum   = new Float_t [AliMUONConstants::NCh()];
  93     fDestepSum = new Float_t [AliMUONConstants::NCh()];
  94     for (Int_t i=0; i<AliMUONConstants::NCh(); i++) {
  95       fStepSum[i] =0.0;
  96       fDestepSum[i]=0.0;
  97     }
  98     // Ratio of particle mean eloss with respect MIP's Khalil Boudjemline, sep 2003, PhD.Thesis and Particle Data Book
  99     fElossRatio = new TF1("ElossRatio","[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x",0.5,5.);
 100     fElossRatio->SetParameter(0,1.02138);
 101     fElossRatio->SetParameter(1,-9.54149e-02);
 102     fElossRatio->SetParameter(2,+7.83433e-02);
 103     fElossRatio->SetParameter(3,-9.98208e-03);
 104     fElossRatio->SetParameter(4,+3.83279e-04);
 105
 106     // Angle effect in tracking chambers at theta =10 degres as a function of ElossRatio (Khalil BOUDJEMLINE sep 2003 Ph.D Thesis) (in micrometers)
 107     fAngleEffect10 = new TF1("AngleEffect10","[0]+[1]*x+[2]*x*x",0.5,3.0);
 108     fAngleEffect10->SetParameter(0, 1.90691e+02);
 109     fAngleEffect10->SetParameter(1,-6.62258e+01);
 110     fAngleEffect10->SetParameter(2,+1.28247e+01);
 111     // Angle effect: Normalisation form theta=10 degres to theta between 0 and 10 (Khalil BOUDJEMLINE sep 2003 Ph.D Thesis)
 112     // Angle with respect to the wires assuming that chambers are perpendicular to the z axis.
 113     fAngleEffectNorma = new TF1("AngleEffectNorma","[0]+[1]*x+[2]*x*x+[3]*x*x*x",0.0,10.0);
 114     fAngleEffectNorma->SetParameter(0,4.148);
 115     fAngleEffectNorma->SetParameter(1,-6.809e-01);
 116     fAngleEffectNorma->SetParameter(2,5.151e-02);
 117     fAngleEffectNorma->SetParameter(3,-1.490e-03);
 118 }
 119
 120 //_____________________________________________________________________________
 121 AliMUONv1::AliMUONv1(const AliMUONv1& right)
 122   : AliMUON(right)
 123 {
 124   // copy constructor (not implemented)
 125
 126   AliFatal("Copy constructor not provided.");
 127 }
 128
 129 //___________________________________________
 130 AliMUONv1::~AliMUONv1()
 131 {
 132 // Destructor
 133         AliDebug(1,Form("dtor this=%p",this));
 134   delete [] fStepSum;
 135   delete [] fDestepSum;
 136   delete fElossRatio;
 137   delete fAngleEffect10;
 138   delete fAngleEffectNorma;
 139 }
 140
 141 //_____________________________________________________________________________
 142 AliMUONv1& AliMUONv1::operator=(const AliMUONv1& right)
 143 {
 144   // assignement operator (not implemented)
 145
 146   // check assignement to self
 147   if (this == &right) return *this;
 148
 149   AliFatal("Assignement operator not provided.");
 150
 151   return *this;
 152 }
 153
 154 //__________________________________________________
 155 void AliMUONv1::CreateGeometry()
 156 {
 157 //
 158 // Construct geometry using geometry builder
 159 //
 160
 161   fGeometryBuilder->CreateGeometry();
 162 }
 163
 164 //________________________________________________________________
 165 void AliMUONv1::CreateMaterials()
 166 {
 167 //
 168 // Construct materials using geometry builder
 169 //
 170
 171   fGeometryBuilder->CreateMaterials();
 172 }
 173
 174 //___________________________________________
 175 void AliMUONv1::Init()
 176 {
 177   AliDebug(1,"Start Init for version 1 - CPC chamber type");
 178
 179   //
 180   // Initialize geometry
 181   //
 182   fGeometryBuilder->InitGeometry();
 183   AliDebug(1,"Finished Init for version 1 - CPC chamber type");
 184
 185
 186   std::string ftype(GetTitle());
 187
 188   // Build segmentation
 189   // using geometry parametrisation
 190   //
 191   AliMUONSegFactory segFactory(GetGeometryTransformer());
 192   fSegmentation = segFactory.CreateSegmentation(ftype);
 193
 194   if (!fSegmentation) {
 195     AliFatal(Form("Wrong factory type : %s",ftype.c_str()));
 196   }
 197
 198   // Build response
 199   //
 200   AliMUONResponseFactory respFactory("default");
 201   respFactory.Build(this);
 202
 203
 204   // Initialize segmentation
 205   //
 206   fSegmentation->Init();
 207
 208   // Initialize trigger circuits
 209   //
 210   for (Int_t i=0; i<AliMUONConstants::NTriggerCircuit(); i++)  {
 211     AliMUONTriggerCircuit* c = (AliMUONTriggerCircuit*)(fTriggerCircuits->At(i));
 212     c->Init(i);
 213 //    c->Print();
 214   }
 215
 216 }
 217
 218 //__________________________________________________________________
 219 Int_t  AliMUONv1::GetChamberId(Int_t volId) const
 220 {
 221 // Check if the volume with specified  volId is a sensitive volume (gas)
 222 // of some chamber and returns the chamber number;
 223 // if not sensitive volume - return 0.
 224 // ---
 225
 226   for (Int_t i = 0; i < AliMUONConstants::NCh(); i++) {
 227     if ( GetGeometry()->GetModule(i)->IsSensitiveVolume(volId) )
 228       return i+1;
 229   }
 230
 231   return 0;
 232 }
 233
 234 //_______________________________________________________________________________
 235 TString  AliMUONv1::CurrentVolumePath() const
 236 {
 237 // Returns current volume path
 238 // (Could be removed when this function is available via gMC)
 239 // ---
 240
 241   TString path = "";
 242   TString name;
 243   Int_t copyNo;
 244   Int_t imother = 0;
 245   do {
 246     name = gMC->CurrentVolOffName(imother);
 247     gMC->CurrentVolOffID(imother++, copyNo);
 248     TString add = "/";
 249     add += name;
 250     add += ".";
 251     add += copyNo;
 252     path.Insert(0,add);
 253   }
 254   while ( name != TString("ALIC") );
 255
 256   return path;
 257 }
 258
 259 //_______________________________________________________________________________
 260 void AliMUONv1::StepManager()
 261 {
 262   // Stepmanager for the chambers
 263   // TBR
 264
 265  if (fStepManagerVersionDE) {
 266     StepManager2();
 267     return;
 268   }
 269
 270
 271   // Only charged tracks
 272   if( !(gMC->TrackCharge()) ) return;
 273   // Only charged tracks
 274
 275   // Only gas gap inside chamber
 276   // Tag chambers and record hits when track enters
 277   static Int_t   idvol=-1;
 278   Int_t   iChamber=0;
 279   Int_t   id=0;
 280   Int_t   copy;
 281   const  Float_t kBig = 1.e10;
 282
 283
 284   //
 285   // Only gas gap inside chamber
 286   // Tag chambers and record hits when track enters
 287   id=gMC->CurrentVolID(copy);
 288   iChamber = GetChamberId(id);
 289   idvol = iChamber -1;
 290
 291   if (idvol == -1) return;
 292
 293   // Filling TrackRefs file for MUON. Our Track references are the active volume of the chambers
 294   if ( (gMC->IsTrackEntering() || gMC->IsTrackExiting() ) )
 295     AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber());
 296
 297    if( gMC->IsTrackEntering() ) {
 298      Float_t theta = fTrackMomentum.Theta();
 299      if ((TMath::Pi()-theta)*kRaddeg>=15.) gMC->SetMaxStep(fStepMaxInActiveGas); // We use Pi-theta because z is negative
 300   }
 301
 302    //  AliDebug(1,
 303    //    Form("Active volume found %d chamber %d Z chamber is %f ",idvol,iChamber,
 304    //    ( (AliMUONChamber*)(*fChambers)[idvol])->Z()));
 305    // Particule id and mass,
 306   Int_t     ipart = gMC->TrackPid();
 307   Float_t   mass  = gMC->TrackMass();
 308
 309   fDestepSum[idvol]+=gMC->Edep();
 310   // Get current particle id (ipart), track position (pos)  and momentum (mom)
 311   if ( fStepSum[idvol]==0.0 )  gMC->TrackMomentum(fTrackMomentum);
 312   fStepSum[idvol]+=gMC->TrackStep();
 313
 314   //  if(AliDebugLevel()) {
 315   //  AliDebug(1,
 316   //         Form("iChamber %d, Particle %d, theta %f phi %f mass %f StepSum %f eloss %g",
 317   //              iChamber,ipart, fTrackMomentum.Theta()*kRaddeg, fTrackMomentum.Phi()*kRaddeg,
 318   //              mass, fStepSum[idvol], gMC->Edep()));
 319   //  AliDebug(1,
 320   //         Form("Track Momentum %f %f %f", fTrackMomentum.X(), fTrackMomentum.Y(),
 321   //              fTrackMomentum.Z()));
 322   //  gMC->TrackPosition(fTrackPosition);
 323   //  AliDebug(1,
 324   //         Form("Track Position %f %f %f",fTrackPosition.X(),fTrackPosition.Y(),
 325   //      fTrackPosition.Z())) ;
 326   //   }
 327
 328   // Track left chamber or StepSum larger than fStepMaxInActiveGas
 329   if ( gMC->IsTrackExiting() ||
 330        gMC->IsTrackStop() ||
 331        gMC->IsTrackDisappeared()||
 332        (fStepSum[idvol]>fStepMaxInActiveGas) ) {
 333
 334     if   ( gMC->IsTrackExiting() ||
 335            gMC->IsTrackStop() ||
 336            gMC->IsTrackDisappeared() ) gMC->SetMaxStep(kBig);
 337
 338     gMC->TrackPosition(fTrackPosition);
 339     Float_t theta = fTrackMomentum.Theta();
 340     Float_t phi   = fTrackMomentum.Phi();
 341
 342     TLorentzVector backToWire( fStepSum[idvol]/2.*sin(theta)*cos(phi),
 343                                fStepSum[idvol]/2.*sin(theta)*sin(phi),
 344                                fStepSum[idvol]/2.*cos(theta),0.0       );
 345     //     AliDebug(1,Form("Exit: Track Position %f %f %f",fTrackPosition.X(),
 346     //                     fTrackPosition.Y(),fTrackPosition.Z())) ;
 347     //     AliDebug(1,Form("Exit: Track backToWire %f %f %f",backToWire.X(),
 348     //                     backToWire.Y(),backToWire.Z()) ;
 349     fTrackPosition-=backToWire;
 350
 351     //-------------- Angle effect
 352     // Ratio between energy loss of particle and Mip as a function of BetaGamma of particle (Energy/Mass)
 353
 354     Float_t betaxGamma    = fTrackMomentum.P()/mass;//  pc/mc2
 355     Float_t sigmaEffect10degrees;
 356     Float_t sigmaEffectThetadegrees;
 357     Float_t eLossParticleELossMip;
 358     Float_t yAngleEffect=0.;
 359     Float_t thetawires      =  TMath::Abs( TMath::ASin( TMath::Sin(TMath::Pi()-theta) * TMath::Sin(phi) ) );// We use Pi-theta because z is negative
 360
 361
 362     if (fAngleEffect){
 363     if ( (betaxGamma >3.2)   &&  (thetawires*kRaddeg<=15.) ) {
 364       betaxGamma=TMath::Log(betaxGamma);
 365       eLossParticleELossMip = fElossRatio->Eval(betaxGamma);
 366       // 10 degrees is a reference for a model (arbitrary)
 367       sigmaEffect10degrees=fAngleEffect10->Eval(eLossParticleELossMip);// in micrometers
 368       // Angle with respect to the wires assuming that chambers are perpendicular to the z axis.
 369       sigmaEffectThetadegrees =  sigmaEffect10degrees/fAngleEffectNorma->Eval(thetawires*kRaddeg);  // For 5mm gap
 370       if ( (iChamber==1)  ||  (iChamber==2) )
 371         sigmaEffectThetadegrees/=(1.09833e+00+1.70000e-02*(thetawires*kRaddeg)); // The gap is different (4mm)
 372       yAngleEffect=1.e-04*gRandom->Gaus(0,sigmaEffectThetadegrees); // Error due to the angle effect in cm
 373     }
 374     }
 375
 376     // Detection elements ids
 377     const AliMUONGeometryModule* kGeometryModule
 378       = GetGeometry()->GetModule(iChamber-1);
 379
 380     AliMUONGeometryDetElement* detElement
 381       = kGeometryModule->FindBySensitiveVolume(CurrentVolumePath());
 382
 383     Int_t detElemId = 0;
 384     if (detElement) detElemId = detElement->GetUniqueID();
 385
 386     if (!detElemId) {
 387       cerr << "Chamber id: "
 388            << setw(3) << iChamber << "  "
 389            << "Current SV: "
 390            <<  CurrentVolumePath()
 391            << "  detElemId: "
 392            << setw(5) << detElemId
 393            << endl;
 394       Double_t x, y, z;
 395       gMC->TrackPosition(x, y, z);
 396       cerr << "   global position: "
 397            << x << ", " << y << ", " << z
 398            << endl;
 399       AliWarning("DetElemId not identified.");
 400     }
 401
 402     // One hit per chamber
 403     GetMUONData()->AddHit(fIshunt,
 404                           gAlice->GetMCApp()->GetCurrentTrackNumber(),
 405                           iChamber, ipart,
 406                           fTrackPosition.X(),
 407                           fTrackPosition.Y()+yAngleEffect,
 408                           fTrackPosition.Z(),
 409                           gMC->TrackTime(),
 410                           fTrackMomentum.P(),
 411                           theta,
 412                           phi,
 413                           fStepSum[idvol],
 414                           fDestepSum[idvol],
 415                           fTrackPosition.X(),
 416                           fTrackPosition.Y(),
 417                           fTrackPosition.Z());
 418
 419 //     if (AliDebugLevel()){
 420 //       AliDebug(1,Form("Exit: Particle exiting from chamber %d",iChamber));
 421 //       AliDebug(1,Form("Exit: StepSum %f eloss geant %g ",fStepSum[idvol],fDestepSum[idvol]));
 422 //       AliDebug(1,Form("Exit: Track Position %f %f %f",fTrackPosition.X(),fTrackPosition.Y(),fTrackPosition.Z())) ;
 423 //     }
 424     fStepSum[idvol]  =0; // Reset for the next event
 425     fDestepSum[idvol]=0; // Reset for the next event
 426   }
 427 }
 428
 429 //_______________________________________________________________________________
 430 void AliMUONv1::StepManager2()
 431 {
 432   // Stepmanager for the chambers
 433
 434   // Only charged tracks
 435   if( !(gMC->TrackCharge()) ) return;
 436   // Only charged tracks
 437
 438   // Only gas gap inside chamber
 439   // Tag chambers and record hits when track enters
 440   static Int_t   idvol=-1;
 441   Int_t   iChamber=0;
 442   Int_t   id=0;
 443   Int_t   copy;
 444   const  Float_t kBig = 1.e10;
 445
 446
 447   //
 448   // Only gas gap inside chamber
 449   // Tag chambers and record hits when track enters
 450   id=gMC->CurrentVolID(copy);
 451   iChamber = GetChamberId(id);
 452   idvol = iChamber -1;
 453
 454   if (idvol == -1) return;
 455
 456   // Filling TrackRefs file for MUON. Our Track references are the active volume of the chambers
 457   if ( (gMC->IsTrackEntering() || gMC->IsTrackExiting() ) )
 458     AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber());
 459
 460    if( gMC->IsTrackEntering() ) {
 461      Float_t theta = fTrackMomentum.Theta();
 462      if ((TMath::Pi()-theta)*kRaddeg>=15.) gMC->SetMaxStep(fStepMaxInActiveGas); // We use Pi-theta because z is negative
 463   }
 464
 465    //   AliDebug(1,
 466    //       Form("Active volume found %d chamber %d Z chamber is %f ",idvol,iChamber,
 467    //            ( (AliMUONChamber*)(*fChambers)[idvol])->Z())) ;
 468   // Particule id and mass,
 469   Int_t     ipart = gMC->TrackPid();
 470   Float_t   mass  = gMC->TrackMass();
 471
 472   fDestepSum[idvol]+=gMC->Edep();
 473   // Get current particle id (ipart), track position (pos)  and momentum (mom)
 474   if ( fStepSum[idvol]==0.0 )  gMC->TrackMomentum(fTrackMomentum);
 475   fStepSum[idvol]+=gMC->TrackStep();
 476
 477   //  if (AliDebugLevel()) {
 478   //   AliDebug(1,Form("Step, iChamber %d, Particle %d, theta %f phi %f mass %f StepSum %f eloss %g",
 479   //                 iChamber,ipart, fTrackMomentum.Theta()*kRaddeg, fTrackMomentum.Phi()*kRaddeg,
 480   //         mass, fStepSum[idvol], gMC->Edep()));
 481   // AliDebug(1,Form("Step:Track Momentum %f %f %f", fTrackMomentum.X(), fTrackMomentum.Y(),
 482   //         fTrackMomentum.Z()));
 483   // gMC->TrackPosition(fTrackPosition);
 484   // AliDebug(1,Form("Step: Track Position %f %f %f",fTrackPosition.X(),
 485   //         fTrackPosition.Y(),fTrackPosition.Z())) ;
 486   //}
 487
 488   // Track left chamber or StepSum larger than fStepMaxInActiveGas
 489   if ( gMC->IsTrackExiting() ||
 490        gMC->IsTrackStop() ||
 491        gMC->IsTrackDisappeared()||
 492        (fStepSum[idvol]>fStepMaxInActiveGas) ) {
 493
 494     if   ( gMC->IsTrackExiting() ||
 495            gMC->IsTrackStop() ||
 496            gMC->IsTrackDisappeared() ) gMC->SetMaxStep(kBig);
 497
 498     gMC->TrackPosition(fTrackPosition);
 499     Float_t theta = fTrackMomentum.Theta();
 500     Float_t phi   = fTrackMomentum.Phi();
 501
 502     TLorentzVector backToWire( fStepSum[idvol]/2.*sin(theta)*cos(phi),
 503                                fStepSum[idvol]/2.*sin(theta)*sin(phi),
 504                                fStepSum[idvol]/2.*cos(theta),0.0       );
 505     //    AliDebug(1,
 506     //       Form("Track Position %f %f %f",fTrackPosition.X(),fTrackPosition.Y(),fTrackPosition.Z()));
 507     // AliDebug(1,
 508     //       Form("Exit: Track backToWire %f %f %f",backToWire.X(),backToWire.Y(),backToWire.Z())) ;
 509     fTrackPosition-=backToWire;
 510
 511     //-------------- Angle effect
 512     // Ratio between energy loss of particle and Mip as a function of BetaGamma of particle (Energy/Mass)
 513
 514     Float_t betaxGamma    = fTrackMomentum.P()/mass;//  pc/mc2
 515     Float_t sigmaEffect10degrees;
 516     Float_t sigmaEffectThetadegrees;
 517     Float_t eLossParticleELossMip;
 518     Float_t yAngleEffect=0.;
 519     Float_t thetawires      =  TMath::Abs( TMath::ASin( TMath::Sin(TMath::Pi()-theta) * TMath::Sin(phi) ) );// We use Pi-theta because z is negative
 520
 521
 522     if (fAngleEffect){
 523     if ( (betaxGamma >3.2)   &&  (thetawires*kRaddeg<=15.) ) {
 524       betaxGamma=TMath::Log(betaxGamma);
 525       eLossParticleELossMip = fElossRatio->Eval(betaxGamma);
 526       // 10 degrees is a reference for a model (arbitrary)
 527       sigmaEffect10degrees=fAngleEffect10->Eval(eLossParticleELossMip);// in micrometers
 528       // Angle with respect to the wires assuming that chambers are perpendicular to the z axis.
 529       sigmaEffectThetadegrees =  sigmaEffect10degrees/fAngleEffectNorma->Eval(thetawires*kRaddeg);  // For 5mm gap
 530       if ( (iChamber==1)  ||  (iChamber==2) )
 531         sigmaEffectThetadegrees/=(1.09833e+00+1.70000e-02*(thetawires*kRaddeg)); // The gap is different (4mm)
 532       yAngleEffect=1.e-04*gRandom->Gaus(0,sigmaEffectThetadegrees); // Error due to the angle effect in cm
 533     }
 534     }
 535
 536     // Detection elements ids
 537     const AliMUONGeometryModule* kGeometryModule
 538       = GetGeometry()->GetModule(iChamber-1);
 539
 540     AliMUONGeometryDetElement* detElement
 541       = kGeometryModule->FindBySensitiveVolume(CurrentVolumePath());
 542
 543     Int_t detElemId = 0;
 544     if (detElement) detElemId = detElement->GetUniqueID();
 545
 546     if (!detElemId) {
 547       cerr << "Chamber id: "
 548            << setw(3) << iChamber << "  "
 549            << "Current SV: "
 550            <<  CurrentVolumePath()
 551            << "  detElemId: "
 552            << setw(5) << detElemId
 553            << endl;
 554       Double_t x, y, z;
 555       gMC->TrackPosition(x, y, z);
 556       cerr << "   global position: "
 557            << x << ", " << y << ", " << z
 558            << endl;
 559       AliError("DetElemId not identified.");
 560     }
 561
 562     // One hit per chamber
 563     GetMUONData()->AddHit2(fIshunt,
 564                           gAlice->GetMCApp()->GetCurrentTrackNumber(),
 565                           detElemId, ipart,
 566                           fTrackPosition.X(),
 567                           fTrackPosition.Y()+yAngleEffect,
 568                           fTrackPosition.Z(),
 569                           gMC->TrackTime(),
 570                           fTrackMomentum.P(),
 571                           theta,
 572                           phi,
 573                           fStepSum[idvol],
 574                           fDestepSum[idvol],
 575                           fTrackPosition.X(),
 576                           fTrackPosition.Y(),
 577                           fTrackPosition.Z());
 578
 579     //       AliDebug(1,Form("Exit: Particle exiting from chamber %d",iChamber));
 580     //       AliDebug(1,Form("Exit: StepSum %f eloss geant %g ",fStepSum[idvol],fDestepSum[idvol]));
 581     //       AliDebug(1,Form("Exit: Track Position %f %f %f",fTrackPosition.X(),fTrackPosition.Y(),fTrackPosition.Z()) ;
 582
 583     fStepSum[idvol]  =0; // Reset for the next event
 584     fDestepSum[idvol]=0; // Reset for the next event
 585   }
 586 }