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