/[PAMELA software]/PamVMC/cal/src/PamVMCDetCalo.cxx
ViewVC logotype

Contents of /PamVMC/cal/src/PamVMCDetCalo.cxx

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.5 - (show annotations) (download)
Fri Jun 12 18:39:11 2009 UTC (15 years, 5 months ago) by pam-rm2
Branch: MAIN
CVS Tags: v1r0, HEAD
Changes since 1.1: +0 -0 lines 
- Introduced user-defined names of output files and random seeds number.
Users can do it use options of PamVMCApplication constructor:
PamVMCApplication(const char* name,  const char *title, const char*
filename="pamtest", Int_t seed=0).
The Random object that I use is TRandom3 object which has astronomical
large period (in case of default initialization 0). All random generators
in the code use this object by calling of gRandom singleton which keeps
it.

- Corrected TOF digitization routine. No problems with TDC hits due to
hadronic interactions anymore.

- Some small changes was done to compile code under Root 5.23. +
geant4_vmc v. 2.6 without any warnings

- Some classes of PamG4RunConfiguartion was changed for geant4_vmc v.
2.6.Some obsolete classes was deleted as soon as developers implemented
regions.

- Navigation was changed from "geomRootToGeant4" to "geomRoot", because on
VMC web page written that as soon as Geant4 has no option ONLY/MANY
translation of overlapped geometry to Geant4 through VGM could be wrong.
I'd like to stay with Root navigation:
http://root.cern.ch/root/vmc/Geant4VMC.html. This should be default
option.

- New Tracker digitization routine written by Sergio was implemented

- PamVMC again became compatible with geant4_vmc v.2.5 and ROOT 5.20.
 The problem was that ROOT developers introduced in TVirtualMC class a new
method SetMagField and new base class:TVirtualMagField from which
user-defined classes shoukd be derived
1 #include "PamVMCDetCalo.h"
2
3 ClassImp(PamVMCDetCalo)
4
5 void PamVMCDetCalo::DefineGeometry(){
6
7
8 TGeoMaterial *m=0;
9 TGeoMixture *x=0;
10 TString s = "ALUMINIUM";
11 m = GetMat(s);
12 if(!m){
13 m = new TGeoMaterial(s,26.98,13,2.7);
14 m->SetUniqueID( 9);
15 }
16
17 s = "SCINT";
18 x = (TGeoMixture*)GetMat(s); //!!!! CHECK IT NOT SCINTILLATOR !!!
19 if(!x){
20 x = new TGeoMixture(s,2, 1.03200 );
21 x->SetUniqueID( 17);
22 x->DefineElement(0,12.01,6,0.922427);
23 x->DefineElement(1,1.01,1,0.7757296E-01);
24 }
25
26 s = "N2GAS";
27 m = GetMat(s);
28 if (!m){
29 m = new TGeoMaterial(s,14.01,7,0.1250000E-02);
30 m->SetUniqueID( 19);
31 }
32
33 s = "SILICON CA";
34 m = GetMat(s);
35 if(!m){
36 m = new TGeoMaterial(s,28.09,14,2.33);
37 m->SetUniqueID( 21);
38 }
39
40 s = "TUNGA";
41 x = (TGeoMixture*)GetMat(s);
42 if(!x){
43 x = new TGeoMixture(s,3, 18.1000 );
44 x->SetUniqueID( 25);
45 x->DefineElement(0,183.5,74,0.9827841);
46 x->DefineElement(1,58.6934,28,0.8272327E-02);
47 x->DefineElement(2,63.456,29,0.8943575E-02);
48 }
49
50 s = "CERAMIC";
51 x = (TGeoMixture*)GetMat(s);
52 if(!x){
53 x = new TGeoMixture(s,4, 2.59400 );
54 x->SetUniqueID( 26);
55 x->DefineElement(0,26.98,13,0.21);
56 x->DefineElement(1,28.09,14,0.217);
57 x->DefineElement(2,16.,8,0.558);
58 x->DefineElement(3,1.01,1,0.1500000E-01);
59 }
60
61 s = "G10C";
62 x = (TGeoMixture*)GetMat(s);
63 if(!x){
64 x = new TGeoMixture(s,4, 1.70000 );
65 x->SetUniqueID( 27);
66 x->DefineElement(0,28.09,14,0.53);
67 x->DefineElement(1,16.,8,0.3);
68 x->DefineElement(2,12.01,6,0.15);
69 x->DefineElement(3,1.01,1,0.2000000E-01);
70 }
71
72 TGeoMedium *n=0;
73 s="ALUMINIUM";
74 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03);
75
76 s = "N2 GAS";
77 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,3,19,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.1000000E-02,0.5000000E-03);
78
79 s = "KAPTON";
80 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,6,17,1,1,40,1,100,0.5000000E-01,0.6000012E-04,0.5000000E-03);
81
82 s = "SICALO";
83 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,11,21,1,1,40,1,100,0.5000000E-01,0.3800000E-03,0.5000000E-03);
84
85 s = "TUNGA";
86 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,15,25,-1,1,40,1,100,0.5000000E-01,0.2630000E-02,0.5000000E-03);
87
88 s ="KAOLINITE";
89 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,16,26,-1,1,40,1,100,0.5000000E-01,0.5000000E-04,0.5000000E-03);
90
91 s = "G10C";
92 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,17,27,-1,1,40,1,100,0.5000000E-01,0.1200000E-02,0.5000000E-03);
93
94
95 TGeoVolume *v=0;
96
97 s = "CALB";
98 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),23.2,23.2,10.589);
99 s = "CALS";
100 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),24.2,24.2,10.639);
101 s = "CAPL";
102 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),12.3,12.3,0.19E-01);
103 s = "CANS";
104 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),12.3,12.3,0.7E-01);
105 s = "CASI";
106 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SICALO"),3.904,3.904,0.19E-01);
107 s = "C10C";
108 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("G10C"),12.3,12.3,0.6E-01);
109 s = "CAAB";
110 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("TUNGA"),12.3,12.3,0.1315);
111 s = "CAKA";
112 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("KAOLINITE"),12.3,12.3,0.25E-02);
113 s = "CAKP";
114 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("KAPTON"),12.3,12.3,0.25E-02);
115 s = "CAGL";
116 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SICALO"),12.3,12.3,0.5E-02);
117 //=====//
118 s = "CAST";
119 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SICALO"),3.904,0.122,0.19E-01);
120 //=====//
121
122 GetVol("CALS")->AddNode(GetVol("CALB"),1,new TGeoTranslation(0.,0.,0.5E-01));
123 GetVol("CALB")->AddNode(GetVol("CAPL"),1,new TGeoTranslation(0.1,-0.5E-01,9.876));
124 GetVol("CAPL")->AddNode(GetVol("CASI"),1,new TGeoTranslation(-8.05,-8.05,0.));
125 //====//
126 for (Int_t i=0; i<32; i++){
127 GetVol("CASI")->AddNode(GetVol("CAST"),i, new TGeoTranslation(0.,-3.782+0.244*i,0.));
128 }
129 //====//
130 GetVol("CAPL")->AddNode(GetVol("CASI"),2,new TGeoTranslation(0.,-8.05,0.));
131 GetVol("CAPL")->AddNode(GetVol("CASI"),3,new TGeoTranslation(8.05,-8.05,0.));
132 GetVol("CAPL")->AddNode(GetVol("CASI"),4,new TGeoTranslation(-8.05,0.,0.));
133 GetVol("CAPL")->AddNode(GetVol("CASI"),5,gGeoIdentity);
134 GetVol("CAPL")->AddNode(GetVol("CASI"),6,new TGeoTranslation(8.05,0.,0.));
135 GetVol("CAPL")->AddNode(GetVol("CASI"),7,new TGeoTranslation(-8.05,8.05,0.));
136 GetVol("CAPL")->AddNode(GetVol("CASI"),8,new TGeoTranslation(0.,8.05,0.));
137 GetVol("CAPL")->AddNode(GetVol("CASI"),9,new TGeoTranslation(8.05,8.05,0.));
138 GetVol("CALB")->AddNode(GetVol("CANS"),1,new TGeoTranslation(0.,-0.1,9.787));
139 GetVol("CANS")->AddNode(GetVol("CAKA"),1,new TGeoTranslation(0.,0.,-0.675E-01));
140 GetVol("CANS")->AddNode(GetVol("C10C"),1,new TGeoTranslation(0.,0.,-0.5E-02));
141 GetVol("CANS")->AddNode(GetVol("CAGL"),1,new TGeoTranslation(0.,0.,0.6E-01));
142 GetVol("CANS")->AddNode(GetVol("CAKP"),1,new TGeoTranslation(0.,0.,0.675E-01));
143
144 GetVol("CALB")->AddNode(GetVol("CAAB"),1,new TGeoTranslation(0.,-0.1,9.5855));
145 GetVol("CALB")->AddNode(GetVol("CANS"),2,new TGeoCombiTrans(0.,-0.1,9.384,GetRot("rot3")));
146 GetVol("CALB")->AddNode(GetVol("CAPL"),2,new TGeoCombiTrans(0.5E-01,0.,9.295,GetRot("rot3")));
147 GetVol("CALB")->AddNode(GetVol("CAPL"),3,new TGeoTranslation(-0.1,-0.15,9.067));
148 GetVol("CALB")->AddNode(GetVol("CANS"),3,new TGeoTranslation(0.,-0.1,8.978));
149 GetVol("CALB")->AddNode(GetVol("CAAB"),2,new TGeoTranslation(0.,-0.1,8.7765));
150 GetVol("CALB")->AddNode(GetVol("CANS"),4,new TGeoCombiTrans(0.,-0.1,8.575,GetRot("rot3")));
151 GetVol("CALB")->AddNode(GetVol("CAPL"),4,new TGeoCombiTrans(-0.5E-01,-0.2,8.486,GetRot("rot3")));
152
153 GetVol("CALB")->AddNode(GetVol("CAPL"),5,new TGeoTranslation(0.1,-0.5E-01,8.058));
154 GetVol("CALB")->AddNode(GetVol("CANS"),5,new TGeoTranslation(0.,-0.1,7.969));
155 GetVol("CALB")->AddNode(GetVol("CAAB"),3,new TGeoTranslation(0.,-0.1,7.7675));
156 GetVol("CALB")->AddNode(GetVol("CANS"),6,new TGeoCombiTrans(0.,-0.1,7.566,GetRot("rot3")));
157 GetVol("CALB")->AddNode(GetVol("CAPL"),6,new TGeoCombiTrans(0.5E-01,0.,7.477,GetRot("rot3")));
158 GetVol("CALB")->AddNode(GetVol("CAPL"),7,new TGeoTranslation(-0.1,-0.15,7.249));
159 GetVol("CALB")->AddNode(GetVol("CANS"),7,new TGeoTranslation(0.,-0.1,7.16));
160 GetVol("CALB")->AddNode(GetVol("CAAB"),4,new TGeoTranslation(0.,-0.1,6.9585));
161 GetVol("CALB")->AddNode(GetVol("CANS"),8,new TGeoCombiTrans(0.,-0.1,6.757,GetRot("rot3")));
162 GetVol("CALB")->AddNode(GetVol("CAPL"),8,new TGeoCombiTrans(-0.5E-01,-0.2,6.668,GetRot("rot3")));
163
164 GetVol("CALB")->AddNode(GetVol("CAPL"),9,new TGeoTranslation(0.1,-0.5E-01,6.24));
165 GetVol("CALB")->AddNode(GetVol("CANS"),9,new TGeoTranslation(0.,-0.1,6.151));
166 GetVol("CALB")->AddNode(GetVol("CAAB"),5,new TGeoTranslation(0.,-0.1,5.9495));
167 GetVol("CALB")->AddNode(GetVol("CANS"),10,new TGeoCombiTrans(0.,-0.1,5.748,GetRot("rot3")));
168 GetVol("CALB")->AddNode(GetVol("CAPL"),10,new TGeoCombiTrans(0.5E-01,0.,5.659,GetRot("rot3")));
169 GetVol("CALB")->AddNode(GetVol("CAPL"),11,new TGeoTranslation(-0.1,-0.15,5.431));
170 GetVol("CALB")->AddNode(GetVol("CANS"),11,new TGeoTranslation(0.,-0.1,5.342));
171 GetVol("CALB")->AddNode(GetVol("CAAB"),6,new TGeoTranslation(0.,-0.1,5.1405));
172 GetVol("CALB")->AddNode(GetVol("CANS"),12,new TGeoCombiTrans(0.,-0.1,4.939,GetRot("rot3")));
173 GetVol("CALB")->AddNode(GetVol("CAPL"),12,new TGeoCombiTrans(-0.5E-01,-0.2,4.85,GetRot("rot3")));
174
175 GetVol("CALB")->AddNode(GetVol("CAPL"),13,new TGeoTranslation(0.1,-0.5E-01,4.422));
176 GetVol("CALB")->AddNode(GetVol("CANS"),13,new TGeoTranslation(0.,-0.1,4.333));
177 GetVol("CALB")->AddNode(GetVol("CAAB"),7,new TGeoTranslation(0.,-0.1,4.1315));
178 GetVol("CALB")->AddNode(GetVol("CANS"),14,new TGeoCombiTrans(0.,-0.1,3.93,GetRot("rot3")));
179 GetVol("CALB")->AddNode(GetVol("CAPL"),14,new TGeoCombiTrans(0.5E-01,0.,3.841,GetRot("rot3")));
180 GetVol("CALB")->AddNode(GetVol("CAPL"),15,new TGeoTranslation(-0.1,-0.15,3.613));
181 GetVol("CALB")->AddNode(GetVol("CANS"),15,new TGeoTranslation(0.,-0.1,3.524));
182 GetVol("CALB")->AddNode(GetVol("CAAB"),8,new TGeoTranslation(0.,-0.1,3.3225));
183 GetVol("CALB")->AddNode(GetVol("CANS"),16,new TGeoCombiTrans(0.,-0.1,3.121,GetRot("rot3")));
184 GetVol("CALB")->AddNode(GetVol("CAPL"),16,new TGeoCombiTrans(-0.5E-01,-0.2,3.032,GetRot("rot3")));
185
186 GetVol("CALB")->AddNode(GetVol("CAPL"),17,new TGeoTranslation(0.1,-0.5E-01,2.604));
187 GetVol("CALB")->AddNode(GetVol("CANS"),17,new TGeoTranslation(0.,-0.1,2.515));
188 GetVol("CALB")->AddNode(GetVol("CAAB"),9,new TGeoTranslation(0.,-0.1,2.3135));
189 GetVol("CALB")->AddNode(GetVol("CANS"),18,new TGeoCombiTrans(0.,-0.1,2.112,GetRot("rot3")));
190 GetVol("CALB")->AddNode(GetVol("CAPL"),18,new TGeoCombiTrans(0.5E-01,0.,2.023,GetRot("rot3")));
191 GetVol("CALB")->AddNode(GetVol("CAPL"),19,new TGeoTranslation(-0.1,-0.15,1.795));
192 GetVol("CALB")->AddNode(GetVol("CANS"),19,new TGeoTranslation(0.,-0.1,1.706));
193 GetVol("CALB")->AddNode(GetVol("CAAB"),10,new TGeoTranslation(0.,-0.1,1.5045));
194 GetVol("CALB")->AddNode(GetVol("CANS"),20,new TGeoCombiTrans(0.,-0.1,1.303,GetRot("rot3")));
195 GetVol("CALB")->AddNode(GetVol("CAPL"),20,new TGeoCombiTrans(-0.5E-01,-0.2,1.214,GetRot("rot3")));
196
197 GetVol("CALB")->AddNode(GetVol("CAPL"),21,new TGeoTranslation(0.1,-0.5E-01,0.786));
198 GetVol("CALB")->AddNode(GetVol("CANS"),21,new TGeoTranslation(0.,-0.1,0.697));
199 GetVol("CALB")->AddNode(GetVol("CAAB"),11,new TGeoTranslation(0.,-0.1,0.4955));
200 GetVol("CALB")->AddNode(GetVol("CANS"),22,new TGeoCombiTrans(0.,-0.1,0.294,GetRot("rot3")));
201 GetVol("CALB")->AddNode(GetVol("CAPL"),22,new TGeoCombiTrans(0.5E-01,0.,0.205,GetRot("rot3")));
202 GetVol("CALB")->AddNode(GetVol("CAPL"),23,new TGeoTranslation(-0.1,-0.15,-0.23E-01));
203 GetVol("CALB")->AddNode(GetVol("CANS"),23,new TGeoTranslation(0.,-0.1,-0.112));
204 GetVol("CALB")->AddNode(GetVol("CAAB"),12,new TGeoTranslation(0.,-0.1,-0.3135));
205 GetVol("CALB")->AddNode(GetVol("CANS"),24,new TGeoCombiTrans(0.,-0.1,-0.515,GetRot("rot3")));
206 GetVol("CALB")->AddNode(GetVol("CAPL"),24,new TGeoCombiTrans(-0.5E-01,-0.2,-0.604,GetRot("rot3")));
207
208 GetVol("CALB")->AddNode(GetVol("CAPL"),25,new TGeoTranslation(0.1,-0.5E-01,-1.032));
209 GetVol("CALB")->AddNode(GetVol("CANS"),25,new TGeoTranslation(0.,-0.1,-1.121));
210 GetVol("CALB")->AddNode(GetVol("CAAB"),13,new TGeoTranslation(0.,-0.1,-1.3225));
211 GetVol("CALB")->AddNode(GetVol("CANS"),26,new TGeoCombiTrans(0.,-0.1,-1.524,GetRot("rot3")));
212 GetVol("CALB")->AddNode(GetVol("CAPL"),26,new TGeoCombiTrans(0.5E-01,0.,-1.613,GetRot("rot3")));
213 GetVol("CALB")->AddNode(GetVol("CAPL"),27,new TGeoTranslation(-0.1,-0.15,-1.841));
214 GetVol("CALB")->AddNode(GetVol("CANS"),27,new TGeoTranslation(0.,-0.1,-1.93));
215 GetVol("CALB")->AddNode(GetVol("CAAB"),14,new TGeoTranslation(0.,-0.1,-2.1315));
216 GetVol("CALB")->AddNode(GetVol("CANS"),28,new TGeoCombiTrans(0.,-0.1,-2.333,GetRot("rot3")));
217 GetVol("CALB")->AddNode(GetVol("CAPL"),28,new TGeoCombiTrans(-0.5E-01,-0.2,-2.422,GetRot("rot3")));
218
219 GetVol("CALB")->AddNode(GetVol("CAPL"),29,new TGeoTranslation(0.1,-0.5E-01,-2.85));
220 GetVol("CALB")->AddNode(GetVol("CANS"),29,new TGeoTranslation(0.,-0.1,-2.939));
221 GetVol("CALB")->AddNode(GetVol("CAAB"),15,new TGeoTranslation(0.,-0.1,-3.1405));
222 GetVol("CALB")->AddNode(GetVol("CANS"),30,new TGeoCombiTrans(0.,-0.1,-3.342,GetRot("rot3")));
223 GetVol("CALB")->AddNode(GetVol("CAPL"),30,new TGeoCombiTrans(0.5E-01,0.,-3.431,GetRot("rot3")));
224 GetVol("CALB")->AddNode(GetVol("CAPL"),31,new TGeoTranslation(-0.1,-0.15,-3.659));
225 GetVol("CALB")->AddNode(GetVol("CANS"),31,new TGeoTranslation(0.,-0.1,-3.748));
226 GetVol("CALB")->AddNode(GetVol("CAAB"),16,new TGeoTranslation(0.,-0.1,-3.9495));
227 GetVol("CALB")->AddNode(GetVol("CANS"),32,new TGeoCombiTrans(0.,-0.1,-4.151,GetRot("rot3")));
228 GetVol("CALB")->AddNode(GetVol("CAPL"),32,new TGeoCombiTrans(-0.5E-01,-0.2,-4.24,GetRot("rot3")));
229
230 GetVol("CALB")->AddNode(GetVol("CAPL"),33,new TGeoTranslation(0.1,-0.5E-01,-4.668));
231 GetVol("CALB")->AddNode(GetVol("CANS"),33,new TGeoTranslation(0.,-0.1,-4.757));
232 GetVol("CALB")->AddNode(GetVol("CAAB"),17,new TGeoTranslation(0.,-0.1,-4.9585));
233 GetVol("CALB")->AddNode(GetVol("CANS"),34,new TGeoCombiTrans(0.,-0.1,-5.16,GetRot("rot3")));
234 GetVol("CALB")->AddNode(GetVol("CAPL"),34,new TGeoCombiTrans(0.5E-01,0.,-5.249,GetRot("rot3")));
235 GetVol("CALB")->AddNode(GetVol("CAPL"),35,new TGeoTranslation(-0.1,-0.15,-5.477));
236 GetVol("CALB")->AddNode(GetVol("CANS"),35,new TGeoTranslation(0.,-0.1,-5.566));
237 GetVol("CALB")->AddNode(GetVol("CAAB"),18,new TGeoTranslation(0.,-0.1,-5.7675));
238 GetVol("CALB")->AddNode(GetVol("CANS"),36,new TGeoCombiTrans(0.,-0.1,-5.969,GetRot("rot3")));
239 GetVol("CALB")->AddNode(GetVol("CAPL"),36,new TGeoCombiTrans(-0.5E-01,-0.2,-6.058,GetRot("rot3")));
240
241 GetVol("CALB")->AddNode(GetVol("CAPL"),37,new TGeoTranslation(0.1,-0.5E-01,-6.486));
242 GetVol("CALB")->AddNode(GetVol("CANS"),37,new TGeoTranslation(0.,-0.1,-6.575));
243 GetVol("CALB")->AddNode(GetVol("CAAB"),19,new TGeoTranslation(0.,-0.1,-6.7765));
244 GetVol("CALB")->AddNode(GetVol("CANS"),38,new TGeoCombiTrans(0.,-0.1,-6.978,GetRot("rot3")));
245 GetVol("CALB")->AddNode(GetVol("CAPL"),38,new TGeoCombiTrans(0.5E-01,0.,-7.067,GetRot("rot3")));
246 GetVol("CALB")->AddNode(GetVol("CAPL"),39,new TGeoTranslation(-0.1,-0.15,-7.295));
247 GetVol("CALB")->AddNode(GetVol("CANS"),39,new TGeoTranslation(0.,-0.1,-7.384));
248 GetVol("CALB")->AddNode(GetVol("CAAB"),20,new TGeoTranslation(0.,-0.1,-7.5855));
249 GetVol("CALB")->AddNode(GetVol("CANS"),40,new TGeoCombiTrans(0.,-0.1,-7.787,GetRot("rot3")));
250 GetVol("CALB")->AddNode(GetVol("CAPL"),40,new TGeoCombiTrans(-0.5E-01,-0.2,-7.876,GetRot("rot3")));
251
252 GetVol("CALB")->AddNode(GetVol("CAPL"),41,new TGeoTranslation(0.1,-0.5E-01,-8.304));
253 GetVol("CALB")->AddNode(GetVol("CANS"),41,new TGeoTranslation(0.,-0.1,-8.393));
254 GetVol("CALB")->AddNode(GetVol("CAAB"),21,new TGeoTranslation(0.,-0.1,-8.5945));
255 GetVol("CALB")->AddNode(GetVol("CANS"),42,new TGeoCombiTrans(0.,-0.1,-8.796,GetRot("rot3")));
256 GetVol("CALB")->AddNode(GetVol("CAPL"),42,new TGeoCombiTrans(0.5E-01,0.,-8.885,GetRot("rot3")));
257 GetVol("CALB")->AddNode(GetVol("CAPL"),43,new TGeoTranslation(-0.1,-0.15,-9.113));
258 GetVol("CALB")->AddNode(GetVol("CANS"),43,new TGeoTranslation(0.,-0.1,-9.202));
259 GetVol("CALB")->AddNode(GetVol("CAAB"),22,new TGeoTranslation(0.,-0.1,-9.4035));
260 GetVol("CALB")->AddNode(GetVol("CANS"),44,new TGeoCombiTrans(0.,-0.1,-9.605,GetRot("rot3")));
261 GetVol("CALB")->AddNode(GetVol("CAPL") ,44,new TGeoCombiTrans(-0.5E-01,-0.2,-9.694,GetRot("rot3")));
262
263
264
265 SetMotherProp(GetVol("CALS"),1,new TGeoTranslation(0.,0.,13.05));
266
267 }
268
269
270 void PamVMCDetCalo::DefineCuts(){
271
272 TString s ="ALUMINIUM"; // default GPAMELA parameters
273 if (GetMed(s) && !GetCC(s))
274 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
275 0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
276 s ="N2 GAS";
277 if (GetMed(s) && !GetCC(s))
278 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
279 0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.));
280 s = "IRON"; //default GPAMELA parameters
281 if (GetMed(s) && !GetCC(s))
282 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
283 0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
284 s = "KAPTON"; // default GPAMELA parameters scintillator?
285 if (GetMed(s) && !GetCC(s))
286 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
287 0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
288 s = "SICALO";
289 if (GetMed(s) && !GetCC(s))
290 SetCC(s, new pCutControl(GetMedID(s), 0.00001, 0.00001, 0.00001, 0.00001, 0.001, 0.001, 0.01, 0.00001,
291 0.00001, 0.01, 1., 1., 1., 0., 1., 1., 1., 1., 1., 1., 1., 2.));
292 s = "TUNGA";
293 if (GetMed(s) && !GetCC(s))
294 SetCC(s, new pCutControl(GetMedID(s), 0.00001, 0.00001, 0.00001, 0.00001, 0.001, 0.001, 0.01, 0.00001,
295 0.00001, 0.01, 1., 1., 1., 0., 1., 1., 1., 1., 1., 1., 1., 2.));
296 s = "KAOLINITE";
297 if (GetMed(s) && !GetCC(s))
298 SetCC(s, new pCutControl(GetMedID(s), 0.00001, 0.00001, 0.00001, 0.00001, 0.001, 0.001, 0.01, 0.00001,
299 0.00001, 0.01, 1., 1., 1., 0., 1., 1., 1., 1., 1., 1., 1., 2.));
300 s = "G10C";
301 if (GetMed(s) && !GetCC(s))
302 SetCC(s, new pCutControl(GetMedID(s), 0.00001, 0.00001, 0.00001, 0.00001, 0.001, 0.001, 0.01, 0.00001,
303 0.00001, 0.01, 1., 1., 1., 0., 1., 1., 1., 1., 1., 1., 1., 2.));
304
305 }
306
  ViewVC Help
Powered by ViewVC 1.1.23