--- gpamela/history/v_100.txt 2002/07/11 16:01:59 3.1.1.1 +++ gpamela/history/v_100.txt 2006/11/30 13:22:35 3.23 @@ -1,7 +1,73 @@ # -# $Id$ +# $Id: v_100.txt,v 3.22 2006/11/28 11:43:41 pam-ba Exp $ +# +# $Log: v_100.txt,v $ +# Revision 3.22 2006/11/28 11:43:41 pam-ba +# History modified +# +# Revision 3.21 2006/11/28 10:26:15 pam-ba +# S3 positioning completed +# +# 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. +# +# Revision 3.8 2005/12/14 03:16:08 cafagna +# Neutron detector added. Geometry and GPCALOR package +# +# Revision 3.7 2005/10/18 08:24:35 cafagna +# History updated +# +# Revision 3.6 2005/07/25 11:53:21 cafagna +# Several updates. See history for details +# +# Revision 3.5 2004/04/06 10:33:46 pamela +# NON-REPRODUCIBILITY problem of a GPAMELA RUN fixed; bug found and fixed filling in the hit structure of the calorimeter +# +# Revision 3.4 2003/12/17 11:32:50 pamela +# CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation +# +# Revision 3.3 2002/12/05 17:27:59 pamela +# New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated +# +# Revision 3.2 2002/12/05 10:17:42 pamela +# Update CAS and CALO geometries and positions. Makefile updated as well +# +# Revision 3.1.1.1 2002/07/11 16:01:59 cafagna +# First GPAMELA release on CVS # -# $Log$ # #CMZ : 3.00/00 11/02/2002 20.05.23 by Unknown #CMZ : 2.03/00 06/11/2000 02.14.56 by Francesco Cafagna @@ -14,7 +80,489 @@ #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 + + + +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 by hand 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. + +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 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. + +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 + +GEN DATA CARD UPDATED + + To enable generation on a surface perpendicular to the XY plane, + GEN gata card has been updated addingh a new parameter: ZDGEN. This is + the dimension, along Z axis , of the generation surface. The Z + position will be randomply chosen according to: Z= ZDGEN*RNDM_NUMBER + + ZGEN, i.e. Z= GEN(6)*RNDM_NOMBER + GEN(3). + +Nov 2005, Bari + +GUHADR AND GUPHAD UPDATED + + To use GCALOR package the hadronic routines have been updated. The + inizialization routine call CALSIG, while the other calls GCALOR. + +NEW GPKEY ADDED: GPCALOR + + This logical has been added to enable the GCALOR package. This flag + 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 + geometry. The directory structure of the repository has been + updated as well. Dimensions has been taken from picture and + literature. A full upgrade to the drawing is needed. + +GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR, GPCALORDES + + GCALOR package contins the CALOR simulation code and an interface + to use it in GEANT. The important feature for us is the usage of + the MICAP code. This is facused on the low energy neutron + simulation. for details see: + 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. + 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 + did store the include files into a newly created directory: + gpcalordes. + The Makefile has been updated accordingly. + Please note that to avoid conflict with CRENLIB distribution the gcalor source file has been named gpcalor.F + 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. + 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 + look for CERN_ROOT environment variable. If found files are + 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 + +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. + +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 + defined file or the default file: INPUTSEED.DAT + + - 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 + random chain of a previous run. In this case the user can save the + initial seeds, used in the run he would like to reproduce, in a + 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 + specify an input binary file using the *FLSF data card, otherwise + the INPUTSEED.DAT file will be used. + +NEW DATA CARD ADDED: *FSFI + + Using this card the user can specify the logical unit and name of + the file storing the initial seeds to be used to initialize the + random number generator. This file must be a FORTRAN binary one + storing four integer numbers. The first two are the number to be + used in the case: HFSF=1, the other two will be used in the case: + HFSF=2. This file can be one created by GPAMELA or by the user + 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. + +NEW DATA CARD ADDED: *LSFI + + Using this card the user can specify the logical unit and name of + the file storing the first and last seeds used in the GPAMELA + 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 file name respectively. + +NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f + + These new programs have been added in the aux directory. Using these a + user defined seed file can be created and re-read. + +NEW VOLUMES ADDED: MSHE, BSPH; PRESSURIZED CONTAINER ADDED + + Alexey Bakaldin, in MEPHI, did add the PAMELA pressurized container to + the simulation. He did defined new volumes filled with aluminum and + placed inside the mother volume. Positions have been fine tuned by + Marialuigia Ambriola and compared to the CAD drawings. + Two new volumes have been added to simulate the container: + - MSHE, a tube simulating the middle part of the container + - BSPH, the spherical bottom part of the container + + To better simulate the upper part the SHEL volume has been modified + into a cone. Dimentions of the top cover: TSPH, have been modified + accordingly. + +DETECTOR POSITIONS REVIEWED + + All detector Z positions have been reviewd to fit into the + simulated pressurized container. + +TRD GEOMETRY AND CALIBRATION REVIEWD + + The TRD geometry has been deeply reviewed. Using the CAD drawings + the carbon fiber frames have been simulated and radiator dimentions + corrected. For this reason the calibration done on the beam tests + 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 + + 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. + + The tracker ladder is now simulated inside a nitrogen box: TPAS, + placed inside an aluminum frame: TRPB. Each silicon ladder has been + simulated using two silicon blocks: TRSL, into each of this block a + 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. + +TRACKER HIT STRUCTURE REVIEWED + + Taking into account the new version of the tracker geometry, the hit + structure for this detector has been revied. + +CALORIMETER GEOMETRY REVIEWED + + Marco Albi reviewed the calorimeter dimentions and positioning. + + +29 March 2004, Bari + +NON-REPRODUCIBILITY PROBLEM OF A GPAMELA RUN FIXED. + The non-reproducibility of a GPAMELA run was due to the random number + initialization in the GARFIELD code. In GARFIELD by default, the initial + seeds of the random number generators are always the same while the random + number generators are called a given number of times (determined by the + hour of the day) during the initialization phase (see init.f subroutine in + 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 + 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 + in GARFIELD will be always the same even if the RNDM GEANT data card is + 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. + + +RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS. + The obsolete RNDM random number generator has been replaced by the GEANT + GRNDN routine in the gpxtr.F subroutine and in the npoiss.F function. + +BUG FOUND AND FIXED: the set and detector calorimeter addresses (ISCAL + and IDCASI variables) used in GUTREV were respectively set to a fixed + values of 12 and 1. The correct values of these variables are stored in + 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. + +9 December 2003, Bari + + CALORIMETER SIMULATION completed! The update of the geometry and of the + special tracking parameters and the tuning of the calorimeter have been + successfully done. A great quantity of simulated data have been produced + in the calorimeter for different particles (muons, electrons and pions) + and momenta (5 and 40 GeV/c) and the output data have been analyzed. The + distributions of the total energy deposited in the calorimeter and the + total number of strips hit have been compared with the respective + distributions produced by the Trieste's tuned standalone Monte Carlo + simulation program of the PAMELA calorimeter. The accord between the + two simulations is excellent. Many thanks to Mirko for his collaboration. + + Working in progress on TRD. The GARFIELD interface to the HEED program is not + optimized to track particle with a charge greater than one and photons. The + program print a warning message to advise the user when it is the case. + +18 April 2003, Bari + + The buffer size of each column of the GPAMELA Ntuple has been increased to + 4096 and set equal to the record length, defined by a call to the HROPEN + routine. + Also the length of the common /PAWC/ (parameter NWPAW) has been increased + to 1.34E8, according to the rule that it has to be larger than the number + of columns times the buffer size. + +10 April 2003, Bari + + The variables in the HIT STRUCTURE of the CALORIMETER and their way to be + filled have been changed according to the electronics system of the real + detector. In fact, because each silicon detector (module) consists of + 32 strips and each strip is connected to those belonging to the two detectors + of the same row (or column) for forming 24 cm long strips, the sum of the + deposited energies in the strips forming a `long strip' is now calculated for + each event (gpucal.F subroutine) and it is stored in a hit only at the + end of the event (gutrev.F subroutine). + The output variables of the GPAMELA en-tuple are then filled in the vectors + ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL), + by a call to the GPDCAL subroutine: + -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. + -ENESTRIP(i) contains the deposited energy in the hit strip; + where i is the number of hit (1