1 |
* |
* |
2 |
* $Id: gpdhit.inc,v 3.3 2005/06/21 02:42:09 cafagna Exp $ |
* $Id: gpdhit.inc,v 3.6 2006/04/10 00:07:21 cafagna Exp $ |
3 |
* |
* |
4 |
* $Log: gpdhit.inc,v $ |
* $Log: gpdhit.inc,v $ |
5 |
|
* Revision 3.6 2006/04/10 00:07:21 cafagna |
6 |
|
* The P0 range in ND hit changed |
7 |
|
* |
8 |
|
* Revision 3.5 2006/04/09 23:28:55 cafagna |
9 |
|
* Several new things, among this: ND and CARD |
10 |
|
* |
11 |
|
* Revision 3.4 2005/12/20 12:21:04 cafagna |
12 |
|
* gpnd directory added along with ND files |
13 |
|
* |
14 |
* Revision 3.3 2005/06/21 02:42:09 cafagna |
* Revision 3.3 2005/06/21 02:42:09 cafagna |
15 |
* Major modification to the geometry and to the random number chain |
* Major modification to the geometry and to the random number chain |
16 |
* |
* |
17 |
* Revision 3.2 2003/12/17 11:32:49 pamela |
* Revision 3.2 2003/12/17 11:32:49 pamela CALO SIMULATION COMPLETED: |
18 |
* CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation |
* geometry and special tracking parameters updated and simulation |
19 |
|
* checked by a comparison with the Trieste's standalone Monte Carlo |
20 |
|
* simulation |
21 |
* |
* |
22 |
* Revision 3.1.1.1 2002/07/11 16:01:59 cafagna |
* Revision 3.1.1.1 2002/07/11 16:01:59 cafagna |
23 |
* First GPAMELA release on CVS |
* First GPAMELA release on CVS |
39 |
#if defined(GPAMELA_INTER) |
#if defined(GPAMELA_INTER) |
40 |
* Interactive case. The average position inside the set has been stored |
* Interactive case. The average position inside the set has been stored |
41 |
DATA CHTOF/'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
DATA CHTOF/'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
42 |
+ 'P0 ','TIME','PATH','IPAR','EREL'/ |
+ 'EREL','TIME','PATH','IPAR','P0 '/ |
43 |
C ML: |
C ML: |
44 |
C DATA CHTRD/'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
C DATA CHTRD/'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
45 |
C + 'EREL','PATH','IPAR','P0 '/ |
C + 'EREL','PATH','IPAR','P0 '/ |
53 |
DATA CHSPE/'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
DATA CHSPE/'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
54 |
+ 'EREL','PATH','IPAR','P0 '/ |
+ 'EREL','PATH','IPAR','P0 '/ |
55 |
DATA CHND/ 'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
DATA CHND/ 'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
56 |
+ 'EREL','TIME','IPAR','P0 '/ |
+ 'EREL','TIME','PATH','IPAR','P0 '/ |
57 |
|
DATA CHCARD/'XAVE','YAVE','ZAVE','XOUT','YOUT','ZOUT', |
58 |
|
+ 'EREL','TIME','PATH','IPAR','P0 '/ |
59 |
#endif |
#endif |
60 |
#if !defined(GPAMELA_INTER) |
#if !defined(GPAMELA_INTER) |
61 |
DATA CHTOF/'XIN ','YIN ','ZIN ','XOUT','YOUT','ZOUT', |
DATA CHTOF/'XIN ','YIN ','ZIN ','XOUT','YOUT','ZOUT', |
62 |
+ 'P0 ','TIME','PATH','IPAR','EREL'/ |
+ 'EREL','TIME','PATH','IPAR','P0 '/ |
63 |
C ML: |
C ML: |
64 |
C DATA CHTRD/'XIN ','YIN ','ZIN ','XOUT','YOUT','ZOUT', |
C DATA CHTRD/'XIN ','YIN ','ZIN ','XOUT','YOUT','ZOUT', |
65 |
C + 'EREL','PATH','IPAR','P0 '/ |
C + 'EREL','PATH','IPAR','P0 '/ |
74 |
+ 'EREL','PATH','IPAR','P0 '/ |
+ 'EREL','PATH','IPAR','P0 '/ |
75 |
C ML: |
C ML: |
76 |
DATA CHND/ 'XIN','YIN','ZIN','XOUT','YOUT','ZOUT', |
DATA CHND/ 'XIN','YIN','ZIN','XOUT','YOUT','ZOUT', |
77 |
+ 'EREL','TIME','IPAR','P0 '/ |
+ 'EREL','TIME','PATH','IPAR','P0 '/ |
78 |
|
DATA CHCARD/'XIN ','YIN ','ZIN ','XOUT','YOUT','ZOUT', |
79 |
|
+ 'EREL','TIME','PATH','IPAR','P0 '/ |
80 |
|
|
81 |
C +SELF. |
C +SELF. |
82 |
#endif |
#endif |
83 |
#if defined(GPAMELA_GARFIELD) |
#if defined(GPAMELA_GARFIELD) |
89 |
* No packing required for the last CAL hit to use the GSCHIT routine |
* No packing required for the last CAL hit to use the GSCHIT routine |
90 |
DATA NBHTOF/NHTOF*32/,NBHTRD/NHTRD*32/, |
DATA NBHTOF/NHTOF*32/,NBHTRD/NHTRD*32/, |
91 |
+ NBHCAT/NHCAT*32/,NBHCAS/NHCAS*32/,NBHSPE/NHSPE*32/, |
+ NBHCAT/NHCAT*32/,NBHCAS/NHCAS*32/,NBHSPE/NHSPE*32/, |
92 |
+ NBHCAL/NHCAL*32/,NBHND/NHND*32/ |
+ NBHCAL/NHCAL*32/,NBHND/NHND*32/,NBHCARD/NHCARD*32/ |
93 |
* The coordinates have to be always positive. We need microns for XYZ and eV. |
* The coordinates have to be always positive. We need microns for XYZ and eV. |
94 |
* Time is in second but we need psec. |
* Time is in second but we need psec. |
95 |
DATA OHTOF/6*100.,0.,0.,0.,0.,0./, |
DATA OHTOF/6*100.,0.,0.,0.,0.,0./, |
96 |
*ML: |
*ML: |
97 |
* + FHTOF/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
* + FHTOF/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
98 |
+ FHTOF/6*10000.,1.E3,1.E12,10000.,1.,1.E9/ |
+ FHTOF/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
99 |
*END ML. |
*END ML. |
100 |
C ML: |
C ML: |
101 |
C DATA OHTRD/6*100.,0.,0.,0.,0./, |
C DATA OHTRD/6*100.,0.,0.,0.,0./, |
113 |
+ FHCAT/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
+ FHCAT/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
114 |
DATA OHCAS/6*100.,0.,0.,0.,0.,0./, |
DATA OHCAS/6*100.,0.,0.,0.,0.,0./, |
115 |
+ FHCAS/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
+ FHCAS/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
116 |
|
DATA OHCARD/6*100.,0.,0.,0.,0.,0./, |
117 |
|
+ FHCARD/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
118 |
DATA OHSPE/6*100.,0.,0.,0.,0./, |
DATA OHSPE/6*100.,0.,0.,0.,0./, |
119 |
+ FHSPE/6*10000.,1.E9,10000.,1.,1.E3/ |
+ FHSPE/6*10000.,1.E9,10000.,1.,1.E3/ |
120 |
*EM: |
*EM: |
126 |
DATA OHCALI/0.,0./,FHCALI/1.,1.E9/ |
DATA OHCALI/0.,0./,FHCALI/1.,1.E9/ |
127 |
*end ml. |
*end ml. |
128 |
*END EM. |
*END EM. |
129 |
DATA OHND/6*100.,0.,0.,0.,0./, |
DATA OHND/6*100.,0.,0.,0.,0.,0./, |
130 |
+ FHND/6*10000.,1.E9,1.E12,1.,1.E9/ |
+ FHND/6*10000.,1.E9,1.E12,10000.,1.,1.E3/ |
131 |
|
|
132 |
|
|
133 |
|
|