changes in the geometry
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALv0.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 /* $Id$ */
17
18 //_________________________________________________________________________
19 // Implementation version v0 of EMCAL Manager class 
20 // An object of this class does not produce hits nor digits
21 // It is the one to use if you do not want to produce outputs in TREEH or TREED
22 //                  
23 //*-- Author: Yves Schutz (SUBATECH)
24
25
26 // --- ROOT system ---
27
28 #include "TTUBS.h"
29 #include "TNode.h"
30 #include "TRandom.h"
31 #include "TGeometry.h"
32
33
34 // --- Standard library ---
35
36 #include <iostream.h> 
37
38 // --- AliRoot header files ---
39
40 #include "AliEMCALv0.h"
41 #include "AliEMCALGeometry.h"
42 #include "AliRun.h"
43 #include "AliMC.h"
44
45 ClassImp(AliEMCALv0)
46
47 //____________________________________________________________________________
48 AliEMCALv0::AliEMCALv0(const char *name, const char *title):
49   AliEMCAL(name,title)
50 {
51   // ctor
52
53 }
54
55 //____________________________________________________________________________
56 void AliEMCALv0::BuildGeometry()
57 {
58
59   const Int_t kColorArm1   = kBlue ;
60   const Int_t kColorArm2   = kBlue ;
61   const Int_t kColorArm1Active   = kRed ;
62   const Int_t kColorArm2Active   = kRed ;
63
64   // make the container of  Arm1
65   
66   new TTUBS("Envelop1", "Tubs that contains arm 1", "void", 
67             fGeom->GetEnvelop(0),     // rmin 
68             fGeom->GetEnvelop(1),     // rmax
69             fGeom->GetEnvelop(2)/2.0, // half length in Z
70             fGeom->GetArm1PhiMin(),   // minimun phi angle
71             fGeom->GetArm1PhiMax()   // maximun phi angle
72             ) ; 
73    // Active material of  Arm1
74  
75   new TTUBS("Arm1", "Active material of  arm 1", "void", 
76             fGeom->GetEnvelop(0),     // rmin 
77             fGeom->GetEnvelop(1),     // rmax
78             fGeom->GetLmat()/2.0,     // half length in Z
79             fGeom->GetArm1PhiMin(),   // minimun phi angle
80             fGeom->GetArm1PhiMax()   // maximun phi angle
81             ) ; 
82   // make the container of  Arm2
83  
84   new TTUBS("Envelop2", "Tubs that contains arm 2", "void", 
85             fGeom->GetEnvelop(0),     // rmin 
86             fGeom->GetEnvelop(1),     // rmax
87             fGeom->GetEnvelop(2)/2.0, // half length in Z
88             fGeom->GetArm2PhiMin(),   // minimun phi angle
89             fGeom->GetArm2PhiMax()   // maximun phi angle
90             ) ;
91            
92   // Active material of  Arm2
93  
94   new TTUBS("Arm2", "Active material of  arm 2", "void", 
95             fGeom->GetEnvelop(0),     // rmin 
96             fGeom->GetEnvelop(1),     // rmax
97             fGeom->GetLmat()/2.0,     // half length in Z           ) ; 
98             fGeom->GetArm2PhiMin(),   // minimun phi angle
99             fGeom->GetArm2PhiMax()   // maximun phi angle
100             ) ;
101
102   TNode * top = gAlice->GetGeometry()->GetNode("alice") ;
103   top->cd();
104   
105   // Arm 1 inside alice
106   TNode * envelop1node = new TNode("Arm1 Envelop", "Arm1 Envelop", "Envelop1") ;
107   envelop1node->SetLineColor(kColorArm1) ;
108   fNodes->Add(envelop1node) ;
109
110   // Arm 2 inside alice
111   TNode * envelop2node = new TNode("Arm2 Envelop", "Arm2 Envelop", "Envelop2") ;
112   envelop2node->SetLineColor(kColorArm2) ;
113   fNodes->Add(envelop2node) ;
114
115   // active material inside Arm 1
116   envelop1node->cd() ; 
117   TNode * arm1node = new TNode("Arm1 Mat", "Arm1 Mat", "Arm1Mat") ;
118   arm1node->SetLineColor(kColorArm1Active) ;
119   fNodes->Add(arm1node) ; 
120
121   // active material inside Arm 2
122   envelop2node->cd() ; 
123   TNode * arm2node = new TNode("Arm2 Mat", "Arm2 Mat", "Arm2Mat") ;
124   arm2node->SetLineColor(kColorArm2Active) ;
125   fNodes->Add(arm2node) ; 
126   
127 }
128
129
130 //____________________________________________________________________________
131 void AliEMCALv0::CreateGeometry()
132 {
133   // Create the EMCAL geometry for Geant
134
135   AliEMCALv0 *emcaltmp = (AliEMCALv0*)gAlice->GetModule("EMCAL") ;
136
137   if ( emcaltmp == NULL ) {
138     
139     fprintf(stderr, "EMCAL detector not found!\n") ;
140     return;
141     
142   }
143   // Get pointer to the array containing media indices
144   Int_t *idtmed = fIdtmed->GetArray() - 1599 ;
145
146
147   // Create a tube sector that contains Arm 1 
148  
149   Float_t envelopA[5] ; 
150   envelopA[0] = fGeom->GetEnvelop(0) ;         // rmin
151   envelopA[1] = fGeom->GetEnvelop(1) ;         // rmax
152   envelopA[2] = fGeom->GetEnvelop(2) / 2.0 ;   // dz
153   envelopA[3] = fGeom->GetArm1PhiMin() ;       // minimun phi angle
154   envelopA[4] = fGeom->GetArm1PhiMax() ;       // maximun phi angle
155
156   gMC->Gsvolu("XEN1", "TUBS ", idtmed[1599], envelopA, 5) ; // filled with air
157
158   // Create a tube sector that contains active material Arm 1 
159  
160   envelopA[2] = fGeom->GetLmat() / 2.0 ;       // dz
161
162   gMC->Gsvolu("XAR1", "TUBS ", idtmed[1601], envelopA, 5) ; // filled with active material (average)
163
164
165   // Create a tube sector that contains Arm 2 
166  
167   envelopA[3] = fGeom->GetArm2PhiMin() ;       // minimun phi angle
168   envelopA[4] = fGeom->GetArm2PhiMax() ;       // maximun phi angle
169
170   gMC->Gsvolu("XEN2", "TUBS ", idtmed[1599], envelopA, 5) ; // filled with air
171
172   // Create a tube sector that contains active material Arm 1 
173  
174   envelopA[2] = fGeom->GetLmat() / 2.0 ;       // dz
175
176   gMC->Gsvolu("XAR2", "TUBS ", idtmed[1601], envelopA, 5) ; // filled with active material (average)
177   
178   Int_t idrotm =1;
179   AliMatrix(idrotm, 90.0, 0., 90.0, 90.0, 0.0, 0.0) ;
180
181   // Position  ENV1 container in ALIC
182   gMC->Gspos("XEN1", 1, "ALIC", 0.0, 0.0, 0.0, idrotm, "ONLY") ;
183   
184   // Position  ARM1  into ENV1
185   gMC->Gspos("XAR1", 1, "XEN1", 0.0, 0.0, 0.0, idrotm, "ONLY") ;
186
187   // Position  ENV2 container in ALIC  
188   gMC->Gspos("XEN2", 1, "ALIC", 0.0, 0.0, 0.0, idrotm, "ONLY") ;
189
190   // Position  ARM2  into ENV2
191   gMC->Gspos("XAR2", 1, "XEN2", 0.0, 0.0, 0.0, idrotm, "ONLY") ;
192   
193 }
194
195 //____________________________________________________________________________
196 void AliEMCALv0::Init(void)
197 {
198   // Just prints an information message
199   
200   Int_t i;
201
202   printf("\n");
203   for(i=0;i<35;i++) printf("*");
204   printf(" EMCAL_INIT ");
205   for(i=0;i<35;i++) printf("*");
206   printf("\n");
207
208   // Here the EMCAL initialisation code (if any!)
209
210   if (fGeom!=0)  
211     cout << "AliEMCAL" << Version() << " : EMCAL geometry intialized for " << fGeom->GetName() << endl ;
212   else
213     cout << "AliEMCAL" << Version() << " : EMCAL geometry initialization failed !" << endl ;   
214   
215   for(i=0;i<80;i++) printf("*");
216   printf("\n");
217   
218 }
219