--- gpamela/history/v_100.txt 2002/12/05 17:27:59 3.3 +++ gpamela/history/v_100.txt 2005/12/14 03:16:08 3.8 @@ -1,7 +1,22 @@ # -# $Id: v_100.txt,v 3.2 2002/12/05 10:17:42 pamela Exp $ +# $Id: v_100.txt,v 3.7 2005/10/18 08:24:35 cafagna Exp $ # # $Log: v_100.txt,v $ +# 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 # @@ -20,31 +35,326 @@ #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 + +Nov 2005, Bari + +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 +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