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