/[PAMELA software]/gpamela/history/v_100.txt
ViewVC logotype

Diff of /gpamela/history/v_100.txt

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.1 by cafagna, Thu Jul 11 16:01:59 2002 UTC revision 3.10 by cafagna, Mon Apr 10 11:07:43 2006 UTC
# Line 1  Line 1 
1    
2    #
3    # $Id: v_100.txt,v 3.9 2005/12/14 03:34:40 cafagna Exp $
4    #
5    # $Log: v_100.txt,v $
6    # Revision 3.9  2005/12/14 03:34:40  cafagna
7    # An update of the history and inform readme files.
8    #
9    # Revision 3.8  2005/12/14 03:16:08  cafagna
10    # Neutron detector added. Geometry and GPCALOR package
11    #
12    # Revision 3.7  2005/10/18 08:24:35  cafagna
13    # History updated
14    #
15    # Revision 3.6  2005/07/25 11:53:21  cafagna
16    # Several updates. See history for details
17    #
18    # Revision 3.5  2004/04/06 10:33:46  pamela
19    # NON-REPRODUCIBILITY problem of a GPAMELA RUN fixed; bug found and fixed filling in the hit structure of the calorimeter
20    #
21    # Revision 3.4  2003/12/17 11:32:50  pamela
22    # CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation
23  #  #
24  # $Id$  # Revision 3.3  2002/12/05 17:27:59  pamela
25    # New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated
26    #
27    # Revision 3.2  2002/12/05 10:17:42  pamela
28    # Update CAS and CALO geometries and positions. Makefile updated as well
29    #
30    # Revision 3.1.1.1  2002/07/11 16:01:59  cafagna
31    # First GPAMELA release on CVS
32  #  #
 # $Log$  
33  #  #
34  #CMZ :  3.00/00 11/02/2002  20.05.23  by  Unknown  #CMZ :  3.00/00 11/02/2002  20.05.23  by  Unknown
35  #CMZ :  2.03/00 06/11/2000  02.14.56  by  Francesco Cafagna  #CMZ :  2.03/00 06/11/2000  02.14.56  by  Francesco Cafagna
# Line 14  Line 42 
42  #CMZ :  1.00/03 30/04/96  12.23.59  by  Francesco Cafagna  #CMZ :  1.00/03 30/04/96  12.23.59  by  Francesco Cafagna
43  #CMZ :  1.00/02 05/04/96  15.31.25  by  Francesco Cafagna  #CMZ :  1.00/02 05/04/96  15.31.25  by  Francesco Cafagna
44  #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
45  #-- Author :    Francesco Cafagna   28/11/95  #-- Author :    Francesco Cafagna   28/11/95
46    
47    Mar 2006, Bari
48    
49    GEN DATA CARD UPDATED
50    
51       To enable generation on a surface perpendicular to the XY plane,
52       GEN gata card has been updated addingh a new parameter: ZDGEN. This is
53       the dimension, along Z axis , of the generation surface. The Z
54       position will be randomply chosen according to: Z= ZDGEN*RNDM_NUMBER +
55       ZGEN, i.e. Z= GEN(6)*RNDM_NOMBER + GEN(3).
56    
57    Nov 2005, Bari
58    
59    GUHADR AND GUPHAD UPDATED
60    
61       To use GCALOR package the hadronic routines have been updated. The
62       inizialization routine call CALSIG, while the other calls GCALOR.
63    
64    NEW GPKEY ADDED: GPCALOR
65    
66       This logical has been added to enable the GCALOR package. This flag
67       is set to true in GPDAT if the data card: HPAK, is set to
68       'GCAL'. The gpkey.inc has been update accordingly.
69    
70      
71    NEUTRON DETECTOR ADDED. NEW DIR: GPND
72    
73       The neutron detector has been added. At the moment it is just the
74       geometry. The directory structure of the repository has been
75       updated as well. Dimensions has been taken from picture and
76       literature. A full upgrade to the drawing is needed.
77    
78    GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR, GPCALORDES
79    
80       GCALOR package contins the CALOR simulation code and an interface
81       to use it in GEANT. The important feature for us is the usage of
82       the MICAP code. This is facused on the low energy neutron
83       simulation. for details see:
84       http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html
85       This package should be distributed with the GEANT library but is
86       not up to date. I did download the latest release and stored into
87       gpcalor directory of the gpamela tree.
88       Then I did clean up the code substituting the explicit inclusion of
89       the commons with a #include cpp directive. In parallel I did
90       extract the commons to include files having the same common name. I
91       did store the include files into a newly created directory:
92       gpcalordes.
93       The Makefile has been updated accordingly.
94       Please note that to avoid conflict with CRENLIB distribution the gcalor source file has been named gpcalor.F
95       NOTE: There are still problem due to different common sizes. In
96       particular the common MICFIL is maller in the geant library
97       libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are
98       present and linked using a wrong version of the common. This still needs to be debuged.
99       NOTE2: The auxiliary files with the cross sections: chetc.dat.gz
100       and xsneut.dat.gz, have been added to the aux directory and moved
101       to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will
102       look for CERN_ROOT environment variable. If found files are
103       searched there at first, then in the working directory. A fool
104       proof policy has to be implemented to avoid problem with
105       synchronization fo these files.
106      
107    
108    The GCALOR package
109    May 2005, Bari
110    
111    Some updates on the latest modification done in the past year.
112    
113    NEW DATA CARD ADDED: HFSF
114    
115       To define a policy for the random number initial seeds
116       definition. Using this card is possible to override GEANT seeds
117       defined via NRDM card. The policy is selected according to the
118       values:
119    
120       - 1: The seeds are initialized to the initial values found in a user
121            defined file or the default file: INPUTSEED.DAT
122      
123       - 2: The seeds are initialized to the final values found in a user defined
124            file or the default file: INPUTSEED.DAT
125    
126       The case 1 must be used in case the user needs to reproduce the
127       random chain of a previous run. In this case the user can save the
128       initial seeds, used in the run he would like to reproduce, in a
129       binary file and pass the filename to the program using the *FLSF
130       data card. In case the user file is not specified the default
131       INPUTSEED.DAT will be used.
132      
133       The case 2 must be used in case the user needs to chain several
134       GPAMELA run and likes to be sure he is starting the random
135       generator using the right sequence. In this case the user must
136       specify an input binary file using the *FLSF data card, otherwise
137       the INPUTSEED.DAT file will be used.
138    
139    NEW DATA CARD ADDED: *FSFI
140    
141       Using this card the user can specify the logical unit and name of
142       the file storing the initial seeds to be used to initialize the
143       random number generator. This file must be a FORTRAN binary one
144       storing four integer numbers. The first two are the number to be
145       used in the case: HFSF=1, the other two will be used in the case:
146       HFSF=2. This file can be one created by GPAMELA or by the user
147       filled with his own seeds. For this purpose an utility program:
148       writeseeds.f, has been added in the aux directory.  In case the
149       *FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will
150       be used as LUN and file name respectively.
151      
152    NEW DATA CARD ADDED: *LSFI
153    
154       Using this card the user can specify the logical unit and name of
155       the file storing the first and last seeds used in the GPAMELA
156       run. This file is a FORTRAN binary one. This file can be used as
157       input one specifying it in the *FSFI data card of the next GPAMELA
158       run.  In case the *LSFI card is not specified the default values: 26
159       and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN
160       and file name respectively.
161      
162    NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f
163    
164       These new programs have been added in the aux directory. Using these a
165       user defined seed file can be created and re-read.
166    
167    NEW VOLUMES ADDED: MSHE, BSPH; PRESSURIZED CONTAINER ADDED
168    
169       Alexey Bakaldin, in MEPHI, did add the PAMELA pressurized container to
170       the simulation. He did defined new volumes filled with aluminum and
171       placed inside the mother volume. Positions have been fine tuned by
172       Marialuigia Ambriola and compared to the CAD drawings.
173       Two new volumes have been added to simulate the container:
174       - MSHE, a tube simulating the middle part of the container
175       - BSPH, the spherical bottom part of the container
176    
177       To better simulate the upper part the SHEL volume has been modified
178       into a cone. Dimentions of the top cover: TSPH, have been modified
179       accordingly.
180    
181    DETECTOR POSITIONS REVIEWED
182    
183       All detector Z positions have been reviewd to fit into the
184       simulated pressurized container.
185    
186    TRD GEOMETRY AND CALIBRATION REVIEWD
187    
188       The TRD geometry has been deeply reviewed. Using the CAD drawings
189       the carbon fiber frames have been simulated and radiator dimentions
190       corrected. For this reason the calibration done on the beam tests
191       has been revied and new sets of calibration constants calculated
192       comparing the beam test data with the GPAMELA results. The new
193       constants are about 3% larger than the previous ones.
194      
195    TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI
196      
197       Thanks to Lorenzo Bonechi for the drawings and explanations. Now the
198       hybrd cards have been put into the simulation and the geometry updated
199       considering the dead zones in the silicon detectors. The hybrid zone
200       has been simulated as well. At the moment the hybrid is simulated as
201       a G10 plates. The full height of the tracker magnet has been
202       reviewed as well.
203    
204       The tracker ladder is now simulated inside a nitrogen box: TPAS,
205       placed inside an aluminum frame: TRPB. Each silicon ladder has been
206       simulated using two silicon blocks: TRSL, into each of this block a
207       smaller silicon detector: TPAI, has been placed inside the larger
208       silicon block TRSL. In this way the subdivided silicon ladder can
209       be upgraded with an indipendend roto-translation for each sensor.
210      
211       The TRPB aluminum frame has been enlarged to fit the external
212       magnet canister frame.
213      
214       The last plane has been flipped with a 180 degree rotation around
215       the X axis.
216      
217    TRACKER HIT STRUCTURE REVIEWED
218    
219       Taking into account the new version of the tracker geometry, the hit
220       structure for this detector has been revied.
221    
222    CALORIMETER GEOMETRY REVIEWED
223    
224       Marco Albi reviewed the calorimeter dimentions and positioning.
225    
226    
227    29 March 2004, Bari
228    
229    NON-REPRODUCIBILITY PROBLEM OF A GPAMELA RUN FIXED.
230       The non-reproducibility of a GPAMELA run was due to the random number
231       initialization in the GARFIELD code. In GARFIELD by default, the initial
232       seeds of the random number generators are always the same while the random
233       number generators are called a given number of times (determined by the
234       hour of the day) during the initialization phase (see init.f subroutine in
235       the GARFIELD code for details). Follows that different runs produce
236       different results without changing the initial seeds. To have identical
237       results in different runs, the GARFIELD program has to start typing the
238       noRNDM_initialisation switch. To avoid of specifying this switch
239       by the user,
240       the GARFIELD package has been upgraded with a patch. In this way the problem
241       is partially solved because, now, the initial seeds of the random generators
242       in GARFIELD will be always the same even if the RNDM GEANT data card is
243       activated by the user for changing the initial seeds in the GPAMELA program.
244       Work is in progress for a more general correction of this problem.
245       Please, use the updated GARFIELD code released with the CVS version v4r1
246       to fix this problem.  
247    
248    
249    RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS.
250       The obsolete RNDM random number generator has been replaced by the GEANT
251       GRNDN routine in the gpxtr.F subroutine and in the npoiss.F function.
252    
253    BUG FOUND AND FIXED: the set and detector calorimeter addresses (ISCAL
254       and IDCASI variables) used in GUTREV were respectively set to a fixed
255       values of 12 and 1. The correct values of these variables are stored in
256       the GPSED common when the set and the detector ZEBRA banks are filled
257       during a run. In general the values of the set and detector addresses
258       depend on the number of active detectors in a given run. ISCAL=12 and
259       IDCASI=1 are only right when all the detectors of GPAMELA are active.
260    
261    9 December 2003, Bari
262    
263       CALORIMETER SIMULATION completed! The update of the geometry and of the
264       special tracking parameters and the tuning of the calorimeter have been
265       successfully done. A great quantity of simulated data have been produced
266       in the calorimeter for different particles (muons, electrons and pions)
267       and momenta (5 and 40 GeV/c) and the output data have been analyzed. The
268       distributions of the total energy deposited in the calorimeter and the
269       total number of strips hit have been compared with the respective
270       distributions produced by the Trieste's tuned standalone Monte Carlo
271       simulation program of the PAMELA calorimeter. The accord between the
272       two simulations is excellent. Many thanks to Mirko for his collaboration.
273    
274       Working in progress on TRD. The GARFIELD interface to the HEED program is not
275       optimized to track particle with a charge greater than one and photons. The
276       program print a warning message to advise the user when it is the case.
277    
278    18 April 2003, Bari
279    
280       The buffer size of each column of the GPAMELA Ntuple has been increased to
281       4096 and set equal to the record length, defined by a call to the HROPEN
282       routine.
283       Also the length of the common /PAWC/ (parameter NWPAW) has been increased
284       to 1.34E8, according to the rule that it has to be larger than the number
285       of columns times the buffer size.
286    
287    10 April 2003, Bari
288    
289       The variables in the HIT STRUCTURE of the CALORIMETER and their way to be
290       filled have been changed according to the electronics system of the real
291       detector. In fact, because each silicon detector (module) consists of
292       32 strips and each strip is connected to those belonging to the two detectors
293       of the same row (or column) for forming 24 cm long strips, the sum of the
294       deposited energies in the strips forming a `long strip' is now calculated for
295       each event (gpucal.F subroutine) and it is stored in a hit only at the
296       end of the event (gutrev.F subroutine).
297       The output variables of the GPAMELA en-tuple are then filled in the vectors
298       ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL),
299       by a call to the GPDCAL subroutine:
300       -ICAPLANE(i) contains the number of hit plane;
301       -ICASTRIP(i) contains the number of hit strip;
302       -ICAMOD(i) can assume different values based on the number of times and
303                  positions in which a `long strip' has been hit.
304       -ENESTRIP(i) contains the deposited energy in the hit strip;
305       where i is the number of hit (1<i<4224).
306       Note that in the calorimeter each hit is filled at the end of the event and
307       that there is a hit for each `long strip' hit from
308       the particle. This use of the hit structure is different for the other
309       detectors and it has been considered to avoid a too big number of hit in the
310       calorimeter due to the showers. Follows that NTHCAL, which is the
311       max number of hit in the calorimeter, is equal to 4224, the total
312       number of `long strips'. So, for each event, the real number of hit will
313       be less or equal to 4224.
314       ICAMOD(i) is an additional information that does not exist in the real
315       detector: if the strip i (i=1,32) of the module 1 or 2 or 3
316       is hit, the value of ICAMOD(i) is respectively incremented of 1, 100, 10000.
317       Analogously it is done, if it is the strip j (j=33,64) of the modules 4, 5
318       and 6 or if it is the strip k (k=65,96) of the modules 7, 8 and 9.
319       For example if we consider the hit 1 of an event, we could read:
320       ICASTRIP(1)=30, ICAPLANE(1)=21, ENESTRIP(1)=0.5E-03 and ICAMOD(1)=10001.
321       It means that the hit 1 contains the information that in the strip 30 of the
322       plane 21 has been deposited a total energy of 0.5E-03 GeV. In addition the
323       `long strip 30' has been hit two times, one in the first module and the
324       other in the third one.
325    
326       The energy deposited in the calorimeter is calculated in GeV.
327    
328       To store the hits in the calorimeter the subroutine GSAHIT is used instead of
329       GSCHIT.
330    
331       To retrieve the hit structure the call to the routine GPRHIT is done instead
332       of a call to the GFHITS subroutine.
333    
334    25 February 2003, Bari
335    
336    BUG found:
337       DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP,
338       DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2,
339       DCUTEW, DCUTEXE variables missed in the commons: gpaer.inc, gpal.inc, gpce.inc,
340       gpcp.inc, gpfe.inc, gpg10c.inc, gpg10.inc, gpkap.inc, gpn2g.inc, gproa.inc,
341       gpscin.inc (obsolete), gpscint.inc, gpsica.inc, gptrad.inc, gpw2.inc, gpw.inc,
342       gpxe.inc, gpdaer.inc, gpdal.inc, gpdce.inc, gpdcp.inc, gpdfe.inc, gpdg10c.inc,
343       gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc,
344       gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc.
345       They have been added in these commons and they have been initialized in the
346       GPSTM subroutine.
347    
348       Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C
349       in the subroutines gpsica.F, gpw2.F, gpw.F, gpce.F and gpg10c.F. They were
350       suggested by Mirko Boezio.
351    
352       Updated the value of the absorption length for silicon in the calorimeter
353       and tracker although this parameter is ignored by GEANT. For this reason
354       it was equal to the radiation length.
355    
356       Updated the relative positions of the calorimeter planes. The corrected
357       shifting are:
358    
359       first view: (Dxo,Dyo)=(0.10,0.05) cm
360       second view: (Dxo,Dyo)=(-0.05,0.10) cm
361       third view: (Dxo,Dyo)=(-0.10,-0.05) cm
362       fourth view: (Dxo,Dyo)=(0.05,-0.10) cm
363    
364    4 November 2002, Bari
365    
366    CAS detectors distances modified
367    
368       The distances between the CAS detectors have been modified based on the
369       latest CAD drawings.
370    
371    2 November 2002, Bari
372    
373    CALORIMETER geometry upgrade
374    
375       The volumes CAPD and CAAD have been taken off from the calorimeter.
376       In addition the logical tree has been slightly changed to make the shifts of
377       the silicon planes into the calorimeter box easier, i.e. the CAPL volume,
378       which was made of the CASI, CAKP, CAGL, C10C and CAKA volumes, has
379       been split up in the volumes CANS and CAPL. Now CANS is made of the CAKP,
380       CAGL, C10C and CAKA volumes while CAPL contains the CASI volume, that has to
381       be shifted as a function of the vertical position in the calorimeter. Also the
382       dimensions of some volumes have been upgraded, including the external ones:
383       CALB and CALS. CALS is an aluminum box of dimensions: 48.4*48.4*21.278 cm^3,
384       having side-walls 1 cm thick and a bottom of 1 mm. The real box is more
385       complicated and the configuration of the bottom should be upgraded if we want
386       a reliable description of the event in the S4 scintillator.
387    
388    22 October 2002, Stockholm
389    
390    ANTICOINC. GEOMETRY UPGRADE
391    
392       The AC geometry has been updated. The top AC scintillator (CAT) now
393       consists of 1 single sheet of scintillator with a hole in the middle
394       and the correct geometry(*). The side AC scintillators (CAS) also
395       have the correct shape. The AC scintillators are placed in aluminum
396       boxes with plastic rims inside. For these rims a 'new' material, PLAS,
397       was defined. PLAS has all the characteristics of SCIN but is
398       non-sensitive. No PMTs or PMT holders have been modelled.
399       (*)-The interfaces on CAT where the PMTs should be located are
400           slightly different from the real case.
401    
402  11 February 2002, Bari  11 February 2002, Bari
403    
404  MACRO CLEAN-UP  MACRO CLEAN-UP
# Line 116  TRACK COMMAND CALLED BY GPGARIN Line 500  TRACK COMMAND CALLED BY GPGARIN
500    
501  TRD IONIZATION ENERGY LOSS GENERATED NOW BY GARFIELD  TRD IONIZATION ENERGY LOSS GENERATED NOW BY GARFIELD
502     To generate the ionization in the TRD straw tubes the HEED program     To generate the ionization in the TRD straw tubes the HEED program
503     interfaced by GARFIELD is used (GEANT does not simulate the ionization     interfaced by GARFIELD is used (GEANT does not correctly simulate
504     in thin layer and in the gas, correctly). The idea is that GEANT tracks     the ionization in thin layer and in the gas). The idea is that GEANT
505     the particle in the gas and then passes the coordinates, translated in     tracks the particle in the gas and then passes the coordinates,
506     the DRS, to GARFIELD. The GARFIELD subroutines are called by GPUTRD.     translated in the DRS, to GARFIELD. The GARFIELD subroutines are
507     The energy loss and the number of clusters in TRD are stored in the     called by GPUTRD. The energy loss and the number of clusters in TRD
508     variables EGARTRD and NGARTRD of the CWN-tplu.     are stored in the variables EGARTRD and NGARTRD of the CWN-tplu.
509    
510   1 May 2001, Bari   1 May 2001, Bari
511    
# Line 200  NEW SEQUENCES ADDED: $XPRINTPLOT,$PRINTP Line 584  NEW SEQUENCES ADDED: $XPRINTPLOT,$PRINTP
584     The definition of the ITRSO detector has been changed in the GPSED routine:     The definition of the ITRSO detector has been changed in the GPSED routine:
585     NVTRD has been forced to 2 for compatibility with GPDTRD.     NVTRD has been forced to 2 for compatibility with GPDTRD.
586    
 3 april 2001, Bari  
   
   
587  28 march 2001, Bari  28 march 2001, Bari
588    
589     ITRSO has been defined as a sensitive detector in GSTMED routine and it has     ITRSO has been defined as a sensitive detector in GSTMED routine and it has
Legend:
Removed from v.3.1  
changed lines
  Added in v.3.10
  ViewVC Help
Powered by ViewVC 1.1.23