1 ///////////////////////////////////////////////////////////////////////////////
3 // Time Of Flight: as for version 1 but full coverage //
4 // This class contains the functions for version 3 of the Time Of Flight //
9 <img src="picts/AliTOFv3Class.gif">
13 ///////////////////////////////////////////////////////////////////////////////
21 //_____________________________________________________________________________
25 // Default constructor
29 //_____________________________________________________________________________
30 AliTOFv3::AliTOFv3(const char *name, const char *title)
34 // Standard constructor
38 //_____________________________________________________________________________
39 void AliTOFv3::CreateGeometry()
42 // Create geometry for Time Of Flight version 0
46 <img src="picts/AliTOFv3.gif">
51 // Create common geometry
53 AliTOF::CreateGeometry();
56 //_____________________________________________________________________________
57 void AliTOFv3::TOFpc(Float_t xm, Float_t ym, Float_t zm0,
58 Float_t zm1, Float_t zm2)
61 // Definition of the Time Of Fligh Resistive Plate Chambers
62 // xm, ym, zm - sizes of TOF modules (large)
65 Float_t zazor, xp, yp, zp;
68 Int_t *idtmed = fIdtmed->GetArray()-499;
72 // Sizes of RPC chamber
73 xp = 3.0; //small pixel
74 //xp = 3.9; //large pixel
75 yp = 12.3*0.05; // 5% X0 of glass
76 zp = 3.0; //small pixel
77 //zp = 4.1; //large pixel
78 // Large not sensitive volumes with CO2
82 gMC->Gsvolu("FBT1", "BOX ", idtmed[506], par, 3); // CO2
83 gMC->Gspos("FBT1", 0, "FTO1", 0., 0., 0., 0, "ONLY");
84 gMC->Gsdvn("FDT1", "FBT1", 2, 3); // 2 large modules along Z
86 gMC->Gsvolu("FBT2", "BOX ", idtmed[506], par, 3); // CO2
87 gMC->Gspos("FBT2", 1, "FTO2", 0., 0., 0., 0, "ONLY");
88 gMC->Gsdvn("FDT2", "FBT2", 2, 3); // 2 (PHOS) modules along Z
90 gMC->Gsvolu("FBT3", "BOX ", idtmed[506], par, 3); // CO2
91 gMC->Gspos("FBT3", 2, "FTO3", 0., 0., 0., 0, "ONLY");
92 gMC->Gsdvn("FDT3", "FBT3", 1, 3); // 1 (RICH) module along Z
94 // subtraction of dead boundaries in X=2 cm and Z=7/2 cm
96 Int_t nz0, nz1, nz2, nx; //- numbers of pixels
97 nx = Int_t (par[0]*2/xp);
98 cout <<"************************* TOF geometry **************************"<<endl;
99 cout<< "nx = "<< nx << " x size = "<< par[0]*2/nx << endl;
101 par[2] = (zm0 / 2.)/2.; //this is half size of module after division by 2
103 nz0 = Int_t (par[2]*2/zp);
104 cout<< "nz0 = "<< nz0 << " z0 size = "<< par[2]*2/nz0 << endl;
105 gMC->Gsvolu("FLT1", "BOX ", idtmed[506], par, 3); // CO2
106 gMC->Gspos("FLT1", 0, "FDT1", 0., 0., 0., 0, "ONLY");
107 par[2] = (zm1 / 2.)/2.; //this is half size of module after division by 2
109 nz1 = Int_t (par[2]*2/zp);
110 cout<< "nz1 = "<< nz1 << " z1 size = "<< par[2]*2/nz1 << endl;
111 gMC->Gsvolu("FLT2", "BOX ", idtmed[506], par, 3); // CO2
112 gMC->Gspos("FLT2", 0, "FDT2", 0., 0., 0., 0, "ONLY");
113 par[2] = (zm2 / 2.); //this is half size of module after division by 1
115 nz2 = Int_t (par[2]*2/zp);
116 cout<< "nz2 = "<< nz2 << " z2 size = "<< par[2]*2/nz2 << endl;
117 gMC->Gsvolu("FLT3", "BOX ", idtmed[506], par, 3); // CO2
118 gMC->Gspos("FLT3", 0, "FDT3", 0., 0., 0., 0, "ONLY");
120 ////////// Layers before detector ////////////////////
121 // Mylar layer in front 0.5mm thick at the beginning
125 ycoor = -ym/2 + par[1];
126 gMC->Gsvolu("FMY1", "BOX ", idtmed[511], par, 3); // Mylar
127 gMC->Gspos("FMY1", 0, "FLT1", 0., ycoor, 0., 0, "ONLY");
128 gMC->Gsvolu("FMY2", "BOX ", idtmed[511], par, 3); // Mylar
129 gMC->Gspos("FMY2", 0, "FLT2", 0., ycoor, 0., 0, "ONLY");
130 gMC->Gsvolu("FMY3", "BOX ", idtmed[511], par, 3); // Mylar
131 gMC->Gspos("FMY3", 0, "FLT3", 0., ycoor, 0., 0, "ONLY");
132 // Honeycomb layer (1cm of special!!! polyethilene)
133 ycoor = ycoor + par[1];
137 ycoor = ycoor + par[1];
138 gMC->Gsvolu("FPL1", "BOX ", idtmed[503], par, 3); // Hony
139 gMC->Gspos("FPL1", 0, "FLT1", 0., ycoor, 0., 0, "ONLY");
140 gMC->Gsvolu("FPL2", "BOX ", idtmed[503], par, 3); // Hony
141 gMC->Gspos("FPL2", 0, "FLT2", 0., ycoor, 0., 0, "ONLY");
142 gMC->Gsvolu("FPL3", "BOX ", idtmed[503], par, 3); // Hony
143 gMC->Gspos("FPL3", 0, "FLT3", 0., ycoor, 0., 0, "ONLY");
145 ///////////////// Detector itself //////////////////////
147 par[1] = yp/2; // 5 %X0 thick of glass
150 gMC->Gsvolu("FLD1", "BOX ", idtmed[514], par, 3); // Glass
151 gMC->Gspos("FLD1", 0, "FLT1", 0., ycoor, 0., 0, "ONLY");
152 gMC->Gsvolu("FLD2", "BOX ", idtmed[514], par, 3); // Glass
153 gMC->Gspos("FLD2", 0, "FLT2", 0., ycoor, 0., 0, "ONLY");
154 gMC->Gsvolu("FLD3", "BOX ", idtmed[514], par, 3); // Glass
155 gMC->Gspos("FLD3", 0, "FLT3", 0., ycoor, 0., 0, "ONLY");
157 gMC->Gsdvn("FLZ1", "FLD1", nz0, 3); //pixel size xp=zp=3
158 gMC->Gsdvn("FLZ2", "FLD2", nz1, 3);
159 gMC->Gsdvn("FLZ3", "FLD3", nz2, 3);
160 gMC->Gsdvn("FLX1", "FLZ1", nx, 1);
161 gMC->Gsdvn("FLX2", "FLZ2", nx, 1);
162 gMC->Gsdvn("FLX3", "FLZ3", nx, 1);
165 par[1] = -1;//yp/2; // 5 %X0 thick of glass
167 gMC->Gsvolu("FPA0", "BOX ", idtmed[514], par, 3);// Glass
168 gMC->Gspos("FPA0", 1, "FLX1", 0., 0., 0., 0, "ONLY");
169 gMC->Gspos("FPA0", 2, "FLX2", 0., 0., 0., 0, "ONLY");
170 gMC->Gspos("FPA0", 3, "FLX3", 0., 0., 0., 0, "ONLY");
171 // Freon gas sencitive volume
175 gMC->Gsvolu("FPG0", "BOX ", idtmed[513], par, 3);// Freon
176 gMC->Gspos("FPG0", 0, "FPA0", 0., 0., 0., 0, "ONLY");
178 ////////// Layers after detector ////////////////////
179 // Honeycomb layer after (3cm)
183 ycoor = -ym/2 + 6. - par[1];
184 gMC->Gsvolu("FPE1", "BOX ", idtmed[503], par, 3); // Hony
185 gMC->Gspos("FPE1", 0, "FLT1", 0., ycoor, 0., 0, "ONLY");
186 gMC->Gsvolu("FPE2", "BOX ", idtmed[503], par, 3); // Hony
187 gMC->Gspos("FPE2", 0, "FLT2", 0., ycoor, 0., 0, "ONLY");
188 gMC->Gsvolu("FPE3", "BOX ", idtmed[503], par, 3); // Hony
189 gMC->Gspos("FPE3", 0, "FLT3", 0., ycoor, 0., 0, "ONLY");
190 // Electronics (Cu) after
192 par[1] = 1.43*0.05 / 2.; // 5% of X0
194 ycoor = -ym/2 + 6.+par[1];
195 gMC->Gsvolu("FEC1", "BOX ", idtmed[501], par, 3); // Cu
196 gMC->Gspos("FEC1", 0, "FLT1", 0., ycoor, 0., 0, "ONLY");
197 gMC->Gsvolu("FEC2", "BOX ", idtmed[501], par, 3); // Cu
198 gMC->Gspos("FEC2", 0, "FLT2", 0., ycoor, 0., 0, "ONLY");
199 gMC->Gsvolu("FEC3", "BOX ", idtmed[501], par, 3); // Cu
200 gMC->Gspos("FEC3", 0, "FLT3", 0., ycoor, 0., 0, "ONLY");
201 // Cooling water after
202 ycoor = ycoor+par[1];
204 par[1] = 36.1*0.02 / 2.; // 2% of X0
206 ycoor = ycoor+par[1];
207 gMC->Gsvolu("FWA1", "BOX ", idtmed[515], par, 3); // Water
208 gMC->Gspos("FWA1", 0, "FLT1", 0., ycoor, 0., 0, "ONLY");
209 gMC->Gsvolu("FWA2", "BOX ", idtmed[515], par, 3); // Water
210 gMC->Gspos("FWA2", 0, "FLT2", 0., ycoor, 0., 0, "ONLY");
211 gMC->Gsvolu("FWA3", "BOX ", idtmed[515], par, 3); // Water
212 gMC->Gspos("FWA3", 0, "FLT3", 0., ycoor, 0., 0, "ONLY");
213 //back plate honycomb (2cm)
217 ycoor = ym/2 - par[1];
218 gMC->Gsvolu("FEG1", "BOX ", idtmed[503], par, 3); // Hony
219 gMC->Gspos("FEG1", 0, "FLT1", 0., ycoor, 0., 0, "ONLY");
220 gMC->Gsvolu("FEG2", "BOX ", idtmed[503], par, 3); // Hony
221 gMC->Gspos("FEG2", 0, "FLT2", 0., ycoor, 0., 0, "ONLY");
222 gMC->Gsvolu("FEG3", "BOX ", idtmed[503], par, 3); // Hony
223 gMC->Gspos("FEG3", 0, "FLT3", 0., ycoor, 0., 0, "ONLY");
226 //_____________________________________________________________________________
227 void AliTOFv3::DrawModule()
230 // Draw a shaded view of the Time Of Flight version 0
232 // Set everything unseen
233 gMC->Gsatt("*", "seen", -1);
235 // Set ALIC mother transparent
236 gMC->Gsatt("ALIC","SEEN",0);
238 // Set the volumes visible
239 gMC->Gsatt("ALIC","SEEN",0);
240 gMC->Gsatt("FBAR","SEEN",1);
241 gMC->Gsatt("FTO1","SEEN",1);
242 gMC->Gsatt("FTO2","SEEN",1);
243 gMC->Gsatt("FTO3","SEEN",1);
244 gMC->Gsatt("FBT1","SEEN",1);
245 gMC->Gsatt("FBT2","SEEN",1);
246 gMC->Gsatt("FBT3","SEEN",1);
247 gMC->Gsatt("FDT1","SEEN",1);
248 gMC->Gsatt("FDT2","SEEN",1);
249 gMC->Gsatt("FDT3","SEEN",1);
250 gMC->Gsatt("FLT1","SEEN",1);
251 gMC->Gsatt("FLT2","SEEN",1);
252 gMC->Gsatt("FLT3","SEEN",1);
254 gMC->Gdopt("hide", "on");
255 gMC->Gdopt("shad", "on");
256 gMC->Gsatt("*", "fill", 7);
257 gMC->SetClipBox(".");
258 gMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
260 gMC->Gdraw("alic", 40, 30, 0, 12, 9.5, .02, .02);
261 gMC->Gdhead(1111, "Time Of Flight");
262 gMC->Gdman(18, 4, "MAN");
263 gMC->Gdopt("hide","off");
266 //_____________________________________________________________________________
267 void AliTOFv3::CreateMaterials()
270 // Define materials for the Time Of Flight
272 AliTOF::CreateMaterials();
275 //_____________________________________________________________________________
276 void AliTOFv3::Init()
279 // Initialise the detector after the geometry has been defined
282 fIdFTO2=gMC->VolId("FTO2");
283 fIdFTO3=gMC->VolId("FTO3");
284 fIdFLT1=gMC->VolId("FLT1");
285 fIdFLT2=gMC->VolId("FLT2");
286 fIdFLT3=gMC->VolId("FLT3");
289 //_____________________________________________________________________________
290 void AliTOFv3::StepManager()
293 // Procedure called at each step in the Time Of Flight
295 TLorentzVector mom, pos;
299 Int_t *idtmed = fIdtmed->GetArray()-499;
300 if(gMC->GetMedium()==idtmed[514-1] &&
301 gMC->IsTrackEntering() && gMC->TrackCharge()
302 && gMC->CurrentVolID(copy)==fIdSens) {
303 TClonesArray &lhits = *fHits;
305 // Record only charged tracks at entrance
306 gMC->CurrentVolOffID(1,copy);
308 gMC->CurrentVolOffID(3,copy);
310 id=gMC->CurrentVolOffID(8,copy);
314 id=gMC->CurrentVolOffID(5,copy);
315 if(id==fIdFLT3) vol[1]+=6;
316 } else if (id==fIdFTO2) {
318 id=gMC->CurrentVolOffID(5,copy);
319 if(id==fIdFLT2) vol[1]+=8;
321 id=gMC->CurrentVolOffID(5,copy);
322 if(id==fIdFLT1) vol[1]+=14;
324 gMC->TrackPosition(pos);
325 gMC->TrackMomentum(mom);
327 Double_t ptot=mom.Rho();
328 Double_t norm=1/ptot;
331 hits[i+3]=mom[i]*norm;
335 new(lhits[fNhits++]) AliTOFhit(fIshunt,gAlice->CurrentTrack(),vol,hits);