First version, stand alone detector
[u/mrichter/AliRoot.git] / CRT / AliCRTv1.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
19 */
20
21 ///////////////////////////////////////////////////////////////////////////////
22 //                                                                           //
23 // ALICE Cosmic Ray Trigger                                                  //
24 //                                                                           //
25 //  This class contains the functions for version 0 of the ALICE Cosmic Ray  //
26 //  Trigger. This vesion is suposed to work as standalone module             //
27 //                                                                           //
28 //
29 //   Authors:
30 //
31 //   Arturo Fernandez <afernand@fcfm.buap.mx>
32 //   Enrique Gamez    <egamez@fcfm.buap.mx>
33 //
34 //   Universidad Autonoma de Puebla
35 //
36 //
37 //Begin_Html
38 /*
39 <img src="picts/AliCRTv1Class.gif">
40 </pre>
41 <br clear=left>
42 <p>The responsible person for this module is
43 <a href="mailto:egamez@fcfm.buap.mx">Enrique Gamez</a>.
44 </font>
45 <pre>
46 */
47 //End_Html
48 //                                                                           //
49 ///////////////////////////////////////////////////////////////////////////////
50
51 #include <iostream.h>
52
53 #include <TGeometry.h>
54 #include <TBRIK.h>
55 #include <TNode.h>
56 #include <TLorentzVector.h>
57
58 #include "AliRun.h"
59 #include "AliMC.h"
60 #include "AliMagF.h"
61 #include "AliConst.h"
62 #include "AliPDG.h"
63
64 #include "AliCRTv1.h"
65 #include "AliCRTConstants.h"
66
67 ClassImp(AliCRTv1)
68  
69 //_____________________________________________________________________________
70 AliCRTv1::AliCRTv1() : AliCRTv0()
71 {
72   //
73   // Default constructor for CRT
74   //
75   fCRTStatus = kTRUE;
76   fRICHStatus = kFALSE;
77   fTPCStatus = kFALSE;
78   fMagnetStatus = kTRUE;
79
80   fCRTModule = kFALSE;
81 }
82  
83 //_____________________________________________________________________________
84 AliCRTv1::AliCRTv1(const char *name, const char *title)
85   : AliCRTv0(name,title)
86 {
87   //
88   // Standard constructor for CRT
89   //
90   //Begin_Html
91   /*
92     <img src="picts/AliCRTv1.gif">
93   */
94   //End_Html
95   fCRTStatus = kTRUE;
96   fCRTModule = kFALSE;
97
98   fRICHStatus = kFALSE;
99   fTPCStatus = kFALSE;
100   fMagnetStatus = kFALSE;
101 }
102
103 //_____________________________________________________________________________
104 AliCRTv1::AliCRTv1(const AliCRTv1& crt)
105 {
106   //
107   // Copy ctor.
108   //
109   crt.Copy(*this);
110 }
111
112 //_____________________________________________________________________________
113 AliCRTv1& AliCRTv1::operator= (const AliCRTv1& crt)
114 {
115   //
116   // Asingment operator
117   //
118   crt.Copy(*this);
119   return *this;
120 }
121
122 //_____________________________________________________________________________
123 void AliCRTv1::CreateGeometry()
124 {
125   //
126   // Create geometry for the CRT array
127   //
128
129   Int_t  idrotm[2499];    // The rotation matrix.
130
131   Int_t * idtmed = fIdtmed->GetArray() - 1099 ;
132
133   //
134   // Shafts.
135   this->CreateShafts();
136
137   //
138   // Molasse.
139   this->CreateMolasse();
140
141
142   //
143   // Scintillators
144
145   Float_t box[3];
146   box[0] = AliCRTConstants::fgCageLenght/2.; // Half Length of the box along the X axis, cm.
147   box[1] = AliCRTConstants::fgCageHeight/2.; // Half Length of the box along the Y axis, cm.
148   box[2] = AliCRTConstants::fgCageWidth/2.;  // Half Length of the box along the Z axis, cm.
149
150   //
151   // Create a big voluem with air barrel above the magnet
152   Float_t barrel[10];
153   Float_t magnetSides = 3.;
154   Float_t planesPerpendicularToZ = 2.;
155   Float_t rMin = 790.;
156   Float_t rMax = rMin + 20.; // 20 cm width
157   barrel[0] = 22.5;
158   barrel[1] = 45*magnetSides;
159   barrel[2] = magnetSides;
160   barrel[3] = planesPerpendicularToZ;
161   barrel[4] = -600.;
162   barrel[5] = rMin;
163   barrel[6] = rMax;
164   barrel[7] = 600.;
165   barrel[8] = rMin;
166   barrel[9] = rMax;
167   gMC->Gsvolu("CRT4", "PGON", idtmed[1114], barrel, 10);
168   gMC->Gspos("CRT4", 1 , "CRT", 0., -30., 0., 0, "ONLY");
169   
170
171   // Create  the current sicuiitllator arry
172   // Define the Scintillators. as a big box.
173   Float_t scint[3];
174   scint[0] = AliCRTConstants::fgActiveAreaLenght/2.;   // Half Length in X
175   scint[1] = AliCRTConstants::fgActiveAreaHeight/2.;   // Half Length in Y
176   scint[2] = AliCRTConstants::fgActiveAreaWidth/2.;    // Half Length in Z
177   gMC->Gsvolu("CRT1", "BOX ", idtmed[1112], scint, 3); // Scintillators
178   //
179   // -- X axis.
180   // we'll start dawing from the center.
181   Float_t initX = 0.;
182   
183   //
184   // -- Y axis
185   Float_t gapY   = 30.;        // 30 cms. above the barrel.
186   // For the height we staimate the from the center of the ceiling,
187   // if were a cilinder, must be about 280cm.
188   Float_t barrelc = 790.; // Barrel radius.
189   Float_t height  = barrelc + gapY - 30.;
190   Float_t initY = height;
191   
192   //
193   // -- Z axis.
194   // we'll start dawing from the center.
195   
196   //
197   // Put 4 modules on the top of the magnet
198   Int_t step = 4;
199   for ( Int_t i = 1 ; i <= 4 ; i++ ) {
200     gMC->Gspos("CRT1", i, "CRT", initX, initY, (i-step)*box[2], 0, "ONLY");
201     step--;
202   }
203   
204   // Modules on the barrel sides.
205   // Because the openenig angle for each face is 22.5, and if we want to
206   //    put the modules right in the middle
207   Float_t xtragap = 10.;
208   Float_t initXside = (height+xtragap)*TMath::Sin(2*22.5*kDegrad);//rigthside
209   Float_t initYside = (height+xtragap)*TMath::Cos(2*22.5*kDegrad);
210   
211   // Put 4 modules on the left side of the magnet
212   // The rotation matrix parameters, for the left side.
213   AliMatrix(idrotm[232], 90., 315., 90., 45., 0., 337.5);
214   Int_t stepl = 4;
215   for ( Int_t i = 1 ; i <= 4 ; i++ ) {
216     gMC->Gspos("CRT1", i+4, "CRT", initXside, initYside, (i-stepl)*box[2],
217                idrotm[232], "ONLY");
218     stepl--;
219   }
220   
221   // Put 4 modules on the right side of the magnet
222   // The rotation matrix parameters for the right side.
223   AliMatrix(idrotm[231], 90., 45., 90., 315., 180., 202.5);
224   Int_t stepr = 4;
225   for ( Int_t i = 1 ; i <= 4 ; i++ ) {
226     gMC->Gspos("CRT1", i+8, "CRT", -initXside, initYside, (i-stepr)*box[2],
227                idrotm[231], "ONLY");
228     stepr--;
229   }
230   
231   this->CreateMagnetGeometry();
232   this->CreateRICHGeometry();
233   this->CreateTPCGeometry();
234   
235 }
236
237 //_____________________________________________________________________________
238 void AliCRTv1::CreateMagnetGeometry()
239 {
240   
241   cout<<"\n\n\tYou are requiring the CRT with the Magnet Activated!\n\n";
242   
243   Int_t  idrotm[2499];    // The rotation matrix.
244
245   Int_t * idtmed = fIdtmed->GetArray() - 1099 ;
246
247   // Disable the CRT StepManager method.
248   fCRTStatus = kFALSE;
249
250   Float_t barrel[10];
251   Float_t magnetSides = 3.;
252   Float_t planesPerpendicularToZ = 2.;
253   //Float_t rMin = 790.;
254   //Float_t rMax = rMin + 20.; // 20 cm width
255
256   // MAgnet
257   // Create the upper faces of the magnet.
258   barrel[0] = 22.5;
259   barrel[1] = 360.;
260   barrel[2] = 8.;
261   barrel[3] = 2.;
262   barrel[4] = -600.;
263   barrel[5] = 580.;
264   barrel[6] = 790.;
265   barrel[7] = 600.;
266   barrel[8] = 580.;
267   barrel[9] = 790.;
268   gMC->Gsvolu("C3MO", "PGON", idtmed[1114], barrel, 10);
269   gMC->Gspos("C3MO", 1, "CRT", 0., -30., 0., 0, "ONLY");
270
271   // Define coils 
272   
273   barrel[5] = 585.;
274   barrel[6] = 690.;
275   barrel[8] = 585.;
276   barrel[9] = 690.;
277   gMC->Gsvolu("C3CO", "PGON", idtmed[1108], barrel, 10); //Aluminium
278   gMC->Gspos("C3CO", 1, "C3MO", 0., 0., 0., 0, "ONLY");
279   
280   barrel[5] = 580.;
281   barrel[6] = 585.;
282   barrel[8] = 580.;
283   barrel[9] = 585.;
284   gMC->Gsvolu("C3C1", "PGON", idtmed[1128], barrel, 10);// Aluminium
285   gMC->Gspos("C3C1", 1, "C3MO", 0., 0., 0., 0, "ONLY");
286
287   // Define yoke 
288   
289   barrel[5] = 690.;
290   barrel[6] = 790.;
291   barrel[8] = 690.;
292   barrel[9] = 790.;
293   gMC->Gsvolu("C3YO", "PGON", idtmed[1109], barrel, 10); // Iron
294   gMC->Gspos("C3YO", 1, "C3MO", 0., 0., 0., 0, "ONLY");
295
296
297   // Now create one inside the magnet as L3C1
298   // voulme for tracking.
299   barrel[0] = 22.5;
300   barrel[1] = 45*magnetSides;
301   barrel[2] = magnetSides;
302   barrel[3] = planesPerpendicularToZ;
303   barrel[4] = -600.;
304   barrel[5] = 575.;
305   barrel[6] = 580.;
306   barrel[7] = 600.;
307   barrel[8] = 575.;
308   barrel[9] = 580.;
309   gMC->Gsvolu("C3CI", "PGON", idtmed[1134], barrel, 10);
310   gMC->Gspos("C3CI", 1 , "CRT", 0., -30., 0., 0, "ONLY");
311
312   // And a detector layer in the door 10 cm thick
313   // Volume for tracking.
314   barrel[0] = 22.5;
315   barrel[1] = 360.;
316   barrel[2] = 8.;
317   barrel[3] = 2.;
318   barrel[4] = 590.;
319   barrel[5] = 0.;
320   barrel[6] = 580.;
321   barrel[7] = 600.;
322   barrel[8] = barrel[5];
323   barrel[9] = barrel[6];
324   gMC->Gsvolu("C3C2", "PGON", idtmed[1154], barrel, 10); // Air
325   gMC->Gspos("C3C2", 1, "CRT",  0., -30., 0., 0, "ONLY");
326   AliMatrix(idrotm[1010], 90., 0., 90., 90., 180., 0.);
327   gMC->Gspos("C3C2", 2, "CRT",  0., -30., 0., idrotm[1010], "ONLY");
328
329
330
331   barrel[4] = 600.;
332   barrel[5] = 0.;
333   barrel[6] = 790.;
334   barrel[7] = 700.;
335   barrel[8] = barrel[5];
336   barrel[9] = barrel[6];
337   gMC->Gsvolu("C3DO", "PGON", idtmed[1174], barrel, 10); // Air
338   gMC->Gspos("C3DO", 1, "CRT", 0., -30., 0., 0, "ONLY");
339   AliMatrix(idrotm[1010], 90., 0., 90., 90., 180., 0.);
340   gMC->Gspos("C3DO", 2, "CRT", 0., -30., 0., idrotm[1010], "ONLY");
341
342   barrel[4] = 610.;
343   barrel[5] = 0.;
344   barrel[6] = 790.;
345   barrel[7] = 700.;
346   barrel[8] = barrel[5];
347   barrel[9] = barrel[6];
348   gMC->Gsvolu("C3FR", "PGON", idtmed[1149], barrel, 10); // Iron
349   gMC->Gspos("C3FR", 1, "C3DO", 0., 0., 0., 0, "ONLY");
350   // INNER LAYER 
351   
352   barrel[4] = 600.;
353   barrel[7] = 610.;
354   gMC->Gsvolu("C3IR", "PGON", idtmed[1149], barrel, 10); //Iron
355   gMC->Gspos("C3IR", 1, "C3DO", 0., 0., 0., 0, "ONLY");
356
357 }
358
359 //_____________________________________________________________________________
360 void AliCRTv1::CreateTPCGeometry()
361 {
362   cout<<"\n\n\tYou are requiring the CRT with the TPC Activated!\n\n";
363   Int_t * idtmed = fIdtmed->GetArray() - 1099 ;
364
365   // Disable the CRT StepManager method.
366   fCRTStatus = kFALSE;
367   // Disable the MAgnet
368   fMagnetStatus = kFALSE;
369   // Disable th RICH
370   fRICHStatus = kFALSE;
371
372   // TPC
373   // Tpc SAndwich 1 - Al
374   // TSA1
375   Float_t tube[5];
376   tube[0]=274.8124;
377   tube[1]=278.;
378   tube[2]=252.1;
379   tube[3] = 0.;
380   tube[4] = 180.;
381   gMC->Gsvolu("CSA1","TUBS",idtmed[1154],tube,5);
382   // TSA1->TOCV (0.,0.,3.) ->TOIN (0.,0.,0.)->TPC (0.,0.,0.)->ALIC(0.,0.,0.)
383   gMC->Gspos("CSA1 ",1,"CRT",0.,0.,0.,0,"ONLY");
384
385 }
386
387 //_____________________________________________________________________________
388 void AliCRTv1::CreateRICHGeometry()
389 {
390
391   cout<<"\n\n\tYou are requiring the CRT with the RICH Activated!\n\n";
392
393   Int_t  idrotm[2499];    // The rotation matrix.
394
395   Int_t * idtmed = fIdtmed->GetArray() - 1099 ;
396
397   // Disable the CRT StepManager method.
398   fCRTStatus = kFALSE;
399   // Disable the MAgnet
400   fMagnetStatus = kFALSE;
401
402
403   // now create  volume to simulate the HMPID volume. CSI
404   Float_t csi_length = 160*.8 + 2.6;
405   Float_t csi_width = 144*.84 + 2*2.6;
406   Float_t tbox[3];
407   tbox[0] = csi_width/2;
408   tbox[1] = 11.5;
409   tbox[2] = csi_length/2;
410   gMC->Gsvolu("CRIC ", "BOX ", idtmed[1174], tbox, 3);
411
412   Double_t dOffset = 490+1.267 - 8/2;  // distance from center of mother volume ALIC to methane
413   
414   Double_t dAlpha = 19.5; // angle between centers of chambers - y-z plane
415   Double_t dAlphaRad = dAlpha*kDegrad;
416   
417   Double_t dBeta = 20.;   // angle between center of chambers - y-x plane
418   Double_t dBetaRad = dBeta*kDegrad;
419    
420   Double_t dRotAngle = 60.;     // the whole RICH is to be rotated in x-y plane + means clockwise rotation 
421   Double_t dRotAngleRad = dRotAngle*kDegrad;
422    
423     
424    TRotMatrix *pRotMatrix; // tmp pointer
425    
426    TVector3 vector(0,dOffset,0); // Position of chamber 2 without rotation
427
428    // Chamber 0  standalone (no other chambers in this row) 
429    AliMatrix(idrotm[1000],90, -dRotAngle+360,90-dAlpha, 90-dRotAngle, dAlpha, -90+300);
430    pRotMatrix=new TRotMatrix("rot993","rot993",90,-dRotAngle, 90-dAlpha,90-dRotAngle,dAlpha, -90);
431    
432    vector.SetXYZ(0,dOffset,0);  vector.RotateX(dAlphaRad); 
433    vector.RotateZ(-dRotAngleRad);
434    
435    gMC->Gspos("CRIC",1,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1000], "ONLY");
436    
437    // Chamber 1   
438    AliMatrix(idrotm[1001],90,-dBeta-dRotAngle,90,90-dBeta-dRotAngle, 0,0);
439
440    pRotMatrix=new TRotMatrix("rot994","rot994",90,-dBeta-dRotAngle,90,90-dBeta-dRotAngle,0,0);  
441    
442    vector.SetXYZ(0,dOffset,0);  vector.RotateZ(-dBetaRad); 
443    vector.RotateZ(-dRotAngleRad);
444    
445    gMC->Gspos("CRIC",2,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1001], "ONLY");           
446    
447    // Chamber 2   the top one with no Alpha-Beta rotation
448    AliMatrix(idrotm[1002],90,-dRotAngle,90,90-dRotAngle,0,0);
449
450    pRotMatrix=new TRotMatrix("rot995","rot995",90,-dRotAngle,90,90-dRotAngle,0,0);
451    
452    vector.SetXYZ(0,dOffset,0);
453    vector.RotateZ(-dRotAngleRad);
454    
455    gMC->Gspos("CRIC",3,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1002], "ONLY");           
456    
457    // Chamber 3
458    AliMatrix(idrotm[1003],90,dBeta-dRotAngle,90.,90+dBeta-dRotAngle,0,0);
459    pRotMatrix=new TRotMatrix("rot996","rot996", 90,dBeta-dRotAngle,90.,90+dBeta-dRotAngle,0,0);
460    
461    vector.SetXYZ(0,dOffset,0);  vector.RotateZ(dBetaRad); 
462    vector.RotateZ(-dRotAngleRad);
463    
464    gMC->Gspos("CRIC",4,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1003], "ONLY");
465
466    // Chamber 4
467    AliMatrix(idrotm[1004],90,360-dBeta-dRotAngle,108.2,90-dBeta-dRotAngle,18.2,90-dBeta-60);
468    pRotMatrix=new TRotMatrix("rot997","rot997",90,360-dBeta-dRotAngle,108.2,90-dBeta-dRotAngle,18.2,90-dBeta);
469    
470    vector.SetXYZ(0,dOffset,0);  vector.RotateZ(-dBetaRad); vector.RotateX(-dAlphaRad); 
471    vector.RotateZ(-dRotAngleRad);
472    
473    gMC->Gspos("CRIC",5,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1004], "ONLY");
474
475    // Chamber 5
476    AliMatrix(idrotm[1005],90,-dRotAngle+360,90+dAlpha,90-dRotAngle,dAlpha,90-60);     
477
478    pRotMatrix=new TRotMatrix("rot998","rot998",90,-dRotAngle,90+dAlpha,90-dRotAngle,dAlpha,90);     
479    
480    vector.SetXYZ(0,dOffset,0); vector.RotateX(-dAlphaRad); 
481    vector.RotateZ(-dRotAngleRad);
482    
483    gMC->Gspos("CRIC",6,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1005], "ONLY");           
484    
485    // Chamber 6           
486    AliMatrix(idrotm[1006],90,dBeta-dRotAngle+360,108.2,90+dBeta-dRotAngle,18.2,90+dBeta-60);
487
488    pRotMatrix=new TRotMatrix("rot999","rot999",90,dBeta-dRotAngle,108.2,90+dBeta-dRotAngle,18.2,90+dBeta);    
489    
490    vector.SetXYZ(0,dOffset,0);  vector.RotateZ(dBetaRad); vector.RotateX(-dAlphaRad); 
491    vector.RotateZ(-dRotAngleRad);
492    
493    gMC->Gspos("CRIC",7,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1006], "ONLY");
494    
495 }
496
497 //_____________________________________________________________________________
498 void AliCRTv1::CreateMolasse()
499 {
500   //
501   //
502   //
503
504   Int_t  idrotm[2499];    // The rotation matrix.
505
506   Int_t * idtmed = fIdtmed->GetArray() - 1099 ;
507
508   //
509   // Molasse
510   //
511
512   // Exactly above the hall
513   Float_t tspar[5];
514   tspar[0] = 1170.;
515   tspar[1] = 1170. + 375.;
516   tspar[2] = (1900.+1150.)/2.+100.;
517   tspar[3] = 0.;
518   tspar[4] = 180.;
519   gMC->Gsvolu("CMO1", "TUBS", idtmed[1123], tspar, 5);
520   gMC->Gspos("CMO1", 1, "CRT", 0., 500., 1900.-tspar[2]+400., 0, "MANY");
521
522   Float_t tbox[3];
523   tbox[0] = 1250.;
524   tbox[1] = (4420. - 1670.)/2.;
525   tbox[2] = (1900.+1150.)/2. + 200.;
526   gMC->Gsvolu("CM12", "BOX", idtmed[1123], tbox, 3);
527   gMC->Gspos("CM12", 1, "CRT", 0., 4420. -tbox[1], 1900.-tbox[2]+400., 0, "MANY");
528
529   AliMatrix(idrotm[2003], 0., 0., 90., 0., 90., 90.);
530   // Along the PM25
531   Float_t tube[3];
532   tube[0] = 455. + 100.;
533   tube[1] = 555. + 375.;
534   tube[2] = (5150. - 1166.)/2.;
535   gMC->Gsvolu("CMO2", "TUBE", idtmed[1123], tube, 3);
536   gMC->Gspos("CMO2", 1, "CRT", -2100., 4420.-tube[2], 0., idrotm[2003], "MANY");
537
538
539   // Along the PGC2
540   tube[0] = 650.;
541   tube[1] = 2987.7;
542   tube[2] = (5150. - 690.)/2.;
543   gMC->Gsvolu("CMO3", "TUBE", idtmed[1123], tube, 3);
544   gMC->Gspos("CMO3", 1, "CRT", 375., 4420.-tube[2], 1900.+2987.7, idrotm[2003], "MANY");
545   // Behind the PGC2 up to the end of the M. volume.
546   tbox[0] = 12073.;
547   tbox[1] = 2575. + 95.;
548   tbox[2] = (12073. - 1900.-2987.7-650.)/2.;
549   gMC->Gsvolu("CMO7", "BOX", idtmed[1123], tbox, 3);
550   gMC->Gspos("CMO7", 1, "CRT", 0., 4420.-tbox[1], 1900.+2987.7+650.+tbox[2], 0, "MANY");
551
552   // Along the PX24 , upper part.
553   tube[0] = 1250.;
554   tube[1] = 2300;
555   tube[2] = 2575. - 1300. + 95.;
556   gMC->Gsvolu("CMO4", "TUBE", idtmed[1123], tube, 3);
557   gMC->Gspos("CMO4", 1, "CRT", 0., 404.+1300.+tube[2], -2300., idrotm[2003], "MANY");
558
559   // Along the PX24 , lower part
560   tspar[0] = 1250.;
561   tspar[1] = 2300;
562   tspar[2] = 1300.;
563   tspar[3] = kRaddeg*TMath::ASin(1070./1150.);
564   tspar[4] = 360. - tspar[3];
565   gMC->Gsvolu("CMO5", "TUBS", idtmed[1123], tspar, 5);
566   gMC->Gspos("CMO5", 1, "CRT", 0., 404., -2300., idrotm[2003], "MANY");
567   // behind the PX24
568   tbox[0] = 12073.;
569   tbox[1] = 2575. + 95.;
570   tbox[2] = 8523./2.;
571   gMC->Gsvolu("CMO6", "BOX", idtmed[1123], tbox, 3);
572   gMC->Gspos("CMO6", 1, "CRT", 0., 4420.-tbox[1], -3550.-tbox[2], 0, "MANY");
573
574
575   // On the right side of th hall
576   tbox[0] = (12073. - 1250.)/2.;
577   tbox[1] = 2575. + 95.;
578   tbox[2] = (8437.7+650.)/2.;
579   gMC->Gsvolu("CMO8", "BOX", idtmed[1123], tbox, 3);
580   gMC->Gspos("CMO8", 1, "CRT", 1250.+tbox[0], 4420.-tbox[1], -3550.+tbox[2], 0, "MANY");
581
582   // on the left side of the hall, behind 
583   tbox[0] = (12073. - 2755.)/2.;
584   tbox[1] = 2575. + 95.;
585   tbox[2] = (8437.7+650.)/2.;
586   gMC->Gsvolu("CMO9", "BOX", idtmed[1123], tbox, 3);
587   gMC->Gspos("CMO9", 1, "CRT", -2755.-tbox[0], 4420.-tbox[1], -3550.+tbox[2], 0, "MANY");
588
589
590   // Molasse betwen the PX24 & PM25 on the left side.
591   tbox[0] = (2755. - 1250.)/2.;
592   tbox[1] = 2575. + 95.;
593   tbox[2] = (3550. - 555.)/2.;
594   gMC->Gsvolu("CM10", "BOX", idtmed[1123], tbox, 3);
595   gMC->Gspos("CM10", 1, "CRT", -1250.-tbox[0], 4420.-tbox[1], -tbox[2]-555., 0, "MANY");
596
597
598   // Molasse betwen the PGC2 & PM25 on the left side.
599   tbox[0] = (2755. - 1250.)/2.;
600   tbox[1] = 2575. + 95.;
601   tbox[2] = (1900.+2987.7 - 555. + 650.)/2.;
602   gMC->Gsvolu("CM11", "BOX", idtmed[1123], tbox, 3);
603   gMC->Gspos("CM11", 1, "CRT", -1250.-tbox[0], 4420.-tbox[1], 555.+tbox[2], 0, "MANY");
604
605
606 }
607
608 //_____________________________________________________________________________
609 void AliCRTv1::CreateShafts()
610 {
611   //
612   //
613   //
614   Int_t  idrotm[2499];    // The rotation matrix.
615
616   Int_t * idtmed = fIdtmed->GetArray() - 1099 ;
617
618   // Create a mother volume.
619   Float_t pbox[3];
620   //pbox[0] = AliCRTConstants::fgDepth*TMath::Tan(67.5*kDegrad);
621   pbox[0] = 12073.;
622   pbox[1] = AliCRTConstants::fgDepth;
623   pbox[2] = pbox[0];
624   gMC->Gsvolu("CRT", "BOX", idtmed[1114], pbox, 3);
625   gMC->Gspos("CRT", 1, "ALIC", 0., 0., 0., 0, "ONLY");
626
627   // HAll ceiling
628   Float_t ptubs[5];
629   ptubs[0] = 1070.;
630   ptubs[1] = 1170.;
631   ptubs[2] = 1900.;
632   ptubs[3] = 0.;
633   ptubs[4] = 180.;
634   gMC->Gsvolu("CHC1", "TUBS", idtmed[1116], ptubs, 5);
635   gMC->Gspos("CHC1", 1, "CRT", 0., 500., 0., 0, "ONLY");
636
637
638   //
639   // Acces shafts
640   //
641   AliMatrix(idrotm[2001], 0., 0., 90., 0., 90., 90.);
642   
643   // PX24
644   ptubs[0] = 1150.;
645   ptubs[1] = 1250.;
646   ptubs[2] = 1300.;
647   ptubs[3] = kRaddeg*TMath::ASin(1070./ptubs[0]);
648   ptubs[4] = 360 - ptubs[3];
649   gMC->Gsvolu("CSF1", "TUBS", idtmed[1116], ptubs, 5);
650   gMC->Gspos("CSF1", 1, "CRT", 0., 404., -2300., idrotm[2001], "MANY");
651
652   Float_t ptube[3];
653   ptube[0] = ptubs[0];
654   ptube[1] = ptubs[1];
655   ptube[2] = 2575. - ptubs[2] + 95.;
656   gMC->Gsvolu("CSF2", "TUBE", idtmed[1116], ptube, 3);
657   gMC->Gspos("CSF2", 1, "CRT", 0., 404.+ptubs[2]+ptube[2], -2300., idrotm[2001], "MANY");
658   
659   // Concrete walls along the shaft
660   pbox[0] = 585./2.;
661   pbox[1] = 2575. + 95.;
662   pbox[2] = 20.;
663   gMC->Gsvolu("CSW1", "BOX", idtmed[1116], pbox, 3);
664   gMC->Gspos("CSW1", 1, "CRT", -290-pbox[0], 404.-1300.+pbox[1], -3450.+210.*2, 0, "MANY");
665   
666   //
667   pbox[0] = 750./2.;
668   pbox[1] = 2575. + 95.;
669   pbox[2] = 20.;
670   gMC->Gsvolu("CSW3", "BOX", idtmed[1116], pbox, 3);
671   gMC->Gspos("CSW3", 1, "CRT", 420.-290.+pbox[0], 404.-1300.+pbox[1], -3450.+210.*2, 0, "MANY");
672   
673   //
674   pbox[0] = 60.;
675   pbox[1] = 2575. + 95.;
676   pbox[2] = 210.;
677   gMC->Gsvolu("CSW2", "BOX", idtmed[1116], pbox, 3);
678   gMC->Gspos("CSW2", 1, "CRT", -290-pbox[0], 404.-1300.+pbox[1], -3450.+pbox[2], 0, "MANY");
679   gMC->Gspos("CSW2", 2, "CRT", 420.-290.+pbox[0], 404.-1300.+pbox[1], -3450.+pbox[2], 0, "MANY");
680   
681   
682   // 
683   pbox[0] = 1000.;
684   pbox[1] = 80.;
685   pbox[2] = 200.;
686   gMC->Gsvolu("CSP1", "BOX", idtmed[1116], pbox, 3);
687   gMC->Gspos("CSP1", 1, "CRT", 0., 2600.-700., -1150-pbox[2], 0, "MANY");
688   
689   //
690   pbox[0] = 340.8;
691   pbox[1] = 300./2.;
692   pbox[2] = 460./2.;
693   gMC->Gsvolu("CSP2", "BOX", idtmed[1116], pbox, 3);
694   gMC->Gspos("CSP2", 1, "CRT", 0., 2950.-700., -3450+pbox[2], 0, "MANY");
695   
696   //
697   pbox[0] = 600.;
698   pbox[1] = 150.;
699   pbox[2] = 75.;
700   gMC->Gsvolu("CSP3", "BOX", idtmed[1116], pbox, 3);
701   gMC->Gspos("CSP3", 1, "CRT", 0., 2950.-700., -1150.-210.-pbox[2], 0, "MANY");
702   
703   //
704   pbox[0] = 600.;
705   pbox[1] = 250.;
706   pbox[2] = 38.;
707   gMC->Gsvolu("CSP4", "BOX", idtmed[1116], pbox, 3);
708   gMC->Gspos("CSP4", 1, "CRT", 0., 2950.-700.+155.+pbox[1], -1150.-210.-pbox[2], 0, "MANY");
709   
710   
711   // Shielding plug
712   pbox[0] = 850.;
713   pbox[1] = 90.;
714   pbox[2] = 720.;
715   gMC->Gsvolu("CSP5", "BOX", idtmed[1116], pbox, 3);
716   gMC->Gspos("CSP5", 1, "CRT", 0., 2950.-700., -3450.+460.+pbox[2], 0, "MANY");
717   
718   //
719   pbox[0] = 80.;
720   pbox[1] = 150.;
721   pbox[2] = 720.;
722   gMC->Gsvolu("CSP6", "BOX", idtmed[1116], pbox, 3);
723   gMC->Gspos("CSP6", 1, "CRT", 1150.-600., 2950.-700., -3450.+460.+pbox[2], 0, "MANY");
724   gMC->Gspos("CSP6", 2, "CRT", -1150.+600., 2950.-700., -3450.+460.+pbox[2], 0, "MANY");
725   
726   
727   //
728   pbox[0] = 130.;
729   pbox[1] = 60.;
730   pbox[2] = 750.;
731   gMC->Gsvolu("CSP7", "BOX", idtmed[1116], pbox, 3);
732   gMC->Gspos("CSP7", 1, "CRT", 850.+pbox[0], 2950.-700.+100., -3450.+460.+pbox[2], 0, "MANY");
733   gMC->Gspos("CSP7", 2, "CRT", -850.-pbox[0], 2950.-700.+100., -3450.+460.+pbox[2], 0, "MANY");
734   
735   
736   // PM25 Acces Shaft
737   ptube[0] = 910./2.;
738   ptube[1] = ptube[0] + 100.;
739   ptube[2] = (5150. - 1166.)/2.;
740   gMC->Gsvolu("CSF3", "TUBE", idtmed[1116], ptube, 3);
741   gMC->Gspos("CSF3", 1, "CRT", -2100., AliCRTConstants::fgDepth-ptube[2], 0., idrotm[2001], "MANY");
742   
743   // PGC2 Access Shaft
744   ptube[0] = 1100./2.;
745   ptube[1] = ptube[0] + 100.;
746   ptube[2] = (5150. - 690.)/2.;
747   gMC->Gsvolu("CSF4", "TUBE", idtmed[1116], ptube, 3);
748   gMC->Gspos("CSF4", 1, "CRT", 375., AliCRTConstants::fgDepth-ptube[2], 1900.+2987.7, idrotm[2001], "MANY");
749
750 }
751
752 //_____________________________________________________________________________
753 void AliCRTv1::DrawDetector()
754 {
755   //
756   // Draw a shaded view of the L3 magnet
757   //
758   cout << "AliCRTv1::DrawModule() : Drawing the module" << endl;
759   
760   
761   Int_t able = 1;
762   Int_t enable = 0;
763   gMC->Gsatt("*", "seen", -1);
764   gMC->Gsatt("alic", "seen", 0);
765   
766   gMC->Gsatt("ALIC","seen",enable);
767   gMC->Gsatt("CRT", "seen",enable);
768   gMC->Gsatt("L3MO","seen",  able); // L3 Magnet
769   //gMC->Gsatt("CRT1","seen",  able); // Scintillators
770   gMC->Gsatt("CRT4","seen",  able); // Scintillators barrel
771   
772   // Draw the molasse volumes
773   gMC->Gsatt("CMO1","seen",enable); // Exactly above the HALL
774   gMC->Gsatt("CMO2","seen",enable); // Molasse, along the PM25
775   gMC->Gsatt("CMO3","seen",enable); // molasse along the PGC2
776   gMC->Gsatt("CMO4","seen",enable); // Molasse, behind the PX24 upper part
777   gMC->Gsatt("CMO5","seen",enable); // molasse behind px24, lower part
778   gMC->Gsatt("CMO6","seen",enable); // behind the PX24
779   gMC->Gsatt("CMO7","seen",enable); // behind the PGC2
780   gMC->Gsatt("CMO8","seen",enable); // on the right side.
781   gMC->Gsatt("CMO9","seen",enable); // on the left side.
782   gMC->Gsatt("CM10","seen",enable); // betwen PX24 & PM25.
783   gMC->Gsatt("CM11","seen",enable); // betwen PGC2 & PM25.
784   gMC->Gsatt("CM12","seen",enable); // box above the hall.
785   
786   gMC->Gdopt("hide", "on");
787   gMC->Gdopt("edge","off");
788   gMC->Gdopt("shad", "on");
789   gMC->Gsatt("*", "fill", 7);
790   gMC->SetClipBox("ALIC", 0, 3000, -3000, 3000, -6000, 6000);
791   gMC->DefaultRange();
792   gMC->Gdraw("alic", 40, 30, 0, 10, 9.5, .009, .009);
793   gMC->Gdhead(1111, "View of CRT(ACORDE)");
794   gMC->Gdman(18, 4, "MAN");
795   
796   
797 }
798
799 //_____________________________________________________________________________
800 void AliCRTv1::Init()
801 {
802   //
803   // Initialise L3 magnet after it has been built
804   Int_t i;
805   //
806   if(fDebug) {
807     printf("\n%s: ",ClassName());
808     for(i=0;i<35;i++) printf("*");
809     printf(" CRTv1_INIT ");
810     for(i=0;i<35;i++) printf("*");
811     printf("\n%s: ",ClassName());
812     //
813     // Here the CRTv1 initialisation code (if any!)
814     for(i=0;i<80;i++) printf("*");
815     printf("\n");
816   }
817
818 }
819
820 //____________________________________________________________________________
821 void AliCRTv1::StepManager()
822 {
823   //
824   // Called for every step in the Cosmic Ray Trigger
825   //
826   static Int_t   vol[5];
827   Int_t          ipart;
828   TLorentzVector pos;
829   TLorentzVector mom;
830
831   static Float_t hits[14];
832   static Float_t eloss;
833   static Float_t elossMag;
834
835   if ( !gMC->IsTrackAlive() ) return;
836
837   if (gMC->IsNewTrack()) {
838     // Reset the deposited energy
839     eloss = 0.;
840     elossMag = 0.;
841   }
842
843   // Add th energy loss in each step.
844   eloss += gMC->Edep();
845
846   gMC->TrackPosition(pos);
847
848   //
849   // CRT
850   //
851
852   if ( gMC->IsTrackEntering() && (strcmp(gMC->CurrentVolName(),"CRT4") == 0)
853        &&(gMC->TrackPid() == kMuonMinus || gMC->TrackPid() == kMuonPlus) ) {
854     
855     // Get current particle id(ipart),track position (pos) and momentum (mom)
856     gMC->TrackPosition(pos);
857     gMC->TrackMomentum(mom);
858     ipart = gMC->TrackPid();
859     
860     vol[0] = 1;
861     vol[1] = 0;
862     vol[2] = 0;
863     vol[3] = 0;
864     vol[4] = 0;
865
866     ipart = gMC->TrackPid();
867     hits[0]  = (Float_t)ipart; //                 (fId)
868     
869     hits[1]  = pos[0]; // X coordinate (fX)
870     hits[2]  = pos[1]; // Y coordinate (fY)
871     hits[3]  = pos[2]; // Z coordinate (fZ)
872     hits[4]  = mom[0]; // Px           (fpxug)
873     hits[5]  = mom[1]; // Py           (fpyug)
874     hits[6]  = mom[2]; // Pz           (fpzug)
875     
876     hits[7]  = gMC->GetMedium();  //layer(flay)
877     hits[8] = eloss;              // Energy loss
878     
879     hits[9] = 1; // CRT mother activated.
880     hits[10] = 0;
881     hits[11] = 0;
882     hits[12] = 0;
883     hits[13] = 0;
884
885     //hits[9] = gAlice->CurrentTrack();
886     
887     AddHit(gAlice->CurrentTrack(),vol, hits);
888     
889     eloss = 0.;
890
891   } else if (gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"CRT1")==0)
892              &&(gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus)) {
893     
894     vol[0] = 0;
895     vol[1] = 1;
896     vol[2] = 0;
897     vol[3] = 0;
898     vol[4] = 0;
899
900     hits[9] = 0; // CRT mother activated.
901     hits[10] = 1;
902     hits[11] = 0;
903     hits[12] = 0;
904     hits[13] = 0;
905
906     //hits[10] = 1;
907     
908     //AddHit(gAlice->CurrentTrack(),vol, hits);
909     
910     //eloss = 0.;
911
912
913   } else if (gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"C3CI")==0)
914       &&(gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus)) {
915
916     //
917     // Inside the magnet, upper part.
918     //
919   
920     // Get current particle id(ipart),track position (pos) and momentum (mom)
921
922     vol[0] = 0;
923     vol[1] = 0;
924     vol[2] = 1;
925     vol[3] = 0;
926     vol[4] = 0;
927
928     hits[9] = 0; // CRT mother activated.
929     hits[10] = 0;
930     hits[11] = 1;
931     hits[12] = 0;
932     hits[13] = 0;
933       
934     AddHit(gAlice->CurrentTrack(),vol, hits);
935     
936     //eloss = 0.;
937
938   } else if ( gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"CRIC")==0)
939        && (gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus) ) {
940
941     //
942     // HMPID
943     //
944     
945     // Get current particle id(ipart),track position (pos) and momentum (mom)
946
947     vol[0] = 0;
948     vol[1] = 0;
949     vol[2] = 0;
950     vol[3] = 1;
951     vol[4] = 0;
952
953     hits[9] = 0;
954     hits[10] = 0;
955     hits[11] = 0;
956     hits[12] = 1;
957     hits[13] = 0;
958     
959     AddHit(gAlice->CurrentTrack(),vol, hits);
960     
961     //eloss = 0.;
962
963
964   } else if (gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"CSA1")==0)
965              &&(gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus)) {
966
967     //
968     // TPC
969     //
970     
971     // Get current particle id(ipart),track position (pos) and momentum (mom)
972     
973     vol[0] = 0;
974     vol[1] = 0;
975     vol[2] = 0;
976     vol[3] = 0;
977     vol[4] = 1;
978
979     hits[9] = 0;
980     hits[10] = 0;
981     hits[11] = 0;
982     hits[12] = 0;
983     hits[13] = 1;
984
985     
986     AddHit(gAlice->CurrentTrack(),vol, hits);
987     
988     //eloss = 0.;
989
990   } else {
991     return;
992   }
993
994
995 }