MUON/AliMUONHit.cxx

   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  **************************************************************************/
  15
  16 /* $Id$ */
  17
  18 // MUON class for MonteCarlo Hits, inherited from AliHit for the
  19 // In addition to the ALiHit data member fX, fY, fZ and fTrack, AliMUONHit contains some info about the particle crossing the chamber:
  20 // Impulsion: fPtot, fPx, fPy and fPz
  21 // Reference position at the center of the chamber (wire plane) fXref, fYref and fZref
  22 // Cumulated path along the active volume fTlength for spliting of hits for very inclined tracks
  23 // Energy loss of the particle inside the gas active volume.
  24 // Incident fTheta and fPhi angle with respect of the wire plane of the chamber.
  25 //
  26
  27 #include "AliMUONHit.h"
  28 #include "AliMpDEManager.h"
  29
  30 #include "AliLog.h"
  31
  32 #include <Riostream.h>
  33 #include <TMath.h>
  34 #include <TString.h>
  35
  36 /// \cond CLASSIMP
  37 ClassImp(AliMUONHit)
  38 /// \endcond
  39
  40 //___________________________________________
  41 AliMUONHit::AliMUONHit()
  42   : AliHit(),
  43     fDetElemId(0),
  44     fParticle(0),
  45     fTheta(0),
  46     fPhi(0),
  47     fTlength(0),
  48     fEloss(0),
  49     fAge(0),
  50     fPHfirst(0),
  51     fPHlast(0),
  52     fPTot(0),
  53     fPx(0),
  54     fPy(0),
  55     fPz(0),
  56     fXref(0),
  57     fYref(0),
  58     fZref(0)
  59 {
  60 /// Default constructor
  61 }
  62
  63 //___________________________________________
  64 AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits)
  65   : AliHit(shunt, track),
  66     fDetElemId(vol[0]),
  67     fParticle(hits[0]),
  68     fTheta(hits[4]),
  69     fPhi(hits[5]),
  70     fTlength(hits[6]),
  71     fEloss(hits[7]),
  72     fAge(hits[14]),
  73     fPHfirst((Int_t)hits[8]),
  74     fPHlast((Int_t)hits[9]),
  75     fPTot(hits[10]),
  76     fPx(hits[11]),
  77     fPy(hits[12]),
  78     fPz(hits[13]),
  79     fXref(0),
  80     fYref(0),
  81     fZref(0)
  82 {
  83 /// Constructor
  84
  85     fX         = hits[1];
  86     fY         = hits[2];
  87     fZ         = hits[3];
  88 }
  89
  90 //___________________________________________
  91 AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart,
  92                        Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
  93                        Float_t theta, Float_t phi, Float_t length, Float_t destep)
  94   : AliHit(shunt, track),
  95     fDetElemId(detElemId),
  96     fParticle(idpart),
  97     fTheta(theta),
  98     fPhi(phi),
  99     fTlength(length),
 100     fEloss(destep),
 101     fAge(tof),
 102     fPHfirst(0),
 103     fPHlast(0),
 104     fPTot(momentum),
 105     fPx(momentum * TMath::Sin(theta) * TMath::Cos(phi)),
 106     fPy(momentum * TMath::Sin(theta) * TMath::Sin(phi)),
 107     fPz(momentum * TMath::Cos(theta)),
 108     fXref(0),
 109     fYref(0),
 110     fZref(0)
 111 {
 112 /// Constructor
 113     fX         = X;
 114     fY         = Y;
 115     fZ         = Z;
 116 }
 117
 118 //-----------------------------------------------------------------------------------------------
 119 AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart,
 120                        Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
 121                        Float_t theta, Float_t phi, Float_t length, Float_t destep,
 122                        Float_t Xref,Float_t Yref,Float_t Zref)
 123   : AliHit(shunt, track),
 124     fDetElemId(detElemId),
 125     fParticle(idpart),
 126     fTheta(theta),
 127     fPhi(phi),
 128     fTlength(length),
 129     fEloss(destep),
 130     fAge(tof),
 131     fPHfirst(0),
 132     fPHlast(0),
 133     fPTot(momentum),
 134     fPx(momentum * TMath::Sin(theta) * TMath::Cos(phi)),
 135     fPy(momentum * TMath::Sin(theta) * TMath::Sin(phi)),
 136     fPz(momentum * TMath::Cos(theta)),
 137     fXref(Xref),
 138     fYref(Yref),
 139     fZref(Zref)
 140 {
 141 /// Constructor
 142
 143     fX         = X;
 144     fY         = Y;
 145     fZ         = Z;
 146 }
 147
 148 //-----------------------------------------------------------------------------------------------
 149 AliMUONHit::~AliMUONHit()
 150 {
 151 /// Dectructor
 152 }
 153
 154 //-----------------------------------------------------------------------------------------------
 155 Int_t  AliMUONHit::Chamber()  const
 156 {
 157 /// Return chamber ID
 158
 159   return AliMpDEManager::GetChamberId(fDetElemId) + 1;
 160 }
 161
 162 //-----------------------------------------------------------------------------------------------
 163 void AliMUONHit::Print(Option_t* opt) const
 164 {
 165 /// Printing hit information
 166 /// "full" option for printing all the information about the hit
 167
 168   TString sopt(opt);
 169   sopt.ToUpper();
 170
 171   if ( sopt.Contains("FULL") ) {
 172     cout <<"<AliMUONHit>: Geant track="   << setw(4)  << Track() <<
 173       ", DetEle="        << setw(4)  << DetElemId() <<
 174       ", (x,y,z)=(" << setw(7) << setprecision(5) << X() << "," << setw(7) << setprecision(5) << Y() <<  "," << setw(7) << setprecision(5) << Z() <<
 175       " )cm, Delta E=" << setw(8) << setprecision(3) << Eloss() << " GeV" << endl;
 176   }
 177   else {
 178     cout << "<AliMUONHit>: DetEle="        << setw(4)  << DetElemId() <<
 179       ", (x,y,z)=(" << setw(7) << setprecision(5) << X() << "," << setw(7) << setprecision(5) << Y() <<  "," << setw(7) << setprecision(5) << Z() <<
 180       " ) cm" <<endl;
 181   }
 182
 183 }