]>
Commit | Line | Data |
---|---|---|
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 | } |