]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - ZDC/AliZDCv3.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (non-incremental) | history | HEAD
Reverting...I was too fast
[u/mrichter/AliRoot.git] / ZDC / AliZDCv3.cxx
... / ...
CommitLineData
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
17///////////////////////////////////////////////////////////////////////
18// //
19// AliZDCv3 --- new ZDC geometry //
20// with both ZDC arms geometry implemented //
21// //
22///////////////////////////////////////////////////////////////////////
23
24// --- Standard libraries
25#include "stdio.h"
26
27// --- ROOT system
28#include <TMath.h>
29#include <TRandom.h>
30#include <TSystem.h>
31#include <TTree.h>
32#include <TVirtualMC.h>
33#include <TGeoManager.h>
34#include <TGeoMatrix.h>
35#include <TGeoTube.h>
36#include <TGeoCone.h>
37#include <TGeoShape.h>
38#include <TGeoScaledShape.h>
39#include <TGeoCompositeShape.h>
40#include <TParticle.h>
41
42// --- AliRoot classes
43#include "AliLog.h"
44#include "AliConst.h"
45#include "AliMagF.h"
46#include "AliRun.h"
47#include "AliZDCv3.h"
48#include "AliMC.h"
49
50class AliZDCHit;
51class AliPDG;
52class AliDetector;
53
54
55ClassImp(AliZDCv3)
56
57//_____________________________________________________________________________
58AliZDCv3::AliZDCv3() :
59 AliZDC(),
60 fMedSensF1(0),
61 fMedSensF2(0),
62 fMedSensZP(0),
63 fMedSensZN(0),
64 fMedSensZEM(0),
65 fMedSensGR(0),
66 fMedSensPI(0),
67 fMedSensTDI(0),
68 fMedSensVColl(0),
69 fMedSensLumi(0),
70 fNalfan(0),
71 fNalfap(0),
72 fNben(0),
73 fNbep(0),
74 fZEMLength(0),
75 fpLostITC(0),
76 fpLostD1C(0),
77 fpcVCollC(0),
78 fpDetectedC(0),
79 fnDetectedC(0),
80 fpLostITA(0),
81 fpLostD1A(0),
82 fpLostTDI(0),
83 fpcVCollA(0),
84 fpDetectedA(0),
85 fnDetectedA(0),
86 fVCollSideCAperture(7./2.),
87 fVCollSideCApertureNeg(7./2.),
88 fVCollSideCCentreY(0.),
89 fVCollSideAAperture(7./2.),
90 fVCollSideAApertureNeg(7./2.),
91 fVCollSideACentreY(0.),
92 fTCDDAperturePos(2.0),
93 fTCDDApertureNeg(2.2),
94 fTDIAperturePos(5.5),
95 fTDIApertureNeg(5.5),
96 fLumiLength(15.)
97{
98 //
99 // Default constructor for Zero Degree Calorimeter
100 //
101 for(Int_t i=0; i<3; i++){
102 fDimZN[i] = fDimZP[i] = 0.;
103 fPosZNC[i] = fPosZNA[i] = fPosZPC[i]= fPosZPA[i] = fPosZEM[i] = 0.;
104 fFibZN[i] = fFibZP[i] = 0.;
105 }
106}
107
108//_____________________________________________________________________________
109AliZDCv3::AliZDCv3(const char *name, const char *title) :
110 AliZDC(name,title),
111 fMedSensF1(0),
112 fMedSensF2(0),
113 fMedSensZP(0),
114 fMedSensZN(0),
115 fMedSensZEM(0),
116 fMedSensGR(0),
117 fMedSensPI(0),
118 fMedSensTDI(0),
119 fMedSensVColl(0),
120 fMedSensLumi(0),
121 fNalfan(90),
122 fNalfap(90),
123 fNben(18),
124 fNbep(28),
125 fZEMLength(0),
126 fpLostITC(0),
127 fpLostD1C(0),
128 fpcVCollC(0),
129 fpDetectedC(0),
130 fnDetectedC(0),
131 fpLostITA(0),
132 fpLostD1A(0),
133 fpLostTDI(0),
134 fpcVCollA(0),
135 fpDetectedA(0),
136 fnDetectedA(0),
137 fVCollSideCAperture(7./2.),
138 fVCollSideCApertureNeg(7./2.),
139 fVCollSideCCentreY(0.),
140 fVCollSideAAperture(7./2.),
141 fVCollSideAApertureNeg(7./2.),
142 fVCollSideACentreY(0.),
143 fTCDDAperturePos(2.0),
144 fTCDDApertureNeg(2.2),
145 fTDIAperturePos(5.5),
146 fTDIApertureNeg(5.5),
147 fLumiLength(15.)
148{
149 //
150 // Standard constructor for Zero Degree Calorimeter
151 //
152 //
153 // Check that DIPO, ABSO, DIPO and SHIL is there (otherwise tracking is wrong!!!)
154
155 AliModule* pipe=gAlice->GetModule("PIPE");
156 AliModule* abso=gAlice->GetModule("ABSO");
157 AliModule* dipo=gAlice->GetModule("DIPO");
158 AliModule* shil=gAlice->GetModule("SHIL");
159 if((!pipe) || (!abso) || (!dipo) || (!shil)) {
160 Error("Constructor","ZDC needs PIPE, ABSO, DIPO and SHIL!!!\n");
161 exit(1);
162 }
163 //
164 Int_t ip,jp,kp;
165 for(ip=0; ip<4; ip++){
166 for(kp=0; kp<fNalfap; kp++){
167 for(jp=0; jp<fNbep; jp++){
168 fTablep[ip][kp][jp] = 0;
169 }
170 }
171 }
172 Int_t in,jn,kn;
173 for(in=0; in<4; in++){
174 for(kn=0; kn<fNalfan; kn++){
175 for(jn=0; jn<fNben; jn++){
176 fTablen[in][kn][jn] = 0;
177 }
178 }
179 }
180 //
181 // Parameters for hadronic calorimeters geometry
182 // Positions updated after post-installation measurements
183 fDimZN[0] = 3.52;
184 fDimZN[1] = 3.52;
185 fDimZN[2] = 50.;
186 fDimZP[0] = 11.2;
187 fDimZP[1] = 6.;
188 fDimZP[2] = 75.;
189 fPosZNC[0] = 0.;
190 fPosZNC[1] = 0.;
191 fPosZNC[2] = -11397.3;
192 fPosZPC[0] = 24.35;
193 fPosZPC[1] = 0.;
194 fPosZPC[2] = -11389.3;
195 fPosZNA[0] = 0.;
196 fPosZNA[1] = 0.;
197 fPosZNA[2] = 11395.8;
198 fPosZPA[0] = 24.35;
199 fPosZPA[1] = 0.;
200 fPosZPA[2] = 11387.8;
201 fFibZN[0] = 0.;
202 fFibZN[1] = 0.01825;
203 fFibZN[2] = 50.;
204 fFibZP[0] = 0.;
205 fFibZP[1] = 0.0275;
206 fFibZP[2] = 75.;
207 // Parameters for EM calorimeter geometry
208 fPosZEM[0] = 8.5;
209 fPosZEM[1] = 0.;
210 fPosZEM[2] = 735.;
211 Float_t kDimZEMPb = 0.15*(TMath::Sqrt(2.)); // z-dimension of the Pb slice
212 Float_t kDimZEMAir = 0.001; // scotch
213 Float_t kFibRadZEM = 0.0315; // External fiber radius (including cladding)
214 Int_t kDivZEM[3] = {92, 0, 20}; // Divisions for EM detector
215 Float_t kDimZEM0 = 2*kDivZEM[2]*(kDimZEMPb+kDimZEMAir+kFibRadZEM*(TMath::Sqrt(2.)));
216 fZEMLength = kDimZEM0;
217
218}
219
220//_____________________________________________________________________________
221void AliZDCv3::CreateGeometry()
222{
223 //
224 // Create the geometry for the Zero Degree Calorimeter version 2
225 //* Initialize COMMON block ZDC_CGEOM
226 //*
227
228 CreateBeamLine();
229 CreateZDC();
230}
231
232//_____________________________________________________________________________
233void AliZDCv3::CreateBeamLine()
234{
235 //
236 // Create the beam line elements
237 //
238
239 Double_t zd1, zd2, zCorrDip, zInnTrip, zD1, zD2;
240 Double_t conpar[9], tubpar[3], tubspar[5], boxpar[3];
241
242 //-- rotation matrices for the legs
243 Int_t irotpipe1, irotpipe2;
244 gMC->Matrix(irotpipe1,90.-1.0027,0.,90.,90.,1.0027,180.);
245 gMC->Matrix(irotpipe2,90.+1.0027,0.,90.,90.,1.0027,0.);
246
247 //
248 Int_t *idtmed = fIdtmed->GetArray();
249
250 ////////////////////////////////////////////////////////////////
251 // //
252 // SIDE C - RB26 (dimuon side) //
253 // //
254 ///////////////////////////////////////////////////////////////
255
256
257 // -- Mother of the ZDCs (Vacuum PCON)
258 zd1 = 1921.6;
259
260 conpar[0] = 0.;
261 conpar[1] = 360.;
262 conpar[2] = 2.;
263 conpar[3] = -13500.;
264 conpar[4] = 0.;
265 conpar[5] = 55.;
266 conpar[6] = -zd1;
267 conpar[7] = 0.;
268 conpar[8] = 55.;
269 gMC->Gsvolu("ZDCC", "PCON", idtmed[10], conpar, 9);
270 gMC->Gspos("ZDCC", 1, "ALIC", 0., 0., 0., 0, "ONLY");
271
272
273 // -- BEAM PIPE from compensator dipole to the beginning of D1)
274 tubpar[0] = 6.3/2.;
275 tubpar[1] = 6.7/2.;
276 // From beginning of ZDC volumes to beginning of D1
277 tubpar[2] = (5838.3-zd1)/2.;
278 gMC->Gsvolu("QT01", "TUBE", idtmed[7], tubpar, 3);
279 gMC->Gspos("QT01", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
280 // Ch.debug
281 //printf(" QT01 TUBE pipe from z = %1.2f to z= %1.2f (D1 beg.)\n",-zd1,-2*tubpar[2]-zd1);
282
283 //-- BEAM PIPE from the end of D1 to the beginning of D2)
284
285 //-- FROM MAGNETIC BEGINNING OF D1 TO MAGNETIC END OF D1
286 //-- Cylindrical pipe (r = 3.47) + conical flare
287 // -> Beginning of D1
288 zd1 += 2.*tubpar[2];
289
290 tubpar[0] = 6.94/2.;
291 tubpar[1] = 7.34/2.;
292 tubpar[2] = (6909.8-zd1)/2.;
293 gMC->Gsvolu("QT02", "TUBE", idtmed[7], tubpar, 3);
294 gMC->Gspos("QT02", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
295 // Ch.debug
296 //printf(" QT02 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
297
298 zd1 += 2.*tubpar[2];
299
300 tubpar[0] = 9./2.;
301 tubpar[1] = 9.6/2.;
302 tubpar[2] = (7022.8-zd1)/2.;
303 gMC->Gsvolu("QT03", "TUBE", idtmed[7], tubpar, 3);
304 gMC->Gspos("QT03", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
305 // Ch.debug
306 //printf(" QT03 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
307
308 zd1 += 2.*tubpar[2];
309
310 conpar[0] = 39.2/2.;
311 conpar[1] = 18./2.;
312 conpar[2] = 18.6/2.;
313 conpar[3] = 9./2.;
314 conpar[4] = 9.6/2.;
315 gMC->Gsvolu("QC01", "CONE", idtmed[7], conpar, 5);
316 gMC->Gspos("QC01", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
317 // Ch.debug
318 //printf(" QC01 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
319
320 zd1 += conpar[0] * 2.;
321
322 // ******************************************************
323 // N.B.-> according to last vacuum layout
324 // private communication by D. Macina, mail 27/1/2009
325 // ******************************************************
326 // 2nd section of VCTCQ+VAMTF+TCTVB+VAMTF+TCLIA+VAMTF+1st part of VCTCP
327 Float_t totLength1 = 160.8 + 78. + 148. + 78. + 148. + 78. + 9.3;
328 //
329 tubpar[0] = 18.6/2.;
330 tubpar[1] = 7.6/2.;
331 tubpar[2] = totLength1/2.;
332// gMC->Gsvolu("QE01", "ELTU", idtmed[7], tubpar, 3);
333 // temporary replace with a scaled tube (AG)
334 TGeoTube *tubeQE01 = new TGeoTube(0.,tubpar[0],tubpar[2]);
335 TGeoScale *scaleQE01 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
336 TGeoScaledShape *sshapeQE01 = new TGeoScaledShape(tubeQE01, scaleQE01);
337 new TGeoVolume("QE01", sshapeQE01, gGeoManager->GetMedium(idtmed[7]));
338
339 tubpar[0] = 18.0/2.;
340 tubpar[1] = 7.0/2.;
341 tubpar[2] = totLength1/2.;
342// gMC->Gsvolu("QE02", "ELTU", idtmed[10], tubpar, 3);
343 // temporary replace with a scaled tube (AG)
344 TGeoTube *tubeQE02 = new TGeoTube(0.,tubpar[0],tubpar[2]);
345 TGeoScale *scaleQE02 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
346 TGeoScaledShape *sshapeQE02 = new TGeoScaledShape(tubeQE02, scaleQE02);
347 new TGeoVolume("QE02", sshapeQE02, gGeoManager->GetMedium(idtmed[10]));
348
349 gMC->Gspos("QE01", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
350 gMC->Gspos("QE02", 1, "QE01", 0., 0., 0., 0, "ONLY");
351 // Ch.debug
352 //printf(" QE01 ELTU from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
353
354 // Vertical collimator jaws (defined ONLY if fVCollAperture<3.5!)
355 if(fVCollSideCAperture<3.5){
356 boxpar[0] = 5.4/2.;
357 boxpar[1] = (3.5-fVCollSideCAperture-fVCollSideCCentreY-0.7)/2.;
358 if(boxpar[1]<0.) boxpar[1]=0.;
359 boxpar[2] = 124.4/2.;
360 printf(" AliZDCv3 -> sideC VCollimator jaws: apertures +%1.2f/-%1.2f center %1.2f [cm]\n",
361 fVCollSideCAperture, fVCollSideCApertureNeg,fVCollSideCCentreY);
362 gMC->Gsvolu("QCVC" , "BOX ", idtmed[13], boxpar, 3);
363 gMC->Gspos("QCVC", 1, "QE02", -boxpar[0], fVCollSideCAperture+fVCollSideCCentreY+boxpar[1], -totLength1/2.+160.8+78.+148./2., 0, "ONLY");
364 gMC->Gspos("QCVC", 2, "QE02", -boxpar[0], -fVCollSideCApertureNeg+fVCollSideCCentreY-boxpar[1], -totLength1/2.+160.8+78.+148./2., 0, "ONLY");
365 }
366
367 zd1 += tubpar[2] * 2.;
368
369 // 2nd part of VCTCP
370 conpar[0] = 31.5/2.;
371 conpar[1] = 21.27/2.;
372 conpar[2] = 21.87/2.;
373 conpar[3] = 18.0/2.;
374 conpar[4] = 18.6/2.;
375 gMC->Gsvolu("QC02", "CONE", idtmed[7], conpar, 5);
376 gMC->Gspos("QC02", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
377 // Ch.debug
378 //printf(" QC02 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
379
380 zd1 += conpar[0] * 2.;
381
382 // 3rd section of VCTCP+VCDWC+VMLGB
383 Float_t totLenght2 = 9.2 + 530.5+40.;
384 tubpar[0] = 21.2/2.;
385 tubpar[1] = 21.9/2.;
386 tubpar[2] = totLenght2/2.;
387 gMC->Gsvolu("QT04", "TUBE", idtmed[7], tubpar, 3);
388 gMC->Gspos("QT04", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
389 // Ch.debug
390 //printf(" QT04 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
391
392 zd1 += tubpar[2] * 2.;
393
394 // First part of VCTCD
395 // skewed transition cone from ID=212.7 mm to ID=797 mm
396 conpar[0] = 121./2.;
397 conpar[1] = 79.7/2.;
398 conpar[2] = 81.3/2.;
399 conpar[3] = 21.27/2.;
400 conpar[4] = 21.87/2.;
401 gMC->Gsvolu("QC03", "CONE", idtmed[7], conpar, 5);
402 gMC->Gspos("QC03", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
403 // Ch.debug
404 //printf(" QC03 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
405
406 zd1 += 2.*conpar[0];
407
408 // VCDGB + 1st part of VCTCH
409 tubpar[0] = 79.7/2.;
410 tubpar[1] = 81.3/2.;
411 tubpar[2] = (5*475.2+97.)/2.;
412 gMC->Gsvolu("QT05", "TUBE", idtmed[7], tubpar, 3);
413 gMC->Gspos("QT05", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
414 // Ch.debug
415 //printf(" QT05 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
416
417 zd1 += 2.*tubpar[2];
418
419 // 2nd part of VCTCH
420 // Transition from ID=797 mm to ID=196 mm:
421 // in order to simulate the thin window opened in the transition cone
422 // we divide the transition cone in three cones:
423 // (1) 8 mm thick (2) 3 mm thick (3) the third 8 mm thick
424
425 // (1) 8 mm thick
426 conpar[0] = 9.09/2.; // 15 degree
427 conpar[1] = 74.82868/2.;
428 conpar[2] = 76.42868/2.; // thickness 8 mm
429 conpar[3] = 79.7/2.;
430 conpar[4] = 81.3/2.; // thickness 8 mm
431 gMC->Gsvolu("QC04", "CONE", idtmed[7], conpar, 5);
432 gMC->Gspos("QC04", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
433 // Ch.debug
434 //printf(" QC04 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
435
436 zd1 += 2.*conpar[0];
437
438 // (2) 3 mm thick
439 conpar[0] = 96.2/2.; // 15 degree
440 conpar[1] = 23.19588/2.;
441 conpar[2] = 23.79588/2.; // thickness 3 mm
442 conpar[3] = 74.82868/2.;
443 conpar[4] = 75.42868/2.; // thickness 3 mm
444 gMC->Gsvolu("QC05", "CONE", idtmed[7], conpar, 5);
445 gMC->Gspos("QC05", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
446 // Ch.debug
447 //printf(" QC05 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
448
449 zd1 += 2.*conpar[0];
450
451 // (3) 8 mm thick
452 conpar[0] = 6.71/2.; // 15 degree
453 conpar[1] = 19.6/2.;
454 conpar[2] = 21.2/2.;// thickness 8 mm
455 conpar[3] = 23.19588/2.;
456 conpar[4] = 24.79588/2.;// thickness 8 mm
457 gMC->Gsvolu("QC06", "CONE", idtmed[7], conpar, 5);
458 gMC->Gspos("QC06", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
459 // Ch.debug
460 //printf(" QC06 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
461
462 zd1 += 2.*conpar[0];
463
464 // VMZAR (5 volumes)
465 tubpar[0] = 20.2/2.;
466 tubpar[1] = 20.6/2.;
467 tubpar[2] = 2.15/2.;
468 gMC->Gsvolu("QT06", "TUBE", idtmed[7], tubpar, 3);
469 gMC->Gspos("QT06", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
470 // Ch.debug
471 //printf(" QT06 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
472
473 zd1 += 2.*tubpar[2];
474
475 conpar[0] = 6.9/2.;
476 conpar[1] = 23.9/2.;
477 conpar[2] = 24.3/2.;
478 conpar[3] = 20.2/2.;
479 conpar[4] = 20.6/2.;
480 gMC->Gsvolu("QC07", "CONE", idtmed[7], conpar, 5);
481 gMC->Gspos("QC07", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
482 // Ch.debug
483 //printf(" QC07 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
484
485 zd1 += 2.*conpar[0];
486
487 tubpar[0] = 23.9/2.;
488 tubpar[1] = 25.5/2.;
489 tubpar[2] = 17.0/2.;
490 gMC->Gsvolu("QT07", "TUBE", idtmed[7], tubpar, 3);
491 gMC->Gspos("QT07", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
492 // Ch.debug
493 //printf(" QT07 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
494
495 zd1 += 2.*tubpar[2];
496
497 conpar[0] = 6.9/2.;
498 conpar[1] = 20.2/2.;
499 conpar[2] = 20.6/2.;
500 conpar[3] = 23.9/2.;
501 conpar[4] = 24.3/2.;
502 gMC->Gsvolu("QC08", "CONE", idtmed[7], conpar, 5);
503 gMC->Gspos("QC08", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
504 // Ch.debug
505 //printf(" QC08 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
506
507 zd1 += 2.*conpar[0];
508
509 tubpar[0] = 20.2/2.;
510 tubpar[1] = 20.6/2.;
511 tubpar[2] = 2.15/2.;
512 gMC->Gsvolu("QT08", "TUBE", idtmed[7], tubpar, 3);
513 gMC->Gspos("QT08", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
514 // Ch.debug
515 //printf(" QT08 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
516
517 zd1 += 2.*tubpar[2];
518
519 // Flange (ID=196 mm)(last part of VMZAR and first part of VCTYB)
520 tubpar[0] = 19.6/2.;
521 tubpar[1] = 25.3/2.;
522 tubpar[2] = 4.9/2.;
523 gMC->Gsvolu("QT09", "TUBE", idtmed[7], tubpar, 3);
524 gMC->Gspos("QT09", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
525 // Ch.debug
526 //printf(" QT09 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
527
528 zd1 += 2.*tubpar[2];
529 // Ch.debug
530 //printf(" Beginning of VCTYB volume @ z = %1.2f \n",-zd1);
531
532 // simulation of the trousers (VCTYB)
533 tubpar[0] = 19.6/2.;
534 tubpar[1] = 20.0/2.;
535 tubpar[2] = 3.9/2.;
536 gMC->Gsvolu("QT10", "TUBE", idtmed[7], tubpar, 3);
537 gMC->Gspos("QT10", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
538 // Ch.debug
539 //printf(" QT10 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
540
541 zd1 += 2.*tubpar[2];
542
543 // transition cone from ID=196. to ID=216.6
544 conpar[0] = 32.55/2.;
545 conpar[1] = 21.66/2.;
546 conpar[2] = 22.06/2.;
547 conpar[3] = 19.6/2.;
548 conpar[4] = 20.0/2.;
549 gMC->Gsvolu("QC09", "CONE", idtmed[7], conpar, 5);
550 gMC->Gspos("QC09", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
551 // Ch.debug
552 //printf(" QC09 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
553
554 zd1 += 2.*conpar[0];
555
556 // tube
557 tubpar[0] = 21.66/2.;
558 tubpar[1] = 22.06/2.;
559 tubpar[2] = 28.6/2.;
560 gMC->Gsvolu("QT11", "TUBE", idtmed[7], tubpar, 3);
561 gMC->Gspos("QT11", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
562 // Ch.debug
563 //printf(" QT11 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
564
565 zd1 += 2.*tubpar[2];
566 // Ch.debug
567 //printf(" Beginning of recombination chamber @ z = %f \n",-zd1);
568
569 // --------------------------------------------------------
570 // RECOMBINATION CHAMBER IMPLEMENTED USING TGeo CLASSES!!!!
571 // author: Chiara (August 2008)
572 // --------------------------------------------------------
573 // TRANSFORMATION MATRICES
574 // Combi transformation:
575 Double_t dx = -3.970000;
576 Double_t dy = 0.000000;
577 Double_t dz = 0.0;
578 // Rotation:
579 Double_t thx = 84.989100; Double_t phx = 180.000000;
580 Double_t thy = 90.000000; Double_t phy = 90.000000;
581 Double_t thz = 185.010900; Double_t phz = 0.000000;
582 TGeoRotation *rotMatrix1c = new TGeoRotation("c",thx,phx,thy,phy,thz,phz);
583 // Combi transformation:
584 dx = -3.970000;
585 dy = 0.000000;
586 dz = 0.0;
587 TGeoCombiTrans *rotMatrix2c = new TGeoCombiTrans("ZDCC_c1", dx,dy,dz,rotMatrix1c);
588 rotMatrix2c->RegisterYourself();
589 // Combi transformation:
590 dx = 3.970000;
591 dy = 0.000000;
592 dz = 0.0;
593 // Rotation:
594 thx = 95.010900; phx = 180.000000;
595 thy = 90.000000; phy = 90.000000;
596 thz = 180.-5.010900; phz = 0.000000;
597 TGeoRotation *rotMatrix3c = new TGeoRotation("",thx,phx,thy,phy,thz,phz);
598 TGeoCombiTrans *rotMatrix4c = new TGeoCombiTrans("ZDCC_c2", dx,dy,dz,rotMatrix3c);
599 rotMatrix4c->RegisterYourself();
600
601 // VOLUMES DEFINITION
602 // Volume: ZDCC
603 TGeoVolume *pZDCC = gGeoManager->GetVolume("ZDCC");
604
605 conpar[0] = (90.1-0.95-0.26-0.0085)/2.;
606 conpar[1] = 0.0/2.;
607 conpar[2] = 21.6/2.;
608 conpar[3] = 0.0/2.;
609 conpar[4] = 5.8/2.;
610 new TGeoCone("QCLext", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
611
612 conpar[0] = (90.1-0.95-0.26-0.0085)/2.;
613 conpar[1] = 0.0/2.;
614 conpar[2] = 21.2/2.;
615 conpar[3] = 0.0/2.;
616 conpar[4] = 5.4/2.;
617 new TGeoCone("QCLint", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
618
619 // Outer trousers
620 TGeoCompositeShape *pOutTrousersC = new TGeoCompositeShape("outTrousersC", "QCLext:ZDCC_c1+QCLext:ZDCC_c2");
621
622 // Volume: QCLext
623 TGeoMedium *medZDCFe = gGeoManager->GetMedium("ZDC_ZIRONT");
624 TGeoVolume *pQCLext = new TGeoVolume("QCLext",pOutTrousersC, medZDCFe);
625 pQCLext->SetLineColor(kGreen);
626 pQCLext->SetVisLeaves(kTRUE);
627 //
628 TGeoTranslation *tr1c = new TGeoTranslation(0., 0., (Double_t) -conpar[0]-0.95-zd1);
629 //printf(" Recombination chamber from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-0.95-zd1);
630 //
631 pZDCC->AddNode(pQCLext, 1, tr1c);
632 // Inner trousers
633 TGeoCompositeShape *pIntTrousersC = new TGeoCompositeShape("intTrousersC", "QCLint:ZDCC_c1+QCLint:ZDCC_c2");
634 // Volume: QCLint
635 TGeoMedium *medZDCvoid = gGeoManager->GetMedium("ZDC_ZVOID");
636 TGeoVolume *pQCLint = new TGeoVolume("QCLint",pIntTrousersC, medZDCvoid);
637 pQCLint->SetLineColor(kTeal);
638 pQCLint->SetVisLeaves(kTRUE);
639 pQCLext->AddNode(pQCLint, 1);
640
641 zd1 += 90.1;
642 Double_t offset = 0.5;
643 zd1 = zd1+offset;
644
645 // second section : 2 tubes (ID = 54. OD = 58.)
646 tubpar[0] = 5.4/2.;
647 tubpar[1] = 5.8/2.;
648 tubpar[2] = 40.0/2.;
649 gMC->Gsvolu("QT12", "TUBE", idtmed[7], tubpar, 3);
650 gMC->Gspos("QT12", 1, "ZDCC", -15.8/2., 0., -tubpar[2]-zd1, 0, "ONLY");
651 gMC->Gspos("QT12", 2, "ZDCC", 15.8/2., 0., -tubpar[2]-zd1, 0, "ONLY");
652 // Ch.debug
653 //printf(" QT12 TUBE from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
654
655 zd1 += 2.*tubpar[2];
656
657 // transition x2zdc to recombination chamber : skewed cone
658 conpar[0] = (10.-0.2-offset)/2.;
659 conpar[1] = 6.3/2.;
660 conpar[2] = 7.0/2.;
661 conpar[3] = 5.4/2.;
662 conpar[4] = 5.8/2.;
663 gMC->Gsvolu("QC10", "CONE", idtmed[7], conpar, 5);
664 gMC->Gspos("QC10", 1, "ZDCC", -7.9-0.175, 0., -conpar[0]-0.1-zd1, irotpipe1, "ONLY");
665 gMC->Gspos("QC10", 2, "ZDCC", 7.9+0.175, 0., -conpar[0]-0.1-zd1, irotpipe2, "ONLY");
666 //printf(" QC10 CONE from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-0.2-zd1);
667
668 zd1 += 2.*conpar[0]+0.2;
669
670 // 2 tubes (ID = 63 mm OD=70 mm)
671 tubpar[0] = 6.3/2.;
672 tubpar[1] = 7.0/2.;
673 tubpar[2] = 639.8/2.;
674 gMC->Gsvolu("QT13", "TUBE", idtmed[7], tubpar, 3);
675 gMC->Gspos("QT13", 1, "ZDCC", -16.5/2., 0., -tubpar[2]-zd1, 0, "ONLY");
676 gMC->Gspos("QT13", 2, "ZDCC", 16.5/2., 0., -tubpar[2]-zd1, 0, "ONLY");
677 //printf(" QT13 TUBE from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
678
679 zd1 += 2.*tubpar[2];
680 //printf(" END OF SIDE C BEAM PIPE DEFINITION @ z = %f\n",-zd1);
681
682
683 // -- Luminometer (Cu box) in front of ZN - side C
684 boxpar[0] = 8.0/2.;
685 boxpar[1] = 8.0/2.;
686 boxpar[2] = fLumiLength/2.;
687 gMC->Gsvolu("QLUC", "BOX ", idtmed[9], boxpar, 3);
688 gMC->Gspos("QLUC", 1, "ZDCC", 0., 0., fPosZNC[2]+66.+boxpar[2], 0, "ONLY");
689 //printf(" QLUC LUMINOMETER from z = %1.2f to z= %1.2f\n", fPosZNC[2]+66., fPosZNC[2]+66.+2*boxpar[2]);
690
691 // -- END OF BEAM PIPE VOLUME DEFINITION FOR SIDE C (RB26 SIDE)
692 // ----------------------------------------------------------------
693
694 ////////////////////////////////////////////////////////////////
695 // //
696 // SIDE A - RB24 //
697 // //
698 ///////////////////////////////////////////////////////////////
699
700 // Rotation Matrices definition
701 Int_t irotpipe3, irotpipe4, irotpipe5;
702 //-- rotation matrices for the tilted cone after the TDI to recenter vacuum chamber
703 gMC->Matrix(irotpipe3,90.-1.8934,0.,90.,90.,1.8934,180.);
704 //-- rotation matrices for the tilted tube before and after the TDI
705 gMC->Matrix(irotpipe4,90.-3.8,0.,90.,90.,3.8,180.);
706 //-- rotation matrix for the tilted cone after the TDI
707 gMC->Matrix(irotpipe5,90.+9.8,0.,90.,90.,9.8,0.);
708
709 // -- Mother of the ZDCs (Vacuum PCON)
710 zd2 = 1910.22;// zd2 initial value
711
712 conpar[0] = 0.;
713 conpar[1] = 360.;
714 conpar[2] = 2.;
715 conpar[3] = zd2;
716 conpar[4] = 0.;
717 conpar[5] = 55.;
718 conpar[6] = 13500.;
719 conpar[7] = 0.;
720 conpar[8] = 55.;
721 gMC->Gsvolu("ZDCA", "PCON", idtmed[10], conpar, 9);
722 gMC->Gspos("ZDCA", 1, "ALIC", 0., 0., 0., 0, "ONLY");
723
724 // To avoid overlaps 1 micron are left between certain volumes!
725 Double_t dxNoOverlap = 0.0;
726 //zd2 += dxNoOverlap;
727
728 // BEAM PIPE from 19.10 m to inner triplet beginning (22.965 m)
729 tubpar[0] = 6.0/2.;
730 tubpar[1] = 6.4/2.;
731 tubpar[2] = 386.28/2. - dxNoOverlap;
732 gMC->Gsvolu("QA01", "TUBE", idtmed[7], tubpar, 3);
733 gMC->Gspos("QA01", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
734 // Ch.debug
735 //printf(" QA01 TUBE centred in %f from z = %1.2f to z= %1.2f (IT begin)\n",tubpar[2]+zd2,zd2,2*tubpar[2]+zd2);
736
737 zd2 += 2.*tubpar[2];
738
739 // -- FIRST SECTION OF THE BEAM PIPE (from beginning of inner triplet to
740 // beginning of D1)
741 tubpar[0] = 6.3/2.;
742 tubpar[1] = 6.7/2.;
743 tubpar[2] = 3541.8/2. - dxNoOverlap;
744 gMC->Gsvolu("QA02", "TUBE", idtmed[7], tubpar, 3);
745 gMC->Gspos("QA02", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
746 // Ch.debug
747 //printf(" QA02 TUBE from z = %1.2f to z= %1.2f (D1 begin)\n",zd2,2*tubpar[2]+zd2);
748
749 zd2 += 2.*tubpar[2];
750
751
752 // -- SECOND SECTION OF THE BEAM PIPE (from the beginning of D1 to the beginning of D2)
753 //
754 // FROM (MAGNETIC) BEGINNING OF D1 TO THE (MAGNETIC) END OF D1 + 126.5 cm
755 // CYLINDRICAL PIPE of diameter increasing from 6.75 cm up to 8.0 cm
756 // from magnetic end :
757 // 1) 80.1 cm still with ID = 6.75 radial beam screen
758 // 2) 2.5 cm conical section from ID = 6.75 to ID = 8.0 cm
759 // 3) 43.9 cm straight section (tube) with ID = 8.0 cm
760
761 tubpar[0] = 6.75/2.;
762 tubpar[1] = 7.15/2.;
763 tubpar[2] = (945.0+80.1)/2.;
764 gMC->Gsvolu("QA03", "TUBE", idtmed[7], tubpar, 3);
765 gMC->Gspos("QA03", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
766 // Ch.debug
767 //printf(" QA03 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
768
769 zd2 += 2.*tubpar[2];
770
771 // Transition Cone from ID=67.5 mm to ID=80 mm
772 conpar[0] = 2.5/2.;
773 conpar[1] = 6.75/2.;
774 conpar[2] = 7.15/2.;
775 conpar[3] = 8.0/2.;
776 conpar[4] = 8.4/2.;
777 gMC->Gsvolu("QA04", "CONE", idtmed[7], conpar, 5);
778 gMC->Gspos("QA04", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
779 //printf(" QA04 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
780
781 zd2 += 2.*conpar[0];
782
783 tubpar[0] = 8.0/2.;
784 tubpar[1] = 8.4/2.;
785 tubpar[2] = (43.9+20.+28.5+28.5)/2.;
786 gMC->Gsvolu("QA05", "TUBE", idtmed[7], tubpar, 3);
787 gMC->Gspos("QA05", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
788 // Ch.debug
789 //printf(" QA05 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
790
791 zd2 += 2.*tubpar[2];
792
793 // Second section of VAEHI (transition cone from ID=80mm to ID=98mm)
794 conpar[0] = 4.0/2.;
795 conpar[1] = 8.0/2.;
796 conpar[2] = 8.4/2.;
797 conpar[3] = 9.8/2.;
798 conpar[4] = 10.2/2.;
799 gMC->Gsvolu("QAV1", "CONE", idtmed[7], conpar, 5);
800 gMC->Gspos("QAV1", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
801 //printf(" QAV1 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
802
803 zd2 += 2.*conpar[0];
804
805 //Third section of VAEHI (transition cone from ID=98mm to ID=90mm)
806 conpar[0] = 1.0/2.;
807 conpar[1] = 9.8/2.;
808 conpar[2] = 10.2/2.;
809 conpar[3] = 9.0/2.;
810 conpar[4] = 9.4/2.;
811 gMC->Gsvolu("QAV2", "CONE", idtmed[7], conpar, 5);
812 gMC->Gspos("QAV2", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
813 //printf(" QAV2 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
814
815 zd2 += 2.*conpar[0];
816
817 // Fourth section of VAEHI (tube ID=90mm)
818 tubpar[0] = 9.0/2.;
819 tubpar[1] = 9.4/2.;
820 tubpar[2] = 31.0/2.;
821 gMC->Gsvolu("QAV3", "TUBE", idtmed[7], tubpar, 3);
822 gMC->Gspos("QAV3", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
823 // Ch.debug
824 //printf(" QAV3 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
825
826 zd2 += 2.*tubpar[2];
827
828 //---------------------------- TCDD beginning ----------------------------------
829 // space for the insertion of the collimator TCDD (2 m)
830 // TCDD ZONE - 1st volume
831 conpar[0] = 1.3/2.;
832 conpar[1] = 9.0/2.;
833 conpar[2] = 13.0/2.;
834 conpar[3] = 9.6/2.;
835 conpar[4] = 13.0/2.;
836 gMC->Gsvolu("Q01T", "CONE", idtmed[7], conpar, 5);
837 gMC->Gspos("Q01T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
838 //printf(" Q01T CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
839
840 zd2 += 2.*conpar[0];
841
842 // TCDD ZONE - 2nd volume
843 tubpar[0] = 9.6/2.;
844 tubpar[1] = 10.0/2.;
845 tubpar[2] = 1.0/2.;
846 gMC->Gsvolu("Q02T", "TUBE", idtmed[7], tubpar, 3);
847 gMC->Gspos("Q02T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
848 // Ch.debug
849 //printf(" Q02T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
850
851 zd2 += 2.*tubpar[2];
852
853 // TCDD ZONE - third volume
854 conpar[0] = 9.04/2.;
855 conpar[1] = 9.6/2.;
856 conpar[2] = 10.0/2.;
857 conpar[3] = 13.8/2.;
858 conpar[4] = 14.2/2.;
859 gMC->Gsvolu("Q03T", "CONE", idtmed[7], conpar, 5);
860 gMC->Gspos("Q03T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
861 //printf(" Q03T CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
862
863 zd2 += 2.*conpar[0];
864
865 // TCDD ZONE - 4th volume
866 tubpar[0] = 13.8/2.;
867 tubpar[1] = 14.2/2.;
868 tubpar[2] = 38.6/2.;
869 gMC->Gsvolu("Q04T", "TUBE", idtmed[7], tubpar, 3);
870 gMC->Gspos("Q04T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
871 // Ch.debug
872 //printf(" Q04T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
873
874 zd2 += 2.*tubpar[2];
875
876 // TCDD ZONE - 5th volume
877 tubpar[0] = 21.0/2.;
878 tubpar[1] = 21.4/2.;
879 tubpar[2] = 100.12/2.;
880 gMC->Gsvolu("Q05T", "TUBE", idtmed[7], tubpar, 3);
881 gMC->Gspos("Q05T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
882 // Ch.debug
883 //printf(" Q05T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
884
885 zd2 += 2.*tubpar[2];
886
887 // TCDD ZONE - 6th volume
888 tubpar[0] = 13.8/2.;
889 tubpar[1] = 14.2/2.;
890 tubpar[2] = 38.6/2.;
891 gMC->Gsvolu("Q06T", "TUBE", idtmed[7], tubpar, 3);
892 gMC->Gspos("Q06T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
893 // Ch.debug
894 //printf(" Q06T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
895
896 zd2 += 2.*tubpar[2];
897
898 // TCDD ZONE - 7th volume
899 conpar[0] = 11.34/2.;
900 conpar[1] = 13.8/2.;
901 conpar[2] = 14.2/2.;
902 conpar[3] = 18.0/2.;
903 conpar[4] = 18.4/2.;
904 gMC->Gsvolu("Q07T", "CONE", idtmed[7], conpar, 5);
905 gMC->Gspos("Q07T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
906 //printf(" Q07T CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
907
908 zd2 += 2.*conpar[0];
909
910 // Upper section : one single phi segment of a tube
911 // 5 parameters for tubs: inner radius = 0.,
912 // outer radius = 7. cm, half length = 50 cm
913 // phi1 = 0., phi2 = 180.
914 tubspar[0] = 0.0/2.;
915 tubspar[1] = 14.0/2.;
916 tubspar[2] = 100.0/2.;
917 tubspar[3] = 0.;
918 tubspar[4] = 180.;
919 gMC->Gsvolu("Q08T", "TUBS", idtmed[7], tubspar, 5);
920 // Ch.debug
921 //printf(" upper part : one single phi segment of a tube (Q08T)\n");
922
923 // rectangular beam pipe inside TCDD upper section (Vacuum)
924 boxpar[0] = 7.0/2.;
925 boxpar[1] = 2.2/2.;
926 boxpar[2] = 100./2.;
927 gMC->Gsvolu("Q09T", "BOX ", idtmed[10], boxpar, 3);
928 // positioning vacuum box in the upper section of TCDD
929 gMC->Gspos("Q09T", 1, "Q08T", 0., 1.1, 0., 0, "ONLY");
930
931 // lower section : one single phi segment of a tube
932 tubspar[0] = 0.0/2.;
933 tubspar[1] = 14.0/2.;
934 tubspar[2] = 100.0/2.;
935 tubspar[3] = 180.;
936 tubspar[4] = 360.;
937 gMC->Gsvolu("Q10T", "TUBS", idtmed[7], tubspar, 5);
938 // rectangular beam pipe inside TCDD lower section (Vacuum)
939 boxpar[0] = 7.0/2.;
940 boxpar[1] = 2.2/2.;
941 boxpar[2] = 100./2.;
942 gMC->Gsvolu("Q11T", "BOX ", idtmed[10], boxpar, 3);
943 // positioning vacuum box in the lower section of TCDD
944 gMC->Gspos("Q11T", 1, "Q10T", 0., -1.1, 0., 0, "ONLY");
945
946 // positioning TCDD elements in ZDCA, (inside TCDD volume)
947 gMC->Gspos("Q08T", 1, "ZDCA", 0., fTCDDAperturePos, -100.+zd2, 0, "ONLY");
948 gMC->Gspos("Q10T", 1, "ZDCA", 0., -fTCDDApertureNeg, -100.+zd2, 0, "ONLY");
949 printf(" AliZDCv3 -> TCDD apertures +%1.2f/-%1.2f cm\n",
950 fTCDDAperturePos, fTCDDApertureNeg);
951
952 // RF screen
953 boxpar[0] = 0.2/2.;
954 boxpar[1] = 4.0/2.;
955 boxpar[2] = 100./2.;
956 gMC->Gsvolu("Q12T", "BOX ", idtmed[7], boxpar, 3);
957 // positioning RF screen at both sides of TCDD
958 gMC->Gspos("Q12T", 1, "ZDCA", tubspar[1]+boxpar[0], 0., -100.+zd2, 0, "ONLY");
959 gMC->Gspos("Q12T", 2, "ZDCA", -tubspar[1]-boxpar[0], 0., -100.+zd2, 0, "ONLY");
960 //---------------------------- TCDD end ---------------------------------------
961
962 // The following elliptical tube 180 mm x 70 mm
963 // (obtained positioning the void QA09 in QA08)
964 // represents VMTSA (780 mm) + space reserved to the TCTVB (1480 mm)+
965 // VMTSA (780 mm) + first part of VCTCP (93 mm)
966
967 tubpar[0] = 18.4/2.;
968 tubpar[1] = 7.4/2.;
969 tubpar[2] = 313.3/2.;
970// gMC->Gsvolu("QA06", "ELTU", idtmed[7], tubpar, 3);
971 // temporary replace with a scaled tube (AG)
972 TGeoTube *tubeQA06 = new TGeoTube(0.,tubpar[0],tubpar[2]);
973 TGeoScale *scaleQA06 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
974 TGeoScaledShape *sshapeQA06 = new TGeoScaledShape(tubeQA06, scaleQA06);
975 new TGeoVolume("QA06", sshapeQA06, gGeoManager->GetMedium(idtmed[7]));
976 //printf(" QA06 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
977
978 tubpar[0] = 18.0/2.;
979 tubpar[1] = 7.0/2.;
980 tubpar[2] = 313.3/2.;
981// gMC->Gsvolu("QA07", "ELTU", idtmed[10], tubpar, 3);
982 // temporary replace with a scaled tube (AG)
983 TGeoTube *tubeQA07 = new TGeoTube(0.,tubpar[0],tubpar[2]);
984 TGeoScale *scaleQA07 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
985 TGeoScaledShape *sshapeQA07 = new TGeoScaledShape(tubeQA07, scaleQA07);
986 new TGeoVolume("QA07", sshapeQA07, gGeoManager->GetMedium(idtmed[10]));
987 //printf(" QA07 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
988 gMC->Gspos("QA06", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
989 gMC->Gspos("QA07", 1, "QA06", 0., 0., 0., 0, "ONLY");
990
991 // Vertical collimator jaws (defined ONLY if fVCollAperture<3.5!)
992 if(fVCollSideAAperture<3.5){
993 boxpar[0] = 5.4/2.;
994 boxpar[1] = (3.5-fVCollSideAAperture-fVCollSideACentreY-0.7)/2.;
995 if(boxpar[1]<0.) boxpar[1]=0.;
996 boxpar[2] = 124.4/2.;
997 gMC->Gsvolu("QCVA" , "BOX ", idtmed[13], boxpar, 3);
998 gMC->Gspos("QCVA", 1, "QA07", -boxpar[0], fVCollSideAAperture+fVCollSideACentreY+boxpar[1], -313.3/2.+78.+148./2., 0, "ONLY");
999 gMC->Gspos("QCVA", 2, "QA07", -boxpar[0], -fVCollSideAApertureNeg+fVCollSideACentreY-boxpar[1], -313.3/2.+78.+148./2., 0, "ONLY");
1000 printf(" AliZDCv3 -> sideA VCollimator jaws: apertures +%1.2f/-%1.2f center %1.2f [cm]\n",
1001 fVCollSideAAperture, fVCollSideAApertureNeg,fVCollSideACentreY);
1002 }
1003
1004 zd2 += 2.*tubpar[2];
1005
1006 // VCTCP second part: transition cone from ID=180 to ID=212.7
1007 conpar[0] = 31.5/2.;
1008 conpar[1] = 18.0/2.;
1009 conpar[2] = 18.6/2.;
1010 conpar[3] = 21.27/2.;
1011 conpar[4] = 21.87/2.;
1012 gMC->Gsvolu("QA08", "CONE", idtmed[7], conpar, 5);
1013 gMC->Gspos("QA08", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1014 // Ch.debug
1015 //printf(" QA08 CONE from z = %Third part of VCTCR: tube (ID=196 mm) f to z = %f\n",zd2,2*conpar[0]+zd2);
1016
1017 zd2 += 2.*conpar[0];
1018
1019 // Tube ID 212.7 mm
1020 // Represents VCTCP third part (92 mm) + VCDWB (765 mm) + VMBGA (400 mm) +
1021 // VCDWE (300 mm) + VMBGA (400 mm)
1022 tubpar[0] = 21.27/2.;
1023 tubpar[1] = 21.87/2.;
1024 tubpar[2] = 195.7/2.;
1025 gMC->Gsvolu("QA09", "TUBE", idtmed[7], tubpar, 3);
1026 gMC->Gspos("QA09", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1027 //printf(" QA09 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1028
1029 zd2 += 2.*tubpar[2];
1030
1031 // skewed transition piece (ID=212.7 mm to 332 mm) (before TDI)
1032 conpar[0] = (50.0-0.73-1.13)/2.;
1033 conpar[1] = 21.27/2.;
1034 conpar[2] = 21.87/2.;
1035 conpar[3] = 33.2/2.;
1036 conpar[4] = 33.8/2.;
1037 gMC->Gsvolu("QA10", "CONE", idtmed[7], conpar, 5);
1038 gMC->Gspos("QA10", 1, "ZDCA", -1.66, 0., conpar[0]+0.73+zd2, irotpipe4, "ONLY");
1039 // Ch.debug
1040 //printf(" QA10 skewed CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+0.73+1.13+zd2);
1041
1042 zd2 += 2.*conpar[0]+0.73+1.13;
1043
1044 // Vacuum chamber containing TDI
1045 tubpar[0] = 0.;
1046 tubpar[1] = 54.6/2.;
1047 tubpar[2] = 540.0/2.;
1048 gMC->Gsvolu("Q13TM", "TUBE", idtmed[10], tubpar, 3);
1049 gMC->Gspos("Q13TM", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1050 tubpar[0] = 54.0/2.;
1051 tubpar[1] = 54.6/2.;
1052 tubpar[2] = 540.0/2.;
1053 gMC->Gsvolu("Q13T", "TUBE", idtmed[7], tubpar, 3);
1054 gMC->Gspos("Q13T", 1, "Q13TM", 0., 0., 0., 0, "ONLY");
1055 // Ch.debug
1056 //printf(" Q13T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1057
1058 zd2 += 2.*tubpar[2];
1059
1060 //---------------- INSERT TDI INSIDE Q13T -----------------------------------
1061 boxpar[0] = 11.0/2.;
1062 boxpar[1] = 9.0/2.;
1063 boxpar[2] = 540.0/2.;
1064 gMC->Gsvolu("QTD1", "BOX ", idtmed[7], boxpar, 3);
1065 gMC->Gspos("QTD1", 1, "Q13TM", -3.8, boxpar[1]+fTDIAperturePos, 0., 0, "ONLY");
1066 boxpar[0] = 11.0/2.;
1067 boxpar[1] = 9.0/2.;
1068 boxpar[2] = 540.0/2.;
1069 gMC->Gsvolu("QTD2", "BOX ", idtmed[7], boxpar, 3);
1070 gMC->Gspos("QTD2", 1, "Q13TM", -3.8, -boxpar[1]-fTDIApertureNeg, 0., 0, "ONLY");
1071 boxpar[0] = 5.1/2.;
1072 boxpar[1] = 0.2/2.;
1073 boxpar[2] = 540.0/2.;
1074 gMC->Gsvolu("QTD3", "BOX ", idtmed[7], boxpar, 3);
1075 gMC->Gspos("QTD3", 1, "Q13TM", -3.8+5.5+boxpar[0], fTDIAperturePos, 0., 0, "ONLY");
1076 gMC->Gspos("QTD3", 2, "Q13TM", -3.8+5.5+boxpar[0], -fTDIApertureNeg, 0., 0, "ONLY");
1077 gMC->Gspos("QTD3", 3, "Q13TM", -3.8-5.5-boxpar[0], fTDIAperturePos, 0., 0, "ONLY");
1078 gMC->Gspos("QTD3", 4, "Q13TM", -3.8-5.5-boxpar[0], -fTDIApertureNeg, 0., 0, "ONLY");
1079 printf(" AliZDCv3 -> TDI apertures +%1.2f/-%1.2f cm\n\n",
1080 fTDIAperturePos, fTDIApertureNeg);
1081 //
1082 tubspar[0] = 12.0/2.;
1083 tubspar[1] = 12.4/2.;
1084 tubspar[2] = 540.0/2.;
1085 tubspar[3] = 90.;
1086 tubspar[4] = 270.;
1087 gMC->Gsvolu("QTD4", "TUBS", idtmed[7], tubspar, 5);
1088 gMC->Gspos("QTD4", 1, "Q13TM", -3.8-10.6, 0., 0., 0, "ONLY");
1089 tubspar[0] = 12.0/2.;
1090 tubspar[1] = 12.4/2.;
1091 tubspar[2] = 540.0/2.;
1092 tubspar[3] = -90.;
1093 tubspar[4] = 90.;
1094 gMC->Gsvolu("QTD5", "TUBS", idtmed[7], tubspar, 5);
1095 gMC->Gspos("QTD5", 1, "Q13TM", -3.8+10.6, 0., 0., 0, "ONLY");
1096 //---------------- END DEFINING TDI INSIDE Q13T -------------------------------
1097
1098 // VCTCG skewed transition piece (ID=332 mm to 212.7 mm) (after TDI)
1099 conpar[0] = (50.0-2.92-1.89)/2.;
1100 conpar[1] = 33.2/2.;
1101 conpar[2] = 33.8/2.;
1102 conpar[3] = 21.27/2.;
1103 conpar[4] = 21.87/2.;
1104 gMC->Gsvolu("QA11", "CONE", idtmed[7], conpar, 5);
1105 gMC->Gspos("QA11", 1, "ZDCA", 4.32-3.8, 0., conpar[0]+2.92+zd2, irotpipe5, "ONLY");
1106 // Ch.debug
1107 //printf(" QA11 skewed CONE from z = %f to z =%f\n",zd2,2*conpar[0]+2.92+1.89+zd2);
1108
1109 zd2 += 2.*conpar[0]+2.92+1.89;
1110
1111 // The following tube ID 212.7 mm
1112 // represents VMBGA (400 mm) + VCDWE (300 mm) + VMBGA (400 mm) +
1113 // BTVTS (600 mm) + VMLGB (400 mm)
1114 tubpar[0] = 21.27/2.;
1115 tubpar[1] = 21.87/2.;
1116 tubpar[2] = 210.0/2.;
1117 gMC->Gsvolu("QA12", "TUBE", idtmed[7], tubpar, 3);
1118 gMC->Gspos("QA12", 1, "ZDCA", 4., 0., tubpar[2]+zd2, 0, "ONLY");
1119 // Ch.debug
1120 //printf(" QA12 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1121
1122 zd2 += 2.*tubpar[2];
1123
1124 // First part of VCTCC
1125 // skewed transition cone from ID=212.7 mm to ID=797 mm
1126 conpar[0] = (121.0-0.37-1.35)/2.;
1127 conpar[1] = 21.27/2.;
1128 conpar[2] = 21.87/2.;
1129 conpar[3] = 79.7/2.;
1130 conpar[4] = 81.3/2.;
1131 gMC->Gsvolu("QA13", "CONE", idtmed[7], conpar, 5);
1132 gMC->Gspos("QA13", 1, "ZDCA", 4.-2., 0., conpar[0]+0.37+zd2, irotpipe3, "ONLY");
1133 // Ch.debug
1134 //printf(" QA13 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+0.37+1.35+zd2);
1135
1136 zd2 += 2.*conpar[0]+0.37+1.35;
1137
1138 // The following tube ID 797 mm --- (volume QA16)
1139 // represents the second part of VCTCC (4272 mm) +
1140 // 4 x VCDGA (4 x 4272 mm) +
1141 // the first part of VCTCR (850 mm)
1142 tubpar[0] = 79.7/2.;
1143 tubpar[1] = 81.3/2.;
1144 tubpar[2] = 2221./2.;
1145 gMC->Gsvolu("QA14", "TUBE", idtmed[7], tubpar, 3);
1146 gMC->Gspos("QA14", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1147 // Ch.debug
1148 //printf(" QA14 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1149
1150 zd2 += 2.*tubpar[2];
1151
1152 // Second part of VCTCR
1153 // Transition from ID=797 mm to ID=196 mm:
1154 // in order to simulate the thin window opened in the transition cone
1155 // we divide the transition cone in three cones:
1156 // (1) 8 mm thick (2) 3 mm thick (3) the third 8 mm thick
1157
1158 // (1) 8 mm thick
1159 conpar[0] = 9.09/2.; // 15 degree
1160 conpar[1] = 79.7/2.;
1161 conpar[2] = 81.3/2.; // thickness 8 mm
1162 conpar[3] = 74.82868/2.;
1163 conpar[4] = 76.42868/2.; // thickness 8 mm
1164 gMC->Gsvolu("QA15", "CONE", idtmed[7], conpar, 5);
1165 gMC->Gspos("QA15", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1166 //printf(" QA15 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1167
1168 zd2 += 2.*conpar[0];
1169
1170 // (2) 3 mm thick
1171 conpar[0] = 96.2/2.; // 15 degree
1172 conpar[1] = 74.82868/2.;
1173 conpar[2] = 75.42868/2.; // thickness 3 mm
1174 conpar[3] = 23.19588/2.;
1175 conpar[4] = 23.79588/2.; // thickness 3 mm
1176 gMC->Gsvolu("QA16", "CONE", idtmed[7], conpar, 5);
1177 gMC->Gspos("QA16", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1178 //printf(" QA16 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1179
1180 zd2 += 2.*conpar[0];
1181
1182 // (3) 8 mm thick
1183 conpar[0] = 6.71/2.; // 15 degree
1184 conpar[1] = 23.19588/2.;
1185 conpar[2] = 24.79588/2.;// thickness 8 mm
1186 conpar[3] = 19.6/2.;
1187 conpar[4] = 21.2/2.;// thickness 8 mm
1188 gMC->Gsvolu("QA17", "CONE", idtmed[7], conpar, 5);
1189 gMC->Gspos("QA17", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1190 //printf(" QA19 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1191
1192 zd2 += 2.*conpar[0];
1193
1194 // Third part of VCTCR: tube (ID=196 mm)
1195 tubpar[0] = 19.6/2.;
1196 tubpar[1] = 21.2/2.;
1197 tubpar[2] = 9.55/2.;
1198 gMC->Gsvolu("QA18", "TUBE", idtmed[7], tubpar, 3);
1199 gMC->Gspos("QA18", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1200 // Ch.debug
1201 //printf(" QA18 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1202
1203 zd2 += 2.*tubpar[2];
1204
1205 // Flange (ID=196 mm) (last part of VCTCR and first part of VMZAR)
1206 tubpar[0] = 19.6/2.;
1207 tubpar[1] = 25.3/2.;
1208 tubpar[2] = 4.9/2.;
1209 gMC->Gsvolu("QF01", "TUBE", idtmed[7], tubpar, 3);
1210 gMC->Gspos("QF01", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1211 // Ch.debug
1212 //printf(" QF01 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1213
1214 zd2 += 2.*tubpar[2];
1215
1216 // VMZAR (5 volumes)
1217 tubpar[0] = 20.2/2.;
1218 tubpar[1] = 20.6/2.;
1219 tubpar[2] = 2.15/2.;
1220 gMC->Gsvolu("QA19", "TUBE", idtmed[7], tubpar, 3);
1221 gMC->Gspos("QA19", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1222 // Ch.debug
1223 //printf(" QA19 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1224
1225 zd2 += 2.*tubpar[2];
1226
1227 conpar[0] = 6.9/2.;
1228 conpar[1] = 20.2/2.;
1229 conpar[2] = 20.6/2.;
1230 conpar[3] = 23.9/2.;
1231 conpar[4] = 24.3/2.;
1232 gMC->Gsvolu("QA20", "CONE", idtmed[7], conpar, 5);
1233 gMC->Gspos("QA20", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1234 // Ch.debug
1235 //printf(" QA20 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1236
1237 zd2 += 2.*conpar[0];
1238
1239 tubpar[0] = 23.9/2.;
1240 tubpar[1] = 25.5/2.;
1241 tubpar[2] = 17.0/2.;
1242 gMC->Gsvolu("QA21", "TUBE", idtmed[7], tubpar, 3);
1243 gMC->Gspos("QA21", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1244 // Ch.debug
1245 //printf(" QA21 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1246
1247 zd2 += 2.*tubpar[2];
1248
1249 conpar[0] = 6.9/2.;
1250 conpar[1] = 23.9/2.;
1251 conpar[2] = 24.3/2.;
1252 conpar[3] = 20.2/2.;
1253 conpar[4] = 20.6/2.;
1254 gMC->Gsvolu("QA22", "CONE", idtmed[7], conpar, 5);
1255 gMC->Gspos("QA22", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1256 // Ch.debug
1257 //printf(" QA22 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1258
1259 zd2 += 2.*conpar[0];
1260
1261 tubpar[0] = 20.2/2.;
1262 tubpar[1] = 20.6/2.;
1263 tubpar[2] = 2.15/2.;
1264 gMC->Gsvolu("QA23", "TUBE", idtmed[7], tubpar, 3);
1265 gMC->Gspos("QA23", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1266 // Ch.debug
1267 //printf(" QA23 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1268
1269 zd2 += 2.*tubpar[2];
1270
1271 // Flange (ID=196 mm)(last part of VMZAR and first part of VCTYD)
1272 tubpar[0] = 19.6/2.;
1273 tubpar[1] = 25.3/2.;
1274 tubpar[2] = 4.9/2.;
1275 gMC->Gsvolu("QF02", "TUBE", idtmed[7], tubpar, 3);
1276 gMC->Gspos("QF02", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1277 // Ch.debug
1278 //printf(" QF02 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1279
1280 zd2 += 2.*tubpar[2];
1281
1282 // simulation of the trousers (VCTYB)
1283 tubpar[0] = 19.6/2.;
1284 tubpar[1] = 20.0/2.;
1285 tubpar[2] = 3.9/2.;
1286 gMC->Gsvolu("QA24", "TUBE", idtmed[7], tubpar, 3);
1287 gMC->Gspos("QA24", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1288 // Ch.debug
1289 //printf(" QA24 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1290
1291 zd2 += 2.*tubpar[2];
1292
1293 // transition cone from ID=196. to ID=216.6
1294 conpar[0] = 32.55/2.;
1295 conpar[1] = 19.6/2.;
1296 conpar[2] = 20.0/2.;
1297 conpar[3] = 21.66/2.;
1298 conpar[4] = 22.06/2.;
1299 gMC->Gsvolu("QA25", "CONE", idtmed[7], conpar, 5);
1300 gMC->Gspos("QA25", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1301 // Ch.debug
1302 //printf(" QA25 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1303
1304 zd2 += 2.*conpar[0];
1305
1306 // tube
1307 tubpar[0] = 21.66/2.;
1308 tubpar[1] = 22.06/2.;
1309 tubpar[2] = 28.6/2.;
1310 gMC->Gsvolu("QA26", "TUBE", idtmed[7], tubpar, 3);
1311 gMC->Gspos("QA26", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1312 // Ch.debug
1313 //printf(" QA26 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1314
1315 zd2 += 2.*tubpar[2];
1316
1317 // --------------------------------------------------------
1318 // RECOMBINATION CHAMBER IMPLEMENTED USING TGeo CLASSES!!!!
1319 // author: Chiara (June 2008)
1320 // --------------------------------------------------------
1321 // TRANSFORMATION MATRICES
1322 // Combi transformation:
1323 dx = -3.970000;
1324 dy = 0.000000;
1325 dz = 0.0;
1326 // Rotation:
1327 thx = 84.989100; phx = 0.000000;
1328 thy = 90.000000; phy = 90.000000;
1329 thz = 5.010900; phz = 180.000000;
1330 TGeoRotation *rotMatrix1 = new TGeoRotation("",thx,phx,thy,phy,thz,phz);
1331 // Combi transformation:
1332 dx = -3.970000;
1333 dy = 0.000000;
1334 dz = 0.0;
1335 TGeoCombiTrans *rotMatrix2 = new TGeoCombiTrans("ZDC_c1", dx,dy,dz,rotMatrix1);
1336 rotMatrix2->RegisterYourself();
1337 // Combi transformation:
1338 dx = 3.970000;
1339 dy = 0.000000;
1340 dz = 0.0;
1341 // Rotation:
1342 thx = 95.010900; phx = 0.000000;
1343 thy = 90.000000; phy = 90.000000;
1344 thz = 5.010900; phz = 0.000000;
1345 TGeoRotation *rotMatrix3 = new TGeoRotation("",thx,phx,thy,phy,thz,phz);
1346 TGeoCombiTrans *rotMatrix4 = new TGeoCombiTrans("ZDC_c2", dx,dy,dz,rotMatrix3);
1347 rotMatrix4->RegisterYourself();
1348
1349
1350 // VOLUMES DEFINITION
1351 // Volume: ZDCA
1352 TGeoVolume *pZDCA = gGeoManager->GetVolume("ZDCA");
1353
1354 conpar[0] = (90.1-0.95-0.26)/2.;
1355 conpar[1] = 0.0/2.;
1356 conpar[2] = 21.6/2.;
1357 conpar[3] = 0.0/2.;
1358 conpar[4] = 5.8/2.;
1359 new TGeoCone("QALext", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
1360
1361 conpar[0] = (90.1-0.95-0.26)/2.;
1362 conpar[1] = 0.0/2.;
1363 conpar[2] = 21.2/2.;
1364 conpar[3] = 0.0/2.;
1365 conpar[4] = 5.4/2.;
1366 new TGeoCone("QALint", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
1367
1368 // Outer trousers
1369 TGeoCompositeShape *pOutTrousers = new TGeoCompositeShape("outTrousers", "QALext:ZDC_c1+QALext:ZDC_c2");
1370
1371 // Volume: QALext
1372 //TGeoMedium *medZDCFe = gGeoManager->GetMedium("ZDC_ZIRON");
1373 TGeoVolume *pQALext = new TGeoVolume("QALext",pOutTrousers, medZDCFe);
1374 pQALext->SetLineColor(kBlue);
1375 pQALext->SetVisLeaves(kTRUE);
1376 //
1377 TGeoTranslation *tr1 = new TGeoTranslation(0., 0., (Double_t) conpar[0]+0.95+zd2);
1378 pZDCA->AddNode(pQALext, 1, tr1);
1379 // Inner trousers
1380 TGeoCompositeShape *pIntTrousers = new TGeoCompositeShape("intTrousers", "QALint:ZDC_c1+QALint:ZDC_c2");
1381 // Volume: QALint
1382 //TGeoMedium *medZDCvoid = gGeoManager->GetMedium("ZDC_ZVOID");
1383 TGeoVolume *pQALint = new TGeoVolume("QALint",pIntTrousers, medZDCvoid);
1384 pQALint->SetLineColor(kAzure);
1385 pQALint->SetVisLeaves(kTRUE);
1386 pQALext->AddNode(pQALint, 1);
1387
1388 zd2 += 90.1;
1389
1390 // second section : 2 tubes (ID = 54. OD = 58.)
1391 tubpar[0] = 5.4/2.;
1392 tubpar[1] = 5.8/2.;
1393 tubpar[2] = 40.0/2.;
1394 gMC->Gsvolu("QA27", "TUBE", idtmed[7], tubpar, 3);
1395 gMC->Gspos("QA27", 1, "ZDCA", -15.8/2., 0., tubpar[2]+zd2, 0, "ONLY");
1396 gMC->Gspos("QA27", 2, "ZDCA", 15.8/2., 0., tubpar[2]+zd2, 0, "ONLY");
1397 // Ch.debug
1398 //printf(" QA27 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1399
1400 zd2 += 2.*tubpar[2];
1401
1402 // transition x2zdc to recombination chamber : skewed cone
1403 conpar[0] = (10.-1.)/2.;
1404 conpar[1] = 5.4/2.;
1405 conpar[2] = 5.8/2.;
1406 conpar[3] = 6.3/2.;
1407 conpar[4] = 7.0/2.;
1408 gMC->Gsvolu("QA28", "CONE", idtmed[7], conpar, 5);
1409 gMC->Gspos("QA28", 1, "ZDCA", -7.9-0.175, 0., conpar[0]+0.5+zd2, irotpipe1, "ONLY");
1410 gMC->Gspos("QA28", 2, "ZDCA", 7.9+0.175, 0., conpar[0]+0.5+zd2, irotpipe2, "ONLY");
1411 //printf(" QA28 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+0.2+zd2);
1412
1413 zd2 += 2.*conpar[0]+1.;
1414
1415 // 2 tubes (ID = 63 mm OD=70 mm)
1416 tubpar[0] = 6.3/2.;
1417 tubpar[1] = 7.0/2.;
1418 tubpar[2] = (342.5+498.3)/2.;
1419 gMC->Gsvolu("QA29", "TUBE", idtmed[7], tubpar, 3);
1420 gMC->Gspos("QA29", 1, "ZDCA", -16.5/2., 0., tubpar[2]+zd2, 0, "ONLY");
1421 gMC->Gspos("QA29", 2, "ZDCA", 16.5/2., 0., tubpar[2]+zd2, 0, "ONLY");
1422 //printf(" QA29 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1423
1424 zd2 += 2.*tubpar[2];
1425
1426 // -- Luminometer (Cu box) in front of ZN - side A
1427 boxpar[0] = 8.0/2.;
1428 boxpar[1] = 8.0/2.;
1429 boxpar[2] = fLumiLength/2.;
1430 gMC->Gsvolu("QLUA", "BOX ", idtmed[9], boxpar, 3);
1431 gMC->Gspos("QLUA", 1, "ZDCA", 0., 0., fPosZNA[2]-66.-boxpar[2], 0, "ONLY");
1432 //printf(" QLUA LUMINOMETER from z = %1.2f to z= %1.2f\n\n", fPosZNA[2]-66., fPosZNA[2]-66.-2*boxpar[2]);
1433
1434 //printf(" END OF BEAM PIPE VOLUME DEFINITION AT z = %f\n",zd2);
1435
1436
1437 // ----------------------------------------------------------------
1438 // -- MAGNET DEFINITION -> LHC OPTICS 6.5
1439 // ----------------------------------------------------------------
1440 // ***************************************************************
1441 // SIDE C - RB26 (dimuon side)
1442 // ***************************************************************
1443 // -- COMPENSATOR DIPOLE (MBXW)
1444 zCorrDip = 1972.5;
1445
1446 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1447 tubpar[0] = 0.;
1448 tubpar[1] = 3.14;
1449 tubpar[2] = 153./2.;
1450 gMC->Gsvolu("MBXW", "TUBE", idtmed[11], tubpar, 3);
1451
1452 // -- YOKE
1453 tubpar[0] = 4.5;
1454 tubpar[1] = 55.;
1455 tubpar[2] = 153./2.;
1456 gMC->Gsvolu("YMBX", "TUBE", idtmed[7], tubpar, 3);
1457
1458 gMC->Gspos("MBXW", 1, "ZDCC", 0., 0., -tubpar[2]-zCorrDip, 0, "ONLY");
1459 gMC->Gspos("YMBX", 1, "ZDCC", 0., 0., -tubpar[2]-zCorrDip, 0, "ONLY");
1460
1461
1462 // -- INNER TRIPLET
1463 zInnTrip = 2296.5;
1464
1465 // -- DEFINE MQXL AND MQX QUADRUPOLE ELEMENT
1466 // -- MQXL
1467 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1468 tubpar[0] = 0.;
1469 tubpar[1] = 3.14;
1470 tubpar[2] = 637./2.;
1471 gMC->Gsvolu("MQXL", "TUBE", idtmed[11], tubpar, 3);
1472
1473 // -- YOKE
1474 tubpar[0] = 3.5;
1475 tubpar[1] = 22.;
1476 tubpar[2] = 637./2.;
1477 gMC->Gsvolu("YMQL", "TUBE", idtmed[7], tubpar, 3);
1478
1479 gMC->Gspos("MQXL", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip, 0, "ONLY");
1480 gMC->Gspos("YMQL", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip, 0, "ONLY");
1481
1482 gMC->Gspos("MQXL", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-2400., 0, "ONLY");
1483 gMC->Gspos("YMQL", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-2400., 0, "ONLY");
1484
1485 // -- MQX
1486 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1487 tubpar[0] = 0.;
1488 tubpar[1] = 3.14;
1489 tubpar[2] = 550./2.;
1490 gMC->Gsvolu("MQX ", "TUBE", idtmed[11], tubpar, 3);
1491
1492 // -- YOKE
1493 tubpar[0] = 3.5;
1494 tubpar[1] = 22.;
1495 tubpar[2] = 550./2.;
1496 gMC->Gsvolu("YMQ ", "TUBE", idtmed[7], tubpar, 3);
1497
1498 gMC->Gspos("MQX ", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-908.5, 0, "ONLY");
1499 gMC->Gspos("YMQ ", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-908.5, 0, "ONLY");
1500
1501 gMC->Gspos("MQX ", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-1558.5, 0, "ONLY");
1502 gMC->Gspos("YMQ ", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-1558.5, 0, "ONLY");
1503
1504 // -- SEPARATOR DIPOLE D1
1505 zD1 = 5838.3001;
1506
1507 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1508 tubpar[0] = 0.;
1509 tubpar[1] = 3.46;
1510 tubpar[2] = 945./2.;
1511 gMC->Gsvolu("MD1 ", "TUBE", idtmed[11], tubpar, 3);
1512
1513 // -- Insert horizontal Cu plates inside D1
1514 // -- (to simulate the vacuum chamber)
1515 boxpar[0] = TMath::Sqrt(tubpar[1]*tubpar[1]-(2.98+0.2)*(2.98+0.2)) - 0.05;
1516 boxpar[1] = 0.2/2.;
1517 boxpar[2] = 945./2.;
1518 gMC->Gsvolu("MD1V", "BOX ", idtmed[6], boxpar, 3);
1519 gMC->Gspos("MD1V", 1, "MD1 ", 0., 2.98+boxpar[1], 0., 0, "ONLY");
1520 gMC->Gspos("MD1V", 2, "MD1 ", 0., -2.98-boxpar[1], 0., 0, "ONLY");
1521
1522 // -- YOKE
1523 tubpar[0] = 3.68;
1524 tubpar[1] = 110./2.;
1525 tubpar[2] = 945./2.;
1526 gMC->Gsvolu("YD1 ", "TUBE", idtmed[7], tubpar, 3);
1527
1528 gMC->Gspos("YD1 ", 1, "ZDCC", 0., 0., -tubpar[2]-zD1, 0, "ONLY");
1529 gMC->Gspos("MD1 ", 1, "ZDCC", 0., 0., -tubpar[2]-zD1, 0, "ONLY");
1530 // Ch debug
1531 //printf(" MD1 from z = %1.2f to z= %1.2f cm\n",-zD1, -zD1-2*tubpar[2]);
1532
1533 // -- DIPOLE D2
1534 zD2 = 12167.8;
1535 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1536 tubpar[0] = 0.;
1537 tubpar[1] = 7.5/2.;
1538 tubpar[2] = 945./2.;
1539 gMC->Gsvolu("MD2 ", "TUBE", idtmed[11], tubpar, 3);
1540
1541 // -- YOKE
1542 tubpar[0] = 0.;
1543 tubpar[1] = 55.;
1544 tubpar[2] = 945./2.;
1545 gMC->Gsvolu("YD2 ", "TUBE", idtmed[7], tubpar, 3);
1546
1547 gMC->Gspos("YD2 ", 1, "ZDCC", 0., 0., -tubpar[2]-zD2, 0, "ONLY");
1548 // Ch debug
1549 //printf(" YD2 from z = %1.2f to z= %1.2f cm\n",-zD2, -zD2-2*tubpar[2]);
1550
1551 gMC->Gspos("MD2 ", 1, "YD2 ", -9.4, 0., 0., 0, "ONLY");
1552 gMC->Gspos("MD2 ", 2, "YD2 ", 9.4, 0., 0., 0, "ONLY");
1553
1554 // ***************************************************************
1555 // SIDE A - RB24
1556 // ***************************************************************
1557
1558 // COMPENSATOR DIPOLE (MCBWA) (2nd compensator)
1559 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1560 tubpar[0] = 0.;
1561 tubpar[1] = 3.;
1562 tubpar[2] = 153./2.;
1563 gMC->Gsvolu("MCBW", "TUBE", idtmed[11], tubpar, 3);
1564 gMC->Gspos("MCBW", 1, "ZDCA", 0., 0., tubpar[2]+zCorrDip, 0, "ONLY");
1565
1566 // -- YOKE
1567 tubpar[0] = 4.5;
1568 tubpar[1] = 55.;
1569 tubpar[2] = 153./2.;
1570 gMC->Gsvolu("YMCB", "TUBE", idtmed[7], tubpar, 3);
1571 gMC->Gspos("YMCB", 1, "ZDCA", 0., 0., tubpar[2]+zCorrDip, 0, "ONLY");
1572
1573 // -- INNER TRIPLET
1574 // -- DEFINE MQX1 AND MQX2 QUADRUPOLE ELEMENT
1575 // -- MQX1
1576 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1577 tubpar[0] = 0.;
1578 tubpar[1] = 3.14;
1579 tubpar[2] = 637./2.;
1580 gMC->Gsvolu("MQX1", "TUBE", idtmed[11], tubpar, 3);
1581 gMC->Gsvolu("MQX4", "TUBE", idtmed[11], tubpar, 3);
1582
1583 // -- YOKE
1584 tubpar[0] = 3.5;
1585 tubpar[1] = 22.;
1586 tubpar[2] = 637./2.;
1587 gMC->Gsvolu("YMQ1", "TUBE", idtmed[7], tubpar, 3);
1588
1589 // -- Q1
1590 gMC->Gspos("MQX1", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip, 0, "ONLY");
1591 gMC->Gspos("YMQ1", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip, 0, "ONLY");
1592
1593 // -- BEAM SCREEN FOR Q1
1594 tubpar[0] = 4.78/2.;
1595 tubpar[1] = 5.18/2.;
1596 tubpar[2] = 637./2.;
1597 gMC->Gsvolu("QBS1", "TUBE", idtmed[6], tubpar, 3);
1598 gMC->Gspos("QBS1", 1, "MQX1", 0., 0., 0., 0, "ONLY");
1599 // INSERT VERTICAL PLATE INSIDE Q1
1600 boxpar[0] = 0.2/2.0;
1601 boxpar[1] = TMath::Sqrt(tubpar[0]*tubpar[0]-(1.9+0.2)*(1.9+0.2));
1602 boxpar[2] =637./2.;
1603 gMC->Gsvolu("QBS2", "BOX ", idtmed[6], boxpar, 3);
1604 gMC->Gspos("QBS2", 1, "MQX1", 1.9+boxpar[0], 0., 0., 0, "ONLY");
1605 gMC->Gspos("QBS2", 2, "MQX1", -1.9-boxpar[0], 0., 0., 0, "ONLY");
1606
1607 // -- Q3
1608 gMC->Gspos("MQX4", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+2400., 0, "ONLY");
1609 gMC->Gspos("YMQ1", 2, "ZDCA", 0., 0., tubpar[2]+zInnTrip+2400., 0, "ONLY");
1610
1611 // -- BEAM SCREEN FOR Q3
1612 tubpar[0] = 5.79/2.;
1613 tubpar[1] = 6.14/2.;
1614 tubpar[2] = 637./2.;
1615 gMC->Gsvolu("QBS3", "TUBE", idtmed[6], tubpar, 3);
1616 gMC->Gspos("QBS3", 1, "MQX4", 0., 0., 0., 0, "ONLY");
1617 // INSERT VERTICAL PLATE INSIDE Q3
1618 boxpar[0] = 0.2/2.0;
1619 boxpar[1] = TMath::Sqrt(tubpar[0]*tubpar[0]-(2.405+0.2)*(2.405+0.2));
1620 boxpar[2] =637./2.;
1621 gMC->Gsvolu("QBS4", "BOX ", idtmed[6], boxpar, 3);
1622 gMC->Gspos("QBS4", 1, "MQX4", 2.405+boxpar[0], 0., 0., 0, "ONLY");
1623 gMC->Gspos("QBS4", 2, "MQX4", -2.405-boxpar[0], 0., 0., 0, "ONLY");
1624
1625
1626
1627 // -- MQX2
1628 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1629 tubpar[0] = 0.;
1630 tubpar[1] = 3.14;
1631 tubpar[2] = 550./2.;
1632 gMC->Gsvolu("MQX2", "TUBE", idtmed[11], tubpar, 3);
1633 gMC->Gsvolu("MQX3", "TUBE", idtmed[11], tubpar, 3);
1634
1635 // -- YOKE
1636 tubpar[0] = 3.5;
1637 tubpar[1] = 22.;
1638 tubpar[2] = 550./2.;
1639 gMC->Gsvolu("YMQ2", "TUBE", idtmed[7], tubpar, 3);
1640
1641 // -- BEAM SCREEN FOR Q2
1642 tubpar[0] = 5.79/2.;
1643 tubpar[1] = 6.14/2.;
1644 tubpar[2] = 550./2.;
1645 gMC->Gsvolu("QBS5", "TUBE", idtmed[6], tubpar, 3);
1646 // VERTICAL PLATE INSIDE Q2
1647 boxpar[0] = 0.2/2.0;
1648 boxpar[1] = TMath::Sqrt(tubpar[0]*tubpar[0]-(2.405+0.2)*(2.405+0.2));
1649 boxpar[2] =550./2.;
1650 gMC->Gsvolu("QBS6", "BOX ", idtmed[6], boxpar, 3);
1651
1652 // -- Q2A
1653 gMC->Gspos("MQX2", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+908.5, 0, "ONLY");
1654 gMC->Gspos("QBS5", 1, "MQX2", 0., 0., 0., 0, "ONLY");
1655 gMC->Gspos("QBS6", 1, "MQX2", 2.405+boxpar[0], 0., 0., 0, "ONLY");
1656 gMC->Gspos("QBS6", 2, "MQX2", -2.405-boxpar[0], 0., 0., 0, "ONLY");
1657 gMC->Gspos("YMQ2", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+908.5, 0, "ONLY");
1658
1659
1660 // -- Q2B
1661 gMC->Gspos("MQX3", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+1558.5, 0, "ONLY");
1662 gMC->Gspos("QBS5", 2, "MQX3", 0., 0., 0., 0, "ONLY");
1663 gMC->Gspos("QBS6", 3, "MQX3", 2.405+boxpar[0], 0., 0., 0, "ONLY");
1664 gMC->Gspos("QBS6", 4, "MQX3", -2.405-boxpar[0], 0., 0., 0, "ONLY");
1665 gMC->Gspos("YMQ2", 2, "ZDCA", 0., 0., tubpar[2]+zInnTrip+1558.5, 0, "ONLY");
1666
1667 // -- SEPARATOR DIPOLE D1
1668 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1669 tubpar[0] = 0.;
1670 tubpar[1] = 6.75/2.;//3.375
1671 tubpar[2] = 945./2.;
1672 gMC->Gsvolu("MD1L", "TUBE", idtmed[11], tubpar, 3);
1673
1674 // -- The beam screen tube is provided by the beam pipe in D1 (QA03 volume)
1675 // -- Insert the beam screen horizontal Cu plates inside D1
1676 // -- (to simulate the vacuum chamber)
1677 boxpar[0] = TMath::Sqrt(tubpar[1]*tubpar[1]-(2.885+0.2)*(2.885+0.2));
1678 boxpar[1] = 0.2/2.;
1679 boxpar[2] =945./2.;
1680 gMC->Gsvolu("QBS7", "BOX ", idtmed[6], boxpar, 3);
1681 gMC->Gspos("QBS7", 1, "MD1L", 0., 2.885+boxpar[1],0., 0, "ONLY");
1682 gMC->Gspos("QBS7", 2, "MD1L", 0., -2.885-boxpar[1],0., 0, "ONLY");
1683
1684 // -- YOKE
1685 tubpar[0] = 3.68;
1686 tubpar[1] = 110./2;
1687 tubpar[2] = 945./2.;
1688 gMC->Gsvolu("YD1L", "TUBE", idtmed[7], tubpar, 3);
1689
1690 gMC->Gspos("YD1L", 1, "ZDCA", 0., 0., tubpar[2]+zD1, 0, "ONLY");
1691 gMC->Gspos("MD1L", 1, "ZDCA", 0., 0., tubpar[2]+zD1, 0, "ONLY");
1692
1693 // -- DIPOLE D2
1694 // -- GAP (VACUUM WITH MAGNETIC FIELD)
1695 tubpar[0] = 0.;
1696 tubpar[1] = 7.5/2.; // this has to be checked
1697 tubpar[2] = 945./2.;
1698 gMC->Gsvolu("MD2L", "TUBE", idtmed[11], tubpar, 3);
1699
1700 // -- YOKE
1701 tubpar[0] = 0.;
1702 tubpar[1] = 55.;
1703 tubpar[2] = 945./2.;
1704 gMC->Gsvolu("YD2L", "TUBE", idtmed[7], tubpar, 3);
1705
1706 gMC->Gspos("YD2L", 1, "ZDCA", 0., 0., tubpar[2]+zD2, 0, "ONLY");
1707
1708 gMC->Gspos("MD2L", 1, "YD2L", -9.4, 0., 0., 0, "ONLY");
1709 gMC->Gspos("MD2L", 2, "YD2L", 9.4, 0., 0., 0, "ONLY");
1710
1711 // -- END OF MAGNET DEFINITION
1712}
1713
1714//_____________________________________________________________________________
1715void AliZDCv3::CreateZDC()
1716{
1717 //
1718 // Create the various ZDCs (ZN + ZP)
1719 //
1720
1721 Float_t dimPb[6], dimVoid[6];
1722
1723 Int_t *idtmed = fIdtmed->GetArray();
1724
1725 // Parameters for hadronic calorimeters geometry
1726 // NB -> parameters used ONLY in CreateZDC()
1727 Float_t fGrvZN[3] = {0.03, 0.03, 50.}; // Grooves for neutron detector
1728 Float_t fGrvZP[3] = {0.04, 0.04, 75.}; // Grooves for proton detector
1729 Int_t fDivZN[3] = {11, 11, 0}; // Division for neutron detector
1730 Int_t fDivZP[3] = {7, 15, 0}; // Division for proton detector
1731 Int_t fTowZN[2] = {2, 2}; // Tower for neutron detector
1732 Int_t fTowZP[2] = {4, 1}; // Tower for proton detector
1733
1734 // Parameters for EM calorimeter geometry
1735 // NB -> parameters used ONLY in CreateZDC()
1736 Float_t kDimZEMPb = 0.15*(TMath::Sqrt(2.)); // z-dimension of the Pb slice
1737 Float_t kFibRadZEM = 0.0315; // External fiber radius (including cladding)
1738 Int_t fDivZEM[3] = {92, 0, 20}; // Divisions for EM detector
1739 Float_t fDimZEM[6] = {fZEMLength, 3.5, 3.5, 45., 0., 0.}; // Dimensions of EM detector
1740 Float_t fFibZEM2 = fDimZEM[2]/TMath::Sin(fDimZEM[3]*kDegrad)-kFibRadZEM;
1741 Float_t fFibZEM[3] = {0., 0.0275, fFibZEM2}; // Fibers for EM calorimeter
1742
1743
1744 //-- Create calorimeters geometry
1745
1746 // -------------------------------------------------------------------------------
1747 //--> Neutron calorimeter (ZN)
1748
1749 gMC->Gsvolu("ZNEU", "BOX ", idtmed[1], fDimZN, 3); // Passive material
1750 gMC->Gsvolu("ZNF1", "TUBE", idtmed[3], fFibZN, 3); // Active material
1751 gMC->Gsvolu("ZNF2", "TUBE", idtmed[4], fFibZN, 3);
1752 gMC->Gsvolu("ZNF3", "TUBE", idtmed[4], fFibZN, 3);
1753 gMC->Gsvolu("ZNF4", "TUBE", idtmed[3], fFibZN, 3);
1754 gMC->Gsvolu("ZNG1", "BOX ", idtmed[12], fGrvZN, 3); // Empty grooves
1755 gMC->Gsvolu("ZNG2", "BOX ", idtmed[12], fGrvZN, 3);
1756 gMC->Gsvolu("ZNG3", "BOX ", idtmed[12], fGrvZN, 3);
1757 gMC->Gsvolu("ZNG4", "BOX ", idtmed[12], fGrvZN, 3);
1758
1759 // Divide ZNEU in towers (for hits purposes)
1760
1761 gMC->Gsdvn("ZNTX", "ZNEU", fTowZN[0], 1); // x-tower
1762 gMC->Gsdvn("ZN1 ", "ZNTX", fTowZN[1], 2); // y-tower
1763
1764 //-- Divide ZN1 in minitowers
1765 // fDivZN[0]= NUMBER OF FIBERS PER TOWER ALONG X-AXIS,
1766 // fDivZN[1]= NUMBER OF FIBERS PER TOWER ALONG Y-AXIS
1767 // (4 fibres per minitower)
1768
1769 gMC->Gsdvn("ZNSL", "ZN1 ", fDivZN[1], 2); // Slices
1770 gMC->Gsdvn("ZNST", "ZNSL", fDivZN[0], 1); // Sticks
1771
1772 // --- Position the empty grooves in the sticks (4 grooves per stick)
1773 Float_t dx = fDimZN[0] / fDivZN[0] / 4.;
1774 Float_t dy = fDimZN[1] / fDivZN[1] / 4.;
1775
1776 gMC->Gspos("ZNG1", 1, "ZNST", 0.-dx, 0.+dy, 0., 0, "ONLY");
1777 gMC->Gspos("ZNG2", 1, "ZNST", 0.+dx, 0.+dy, 0., 0, "ONLY");
1778 gMC->Gspos("ZNG3", 1, "ZNST", 0.-dx, 0.-dy, 0., 0, "ONLY");
1779 gMC->Gspos("ZNG4", 1, "ZNST", 0.+dx, 0.-dy, 0., 0, "ONLY");
1780
1781 // --- Position the fibers in the grooves
1782 gMC->Gspos("ZNF1", 1, "ZNG1", 0., 0., 0., 0, "ONLY");
1783 gMC->Gspos("ZNF2", 1, "ZNG2", 0., 0., 0., 0, "ONLY");
1784 gMC->Gspos("ZNF3", 1, "ZNG3", 0., 0., 0., 0, "ONLY");
1785 gMC->Gspos("ZNF4", 1, "ZNG4", 0., 0., 0., 0, "ONLY");
1786
1787 // --- Position the neutron calorimeter in ZDC
1788 // -- Rotation of ZDCs
1789 Int_t irotzdc;
1790 gMC->Matrix(irotzdc, 90., 180., 90., 90., 180., 0.);
1791 //
1792 gMC->Gspos("ZNEU", 1, "ZDCC", fPosZNC[0], fPosZNC[1], fPosZNC[2]-fDimZN[2], irotzdc, "ONLY");
1793 //Ch debug
1794 //printf("\n ZN -> %f < z < %f cm\n",fPosZN[2],fPosZN[2]-2*fDimZN[2]);
1795
1796 // --- Position the neutron calorimeter in ZDC2 (left line)
1797 // -- No Rotation of ZDCs
1798 gMC->Gspos("ZNEU", 2, "ZDCA", fPosZNA[0], fPosZNA[1], fPosZNA[2]+fDimZN[2], 0, "ONLY");
1799 //Ch debug
1800 //printf("\n ZN left -> %f < z < %f cm\n",fPosZNl[2],fPosZNl[2]+2*fDimZN[2]);
1801
1802
1803 // -------------------------------------------------------------------------------
1804 //--> Proton calorimeter (ZP)
1805
1806 gMC->Gsvolu("ZPRO", "BOX ", idtmed[2], fDimZP, 3); // Passive material
1807 gMC->Gsvolu("ZPF1", "TUBE", idtmed[3], fFibZP, 3); // Active material
1808 gMC->Gsvolu("ZPF2", "TUBE", idtmed[4], fFibZP, 3);
1809 gMC->Gsvolu("ZPF3", "TUBE", idtmed[4], fFibZP, 3);
1810 gMC->Gsvolu("ZPF4", "TUBE", idtmed[3], fFibZP, 3);
1811 gMC->Gsvolu("ZPG1", "BOX ", idtmed[12], fGrvZP, 3); // Empty grooves
1812 gMC->Gsvolu("ZPG2", "BOX ", idtmed[12], fGrvZP, 3);
1813 gMC->Gsvolu("ZPG3", "BOX ", idtmed[12], fGrvZP, 3);
1814 gMC->Gsvolu("ZPG4", "BOX ", idtmed[12], fGrvZP, 3);
1815
1816 //-- Divide ZPRO in towers(for hits purposes)
1817
1818 gMC->Gsdvn("ZPTX", "ZPRO", fTowZP[0], 1); // x-tower
1819 gMC->Gsdvn("ZP1 ", "ZPTX", fTowZP[1], 2); // y-tower
1820
1821
1822 //-- Divide ZP1 in minitowers
1823 // fDivZP[0]= NUMBER OF FIBERS ALONG X-AXIS PER MINITOWER,
1824 // fDivZP[1]= NUMBER OF FIBERS ALONG Y-AXIS PER MINITOWER
1825 // (4 fiber per minitower)
1826
1827 gMC->Gsdvn("ZPSL", "ZP1 ", fDivZP[1], 2); // Slices
1828 gMC->Gsdvn("ZPST", "ZPSL", fDivZP[0], 1); // Sticks
1829
1830 // --- Position the empty grooves in the sticks (4 grooves per stick)
1831 dx = fDimZP[0] / fTowZP[0] / fDivZP[0] / 2.;
1832 dy = fDimZP[1] / fTowZP[1] / fDivZP[1] / 2.;
1833
1834 gMC->Gspos("ZPG1", 1, "ZPST", 0.-dx, 0.+dy, 0., 0, "ONLY");
1835 gMC->Gspos("ZPG2", 1, "ZPST", 0.+dx, 0.+dy, 0., 0, "ONLY");
1836 gMC->Gspos("ZPG3", 1, "ZPST", 0.-dx, 0.-dy, 0., 0, "ONLY");
1837 gMC->Gspos("ZPG4", 1, "ZPST", 0.+dx, 0.-dy, 0., 0, "ONLY");
1838
1839 // --- Position the fibers in the grooves
1840 gMC->Gspos("ZPF1", 1, "ZPG1", 0., 0., 0., 0, "ONLY");
1841 gMC->Gspos("ZPF2", 1, "ZPG2", 0., 0., 0., 0, "ONLY");
1842 gMC->Gspos("ZPF3", 1, "ZPG3", 0., 0., 0., 0, "ONLY");
1843 gMC->Gspos("ZPF4", 1, "ZPG4", 0., 0., 0., 0, "ONLY");
1844
1845
1846 // --- Position the proton calorimeter in ZDCC
1847 gMC->Gspos("ZPRO", 1, "ZDCC", fPosZPC[0], fPosZPC[1], fPosZPC[2]-fDimZP[2], irotzdc, "ONLY");
1848 //Ch debug
1849 //printf("\n ZP -> %f < z < %f cm\n",fPosZP[2],fPosZP[2]-2*fDimZP[2]);
1850
1851 // --- Position the proton calorimeter in ZDCA
1852 // --- No rotation
1853 gMC->Gspos("ZPRO", 2, "ZDCA", fPosZPA[0], fPosZPA[1], fPosZPA[2]+fDimZP[2], 0, "ONLY");
1854 //Ch debug
1855 //printf("\n ZP left -> %f < z < %f cm\n",fPosZPl[2],fPosZPl[2]+2*fDimZP[2]);
1856
1857
1858 // -------------------------------------------------------------------------------
1859 // -> EM calorimeter (ZEM)
1860
1861 gMC->Gsvolu("ZEM ", "PARA", idtmed[10], fDimZEM, 6);
1862
1863 Int_t irot1, irot2;
1864 gMC->Matrix(irot1,0.,0.,90.,90.,-90.,0.); // Rotation matrix 1
1865 gMC->Matrix(irot2,180.,0.,90.,fDimZEM[3]+90.,90.,fDimZEM[3]);// Rotation matrix 2
1866 //printf("irot1 = %d, irot2 = %d \n", irot1, irot2);
1867
1868 gMC->Gsvolu("ZEMF", "TUBE", idtmed[3], fFibZEM, 3); // Active material
1869
1870 gMC->Gsdvn("ZETR", "ZEM ", fDivZEM[2], 1); // Tranches
1871
1872 dimPb[0] = kDimZEMPb; // Lead slices
1873 dimPb[1] = fDimZEM[2];
1874 dimPb[2] = fDimZEM[1];
1875 //dimPb[3] = fDimZEM[3]; //controllare
1876 dimPb[3] = 90.-fDimZEM[3]; //originale
1877 dimPb[4] = 0.;
1878 dimPb[5] = 0.;
1879 gMC->Gsvolu("ZEL0", "PARA", idtmed[5], dimPb, 6);
1880 gMC->Gsvolu("ZEL1", "PARA", idtmed[5], dimPb, 6);
1881 gMC->Gsvolu("ZEL2", "PARA", idtmed[5], dimPb, 6);
1882
1883 // --- Position the lead slices in the tranche
1884 Float_t zTran = fDimZEM[0]/fDivZEM[2];
1885 Float_t zTrPb = -zTran+kDimZEMPb;
1886 gMC->Gspos("ZEL0", 1, "ZETR", zTrPb, 0., 0., 0, "ONLY");
1887 gMC->Gspos("ZEL1", 1, "ZETR", kDimZEMPb, 0., 0., 0, "ONLY");
1888
1889 // --- Vacuum zone (to be filled with fibres)
1890 dimVoid[0] = (zTran-2*kDimZEMPb)/2.;
1891 dimVoid[1] = fDimZEM[2];
1892 dimVoid[2] = fDimZEM[1];
1893 dimVoid[3] = 90.-fDimZEM[3];
1894 dimVoid[4] = 0.;
1895 dimVoid[5] = 0.;
1896 gMC->Gsvolu("ZEV0", "PARA", idtmed[10], dimVoid,6);
1897 gMC->Gsvolu("ZEV1", "PARA", idtmed[10], dimVoid,6);
1898
1899 // --- Divide the vacuum slice into sticks along x axis
1900 gMC->Gsdvn("ZES0", "ZEV0", fDivZEM[0], 3);
1901 gMC->Gsdvn("ZES1", "ZEV1", fDivZEM[0], 3);
1902
1903 // --- Positioning the fibers into the sticks
1904 gMC->Gspos("ZEMF", 1,"ZES0", 0., 0., 0., irot2, "ONLY");
1905 gMC->Gspos("ZEMF", 1,"ZES1", 0., 0., 0., irot2, "ONLY");
1906
1907 // --- Positioning the vacuum slice into the tranche
1908 //Float_t displFib = fDimZEM[1]/fDivZEM[0];
1909 gMC->Gspos("ZEV0", 1,"ZETR", -dimVoid[0], 0., 0., 0, "ONLY");
1910 gMC->Gspos("ZEV1", 1,"ZETR", -dimVoid[0]+zTran, 0., 0., 0, "ONLY");
1911
1912 // --- Positioning the ZEM into the ZDC - rotation for 90 degrees
1913 // NB -> ZEM is positioned in ALIC (instead of in ZDC) volume
1914 gMC->Gspos("ZEM ", 1,"ALIC", -fPosZEM[0], fPosZEM[1], fPosZEM[2]+fDimZEM[0], irot1, "ONLY");
1915
1916 // Second EM ZDC (same side w.r.t. IP, just on the other side w.r.t. beam pipe)
1917 gMC->Gspos("ZEM ", 2,"ALIC", fPosZEM[0], fPosZEM[1], fPosZEM[2]+fDimZEM[0], irot1, "ONLY");
1918
1919 // --- Adding last slice at the end of the EM calorimeter
1920 Float_t zLastSlice = fPosZEM[2]+kDimZEMPb+2*fDimZEM[0];
1921 gMC->Gspos("ZEL2", 1,"ALIC", fPosZEM[0], fPosZEM[1], zLastSlice, irot1, "ONLY");
1922 //Ch debug
1923 //printf("\n ZEM lenght = %f cm\n",2*fZEMLength);
1924 //printf("\n ZEM -> %f < z < %f cm\n",fPosZEM[2],fPosZEM[2]+2*fZEMLength+zLastSlice+kDimZEMPb);
1925
1926}
1927
1928//_____________________________________________________________________________
1929void AliZDCv3::CreateMaterials()
1930{
1931 //
1932 // Create Materials for the Zero Degree Calorimeter
1933 //
1934 Float_t dens, ubuf[1], wmat[3], a[3], z[3];
1935
1936 // --- W alloy -> ZN passive material
1937 dens = 17.6;
1938 a[0] = 183.85;
1939 a[1] = 55.85;
1940 a[2] = 58.71;
1941 z[0] = 74.;
1942 z[1] = 26.;
1943 z[2] = 28.;
1944 wmat[0] = .93;
1945 wmat[1] = .03;
1946 wmat[2] = .04;
1947 AliMixture(1, "WALL", a, z, dens, 3, wmat);
1948
1949 // --- Brass (CuZn) -> ZP passive material
1950 dens = 8.48;
1951 a[0] = 63.546;
1952 a[1] = 65.39;
1953 z[0] = 29.;
1954 z[1] = 30.;
1955 wmat[0] = .63;
1956 wmat[1] = .37;
1957 AliMixture(2, "BRASS", a, z, dens, 2, wmat);
1958
1959 // --- SiO2
1960 dens = 2.64;
1961 a[0] = 28.086;
1962 a[1] = 15.9994;
1963 z[0] = 14.;
1964 z[1] = 8.;
1965 wmat[0] = 1.;
1966 wmat[1] = 2.;
1967 AliMixture(3, "SIO2", a, z, dens, -2, wmat);
1968
1969 // --- Lead
1970 ubuf[0] = 1.12;
1971 AliMaterial(5, "LEAD", 207.19, 82., 11.35, .56, 0., ubuf, 1);
1972
1973 // --- Copper (energy loss taken into account)
1974 ubuf[0] = 1.10;
1975 AliMaterial(6, "COPP0", 63.54, 29., 8.96, 1.4, 0., ubuf, 1);
1976
1977 // --- Copper
1978 ubuf[0] = 1.10;
1979 AliMaterial(9, "COPP1", 63.54, 29., 8.96, 1.4, 0., ubuf, 1);
1980
1981 // --- Iron (energy loss taken into account)
1982 ubuf[0] = 1.1;
1983 AliMaterial(7, "IRON0", 55.85, 26., 7.87, 1.76, 0., ubuf, 1);
1984
1985 // --- Iron (no energy loss)
1986 ubuf[0] = 1.1;
1987 AliMaterial(8, "IRON1", 55.85, 26., 7.87, 1.76, 0., ubuf, 1);
1988
1989 // --- Tatalum
1990 ubuf[0] = 1.1;
1991 AliMaterial(13, "TANT", 183.84, 74., 19.3, 0.35, 0., ubuf, 1);
1992
1993 // ---------------------------------------------------------
1994 Float_t aResGas[3]={1.008,12.0107,15.9994};
1995 Float_t zResGas[3]={1.,6.,8.};
1996 Float_t wResGas[3]={0.28,0.28,0.44};
1997 Float_t dResGas = 3.2E-14;
1998
1999 // --- Vacuum (no magnetic field)
2000 AliMixture(10, "VOID", aResGas, zResGas, dResGas, 3, wResGas);
2001
2002 // --- Vacuum (with magnetic field)
2003 AliMixture(11, "VOIM", aResGas, zResGas, dResGas, 3, wResGas);
2004
2005 // --- Air (no magnetic field)
2006 Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
2007 Float_t zAir[4]={6.,7.,8.,18.};
2008 Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
2009 Float_t dAir = 1.20479E-3;
2010 //
2011 AliMixture(12, "Air $", aAir, zAir, dAir, 4, wAir);
2012
2013 // --- Definition of tracking media:
2014
2015 // --- Tantalum = 1 ;
2016 // --- Brass = 2 ;
2017 // --- Fibers (SiO2) = 3 ;
2018 // --- Fibers (SiO2) = 4 ;
2019 // --- Lead = 5 ;
2020 // --- Copper (with high thr.)= 6 ;
2021 // --- Copper (with low thr.)= 9;
2022 // --- Iron (with energy loss) = 7 ;
2023 // --- Iron (without energy loss) = 8 ;
2024 // --- Vacuum (no field) = 10
2025 // --- Vacuum (with field) = 11
2026 // --- Air (no field) = 12
2027
2028 // ****************************************************
2029 // Tracking media parameters
2030 //
2031 Float_t epsil = 0.01; // Tracking precision,
2032 Float_t stmin = 0.01; // Min. value 4 max. step (cm)
2033 Float_t stemax = 1.; // Max. step permitted (cm)
2034 Float_t tmaxfd = 0.; // Maximum angle due to field (degrees)
2035 Float_t tmaxfdv = 0.1; // Maximum angle due to field (degrees)
2036 Float_t deemax = -1.; // Maximum fractional energy loss
2037 Float_t nofieldm = 0.; // Max. field value (no field)
2038 Float_t fieldm = 45.; // Max. field value (with field)
2039 Int_t isvol = 0; // ISVOL =0 -> not sensitive volume
2040 Int_t isvolActive = 1; // ISVOL =1 -> sensitive volume
2041 Int_t inofld = 0; // IFIELD=0 -> no magnetic field
2042 Int_t ifield =2; // IFIELD=2 -> magnetic field defined in AliMagFC.h
2043 // *****************************************************
2044
2045 AliMedium(1, "ZWALL", 1, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2046 AliMedium(2, "ZBRASS",2, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2047 AliMedium(3, "ZSIO2", 3, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2048 AliMedium(4, "ZQUAR", 3, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2049 AliMedium(5, "ZLEAD", 5, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2050 AliMedium(6, "ZCOPP", 6, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2051 AliMedium(7, "ZIRON", 7, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2052 AliMedium(8, "ZIRONN",8, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2053 AliMedium(9, "ZCOPL", 6, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2054 AliMedium(10,"ZVOID",10, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2055 AliMedium(11,"ZVOIM",11, isvol, ifield, fieldm, tmaxfdv, stemax, deemax, epsil, stmin);
2056 AliMedium(12,"ZAIR", 12, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2057 AliMedium(13,"ZTANT",13, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2058 AliMedium(14, "ZIRONT", 7, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2059
2060}
2061
2062//_____________________________________________________________________________
2063void AliZDCv3::AddAlignableVolumes() const
2064{
2065 //
2066 // Create entries for alignable volumes associating the symbolic volume
2067 // name with the corresponding volume path. Needs to be syncronized with
2068 // eventual changes in the geometry.
2069 //
2070 TString volpath1 = "ALIC_1/ZDCC_1/ZNEU_1";
2071 TString volpath2 = "ALIC_1/ZDCC_1/ZPRO_1";
2072 TString volpath3 = "ALIC_1/ZDCA_1/ZNEU_2";
2073 TString volpath4 = "ALIC_1/ZDCA_1/ZPRO_2";
2074
2075 TString symname1="ZDC/NeutronZDC_C";
2076 TString symname2="ZDC/ProtonZDC_C";
2077 TString symname3="ZDC/NeutronZDC_A";
2078 TString symname4="ZDC/ProtonZDC_A";
2079
2080 if(!gGeoManager->SetAlignableEntry(symname1.Data(),volpath1.Data()))
2081 AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname1.Data(),volpath1.Data()));
2082
2083 if(!gGeoManager->SetAlignableEntry(symname2.Data(),volpath2.Data()))
2084 AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname2.Data(),volpath2.Data()));
2085
2086 if(!gGeoManager->SetAlignableEntry(symname3.Data(),volpath3.Data()))
2087 AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname1.Data(),volpath1.Data()));
2088
2089 if(!gGeoManager->SetAlignableEntry(symname4.Data(),volpath4.Data()))
2090 AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname2.Data(),volpath2.Data()));
2091
2092}
2093
2094
2095//_____________________________________________________________________________
2096void AliZDCv3::Init()
2097{
2098 InitTables();
2099 Int_t *idtmed = fIdtmed->GetArray();
2100 //
2101 fMedSensZN = idtmed[1]; // Sensitive volume: ZN passive material
2102 fMedSensZP = idtmed[2]; // Sensitive volume: ZP passive material
2103 fMedSensF1 = idtmed[3]; // Sensitive volume: fibres type 1
2104 fMedSensF2 = idtmed[4]; // Sensitive volume: fibres type 2
2105 fMedSensZEM = idtmed[5]; // Sensitive volume: ZEM passive material
2106 fMedSensTDI = idtmed[6]; // Sensitive volume: TDI Cu shield
2107 fMedSensPI = idtmed[7]; // Sensitive volume: beam pipes
2108 fMedSensLumi = idtmed[9]; // Sensitive volume: luminometer
2109 fMedSensGR = idtmed[12]; // Sensitive volume: air into the grooves
2110 fMedSensVColl = idtmed[13]; // Sensitive volume: collimator jaws
2111}
2112
2113//_____________________________________________________________________________
2114void AliZDCv3::InitTables()
2115{
2116 //
2117 // Read light tables for Cerenkov light production parameterization
2118 //
2119
2120 Int_t k, j;
2121 int read=1;
2122
2123 // --- Reading light tables for ZN
2124 char *lightfName1 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362207s");
2125 FILE *fp1 = fopen(lightfName1,"r");
2126 if(fp1 == NULL){
2127 printf("Cannot open file fp1 \n");
2128 return;
2129 }
2130 else{
2131 for(k=0; k<fNalfan; k++){
2132 for(j=0; j<fNben; j++){
2133 read = fscanf(fp1,"%f",&fTablen[0][k][j]);
2134 if(read==0) AliDebug(3, " Error in reading light table 1");
2135 }
2136 }
2137 fclose(fp1);
2138 }
2139 char *lightfName2 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362208s");
2140 FILE *fp2 = fopen(lightfName2,"r");
2141 if(fp2 == NULL){
2142 printf("Cannot open file fp2 \n");
2143 return;
2144 }
2145 else{
2146 for(k=0; k<fNalfan; k++){
2147 for(j=0; j<fNben; j++){
2148 read = fscanf(fp2,"%f",&fTablen[1][k][j]);
2149 if(read==0) AliDebug(3, " Error in reading light table 2");
2150 }
2151 }
2152 fclose(fp2);
2153 }
2154 char *lightfName3 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362209s");
2155 FILE *fp3 = fopen(lightfName3,"r");
2156 if(fp3 == NULL){
2157 printf("Cannot open file fp3 \n");
2158 return;
2159 }
2160 else{
2161 for(k=0; k<fNalfan; k++){
2162 for(j=0; j<fNben; j++){
2163 read = fscanf(fp3,"%f",&fTablen[2][k][j]);
2164 if(read==0) AliDebug(3, " Error in reading light table 3");
2165 }
2166 }
2167 fclose(fp3);
2168 }
2169 char *lightfName4 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362210s");
2170 FILE *fp4 = fopen(lightfName4,"r");
2171 if(fp4 == NULL){
2172 printf("Cannot open file fp4 \n");
2173 return;
2174 }
2175 else{
2176 for(k=0; k<fNalfan; k++){
2177 for(j=0; j<fNben; j++){
2178 read = fscanf(fp4,"%f",&fTablen[3][k][j]);
2179 if(read==0) AliDebug(3, " Error in reading light table 4");
2180 }
2181 }
2182 fclose(fp4);
2183 }
2184
2185 // --- Reading light tables for ZP and ZEM
2186 char *lightfName5 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552207s");
2187 FILE *fp5 = fopen(lightfName5,"r");
2188 if(fp5 == NULL){
2189 printf("Cannot open file fp5 \n");
2190 return;
2191 }
2192 else{
2193 for(k=0; k<fNalfap; k++){
2194 for(j=0; j<fNbep; j++){
2195 read = fscanf(fp5,"%f",&fTablep[0][k][j]);
2196 if(read==0) AliDebug(3, " Error in reading light table 5");
2197 }
2198 }
2199 fclose(fp5);
2200 }
2201 char *lightfName6 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552208s");
2202 FILE *fp6 = fopen(lightfName6,"r");
2203 if(fp6 == NULL){
2204 printf("Cannot open file fp6 \n");
2205 return;
2206 }
2207 else{
2208 for(k=0; k<fNalfap; k++){
2209 for(j=0; j<fNbep; j++){
2210 read = fscanf(fp6,"%f",&fTablep[1][k][j]);
2211 if(read==0) AliDebug(3, " Error in reading light table 6");
2212 }
2213 }
2214 fclose(fp6);
2215 }
2216 char *lightfName7 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552209s");
2217 FILE *fp7 = fopen(lightfName7,"r");
2218 if(fp7 == NULL){
2219 printf("Cannot open file fp7 \n");
2220 return;
2221 }
2222 else{
2223 for(k=0; k<fNalfap; k++){
2224 for(j=0; j<fNbep; j++){
2225 read = fscanf(fp7,"%f",&fTablep[2][k][j]);
2226 if(read==0) AliDebug(3, " Error in reading light table 7");
2227 }
2228 }
2229 fclose(fp7);
2230 }
2231 char *lightfName8 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552210s");
2232 FILE *fp8 = fopen(lightfName8,"r");
2233 if(fp8 == NULL){
2234 printf("Cannot open file fp8 \n");
2235 return;
2236 }
2237 else{
2238 for(k=0; k<fNalfap; k++){
2239 for(j=0; j<fNbep; j++){
2240 read = fscanf(fp8,"%f",&fTablep[3][k][j]);
2241 if(read==0) AliDebug(3, " Error in reading light table 8");
2242 }
2243 }
2244 fclose(fp8);
2245 }
2246
2247}
2248//_____________________________________________________________________________
2249void AliZDCv3::StepManager()
2250{
2251 //
2252 // Routine called at every step in the Zero Degree Calorimeters
2253 //
2254 Int_t j, vol[2]={0,0}, ibeta=0, ialfa=0, ibe=0, nphe=0;
2255 Float_t hits[13], x[3], xdet[3]={999.,999.,999.}, um[3], ud[3];
2256 Float_t destep=0., be=0., out=0.;
2257 Double_t s[3], p[4];
2258 const char *knamed;
2259 //
2260 for(j=0;j<13;j++) hits[j]=-999.;
2261 //
2262 // --- This part is for no shower developement in beam pipe, TDI, VColl
2263 // If particle interacts with beam pipe, TDI, VColl -> return
2264 if(fNoShower==1 && ((gMC->CurrentMedium() == fMedSensPI) || (gMC->CurrentMedium() == fMedSensTDI) ||
2265 (gMC->CurrentMedium() == fMedSensVColl || (gMC->CurrentMedium() == fMedSensLumi)))){
2266
2267 // If option NoShower is set -> StopTrack
2268
2269 Int_t ipr = 0;
2270 gMC->TrackPosition(s[0],s[1],s[2]);
2271 if(gMC->CurrentMedium() == fMedSensPI){
2272 knamed = gMC->CurrentVolName();
2273 if(!strncmp(knamed,"YMQ",3)){
2274 if(s[2]<0) fpLostITC += 1;
2275 else fpLostITA += 1;
2276 ipr=1;
2277 }
2278 else if(!strncmp(knamed,"YD1",3)){
2279 if(s[2]<0) fpLostD1C += 1;
2280 else fpLostD1A += 1;
2281 ipr=1;
2282 }
2283 }
2284 else if(gMC->CurrentMedium() == fMedSensTDI){
2285 knamed = gMC->CurrentVolName();
2286 if(!strncmp(knamed,"MD1",3)){
2287 if(s[2]<0) fpLostD1C += 1;
2288 else fpLostD1A += 1;
2289 ipr=1;
2290 }
2291 else if(!strncmp(knamed,"QTD",3)) fpLostTDI += 1;
2292 }
2293 else if(gMC->CurrentMedium() == fMedSensVColl){
2294 knamed = gMC->CurrentVolName();
2295 if(!strncmp(knamed,"QCVC",4)) fpcVCollC++;
2296 else if(!strncmp(knamed,"QCVA",4)) fpcVCollA++;
2297 ipr=1;
2298 }
2299 //
2300 //gMC->TrackMomentum(p[0], p[1], p[2], p[3]);
2301 //printf("\t Particle: mass = %1.3f, E = %1.3f GeV, pz = %1.2f GeV -> stopped in volume %s\n",
2302 // gMC->TrackMass(), p[3], p[2], gMC->CurrentVolName());
2303 //
2304 if(ipr!=0){
2305 printf("\n\t **********************************\n");
2306 printf("\t ********** Side C **********\n");
2307 printf("\t # of particles in IT = %d\n",fpLostITC);
2308 printf("\t # of particles in D1 = %d\n",fpLostD1C);
2309 printf("\t # of particles in VColl = %d\n",fpcVCollC);
2310 printf("\t ********** Side A **********\n");
2311 printf("\t # of particles in IT = %d\n",fpLostITA);
2312 printf("\t # of particles in D1 = %d\n",fpLostD1A);
2313 printf("\t # of particles in TDI = %d\n",fpLostTDI);
2314 printf("\t # of particles in VColl = %d\n",fpcVCollA);
2315 printf("\t **********************************\n");
2316 }
2317 gMC->StopTrack();
2318 return;
2319 }
2320
2321 if((gMC->CurrentMedium() == fMedSensZN) || (gMC->CurrentMedium() == fMedSensZP) ||
2322 (gMC->CurrentMedium() == fMedSensGR) || (gMC->CurrentMedium() == fMedSensF1) ||
2323 (gMC->CurrentMedium() == fMedSensF2) || (gMC->CurrentMedium() == fMedSensZEM)){
2324
2325
2326 //Particle coordinates
2327 gMC->TrackPosition(s[0],s[1],s[2]);
2328 for(j=0; j<=2; j++) x[j] = s[j];
2329 hits[0] = x[0];
2330 hits[1] = x[1];
2331 hits[2] = x[2];
2332
2333 // Determine in which ZDC the particle is
2334 knamed = gMC->CurrentVolName();
2335 if(!strncmp(knamed,"ZN",2)){
2336 if(x[2]<0.) vol[0]=1; // ZNC (dimuon side)
2337 else if(x[2]>0.) vol[0]=4; //ZNA
2338 }
2339 else if(!strncmp(knamed,"ZP",2)){
2340 if(x[2]<0.) vol[0]=2; //ZPC (dimuon side)
2341 else if(x[2]>0.) vol[0]=5; //ZPA
2342 }
2343 else if(!strncmp(knamed,"ZE",2)) vol[0]=3; //ZEM
2344
2345 // Determine in which quadrant the particle is
2346 if(vol[0]==1){ //Quadrant in ZNC
2347 // Calculating particle coordinates inside ZNC
2348 xdet[0] = x[0]-fPosZNC[0];
2349 xdet[1] = x[1]-fPosZNC[1];
2350 // Calculating quadrant in ZN
2351 if(xdet[0]<=0.){
2352 if(xdet[1]<=0.) vol[1]=1;
2353 else vol[1]=3;
2354 }
2355 else if(xdet[0]>0.){
2356 if(xdet[1]<=0.) vol[1]=2;
2357 else vol[1]=4;
2358 }
2359 }
2360
2361 else if(vol[0]==2){ //Quadrant in ZPC
2362 // Calculating particle coordinates inside ZPC
2363 xdet[0] = x[0]-fPosZPC[0];
2364 xdet[1] = x[1]-fPosZPC[1];
2365 if(xdet[0]>=fDimZP[0]) xdet[0]=fDimZP[0]-0.01;
2366 if(xdet[0]<=-fDimZP[0]) xdet[0]=-fDimZP[0]+0.01;
2367 // Calculating tower in ZP
2368 Float_t xqZP = xdet[0]/(fDimZP[0]/2.);
2369 for(int i=1; i<=4; i++){
2370 if(xqZP>=(i-3) && xqZP<(i-2)){
2371 vol[1] = i;
2372 break;
2373 }
2374 }
2375 }
2376 //
2377 // Quadrant in ZEM: vol[1] = 1 -> particle in 1st ZEM (placed at x = 8.5 cm)
2378 // vol[1] = 2 -> particle in 2nd ZEM (placed at x = -8.5 cm)
2379 else if(vol[0] == 3){
2380 if(x[0]>0.){
2381 vol[1] = 1;
2382 // Particle x-coordinate inside ZEM1
2383 xdet[0] = x[0]-fPosZEM[0];
2384 }
2385 else{
2386 vol[1] = 2;
2387 // Particle x-coordinate inside ZEM2
2388 xdet[0] = x[0]+fPosZEM[0];
2389 }
2390 xdet[1] = x[1]-fPosZEM[1];
2391 }
2392 //
2393 else if(vol[0]==4){ //Quadrant in ZNA
2394 // Calculating particle coordinates inside ZNA
2395 xdet[0] = x[0]-fPosZNA[0];
2396 xdet[1] = x[1]-fPosZNA[1];
2397 // Calculating quadrant in ZNA
2398 if(xdet[0]>=0.){
2399 if(xdet[1]<=0.) vol[1]=1;
2400 else vol[1]=3;
2401 }
2402 else if(xdet[0]<0.){
2403 if(xdet[1]<=0.) vol[1]=2;
2404 else vol[1]=4;
2405 }
2406 }
2407 //
2408 else if(vol[0]==5){ //Quadrant in ZPA
2409 // Calculating particle coordinates inside ZPA
2410 xdet[0] = x[0]-fPosZPA[0];
2411 xdet[1] = x[1]-fPosZPA[1];
2412 if(xdet[0]>=fDimZP[0]) xdet[0]=fDimZP[0]-0.01;
2413 if(xdet[0]<=-fDimZP[0]) xdet[0]=-fDimZP[0]+0.01;
2414 // Calculating tower in ZP
2415 Float_t xqZP = -xdet[0]/(fDimZP[0]/2.);
2416 for(int i=1; i<=4; i++){
2417 if(xqZP>=(i-3) && xqZP<(i-2)){
2418 vol[1] = i;
2419 break;
2420 }
2421 }
2422 }
2423 if((vol[1]!=1) && (vol[1]!=2) && (vol[1]!=3) && (vol[1]!=4))
2424 AliError(Form(" WRONG tower for det %d: tow %d with xdet=(%f, %f)\n",
2425 vol[0], vol[1], xdet[0], xdet[1]));
2426 // Ch. debug
2427 //printf("\t *** det %d vol %d xdet(%f, %f)\n",vol[0], vol[1], xdet[0], xdet[1]);
2428
2429
2430 // Store impact point and kinetic energy of the ENTERING particle
2431
2432 if(gMC->IsTrackEntering()){
2433 //Particle energy
2434 gMC->TrackMomentum(p[0],p[1],p[2],p[3]);
2435 hits[3] = p[3];
2436
2437 // Impact point on ZDC
2438 // X takes into account the LHC x-axis sign
2439 // which is opposite to positive x on detector front face
2440 // for side A detectors (ZNA and ZPA)
2441 if(vol[0]==4 || vol[0]==5){
2442 hits[4] = -xdet[0];
2443 }
2444 else{
2445 hits[4] = xdet[0];
2446 }
2447 hits[5] = xdet[1];
2448 hits[6] = 0;
2449 hits[7] = 0;
2450 hits[8] = 0;
2451 hits[9] = 0;
2452 //
2453 Int_t curTrackN = gAlice->GetMCApp()->GetCurrentTrackNumber();
2454 TParticle *part = gAlice->GetMCApp()->Particle(curTrackN);
2455 hits[10] = part->GetPdgCode();
2456 //printf("\t PDGCode = %d\n", part->GetPdgCode());
2457 //
2458 Int_t imo = part->GetFirstMother();
2459 if(imo>0){
2460 TParticle * pmot = gAlice->GetMCApp()->Particle(imo);
2461 hits[11] = pmot->GetPdgCode();
2462 }
2463 else hits[11]=0;
2464 //
2465 hits[12] = 1.0e09*gMC->TrackTime(); // in ns!
2466 //printf("\t TrackTime = %f\n", hits[12]);
2467
2468 AddHit(curTrackN, vol, hits);
2469
2470 if(fNoShower==1){
2471 if(vol[0]==1){
2472 fnDetectedC += 1;
2473 if(fnDetectedC==1) printf(" ### Particle in ZNC\n\n");
2474 }
2475 else if(vol[0]==2){
2476 fpDetectedC += 1;
2477 if(fpDetectedC==1) printf(" ### Particle in ZPC\n\n");
2478 }
2479 else if(vol[0]==4){
2480 fnDetectedA += 1;
2481 if(fnDetectedA==1) printf(" ### Particle in ZNA\n\n");
2482 }
2483 else if(vol[0]==5){
2484 fpDetectedA += 1;
2485 if(fpDetectedA==1) printf(" ### Particle in ZPA\n\n");
2486 }
2487 //
2488 //printf("\t Pc: x %1.2f y %1.2f z %1.2f E %1.2f GeV pz = %1.2f GeV in volume %s\n",
2489 // x[0],x[1],x[3],p[3],p[2],gMC->CurrentVolName());
2490 //
2491 gMC->StopTrack();
2492 return;
2493 }
2494 }
2495
2496 // Particle energy loss
2497 if(gMC->Edep() != 0){
2498 hits[9] = gMC->Edep();
2499 hits[7] = 0.;
2500 hits[8] = 0.;
2501 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2502 }
2503 }
2504
2505
2506 // *** Light production in fibres
2507 if((gMC->CurrentMedium() == fMedSensF1) || (gMC->CurrentMedium() == fMedSensF2)){
2508
2509 //Select charged particles
2510 if((destep=gMC->Edep())){
2511
2512 // Particle velocity
2513 Float_t beta = 0.;
2514 gMC->TrackMomentum(p[0],p[1],p[2],p[3]);
2515 Float_t ptot=TMath::Sqrt(p[0]*p[0]+p[1]*p[1]+p[2]*p[2]);
2516 if(p[3] > 0.00001) beta = ptot/p[3];
2517 else return;
2518 if(beta<0.67)return;
2519 else if((beta>=0.67) && (beta<=0.75)) ibeta = 0;
2520 else if((beta>0.75) && (beta<=0.85)) ibeta = 1;
2521 else if((beta>0.85) && (beta<=0.95)) ibeta = 2;
2522 else if(beta>0.95) ibeta = 3;
2523
2524 // Angle between particle trajectory and fibre axis
2525 // 1 -> Momentum directions
2526 um[0] = p[0]/ptot;
2527 um[1] = p[1]/ptot;
2528 um[2] = p[2]/ptot;
2529 gMC->Gmtod(um,ud,2);
2530 // 2 -> Angle < limit angle
2531 Double_t alfar = TMath::ACos(ud[2]);
2532 Double_t alfa = alfar*kRaddeg;
2533 if(alfa>=110.) return;
2534 //
2535 ialfa = Int_t(1.+alfa/2.);
2536
2537 // Distance between particle trajectory and fibre axis
2538 gMC->TrackPosition(s[0],s[1],s[2]);
2539 for(j=0; j<=2; j++){
2540 x[j] = s[j];
2541 }
2542 gMC->Gmtod(x,xdet,1);
2543 if(TMath::Abs(ud[0])>0.00001){
2544 Float_t dcoeff = ud[1]/ud[0];
2545 be = TMath::Abs((xdet[1]-dcoeff*xdet[0])/TMath::Sqrt(dcoeff*dcoeff+1.));
2546 }
2547 else{
2548 be = TMath::Abs(ud[0]);
2549 }
2550
2551 ibe = Int_t(be*1000.+1);
2552
2553 //Looking into the light tables
2554 Float_t charge = gMC->TrackCharge();
2555
2556 if(vol[0]==1 || vol[0]==4) { // (1) ZN fibres
2557 if(ibe>fNben) ibe=fNben;
2558 out = charge*charge*fTablen[ibeta][ialfa][ibe];
2559 nphe = gRandom->Poisson(out);
2560 // Ch. debug
2561 //if(ibeta==3) printf("\t %f \t %f \t %f\n",alfa, be, out);
2562 //printf("\t ibeta = %d, ialfa = %d, ibe = %d -> nphe = %d\n\n",ibeta,ialfa,ibe,nphe);
2563 if(gMC->CurrentMedium() == fMedSensF1){
2564 hits[7] = nphe; //fLightPMQ
2565 hits[8] = 0;
2566 hits[9] = 0;
2567 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2568 }
2569 else{
2570 hits[7] = 0;
2571 hits[8] = nphe; //fLightPMC
2572 hits[9] = 0;
2573 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2574 }
2575 }
2576 else if(vol[0]==2 || vol[0]==5) {// (2) ZP fibres
2577 if(ibe>fNbep) ibe=fNbep;
2578 out = charge*charge*fTablep[ibeta][ialfa][ibe];
2579 nphe = gRandom->Poisson(out);
2580 if(gMC->CurrentMedium() == fMedSensF1){
2581 hits[7] = nphe; //fLightPMQ
2582 hits[8] = 0;
2583 hits[9] = 0;
2584 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2585 }
2586 else{
2587 hits[7] = 0;
2588 hits[8] = nphe; //fLightPMC
2589 hits[9] = 0;
2590 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2591 }
2592 }
2593 else if(vol[0]==3) { // (3) ZEM fibres
2594 if(ibe>fNbep) ibe=fNbep;
2595 out = charge*charge*fTablep[ibeta][ialfa][ibe];
2596 gMC->TrackPosition(s[0],s[1],s[2]);
2597 Float_t xalic[3];
2598 for(j=0; j<3; j++){
2599 xalic[j] = s[j];
2600 }
2601 // z-coordinate from ZEM front face
2602 // NB-> fPosZEM[2]+fZEMLength = -1000.+2*10.3 = 979.69 cm
2603 Float_t z = -xalic[2]+fPosZEM[2]+2*fZEMLength-xalic[1];
2604 //z = xalic[2]-fPosZEM[2]-fZEMLength-xalic[1]*(TMath::Tan(45.*kDegrad));
2605 //printf(" fPosZEM[2]+2*fZEMLength = %f", fPosZEM[2]+2*fZEMLength);
2606 //
2607 // Parametrization for light guide uniformity
2608 // NEW!!! Light guide tilted @ 51 degrees
2609 Float_t guiPar[4]={0.31,-0.0006305,0.01337,0.8895};
2610 Float_t guiEff = guiPar[0]*(guiPar[1]*z*z+guiPar[2]*z+guiPar[3]);
2611 out = out*guiEff;
2612 nphe = gRandom->Poisson(out);
2613 //printf(" out*guiEff = %f nphe = %d", out, nphe);
2614 if(vol[1] == 1){
2615 hits[7] = 0;
2616 hits[8] = nphe; //fLightPMC (ZEM1)
2617 hits[9] = 0;
2618 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2619 }
2620 else{
2621 hits[7] = nphe; //fLightPMQ (ZEM2)
2622 hits[8] = 0;
2623 hits[9] = 0;
2624 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2625 }
2626 }
2627 }
2628 }
2629}