--- gpamela/history/v_100.txt 2006/11/30 13:22:35 3.23 +++ gpamela/history/v_100.txt 2007/10/02 17:58:32 3.27 @@ -1,12 +1,15 @@ +# 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.22 2006/11/28 11:43:41 pam-ba Exp $ +# $Id: v_100.txt,v 3.26 2007/03/30 15:24:54 cafagna Exp $ # # $Log: v_100.txt,v $ -# Revision 3.22 2006/11/28 11:43:41 pam-ba -# History modified +# 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.21 2006/11/28 10:26:15 pam-ba -# S3 positioning completed +# 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 @@ -80,39 +83,70 @@ #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 -30 November 2006, Bari F. Cafagna - -BUG FIXED - - For some reason there was a conflict using FLUKA and the GARFILD - interface for the TRD. The file TRD_ATT.TXT was treated as flukaerr - file by the compiler. It is not clear why. The code was correct, - the logical units and names were not inconflict, but the compiler - assigned to them the same internal number. It may be a problem - related to the optimization phase. I did close the file just after - the read-out in GPDAT.F and commented out the closing in UGLAST.F +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. +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 by hand based on the positions given by +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 because the thicknesses of the mylar and of the glue were -not considered in that document. Also the height of S3 has been modified as -suggested from the Naples PAMELA collaborators. +lightly approximated. ND GEOMETRY COMPLETED. The He3 and the polyethylene (CH2) have been added to fill @@ -120,14 +154,16 @@ details of these materials. TOP PLATE GEOMETRY ADDED. -The top plate geometry has been simulated as a rectangular window filled -with a mixture of titanium (90%), aluminum (6%) and vanadium (4%). It has been -put in a N2 box having the same planar dimensions of CATA. In this box there -are also S2 and CATA at the right positions. Then the top plate geometry has -been completed with other eight small volumes of titanium positioned at the -corners of the N2 box. Four of them are at the same height of CATA, the -other four are at the same height of S2. The simulated geometry of the top -plate is based on a simplified version of the CAD drawings. +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 @@ -142,10 +178,10 @@ 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. +a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA. -June 2006, Bari +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 @@ -154,7 +190,7 @@ 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. +document it has been neglected. May 2006, Bari & Tor Vergata @@ -180,7 +216,7 @@ 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. + In this case in the ntple up to the maximum number of hits will be stored. April 2006, Bari @@ -223,7 +259,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 @@ -240,7 +276,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 @@ -251,7 +287,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 @@ -259,9 +295,9 @@ 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 @@ -272,19 +308,19 @@ 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 @@ -293,7 +329,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 @@ -311,8 +347,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 @@ -320,9 +356,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 @@ -355,15 +391,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 @@ -371,21 +407,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 @@ -399,7 +435,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 @@ -407,7 +443,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. @@ -420,7 +456,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 @@ -464,7 +500,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