/[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.4 by pamela, Wed Dec 17 11:32:50 2003 UTC revision 3.6 by cafagna, Mon Jul 25 11:53:21 2005 UTC
# Line 1  Line 1 
1  #  #
2  # $Id: v_100.txt,v 3.3 2002/12/05 17:27:59 pamela Exp $  # $Id: v_100.txt,v 3.5 2004/04/06 10:33:46 pamela Exp $
3  #  #
4  # $Log: v_100.txt,v $  # $Log: v_100.txt,v $
5    # Revision 3.5  2004/04/06 10:33:46  pamela
6    # NON-REPRODUCIBILITY problem of a GPAMELA RUN fixed; bug found and fixed filling in the hit structure of the calorimeter
7    #
8    # Revision 3.4  2003/12/17 11:32:50  pamela
9    # CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation
10    #
11  # Revision 3.3  2002/12/05 17:27:59  pamela  # Revision 3.3  2002/12/05 17:27:59  pamela
12  # New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated  # New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated
13  #  #
# Line 23  Line 29 
29  #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
30  #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
31  #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
32  #-- Author :    Francesco Cafagna   28/11/95  #-- Author :    Francesco Cafagna   28/11/95
33    
34    May 2005, Bari
35    
36    Some updates on the latest modification done in the past year.
37    
38    NEW DATA CARD ADDED: HFSF
39    
40       To define a policy for the random number initial seeds
41       definition. Using this card is possible to override GEANT seeds
42       defined via NRDM card. The policy is selected according to the
43       values:
44    
45       - 1: The seeds are initialized to the initial values found in a user
46            defined file or the default file: INPUTSEED.DAT
47      
48       - 2: The seeds are initialized to the final values found in a user defined
49            file or the default file: INPUTSEED.DAT
50    
51       The case 1 must be used in case the user needs to reproduce the
52       random chain of a previous run. In this case the user can save the
53       initial seeds, used in the run he would like to reproduce, in a
54       binary file and pass the filename to the program using the *FLSF
55       data card. In case the user file is not specified the default
56       INPUTSEED.DAT will be used.
57      
58       The case 2 must be used in case the user needs to chain several
59       GPAMELA run and likes to be sure he is starting the random
60       generator using the right sequence. In this case the user must
61       specify an input binary file using the *FLSF data card, otherwise
62       the INPUTSEED.DAT file will be used.
63    
64    NEW DATA CARD ADDED: *FSFI
65    
66       Using this card the user can specify the logical unit and name of
67       the file storing the initial seeds to be used to initialize the
68       random number generator. This file must be a FORTRAN binary one
69       storing four integer numbers. The first two are the number to be
70       used in the case: HFSF=1, the other two will be used in the case:
71       HFSF=2. This file can be one created by GPAMELA or by the user
72       filled with his own seeds. For this purpose an utility program:
73       writeseeds.f, has been added in the aux directory.  In case the
74       *FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will
75       be used as LUN and file name respectively.
76      
77    NEW DATA CARD ADDED: *LSFI
78    
79       Using this card the user can specify the logical unit and name of
80       the file storing the first and last seeds used in the GPAMELA
81       run. This file is a FORTRAN binary one. This file can be used as
82       input one specifying it in the *FSFI data card of the next GPAMELA
83       run.  In case the *LSFI card is not specified the default values: 26
84       and OUTPUTSEEDS.DAT, will be used as LUN and file name respectively.
85      
86    NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f
87    
88       These new programs have been added in the aux directory. Using these a
89       user defined seed file can be created and re-read.
90    
91    NEW VOLUMES ADDED: MSHE, BSPH; PRESSURIZED CONTAINER ADDED
92    
93       Alexey Bakaldin, in MEPHI, did add the PAMELA pressurized container to
94       the simulation. He did defined new volumes filled with aluminum and
95       placed inside the mother volume. Positions have been fine tuned by
96       Marialuigia Ambriola and compared to the CAD drawings.
97       Two new volumes have been added to simulate the container:
98       - MSHE, a tube simulating the middle part of the container
99       - BSPH, the spherical bottom part of the container
100    
101       To better simulate the upper part the SHEL volume has been modified
102       into a cone. Dimentions of the top cover: TSPH, have been modified
103       accordingly.
104    
105    DETECTOR POSITIONS REVIEWED
106    
107       All detector Z positions have been reviewd to fit into the
108       simulated pressurized container.
109    
110    TRD GEOMETRY AND CALIBRATION REVIEWD
111    
112       The TRD geometry has been deeply reviewed. Using the CAD drawings
113       the carbon fiber frames have been simulated and radiator dimentions
114       corrected. For this reason the calibration done on the beam tests
115       has been revied and new sets of calibration constants calculated
116       comparing the beam test data with the GPAMELA results. The new
117       constants are about 3% larger than the previous ones.
118      
119    TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI
120      
121       Thanks to Lorenzo Bonechi for the drawings and explanations. Now the
122       hybrd cards have been put into the simulation and the geometry updated
123       considering the dead zones in the silicon detectors. The hybrid zone
124       has been simulated as well. At the moment the hybrid is simulated as
125       a G10 plates. The full height of the tracker magnet has been
126       reviewed as well.
127    
128       The tracker ladder is now simulated inside a nitrogen box: TPAS,
129       placed inside an aluminum frame: TRPB. Each silicon ladder has been
130       simulated using two silicon blocks: TRSL, into each of this block a
131       smaller silicon detector: TPAI, has been placed inside the larger
132       silicon block TRSL. In this way the subdivided silicon ladder can
133       be upgraded with an indipendend roto-translation for each sensor.
134      
135       The TRPB aluminum frame has been enlarged to fit the external
136       magnet canister frame.
137      
138       The last plane has been flipped with a 180 degree rotation around
139       the X axis.
140      
141    TRACKER HIT STRUCTURE REVIEWED
142    
143       Taking into account the new version of the tracker geometry, the hit
144       structure for this detector has been revied.
145    
146    CALORIMETER GEOMETRY REVIEWED
147    
148       Marco Albi reviewd the calorimeter dimentions and positioning.
149    
150    
151    29 March 2004, Bari
152    
153    NON-REPRODUCIBILITY PROBLEM OF A GPAMELA RUN FIXED.
154       The non-reproducibility of a GPAMELA run was due to the random number
155       initialization in the GARFIELD code. In GARFIELD by default, the initial
156       seeds of the random number generators are always the same while the random
157       number generators are called a given number of times (determined by the
158       hour of the day) during the initialization phase (see init.f subroutine in
159       the GARFIELD code for details). Follows that different runs produce
160       different results without changing the initial seeds. To have identical
161       results in different runs, the GARFIELD program has to start typing the
162       noRNDM_initialisation switch. To avoid of specifying this switch
163       by the user,
164       the GARFIELD package has been upgraded with a patch. In this way the problem
165       is partially solved because, now, the initial seeds of the random generators
166       in GARFIELD will be always the same even if the RNDM GEANT data card is
167       activated by the user for changing the initial seeds in the GPAMELA program.
168       Work is in progress for a more general correction of this problem.
169       Please, use the updated GARFIELD code released with the CVS version v4r1
170       to fix this problem.  
171    
172    
173    RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS.
174       The obsolete RNDM random number generator has been replaced by the GEANT
175       GRNDN routine in the gpxtr.F subroutine and in the npoiss.F function.
176    
177    BUG FOUND AND FIXED: the set and detector calorimeter addresses (ISCAL
178       and IDCASI variables) used in GUTREV were respectively set to a fixed
179       values of 12 and 1. The correct values of these variables are stored in
180       the GPSED common when the set and the detector ZEBRA banks are filled
181       during a run. In general the values of the set and detector addresses
182       depend on the number of active detectors in a given run. ISCAL=12 and
183       IDCASI=1 are only right when all the detectors of GPAMELA are active.
184    
185  9 December 2003, Bari  9 December 2003, Bari
186    
187  CALORIMETER SIMULATION completed! The update of the geometry and of the     CALORIMETER SIMULATION completed! The update of the geometry and of the
188  special tracking parameters and the tuning of the calorimeter have been     special tracking parameters and the tuning of the calorimeter have been
189  successfully done. A great quantity of simulated data have been produced     successfully done. A great quantity of simulated data have been produced
190  in the calorimeter for different particles (muons, electrons and pions)     in the calorimeter for different particles (muons, electrons and pions)
191  and momenta (5 and 40 GeV/c) and the output data have been analyzed. The     and momenta (5 and 40 GeV/c) and the output data have been analyzed. The
192  distributions of the total energy deposited in the calorimeter and the     distributions of the total energy deposited in the calorimeter and the
193  total number of strips hit have been compared with the respective     total number of strips hit have been compared with the respective
194  distributions produced by the Trieste's tuned standalone Monte Carlo     distributions produced by the Trieste's tuned standalone Monte Carlo
195  simulation program of the PAMELA calorimeter. The accord between the     simulation program of the PAMELA calorimeter. The accord between the
196  two simulations is excellent. Many thanks to Mirko for his collaboration.     two simulations is excellent. Many thanks to Mirko for his collaboration.
197    
198  Working in progress on TRD. The GARFIELD interface to the HEED program is not     Working in progress on TRD. The GARFIELD interface to the HEED program is not
199  optimized to track particle with a charge greater than one and photons. The     optimized to track particle with a charge greater than one and photons. The
200  program print a warning message to advise the user when it is the case.     program print a warning message to advise the user when it is the case.
201    
202  18 April 2003, Bari  18 April 2003, Bari
203    
204  The buffer size of each column of the GPAMELA Ntuple has been increased to     The buffer size of each column of the GPAMELA Ntuple has been increased to
205  4096 and set equal to the record length, defined by a call to the HROPEN     4096 and set equal to the record length, defined by a call to the HROPEN
206  routine.     routine.
207  Also the length of the common /PAWC/ (parameter NWPAW) has been increased     Also the length of the common /PAWC/ (parameter NWPAW) has been increased
208  to 1.34E8, according to the rule that it has to be larger than the number     to 1.34E8, according to the rule that it has to be larger than the number
209  of columns times the buffer size.     of columns times the buffer size.
210    
211  10 April 2003, Bari  10 April 2003, Bari
212    
213  The variables in the HIT STRUCTURE of the CALORIMETER and their way to be     The variables in the HIT STRUCTURE of the CALORIMETER and their way to be
214  filled have been changed according to the electronics system of the real     filled have been changed according to the electronics system of the real
215  detector. In fact, because each silicon detector (module) consists of     detector. In fact, because each silicon detector (module) consists of
216  32 strips and each strip is connected to those belonging to the two detectors     32 strips and each strip is connected to those belonging to the two detectors
217  of the same row (or column) for forming 24 cm long strips, the sum of the     of the same row (or column) for forming 24 cm long strips, the sum of the
218  deposited energies in the strips forming a `long strip' is now calculated for     deposited energies in the strips forming a `long strip' is now calculated for
219  each event (gpucal.F subroutine) and it is stored in a hit only at the     each event (gpucal.F subroutine) and it is stored in a hit only at the
220  end of the event (gutrev.F subroutine).     end of the event (gutrev.F subroutine).
221  The output variables of the GPAMELA en-tuple are then filled in the vectors     The output variables of the GPAMELA en-tuple are then filled in the vectors
222  ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL),     ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL),
223  by a call to the GPDCAL subroutine:     by a call to the GPDCAL subroutine:
224  -ICAPLANE(i) contains the number of hit plane;     -ICAPLANE(i) contains the number of hit plane;
225  -ICASTRIP(i) contains the number of hit strip;     -ICASTRIP(i) contains the number of hit strip;
226  -ICAMOD(i) can assume different values based on the number of times and     -ICAMOD(i) can assume different values based on the number of times and
227             positions in which a `long strip' has been hit.                positions in which a `long strip' has been hit.
228  -ENESTRIP(i) contains the deposited energy in the hit strip;     -ENESTRIP(i) contains the deposited energy in the hit strip;
229  where i is the number of hit (1<i<4224).     where i is the number of hit (1<i<4224).
230  Note that in the calorimeter each hit is filled at the end of the event and     Note that in the calorimeter each hit is filled at the end of the event and
231  that there is a hit for each `long strip' hit from     that there is a hit for each `long strip' hit from
232  the particle. This use of the hit structure is different for the other     the particle. This use of the hit structure is different for the other
233  detectors and it has been considered to avoid a too big number of hit in the     detectors and it has been considered to avoid a too big number of hit in the
234  calorimeter due to the showers. Follows that NTHCAL, which is the     calorimeter due to the showers. Follows that NTHCAL, which is the
235  max number of hit in the calorimeter, is equal to 4224, the total     max number of hit in the calorimeter, is equal to 4224, the total
236  number of `long strips'. So, for each event, the real number of hit will     number of `long strips'. So, for each event, the real number of hit will
237  be less or equal to 4224.     be less or equal to 4224.
238  ICAMOD(i) is an additional information that does not exist in the real     ICAMOD(i) is an additional information that does not exist in the real
239  detector: if the strip i (i=1,32) of the module 1 or 2 or 3     detector: if the strip i (i=1,32) of the module 1 or 2 or 3
240  is hit, the value of ICAMOD(i) is respectively incremented of 1, 100, 10000.     is hit, the value of ICAMOD(i) is respectively incremented of 1, 100, 10000.
241  Analogously it is done, if it is the strip j (j=33,64) of the modules 4, 5     Analogously it is done, if it is the strip j (j=33,64) of the modules 4, 5
242  and 6 or if it is the strip k (k=65,96) of the modules 7, 8 and 9.     and 6 or if it is the strip k (k=65,96) of the modules 7, 8 and 9.
243  For example if we consider the hit 1 of an event, we could read:     For example if we consider the hit 1 of an event, we could read:
244  ICASTRIP(1)=30, ICAPLANE(1)=21, ENESTRIP(1)=0.5E-03 and ICAMOD(1)=10001.     ICASTRIP(1)=30, ICAPLANE(1)=21, ENESTRIP(1)=0.5E-03 and ICAMOD(1)=10001.
245  It means that the hit 1 contains the information that in the strip 30 of the     It means that the hit 1 contains the information that in the strip 30 of the
246  plane 21 has been deposited a total energy of 0.5E-03 GeV. In addition the     plane 21 has been deposited a total energy of 0.5E-03 GeV. In addition the
247  `long strip 30' has been hit two times, one in the first module and the     `long strip 30' has been hit two times, one in the first module and the
248  other in the third one.     other in the third one.
249    
250  The energy deposited in the calorimeter is calculated in GeV.     The energy deposited in the calorimeter is calculated in GeV.
251    
252  To store the hits in the calorimeter the subroutine GSAHIT is used instead of     To store the hits in the calorimeter the subroutine GSAHIT is used instead of
253  GSCHIT.     GSCHIT.
254    
255  To retrieve the hit structure the call to the routine GPRHIT is done instead     To retrieve the hit structure the call to the routine GPRHIT is done instead
256  of a call to the GFHITS subroutine.     of a call to the GFHITS subroutine.
257    
258  25 February 2003, Bari  25 February 2003, Bari
259    
260  BUG found:  BUG found:
261  DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP,     DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP,
262  DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2,     DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2,
263  DCUTEW, DCUTEXE variables missed in the commons: gpaer.inc, gpal.inc, gpce.inc,     DCUTEW, DCUTEXE variables missed in the commons: gpaer.inc, gpal.inc, gpce.inc,
264  gpcp.inc, gpfe.inc, gpg10c.inc, gpg10.inc, gpkap.inc, gpn2g.inc, gproa.inc,     gpcp.inc, gpfe.inc, gpg10c.inc, gpg10.inc, gpkap.inc, gpn2g.inc, gproa.inc,
265  gpscin.inc (obsolete), gpscint.inc, gpsica.inc, gptrad.inc, gpw2.inc, gpw.inc,     gpscin.inc (obsolete), gpscint.inc, gpsica.inc, gptrad.inc, gpw2.inc, gpw.inc,
266  gpxe.inc, gpdaer.inc, gpdal.inc, gpdce.inc, gpdcp.inc, gpdfe.inc, gpdg10c.inc,     gpxe.inc, gpdaer.inc, gpdal.inc, gpdce.inc, gpdcp.inc, gpdfe.inc, gpdg10c.inc,
267  gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc,     gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc,
268  gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc.     gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc.
269  They have been added in these commons and they have been initialized in the     They have been added in these commons and they have been initialized in the
270  GPSTM subroutine.     GPSTM subroutine.
271    
272  Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C     Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C
273  in the subroutines gpsica.F, gpw2.F, gpw.F, gpce.F and gpg10c.F. They were     in the subroutines gpsica.F, gpw2.F, gpw.F, gpce.F and gpg10c.F. They were
274  suggested by Mirko Boezio.     suggested by Mirko Boezio.
275    
276  Updated the value of the absorption length for silicon in the calorimeter     Updated the value of the absorption length for silicon in the calorimeter
277  and tracker although this parameter is ignored by GEANT. For this reason     and tracker although this parameter is ignored by GEANT. For this reason
278  it was equal to the radiation length.     it was equal to the radiation length.
279    
280  Updated the relative positions of the calorimeter planes. The corrected     Updated the relative positions of the calorimeter planes. The corrected
281  shifting are:     shifting are:
282    
283  first view: (Dxo,Dyo)=(0.10,0.05) cm     first view: (Dxo,Dyo)=(0.10,0.05) cm
284  second view: (Dxo,Dyo)=(-0.05,0.10) cm     second view: (Dxo,Dyo)=(-0.05,0.10) cm
285  third view: (Dxo,Dyo)=(-0.10,-0.05) cm     third view: (Dxo,Dyo)=(-0.10,-0.05) cm
286  fourth view: (Dxo,Dyo)=(0.05,-0.10) cm     fourth view: (Dxo,Dyo)=(0.05,-0.10) cm
287    
288  4 November 2002, Bari  4 November 2002, Bari
289    
290  CAS detectors distances modified  CAS detectors distances modified
291    
292  The distances between the CAS detectors have been modified based on the     The distances between the CAS detectors have been modified based on the
293  latest CAD drawings.     latest CAD drawings.
294    
295  2 November 2002, Bari  2 November 2002, Bari
296    
297  CALORIMETER geometry upgrade  CALORIMETER geometry upgrade
298    
299  The volumes CAPD and CAAD have been taken off from the calorimeter.     The volumes CAPD and CAAD have been taken off from the calorimeter.
300  In addition the logical tree has been slightly changed to make the shifts of     In addition the logical tree has been slightly changed to make the shifts of
301  the silicon planes into the calorimeter box easier, i.e. the CAPL volume,     the silicon planes into the calorimeter box easier, i.e. the CAPL volume,
302  which was made of the CASI, CAKP, CAGL, C10C and CAKA volumes, has     which was made of the CASI, CAKP, CAGL, C10C and CAKA volumes, has
303  been split up in the volumes CANS and CAPL. Now CANS is made of the CAKP,     been split up in the volumes CANS and CAPL. Now CANS is made of the CAKP,
304  CAGL, C10C and CAKA volumes while CAPL contains the CASI volume, that has to     CAGL, C10C and CAKA volumes while CAPL contains the CASI volume, that has to
305  be shifted as a function of the vertical position in the calorimeter. Also the     be shifted as a function of the vertical position in the calorimeter. Also the
306  dimensions of some volumes have been upgraded, including the external ones:     dimensions of some volumes have been upgraded, including the external ones:
307  CALB and CALS. CALS is an aluminum box of dimensions: 48.4*48.4*21.278 cm^3,     CALB and CALS. CALS is an aluminum box of dimensions: 48.4*48.4*21.278 cm^3,
308  having side-walls 1 cm thick and a bottom of 1 mm. The real box is more     having side-walls 1 cm thick and a bottom of 1 mm. The real box is more
309  complicated and the configuration of the bottom should be upgraded if we want     complicated and the configuration of the bottom should be upgraded if we want
310  a reliable description of the event in the S4 scintillator.     a reliable description of the event in the S4 scintillator.
311    
312  22 October 2002, Stockholm  22 October 2002, Stockholm
313    
Legend:
Removed from v.3.4  
changed lines
  Added in v.3.6
  ViewVC Help
Powered by ViewVC 1.1.23