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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.1.1.1 - (hide annotations) (download) (vendor branch)
Tue Oct 15 15:51:13 2013 UTC (12 years ago) by formato
Branch: MAIN, rel
CVS Tags: reltag, HEAD
Changes since 1.1: +0 -0 lines 
PamVMC update
1 formato 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