1 |
|
# Revision 3.21 2006/12/05 18:45:29 bruno |
2 |
|
# Update of the CARD geometry |
3 |
|
# Simulated also legs, table and magnetic screen |
4 |
# |
# |
5 |
# $Id: v_100.txt,v 3.16 2006/10/12 11:11:21 pam-ba Exp $ |
# $Id: v_100.txt,v 3.20 2006/11/16 18:45:29 pam-ba Exp $ |
6 |
# |
# |
7 |
# $Log: v_100.txt,v $ |
# $Log: v_100.txt,v $ |
8 |
|
# Revision 3.20 2006/11/16 18:45:29 pam-ba |
9 |
|
# Simulated an aluminum container for S4 |
10 |
|
# |
11 |
|
# Revision 3.19 2006/11/16 12:01:13 pam-ba |
12 |
|
# S4 dimensions corrected. |
13 |
|
# |
14 |
|
# Revision 3.18 2006/11/10 11:39:35 pam-ba |
15 |
|
# S2 and S1 z-positions corrected, He3 and plystyrene mixture added, Top Plate geometry simulated and titanium mixture added. |
16 |
|
# |
17 |
|
# Revision 3.17 2006/10/13 16:36:59 pam-ba |
18 |
|
# Added a new material, the cadmium, for ND |
19 |
|
# |
20 |
# Revision 3.16 2006/10/12 11:11:21 pam-ba |
# Revision 3.16 2006/10/12 11:11:21 pam-ba |
21 |
# ND geometry updated. |
# ND geometry updated. |
22 |
# |
# |
77 |
#CMZ : 1.00/03 30/04/96 12.23.59 by Francesco Cafagna |
#CMZ : 1.00/03 30/04/96 12.23.59 by Francesco Cafagna |
78 |
#CMZ : 1.00/02 05/04/96 15.31.25 by Francesco Cafagna |
#CMZ : 1.00/02 05/04/96 15.31.25 by Francesco Cafagna |
79 |
#CMZ : 1.00/01 28/11/95 18.51.23 by Francesco Cafagna |
#CMZ : 1.00/01 28/11/95 18.51.23 by Francesco Cafagna |
80 |
#-- Author : Francesco Cafagna 28/11/95 |
#-- Author : Francesco Cafagna 28/11/95 |
81 |
|
|
82 |
|
December 2006, Bari |
83 |
|
|
84 |
|
CARD GEOMETRY: |
85 |
|
Positions, dimensions and quotes have been updated. Indeed, the aluminium table |
86 |
|
under S1, legs included, and the magnetic screen between the CARD and S2 have |
87 |
|
been added. |
88 |
|
|
89 |
|
#------------------------------------------------------------------------- |
90 |
|
|
91 |
|
November 2006, Bari |
92 |
|
|
93 |
|
|
94 |
|
S4 GEOMETRY UPDATED: |
95 |
|
The x and y dimensions have been corrected because they were wrong. |
96 |
|
A container of aluminum has been added having walls of thickness equal |
97 |
|
to 0.1 cm. |
98 |
|
|
99 |
|
TOF POSITIONS DEFINITELY UPDATED: |
100 |
|
S2 and S1 have been positioned again, after the simulation of the top plate. |
101 |
|
Before now, the positions were put ad hoc based on the positions given by |
102 |
|
Sergio Ricciarini in the document 'Main geometrical parameters of the PAMELA |
103 |
|
sub-detectors' released by O. Adriani, L. Bonechi, E. Mocchiutti and |
104 |
|
S. Ricciarini on the 20th of December 2005.. These positions were |
105 |
|
lightly approximated. |
106 |
|
|
107 |
|
ND GEOMETRY COMPLETED. |
108 |
|
The He3 and the polyethylene (CH2) have been added to fill |
109 |
|
the volumes NDTI and NDPB in the neutron detector. Leonov gave us the |
110 |
|
details of these materials. |
111 |
|
|
112 |
|
TOP PLATE GEOMETRY ADDED. |
113 |
|
The top plate geometry has been simulated as a N2 box that has the same |
114 |
|
dimensions of CATA. Inside it, at the right positions, there are S2 and CATA. |
115 |
|
The real top plate is a rectangular window having a thickness of 0.5 cm |
116 |
|
(volume TPTL emptied by the volume TPCV). This volume is filled with a mixture |
117 |
|
of titanium (90%), aluminum (6%) and vanadium (4%). Other parts of titanium |
118 |
|
are the four volumes TPTU, which are positioned at the corners of TPTL at the |
119 |
|
same height of CATA, and the four volumes TPTM, which are positioned at the |
120 |
|
corners of TPTL, at the same height of S2. |
121 |
|
The geometry of the top plate is based on a simplified version of the CAD |
122 |
|
drawings. |
123 |
|
|
124 |
October 2006, Bari |
October 2006, Bari |
125 |
|
|
126 |
The geometry of the neutron detector has been updated. Some dimensions of some |
The geometry of the neutron detector has been updated. Some dimensions of some |
127 |
volumes have been corrected while the aluminum cover and the aluminum boxes to |
volumes have been corrected and the aluminum cover and the aluminum boxes to |
128 |
put the cables have been added. A new material, the cadmium, has been also |
put the cables have been added. A new material, the cadmium, has been also |
129 |
defined. |
defined. |
130 |
|
|
134 |
missing in the subroutine gpgeo.F, has been added; the meaning of the SPHE |
missing in the subroutine gpgeo.F, has been added; the meaning of the SPHE |
135 |
data card has been changed. Before the correction the data card: |
data card has been changed. Before the correction the data card: |
136 |
NDET 'SPHE' was used to delete the spherical top shell to substitute it with |
NDET 'SPHE' was used to delete the spherical top shell to substitute it with |
137 |
a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA. |
a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA. |
138 |
|
|
139 |
|
|
140 |
June 2006, Bari |
June 2006, Bari |
141 |
|
|
142 |
The center of the scintillator planes S22Y (variable ZPAMS22Y in gpdgeo.inc) |
The center of the scintillator planes S22Y (variable ZPAMS22Y in gpdgeo.inc) |
143 |
and S12X (variable ZPAMS12X in gpdgeo.inc) has been positioned at the |
and S12X (variable ZPAMS12X in gpdgeo.inc) has been positioned at the |
146 |
E. Mocchiutti, S. Ricciarini, 20 December 2005). Follows that the positions |
E. Mocchiutti, S. Ricciarini, 20 December 2005). Follows that the positions |
147 |
of S21Y and S12X are higher than those in the cited document due to the fact |
of S21Y and S12X are higher than those in the cited document due to the fact |
148 |
that in GPAMELA the thickness of the mylar has been considered while in the |
that in GPAMELA the thickness of the mylar has been considered while in the |
149 |
document it has been neglected. |
document it has been neglected. |
150 |
|
|
151 |
|
|
152 |
May 2006, Bari & Tor Vergata |
May 2006, Bari & Tor Vergata |
172 |
ntple booking have been updated. The GPDCAL routine has been fixed |
ntple booking have been updated. The GPDCAL routine has been fixed |
173 |
so to handle the case in wich hits stored are more than the maximum |
so to handle the case in wich hits stored are more than the maximum |
174 |
number of hit. |
number of hit. |
175 |
In this case in the ntple up to the maximum number of hits will be stored. |
In this case in the ntple up to the maximum number of hits will be stored. |
176 |
|
|
177 |
April 2006, Bari |
April 2006, Bari |
178 |
|
|
215 |
is set to true in GPDAT if the data card: HPAK, is set to |
is set to true in GPDAT if the data card: HPAK, is set to |
216 |
'GCAL'. The gpkey.inc has been update accordingly. |
'GCAL'. The gpkey.inc has been update accordingly. |
217 |
|
|
218 |
|
|
219 |
NEUTRON DETECTOR ADDED. NEW DIR: GPND |
NEUTRON DETECTOR ADDED. NEW DIR: GPND |
220 |
|
|
221 |
The neutron detector has been added. At the moment it is just the |
The neutron detector has been added. At the moment it is just the |
232 |
http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html |
http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html |
233 |
This package should be distributed with the GEANT library but is |
This package should be distributed with the GEANT library but is |
234 |
not up to date. I did download the latest release and stored into |
not up to date. I did download the latest release and stored into |
235 |
gpcalor directory of the gpamela tree. |
gpcalor directory of the gpamela tree. |
236 |
Then I did clean up the code substituting the explicit inclusion of |
Then I did clean up the code substituting the explicit inclusion of |
237 |
the commons with a #include cpp directive. In parallel I did |
the commons with a #include cpp directive. In parallel I did |
238 |
extract the commons to include files having the same common name. I |
extract the commons to include files having the same common name. I |
243 |
NOTE: There are still problem due to different common sizes. In |
NOTE: There are still problem due to different common sizes. In |
244 |
particular the common MICFIL is maller in the geant library |
particular the common MICFIL is maller in the geant library |
245 |
libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are |
libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are |
246 |
present and linked using a wrong version of the common. This still needs to be debuged. |
present and linked using a wrong version of the common. This still needs to be debuged. |
247 |
NOTE2: The auxiliary files with the cross sections: chetc.dat.gz |
NOTE2: The auxiliary files with the cross sections: chetc.dat.gz |
248 |
and xsneut.dat.gz, have been added to the aux directory and moved |
and xsneut.dat.gz, have been added to the aux directory and moved |
249 |
to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will |
to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will |
251 |
searched there at first, then in the working directory. A fool |
searched there at first, then in the working directory. A fool |
252 |
proof policy has to be implemented to avoid problem with |
proof policy has to be implemented to avoid problem with |
253 |
synchronization fo these files. |
synchronization fo these files. |
|
|
|
254 |
|
|
255 |
The GCALOR package |
|
256 |
|
The GCALOR package |
257 |
|
|
258 |
June 2005, Bari |
June 2005, Bari |
259 |
|
|
264 |
|
|
265 |
May 2005, Bari |
May 2005, Bari |
266 |
|
|
267 |
Some updates on the latest modification done in the past year. |
Some updates on the latest modification done in the past year. |
268 |
|
|
269 |
NEW DATA CARD ADDED: HFSF |
NEW DATA CARD ADDED: HFSF |
270 |
|
|
271 |
To define a policy for the random number initial seeds |
To define a policy for the random number initial seeds |
272 |
definition. Using this card is possible to override GEANT seeds |
definition. Using this card is possible to override GEANT seeds |
273 |
defined via NRDM card. The policy is selected according to the |
defined via NRDM card. The policy is selected according to the |
274 |
values: |
values: |
275 |
|
|
276 |
- 1: The seeds are initialized to the initial values found in a user |
- 1: The seeds are initialized to the initial values found in a user |
277 |
defined file or the default file: INPUTSEED.DAT |
defined file or the default file: INPUTSEED.DAT |
278 |
|
|
279 |
- 2: The seeds are initialized to the final values found in a user defined |
- 2: The seeds are initialized to the final values found in a user defined |
280 |
file or the default file: INPUTSEED.DAT |
file or the default file: INPUTSEED.DAT |
281 |
|
|
282 |
The case 1 must be used in case the user needs to reproduce the |
The case 1 must be used in case the user needs to reproduce the |
285 |
binary file and pass the filename to the program using the *FLSF |
binary file and pass the filename to the program using the *FLSF |
286 |
data card. In case the user file is not specified the default |
data card. In case the user file is not specified the default |
287 |
INPUTSEED.DAT will be used. |
INPUTSEED.DAT will be used. |
288 |
|
|
289 |
The case 2 must be used in case the user needs to chain several |
The case 2 must be used in case the user needs to chain several |
290 |
GPAMELA run and likes to be sure he is starting the random |
GPAMELA run and likes to be sure he is starting the random |
291 |
generator using the right sequence. In this case the user must |
generator using the right sequence. In this case the user must |
303 |
filled with his own seeds. For this purpose an utility program: |
filled with his own seeds. For this purpose an utility program: |
304 |
writeseeds.f, has been added in the aux directory. In case the |
writeseeds.f, has been added in the aux directory. In case the |
305 |
*FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will |
*FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will |
306 |
be used as LUN and file name respectively. |
be used as LUN and file name respectively. |
307 |
|
|
308 |
NEW DATA CARD ADDED: *LSFI |
NEW DATA CARD ADDED: *LSFI |
309 |
|
|
310 |
Using this card the user can specify the logical unit and name of |
Using this card the user can specify the logical unit and name of |
312 |
run. This file is a FORTRAN binary one. This file can be used as |
run. This file is a FORTRAN binary one. This file can be used as |
313 |
input one specifying it in the *FSFI data card of the next GPAMELA |
input one specifying it in the *FSFI data card of the next GPAMELA |
314 |
run. In case the *LSFI card is not specified the default values: 26 |
run. In case the *LSFI card is not specified the default values: 26 |
315 |
and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN |
and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN |
316 |
and file name respectively. |
and file name respectively. |
317 |
|
|
318 |
NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f |
NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f |
319 |
|
|
320 |
These new programs have been added in the aux directory. Using these a |
These new programs have been added in the aux directory. Using these a |
347 |
has been revied and new sets of calibration constants calculated |
has been revied and new sets of calibration constants calculated |
348 |
comparing the beam test data with the GPAMELA results. The new |
comparing the beam test data with the GPAMELA results. The new |
349 |
constants are about 3% larger than the previous ones. |
constants are about 3% larger than the previous ones. |
350 |
|
|
351 |
TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI |
TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI |
352 |
|
|
353 |
Thanks to Lorenzo Bonechi for the drawings and explanations. Now the |
Thanks to Lorenzo Bonechi for the drawings and explanations. Now the |
354 |
hybrd cards have been put into the simulation and the geometry updated |
hybrd cards have been put into the simulation and the geometry updated |
355 |
considering the dead zones in the silicon detectors. The hybrid zone |
considering the dead zones in the silicon detectors. The hybrid zone |
356 |
has been simulated as well. At the moment the hybrid is simulated as |
has been simulated as well. At the moment the hybrid is simulated as |
357 |
a G10 plates. The full height of the tracker magnet has been |
a G10 plates. The full height of the tracker magnet has been |
358 |
reviewed as well. |
reviewed as well. |
359 |
|
|
360 |
The tracker ladder is now simulated inside a nitrogen box: TPAS, |
The tracker ladder is now simulated inside a nitrogen box: TPAS, |
361 |
placed inside an aluminum frame: TRPB. Each silicon ladder has been |
placed inside an aluminum frame: TRPB. Each silicon ladder has been |
363 |
smaller silicon detector: TPAI, has been placed inside the larger |
smaller silicon detector: TPAI, has been placed inside the larger |
364 |
silicon block TRSL. In this way the subdivided silicon ladder can |
silicon block TRSL. In this way the subdivided silicon ladder can |
365 |
be upgraded with an indipendend roto-translation for each sensor. |
be upgraded with an indipendend roto-translation for each sensor. |
366 |
|
|
367 |
The TRPB aluminum frame has been enlarged to fit the external |
The TRPB aluminum frame has been enlarged to fit the external |
368 |
magnet canister frame. |
magnet canister frame. |
369 |
|
|
370 |
The last plane has been flipped with a 180 degree rotation around |
The last plane has been flipped with a 180 degree rotation around |
371 |
the X axis. |
the X axis. |
372 |
|
|
373 |
TRACKER HIT STRUCTURE REVIEWED |
TRACKER HIT STRUCTURE REVIEWED |
374 |
|
|
375 |
Taking into account the new version of the tracker geometry, the hit |
Taking into account the new version of the tracker geometry, the hit |
376 |
structure for this detector has been revied. |
structure for this detector has been revied. |
377 |
|
|
378 |
CALORIMETER GEOMETRY REVIEWED |
CALORIMETER GEOMETRY REVIEWED |
379 |
|
|
380 |
Marco Albi reviewed the calorimeter dimentions and positioning. |
Marco Albi reviewed the calorimeter dimentions and positioning. |
381 |
|
|
382 |
|
|
383 |
29 March 2004, Bari |
29 March 2004, Bari |
391 |
the GARFIELD code for details). Follows that different runs produce |
the GARFIELD code for details). Follows that different runs produce |
392 |
different results without changing the initial seeds. To have identical |
different results without changing the initial seeds. To have identical |
393 |
results in different runs, the GARFIELD program has to start typing the |
results in different runs, the GARFIELD program has to start typing the |
394 |
noRNDM_initialisation switch. To avoid of specifying this switch |
noRNDM_initialisation switch. To avoid of specifying this switch |
395 |
by the user, |
by the user, |
396 |
the GARFIELD package has been upgraded with a patch. In this way the problem |
the GARFIELD package has been upgraded with a patch. In this way the problem |
397 |
is partially solved because, now, the initial seeds of the random generators |
is partially solved because, now, the initial seeds of the random generators |
399 |
activated by the user for changing the initial seeds in the GPAMELA program. |
activated by the user for changing the initial seeds in the GPAMELA program. |
400 |
Work is in progress for a more general correction of this problem. |
Work is in progress for a more general correction of this problem. |
401 |
Please, use the updated GARFIELD code released with the CVS version v4r1 |
Please, use the updated GARFIELD code released with the CVS version v4r1 |
402 |
to fix this problem. |
to fix this problem. |
403 |
|
|
404 |
|
|
405 |
RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS. |
RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS. |
412 |
the GPSED common when the set and the detector ZEBRA banks are filled |
the GPSED common when the set and the detector ZEBRA banks are filled |
413 |
during a run. In general the values of the set and detector addresses |
during a run. In general the values of the set and detector addresses |
414 |
depend on the number of active detectors in a given run. ISCAL=12 and |
depend on the number of active detectors in a given run. ISCAL=12 and |
415 |
IDCASI=1 are only right when all the detectors of GPAMELA are active. |
IDCASI=1 are only right when all the detectors of GPAMELA are active. |
416 |
|
|
417 |
9 December 2003, Bari |
9 December 2003, Bari |
418 |
|
|
456 |
-ICAPLANE(i) contains the number of hit plane; |
-ICAPLANE(i) contains the number of hit plane; |
457 |
-ICASTRIP(i) contains the number of hit strip; |
-ICASTRIP(i) contains the number of hit strip; |
458 |
-ICAMOD(i) can assume different values based on the number of times and |
-ICAMOD(i) can assume different values based on the number of times and |
459 |
positions in which a `long strip' has been hit. |
positions in which a `long strip' has been hit. |
460 |
-ENESTRIP(i) contains the deposited energy in the hit strip; |
-ENESTRIP(i) contains the deposited energy in the hit strip; |
461 |
where i is the number of hit (1<i<4224). |
where i is the number of hit (1<i<4224). |
462 |
Note that in the calorimeter each hit is filled at the end of the event and |
Note that in the calorimeter each hit is filled at the end of the event and |
482 |
The energy deposited in the calorimeter is calculated in GeV. |
The energy deposited in the calorimeter is calculated in GeV. |
483 |
|
|
484 |
To store the hits in the calorimeter the subroutine GSAHIT is used instead of |
To store the hits in the calorimeter the subroutine GSAHIT is used instead of |
485 |
GSCHIT. |
GSCHIT. |
486 |
|
|
487 |
To retrieve the hit structure the call to the routine GPRHIT is done instead |
To retrieve the hit structure the call to the routine GPRHIT is done instead |
488 |
of a call to the GFHITS subroutine. |
of a call to the GFHITS subroutine. |
489 |
|
|
490 |
25 February 2003, Bari |
25 February 2003, Bari |
491 |
|
|
492 |
BUG found: |
BUG found: |
493 |
DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP, |
DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP, |
494 |
DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2, |
DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2, |
512 |
Updated the relative positions of the calorimeter planes. The corrected |
Updated the relative positions of the calorimeter planes. The corrected |
513 |
shifting are: |
shifting are: |
514 |
|
|
515 |
first view: (Dxo,Dyo)=(0.10,0.05) cm |
first view: (Dxo,Dyo)=(0.10,0.05) cm |
516 |
second view: (Dxo,Dyo)=(-0.05,0.10) cm |
second view: (Dxo,Dyo)=(-0.05,0.10) cm |
517 |
third view: (Dxo,Dyo)=(-0.10,-0.05) cm |
third view: (Dxo,Dyo)=(-0.10,-0.05) cm |
518 |
fourth view: (Dxo,Dyo)=(0.05,-0.10) cm |
fourth view: (Dxo,Dyo)=(0.05,-0.10) cm |
519 |
|
|
520 |
4 November 2002, Bari |
4 November 2002, Bari |
521 |
|
|
522 |
CAS detectors distances modified |
CAS detectors distances modified |
523 |
|
|
524 |
The distances between the CAS detectors have been modified based on the |
The distances between the CAS detectors have been modified based on the |
525 |
latest CAD drawings. |
latest CAD drawings. |
526 |
|
|
527 |
2 November 2002, Bari |
2 November 2002, Bari |
528 |
|
|
539 |
CALB and CALS. CALS is an aluminum box of dimensions: 48.4*48.4*21.278 cm^3, |
CALB and CALS. CALS is an aluminum box of dimensions: 48.4*48.4*21.278 cm^3, |
540 |
having side-walls 1 cm thick and a bottom of 1 mm. The real box is more |
having side-walls 1 cm thick and a bottom of 1 mm. The real box is more |
541 |
complicated and the configuration of the bottom should be upgraded if we want |
complicated and the configuration of the bottom should be upgraded if we want |
542 |
a reliable description of the event in the S4 scintillator. |
a reliable description of the event in the S4 scintillator. |
543 |
|
|
544 |
22 October 2002, Stockholm |
22 October 2002, Stockholm |
545 |
|
|
1203 |
CAT dimension modified |
CAT dimension modified |
1204 |
CAT dimensions modified after PAOLO simulation and Tom design. Here |
CAT dimensions modified after PAOLO simulation and Tom design. Here |
1205 |
are the new dimensions: |
are the new dimensions: |
1206 |
DATA CATB/22.05,21.05,0.55/ |
DATA CATB/22.05,21.05,0.55/ |
1207 |
DATA CATL/15.75,6.675,0.55/ |
DATA CATL/15.75,6.675,0.55/ |
1208 |
DATA CATT/6.3,14.375,0.55/ |
DATA CATT/6.3,14.375,0.55/ |
1209 |
DATA ZCAT/75.75/ |
DATA ZCAT/75.75/ |