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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (hide annotations) (download)
Tue Oct 15 15:51:28 2013 UTC (11 years, 2 months ago) by formato
Branch point for: MAIN, rel
Initial revision

1 formato 1.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