/[PAMELA software]/PamVMC_update/trk/include/f77/gpdgeo.inc
ViewVC logotype

Contents of /PamVMC_update/trk/include/f77/gpdgeo.inc

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1.1.1 - (show annotations) (download) (vendor branch)
Tue Oct 15 15:51:28 2013 UTC (11 years, 1 month ago) by formato
Branch: MAIN, rel
CVS Tags: reltag, HEAD
Changes since 1.1: +0 -0 lines
PamVMC update

1 *
2 * $Id: gpdgeo.inc,v 1.5 2009-06-12 18:39:52 pam-rm2 Exp $
3 *
4 * $Log: gpdgeo.inc,v $
5 * Revision 1.5 2009-06-12 18:39:52 pam-rm2
6 * - Introduced user-defined names of output files and random seeds number.
7 * Users can do it use options of PamVMCApplication constructor:
8 * PamVMCApplication(const char* name, const char *title, const char*
9 * filename="pamtest", Int_t seed=0).
10 * The Random object that I use is TRandom3 object which has astronomical
11 * large period (in case of default initialization 0). All random generators
12 * in the code use this object by calling of gRandom singleton which keeps
13 * it.
14 *
15 * - Corrected TOF digitization routine. No problems with TDC hits due to
16 * hadronic interactions anymore.
17 *
18 * - Some small changes was done to compile code under Root 5.23. +
19 * geant4_vmc v. 2.6 without any warnings
20 *
21 * - Some classes of PamG4RunConfiguartion was changed for geant4_vmc v.
22 * 2.6.Some obsolete classes was deleted as soon as developers implemented
23 * regions.
24 *
25 * - Navigation was changed from "geomRootToGeant4" to "geomRoot", because on
26 * VMC web page written that as soon as Geant4 has no option ONLY/MANY
27 * translation of overlapped geometry to Geant4 through VGM could be wrong.
28 * I'd like to stay with Root navigation:
29 * http://root.cern.ch/root/vmc/Geant4VMC.html. This should be default
30 * option.
31 *
32 * - New Tracker digitization routine written by Sergio was implemented
33 *
34 * - PamVMC again became compatible with geant4_vmc v.2.5 and ROOT 5.20.
35 * The problem was that ROOT developers introduced in TVirtualMC class a new
36 * method SetMagField and new base class:TVirtualMagField from which
37 * user-defined classes shoukd be derived
38 *
39 * Revision 1.1 2009-02-19 17:46:24 nikolas
40 * Cleaning up before releasing
41 *
42 * Revision 3.23 2008/01/29 18:25:16 pamela
43 * Review of the shell and TOF geometries. New materials around the shell
44 * added to simulate the insulation.
45 *
46 * Revision 3.22 2006/12/17 14:53:11 cafagna
47 * Full review of the CARD geometry. S1 table added
48 *
49 * Revision 3.21 2006/11/30 12:29:14 cafagna
50 * CARD geometry updated. New S1 "table", legs included, added. Top magnetic screen has been added as well.
51 *
52 * Revision 3.20 2006/11/28 10:26:14 pam-ba
53 * S3 positioning completed
54 *
55 * Revision 3.19 2006/11/16 18:45:29 pam-ba
56 * Simulated an aluminum container for S4
57 *
58 * Revision 3.18 2006/11/16 12:01:13 pam-ba
59 * S4 dimensions corrected.
60 *
61 * Revision 3.17 2006/11/10 11:39:34 pam-ba
62 * S2 and S1 z-positions corrected, He3 and plystyrene mixture added, Top Plate geometry simulated and titanium mixture added.
63 *
64 * Revision 3.16 2006/10/12 11:11:21 pam-ba
65 * ND geometry updated.
66 *
67 * Revision 3.15 2006/06/30 15:38:16 pam-ba
68 * S22 and S12 heights positioned in GPAMELA at the nominal heights in PAMELA (see document: Main geometrical parameters of the PAMELA sub-detectors, 20 December 2005)
69 *
70 * Revision 3.14 2006/05/18 10:52:31 pam-ba
71 * TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
72 *
73 * Revision 3.13 2006/05/09 12:01:16 cafagna
74 * Calo geometry updated
75 *
76 * Revision 3.12 2006/05/02 12:30:02 bottai
77 * change in trcp dimension in order to not overlap with its mother
78 *
79 * Revision 3.11 2006/05/02 11:32:14 pam-ba
80 * TOF geometry and position updated and a new material, the mylar, added
81 *
82 * Revision 3.10 2006/04/09 23:28:54 cafagna
83 * Several new things, among this: ND and CARD
84 *
85 * Revision 3.9 2005/12/20 12:21:04 cafagna
86 * gpnd directory added along with ND files
87 *
88 * Revision 3.8 2005/12/16 10:20:18 cafagna
89 * Bug fixed in the new update
90 *
91 * Revision 3.7 2005/12/16 09:13:55 cafagna
92 * New small valume added to the tracker frame
93 *
94 * Revision 3.6 2005/12/13 10:31:29 pam-ba
95 * added a new volume, TPGI, in the spectrometer
96 *
97 * Revision 3.5 2005/12/05 12:15:20 pam-ba
98 * new spectrometer geometry and internal magnetic field
99 *
100 * Revision 3.4 2005/07/25 11:53:21 cafagna
101 * Several updates. See history for details
102 *
103 * Revision 3.3 2005/06/21 02:42:08 cafagna
104 * Major modification to the geometry and to the random number chain
105 *
106 * Revision 3.2 2002/12/05 10:17:42 pamela
107 * Update CAS and CALO geometries and positions. Makefile updated as well
108 *
109 * Revision 3.1.1.1 2002/07/11 16:01:59 cafagna
110 * First GPAMELA release on CVS
111 *
112 *
113 *
114 * gpdgeo.inc
115 *
116 * 16/07/2004 16.44.06 by Alexey Bakaldin
117 * 17/10/2002 22.31.45 by Jens Lund
118 *CMZ : 3.00/00 14/11/2000 09.13.05 by Emiliano Mocchiutti
119 *CMZ : 2.02/00 12/10/2000 18.43.21 by Francesco Cafagna
120 *CMZ : 2.01/00 05/04/2000 14.35.17 by Marialuigia Ambriola
121 *CMZ : 2.00/00 03/03/2000 15.39.05 by Francesco Cafagna
122 *CMZ : 1.02/00 04/02/2000 18.40.04 by Francesco Cafagna
123 *CMZ : 1.01/00 06/05/96 11.40.48 by Francesco Cafagna
124 *CMZ : 1.00/02 02/04/96 16.44.16 by Francesco Cafagna
125 *-- Author : Francesco Cafagna 05/12/95
126 DATA PAME/0.,53.,130./
127 * AB:
128 C DATA SHEL/69.8,70.,75./
129 cml DATA SHEL/40.6,52.35,52.5,43.62253,43.77253/
130 cml DATA MSHE/52.25,52.5,20./
131 cml DATA BSPH/52.25,52.5,90.,180.,0.,360./
132 cml:
133 * ab mod DATA SHEL/38.25,52.5,52.65,43.77,43.92/
134 * ab mod DATA MSHE/52.5,52.75,20./
135 * ab mod DATA BSPH/52.5,52.75,90.,180.,0.,360./
136 c end ml
137 * END: AB.
138 C DATA TFLA/126.8,127.,0.,33.45,0.,360./
139 DATA TFLA/0.,43.77253,0.125/
140 * AB:
141 C DATA TSPH/0., 150.,0.,0.,0.,360./
142 cml DATA TSPH/239.83,240.,0.,0.,0.,360./
143 cml:
144 c ab mod DATA TSPH/240.,240.17,0.,10.5062,0.,360./
145 * pressurized container upgrated by ab:
146 DATA SHEL/38.8717,52.5,52.651,43.6423,43.7932/
147 DATA MSHE/52.5,52.75,19.95/
148 DATA BSPH/52.5,52.75,90.,180.,0.,360./
149 DATA TSPH/240.,240.17,0.,10.5062,0.,360./
150 DATA COV1/240.17,240.1715,0.,10.5062,0.,360./
151 DATA COV2/240.1715,240.2165,0.,10.5062,0.,360./
152 DATA COV3/240.2165,240.2315,0.,10.5062,0.,360./
153 DATA ZTSPH/-130.716/
154 DATA ZSHEL/66.5559/
155 DATA ZMSHE/7.7342/
156 DATA ZBSPH/-12.2158/
157 c end ml.
158 * END: AB.
159 DATA GAPTOP/2./
160 c ml:26/10/04
161 c DATA S1/20.4,16.5,.375/
162 c DATA S2/20.4,16.5,.375/
163 c DATA S1A/9.,7.5,.375/
164 c DATA S2A/9.,7.5,.375/
165 DATA S3/0.,0.,0./
166 c DATA S4/18.,16.,.5/
167 DATA S4/24.1,24.1,0.5/
168 DATA S4AL/24.2,24.2,0.6/
169 DATA S11M/0.,0.,0./
170 DATA S12M/0.,0.,0./
171 DATA S21M/0.,0.,0./
172 DATA S22M/0.,0.,0./
173 DATA S31M/0.,0.,0./
174 DATA S32M/0.,0.,0./
175 DATA POL1/0.,0.,0./
176 DATA POL2/0.,0.,0./
177 DATA POLY/0.,0.,0./
178 DATA S11Y/2.55,16.5,0.35/
179 DATA S12X/20.4,2.75,0.35/
180 DATA S21X/9.,3.75,0.25/
181 DATA S22Y/4.5,7.5,0.25/
182 DATA S31Y/3.,7.5,0.35/
183 DATA S32X/9.,2.5,0.35/
184 c DATA HMYLTOF/0.0025/
185 DATA HMYLTOF/0.0023/
186 DATA HGLUTOF/0.01/
187 DATA HPVCMAX/0.9/
188 DATA HPVCMIN/0.4/
189 DATA HALTOF/0.033/
190 DATA HAIRS3/0.15/
191 DATA HTABLE/24.7/
192 DATA DTABLEAL/1.7/
193 DATA XTPLA/23.2/
194 DATA YTPLA/22.1/
195 c ml: 28/10/04:
196 c dimensions of top plate. At the moment, it is a reduced version,
197 c i.e. a gap of gas with a surface equal to S21.
198 c DATA TPLA/9.,7.5,.25/
199 c dimensions of top plate from CAD draws:
200 c DATA TPLA/21.9,21.1,1.6/
201 DATA TPLA/28.,27.165,1.6/
202 C DATA TPCA/18.1,18.2,0.6/
203 C DATA TPS2/10.,8.1,0.75/
204 DATA TPCV/8.5,7.,0.25/
205 DATA TPTU/2.55,1.95,0.6/
206 DATA TPTM/3.45,3.3,0.75/
207 DATA TPTL/21.9,21.1,.25/
208 c end ml.
209 c DATA SCTIC/0.025/
210 c ml:26/10/04:
211 DATA NS11Y/8/
212 DATA NS12X/6/
213 DATA NS21X/2/
214 DATA NS22Y/2/
215 DATA NS31Y/3/
216 DATA NS32X/3/
217 c end ml.
218 C ml 28/03/06: official z coordinates in the PAMELA reference system with
219 c origin in the centre of the central plane of the central magnetic module:
220 c DATA ZPAMS31,ZPAMS32/-23.49,-24.34/
221 DATA DTPS3/0.845/
222 c DATA ZS11,ZS12,ZS21,ZS22,ZS31,ZS32,ZS4
223 c ml 22/10/04:
224 c + /107.8,106.8,74.9,74.1,27.1,26.3,2.3/
225 c ml 28/10/04: introducing the top plate and calculating the zpos of S21, S22
226 c starting from the zpos of CATA:
227 c + /106.66,105.66,74.9,74.1,27.1,26.3,2.3/
228 c + /106.66,105.66,74.875,74.125,27.1,26.3,2.3/
229 c DATA ZTPLA/73.5/
230 C HORMG is the distance beteween the bottom the tracker and the center of
231 c the third magnetic module.
232 DATA HORMG/22.57/
233 c end ml 28/10/04
234 c end ml 22/10/04:
235 C CAT data
236 *EM:
237 *JeL:
238 DATA CATA/28.0,0.6,27.165/
239 *END: JeL.
240 c DATA CATB/22.05,21.05,0.55/
241 DATA CATB/22.2,21.3,0.5/
242 c DATA CATL/15.75,6.675,0.55/
243 c DATA CATT/6.3,14.375,0.55/
244 DATA CATL/2.6,14.9,0.4/
245 DATA CATT/8.45,21.2,0.4/
246 DATA CATH/9.4,7.7,0.5/
247 DATA CATF/2.6,3.15,0.4/
248 DATA CATO/9.3,0.05,0.5/
249 DATA CATP/0.05,7.7,0.5/
250 c DATA ZCAT/75.85/
251 C CAS data
252 *EM:
253 *JeL:
254 DATA CASA/21.45,0.6,20.55/
255 *END: JeL.
256 C DATA CASX/15.25,0.55,23./
257 C DATA CASY/0.55,14.75,23./
258 DATA CSSX/15.7,0.5,21.6/
259 DATA CSSY/0.5,15.7,21.6/
260 DATA CASX/15.6,0.4,21.5/
261 DATA CASY/0.4,15.6,21.5/
262 DATA CASXD/2.4/
263 DATA CASYD/2.2/
264 C DATA ZCAS/51.25/
265 *END: EM.
266 C CARD data
267 DATA CARDB/27.,22.8,12.35/,
268 + CAR/18.15,4.8,1.8,3.6,5.8,0.4,76.,19.45/,
269 + ATZ/1.1/, ! aluminium box thickness
270 + ATY/0.15/, ! aluminium box thickness
271 + PT/0.9/, ! plastic box thickness
272 + CARA/31.98273,17.45,20.,19.123,16.80,47.358,30.4/,
273 + CARB/37.18273,20.95,17.,18.67,16.014,47.661,35.6/,
274 + LEGB/2.15,3.575,11.85/,
275 + LEG1/0.65,0.,0.,0.15,6.85,6.85,0.,0.15,10.85,10.85,0./,
276 + LEG2/0.65,1.65,11.35/,
277 + LEG3/0.1,1.3,0.325/,
278 + LEG4/0.65,1.625,11.85/,
279 + LEG5/1.125,2.575,11.85/,
280 + LEG6/0.1,0.55,11.2/,
281 + LEG7/1.1,0.75,3.25/,
282 + LEG8/1.1,0.75,4.35/,
283 + LEG9/1.1,0.75,3.25/,
284 + LEGP/1.95,0.65,1.125,0.375,11.85,3.8,0.7,0.6/,
285 + TPLT/27.,22.8,0.5/,
286 + TOPC/20.4,16.,0.5/,
287 + TOPT/0.5,0.,0.,2.25,2.25,0.,-26.565,2.25,2.25,0.,-26.565/,
288 + TOHO/0.,0.7,0.5/,
289 + TH11/2.8,3.2,0.5/,
290 + TH12/2.8,5.3,0.5/,
291 + TH21/1.7,3.,0.5/,
292 + TH22/4.9,3.,0.5/,
293 + TOPP/23.6,14.4,5.5,17.2,19.4,10.2,4.5,6.2,4.1/,
294 + MGSC/16.,15.4,0.12456/,
295 + MGSH/9.8,8.,0.12456/
296 C TRD data
297 c ml: 11/11/04
298 c DATA TRFR/19,17.5,.2/
299 DATA TRFR/21.1,19.0,.275/
300 DATA TRFD/19.75,13.9,.2425/
301 DATA TRFU/19.4,13.55,.2425/
302 DATA TRFM/21.1,19.0,0.1375/
303 DATA TRFL/20.175,14.325,1.13/
304 c DATA TRFI/15.5,14.,.2/
305 DATA TRR0/21.1,38.5,.25/
306 DATA TRI0/10.65,8.85,.15/
307 DATA TRRF/16.0,15.4,0.1/
308 DATA TRRI/9.8,8.0,0.1/
309 c end ml.
310 DATA TRBS/3.3,18.3,1./
311 *ml: 10/11/66:
312 DATA TRAL/3.3,2.15,1./
313 DATA TRR2/19.025,13.175,1.51/
314 *end ml.
315 DATA TRDS/24.,24.5,13.5/
316 DATA TRDT/3.,3.,.25/
317 DATA TRDB/17.7,18.5,0/
318 DATA TRSO/0.197,0.2,14./
319 C DATA TRSI/0.,0.197,14./
320 DATA TRSI/0.,0.2,14./
321 *ml: 12/11/04:
322 c DATA TRRA/13.2,14.,0/
323 DATA TRRA/18.25,12.4,0.76/
324 *end ml.
325 *ml: 17/11/04:
326 c DATA TRAN/1.75,1.75,1.0/
327 c DATA TRAI/1.75,1.65,0.9/
328 DATA TRAN/1.8,1.8,1.0/
329 DATA TRAI/1.0,1.0,1.0/
330 *end ml.
331 DATA ZTRD/90.55/
332 DATA GAPTRD/1.5/
333 C Spectrometer data, mod by Caf on 20 Jul 2005
334 DATA SPEB/15.,14.4,22.615/
335 c ml: 31/10/05:
336 c DATA MGFR/14.,13.4,4.0/
337 DATA MGFR/15.,14.4,4.45/
338 c end ml.
339 c ml: 26/10/04:
340 c DATA MGPL/12.,11.4,4.05/
341 DATA MGPL/12.,11.4,4.0/
342 DATA MGPA/8.1,6.6,4.0/
343 c DATA MGPI/8.05,6.55,4.05/
344 c ml: 31/10/05:
345 c DATA MGPI/8.1,6.6,4.0/
346 C DATA MGPI/8.07,6.57,4./
347 C DATA MGPI/8.07,6.57,4.03/
348 DATA MGPI/8.07,6.57,4./
349 C end ml.
350 c ml:31/10/05:
351 c new volumes:
352 C DATA TPGA/10.5,10.65,0.055/
353 DATA TPGA/10.29,10.65,0.04/
354 DATA TPGU/10.49,12.15,0.01/
355 DATA TPGD/12.,11.4,0.005/
356 DATA TPGI/8.07,6.57,0.015/
357 c end ml.
358 c ml:17/12/04:
359 c DATA TRPB/12.,11.4,0.4/
360 c$$$ mod by Caf on 20 Jul 2005
361 c$$$ DATA TRPB/10.49,12.15,0.365/
362 c ml 31/10/05:
363 c DATA TRPB/15.,14.4,0.365/
364 DATA TRPB/10.49,12.15,0.365/
365 c end ml.
366 ccc ml: 20/04/05:
367 cc DATA TRPL/0.,7.,0.4/
368 C # DATA TPAI/8.0995,8.5,0.055/
369 C # DATA TPAS/8.0995,9.75,0.06/
370 C # Mod by CAF on 15 Jul 2005
371 c DATA TPAI/8.0995,1.25,0.055/
372 c DATA TPAS/8.0995,9.75,0.365/
373 c ml: 15/12/05:
374 c ml 19/9/05:
375 c DATA TPAS/8.144,9.9,0.365/
376 DATA TPAS/8.2,9.9,0.365/
377 DATA TPAI/8.144,1.25,0.05/
378 C # End mod
379 c DATA TRSL/0.,7.,0./
380 C # DATA TRSL/2.6665,7.,0.015/
381 DATA TRSL/2.6665,3.5,0.015/
382 c DATA TSPA/2.666,3.5,0.015/
383 DATA TSPA/2.596,3.413,0.015/
384 DATA THBP/2.6665,2.75,0.015/
385 c DATA TRCP/0.005,7.,0.3/
386 c DATA TRCP/0.025,10.125,0.25/
387 c ml: 19/9/05: (based on the measure of the final version of the tracker plane)
388 c DATA TRCP/0.035,10.125,0.25/
389 c change TRCP in order to not overla with mother
390 DATA TRCP/0.035,9.9,0.25/
391 DATA TBAL/1.49,2.25,0.365/
392 c ml:16/11/05: distance between the bottom of the tracker and the bottom of
393 c the base plate.
394 DATA ZSPEBP/2.97/
395 c DATA ZSPEC/50.6/
396 C ml: 15/12/05: xglue is the dimension along x of the glue between each carbon
397 c fiber barrel and the TRSL volume
398 DATA XGLUE/0.0015/
399 C # DATA NTRSL/3/
400 DATA NTRSL/6/
401 DATA NTHBP/3/
402 c end ml.
403 C Calorimeter data
404 *EM:
405 * DATA CASI/3.9,3.9,0.019/
406 DATA CASI/3.904,3.904,0.019/
407 *F.V.&ML:
408 C DATA CG10/12.1,12.1,0.0875/
409 C DATA C10C/12.1,12.1,0.08/
410 C DATA CAAB/12.1,12.1,0.133/
411 C DATA CAPL/12.1,12.1,0./
412 C DATA CAAD/12.1,12.1,0.133/
413 C DATA CAPD/12.1,12.1,0./
414 C DATA CALB/19.7,19.7,10.2/
415 *MA
416 DATA DZM0/0.19/, DZST/0.390/
417 DATA CG10/12.3,12.3,0.0875/
418 C # DATA C10C/12.3,12.3,0.08/
419 C # DATA CAAB/12.3,12.3,0.1315/
420 DATA C10C/12.3,12.3,0.06/
421 DATA CAAB/12.3,12.3,0.1315/
422 * end MA
423 DATA CAPL/12.3,12.3,0./
424 DATA CANS/12.3,12.3,0./
425 DATA CAAD/12.3,12.3,0.133/
426 DATA CAPD/12.3,12.3,0./
427 *ML:
428 * DATA CALB/19.7,19.7,10.45/
429 DATA CALB/23.2,23.2,10.589/
430 *END ML.
431 *END F.V.&ML.
432 *ML:
433 * DATA CALS/20.,20.,10.5/
434 DATA CALS/24.2,24.2,10.639/
435 DATA ZCAL/13.05/
436 C with ZCAL=13.03 the distance between the bottom of the magnet and
437 C the top of the top silicon strips is about 44.24 mm, as it has to be
438 C if 37 mm is the height of the base plate and 7.24 mm is the distance
439 C between the bottom of the end of the base plate and the top silicium
440 C strips. Because of the difference is small and the height of the base
441 C plate could be approximate, I left ZCAL=13.05.
442 C DATA ZCAL/13.03/
443 *END ML.
444 * DATA CALSTWID/0.36/
445 DATA CALSTWID/0.244/
446 *
447 * NOTE: IF you change NCAPL, remember to modify the NCPL, NCSI, NOCAPL,
448 * NOCASI array dimentions in $GPKEY, $GPDKEY common and data block.
449 * The corresponding FFKEY call in GPFFR must be changed as well.
450 * The GPDAT loop on these arrays should be checked too.
451 *
452 DATA NCAPL/22/
453 DATA NCAPLD/1/
454 DATA DCASIX/0.242/
455 DATA DCASIY/0.242/
456 DATA DCASIZ/0.694/
457 DATA NCASTR/32/
458 DATA CAKA/12.3,12.3,0.0025/
459 DATA CAGL/12.3,12.3,0.005/
460 DATA CAKP/12.3,12.3,0.0025/
461 C
462 ***** Neutron detector
463 DATA NDBI/21.2,21.2,6.875/,NDTU/0.,0.925,10./,NDTI/0.,0.92,10./,
464 + NDPB/21.2,21.2,1.0/,NDCD/21.25,21.25,6.9/
465 DATA NDBX/21.45,21.45,7./
466 DATA NDCO/27.5,30.,0.5/
467 DATA NDBO/27.5,30.,7.5/
468 DATA NDBS/15.,4.275,7./
469 DATA NDSS/14.8,4.075,6.8/
470 DATA NDCI/21.2,11.3,0.325/
471 DATA NDCM/21.2,4.25,0.325/
472 DATA NDCE/14.6,2.875,0.325/
473 DATA NDSI/1.,7.26,0.325/
474 DATA NDSM/1.,5.8,0.325/
475 DATA NDSE/4.2,1.05,0.325/
476 DATA TNDBX/0.2/
477 DATA TNDCD/0.05/
478 DATA GNDTU/0.5/

  ViewVC Help
Powered by ViewVC 1.1.23