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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.5 - (hide 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 nikolas 1.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