--- gpamela/history/v_100.txt 2006/04/10 11:07:43 3.10 +++ gpamela/history/v_100.txt 2008/01/29 18:25:18 3.28 @@ -1,8 +1,52 @@ - +# Revision 3.21 2006/12/05 18:45:29 bruno +# Update of the CARD geometry +# Simulated also legs, table and magnetic screen # -# $Id: v_100.txt,v 3.9 2005/12/14 03:34:40 cafagna Exp $ +# $Id: v_100.txt,v 3.27 2007/10/02 17:58:32 cafagna Exp $ # # $Log: v_100.txt,v $ +# Revision 3.27 2007/10/02 17:58:32 cafagna +# Power law spectra introduced +# +# Revision 3.26 2007/03/30 15:24:54 cafagna +# 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 +# +# Revision 3.25 2006/12/17 14:53:18 cafagna +# Full review of the CARD geometry. S1 table added +# +# Revision 3.20 2006/11/16 18:45:29 pam-ba +# Simulated an aluminum container for S4 +# +# Revision 3.19 2006/11/16 12:01:13 pam-ba +# S4 dimensions corrected. +# +# Revision 3.18 2006/11/10 11:39:35 pam-ba +# S2 and S1 z-positions corrected, He3 and plystyrene mixture added, Top Plate geometry simulated and titanium mixture added. +# +# Revision 3.17 2006/10/13 16:36:59 pam-ba +# Added a new material, the cadmium, for ND +# +# Revision 3.16 2006/10/12 11:11:21 pam-ba +# ND geometry updated. +# +# Revision 3.15 2006/10/02 11:17:30 pam-ba +# NDET 'SPHE' data card meaning changed. Now it eliminates the whole PAMELA container. +# +# Revision 3.14 2006/06/30 15:38:16 pam-ba +# 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) +# +# Revision 3.13 2006/06/05 13:56:17 pamela +# Gigantic resonance added for gamma enetering in the calorimeter absorber +# +# Revision 3.12 2006/05/18 10:52:32 pam-ba +# TOF geometry completed and a new material, the polystyrene (density 35 g/l), added +# +# Revision 3.11 2006/05/11 23:53:15 cafagna +# More bugs fixed in the CALO ntple structure filling +# +# Revision 3.10 2006/04/10 11:07:43 cafagna +# GEN data card updated, ZDGEN added +# # Revision 3.9 2005/12/14 03:34:40 cafagna # An update of the history and inform readme files. # @@ -42,7 +86,193 @@ #CMZ : 1.00/03 30/04/96 12.23.59 by Francesco Cafagna #CMZ : 1.00/02 05/04/96 15.31.25 by Francesco Cafagna #CMZ : 1.00/01 28/11/95 18.51.23 by Francesco Cafagna -#-- Author : Francesco Cafagna 28/11/95 +#-- Author : Francesco Cafagna 28/11/95 + +Jan 2008, Bari + +PRESSURIZED CONTAINER GEOMETRY: + Pressurized container geometry has been updated: top container TSPH + now has been placed at the maximum distance, quoted in the ICD, + from S1: 5cm. This distance corresponds to the maximum one in + flight due to the cover deformation due to the internal + pressure. Three thermal cover layers (nylon6, mylar and aluminium, + respectively in the new volumes COV1, COV2 and COV3) have been + added to the top container. + +NEW MATERIALS AND TRACKING MEDIA: + + Two new materials have been implemented: AMg6M mixture for volume + TSPH, and Nylon6 mixture for volume COV1. New tracking medium have + been added for Nylon6 Mixture (TNYL), AMg6M Mixture (TAMG) and Iron + (TIRO). + +TOF GEOMETRY: + TOF geometry has been updated: Scintillators positions in the + GPAMELA reference system have been corrected shifting all + scintillator paddle by 23 microns. + +TRACK NUMBER RECORDED FOR TOF: + The TOF structure now register the stack number into the HIT + structure. In GUSTEP the stack number is tested, if is greater than + 0, i.e. the particle is not the primary, the particle number is set + to 1000+IPART and passed to GPUTOF. Accordingly, now IPARTOF + variable can assume 2000 values instead of 100. + +POWER LAW SPECTRA GENERATION: + PLAW datacard parameters have been increased to 4. The meaning of + the 4th parameter is: extract random kinetic energy if PLAW(4)<0., + or extract random momentum if PLAW(4)>=0.. + +Sept 2007, Bari + +POWER LAW SPECTRA GENERATION: + It is now possible to select a power law spectra. Using the + datacard: PLAW, is possible to pass the spectral index and the + kinetic energy range. In PAWL(1) the user can pass the signed power + law index, i.e. PAWL -2.5, will generate a power law spectrum + according to exp(-2.5). Please note that the range is in kinetic + energy. In gukine the generated particle charateristics will be + retrived and the mass calculated from this infos. The KINE data + card assigned momentum will be overriden as well. + +PATHLENGHT CALCULATION MODIFIED IN TOF: + The pathlength calulation saved in the TOF hit structure has now + been modified. Insted of the pathlenght inside each paddle the + average pathlenght from the generation vertex to the hit in the + paddle, is saved. In this way it will be possible to calculate the + beta using this variable. + +ERROR IN MATERIAL ASSIGNMENT FOR NEUTRON DETECTOR: + The material assigned to the NDTU volume was wrong. It was + scintillator, IRON has been asigned now. + +March 2007, Bari + +CALO GEOMETRY: + An error was found in the original document describing PAMELA + geometry. The shift along X axis of the X view planes are swapped, + i.e. they don't start with X=+0.05 but with X=-0.05. The sequence + is: 0.05m -0.05, 0.05, -0.05 and so on. + +MAKEFILE FIXED: + Option -fdebug_kludge has been discharged + +GASATT WRNING FIXED: + There was a warning message from a graphic routine. The GSATT calls has been commented out. + + +December 2006, Bari + +CARD GEOMETRY: +Positions, dimensions and quotes have been updated. Indeed, the aluminium table +under S1, legs included, and the magnetic screen between the CARD and S2 have +been added. + +#------------------------------------------------------------------------- + +November 2006, Bari + + +S4 GEOMETRY UPDATED: +The x and y dimensions have been corrected because they were wrong. +A container of aluminum has been added having walls of thickness equal +to 0.1 cm. + +TOF POSITIONS DEFINITELY UPDATED: +S2 and S1 have been positioned again, after the simulation of the top plate. +Before now, the positions were put ad hoc based on the positions given by +Sergio Ricciarini in the document 'Main geometrical parameters of the PAMELA +sub-detectors' released by O. Adriani, L. Bonechi, E. Mocchiutti and +S. Ricciarini on the 20th of December 2005.. These positions were +lightly approximated. + +ND GEOMETRY COMPLETED. +The He3 and the polyethylene (CH2) have been added to fill +the volumes NDTI and NDPB in the neutron detector. Leonov gave us the +details of these materials. + +TOP PLATE GEOMETRY ADDED. +The top plate geometry has been simulated as a N2 box that has the same +dimensions of CATA. Inside it, at the right positions, there are S2 and CATA. +The real top plate is a rectangular window having a thickness of 0.5 cm +(volume TPTL emptied by the volume TPCV). This volume is filled with a mixture +of titanium (90%), aluminum (6%) and vanadium (4%). Other parts of titanium +are the four volumes TPTU, which are positioned at the corners of TPTL at the +same height of CATA, and the four volumes TPTM, which are positioned at the +corners of TPTL, at the same height of S2. +The geometry of the top plate is based on a simplified version of the CAD +drawings. + +October 2006, Bari + +The geometry of the neutron detector has been updated. Some dimensions of some +volumes have been corrected and the aluminum cover and the aluminum boxes to +put the cables have been added. A new material, the cadmium, has been also +defined. + +September 2006, Bari + +SPHE and ND data card bugs fixed: the definition of the ND data card, +missing in the subroutine gpgeo.F, has been added; the meaning of the SPHE +data card has been changed. Before the correction the data card: +NDET 'SPHE' was used to delete the spherical top shell to substitute it with +a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA. + + +June 2006, Bari + +The center of the scintillator planes S22Y (variable ZPAMS22Y in gpdgeo.inc) +and S12X (variable ZPAMS12X in gpdgeo.inc) has been positioned at the +nominal height as measured in PAMELA (See the document: "Main geometrical +parameters of the PAMELA sub-detectors" by O. Adriani, L. Bonechi, +E. Mocchiutti, S. Ricciarini, 20 December 2005). Follows that the positions +of S21Y and S12X are higher than those in the cited document due to the fact +that in GPAMELA the thickness of the mylar has been considered while in the +document it has been neglected. + + +May 2006, Bari & Tor Vergata + +GIGANTIC RESONANCE FOR NEUTRON DETECTOR ADDED + + Routines to simulate the gigantic resonance of gammas in Tungsten + have been added. The GPGIG routine is called in GUSTEP if a gamma + enter the calorimeter absorber. This is the steering routine to + simulate the production of neutrons from gigantic resonance. It + does checks on STEP lenght. If the range is smaller than the other + selected for that step, it does generate the neutron and stops the + gamma. Please note that the neutron has now a new particle + number. This is to tag the gigantic resonance neutrons. + + +May 2006, Bari & Florence + +CAL HIT STRUCTURE BUGS FIXED + + The maximum number of hit is now different for the two hit + structures: CALST and CALI. Vectors inizialization and HBOOK + ntple booking have been updated. The GPDCAL routine has been fixed + so to handle the case in wich hits stored are more than the maximum + number of hit. + In this case in the ntple up to the maximum number of hits will be stored. + +April 2006, Bari + +TOF GEOMETRY AND POSITIONS UPDATED AND NEW MIXTURES ADDED + + The TOF geometry has been modified. The following boxes have been + added: POL1, POL2 and POLY made of polystyrene, S11M, S12M, S21M, + S22M, S31M and S32M made of mylar, S1A, S2A and S3 made of air and + S1 and S2 made of aluminum. Each scintillator paddle has been put + in his mylar box and the other materials: air, polystyrene, and + aluminum have been added at their nominal positions. According to + Naples people the araldite glue has been simulated has an air + gap. For this work two new materials: the Mylar (MYLAR) and the + polystyrene (POLYSTYRENE) with a density of 35 g/l have been + defined as a mixture. The positions of the three bottom + scintillator planes that contain respectively the S12X, S22Y and + S32X paddles have been regulated according on their official + positions in PAMELA. Mar 2006, Bari @@ -67,7 +297,7 @@ is set to true in GPDAT if the data card: HPAK, is set to 'GCAL'. The gpkey.inc has been update accordingly. - + NEUTRON DETECTOR ADDED. NEW DIR: GPND The neutron detector has been added. At the moment it is just the @@ -84,7 +314,7 @@ http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html This package should be distributed with the GEANT library but is not up to date. I did download the latest release and stored into - gpcalor directory of the gpamela tree. + gpcalor directory of the gpamela tree. Then I did clean up the code substituting the explicit inclusion of the commons with a #include cpp directive. In parallel I did extract the commons to include files having the same common name. I @@ -95,7 +325,7 @@ NOTE: There are still problem due to different common sizes. In particular the common MICFIL is maller in the geant library libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are - 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. NOTE2: The auxiliary files with the cross sections: chetc.dat.gz and xsneut.dat.gz, have been added to the aux directory and moved to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will @@ -103,24 +333,32 @@ searched there at first, then in the working directory. A fool proof policy has to be implemented to avoid problem with synchronization fo these files. - -The GCALOR package + +The GCALOR package + +June 2005, Bari + +TOF SCINTILLATOR PADDLES UPDATED + + The dimensions and the number of the scintillator paddles for each + TOF planes have been updated. + May 2005, Bari -Some updates on the latest modification done in the past year. +Some updates on the latest modification done in the past year. -NEW DATA CARD ADDED: HFSF +NEW DATA CARD ADDED: HFSF To define a policy for the random number initial seeds definition. Using this card is possible to override GEANT seeds defined via NRDM card. The policy is selected according to the values: - - 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 defined file or the default file: INPUTSEED.DAT - - - 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 file or the default file: INPUTSEED.DAT The case 1 must be used in case the user needs to reproduce the @@ -129,7 +367,7 @@ binary file and pass the filename to the program using the *FLSF data card. In case the user file is not specified the default INPUTSEED.DAT will be used. - + The case 2 must be used in case the user needs to chain several GPAMELA run and likes to be sure he is starting the random generator using the right sequence. In this case the user must @@ -147,8 +385,8 @@ filled with his own seeds. For this purpose an utility program: writeseeds.f, has been added in the aux directory. In case the *FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will - be used as LUN and file name respectively. - + be used as LUN and file name respectively. + NEW DATA CARD ADDED: *LSFI Using this card the user can specify the logical unit and name of @@ -156,9 +394,9 @@ run. This file is a FORTRAN binary one. This file can be used as input one specifying it in the *FSFI data card of the next GPAMELA run. In case the *LSFI card is not specified the default values: 26 - and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN + and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN and file name respectively. - + NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f These new programs have been added in the aux directory. Using these a @@ -191,15 +429,15 @@ has been revied and new sets of calibration constants calculated comparing the beam test data with the GPAMELA results. The new constants are about 3% larger than the previous ones. - -TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI - + +TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI + Thanks to Lorenzo Bonechi for the drawings and explanations. Now the hybrd cards have been put into the simulation and the geometry updated considering the dead zones in the silicon detectors. The hybrid zone has been simulated as well. At the moment the hybrid is simulated as a G10 plates. The full height of the tracker magnet has been - reviewed as well. + reviewed as well. The tracker ladder is now simulated inside a nitrogen box: TPAS, placed inside an aluminum frame: TRPB. Each silicon ladder has been @@ -207,21 +445,21 @@ smaller silicon detector: TPAI, has been placed inside the larger silicon block TRSL. In this way the subdivided silicon ladder can be upgraded with an indipendend roto-translation for each sensor. - + The TRPB aluminum frame has been enlarged to fit the external magnet canister frame. - - The last plane has been flipped with a 180 degree rotation around - the X axis. - + + The last plane has been flipped with a 180 degree rotation around + the X axis. + TRACKER HIT STRUCTURE REVIEWED Taking into account the new version of the tracker geometry, the hit - structure for this detector has been revied. + structure for this detector has been revied. CALORIMETER GEOMETRY REVIEWED - Marco Albi reviewed the calorimeter dimentions and positioning. + Marco Albi reviewed the calorimeter dimentions and positioning. 29 March 2004, Bari @@ -235,7 +473,7 @@ the GARFIELD code for details). Follows that different runs produce different results without changing the initial seeds. To have identical results in different runs, the GARFIELD program has to start typing the - noRNDM_initialisation switch. To avoid of specifying this switch + noRNDM_initialisation switch. To avoid of specifying this switch by the user, the GARFIELD package has been upgraded with a patch. In this way the problem is partially solved because, now, the initial seeds of the random generators @@ -243,7 +481,7 @@ activated by the user for changing the initial seeds in the GPAMELA program. Work is in progress for a more general correction of this problem. Please, use the updated GARFIELD code released with the CVS version v4r1 - to fix this problem. + to fix this problem. RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS. @@ -256,7 +494,7 @@ the GPSED common when the set and the detector ZEBRA banks are filled during a run. In general the values of the set and detector addresses depend on the number of active detectors in a given run. ISCAL=12 and - 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. 9 December 2003, Bari @@ -300,7 +538,7 @@ -ICAPLANE(i) contains the number of hit plane; -ICASTRIP(i) contains the number of hit strip; -ICAMOD(i) can assume different values based on the number of times and - positions in which a `long strip' has been hit. + positions in which a `long strip' has been hit. -ENESTRIP(i) contains the deposited energy in the hit strip; where i is the number of hit (1