hole region in front of PHOS detector: update of sectors' numbers
[u/mrichter/AliRoot.git] / TOF / AliTOFv5T0.cxx
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 $Log$
18 Revision 1.22  2007/10/07 19:40:46  decaro
19 right handling of l2t matrices and alignable entries in case of TOF staging geometry
20
21 Revision 1.21  2007/10/04 13:15:30  arcelli
22 updates to comply with AliTOFGeometryV5 becoming AliTOFGeometry
23
24 Revision 1.20  2007/10/03 18:07:22  arcelli
25 right handling of l2t matrices and alignable entries in case of TOF holes (Annalisa)
26
27 Revision 1.19  2007/10/03 10:41:12  arcelli
28 adding tracking-to-local matrices for new AliTOFcluster
29
30 Revision 1.18  2007/07/27 08:14:48  morsch
31 Write all track references into the same branch.
32
33 Revision 1.17  2007/05/14 14:41:13  decaro
34 Fix a bug in FTOA volume positionig inside BTOF13,14,15,16,17 in case of holes in 11th and 12th sectors
35
36 Revision 1.16  2007/05/04 12:59:26  arcelli
37 Change the TOF SM paths for misalignment (one layer up)
38
39 Revision 1.15  2007/02/19 15:41:55  decaro
40 Coding convention: few corrections
41
42 Revision 1.14  2006/10/17 15:33:14  arcelli
43 Moving some printout from Info to Debug level
44
45 Revision 1.13  2006/10/12 16:35:43  arcelli
46 definition of the alignable volumes symbolic names added
47
48 Revision 1.12  2006/08/22 13:34:46  arcelli
49 removal of effective c++ warnings (C.Zampolli)
50
51 Revision 1.11  2006/07/12 16:03:44  arcelli
52 updates to match the new numbering of the TOF/TRD mother volumes in FRAME (ALICE convention)
53
54 Revision 1.10  2006/05/10 18:40:17  hristov
55 Larger strings for the names
56
57 Revision 1.9  2006/05/04 19:41:42  hristov
58 Possibility for partial TOF geometry (S.Arcelli)
59
60 Revision 1.8  2006/04/20 22:30:50  hristov
61 Coding conventions (Annalisa)
62
63 Revision 1.7  2006/04/16 22:29:05  hristov
64 Coding conventions (Annalisa)
65
66 Revision 1.6  2006/03/20 08:20:35  decaro
67 Al layer: positioning correction
68
69 Revision 1.5  2006/03/20 07:54:20  decaro
70 Correction of some layer thickness
71
72 Revision 1.4  2006/03/13 12:35:44  decaro
73 Suppression of fractional Z warning
74
75 Revision 1.3  2006/02/28 10:38:00  decaro
76 AliTOFGeometry::fAngles, AliTOFGeometry::fHeights,
77 AliTOFGeometry::fDistances arrays: dimension definition in the right
78 location
79
80 Revision 1.2  2006/02/27 18:12:14  decaro
81 Remove in StepManager the dependence of hit indexes from parametrized
82 TOF position
83
84 Revision 1.1  2005/12/15 08:55:33  decaro
85 New TOF geometry description (V5) -G. Cara Romeo and A. De Caro
86
87
88 Revision 0.1 2004 November G. Cara Romeo and A. De Caro
89         Implement new TOF geometry version
90         in order to
91            suppress few volume overlaps
92                 (in the 4th TOF geometry version),
93            insert the realistic strip numbers and positions
94
95 */
96
97 ///////////////////////////////////////////////////////////////////////////////
98 //                                                                           //
99 //  This class contains the functions for version 5 of the Time Of Flight    //
100 //  detector.                                                                //
101 //                                                                           //
102 //  VERSION WITH 5 MODULES AND TILTED STRIPS                                 //
103 //                                                                           //
104 //  FULL COVERAGE VERSION + OPTION for PHOS holes                            //
105 //                                                                           //
106 //                                                                           //
107 //Begin_Html                                                                 //
108 /*                                                                           //
109 <img src="picts/AliTOFv5T0Class.gif">                                        //
110 */                                                                           //
111 //End_Html                                                                   //
112 //                                                                           //
113 ///////////////////////////////////////////////////////////////////////////////
114
115 #include "TBRIK.h"
116 #include "TGeometry.h"
117 #include "TLorentzVector.h"
118 #include "TNode.h"
119 #include "TVirtualMC.h"
120 #include "TGeoManager.h"
121 #include <TGeoMatrix.h>
122 #include <TGeoPhysicalNode.h>
123 #include <TGeoVolume.h>
124
125 #include "AliConst.h"
126 #include "AliLog.h"
127 #include "AliMagF.h"
128 #include "AliMC.h"
129 #include "AliRun.h"
130 #include "AliTrackReference.h"
131
132 #include "AliTOFGeometry.h"
133 #include "AliTOFv5T0.h"
134
135 extern TDirectory *gDirectory;
136 extern TVirtualMC *gMC;
137 extern TGeoManager *gGeoManager;
138
139 extern AliRun *gAlice;
140
141 ClassImp(AliTOFv5T0)
142
143 //_____________________________________________________________________________
144   AliTOFv5T0::AliTOFv5T0():
145   fIdFTOA(-1),
146   fIdFTOB(-1),
147   fIdFTOC(-1),
148   fIdFLTA(-1),
149   fIdFLTB(-1),
150   fIdFLTC(-1),
151   fTOFHoles(kFALSE)
152 {
153   //
154   // Default constructor
155   //
156 }
157  
158 //_____________________________________________________________________________
159 AliTOFv5T0::AliTOFv5T0(const char *name, const char *title):
160   AliTOF(name,title,"tzero"),
161   fIdFTOA(-1),
162   fIdFTOB(-1),
163   fIdFTOC(-1),
164   fIdFLTA(-1),
165   fIdFLTB(-1),
166   fIdFLTC(-1),
167   fTOFHoles(kFALSE)
168 {
169   //
170   // Standard constructor
171   //
172   //
173   // Check that FRAME is there otherwise we have no place where to
174   // put TOF
175
176
177   AliModule* frame = (AliModule*)gAlice->GetModule("FRAME");
178   if(!frame) {
179     AliFatal("TOF needs FRAME to be present");
180   } else{
181     
182     if (fTOFGeometry) delete fTOFGeometry;
183     fTOFGeometry = new AliTOFGeometry();
184
185     if(frame->IsVersion()==1) {
186       AliDebug(1,Form("Frame version %d", frame->IsVersion())); 
187       AliDebug(1,"Full Coverage for TOF");
188       fTOFHoles=false;}    
189     else {
190       AliDebug(1,Form("Frame version %d", frame->IsVersion())); 
191       AliDebug(1,"TOF with Holes for PHOS");
192       fTOFHoles=true;}      
193   }
194   fTOFGeometry->SetHoles(fTOFHoles);
195
196   //AliTOF::fTOFGeometry = fTOFGeometry;
197
198   // Save the geometry
199   TDirectory* saveDir = gDirectory;
200   gAlice->GetRunLoader()->CdGAFile();
201   fTOFGeometry->Write("TOFgeometry");
202   saveDir->cd();
203
204
205
206 //_____________________________________________________________________________
207 void AliTOFv5T0::AddAlignableVolumes() const
208 {
209   //
210   // Create entries for alignable volumes associating the symbolic volume
211   // name with the corresponding volume path. Needs to be syncronized with
212   // eventual changes in the geometry.
213   //
214
215   TString volPath;
216   TString symName;
217
218   TString vpL0  = "ALIC_1/B077_1/BSEGMO";
219   TString vpL1 = "_1/BTOF";
220   TString vpL2 = "_1";
221   TString vpL3 = "/FTOA_0";
222   TString vpL4 = "/FLTA_0/FSTR_";
223
224   TString snSM  = "TOF/sm";
225   TString snSTRIP = "/strip";
226
227   Int_t nSectors=fTOFGeometry->NSectors();
228   Int_t nStrips =fTOFGeometry->NStripA()+
229                  2*fTOFGeometry->NStripB()+
230                  2*fTOFGeometry->NStripC();
231
232   //
233   // The TOF MRPC Strips
234   // The symbolic names are: TOF/sm00/strip01
235   //                           ...
236   //                         TOF/sm17/strip91
237  
238   Int_t imod=0;
239
240   for (Int_t isect = 0; isect < nSectors; isect++) {
241     for (Int_t istr = 1; istr <= nStrips; istr++) {
242
243       //if (fTOFSectors[isect]==-1) continue;
244
245       if (fTOFHoles && (isect==13 || isect==14 || isect==15)) {
246         if (istr<39) {
247           vpL3 = "/FTOB_0";
248           vpL4 = "/FLTB_0/FSTR_";
249         }
250         else if (istr>53) {
251           vpL3 = "/FTOC_0";
252           vpL4 = "/FLTC_0/FSTR_";
253         }
254         else continue;
255       }
256       else {
257         vpL3 = "/FTOA_0";
258         vpL4 = "/FLTA_0/FSTR_";
259       }
260
261       volPath  = vpL0;
262       volPath += isect;
263       volPath += vpL1;
264       volPath += isect;
265       volPath += vpL2;
266       volPath += vpL3;
267       volPath += vpL4;
268       volPath += istr;
269
270       
271       symName  = snSM;
272       symName += Form("%02d",isect);
273       symName += snSTRIP;
274       symName += Form("%02d",istr);
275             
276       AliDebug(2,"--------------------------------------------"); 
277       AliDebug(2,Form("Alignable object %d", imod)); 
278       AliDebug(2,Form("volPath=%s\n",volPath.Data()));
279       AliDebug(2,Form("symName=%s\n",symName.Data()));
280       AliDebug(2,"--------------------------------------------"); 
281               
282       gGeoManager->SetAlignableEntry(symName.Data(),volPath.Data());
283
284       //T2L matrices for alignment
285       TGeoPNEntry *e = gGeoManager->GetAlignableEntry(symName.Data());
286       if (e) {
287         const char *path = e->GetTitle();
288         if (!gGeoManager->cd(path)) {
289           AliFatal(Form("Volume path %s not valid!",path));
290         }
291         TGeoHMatrix *globMatrix = gGeoManager->GetCurrentMatrix();
292         Double_t phi = 20.0 * (isect % 18) + 10.0;
293         TGeoHMatrix *t2l  = new TGeoHMatrix();
294         t2l->RotateZ(phi);
295         t2l->MultiplyLeft(&(globMatrix->Inverse()));
296         e->SetMatrix(t2l);
297       }
298       else {
299         AliError(Form("Alignable entry %s is not valid!",symName.Data()));
300       }
301       imod++;
302     }
303   }
304
305
306   //
307   // The TOF supermodules
308   // The symbolic names are: TOF/sm00
309   //                           ...
310   //                         TOF/sm17
311   //
312   for (Int_t isect = 0; isect < nSectors; isect++) {
313
314     volPath  = vpL0;
315     volPath += isect;
316     volPath += vpL1;
317     volPath += isect;
318     volPath += vpL2;
319
320     symName  = snSM;
321     symName += Form("%02d",isect);
322
323       AliDebug(2,"--------------------------------------------"); 
324       AliDebug(2,Form("Alignable object %d", isect+imod)); 
325       AliDebug(2,Form("volPath=%s\n",volPath.Data()));
326       AliDebug(2,Form("symName=%s\n",symName.Data()));
327       AliDebug(2,"--------------------------------------------"); 
328               
329     gGeoManager->SetAlignableEntry(symName.Data(),volPath.Data());
330
331   }
332   
333 }
334 //____________________________________________________________________________
335 void AliTOFv5T0::BuildGeometry()
336 {
337   //
338   // Build TOF ROOT geometry for the ALICE event display
339   //
340   TNode *node, *top;
341   const int kColorTOF  = 27;
342   
343   TGeometry *globalGeometry = (TGeometry*)gAlice->GetGeometry();
344
345   // Find top TNODE
346   top = globalGeometry->GetNode("alice");
347   
348   // Position the different copies
349   const Float_t krTof  =(fTOFGeometry->Rmax()+fTOFGeometry->Rmin())/2.;
350   const Float_t khTof  = fTOFGeometry->Rmax()-fTOFGeometry->Rmin();
351   const Int_t   kNTof  = fTOFGeometry->NSectors();
352   const Float_t kangle = k2PI/kNTof;
353
354   const Float_t kInterCentrModBorder1 = 49.5;
355   const Float_t kInterCentrModBorder2 = 57.5;
356
357   Float_t ang;
358   
359   // define offset for nodes
360   Float_t zOffsetB = (fTOFGeometry->ZlenA()*0.5 + (kInterCentrModBorder1+kInterCentrModBorder2)*0.5)*0.5;
361   Float_t zOffsetA = 0.;
362   // Define TOF basic volume
363   
364   char nodeName0[16], nodeName1[16], nodeName2[16];
365   char nodeName3[16], nodeName4[16], rotMatNum[16];
366
367   if (fTOFHoles) {
368     new TBRIK("S_TOF_B","TOF box","void",
369               fTOFGeometry->StripLength()*0.5, khTof*0.5, fTOFGeometry->ZlenB()*0.5);
370     new TBRIK("S_TOF_C","TOF box","void",
371               fTOFGeometry->StripLength()*0.5, khTof*0.5, fTOFGeometry->ZlenB()*0.5);
372   }
373   new TBRIK("S_TOF_A","TOF box","void",
374             fTOFGeometry->StripLength()*0.5, khTof*0.5, fTOFGeometry->ZlenA()*0.5);
375   
376   for (Int_t nodeNum=1;nodeNum<19;nodeNum++){
377     
378     if (nodeNum<10) {
379       sprintf(rotMatNum,"rot50%i",nodeNum);
380       sprintf(nodeName0,"FTO00%i",nodeNum);
381       sprintf(nodeName1,"FTO10%i",nodeNum);
382       sprintf(nodeName2,"FTO20%i",nodeNum);
383       sprintf(nodeName3,"FTO30%i",nodeNum);
384       sprintf(nodeName4,"FTO40%i",nodeNum);
385     }
386     if (nodeNum>9) {
387       sprintf(rotMatNum,"rot5%i",nodeNum);
388       sprintf(nodeName0,"FTO0%i",nodeNum);
389       sprintf(nodeName1,"FTO1%i",nodeNum);
390       sprintf(nodeName2,"FTO2%i",nodeNum);
391       sprintf(nodeName3,"FTO3%i",nodeNum);
392       sprintf(nodeName4,"FTO4%i",nodeNum);
393     }
394     
395     new TRotMatrix(rotMatNum,rotMatNum,90,-20*nodeNum,90,90-20*nodeNum,0,0);
396     ang = (4.5-nodeNum) * kangle;
397
398     if (fTOFHoles) {   
399       top->cd();
400       node = new TNode(nodeName2,nodeName2,"S_TOF_B", krTof*TMath::Cos(ang), krTof*TMath::Sin(ang), zOffsetB,rotMatNum);
401       node->SetLineColor(kColorTOF);
402       fNodes->Add(node);
403       
404       top->cd();
405       node = new TNode(nodeName3,nodeName3,"S_TOF_C", krTof*TMath::Cos(ang), krTof*TMath::Sin(ang),-zOffsetB,rotMatNum);
406       node->SetLineColor(kColorTOF);
407       fNodes->Add(node);
408     }
409
410     top->cd();
411     node = new TNode(nodeName4,nodeName4,"S_TOF_A", krTof*TMath::Cos(ang), krTof*TMath::Sin(ang), zOffsetA,rotMatNum);
412     node->SetLineColor(kColorTOF);
413     fNodes->Add(node);
414   } // end loop on nodeNum
415
416 }
417
418 //_____________________________________________________________________________
419 void AliTOFv5T0::CreateGeometry()
420 {
421   //
422   // Create geometry for Time Of Flight version 0
423   //
424   //Begin_Html
425   /*
426     <img src="picts/AliTOFv5T0.gif">
427   */
428   //End_Html
429   //
430   // Creates common geometry
431   //
432   AliTOF::CreateGeometry();
433 }
434  
435
436 //_____________________________________________________________________________
437 void AliTOFv5T0::TOFpc(Float_t xtof,  Float_t ytof, Float_t zlenA,
438                        Float_t zlenB)
439 {
440   //
441   // Definition of the Time Of Fligh Resistive Plate Chambers
442   //
443
444   const Float_t kPi = TMath::Pi();
445
446   const Float_t kInterCentrModBorder1 = 49.5;
447   const Float_t kInterCentrModBorder2 = 57.5;
448   const Float_t kExterInterModBorder1 = 196.0;
449   const Float_t kExterInterModBorder2 = 203.5;
450
451   const Float_t kLengthExInModBorder  = 4.7;
452   const Float_t kLengthInCeModBorder  = 7.0;
453
454   // Al layers over all internal module walls (cm)
455   const Float_t khAlWall = 0.03;
456
457   // module wall thickness (cm)
458   const Float_t kModuleWallThickness = 0.3;
459
460   // Al honeycomb layer between strips and cards (cm)
461   const Float_t kHoneycombLayerThickness = 1.5;
462
463   AliDebug(2,Form("zlenA*0.5 = %d", zlenA*0.5));
464   AliDebug(1, "************************* TOF geometry **************************");
465   
466   // Definition of the Time Of Fligh Resistive Plate Chambers
467   // xFLT, yFLT, zFLT - sizes of TOF modules (large)
468   
469   Float_t  xcoor, ycoor, zcoor;
470   Float_t  par[3];
471   Int_t    *idtmed = fIdtmed->GetArray()-499;
472   Int_t    idrotm[100];
473
474   par[0] =  xtof * 0.5;
475   par[1] =  ytof * 0.5;
476   par[2] = zlenA * 0.5;
477   gMC->Gsvolu("FTOA", "BOX ", idtmed[503], par, 3);  // fibre glass
478   
479   if (fTOFHoles) {
480     par[0] =  xtof * 0.5;
481     par[1] =  ytof * 0.5;
482     par[2] = (zlenA*0.5 - kInterCentrModBorder1)*0.5;
483     gMC->Gsvolu("FTOB", "BOX ", idtmed[503], par, 3);  // fibre glass
484     gMC->Gsvolu("FTOC", "BOX ", idtmed[503], par, 3);  // fibre glass
485   }
486
487   // Positioning of fibre glass modules (FTOA, FTOB and FTOC)
488   
489   //AliMatrix(idrotm[0], 90.,  0., 0., 0., 90.,-90.);
490   AliMatrix(idrotm[0], 90.,  0., 0., 0., 90.,270.);
491
492   xcoor = 0.;
493   ycoor = 0.;
494   zcoor = 0.;
495   for(Int_t isec=0;isec<18;isec++){
496     if(fTOFSectors[isec]==-1)continue;
497     char name[16];
498     sprintf(name, "BTOF%d",isec);
499     if (fTOFHoles && (isec==13 || isec==14 || isec==15)) {
500     //    if (fTOFHoles && (isec==16||isec==17)) { \\Old 6h convention
501       xcoor = 0.;
502       ycoor = (zlenA*0.5 + kInterCentrModBorder1)*0.5;
503       zcoor = 0.;
504       gMC->Gspos("FTOB", 0, name, xcoor, ycoor, zcoor, idrotm[0], "ONLY");
505       gMC->Gspos("FTOC", 0, name, xcoor,-ycoor, zcoor, idrotm[0], "ONLY");
506     }
507     else {
508       xcoor = 0.;
509       ycoor = 0.;
510       zcoor = 0.;
511       gMC->Gspos("FTOA", 0,name, xcoor, ycoor, zcoor, idrotm[0], "ONLY");
512     }
513
514   }
515   // Large not sensitive volumes with Insensitive Freon (FLTA, FLTB and FLTC)
516   
517   Float_t xFLT, yFLT, zFLTA;
518   
519   xFLT  = xtof  - kModuleWallThickness*2.;
520   yFLT  = ytof  - kModuleWallThickness*2.;
521   zFLTA = zlenA - kModuleWallThickness*2.;
522   
523   par[0] = xFLT*0.5;
524   par[1] = yFLT*0.5;
525   par[2] = zFLTA*0.5;
526   gMC->Gsvolu("FLTA", "BOX ", idtmed[507], par, 3); //  Freon mix
527
528   xcoor = 0.;
529   ycoor = 0.;
530   zcoor = 0.;
531   gMC->Gspos ("FLTA", 0, "FTOA", xcoor, ycoor, zcoor, 0, "ONLY");
532
533   if (fTOFHoles) {
534     par[0] = xFLT*0.5;
535     par[1] = yFLT*0.5;
536     par[2] = (zlenA*0.5 - kInterCentrModBorder1-kModuleWallThickness)*0.5;
537     gMC->Gsvolu("FLTB", "BOX ", idtmed[507], par, 3); // Freon mix
538     gMC->Gsvolu("FLTC", "BOX ", idtmed[507], par, 3); // Freon mix
539
540     xcoor = 0.;
541     ycoor = 0.;
542     zcoor = kModuleWallThickness*0.5;
543     gMC->Gspos ("FLTB", 0, "FTOB", xcoor, ycoor, zcoor, 0, "ONLY");
544     gMC->Gspos ("FLTC", 0, "FTOC", xcoor, ycoor,-zcoor, 0, "ONLY");
545   }
546
547   // Layer of Aluminum before detector (FALA, FALB and FALC)
548
549   par[0] = xFLT*0.5;
550   par[1] = khAlWall*0.5;
551   par[2] = kInterCentrModBorder1 - (kModuleWallThickness + khAlWall);
552   gMC->Gsvolu("FALA", "BOX ", idtmed[505], par, 3); // Alluminium
553
554   xcoor = 0.;
555   ycoor = (-yFLT + khAlWall)*0.5;
556   zcoor = 0.;
557   gMC->Gspos ("FALA", 0, "FLTA", xcoor, ycoor, zcoor, 0, "ONLY");
558
559   par[0] = xFLT*0.5;
560   par[1] = khAlWall*0.5;
561   par[2] = (kExterInterModBorder2 - kInterCentrModBorder1 - 2.*(kModuleWallThickness + khAlWall))*0.5;
562   gMC->Gsvolu("FALB", "BOX ", idtmed[505], par, 3); // Alluminium
563
564   xcoor = 0.;
565   ycoor = (-yFLT + khAlWall)*0.5;
566   zcoor = (kExterInterModBorder2 + kInterCentrModBorder1)*0.5;
567   gMC->Gspos ("FALB", 1, "FLTA", xcoor, ycoor, zcoor, 0, "ONLY");
568   gMC->Gspos ("FALB", 2, "FLTA", xcoor, ycoor,-zcoor, 0, "ONLY");
569
570   par[0] = xFLT*0.5;
571   par[1] = khAlWall*0.5;
572   par[2] = (zlenA*0.5 - kExterInterModBorder2 - 2.*(kModuleWallThickness + khAlWall))*0.5;
573   gMC->Gsvolu("FALC", "BOX ", idtmed[505], par, 3); // Alluminium
574
575   xcoor = 0.;
576   ycoor = (-yFLT + khAlWall)*0.5;
577   zcoor = (kExterInterModBorder2+zlenA*0.5)*0.5;
578   gMC->Gspos ("FALC", 1, "FLTA", xcoor, ycoor, zcoor, 0, "ONLY");
579   gMC->Gspos ("FALC", 2, "FLTA", xcoor, ycoor,-zcoor, 0, "ONLY");
580
581   if (fTOFHoles) {
582     xcoor = 0.;
583     ycoor = (-yFLT + khAlWall)*0.5;
584     zcoor = (zlenA*0.5 - kExterInterModBorder2)*0.5 - kModuleWallThickness*0.5;
585     gMC->Gspos ("FALB", 1, "FLTB", xcoor, ycoor, zcoor, 0, "ONLY");
586     gMC->Gspos ("FALB", 2, "FLTC", xcoor, ycoor,-zcoor, 0, "ONLY");
587
588     xcoor = 0.;
589     ycoor = (-yFLT + khAlWall)*0.5;
590     zcoor = (kExterInterModBorder2 - kInterCentrModBorder1)*0.5 + kModuleWallThickness*0.5;
591     gMC->Gspos ("FALC", 1, "FLTB", xcoor, ycoor,-zcoor, 0, "ONLY");
592     gMC->Gspos ("FALC", 2, "FLTC", xcoor, ycoor, zcoor, 0, "ONLY");
593   }
594
595   Float_t y0, alpha, tgal, beta, tgbe, trpa[11];
596
597   // Fibre glass walls between central and intermediate modules (FWZ1 and FWZ2; holes -> FZ1B, FZ1C, FZ2B)
598
599   tgal = (yFLT*0.5 - 2.*kLengthInCeModBorder)/(kInterCentrModBorder2 - kInterCentrModBorder1);
600   alpha = TMath::ATan(tgal);
601   beta = (kPi*0.5 - alpha)*0.5;
602   tgbe = TMath::Tan(beta);
603   trpa[0]  = xFLT*0.5;
604   trpa[1]  = 0.;
605   trpa[2]  = 0.;
606   trpa[3]  = kModuleWallThickness;
607   trpa[4]  = (kLengthInCeModBorder - kModuleWallThickness*tgbe)*0.5;
608   trpa[5]  = (kLengthInCeModBorder + kModuleWallThickness*tgbe)*0.5;
609   trpa[6]  = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
610   trpa[7]  = kModuleWallThickness;
611   trpa[8]  = (kLengthInCeModBorder - kModuleWallThickness*tgbe)*0.5;
612   trpa[9]  = (kLengthInCeModBorder + kModuleWallThickness*tgbe)*0.5;
613   trpa[10] = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
614   gMC->Gsvolu("FWZ1","TRAP", idtmed[503], trpa, 11);   // fibre glass
615
616   AliMatrix (idrotm[1],90., 90.,180.,0.,90.,180.);
617   AliMatrix (idrotm[4],90., 90.,  0.,0.,90.,  0.);
618
619   xcoor = 0.;
620   ycoor = -yFLT*0.5 + kLengthInCeModBorder*0.5;
621   zcoor = kInterCentrModBorder1;
622   gMC->Gspos("FWZ1", 1,"FLTA", xcoor, ycoor, zcoor,idrotm[1],"ONLY");
623   gMC->Gspos("FWZ1", 2,"FLTA", xcoor, ycoor,-zcoor,idrotm[4],"ONLY");
624
625   if (fTOFHoles) {
626     y0 = kLengthInCeModBorder - kModuleWallThickness*0.5*tgbe;
627     trpa[0]  = xFLT*0.5;
628     trpa[1]  = 0.;
629     trpa[2]  = 0.;
630     trpa[3]  = kModuleWallThickness*0.5;
631     trpa[4]  = (y0 - kModuleWallThickness*0.5*tgbe)*0.5;
632     trpa[5]  = (y0 + kModuleWallThickness*0.5*tgbe)*0.5;
633     trpa[6]  = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
634     trpa[7]  = kModuleWallThickness*0.5;
635     trpa[8]  = (y0 - kModuleWallThickness*0.5*tgbe)*0.5;
636     trpa[9]  = (y0 + kModuleWallThickness*0.5*tgbe)*0.5;
637     trpa[10] = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
638     gMC->Gsvolu("FZ1B","TRAP", idtmed[503], trpa, 11);   // fibre glass
639
640     xcoor = 0.;
641     ycoor = -yFLT*0.5 + kLengthInCeModBorder*0.5 - kModuleWallThickness*0.25*tgbe;
642     zcoor = -kInterCentrModBorder1 + (zlenA*0.5 + kInterCentrModBorder1)*0.5 - kModuleWallThickness;
643     gMC->Gspos("FZ1B", 1,"FLTB", xcoor, ycoor, zcoor,idrotm[4],"ONLY");
644     gMC->Gspos("FZ1B", 2,"FLTC", xcoor, ycoor,-zcoor,idrotm[1],"ONLY");
645   }
646
647   AliMatrix (idrotm[2],90.,270.,  0.,0.,90.,180.);
648   AliMatrix (idrotm[5],90.,270.,180.,0.,90.,  0.);
649
650   xcoor = 0.;
651   ycoor = -kLengthInCeModBorder*0.5;
652   zcoor = kInterCentrModBorder2;
653   gMC->Gspos("FWZ1", 3,"FLTA", xcoor, ycoor, zcoor,idrotm[2],"ONLY");
654   gMC->Gspos("FWZ1", 4,"FLTA", xcoor, ycoor,-zcoor,idrotm[5],"ONLY");
655
656   if (fTOFHoles) {
657     y0 = kLengthInCeModBorder + kModuleWallThickness*0.5*tgbe;
658     trpa[0]  = xFLT*0.5;
659     trpa[1]  = 0.;
660     trpa[2]  = 0.;
661     trpa[3]  = kModuleWallThickness*0.5;
662     trpa[4]  = (y0 - kModuleWallThickness*0.5*tgbe)*0.5;
663     trpa[5]  = (y0 + kModuleWallThickness*0.5*tgbe)*0.5;
664     trpa[6]  = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
665     trpa[7]  = kModuleWallThickness*0.5;
666     trpa[8]  = (y0 - kModuleWallThickness*0.5*tgbe)*0.5;
667     trpa[9]  = (y0 + kModuleWallThickness*0.5*tgbe)*0.5;
668     trpa[10] = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
669     gMC->Gsvolu("FZ1C","TRAP", idtmed[503], trpa, 11);   // fibre glass
670
671     xcoor = 0.;
672     ycoor = -kLengthInCeModBorder*0.5 - kModuleWallThickness*0.25*tgbe;
673     zcoor = -kInterCentrModBorder2 + (zlenA*0.5 + kInterCentrModBorder1)*0.5 - kModuleWallThickness;
674     gMC->Gspos("FZ1C", 1,"FLTB", xcoor, ycoor, zcoor,idrotm[5],"ONLY");
675     gMC->Gspos("FZ1C", 2,"FLTC", xcoor, ycoor,-zcoor,idrotm[2],"ONLY");
676   }
677
678   trpa[0] = 0.5*(kInterCentrModBorder2 - kInterCentrModBorder1)/TMath::Cos(alpha);
679   trpa[1] = kModuleWallThickness;
680   trpa[2] = xFLT*0.5;
681   trpa[3] = -beta*kRaddeg;
682   trpa[4] = 0.;
683   trpa[5] = 0.;
684   gMC->Gsvolu("FWZ2","PARA", idtmed[503], trpa, 6);    // fibre glass
685
686   AliMatrix (idrotm[3],     alpha*kRaddeg,90.,90.+alpha*kRaddeg,90.,90.,180.);
687   AliMatrix (idrotm[6],180.-alpha*kRaddeg,90.,90.-alpha*kRaddeg,90.,90.,  0.);
688
689   xcoor = 0.;
690   ycoor = -yFLT*0.25;
691   zcoor = (kInterCentrModBorder2 + kInterCentrModBorder1)*0.5;
692   gMC->Gspos("FWZ2", 1,"FLTA", xcoor, ycoor, zcoor,idrotm[3],"ONLY");
693   gMC->Gspos("FWZ2", 2,"FLTA", xcoor, ycoor,-zcoor,idrotm[6],"ONLY");
694
695   if (fTOFHoles) {
696     trpa[0] = 0.5*(kInterCentrModBorder2 - kInterCentrModBorder1)/TMath::Cos(alpha);
697     trpa[1] = kModuleWallThickness*0.5;
698     trpa[2] = xFLT*0.5;
699     trpa[3] = -beta*kRaddeg;
700     trpa[4] = 0.;
701     trpa[5] = 0.;
702     gMC->Gsvolu("FZ2B","PARA", idtmed[503], trpa, 6);    // fibre glass
703
704     xcoor = 0.;
705     ycoor = -yFLT*0.25 - kModuleWallThickness*0.5*tgbe;
706     zcoor = -(kInterCentrModBorder2 + kInterCentrModBorder1)*0.5 + (zlenA*0.5 + kInterCentrModBorder1)*0.5 - kModuleWallThickness;
707     gMC->Gspos("FZ2B", 1,"FLTB", xcoor, ycoor, zcoor,idrotm[6],"ONLY");
708     gMC->Gspos("FZ2B", 2,"FLTC", xcoor, ycoor,-zcoor,idrotm[3],"ONLY");
709   }
710
711   // Fibre glass walls between intermediate and lateral modules (FWZ3 and FWZ4)
712
713   tgal = (yFLT*0.5 - 2.*kLengthExInModBorder)/(kExterInterModBorder2 - kExterInterModBorder1);
714   alpha = TMath::ATan(tgal);
715   beta = (kPi*0.5 - alpha)*0.5;
716   tgbe = TMath::Tan(beta);
717   trpa[0]  = xFLT*0.5;
718   trpa[1]  = 0.;
719   trpa[2]  = 0.;
720   trpa[3]  = kModuleWallThickness;
721   trpa[4]  = (kLengthExInModBorder - kModuleWallThickness*tgbe)*0.5;
722   trpa[5]  = (kLengthExInModBorder + kModuleWallThickness*tgbe)*0.5;
723   trpa[6]  = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
724   trpa[7]  = kModuleWallThickness;
725   trpa[8]  = (kLengthExInModBorder - kModuleWallThickness*tgbe)*0.5;
726   trpa[9]  = (kLengthExInModBorder + kModuleWallThickness*tgbe)*0.5;
727   trpa[10] = TMath::ATan(tgbe*0.5)*kRaddeg; //TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*kRaddeg;
728   gMC->Gsvolu("FWZ3","TRAP", idtmed[503], trpa, 11);    // fibre glass
729
730   xcoor = 0.;
731   ycoor = -kLengthExInModBorder*0.5;
732   zcoor = kExterInterModBorder1;
733   gMC->Gspos("FWZ3", 1,"FLTA", xcoor, ycoor, zcoor,idrotm[5],"ONLY");
734   gMC->Gspos("FWZ3", 2,"FLTA", xcoor, ycoor,-zcoor,idrotm[2],"ONLY");
735
736   if (fTOFHoles) {
737     xcoor = 0.;
738     ycoor = -kLengthExInModBorder*0.5;
739     zcoor = -kExterInterModBorder1 + (zlenA*0.5 + kInterCentrModBorder1 - kModuleWallThickness)*0.5;
740     gMC->Gspos("FWZ3", 5,"FLTB", xcoor, ycoor, zcoor,idrotm[2],"ONLY");
741     gMC->Gspos("FWZ3", 6,"FLTC", xcoor, ycoor,-zcoor,idrotm[5],"ONLY");
742   }
743
744   xcoor = 0.;
745   ycoor = -yFLT*0.5 + kLengthExInModBorder*0.5;
746   zcoor = kExterInterModBorder2;
747   gMC->Gspos("FWZ3", 3,"FLTA", xcoor, ycoor, zcoor,idrotm[4],"ONLY");
748   gMC->Gspos("FWZ3", 4,"FLTA", xcoor, ycoor,-zcoor,idrotm[1],"ONLY");
749
750   if (fTOFHoles) {
751     xcoor = 0.;
752     ycoor = -yFLT*0.5 + kLengthExInModBorder*0.5;
753     zcoor = -kExterInterModBorder2 + (zlenA*0.5 + kInterCentrModBorder1 - kModuleWallThickness)*0.5;
754     gMC->Gspos("FWZ3", 7,"FLTB", xcoor, ycoor, zcoor,idrotm[1],"ONLY");
755     gMC->Gspos("FWZ3", 8,"FLTC", xcoor, ycoor,-zcoor,idrotm[4],"ONLY");
756   }
757
758   trpa[0] = 0.5*(kExterInterModBorder2 - kExterInterModBorder1)/TMath::Cos(alpha);
759   trpa[1] = kModuleWallThickness;
760   trpa[2] = xFLT*0.5;
761   trpa[3] = -beta*kRaddeg;
762   trpa[4] = 0.;
763   trpa[5] = 0.;
764   gMC->Gsvolu("FWZ4","PARA", idtmed[503], trpa, 6);    // fibre glass
765
766   AliMatrix (idrotm[13],alpha*kRaddeg,90.,90.+alpha*kRaddeg,90.,90.,180.);
767   AliMatrix (idrotm[16],180.-alpha*kRaddeg,90.,90.-alpha*kRaddeg,90.,90.,0.);
768
769   xcoor = 0.;
770   ycoor = -yFLT*0.25;
771   zcoor = (kExterInterModBorder2 + kExterInterModBorder1)*0.5;
772   gMC->Gspos("FWZ4", 1,"FLTA", xcoor, ycoor, zcoor,idrotm[16],"ONLY");
773   gMC->Gspos("FWZ4", 2,"FLTA", xcoor, ycoor,-zcoor,idrotm[13],"ONLY");
774
775   if (fTOFHoles) {
776     xcoor = 0.;
777     ycoor = -yFLT*0.25;
778     zcoor = -(kExterInterModBorder2 + kExterInterModBorder1)*0.5 + (zlenA*0.5 + kInterCentrModBorder1 - kModuleWallThickness)*0.5;
779     gMC->Gspos("FWZ4", 3,"FLTB", xcoor, ycoor, zcoor,idrotm[13],"ONLY");
780     gMC->Gspos("FWZ4", 4,"FLTC", xcoor, ycoor,-zcoor,idrotm[16],"ONLY");
781   }
782
783
784   ///////////////// Detector itself //////////////////////
785
786   const Int_t    knx   = fTOFGeometry->NpadX();  // number of pads along x
787   const Int_t    knz   = fTOFGeometry->NpadZ();  // number of pads along z
788   const Float_t  kPadX = fTOFGeometry->XPad();   // pad length along x
789   const Float_t  kPadZ = fTOFGeometry->ZPad();   // pad length along z
790
791   // new description for strip volume -double stack strip-
792   // -- all constants are expressed in cm
793   // heigth of different layers
794   const Float_t khhony   = 1.0    ;   // heigth of HONY  Layer
795   const Float_t khpcby   = 0.08   ;   // heigth of PCB   Layer
796   const Float_t khrgly   = 0.055  ;   // heigth of RED GLASS  Layer
797
798   const Float_t khfiliy  = 0.125  ;   // heigth of FISHLINE  Layer
799   const Float_t khglassy = 0.160*0.5; // heigth of GLASS  Layer
800   const Float_t khglfy   = khfiliy+2.*khglassy;// heigth of GLASS+FISHLINE  Layer
801
802   const Float_t khcpcby  = 0.16   ;   // heigth of PCB  Central Layer
803   const Float_t kwhonz   = 8.1    ;   // z dimension of HONEY  Layer
804   const Float_t kwpcbz1  = 10.6   ;   // z dimension of PCB  Lower Layer
805   const Float_t kwpcbz2  = 11.6   ;   // z dimension of PCB  Upper Layer
806   const Float_t kwcpcbz  = 12.4   ;   // z dimension of PCB  Central Layer
807   const Float_t kwrglz   = 8.     ;   // z dimension of RED GLASS  Layer
808   const Float_t kwglfz   = 7.     ;   // z dimension of GLASS+FISHLN Layer
809   const Float_t klsensmx = knx*kPadX; // length of Sensitive Layer
810   const Float_t khsensmy = 0.05;//0.11;//0.16;// heigth of Sensitive Layer
811   const Float_t kwsensmz = knz*kPadZ; // width of Sensitive Layer
812   
813   // heigth of the FSTR Volume (the strip volume)
814   const Float_t khstripy = 2.*khhony+2.*khpcby+4.*khrgly+2.*khglfy+khcpcby;
815
816   // width  of the FSTR Volume (the strip volume)
817   const Float_t kwstripz = kwcpcbz;
818   // length of the FSTR Volume (the strip volume)
819   const Float_t klstripx = fTOFGeometry->StripLength();
820   
821   Float_t parfp[3]={klstripx*0.5,khstripy*0.5,kwstripz*0.5};
822   // Coordinates of the strip center in the strip reference frame;
823   // used for positioning internal strip volumes
824   Float_t posfp[3]={0.,0.,0.};  
825
826   // FSTR volume definition-filling this volume with non sensitive Gas Mixture
827   gMC->Gsvolu("FSTR","BOX",idtmed[507],parfp,3);     // Freon mix
828
829   //-- HONY Layer definition
830   parfp[1] = khhony*0.5;
831   parfp[2] = kwhonz*0.5;
832   gMC->Gsvolu("FHON","BOX",idtmed[501],parfp,3);     // honeycomb (Nomex)
833   // positioning 2 HONY Layers on FSTR volume
834   posfp[1] =-khstripy*0.5+parfp[1];
835   gMC->Gspos("FHON",1,"FSTR",0., posfp[1],0.,0,"ONLY");
836   gMC->Gspos("FHON",2,"FSTR",0.,-posfp[1],0.,0,"ONLY");
837   
838   //-- PCB Layer definition
839   parfp[1] = khpcby*0.5;
840   parfp[2] = kwpcbz1*0.5;
841   gMC->Gsvolu("FPC1","BOX",idtmed[502],parfp,3);     // G10
842   parfp[2] = kwpcbz2*0.5;
843   gMC->Gsvolu("FPC2","BOX",idtmed[502],parfp,3);     // G10
844   // positioning 2 PCB Layers on FSTR volume
845   posfp[1] =-khstripy*0.5+khhony+parfp[1];
846   gMC->Gspos("FPC1",1,"FSTR",0.,-posfp[1],0.,0,"ONLY");
847   gMC->Gspos("FPC2",1,"FSTR",0., posfp[1],0.,0,"ONLY");
848
849   //-- central PCB layer definition
850   parfp[1] = khcpcby*0.5;
851   parfp[2] = kwcpcbz*0.5;
852   gMC->Gsvolu("FPCB","BOX",idtmed[502],parfp,3);     // G10
853   // positioning the central PCB layer
854   gMC->Gspos("FPCB",1,"FSTR",0.,0.,0.,0,"ONLY");
855   
856   //      Sensitive volume
857   Float_t parfs[3] = {klsensmx*0.5, khsensmy*0.5, kwsensmz*0.5};
858   gMC->Gsvolu("FSEN","BOX",idtmed[508],parfs,3);     // sensitive ...
859   // dividing FSEN along z in knz=2 and along x in knx=48
860   gMC->Gsdvn("FSEZ","FSEN",knz,3);
861   gMC->Gsdvn("FPAD","FSEZ",knx,1);
862   // positioning a Sensitive layer inside FPCB
863   gMC->Gspos("FSEN",1,"FPCB",0.,0.,0.,0,"ONLY");
864
865   //-- RED GLASS Layer definition
866   parfp[1] = khrgly*0.5;
867   parfp[2] = kwrglz*0.5;
868   gMC->Gsvolu("FRGL","BOX",idtmed[509],parfp,3);     // glass
869   // positioning 4 RED GLASS Layers on FSTR volume
870   posfp[1] = -khstripy*0.5+khhony+khpcby+parfp[1];
871   gMC->Gspos("FRGL",1,"FSTR",0., posfp[1],0.,0,"ONLY");
872   gMC->Gspos("FRGL",4,"FSTR",0.,-posfp[1],0.,0,"ONLY");
873   posfp[1] = (khcpcby+khrgly)*0.5;
874   gMC->Gspos("FRGL",2,"FSTR",0.,-posfp[1],0.,0,"ONLY");
875   gMC->Gspos("FRGL",3,"FSTR",0., posfp[1],0.,0,"ONLY");
876
877   //-- GLASS Layer definition
878   parfp[1] = khglassy*0.5;
879   parfp[2] = kwglfz*0.5;
880   gMC->Gsvolu("FGLA","BOX",idtmed[509],parfp,3);     // glass
881
882   // positioning 4 GLASS Layers on FSTR volume
883   posfp[1] = -khstripy*0.5+khhony+khpcby+khrgly+parfp[1];
884   gMC->Gspos("FGLA",1,"FSTR",0.,-posfp[1],0.,0,"ONLY");
885   gMC->Gspos("FGLA",4,"FSTR",0., posfp[1],0.,0,"ONLY");
886   posfp[1] = khcpcby*0.5+khrgly+khglassy*0.5;
887   gMC->Gspos("FGLA",2,"FSTR",0.,-posfp[1],0.,0,"ONLY");
888   gMC->Gspos("FGLA",3,"FSTR",0., posfp[1],0.,0,"ONLY");
889
890   //-- FREON Layer definition
891   parfp[1] = khfiliy*0.5;
892   gMC->Gsvolu("FFIS","BOX",idtmed[507],parfp,3);     // freon
893
894   // positioning 2 FREON Layers on FSTR volume
895   posfp[1] = -khstripy*0.5+khhony+khpcby+khrgly+khglassy+parfp[1];
896   gMC->Gspos("FFIS",1,"FSTR",0.,-posfp[1],0.,0,"ONLY");
897   gMC->Gspos("FFIS",2,"FSTR",0., posfp[1],0.,0,"ONLY");
898
899   /*
900   //-- GLASS+FISHLINE Layer definition
901   parfp[1] = khglfy*0.5;
902   parfp[2] = kwglfz*0.5;
903   gMC->Gsvolu("FGLF","BOX",idtmed[504],parfp,3);
904
905   // positioning 2 GLASS+FISHLINE Layers on FSTR volume
906   posfp[1] = (khcpcby+khglfy)*0.5+khrgly;
907   gMC->Gspos("FGLF",1,"FSTR",0.,-posfp[1],0.,0,"ONLY");
908   gMC->Gspos("FGLF",2,"FSTR",0., posfp[1],0.,0,"ONLY");
909   */
910
911   //  Positioning the Strips  (FSTR) in the FLT volumes
912   Int_t maxStripNumbers [5] ={fTOFGeometry->NStripC(),
913                               fTOFGeometry->NStripB(),
914                               fTOFGeometry->NStripA(),
915                               fTOFGeometry->NStripB(),
916                               fTOFGeometry->NStripC()};
917
918   Int_t totalStrip = 0;
919   Float_t xpos, zpos, ypos, ang;
920   for(Int_t iplate =0; iplate < fTOFGeometry->NPlates(); iplate++){
921     if (iplate>0) totalStrip += maxStripNumbers[iplate-1];
922     for(Int_t istrip =0; istrip < maxStripNumbers[iplate]; istrip++){
923
924       ang = fTOFGeometry->GetAngles(iplate,istrip);
925       AliDebug(1, Form(" iplate = %1i, istrip = %2i ---> ang = %f", iplate, istrip, ang));
926  
927       if (ang>0.)       AliMatrix (idrotm[istrip+totalStrip+1],90.,0.,90.+ang,90., ang, 90.);
928       else if (ang==0.) AliMatrix (idrotm[istrip+totalStrip+1],90.,0.,90.,90., 0., 0.);
929       else if (ang<0.)  AliMatrix (idrotm[istrip+totalStrip+1],90.,0.,90.+ang,90.,-ang,270.);
930
931       xpos = 0.;
932       zpos = fTOFGeometry->GetDistances(iplate,istrip);
933       ypos = fTOFGeometry->GetHeights(iplate,istrip);
934
935       gMC->Gspos("FSTR",istrip+totalStrip+1,"FLTA", xpos, ypos,-zpos,idrotm[istrip+totalStrip+1],  "ONLY");
936
937       if (fTOFHoles) {
938         if (istrip+totalStrip+1>53) gMC->Gspos("FSTR",istrip+totalStrip+1,"FLTC", xpos, ypos,-zpos-(zlenA*0.5 + kInterCentrModBorder1 - kModuleWallThickness)*0.5,idrotm[istrip+totalStrip+1],"ONLY");
939         if (istrip+totalStrip+1<39) gMC->Gspos("FSTR",istrip+totalStrip+1,"FLTB", xpos, ypos,-zpos+(zlenA*0.5 + kInterCentrModBorder1 - kModuleWallThickness)*0.5,idrotm[istrip+totalStrip+1],"ONLY");
940       }
941     }
942   }
943
944   //  1.5 cm Al honeycomb layer between strips and cards
945   par[0] = xFLT*0.5;
946   par[1] = kHoneycombLayerThickness*0.5;
947   par[2] = zFLTA*0.5;
948   gMC->Gsvolu("FPEA", "BOX ", idtmed[506], par, 3);   // Al honeycomb
949
950   xcoor = 0.;
951   ycoor = kHoneycombLayerThickness*0.5;
952   zcoor = 0.;
953   gMC->Gspos ("FPEA", 0, "FLTA", xcoor, ycoor, zcoor, 0, "ONLY");
954
955   if (fTOFHoles) {
956     par[0] = xFLT*0.5;
957     par[1] = kHoneycombLayerThickness*0.5;
958     par[2] = (zlenA*0.5 - kInterCentrModBorder2-kModuleWallThickness)*0.5;
959     gMC->Gsvolu("FPEB", "BOX ", idtmed[506], par, 3);   // Al honeycomb
960
961     xcoor = 0.;
962     ycoor = kHoneycombLayerThickness*0.5;
963     zcoor = (kInterCentrModBorder2-kInterCentrModBorder1)*0.5;
964     gMC->Gspos ("FPEB", 1, "FLTB", xcoor, ycoor,-zcoor, 0, "ONLY");
965     gMC->Gspos ("FPEB", 2, "FLTC", xcoor, ycoor, zcoor, 0, "ONLY");
966   }
967
968   // frame of Air
969   par[0] = xFLT*0.5;
970   par[1] = (yFLT*0.5 - kHoneycombLayerThickness)*0.5;
971   par[2] = zFLTA *0.5;
972   gMC->Gsvolu("FAIA", "BOX ", idtmed[500], par, 3); // Air
973
974   xcoor = 0.;
975   ycoor = kHoneycombLayerThickness + (yFLT*0.5 - kHoneycombLayerThickness)*0.5;
976   zcoor = 0.;
977   gMC->Gspos ("FAIA", 0, "FLTA", xcoor, ycoor, zcoor, 0, "ONLY");
978
979   if (fTOFHoles) {
980     par[0] = xFLT*0.5;
981     par[1] = (yFLT*0.5 - kHoneycombLayerThickness)*0.5;
982     par[2] = (zlenA*0.5 - kInterCentrModBorder2 - kModuleWallThickness)*0.5;
983     gMC->Gsvolu("FAIB", "BOX ", idtmed[500], par, 3); // Air
984     gMC->Gsvolu("FAIC", "BOX ", idtmed[500], par, 3); // Air
985
986     xcoor = 0.;
987     ycoor = kHoneycombLayerThickness + (yFLT*0.5 - kHoneycombLayerThickness)*0.5;
988     zcoor = (kInterCentrModBorder2-kInterCentrModBorder1)*0.5;
989     gMC->Gspos ("FAIB", 0, "FLTB", xcoor, ycoor,-zcoor, 0, "ONLY");
990     gMC->Gspos ("FAIC", 0, "FLTC", xcoor, ycoor, zcoor, 0, "ONLY");
991   }
992
993   // start with cards and cooling tubes
994   // finally, cards, cooling tubes and layer for thermal dispersion
995   // 3 volumes
996   
997   // see GEOM200 in GEANT manual
998
999   Float_t cardpar[3];
1000
1001   // card volume definition
1002   cardpar[0]= xFLT*0.5;
1003   cardpar[1]= 5.;
1004   cardpar[2]= 0.1;
1005   gMC->Gsvolu("FCAR", "BOX ", idtmed[502], cardpar, 3); // PCB Card 
1006
1007   //alu plate volume definition
1008   cardpar[1]= 3.5;
1009   cardpar[2]= 0.05;
1010   gMC->Gsvolu("FALP", "BOX ", idtmed[505], cardpar, 3); // Alu Plate
1011
1012   // tube volume definition
1013   Float_t tubepar[3];
1014   tubepar[0]= 0.;
1015   tubepar[1]= 0.4;
1016   tubepar[2]= 61.;
1017   gMC->Gsvolu("FTUB", "TUBE", idtmed[511], tubepar, 3); // cooling tubes (steel)
1018
1019   //tubepar[0]= 0.;
1020   tubepar[1]= 0.35;
1021   //tubepar[2]= 61.;
1022   gMC->Gsvolu("FITU", "TUBE", idtmed[510], tubepar, 3); // cooling water
1023   // positioning water tube into the steel one
1024   gMC->Gspos("FITU",1,"FTUB",0.,0.,0.,0,"ONLY");
1025
1026   // rotation matrix
1027   AliMatrix(idrotm[99], 180., 90., 90., 90., 90., 0.);
1028
1029   // central module positioning
1030   Float_t cardpos[3], aplpos2;
1031   Float_t stepforcardA = 6.625;
1032   Float_t tdis = 0.6;
1033   Float_t aplpos1 = -2.;
1034
1035   cardpos[0]= 0.;
1036   cardpos[1]= -0.5;
1037   cardpos[2]= -53.;
1038   //  tubepos= -53.+tdis;
1039   Int_t icard;
1040   for (icard=39; icard<54; ++icard) {
1041     cardpos[2]= cardpos[2]+stepforcardA;
1042     aplpos2 = cardpos[2]+0.15;
1043     gMC->Gspos("FCAR",icard,"FAIA",cardpos[0],cardpos[1],     cardpos[2],         0,"ONLY"); 
1044     gMC->Gspos("FALP",icard,"FAIA",cardpos[0],   aplpos1,        aplpos2,         0,"ONLY");
1045     gMC->Gspos("FTUB",icard,"FAIA",        0.,cardpos[1],cardpos[2]+tdis,idrotm[99],"ONLY");
1046   }
1047   
1048   // intermediate module positioning
1049   Float_t stepforcardB= 7.05;
1050   Float_t offs = 53.;
1051
1052   cardpos[2]= offs;
1053   for (icard=20; icard<39; ++icard) {
1054     cardpos[2]= cardpos[2]+stepforcardB;
1055     aplpos2 = cardpos[2]+0.15;
1056
1057     gMC->Gspos("FCAR",icard+34,"FAIA",cardpos[0],cardpos[1],      cardpos[2],         0,"ONLY");
1058     gMC->Gspos("FALP",icard+34,"FAIA",cardpos[0],   aplpos1,         aplpos2,         0,"ONLY");
1059     gMC->Gspos("FTUB",icard+34,"FAIA",        0.,cardpos[1], cardpos[2]+tdis,idrotm[99],"ONLY");
1060     gMC->Gspos("FCAR",58-icard,"FAIA",cardpos[0],cardpos[1],     -cardpos[2],         0,"ONLY");
1061     gMC->Gspos("FALP",58-icard,"FAIA",cardpos[0],   aplpos1,        -aplpos2,         0,"ONLY");
1062     gMC->Gspos("FTUB",58-icard,"FAIA",        0.,cardpos[1],-cardpos[2]-tdis,idrotm[99],"ONLY");
1063
1064     if (fTOFHoles) {
1065       gMC->Gspos("FCAR",icard+34+182,"FAIC",cardpos[0],cardpos[1],      cardpos[2]-(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1066       gMC->Gspos("FALP",icard+34+182,"FAIC",cardpos[0],   aplpos1,         aplpos2-(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1067       gMC->Gspos("FTUB",icard+34+182,"FAIC",        0.,cardpos[1], cardpos[2]+tdis-(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,idrotm[99],"ONLY");
1068       gMC->Gspos("FCAR",58-icard+ 91,"FAIB",cardpos[0],cardpos[1],     -cardpos[2]+(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1069       gMC->Gspos("FALP",58-icard+ 91,"FAIB",cardpos[0],   aplpos1,        -aplpos2+(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1070       gMC->Gspos("FTUB",58-icard+ 91,"FAIB",        0.,cardpos[1],-cardpos[2]-tdis+(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,idrotm[99],"ONLY");
1071     }
1072     
1073   }
1074
1075   // outer module positioning
1076   Float_t stepforcardC= 8.45238;
1077   offs += zlenB;
1078   cardpos[2]= offs;
1079   for (icard=1; icard<20; ++icard) {
1080     cardpos[2]= cardpos[2]+stepforcardC;
1081     aplpos2 = cardpos[2]+0.15;
1082
1083     gMC->Gspos("FCAR",icard+72,"FAIA",cardpos[0],cardpos[1],      cardpos[2],         0,"ONLY"); 
1084     gMC->Gspos("FALP",icard+72,"FAIA",cardpos[0],   aplpos1,         aplpos2,         0,"ONLY");
1085     gMC->Gspos("FTUB",icard+72,"FAIA",        0.,cardpos[1], cardpos[2]+tdis,idrotm[99],"ONLY");
1086     gMC->Gspos("FCAR",20-icard,"FAIA",cardpos[0],cardpos[1],     -cardpos[2],         0,"ONLY");
1087     gMC->Gspos("FALP",20-icard,"FAIA",cardpos[0],   aplpos1,        -aplpos2,         0,"ONLY");
1088     gMC->Gspos("FTUB",20-icard,"FAIA",        0.,cardpos[1],-cardpos[2]-tdis,idrotm[99],"ONLY");
1089
1090     if (fTOFHoles) {
1091       gMC->Gspos("FCAR",icard+72+182,"FAIC",cardpos[0],cardpos[1],      cardpos[2]-(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1092       gMC->Gspos("FALP",icard+72+182,"FAIC",cardpos[0],   aplpos1,         aplpos2-(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1093       gMC->Gspos("FTUB",icard+72+182,"FAIC",        0.,cardpos[1], cardpos[2]+tdis-(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,idrotm[99],"ONLY");
1094       gMC->Gspos("FCAR",20-icard+ 91,"FAIB",cardpos[0],cardpos[1],     -cardpos[2]+(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1095       gMC->Gspos("FALP",20-icard+ 91,"FAIB",cardpos[0],   aplpos1,        -aplpos2+(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,         0,"ONLY");
1096       gMC->Gspos("FTUB",20-icard+ 91,"FAIB",        0.,cardpos[1],-cardpos[2]-tdis+(zlenA*0.5 + kInterCentrModBorder2 - kModuleWallThickness)*0.5,idrotm[99],"ONLY");
1097     }
1098   }
1099
1100 }
1101 //_____________________________________________________________________________
1102 void AliTOFv5T0::DrawModule() const
1103 {
1104   //
1105   // Draw a shaded view of the Time Of Flight version 4
1106   //
1107
1108   // Set everything unseen
1109   gMC->Gsatt("*", "seen", -1);
1110
1111   //
1112   //Set volumes visible
1113   // 
1114
1115   //Set ALIC mother transparent
1116   gMC->Gsatt("ALIC","SEEN", 0);
1117
1118 //=====> Level 1
1119   // Level 1 for TOF volumes
1120   gMC->Gsatt("B077","seen", 0);
1121
1122 //=====> Level 2
1123   // Level 2 for TOF volumes
1124   gMC->Gsatt("B071","seen", 0);
1125   gMC->Gsatt("B074","seen", 0);
1126   gMC->Gsatt("B075","seen", 0);
1127   gMC->Gsatt("B076","seen",-1); // all B076 sub-levels skipped -
1128   gMC->Gsatt("B080","seen", 0);  // B080 does not has sub-level                
1129
1130   // Level 2 of B071
1131   gMC->Gsatt("B056","seen", 0);  // B056 does not has sub-levels  -
1132   gMC->Gsatt("B063","seen",-1); // all B063 sub-levels skipped   -
1133   gMC->Gsatt("B065","seen",-1); // all B065 sub-levels skipped   -
1134   gMC->Gsatt("B067","seen",-1); // all B067 sub-levels skipped   -
1135   gMC->Gsatt("B072","seen",-1); // all B072 sub-levels skipped   -
1136
1137   gMC->Gsatt("BTR1","seen", 0);  // all BTR1 sub-levels skipped   -
1138   gMC->Gsatt("BTO1","seen", 0);
1139
1140   // Level 2 of B074
1141   gMC->Gsatt("BTR2","seen", 0);  // all BTR1 sub-levels skipped   -
1142   gMC->Gsatt("BTO2","seen", 0);
1143
1144   // Level 2 of B075
1145   gMC->Gsatt("BTR3","seen", 0);  // all BTR1 sub-levels skipped   -
1146   gMC->Gsatt("BTO3","seen", 0);
1147
1148   // Level 3 of B071, B074 and B075
1149   gMC->Gsatt("FTOA","SEEN", 0);
1150   if (fTOFHoles) gMC->Gsatt("FTOB","SEEN", 0);
1151
1152   // Level 4 of B071, B074 and B075
1153   gMC->Gsatt("FLTA","SEEN", 0);
1154   if (fTOFHoles) gMC->Gsatt("FLTB","SEEN",0);
1155   if (fTOFHoles) gMC->Gsatt("FLTC","SEEN",0);
1156
1157   // Level 5 of B071, B074 and B075
1158   gMC->Gsatt("FAIA","SEEN",-1);  // all FAIA sub-levels skipped   -
1159   if (fTOFHoles) gMC->Gsatt("FAIB","SEEN",-1);  // all FAIB sub-levels skipped   -
1160   if (fTOFHoles) gMC->Gsatt("FAIC","SEEN",-1);  // all FAIC sub-levels skipped   -
1161
1162   gMC->Gsatt("FALA","SEEN", 0);
1163   if (fTOFHoles) gMC->Gsatt("FALB","SEEN", 0);
1164
1165   gMC->Gsatt("FPEA","SEEN", 1);
1166   if (fTOFHoles) gMC->Gsatt("FPEB","SEEN", 1);
1167
1168   gMC->Gsatt("FSTR","SEEN",-1);  // all FSTR sub-levels skipped   -
1169
1170   gMC->Gsatt("FWZ1","SEEN", 0);
1171   gMC->Gsatt("FWZ2","SEEN", 0);
1172   gMC->Gsatt("FWZ3","SEEN", 0);
1173   gMC->Gsatt("FWZ4","SEEN", 0);
1174   if (fTOFHoles) {
1175     gMC->Gsatt("FZ1B","SEEN", 0);
1176     gMC->Gsatt("FZ1C","SEEN", 0);
1177     gMC->Gsatt("FZ2B","SEEN", 0);
1178   }
1179
1180   gMC->Gdopt("hide", "on");
1181   gMC->Gdopt("shad", "on");
1182   gMC->Gsatt("*", "fill", 7);
1183   gMC->SetClipBox(".");
1184   gMC->SetClipBox("*", 100, 1000, 100, 1000, 100, 1000);
1185   gMC->DefaultRange();
1186   gMC->Gdraw("alic", 40, 30, 0, 10, 9.5, .018, .018);
1187   gMC->Gdhead(1111, "Time Of Flight");
1188   gMC->Gdman(18, 3, "MAN");
1189   gMC->Gdopt("hide","off");
1190 }
1191 //_____________________________________________________________________________
1192 void AliTOFv5T0::DrawDetectorModules() const
1193 {
1194   //
1195   // Draw a shaded view of the TOF detector version 4
1196   //
1197  
1198   // Set everything unseen
1199   gMC->Gsatt("*", "seen", -1);
1200
1201   //
1202   //Set volumes visible
1203   // 
1204
1205   //Set ALIC mother transparent
1206   gMC->Gsatt("ALIC","SEEN", 0);
1207
1208 //=====> Level 1
1209   // Level 1 for TOF volumes
1210   gMC->Gsatt("B077","seen", 0);
1211
1212 //=====> Level 2
1213   // Level 2 for TOF volumes
1214   gMC->Gsatt("B071","seen", 0);
1215   gMC->Gsatt("B074","seen", 0);
1216   gMC->Gsatt("B075","seen", 0);
1217   gMC->Gsatt("B076","seen",-1); // all B076 sub-levels skipped -
1218   gMC->Gsatt("B080","seen", 0);  // B080 does not has sub-level                
1219
1220   // Level 2 of B071
1221   gMC->Gsatt("B056","seen", 0);  // B056 does not has sub-levels  -
1222   gMC->Gsatt("B063","seen",-1); // all B063 sub-levels skipped   -
1223   gMC->Gsatt("B065","seen",-1); // all B065 sub-levels skipped   -
1224   gMC->Gsatt("B067","seen",-1); // all B067 sub-levels skipped   -
1225   gMC->Gsatt("B072","seen",-1); // all B072 sub-levels skipped   -
1226
1227   gMC->Gsatt("BTR1","seen", 0);  // all BTR1 sub-levels skipped   -
1228   gMC->Gsatt("BTO1","seen", 0);
1229
1230   // Level 2 of B074
1231   gMC->Gsatt("BTR2","seen", 0);  // all BTR1 sub-levels skipped   -
1232   gMC->Gsatt("BTO2","seen", 0);
1233
1234   // Level 2 of B075
1235   gMC->Gsatt("BTR3","seen", 0);  // all BTR1 sub-levels skipped   -
1236   gMC->Gsatt("BTO3","seen", 0);
1237
1238   // Level 3 of B071, B075 and B074
1239   gMC->Gsatt("FTOA","seen",-2);  // all FTOA sub-levels skipped   -
1240   if (fTOFHoles) {
1241     gMC->Gsatt("FTOB","seen",-2);  // all FTOB sub-levels skipped   -
1242     gMC->Gsatt("FTOC","seen",-2);  // all FTOC sub-levels skipped   -
1243   }
1244
1245   gMC->Gdopt("hide","on");
1246   gMC->Gdopt("shad","on");
1247   gMC->Gsatt("*", "fill", 5);
1248   gMC->SetClipBox(".");
1249   gMC->SetClipBox("*", 100, 1000, 100, 1000, 0, 1000);
1250   gMC->DefaultRange();
1251   gMC->Gdraw("alic", 40, 30, 0, 10, 9.5, .018, .018);
1252   gMC->Gdhead(1111,"TOF detector");
1253   gMC->Gdman(18, 3, "MAN");
1254   gMC->Gdopt("hide","off");
1255 }                                 
1256
1257 //_____________________________________________________________________________
1258 void AliTOFv5T0::DrawDetectorStrips() const
1259 {
1260   //
1261   // Draw a shaded view of the TOF strips for version 4
1262   //
1263
1264   // Set everything unseen
1265   gMC->Gsatt("*", "seen", -1);
1266
1267   //
1268   //Set volumes visible
1269   // 
1270   
1271   //Set ALIC mother transparent
1272   gMC->Gsatt("ALIC","SEEN", 0);
1273   
1274 //=====> Level 1
1275   // Level 1 for TOF volumes
1276   gMC->Gsatt("B077","seen", 0);
1277
1278 //=====> Level 2
1279   // Level 2 for TOF volumes
1280   gMC->Gsatt("B071","seen", 0);
1281   gMC->Gsatt("B074","seen", 0);
1282   gMC->Gsatt("B075","seen", 0);
1283   gMC->Gsatt("B076","seen",-1); // all B076 sub-levels skipped -
1284   gMC->Gsatt("B080","seen", 0);  // B080 does not has sub-level                
1285
1286   // Level 2 of B071
1287   gMC->Gsatt("B063","seen",-1); // all B063 sub-levels skipped   -
1288   gMC->Gsatt("B065","seen",-1); // all B065 sub-levels skipped   -
1289   gMC->Gsatt("B067","seen",-1); // all B067 sub-levels skipped   -
1290   gMC->Gsatt("B056","seen", 0);  // B056 does not has sub-levels  -
1291   gMC->Gsatt("B072","seen",-1); // all B072 sub-levels skipped   -
1292
1293   gMC->Gsatt("BTR1","seen", 0);  // all BTR1 sub-levels skipped   -
1294   gMC->Gsatt("BTO1","seen", 0);
1295
1296   // Level 2 of B074
1297   gMC->Gsatt("BTR2","seen", 0);  // all BTR1 sub-levels skipped   -
1298   gMC->Gsatt("BTO2","seen", 0);
1299
1300   // Level 2 of B075
1301   gMC->Gsatt("BTR3","seen", 0);  // all BTR1 sub-levels skipped   -
1302   gMC->Gsatt("BTO3","seen", 0);
1303
1304   // Level 3 of B071, B074 and B075
1305   gMC->Gsatt("FTOA","SEEN", 0);
1306   if (fTOFHoles) {
1307     gMC->Gsatt("FTOB","SEEN", 0);
1308     gMC->Gsatt("FTOC","SEEN", 0);
1309   }
1310
1311   // Level 4 of B071, B074 and B075
1312   gMC->Gsatt("FLTA","SEEN", 0);
1313   if (fTOFHoles) {
1314     gMC->Gsatt("FLTB","SEEN", 0);
1315     gMC->Gsatt("FLTC","SEEN", 0);
1316   }
1317
1318   // Level 5 of B071, B074 and B075
1319   gMC->Gsatt("FAIA","SEEN",-1);  // all FAIA sub-levels skipped   -
1320   if (fTOFHoles) {
1321     gMC->Gsatt("FAIB","SEEN",-1);  // all FAIB sub-levels skipped   -
1322     gMC->Gsatt("FAIC","SEEN",-1);  // all FAIC sub-levels skipped   -
1323   }
1324
1325   gMC->Gsatt("FALA","SEEN", 0);
1326   if (fTOFHoles) gMC->Gsatt("FALB","SEEN", 0);
1327
1328   gMC->Gsatt("FPEA","SEEN", 0);
1329   if (fTOFHoles) gMC->Gsatt("FPEB","SEEN", 0);
1330
1331   gMC->Gsatt("FSTR","SEEN",-2);  // all FSTR sub-levels skipped   -
1332
1333   gMC->Gsatt("FWZ1","SEEN", 0);
1334   gMC->Gsatt("FWZ2","SEEN", 0);
1335   gMC->Gsatt("FWZ3","SEEN", 0);
1336   gMC->Gsatt("FWZ4","SEEN", 0);
1337   if (fTOFHoles){
1338     gMC->Gsatt("FZ1B","SEEN", 0);
1339     gMC->Gsatt("FZ1C","SEEN", 0);
1340     gMC->Gsatt("FZ2B","SEEN", 0);
1341   }
1342
1343   /*
1344   // Level 2 of FAIA
1345   // Level 2 of FAIB
1346   // Level 2 of FAIC
1347   gMC->Gsatt("FALP","SEEN",0);
1348   gMC->Gsatt("FCAR","SEEN",0);
1349   gMC->Gsatt("FTUB","SEEN",-1);  // all FTUB sub-levels skipped   -
1350
1351   // Level 2 of FTUB
1352   gMC->Gsatt("FITU","SEEN",0);
1353   */
1354
1355   /*
1356   // Level 2 of FSTR
1357   gMC->Gsatt("FGLF","SEEN",0);
1358   gMC->Gsatt("FHON","SEEN",0);
1359   gMC->Gsatt("FPC1","SEEN",0);
1360   gMC->Gsatt("FPC2","SEEN",0);
1361   gMC->Gsatt("FPCB","SEEN",0);
1362   gMC->Gsatt("FRGL","SEEN",0);
1363
1364   // Level 2 of FPCB => Level 3 of FSTR
1365   gMC->Gsatt("FSEN","SEEN",0);
1366   gMC->Gsatt("FSEZ","SEEN",0);
1367   gMC->Gsatt("FPAD","SEEN",1);
1368   */
1369
1370   gMC->Gdopt("hide","on");
1371   gMC->Gdopt("shad","on");
1372   gMC->Gsatt("*", "fill", 5);
1373   gMC->SetClipBox(".");
1374   gMC->SetClipBox("*", 0, 1000, 0, 1000, 0, 1000);
1375   gMC->DefaultRange();
1376   gMC->Gdraw("alic", 40, 30, 0, 10, 9.5, .018, .018);
1377   gMC->Gdhead(1111,"TOF Strips");
1378   gMC->Gdman(18, 3, "MAN");
1379   gMC->Gdopt("hide","off");
1380 }
1381
1382 //_____________________________________________________________________________
1383 void AliTOFv5T0::CreateMaterials()
1384 {
1385   //
1386   // Define materials for the Time Of Flight
1387   //
1388
1389   //AliTOF::CreateMaterials();
1390
1391   AliMagF *magneticField = (AliMagF*)gAlice->Field();
1392
1393   Int_t   isxfld = magneticField->Integ();
1394   Float_t sxmgmx = magneticField->Max();
1395
1396   Float_t we[7], ae[7], na[7], fr[7], vl[7];
1397   Int_t i;
1398
1399   //--- Quartz (SiO2) to simulate float glass
1400   //    density tuned to have correct float glass 
1401   //    radiation length
1402   Float_t   aq[2] = { 28.0855,15.9994 };
1403   Float_t   zq[2] = { 14.,8. };
1404   Float_t   wq[2] = { 1.,2. };
1405   Float_t   dq = 2.55; // std value: 2.2
1406   Int_t nq = -2;
1407
1408   // --- Nomex
1409   Float_t anox[4] = {12.01,1.01,16.00,14.01};
1410   Float_t znox[4] = { 6.,  1.,  8.,  7.};
1411   Float_t wnox[4] = {14., 22., 2., 2.};
1412   Float_t dnox  = 0.048;
1413   Int_t nnox   = -4;
1414
1415   //  { Si, C, H, O }
1416   Float_t ag10[4] = {28.09,12.01,1.01,16.00};
1417   Float_t zg10[4] = {14.,  6.,  1.,  8.};
1418   Float_t wmatg10[4];
1419   Int_t nlmatg10  = 4;
1420   for (i = 0; i < nlmatg10; ++i) {
1421     ae[i] = ag10[i];
1422     vl[i] = 1.;
1423   }
1424   ae[4] = 16.00;
1425   vl[4] = 1.;
1426   na[0] = 1.; 
1427   na[1] = 14.;
1428   na[2] = 20.;
1429   na[3] = 2.;
1430   na[4] = 3.;
1431   fr[0] = 0.6;
1432   fr[1] = 0.4;
1433   fr[2] = 0.4;
1434   fr[3] = 0.6;
1435   fr[4] = 0.4;
1436   MaterialMixer(we,ae,na,fr,vl,5);
1437   we[3] += we[4];
1438   wmatg10[0]= we[0];
1439   wmatg10[1]= we[1];
1440   wmatg10[2]= we[2];
1441   wmatg10[3]= we[3];
1442   Float_t densg10 = 1.7;
1443
1444   // -- Water
1445   Float_t awa[2] = {  1., 16. };
1446   Float_t zwa[2] = {  1.,  8. };
1447   Float_t wwa[2] = {  2.,  1. };
1448   Float_t dwa    = 1.0;
1449   Int_t nwa = -2;
1450
1451   // stainless steel
1452   Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
1453   Float_t zsteel[4] = { 26.,24.,28.,14. };
1454   Float_t wsteel[4] = { .715,.18,.1,.005 };
1455
1456   // AIR
1457   Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
1458   Float_t zAir[4]={6.,7.,8.,18.};
1459   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
1460   Float_t dAir = 1.20479E-3;
1461
1462   // --- fibre glass
1463   Float_t afg[4] = {28.09,16.00,12.01,1.01};
1464   Float_t zfg[4] = {14., 8., 6., 1.};
1465   Float_t wfg[4] = {0.12906,0.29405,0.51502,0.06187};
1466   Float_t dfg  = 1.111;
1467   Int_t nfg   = 4;
1468
1469   // --- Freon C2F4H2 + SF6
1470   Float_t afre[4]= {12.01,1.01,19.00,32.07};
1471   Float_t zfre[4]= { 6., 1., 9., 16.};
1472   Float_t wfre[4]= {0.21250,0.01787,0.74827,0.021355};
1473   Float_t densfre= 0.00375;
1474   Int_t nfre  = 4;
1475
1476   //char namat[15] = "            ";
1477   //Float_t ama[2], zma[2], dma, radl, absl, buf[1];
1478   //Int_t nbuf;
1479
1480   AliMixture ( 0, "Air$", aAir, zAir, dAir, 4, wAir);
1481   AliMixture ( 1, "Nomex$", anox, znox, dnox, nnox, wnox);
1482   AliMixture ( 2, "G10$", ag10, zg10, densg10, nlmatg10, wmatg10);
1483   AliMixture ( 3, "fibre glass$", afg, zfg, dfg, nfg, wfg);
1484   AliMaterial( 4, "Al $", 26.98, 13., 2.7, 8.9, 37.2);
1485   AliMaterial( 5, "Al honeycomb$", 26.98, 13., 0.0496, 483., 2483.);
1486   AliMixture ( 6, "Freon$",  afre, zfre, densfre, nfre, wfre);
1487   AliMixture ( 7, "Glass$", aq, zq, dq, nq, wq);
1488   /*
1489   // get freon and glass
1490   gMC->Gfmate((*fIdmate)[6],namat,ama[0],zma[0],dma,radl,absl,buf,nbuf);
1491   gMC->Gfmate((*fIdmate)[7],namat,ama[1],zma[1],dma,radl,absl,buf,nbuf);
1492
1493   // --- glass-freon
1494   Float_t wgfr[2]= {0.0011,0.9989};
1495   Float_t dgfr = 1.434;
1496   Int_t ngfr  = 2;
1497   AliMixture ( 8, "glass-freon$", ama, zma, dgfr, ngfr, wgfr);
1498   */
1499   AliMixture ( 9, "Water$",  awa, zwa, dwa, nwa, wwa);
1500   AliMixture (10, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
1501
1502   Float_t epsil, stmin, deemax, stemax;
1503  
1504   //   Previous data
1505   //       EPSIL  = 0.1   ! Tracking precision,
1506   //       STEMAX = 0.1   ! Maximum displacement for multiple scattering
1507   //       DEEMAX = 0.1   ! Maximum fractional energy loss, DLS
1508   //       STMIN  = 0.1
1509
1510   //   New data
1511   epsil  = .001;  // Tracking precision,
1512   stemax = -1.;   // Maximum displacement for multiple scattering
1513   deemax = -.3;   // Maximum fractional energy loss, DLS
1514   stmin  = -.8;
1515
1516   AliMedium( 1, "Air$",         0, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1517   AliMedium( 2,"Nomex$",        1, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1518   AliMedium( 3,"G10$",          2, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1519   AliMedium( 4,"fibre glass$",  3, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1520   //AliMedium( 5,"glass-freon$",  8, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1521   AliMedium( 6,"Al Frame$",     4, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1522   AliMedium( 7,"Al honeycomb$", 5, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1523   AliMedium( 8,"Fre$",          6, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1524   AliMedium( 9,"PCB-S$",        2, 1, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1525   AliMedium(10,"Glass$",        7, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1526   AliMedium(11,"Water$",        9, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1527   AliMedium(12,"STEEL$",       10, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
1528
1529 }
1530 //_____________________________________________________________________________
1531 void AliTOFv5T0::Init()
1532 {
1533   //
1534   // Initialise the detector after the geometry has been defined
1535   //
1536   AliDebug(1, "**************************************"
1537            "  TOF  "
1538            "**************************************");
1539   AliDebug(1, "  Version 4 of TOF initialing, "
1540            "symmetric TOF - Full Coverage version");
1541   
1542   AliTOF::Init();
1543   
1544   fIdFTOA = gMC->VolId("FTOA");
1545   if (fTOFHoles) {
1546     fIdFTOB = gMC->VolId("FTOB");
1547     fIdFTOC = gMC->VolId("FTOC");
1548   }
1549   fIdFLTA = gMC->VolId("FLTA");
1550   if (fTOFHoles) {
1551     fIdFLTB = gMC->VolId("FLTB");
1552     fIdFLTC = gMC->VolId("FLTC");
1553   }
1554
1555   AliDebug(1, "**************************************"
1556            "  TOF  "
1557            "**************************************");
1558 }
1559  
1560 //_____________________________________________________________________________
1561 void AliTOFv5T0::StepManager()
1562 {
1563
1564   //
1565   // Procedure called at each step in the Time Of Flight
1566   //
1567
1568   TLorentzVector mom, pos;
1569   Float_t xm[3],pm[3],xpad[3],ppad[3];
1570   Float_t hits[14];
1571   Int_t   vol[5];
1572   Int_t   sector, plate, padx, padz, strip;
1573   Int_t   copy, padzid, padxid, stripid, i;
1574   Int_t   *idtmed = fIdtmed->GetArray()-499;
1575   Float_t incidenceAngle;
1576
1577   const char* volpath;
1578
1579   Int_t index = 0;
1580
1581   if(
1582      gMC->IsTrackEntering()
1583      && gMC->TrackCharge()
1584      //&& gMC->GetMedium()==idtmed[508]
1585      && gMC->CurrentMedium()==idtmed[508]
1586      && gMC->CurrentVolID(copy)==fIdSens
1587      )
1588   {
1589
1590     AliMC *mcApplication = (AliMC*)gAlice->GetMCApp();
1591
1592     AddTrackReference(mcApplication->GetCurrentTrackNumber(), AliTrackReference::kTOF);
1593     //AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber());
1594
1595     // getting information about hit volumes
1596     
1597     padzid=gMC->CurrentVolOffID(1,copy);
1598     padz=copy;
1599     padz--;
1600
1601     padxid=gMC->CurrentVolOffID(0,copy);
1602     padx=copy; 
1603     padx--;
1604     
1605     stripid=gMC->CurrentVolOffID(4,copy);
1606     strip=copy; 
1607     strip--;
1608
1609     gMC->TrackPosition(pos);
1610     gMC->TrackMomentum(mom);
1611
1612     Double_t normMom=1./mom.Rho();
1613
1614     //  getting the coordinates in pad ref system
1615
1616     xm[0] = (Float_t)pos.X();
1617     xm[1] = (Float_t)pos.Y();
1618     xm[2] = (Float_t)pos.Z();
1619
1620     pm[0] = (Float_t)mom.X()*normMom;
1621     pm[1] = (Float_t)mom.Y()*normMom;
1622     pm[2] = (Float_t)mom.Z()*normMom;
1623  
1624     gMC->Gmtod(xm,xpad,1); // from MRS to DRS: coordinates convertion
1625     gMC->Gmtod(pm,ppad,2); // from MRS to DRS: direction cosinus convertion
1626
1627
1628     if (TMath::Abs(ppad[1])>1) {
1629       AliWarning("Abs(ppad) > 1");
1630       ppad[1]=TMath::Sign((Float_t)1,ppad[1]);
1631     }
1632     incidenceAngle = TMath::ACos(ppad[1])*kRaddeg;
1633
1634     plate = -1;
1635     if      (strip <  fTOFGeometry->NStripC()) {
1636       plate = 0;
1637       //strip = strip;
1638     }
1639     else if (strip >= fTOFGeometry->NStripC() && 
1640              strip <  fTOFGeometry->NStripC() + fTOFGeometry->NStripB()) {
1641       plate = 1;
1642       strip = strip - fTOFGeometry->NStripC();
1643     }
1644     else if (strip >= fTOFGeometry->NStripC() + fTOFGeometry->NStripB() &&
1645              strip <  fTOFGeometry->NStripC() + fTOFGeometry->NStripB() + fTOFGeometry->NStripA()) {
1646       plate = 2;
1647       strip = strip - fTOFGeometry->NStripC() - fTOFGeometry->NStripB();
1648     }
1649     else if (strip >= fTOFGeometry->NStripC() + fTOFGeometry->NStripB() + fTOFGeometry->NStripA() &&
1650              strip <  fTOFGeometry->NStripC() + fTOFGeometry->NStripB() + fTOFGeometry->NStripA() + fTOFGeometry->NStripB()) {
1651       plate = 3;
1652       strip = strip - fTOFGeometry->NStripC() - fTOFGeometry->NStripB() - fTOFGeometry->NStripA();
1653     }
1654     else                                {
1655       plate = 4;
1656       strip = strip - fTOFGeometry->NStripC() - fTOFGeometry->NStripB() - fTOFGeometry->NStripA() - fTOFGeometry->NStripB();
1657     }
1658
1659     volpath=gMC->CurrentVolOffName(7);
1660     index=atoi(&volpath[4]);
1661     sector=-1;
1662     sector=index;
1663
1664     //Old 6h convention
1665     // if(index<5){
1666     //   sector=index+13;
1667     //  }
1668     // else{
1669     //   sector=index-5;
1670     // } 
1671  
1672     for(i=0;i<3;++i) {
1673       hits[i]   = pos[i];
1674       hits[i+3] = pm[i];
1675     }
1676
1677     hits[6] = mom.Rho();
1678     hits[7] = pos[3];
1679     hits[8] = xpad[0];
1680     hits[9] = xpad[1];
1681     hits[10]= xpad[2];
1682     hits[11]= incidenceAngle;
1683     hits[12]= gMC->Edep();
1684     hits[13]= gMC->TrackLength();
1685     
1686     vol[0]= sector;
1687     vol[1]= plate;
1688     vol[2]= strip;
1689     vol[3]= padx;
1690     vol[4]= padz;    
1691
1692     AddT0Hit(mcApplication->GetCurrentTrackNumber(),vol, hits);
1693     //AddT0Hit(gAlice->GetMCApp()->GetCurrentTrackNumber(),vol, hits);
1694   }
1695 }
1696 //-------------------------------------------------------------------
1697 void AliTOFv5T0::MaterialMixer(Float_t* p,Float_t* a,Float_t* m,Float_t* d,Float_t* s,Int_t n) const
1698 {
1699   // a[] atomic weights vector  (in)
1700   //     (atoms present in more compound appear separately)
1701   // m[] number of corresponding atoms in the mixture  (in)
1702   // d[] fraction of the compound relative to the corresponding atoms  (in)
1703   // s[] further possible weights     "        "        "        "     (in)
1704   Float_t t = 0.;
1705   for (Int_t i = 0; i < n; ++i) {
1706     p[i] = a[i]*m[i]*d[i]*s[i];
1707     t  += p[i];
1708   }
1709   for (Int_t i = 0; i < n; ++i) {
1710     p[i] = p[i]/t;
1711     //    AliInfo(Form((\n weight[%i] = %f (,i,p[i]));
1712   }
1713 }