]>
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 | #include "TVirtualMC.h" | |
19 | #include "TParticle.h" | |
20 | ||
21 | #include "AliTrackReference.h" | |
22 | #include <Riostream.h> | |
23 | ||
24 | // | |
25 | // Track Reference object is created every time particle is | |
26 | // crossing detector bounds. The object is created by Step Manager | |
27 | // | |
28 | // The class stores the following informations: | |
29 | // track label, | |
30 | // track position: X,Y,X | |
31 | // track momentum px, py, pz | |
32 | // track length and time of fligth: both in cm | |
33 | // status bits from Monte Carlo | |
34 | // | |
35 | ||
36 | ||
37 | ClassImp(AliTrackReference) | |
38 | ||
39 | //_______________________________________________________________________ | |
40 | AliTrackReference::AliTrackReference(): | |
41 | TObject(), | |
42 | fTrack(0), | |
43 | fX(0), | |
44 | fY(0), | |
45 | fZ(0), | |
46 | fPx(0), | |
47 | fPy(0), | |
48 | fPz(0), | |
49 | fLength(0), | |
50 | fTime(0), | |
51 | fUserId(0), | |
52 | fDetectorId(-999) | |
53 | { | |
54 | // | |
55 | // Default constructor | |
56 | // Creates empty object | |
57 | ||
58 | for(Int_t i=0; i<16; i++) ResetBit(BIT(i)); | |
59 | } | |
60 | ||
61 | AliTrackReference::AliTrackReference(const AliTrackReference &tr) : | |
62 | TObject(), | |
63 | fTrack(tr.fTrack), | |
64 | fX(tr.fX), | |
65 | fY(tr.fY), | |
66 | fZ(tr.fZ), | |
67 | fPx(tr.fPx), | |
68 | fPy(tr.fPy), | |
69 | fPz(tr.fPz), | |
70 | fLength(tr.fLength), | |
71 | fTime(tr.fTime), | |
72 | fUserId(tr.fUserId), | |
73 | fDetectorId(tr.fDetectorId) | |
74 | { | |
75 | // Copy Constructor | |
76 | } | |
77 | ||
78 | //_______________________________________________________________________ | |
79 | AliTrackReference::AliTrackReference(Int_t label, Int_t id) : | |
80 | TObject(), | |
81 | fTrack(label), | |
82 | fX(0), | |
83 | fY(0), | |
84 | fZ(0), | |
85 | fPx(0), | |
86 | fPy(0), | |
87 | fPz(0), | |
88 | fLength(gMC->TrackLength()), | |
89 | fTime(gMC->TrackTime()), | |
90 | fUserId(0), | |
91 | fDetectorId(id) | |
92 | { | |
93 | // | |
94 | // Create Reference object out of label and | |
95 | // data in TVirtualMC object | |
96 | // | |
97 | // Creates an object and fill all parameters | |
98 | // from data in VirtualMC | |
99 | // | |
100 | // Sylwester Radomski, (S.Radomski@gsi.de) | |
101 | // GSI, Jan 31, 2003 | |
102 | // | |
103 | ||
104 | Double_t vec[4]; | |
105 | ||
106 | gMC->TrackPosition(vec[0],vec[1],vec[2]); | |
107 | ||
108 | fX = vec[0]; | |
109 | fY = vec[1]; | |
110 | fZ = vec[2]; | |
111 | ||
112 | gMC->TrackMomentum(vec[0],vec[1],vec[2],vec[3]); | |
113 | ||
114 | fPx = vec[0]; | |
115 | fPy = vec[1]; | |
116 | fPz = vec[2]; | |
117 | ||
118 | // Set Up status code | |
119 | // Copy Bits from virtual MC | |
120 | ||
121 | for(Int_t i=0; i<16; i++) ResetBit(BIT(i)); | |
122 | ||
123 | SetBit(BIT(0), gMC->IsNewTrack()); | |
124 | SetBit(BIT(1), gMC->IsTrackAlive()); | |
125 | SetBit(BIT(2), gMC->IsTrackDisappeared()); | |
126 | SetBit(BIT(3), gMC->IsTrackEntering()); | |
127 | SetBit(BIT(4), gMC->IsTrackExiting()); | |
128 | SetBit(BIT(5), gMC->IsTrackInside()); | |
129 | SetBit(BIT(6), gMC->IsTrackOut()); | |
130 | SetBit(BIT(7), gMC->IsTrackStop()); | |
131 | // | |
132 | // This particle has to be kept | |
133 | ||
134 | } | |
135 | /* | |
136 | AliExternalTrackParam * AliTrackReference::MakeTrack(const AliTrackReference *ref, Double_t mass) | |
137 | { | |
138 | // | |
139 | // Make dummy track from the track reference | |
140 | // negative mass means opposite charge | |
141 | // | |
142 | Double_t xx[5]; | |
143 | Double_t cc[15]; | |
144 | for (Int_t i=0;i<15;i++) cc[i]=0; | |
145 | Double_t x = ref->X(), y = ref->Y(), z = ref->Z(); | |
146 | Double_t alpha = TMath::ATan2(y,x); | |
147 | Double_t xr = TMath::Sqrt(x*x+y*y); | |
148 | xx[0] = ref->LocalY(); | |
149 | xx[1] = z; | |
150 | xx[3] = ref->Pz()/ref->Pt(); | |
151 | xx[4] = 1./ref->Pt(); | |
152 | if (mass<0) xx[4]*=-1.; // negative mass - negative direction | |
153 | Double_t alphap = TMath::ATan2(ref->Py(),ref->Px())-alpha; | |
154 | if (alphap> TMath::Pi()) alphap-=TMath::Pi(); | |
155 | if (alphap<-TMath::Pi()) alphap+=TMath::Pi(); | |
156 | xx[2] = TMath::Sin(alphap); | |
157 | ||
158 | AliExternalTrackParam * track = new AliExternalTrackParam(xr,alpha,xx,cc); | |
159 | return track; | |
160 | } | |
161 | */ | |
162 | //_______________________________________________________________________ | |
163 | void | |
164 | AliTrackReference::Print(Option_t* /*opt*/) const | |
165 | { | |
166 | cout << Form("Label %d P=%7.2f (PX,PY,PZ)=(%7.2f,%7.2f,%7.2f) (X,Y,Z)=(%7.2f,%7.2f,%7.2f)" | |
167 | " Length=%7.2f Time=%7.2f UserId=%d", | |
168 | Label(),P(),Px(),Py(),Pz(),X(),Y(),Z(),GetLength(),GetTime(),UserId()) << endl; | |
169 | } | |
170 |