| 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.3 2002/12/05 17:27:59 pamela Exp $ |
# $Id: v_100.txt,v 3.27 2007/10/02 17:58:32 cafagna Exp $ |
| 6 |
# |
# |
| 7 |
# $Log: v_100.txt,v $ |
# $Log: v_100.txt,v $ |
| 8 |
|
# Revision 3.27 2007/10/02 17:58:32 cafagna |
| 9 |
|
# Power law spectra introduced |
| 10 |
|
# |
| 11 |
|
# Revision 3.26 2007/03/30 15:24:54 cafagna |
| 12 |
|
# Calo geometry modified. X-view plane shifts in X coordinates have been swapped, i.e. they do start with X=0.05 instead of X=-0.05 |
| 13 |
|
# |
| 14 |
|
# Revision 3.25 2006/12/17 14:53:18 cafagna |
| 15 |
|
# Full review of the CARD geometry. S1 table added |
| 16 |
|
# |
| 17 |
|
# Revision 3.20 2006/11/16 18:45:29 pam-ba |
| 18 |
|
# Simulated an aluminum container for S4 |
| 19 |
|
# |
| 20 |
|
# Revision 3.19 2006/11/16 12:01:13 pam-ba |
| 21 |
|
# S4 dimensions corrected. |
| 22 |
|
# |
| 23 |
|
# Revision 3.18 2006/11/10 11:39:35 pam-ba |
| 24 |
|
# S2 and S1 z-positions corrected, He3 and plystyrene mixture added, Top Plate geometry simulated and titanium mixture added. |
| 25 |
|
# |
| 26 |
|
# Revision 3.17 2006/10/13 16:36:59 pam-ba |
| 27 |
|
# Added a new material, the cadmium, for ND |
| 28 |
|
# |
| 29 |
|
# Revision 3.16 2006/10/12 11:11:21 pam-ba |
| 30 |
|
# ND geometry updated. |
| 31 |
|
# |
| 32 |
|
# Revision 3.15 2006/10/02 11:17:30 pam-ba |
| 33 |
|
# NDET 'SPHE' data card meaning changed. Now it eliminates the whole PAMELA container. |
| 34 |
|
# |
| 35 |
|
# Revision 3.14 2006/06/30 15:38:16 pam-ba |
| 36 |
|
# 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) |
| 37 |
|
# |
| 38 |
|
# Revision 3.13 2006/06/05 13:56:17 pamela |
| 39 |
|
# Gigantic resonance added for gamma enetering in the calorimeter absorber |
| 40 |
|
# |
| 41 |
|
# Revision 3.12 2006/05/18 10:52:32 pam-ba |
| 42 |
|
# TOF geometry completed and a new material, the polystyrene (density 35 g/l), added |
| 43 |
|
# |
| 44 |
|
# Revision 3.11 2006/05/11 23:53:15 cafagna |
| 45 |
|
# More bugs fixed in the CALO ntple structure filling |
| 46 |
|
# |
| 47 |
|
# Revision 3.10 2006/04/10 11:07:43 cafagna |
| 48 |
|
# GEN data card updated, ZDGEN added |
| 49 |
|
# |
| 50 |
|
# Revision 3.9 2005/12/14 03:34:40 cafagna |
| 51 |
|
# An update of the history and inform readme files. |
| 52 |
|
# |
| 53 |
|
# Revision 3.8 2005/12/14 03:16:08 cafagna |
| 54 |
|
# Neutron detector added. Geometry and GPCALOR package |
| 55 |
|
# |
| 56 |
|
# Revision 3.7 2005/10/18 08:24:35 cafagna |
| 57 |
|
# History updated |
| 58 |
|
# |
| 59 |
|
# Revision 3.6 2005/07/25 11:53:21 cafagna |
| 60 |
|
# Several updates. See history for details |
| 61 |
|
# |
| 62 |
|
# Revision 3.5 2004/04/06 10:33:46 pamela |
| 63 |
|
# NON-REPRODUCIBILITY problem of a GPAMELA RUN fixed; bug found and fixed filling in the hit structure of the calorimeter |
| 64 |
|
# |
| 65 |
|
# Revision 3.4 2003/12/17 11:32:50 pamela |
| 66 |
|
# CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation |
| 67 |
|
# |
| 68 |
# Revision 3.3 2002/12/05 17:27:59 pamela |
# Revision 3.3 2002/12/05 17:27:59 pamela |
| 69 |
# New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated |
# New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated |
| 70 |
# |
# |
| 88 |
#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 |
| 89 |
#-- Author : Francesco Cafagna 28/11/95 |
#-- Author : Francesco Cafagna 28/11/95 |
| 90 |
|
|
| 91 |
|
Jan 2008, Bari |
| 92 |
|
|
| 93 |
|
PRESSURIZED CONTAINER GEOMETRY: |
| 94 |
|
Pressurized container geometry has been updated: top container TSPH |
| 95 |
|
now has been placed at the maximum distance, quoted in the ICD, |
| 96 |
|
from S1: 5cm. This distance corresponds to the maximum one in |
| 97 |
|
flight due to the cover deformation due to the internal |
| 98 |
|
pressure. Three thermal cover layers (nylon6, mylar and aluminium, |
| 99 |
|
respectively in the new volumes COV1, COV2 and COV3) have been |
| 100 |
|
added to the top container. |
| 101 |
|
|
| 102 |
|
NEW MATERIALS AND TRACKING MEDIA: |
| 103 |
|
|
| 104 |
|
Two new materials have been implemented: AMg6M mixture for volume |
| 105 |
|
TSPH, and Nylon6 mixture for volume COV1. New tracking medium have |
| 106 |
|
been added for Nylon6 Mixture (TNYL), AMg6M Mixture (TAMG) and Iron |
| 107 |
|
(TIRO). |
| 108 |
|
|
| 109 |
|
TOF GEOMETRY: |
| 110 |
|
TOF geometry has been updated: Scintillators positions in the |
| 111 |
|
GPAMELA reference system have been corrected shifting all |
| 112 |
|
scintillator paddle by 23 microns. |
| 113 |
|
|
| 114 |
|
TRACK NUMBER RECORDED FOR TOF: |
| 115 |
|
The TOF structure now register the stack number into the HIT |
| 116 |
|
structure. In GUSTEP the stack number is tested, if is greater than |
| 117 |
|
0, i.e. the particle is not the primary, the particle number is set |
| 118 |
|
to 1000+IPART and passed to GPUTOF. Accordingly, now IPARTOF |
| 119 |
|
variable can assume 2000 values instead of 100. |
| 120 |
|
|
| 121 |
|
POWER LAW SPECTRA GENERATION: |
| 122 |
|
PLAW datacard parameters have been increased to 4. The meaning of |
| 123 |
|
the 4th parameter is: extract random kinetic energy if PLAW(4)<0., |
| 124 |
|
or extract random momentum if PLAW(4)>=0.. |
| 125 |
|
|
| 126 |
|
Sept 2007, Bari |
| 127 |
|
|
| 128 |
|
POWER LAW SPECTRA GENERATION: |
| 129 |
|
It is now possible to select a power law spectra. Using the |
| 130 |
|
datacard: PLAW, is possible to pass the spectral index and the |
| 131 |
|
kinetic energy range. In PAWL(1) the user can pass the signed power |
| 132 |
|
law index, i.e. PAWL -2.5, will generate a power law spectrum |
| 133 |
|
according to exp(-2.5). Please note that the range is in kinetic |
| 134 |
|
energy. In gukine the generated particle charateristics will be |
| 135 |
|
retrived and the mass calculated from this infos. The KINE data |
| 136 |
|
card assigned momentum will be overriden as well. |
| 137 |
|
|
| 138 |
|
PATHLENGHT CALCULATION MODIFIED IN TOF: |
| 139 |
|
The pathlength calulation saved in the TOF hit structure has now |
| 140 |
|
been modified. Insted of the pathlenght inside each paddle the |
| 141 |
|
average pathlenght from the generation vertex to the hit in the |
| 142 |
|
paddle, is saved. In this way it will be possible to calculate the |
| 143 |
|
beta using this variable. |
| 144 |
|
|
| 145 |
|
ERROR IN MATERIAL ASSIGNMENT FOR NEUTRON DETECTOR: |
| 146 |
|
The material assigned to the NDTU volume was wrong. It was |
| 147 |
|
scintillator, IRON has been asigned now. |
| 148 |
|
|
| 149 |
|
March 2007, Bari |
| 150 |
|
|
| 151 |
|
CALO GEOMETRY: |
| 152 |
|
An error was found in the original document describing PAMELA |
| 153 |
|
geometry. The shift along X axis of the X view planes are swapped, |
| 154 |
|
i.e. they don't start with X=+0.05 but with X=-0.05. The sequence |
| 155 |
|
is: 0.05m -0.05, 0.05, -0.05 and so on. |
| 156 |
|
|
| 157 |
|
MAKEFILE FIXED: |
| 158 |
|
Option -fdebug_kludge has been discharged |
| 159 |
|
|
| 160 |
|
GASATT WRNING FIXED: |
| 161 |
|
There was a warning message from a graphic routine. The GSATT calls has been commented out. |
| 162 |
|
|
| 163 |
|
|
| 164 |
|
December 2006, Bari |
| 165 |
|
|
| 166 |
|
CARD GEOMETRY: |
| 167 |
|
Positions, dimensions and quotes have been updated. Indeed, the aluminium table |
| 168 |
|
under S1, legs included, and the magnetic screen between the CARD and S2 have |
| 169 |
|
been added. |
| 170 |
|
|
| 171 |
|
#------------------------------------------------------------------------- |
| 172 |
|
|
| 173 |
|
November 2006, Bari |
| 174 |
|
|
| 175 |
|
|
| 176 |
|
S4 GEOMETRY UPDATED: |
| 177 |
|
The x and y dimensions have been corrected because they were wrong. |
| 178 |
|
A container of aluminum has been added having walls of thickness equal |
| 179 |
|
to 0.1 cm. |
| 180 |
|
|
| 181 |
|
TOF POSITIONS DEFINITELY UPDATED: |
| 182 |
|
S2 and S1 have been positioned again, after the simulation of the top plate. |
| 183 |
|
Before now, the positions were put ad hoc based on the positions given by |
| 184 |
|
Sergio Ricciarini in the document 'Main geometrical parameters of the PAMELA |
| 185 |
|
sub-detectors' released by O. Adriani, L. Bonechi, E. Mocchiutti and |
| 186 |
|
S. Ricciarini on the 20th of December 2005.. These positions were |
| 187 |
|
lightly approximated. |
| 188 |
|
|
| 189 |
|
ND GEOMETRY COMPLETED. |
| 190 |
|
The He3 and the polyethylene (CH2) have been added to fill |
| 191 |
|
the volumes NDTI and NDPB in the neutron detector. Leonov gave us the |
| 192 |
|
details of these materials. |
| 193 |
|
|
| 194 |
|
TOP PLATE GEOMETRY ADDED. |
| 195 |
|
The top plate geometry has been simulated as a N2 box that has the same |
| 196 |
|
dimensions of CATA. Inside it, at the right positions, there are S2 and CATA. |
| 197 |
|
The real top plate is a rectangular window having a thickness of 0.5 cm |
| 198 |
|
(volume TPTL emptied by the volume TPCV). This volume is filled with a mixture |
| 199 |
|
of titanium (90%), aluminum (6%) and vanadium (4%). Other parts of titanium |
| 200 |
|
are the four volumes TPTU, which are positioned at the corners of TPTL at the |
| 201 |
|
same height of CATA, and the four volumes TPTM, which are positioned at the |
| 202 |
|
corners of TPTL, at the same height of S2. |
| 203 |
|
The geometry of the top plate is based on a simplified version of the CAD |
| 204 |
|
drawings. |
| 205 |
|
|
| 206 |
|
October 2006, Bari |
| 207 |
|
|
| 208 |
|
The geometry of the neutron detector has been updated. Some dimensions of some |
| 209 |
|
volumes have been corrected and the aluminum cover and the aluminum boxes to |
| 210 |
|
put the cables have been added. A new material, the cadmium, has been also |
| 211 |
|
defined. |
| 212 |
|
|
| 213 |
|
September 2006, Bari |
| 214 |
|
|
| 215 |
|
SPHE and ND data card bugs fixed: the definition of the ND data card, |
| 216 |
|
missing in the subroutine gpgeo.F, has been added; the meaning of the SPHE |
| 217 |
|
data card has been changed. Before the correction the data card: |
| 218 |
|
NDET 'SPHE' was used to delete the spherical top shell to substitute it with |
| 219 |
|
a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA. |
| 220 |
|
|
| 221 |
|
|
| 222 |
|
June 2006, Bari |
| 223 |
|
|
| 224 |
|
The center of the scintillator planes S22Y (variable ZPAMS22Y in gpdgeo.inc) |
| 225 |
|
and S12X (variable ZPAMS12X in gpdgeo.inc) has been positioned at the |
| 226 |
|
nominal height as measured in PAMELA (See the document: "Main geometrical |
| 227 |
|
parameters of the PAMELA sub-detectors" by O. Adriani, L. Bonechi, |
| 228 |
|
E. Mocchiutti, S. Ricciarini, 20 December 2005). Follows that the positions |
| 229 |
|
of S21Y and S12X are higher than those in the cited document due to the fact |
| 230 |
|
that in GPAMELA the thickness of the mylar has been considered while in the |
| 231 |
|
document it has been neglected. |
| 232 |
|
|
| 233 |
|
|
| 234 |
|
May 2006, Bari & Tor Vergata |
| 235 |
|
|
| 236 |
|
GIGANTIC RESONANCE FOR NEUTRON DETECTOR ADDED |
| 237 |
|
|
| 238 |
|
Routines to simulate the gigantic resonance of gammas in Tungsten |
| 239 |
|
have been added. The GPGIG routine is called in GUSTEP if a gamma |
| 240 |
|
enter the calorimeter absorber. This is the steering routine to |
| 241 |
|
simulate the production of neutrons from gigantic resonance. It |
| 242 |
|
does checks on STEP lenght. If the range is smaller than the other |
| 243 |
|
selected for that step, it does generate the neutron and stops the |
| 244 |
|
gamma. Please note that the neutron has now a new particle |
| 245 |
|
number. This is to tag the gigantic resonance neutrons. |
| 246 |
|
|
| 247 |
|
|
| 248 |
|
May 2006, Bari & Florence |
| 249 |
|
|
| 250 |
|
CAL HIT STRUCTURE BUGS FIXED |
| 251 |
|
|
| 252 |
|
The maximum number of hit is now different for the two hit |
| 253 |
|
structures: CALST and CALI. Vectors inizialization and HBOOK |
| 254 |
|
ntple booking have been updated. The GPDCAL routine has been fixed |
| 255 |
|
so to handle the case in wich hits stored are more than the maximum |
| 256 |
|
number of hit. |
| 257 |
|
In this case in the ntple up to the maximum number of hits will be stored. |
| 258 |
|
|
| 259 |
|
April 2006, Bari |
| 260 |
|
|
| 261 |
|
TOF GEOMETRY AND POSITIONS UPDATED AND NEW MIXTURES ADDED |
| 262 |
|
|
| 263 |
|
The TOF geometry has been modified. The following boxes have been |
| 264 |
|
added: POL1, POL2 and POLY made of polystyrene, S11M, S12M, S21M, |
| 265 |
|
S22M, S31M and S32M made of mylar, S1A, S2A and S3 made of air and |
| 266 |
|
S1 and S2 made of aluminum. Each scintillator paddle has been put |
| 267 |
|
in his mylar box and the other materials: air, polystyrene, and |
| 268 |
|
aluminum have been added at their nominal positions. According to |
| 269 |
|
Naples people the araldite glue has been simulated has an air |
| 270 |
|
gap. For this work two new materials: the Mylar (MYLAR) and the |
| 271 |
|
polystyrene (POLYSTYRENE) with a density of 35 g/l have been |
| 272 |
|
defined as a mixture. The positions of the three bottom |
| 273 |
|
scintillator planes that contain respectively the S12X, S22Y and |
| 274 |
|
S32X paddles have been regulated according on their official |
| 275 |
|
positions in PAMELA. |
| 276 |
|
|
| 277 |
|
Mar 2006, Bari |
| 278 |
|
|
| 279 |
|
GEN DATA CARD UPDATED |
| 280 |
|
|
| 281 |
|
To enable generation on a surface perpendicular to the XY plane, |
| 282 |
|
GEN gata card has been updated addingh a new parameter: ZDGEN. This is |
| 283 |
|
the dimension, along Z axis , of the generation surface. The Z |
| 284 |
|
position will be randomply chosen according to: Z= ZDGEN*RNDM_NUMBER + |
| 285 |
|
ZGEN, i.e. Z= GEN(6)*RNDM_NOMBER + GEN(3). |
| 286 |
|
|
| 287 |
|
Nov 2005, Bari |
| 288 |
|
|
| 289 |
|
GUHADR AND GUPHAD UPDATED |
| 290 |
|
|
| 291 |
|
To use GCALOR package the hadronic routines have been updated. The |
| 292 |
|
inizialization routine call CALSIG, while the other calls GCALOR. |
| 293 |
|
|
| 294 |
|
NEW GPKEY ADDED: GPCALOR |
| 295 |
|
|
| 296 |
|
This logical has been added to enable the GCALOR package. This flag |
| 297 |
|
is set to true in GPDAT if the data card: HPAK, is set to |
| 298 |
|
'GCAL'. The gpkey.inc has been update accordingly. |
| 299 |
|
|
| 300 |
|
|
| 301 |
|
NEUTRON DETECTOR ADDED. NEW DIR: GPND |
| 302 |
|
|
| 303 |
|
The neutron detector has been added. At the moment it is just the |
| 304 |
|
geometry. The directory structure of the repository has been |
| 305 |
|
updated as well. Dimensions has been taken from picture and |
| 306 |
|
literature. A full upgrade to the drawing is needed. |
| 307 |
|
|
| 308 |
|
GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR, GPCALORDES |
| 309 |
|
|
| 310 |
|
GCALOR package contins the CALOR simulation code and an interface |
| 311 |
|
to use it in GEANT. The important feature for us is the usage of |
| 312 |
|
the MICAP code. This is facused on the low energy neutron |
| 313 |
|
simulation. for details see: |
| 314 |
|
http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html |
| 315 |
|
This package should be distributed with the GEANT library but is |
| 316 |
|
not up to date. I did download the latest release and stored into |
| 317 |
|
gpcalor directory of the gpamela tree. |
| 318 |
|
Then I did clean up the code substituting the explicit inclusion of |
| 319 |
|
the commons with a #include cpp directive. In parallel I did |
| 320 |
|
extract the commons to include files having the same common name. I |
| 321 |
|
did store the include files into a newly created directory: |
| 322 |
|
gpcalordes. |
| 323 |
|
The Makefile has been updated accordingly. |
| 324 |
|
Please note that to avoid conflict with CRENLIB distribution the gcalor source file has been named gpcalor.F |
| 325 |
|
NOTE: There are still problem due to different common sizes. In |
| 326 |
|
particular the common MICFIL is maller in the geant library |
| 327 |
|
libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are |
| 328 |
|
present and linked using a wrong version of the common. This still needs to be debuged. |
| 329 |
|
NOTE2: The auxiliary files with the cross sections: chetc.dat.gz |
| 330 |
|
and xsneut.dat.gz, have been added to the aux directory and moved |
| 331 |
|
to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will |
| 332 |
|
look for CERN_ROOT environment variable. If found files are |
| 333 |
|
searched there at first, then in the working directory. A fool |
| 334 |
|
proof policy has to be implemented to avoid problem with |
| 335 |
|
synchronization fo these files. |
| 336 |
|
|
| 337 |
|
|
| 338 |
|
The GCALOR package |
| 339 |
|
|
| 340 |
|
June 2005, Bari |
| 341 |
|
|
| 342 |
|
TOF SCINTILLATOR PADDLES UPDATED |
| 343 |
|
|
| 344 |
|
The dimensions and the number of the scintillator paddles for each |
| 345 |
|
TOF planes have been updated. |
| 346 |
|
|
| 347 |
|
May 2005, Bari |
| 348 |
|
|
| 349 |
|
Some updates on the latest modification done in the past year. |
| 350 |
|
|
| 351 |
|
NEW DATA CARD ADDED: HFSF |
| 352 |
|
|
| 353 |
|
To define a policy for the random number initial seeds |
| 354 |
|
definition. Using this card is possible to override GEANT seeds |
| 355 |
|
defined via NRDM card. The policy is selected according to the |
| 356 |
|
values: |
| 357 |
|
|
| 358 |
|
- 1: The seeds are initialized to the initial values found in a user |
| 359 |
|
defined file or the default file: INPUTSEED.DAT |
| 360 |
|
|
| 361 |
|
- 2: The seeds are initialized to the final values found in a user defined |
| 362 |
|
file or the default file: INPUTSEED.DAT |
| 363 |
|
|
| 364 |
|
The case 1 must be used in case the user needs to reproduce the |
| 365 |
|
random chain of a previous run. In this case the user can save the |
| 366 |
|
initial seeds, used in the run he would like to reproduce, in a |
| 367 |
|
binary file and pass the filename to the program using the *FLSF |
| 368 |
|
data card. In case the user file is not specified the default |
| 369 |
|
INPUTSEED.DAT will be used. |
| 370 |
|
|
| 371 |
|
The case 2 must be used in case the user needs to chain several |
| 372 |
|
GPAMELA run and likes to be sure he is starting the random |
| 373 |
|
generator using the right sequence. In this case the user must |
| 374 |
|
specify an input binary file using the *FLSF data card, otherwise |
| 375 |
|
the INPUTSEED.DAT file will be used. |
| 376 |
|
|
| 377 |
|
NEW DATA CARD ADDED: *FSFI |
| 378 |
|
|
| 379 |
|
Using this card the user can specify the logical unit and name of |
| 380 |
|
the file storing the initial seeds to be used to initialize the |
| 381 |
|
random number generator. This file must be a FORTRAN binary one |
| 382 |
|
storing four integer numbers. The first two are the number to be |
| 383 |
|
used in the case: HFSF=1, the other two will be used in the case: |
| 384 |
|
HFSF=2. This file can be one created by GPAMELA or by the user |
| 385 |
|
filled with his own seeds. For this purpose an utility program: |
| 386 |
|
writeseeds.f, has been added in the aux directory. In case the |
| 387 |
|
*FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will |
| 388 |
|
be used as LUN and file name respectively. |
| 389 |
|
|
| 390 |
|
NEW DATA CARD ADDED: *LSFI |
| 391 |
|
|
| 392 |
|
Using this card the user can specify the logical unit and name of |
| 393 |
|
the file storing the first and last seeds used in the GPAMELA |
| 394 |
|
run. This file is a FORTRAN binary one. This file can be used as |
| 395 |
|
input one specifying it in the *FSFI data card of the next GPAMELA |
| 396 |
|
run. In case the *LSFI card is not specified the default values: 26 |
| 397 |
|
and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN |
| 398 |
|
and file name respectively. |
| 399 |
|
|
| 400 |
|
NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f |
| 401 |
|
|
| 402 |
|
These new programs have been added in the aux directory. Using these a |
| 403 |
|
user defined seed file can be created and re-read. |
| 404 |
|
|
| 405 |
|
NEW VOLUMES ADDED: MSHE, BSPH; PRESSURIZED CONTAINER ADDED |
| 406 |
|
|
| 407 |
|
Alexey Bakaldin, in MEPHI, did add the PAMELA pressurized container to |
| 408 |
|
the simulation. He did defined new volumes filled with aluminum and |
| 409 |
|
placed inside the mother volume. Positions have been fine tuned by |
| 410 |
|
Marialuigia Ambriola and compared to the CAD drawings. |
| 411 |
|
Two new volumes have been added to simulate the container: |
| 412 |
|
- MSHE, a tube simulating the middle part of the container |
| 413 |
|
- BSPH, the spherical bottom part of the container |
| 414 |
|
|
| 415 |
|
To better simulate the upper part the SHEL volume has been modified |
| 416 |
|
into a cone. Dimentions of the top cover: TSPH, have been modified |
| 417 |
|
accordingly. |
| 418 |
|
|
| 419 |
|
DETECTOR POSITIONS REVIEWED |
| 420 |
|
|
| 421 |
|
All detector Z positions have been reviewd to fit into the |
| 422 |
|
simulated pressurized container. |
| 423 |
|
|
| 424 |
|
TRD GEOMETRY AND CALIBRATION REVIEWD |
| 425 |
|
|
| 426 |
|
The TRD geometry has been deeply reviewed. Using the CAD drawings |
| 427 |
|
the carbon fiber frames have been simulated and radiator dimentions |
| 428 |
|
corrected. For this reason the calibration done on the beam tests |
| 429 |
|
has been revied and new sets of calibration constants calculated |
| 430 |
|
comparing the beam test data with the GPAMELA results. The new |
| 431 |
|
constants are about 3% larger than the previous ones. |
| 432 |
|
|
| 433 |
|
TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI |
| 434 |
|
|
| 435 |
|
Thanks to Lorenzo Bonechi for the drawings and explanations. Now the |
| 436 |
|
hybrd cards have been put into the simulation and the geometry updated |
| 437 |
|
considering the dead zones in the silicon detectors. The hybrid zone |
| 438 |
|
has been simulated as well. At the moment the hybrid is simulated as |
| 439 |
|
a G10 plates. The full height of the tracker magnet has been |
| 440 |
|
reviewed as well. |
| 441 |
|
|
| 442 |
|
The tracker ladder is now simulated inside a nitrogen box: TPAS, |
| 443 |
|
placed inside an aluminum frame: TRPB. Each silicon ladder has been |
| 444 |
|
simulated using two silicon blocks: TRSL, into each of this block a |
| 445 |
|
smaller silicon detector: TPAI, has been placed inside the larger |
| 446 |
|
silicon block TRSL. In this way the subdivided silicon ladder can |
| 447 |
|
be upgraded with an indipendend roto-translation for each sensor. |
| 448 |
|
|
| 449 |
|
The TRPB aluminum frame has been enlarged to fit the external |
| 450 |
|
magnet canister frame. |
| 451 |
|
|
| 452 |
|
The last plane has been flipped with a 180 degree rotation around |
| 453 |
|
the X axis. |
| 454 |
|
|
| 455 |
|
TRACKER HIT STRUCTURE REVIEWED |
| 456 |
|
|
| 457 |
|
Taking into account the new version of the tracker geometry, the hit |
| 458 |
|
structure for this detector has been revied. |
| 459 |
|
|
| 460 |
|
CALORIMETER GEOMETRY REVIEWED |
| 461 |
|
|
| 462 |
|
Marco Albi reviewed the calorimeter dimentions and positioning. |
| 463 |
|
|
| 464 |
|
|
| 465 |
|
29 March 2004, Bari |
| 466 |
|
|
| 467 |
|
NON-REPRODUCIBILITY PROBLEM OF A GPAMELA RUN FIXED. |
| 468 |
|
The non-reproducibility of a GPAMELA run was due to the random number |
| 469 |
|
initialization in the GARFIELD code. In GARFIELD by default, the initial |
| 470 |
|
seeds of the random number generators are always the same while the random |
| 471 |
|
number generators are called a given number of times (determined by the |
| 472 |
|
hour of the day) during the initialization phase (see init.f subroutine in |
| 473 |
|
the GARFIELD code for details). Follows that different runs produce |
| 474 |
|
different results without changing the initial seeds. To have identical |
| 475 |
|
results in different runs, the GARFIELD program has to start typing the |
| 476 |
|
noRNDM_initialisation switch. To avoid of specifying this switch |
| 477 |
|
by the user, |
| 478 |
|
the GARFIELD package has been upgraded with a patch. In this way the problem |
| 479 |
|
is partially solved because, now, the initial seeds of the random generators |
| 480 |
|
in GARFIELD will be always the same even if the RNDM GEANT data card is |
| 481 |
|
activated by the user for changing the initial seeds in the GPAMELA program. |
| 482 |
|
Work is in progress for a more general correction of this problem. |
| 483 |
|
Please, use the updated GARFIELD code released with the CVS version v4r1 |
| 484 |
|
to fix this problem. |
| 485 |
|
|
| 486 |
|
|
| 487 |
|
RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS. |
| 488 |
|
The obsolete RNDM random number generator has been replaced by the GEANT |
| 489 |
|
GRNDN routine in the gpxtr.F subroutine and in the npoiss.F function. |
| 490 |
|
|
| 491 |
|
BUG FOUND AND FIXED: the set and detector calorimeter addresses (ISCAL |
| 492 |
|
and IDCASI variables) used in GUTREV were respectively set to a fixed |
| 493 |
|
values of 12 and 1. The correct values of these variables are stored in |
| 494 |
|
the GPSED common when the set and the detector ZEBRA banks are filled |
| 495 |
|
during a run. In general the values of the set and detector addresses |
| 496 |
|
depend on the number of active detectors in a given run. ISCAL=12 and |
| 497 |
|
IDCASI=1 are only right when all the detectors of GPAMELA are active. |
| 498 |
|
|
| 499 |
9 December 2003, Bari |
9 December 2003, Bari |
| 500 |
|
|
| 501 |
CALORIMETER SIMULATION completed! The update of the geometry and of the |
CALORIMETER SIMULATION completed! The update of the geometry and of the |
| 502 |
special tracking parameters and the tuning of the calorimeter have been |
special tracking parameters and the tuning of the calorimeter have been |
| 503 |
successfully done. A great quantity of simulated data have been produced |
successfully done. A great quantity of simulated data have been produced |
| 504 |
in the calorimeter for different particles (muons, electrons and pions) |
in the calorimeter for different particles (muons, electrons and pions) |
| 505 |
and momenta (5 and 40 GeV/c) and the output data have been analyzed. The |
and momenta (5 and 40 GeV/c) and the output data have been analyzed. The |
| 506 |
distributions of the total energy deposited in the calorimeter and the |
distributions of the total energy deposited in the calorimeter and the |
| 507 |
total number of strips hit have been compared with the respective |
total number of strips hit have been compared with the respective |
| 508 |
distributions produced by the Trieste's tuned standalone Monte Carlo |
distributions produced by the Trieste's tuned standalone Monte Carlo |
| 509 |
simulation program of the PAMELA calorimeter. The accord between the |
simulation program of the PAMELA calorimeter. The accord between the |
| 510 |
two simulations is excellent. Many thanks to Mirko for his collaboration. |
two simulations is excellent. Many thanks to Mirko for his collaboration. |
| 511 |
|
|
| 512 |
Working in progress on TRD. The GARFIELD interface to the HEED program is not |
Working in progress on TRD. The GARFIELD interface to the HEED program is not |
| 513 |
optimized to track particle with a charge greater than one and photons. The |
optimized to track particle with a charge greater than one and photons. The |
| 514 |
program print a warning message to advise the user when it is the case. |
program print a warning message to advise the user when it is the case. |
| 515 |
|
|
| 516 |
18 April 2003, Bari |
18 April 2003, Bari |
| 517 |
|
|
| 518 |
The buffer size of each column of the GPAMELA Ntuple has been increased to |
The buffer size of each column of the GPAMELA Ntuple has been increased to |
| 519 |
4096 and set equal to the record length, defined by a call to the HROPEN |
4096 and set equal to the record length, defined by a call to the HROPEN |
| 520 |
routine. |
routine. |
| 521 |
Also the length of the common /PAWC/ (parameter NWPAW) has been increased |
Also the length of the common /PAWC/ (parameter NWPAW) has been increased |
| 522 |
to 1.34E8, according to the rule that it has to be larger than the number |
to 1.34E8, according to the rule that it has to be larger than the number |
| 523 |
of columns times the buffer size. |
of columns times the buffer size. |
| 524 |
|
|
| 525 |
10 April 2003, Bari |
10 April 2003, Bari |
| 526 |
|
|
| 527 |
The variables in the HIT STRUCTURE of the CALORIMETER and their way to be |
The variables in the HIT STRUCTURE of the CALORIMETER and their way to be |
| 528 |
filled have been changed according to the electronics system of the real |
filled have been changed according to the electronics system of the real |
| 529 |
detector. In fact, because each silicon detector (module) consists of |
detector. In fact, because each silicon detector (module) consists of |
| 530 |
32 strips and each strip is connected to those belonging to the two detectors |
32 strips and each strip is connected to those belonging to the two detectors |
| 531 |
of the same row (or column) for forming 24 cm long strips, the sum of the |
of the same row (or column) for forming 24 cm long strips, the sum of the |
| 532 |
deposited energies in the strips forming a `long strip' is now calculated for |
deposited energies in the strips forming a `long strip' is now calculated for |
| 533 |
each event (gpucal.F subroutine) and it is stored in a hit only at the |
each event (gpucal.F subroutine) and it is stored in a hit only at the |
| 534 |
end of the event (gutrev.F subroutine). |
end of the event (gutrev.F subroutine). |
| 535 |
The output variables of the GPAMELA en-tuple are then filled in the vectors |
The output variables of the GPAMELA en-tuple are then filled in the vectors |
| 536 |
ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL), |
ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL), |
| 537 |
by a call to the GPDCAL subroutine: |
by a call to the GPDCAL subroutine: |
| 538 |
-ICAPLANE(i) contains the number of hit plane; |
-ICAPLANE(i) contains the number of hit plane; |
| 539 |
-ICASTRIP(i) contains the number of hit strip; |
-ICASTRIP(i) contains the number of hit strip; |
| 540 |
-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 |
| 541 |
positions in which a `long strip' has been hit. |
positions in which a `long strip' has been hit. |
| 542 |
-ENESTRIP(i) contains the deposited energy in the hit strip; |
-ENESTRIP(i) contains the deposited energy in the hit strip; |
| 543 |
where i is the number of hit (1<i<4224). |
where i is the number of hit (1<i<4224). |
| 544 |
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 |
| 545 |
that there is a hit for each `long strip' hit from |
that there is a hit for each `long strip' hit from |
| 546 |
the particle. This use of the hit structure is different for the other |
the particle. This use of the hit structure is different for the other |
| 547 |
detectors and it has been considered to avoid a too big number of hit in the |
detectors and it has been considered to avoid a too big number of hit in the |
| 548 |
calorimeter due to the showers. Follows that NTHCAL, which is the |
calorimeter due to the showers. Follows that NTHCAL, which is the |
| 549 |
max number of hit in the calorimeter, is equal to 4224, the total |
max number of hit in the calorimeter, is equal to 4224, the total |
| 550 |
number of `long strips'. So, for each event, the real number of hit will |
number of `long strips'. So, for each event, the real number of hit will |
| 551 |
be less or equal to 4224. |
be less or equal to 4224. |
| 552 |
ICAMOD(i) is an additional information that does not exist in the real |
ICAMOD(i) is an additional information that does not exist in the real |
| 553 |
detector: if the strip i (i=1,32) of the module 1 or 2 or 3 |
detector: if the strip i (i=1,32) of the module 1 or 2 or 3 |
| 554 |
is hit, the value of ICAMOD(i) is respectively incremented of 1, 100, 10000. |
is hit, the value of ICAMOD(i) is respectively incremented of 1, 100, 10000. |
| 555 |
Analogously it is done, if it is the strip j (j=33,64) of the modules 4, 5 |
Analogously it is done, if it is the strip j (j=33,64) of the modules 4, 5 |
| 556 |
and 6 or if it is the strip k (k=65,96) of the modules 7, 8 and 9. |
and 6 or if it is the strip k (k=65,96) of the modules 7, 8 and 9. |
| 557 |
For example if we consider the hit 1 of an event, we could read: |
For example if we consider the hit 1 of an event, we could read: |
| 558 |
ICASTRIP(1)=30, ICAPLANE(1)=21, ENESTRIP(1)=0.5E-03 and ICAMOD(1)=10001. |
ICASTRIP(1)=30, ICAPLANE(1)=21, ENESTRIP(1)=0.5E-03 and ICAMOD(1)=10001. |
| 559 |
It means that the hit 1 contains the information that in the strip 30 of the |
It means that the hit 1 contains the information that in the strip 30 of the |
| 560 |
plane 21 has been deposited a total energy of 0.5E-03 GeV. In addition the |
plane 21 has been deposited a total energy of 0.5E-03 GeV. In addition the |
| 561 |
`long strip 30' has been hit two times, one in the first module and the |
`long strip 30' has been hit two times, one in the first module and the |
| 562 |
other in the third one. |
other in the third one. |
| 563 |
|
|
| 564 |
The energy deposited in the calorimeter is calculated in GeV. |
The energy deposited in the calorimeter is calculated in GeV. |
| 565 |
|
|
| 566 |
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 |
| 567 |
GSCHIT. |
GSCHIT. |
| 568 |
|
|
| 569 |
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 |
| 570 |
of a call to the GFHITS subroutine. |
of a call to the GFHITS subroutine. |
| 571 |
|
|
| 572 |
25 February 2003, Bari |
25 February 2003, Bari |
| 573 |
|
|
| 574 |
BUG found: |
BUG found: |
| 575 |
DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP, |
DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP, |
| 576 |
DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2, |
DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2, |
| 577 |
DCUTEW, DCUTEXE variables missed in the commons: gpaer.inc, gpal.inc, gpce.inc, |
DCUTEW, DCUTEXE variables missed in the commons: gpaer.inc, gpal.inc, gpce.inc, |
| 578 |
gpcp.inc, gpfe.inc, gpg10c.inc, gpg10.inc, gpkap.inc, gpn2g.inc, gproa.inc, |
gpcp.inc, gpfe.inc, gpg10c.inc, gpg10.inc, gpkap.inc, gpn2g.inc, gproa.inc, |
| 579 |
gpscin.inc (obsolete), gpscint.inc, gpsica.inc, gptrad.inc, gpw2.inc, gpw.inc, |
gpscin.inc (obsolete), gpscint.inc, gpsica.inc, gptrad.inc, gpw2.inc, gpw.inc, |
| 580 |
gpxe.inc, gpdaer.inc, gpdal.inc, gpdce.inc, gpdcp.inc, gpdfe.inc, gpdg10c.inc, |
gpxe.inc, gpdaer.inc, gpdal.inc, gpdce.inc, gpdcp.inc, gpdfe.inc, gpdg10c.inc, |
| 581 |
gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc, |
gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc, |
| 582 |
gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc. |
gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc. |
| 583 |
They have been added in these commons and they have been initialized in the |
They have been added in these commons and they have been initialized in the |
| 584 |
GPSTM subroutine. |
GPSTM subroutine. |
| 585 |
|
|
| 586 |
Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C |
Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C |
| 587 |
in the subroutines gpsica.F, gpw2.F, gpw.F, gpce.F and gpg10c.F. They were |
in the subroutines gpsica.F, gpw2.F, gpw.F, gpce.F and gpg10c.F. They were |
| 588 |
suggested by Mirko Boezio. |
suggested by Mirko Boezio. |
| 589 |
|
|
| 590 |
Updated the value of the absorption length for silicon in the calorimeter |
Updated the value of the absorption length for silicon in the calorimeter |
| 591 |
and tracker although this parameter is ignored by GEANT. For this reason |
and tracker although this parameter is ignored by GEANT. For this reason |
| 592 |
it was equal to the radiation length. |
it was equal to the radiation length. |
| 593 |
|
|
| 594 |
Updated the relative positions of the calorimeter planes. The corrected |
Updated the relative positions of the calorimeter planes. The corrected |
| 595 |
shifting are: |
shifting are: |
| 596 |
|
|
| 597 |
first view: (Dxo,Dyo)=(0.10,0.05) cm |
first view: (Dxo,Dyo)=(0.10,0.05) cm |
| 598 |
second view: (Dxo,Dyo)=(-0.05,0.10) cm |
second view: (Dxo,Dyo)=(-0.05,0.10) cm |
| 599 |
third view: (Dxo,Dyo)=(-0.10,-0.05) cm |
third view: (Dxo,Dyo)=(-0.10,-0.05) cm |
| 600 |
fourth view: (Dxo,Dyo)=(0.05,-0.10) cm |
fourth view: (Dxo,Dyo)=(0.05,-0.10) cm |
| 601 |
|
|
| 602 |
4 November 2002, Bari |
4 November 2002, Bari |
| 603 |
|
|
| 604 |
CAS detectors distances modified |
CAS detectors distances modified |
| 605 |
|
|
| 606 |
The distances between the CAS detectors have been modified based on the |
The distances between the CAS detectors have been modified based on the |
| 607 |
latest CAD drawings. |
latest CAD drawings. |
| 608 |
|
|
| 609 |
2 November 2002, Bari |
2 November 2002, Bari |
| 610 |
|
|
| 611 |
CALORIMETER geometry upgrade |
CALORIMETER geometry upgrade |
| 612 |
|
|
| 613 |
The volumes CAPD and CAAD have been taken off from the calorimeter. |
The volumes CAPD and CAAD have been taken off from the calorimeter. |
| 614 |
In addition the logical tree has been slightly changed to make the shifts of |
In addition the logical tree has been slightly changed to make the shifts of |
| 615 |
the silicon planes into the calorimeter box easier, i.e. the CAPL volume, |
the silicon planes into the calorimeter box easier, i.e. the CAPL volume, |
| 616 |
which was made of the CASI, CAKP, CAGL, C10C and CAKA volumes, has |
which was made of the CASI, CAKP, CAGL, C10C and CAKA volumes, has |
| 617 |
been split up in the volumes CANS and CAPL. Now CANS is made of the CAKP, |
been split up in the volumes CANS and CAPL. Now CANS is made of the CAKP, |
| 618 |
CAGL, C10C and CAKA volumes while CAPL contains the CASI volume, that has to |
CAGL, C10C and CAKA volumes while CAPL contains the CASI volume, that has to |
| 619 |
be shifted as a function of the vertical position in the calorimeter. Also the |
be shifted as a function of the vertical position in the calorimeter. Also the |
| 620 |
dimensions of some volumes have been upgraded, including the external ones: |
dimensions of some volumes have been upgraded, including the external ones: |
| 621 |
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, |
| 622 |
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 |
| 623 |
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 |
| 624 |
a reliable description of the event in the S4 scintillator. |
a reliable description of the event in the S4 scintillator. |
| 625 |
|
|
| 626 |
22 October 2002, Stockholm |
22 October 2002, Stockholm |
| 627 |
|
|
| 1285 |
CAT dimension modified |
CAT dimension modified |
| 1286 |
CAT dimensions modified after PAOLO simulation and Tom design. Here |
CAT dimensions modified after PAOLO simulation and Tom design. Here |
| 1287 |
are the new dimensions: |
are the new dimensions: |
| 1288 |
DATA CATB/22.05,21.05,0.55/ |
DATA CATB/22.05,21.05,0.55/ |
| 1289 |
DATA CATL/15.75,6.675,0.55/ |
DATA CATL/15.75,6.675,0.55/ |
| 1290 |
DATA CATT/6.3,14.375,0.55/ |
DATA CATT/6.3,14.375,0.55/ |
| 1291 |
DATA ZCAT/75.75/ |
DATA ZCAT/75.75/ |
Parent Directory
| 
Revision Log
| 
Patch