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