/[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.23 by cafagna, Thu Nov 30 13:22:35 2006 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.22 2006/11/28 11:43:41 pam-ba Exp $
3  #  #
4  # $Log: v_100.txt,v $  # $Log: v_100.txt,v $
5    # Revision 3.22  2006/11/28 11:43:41  pam-ba
6    # History modified
7    #
8    # Revision 3.21  2006/11/28 10:26:15  pam-ba
9    # S3 positioning completed
10    #
11    # Revision 3.20  2006/11/16 18:45:29  pam-ba
12    # Simulated an aluminum container for S4
13    #
14    # Revision 3.19  2006/11/16 12:01:13  pam-ba
15    # S4 dimensions corrected.
16    #
17    # Revision 3.18  2006/11/10 11:39:35  pam-ba
18    # S2 and S1 z-positions corrected, He3 and plystyrene mixture added, Top Plate geometry simulated and titanium mixture added.
19    #
20    # Revision 3.17  2006/10/13 16:36:59  pam-ba
21    # Added a new material, the cadmium, for ND
22    #
23    # Revision 3.16  2006/10/12 11:11:21  pam-ba
24    # ND geometry updated.
25    #
26    # Revision 3.15  2006/10/02 11:17:30  pam-ba
27    # NDET 'SPHE' data card meaning changed. Now it eliminates the whole PAMELA container.
28    #
29    # Revision 3.14  2006/06/30 15:38:16  pam-ba
30    # 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)
31    #
32    # Revision 3.13  2006/06/05 13:56:17  pamela
33    # Gigantic resonance added for gamma enetering in the calorimeter absorber
34    #
35    # Revision 3.12  2006/05/18 10:52:32  pam-ba
36    # TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
37    #
38    # Revision 3.11  2006/05/11 23:53:15  cafagna
39    # More bugs fixed in the CALO ntple structure filling
40    #
41    # Revision 3.10  2006/04/10 11:07:43  cafagna
42    # GEN data card updated, ZDGEN added
43    #
44    # Revision 3.9  2005/12/14 03:34:40  cafagna
45    # An update of the history and inform readme files.
46    #
47    # Revision 3.8  2005/12/14 03:16:08  cafagna
48    # Neutron detector added. Geometry and GPCALOR package
49    #
50    # Revision 3.7  2005/10/18 08:24:35  cafagna
51    # History updated
52    #
53    # Revision 3.6  2005/07/25 11:53:21  cafagna
54    # Several updates. See history for details
55    #
56    # Revision 3.5  2004/04/06 10:33:46  pamela
57    # NON-REPRODUCIBILITY problem of a GPAMELA RUN fixed; bug found and fixed filling in the hit structure of the calorimeter
58    #
59    # Revision 3.4  2003/12/17 11:32:50  pamela
60    # CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation
61    #
62  # Revision 3.3  2002/12/05 17:27:59  pamela  # Revision 3.3  2002/12/05 17:27:59  pamela
63  # New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated  # New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated
64  #  #
# Line 23  Line 80 
80  #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
81  #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
82  #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
83  #-- Author :    Francesco Cafagna   28/11/95  #-- Author :    Francesco Cafagna   28/11/95
84    
85    30 November 2006, Bari F. Cafagna
86    
87    BUG FIXED
88    
89       For some reason there was a conflict using FLUKA and the GARFILD
90       interface for the TRD. The file TRD_ATT.TXT was treated as flukaerr
91       file by the compiler. It is not clear why. The code was correct,
92       the logical units and names were not inconflict, but the compiler
93       assigned to them the same internal number. It may be a problem
94       related to the optimization phase. I did close the file just after
95       the read-out in GPDAT.F and commented out the closing in UGLAST.F
96    
97    
98    
99    November 2006, Bari
100    
101    
102    S4 GEOMETRY UPDATED:
103    The x and y dimensions have been corrected because they were wrong.
104    A container of aluminum has been added having walls of thickness equal
105    to 0.1 cm.  
106    
107    TOF POSITIONS DEFINITELY UPDATED:
108    S2 and S1 have been positioned again, after the simulation of the top plate.
109    Before now, the positions were put by hand based on the positions given by
110    Sergio Ricciarini in the document 'Main geometrical parameters of the PAMELA
111    sub-detectors' released by O. Adriani, L. Bonechi, E. Mocchiutti and
112    S. Ricciarini on the 20th of December 2005.. These positions were
113    lightly approximated because the thicknesses of the mylar and of the glue were
114    not considered in that document. Also the height of S3 has been modified as
115    suggested from the Naples PAMELA collaborators.
116    
117    ND GEOMETRY COMPLETED.
118    The He3 and the polyethylene (CH2) have been added to fill
119    the volumes NDTI and NDPB in the neutron detector. Leonov gave us the
120    details of these materials.
121    
122    TOP PLATE GEOMETRY ADDED.
123    The top plate geometry has been simulated as a rectangular window filled
124    with a mixture of titanium (90%), aluminum (6%) and vanadium (4%). It has been
125    put in a N2 box having the same planar dimensions of CATA. In this box there
126    are also S2 and CATA at the right positions. Then the top plate geometry has
127    been completed with other eight small volumes of titanium positioned at the
128    corners of the N2 box. Four of them are at the same height of CATA, the
129    other four are at the same height of S2. The simulated geometry of the top
130    plate is based on a simplified version of the CAD drawings.
131    
132    October 2006, Bari
133    
134    The geometry of the neutron detector has been updated. Some dimensions of some
135    volumes have been corrected and the aluminum cover and the aluminum boxes to
136    put the cables have been added. A new material, the cadmium, has been also
137    defined.
138    
139    September 2006, Bari
140    
141    SPHE and ND data card bugs fixed: the definition of the ND data card,
142    missing in the subroutine gpgeo.F, has been added; the meaning of the SPHE
143    data card has been changed. Before the correction the data card:
144    NDET 'SPHE' was used to delete the spherical top shell to substitute it with
145    a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA.
146    
147    
148    June 2006, Bari
149    
150    The center of the scintillator planes S22Y (variable ZPAMS22Y in gpdgeo.inc)
151    and S12X (variable ZPAMS12X in gpdgeo.inc) has been positioned at the
152    nominal height as measured in PAMELA (See the document: "Main geometrical
153    parameters of the PAMELA sub-detectors" by O. Adriani, L. Bonechi,
154    E. Mocchiutti, S. Ricciarini, 20 December 2005). Follows that the positions
155    of S21Y and S12X are higher than those in the cited document due to the fact
156    that in GPAMELA the thickness of the mylar has been considered while in the
157    document it has been neglected.
158    
159    
160    May 2006, Bari & Tor Vergata
161    
162    GIGANTIC RESONANCE FOR NEUTRON DETECTOR ADDED
163    
164       Routines to simulate the gigantic resonance of gammas in Tungsten
165       have been added.  The GPGIG routine is called in GUSTEP if a gamma
166       enter the calorimeter absorber.  This is the steering routine to
167       simulate the production of neutrons from gigantic resonance.  It
168       does checks on STEP lenght. If the range is smaller than the other
169       selected for that step, it does generate the neutron and stops the
170       gamma. Please note that the neutron has now a new particle
171       number. This is to tag the gigantic resonance neutrons.
172    
173    
174    May 2006, Bari & Florence
175    
176    CAL HIT STRUCTURE BUGS FIXED
177    
178       The maximum number of hit is now different for the two hit
179       structures: CALST and CALI. Vectors inizialization and HBOOK
180       ntple booking have been updated. The GPDCAL routine has been fixed
181       so to handle the case in wich hits stored are more than the maximum
182       number of hit.
183       In this case in the ntple up to the maximum number of hits will be stored.
184    
185    April 2006, Bari
186    
187    TOF GEOMETRY AND POSITIONS UPDATED AND NEW MIXTURES ADDED
188    
189       The TOF geometry has been modified. The following boxes have been
190       added: POL1, POL2 and POLY made of polystyrene, S11M, S12M, S21M,
191       S22M, S31M and S32M made of mylar, S1A, S2A and S3 made of air and
192       S1 and S2 made of aluminum. Each scintillator paddle has been put
193       in his mylar box and the other materials: air, polystyrene, and
194       aluminum have been added at their nominal positions.  According to
195       Naples people the araldite glue has been simulated has an air
196       gap. For this work two new materials: the Mylar (MYLAR) and the
197       polystyrene (POLYSTYRENE) with a density of 35 g/l have been
198       defined as a mixture.  The positions of the three bottom
199       scintillator planes that contain respectively the S12X, S22Y and
200       S32X paddles have been regulated according on their official
201       positions in PAMELA.
202    
203    Mar 2006, Bari
204    
205    GEN DATA CARD UPDATED
206    
207       To enable generation on a surface perpendicular to the XY plane,
208       GEN gata card has been updated addingh a new parameter: ZDGEN. This is
209       the dimension, along Z axis , of the generation surface. The Z
210       position will be randomply chosen according to: Z= ZDGEN*RNDM_NUMBER +
211       ZGEN, i.e. Z= GEN(6)*RNDM_NOMBER + GEN(3).
212    
213    Nov 2005, Bari
214    
215    GUHADR AND GUPHAD UPDATED
216    
217       To use GCALOR package the hadronic routines have been updated. The
218       inizialization routine call CALSIG, while the other calls GCALOR.
219    
220    NEW GPKEY ADDED: GPCALOR
221    
222       This logical has been added to enable the GCALOR package. This flag
223       is set to true in GPDAT if the data card: HPAK, is set to
224       'GCAL'. The gpkey.inc has been update accordingly.
225    
226      
227    NEUTRON DETECTOR ADDED. NEW DIR: GPND
228    
229       The neutron detector has been added. At the moment it is just the
230       geometry. The directory structure of the repository has been
231       updated as well. Dimensions has been taken from picture and
232       literature. A full upgrade to the drawing is needed.
233    
234    GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR, GPCALORDES
235    
236       GCALOR package contins the CALOR simulation code and an interface
237       to use it in GEANT. The important feature for us is the usage of
238       the MICAP code. This is facused on the low energy neutron
239       simulation. for details see:
240       http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html
241       This package should be distributed with the GEANT library but is
242       not up to date. I did download the latest release and stored into
243       gpcalor directory of the gpamela tree.
244       Then I did clean up the code substituting the explicit inclusion of
245       the commons with a #include cpp directive. In parallel I did
246       extract the commons to include files having the same common name. I
247       did store the include files into a newly created directory:
248       gpcalordes.
249       The Makefile has been updated accordingly.
250       Please note that to avoid conflict with CRENLIB distribution the gcalor source file has been named gpcalor.F
251       NOTE: There are still problem due to different common sizes. In
252       particular the common MICFIL is maller in the geant library
253       libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are
254       present and linked using a wrong version of the common. This still needs to be debuged.
255       NOTE2: The auxiliary files with the cross sections: chetc.dat.gz
256       and xsneut.dat.gz, have been added to the aux directory and moved
257       to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will
258       look for CERN_ROOT environment variable. If found files are
259       searched there at first, then in the working directory. A fool
260       proof policy has to be implemented to avoid problem with
261       synchronization fo these files.
262      
263    
264    The GCALOR package
265    
266    June 2005, Bari
267    
268    TOF SCINTILLATOR PADDLES UPDATED
269    
270       The dimensions and the number of the scintillator paddles for each
271       TOF planes have been updated.
272    
273    May 2005, Bari
274    
275    Some updates on the latest modification done in the past year.
276    
277    NEW DATA CARD ADDED: HFSF
278    
279       To define a policy for the random number initial seeds
280       definition. Using this card is possible to override GEANT seeds
281       defined via NRDM card. The policy is selected according to the
282       values:
283    
284       - 1: The seeds are initialized to the initial values found in a user
285            defined file or the default file: INPUTSEED.DAT
286      
287       - 2: The seeds are initialized to the final values found in a user defined
288            file or the default file: INPUTSEED.DAT
289    
290       The case 1 must be used in case the user needs to reproduce the
291       random chain of a previous run. In this case the user can save the
292       initial seeds, used in the run he would like to reproduce, in a
293       binary file and pass the filename to the program using the *FLSF
294       data card. In case the user file is not specified the default
295       INPUTSEED.DAT will be used.
296      
297       The case 2 must be used in case the user needs to chain several
298       GPAMELA run and likes to be sure he is starting the random
299       generator using the right sequence. In this case the user must
300       specify an input binary file using the *FLSF data card, otherwise
301       the INPUTSEED.DAT file will be used.
302    
303    NEW DATA CARD ADDED: *FSFI
304    
305       Using this card the user can specify the logical unit and name of
306       the file storing the initial seeds to be used to initialize the
307       random number generator. This file must be a FORTRAN binary one
308       storing four integer numbers. The first two are the number to be
309       used in the case: HFSF=1, the other two will be used in the case:
310       HFSF=2. This file can be one created by GPAMELA or by the user
311       filled with his own seeds. For this purpose an utility program:
312       writeseeds.f, has been added in the aux directory.  In case the
313       *FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will
314       be used as LUN and file name respectively.
315      
316    NEW DATA CARD ADDED: *LSFI
317    
318       Using this card the user can specify the logical unit and name of
319       the file storing the first and last seeds used in the GPAMELA
320       run. This file is a FORTRAN binary one. This file can be used as
321       input one specifying it in the *FSFI data card of the next GPAMELA
322       run.  In case the *LSFI card is not specified the default values: 26
323       and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN
324       and file name respectively.
325      
326    NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f
327    
328       These new programs have been added in the aux directory. Using these a
329       user defined seed file can be created and re-read.
330    
331    NEW VOLUMES ADDED: MSHE, BSPH; PRESSURIZED CONTAINER ADDED
332    
333       Alexey Bakaldin, in MEPHI, did add the PAMELA pressurized container to
334       the simulation. He did defined new volumes filled with aluminum and
335       placed inside the mother volume. Positions have been fine tuned by
336       Marialuigia Ambriola and compared to the CAD drawings.
337       Two new volumes have been added to simulate the container:
338       - MSHE, a tube simulating the middle part of the container
339       - BSPH, the spherical bottom part of the container
340    
341       To better simulate the upper part the SHEL volume has been modified
342       into a cone. Dimentions of the top cover: TSPH, have been modified
343       accordingly.
344    
345    DETECTOR POSITIONS REVIEWED
346    
347       All detector Z positions have been reviewd to fit into the
348       simulated pressurized container.
349    
350    TRD GEOMETRY AND CALIBRATION REVIEWD
351    
352       The TRD geometry has been deeply reviewed. Using the CAD drawings
353       the carbon fiber frames have been simulated and radiator dimentions
354       corrected. For this reason the calibration done on the beam tests
355       has been revied and new sets of calibration constants calculated
356       comparing the beam test data with the GPAMELA results. The new
357       constants are about 3% larger than the previous ones.
358      
359    TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI
360      
361       Thanks to Lorenzo Bonechi for the drawings and explanations. Now the
362       hybrd cards have been put into the simulation and the geometry updated
363       considering the dead zones in the silicon detectors. The hybrid zone
364       has been simulated as well. At the moment the hybrid is simulated as
365       a G10 plates. The full height of the tracker magnet has been
366       reviewed as well.
367    
368       The tracker ladder is now simulated inside a nitrogen box: TPAS,
369       placed inside an aluminum frame: TRPB. Each silicon ladder has been
370       simulated using two silicon blocks: TRSL, into each of this block a
371       smaller silicon detector: TPAI, has been placed inside the larger
372       silicon block TRSL. In this way the subdivided silicon ladder can
373       be upgraded with an indipendend roto-translation for each sensor.
374      
375       The TRPB aluminum frame has been enlarged to fit the external
376       magnet canister frame.
377      
378       The last plane has been flipped with a 180 degree rotation around
379       the X axis.
380      
381    TRACKER HIT STRUCTURE REVIEWED
382    
383       Taking into account the new version of the tracker geometry, the hit
384       structure for this detector has been revied.
385    
386    CALORIMETER GEOMETRY REVIEWED
387    
388       Marco Albi reviewed the calorimeter dimentions and positioning.
389    
390    
391    29 March 2004, Bari
392    
393    NON-REPRODUCIBILITY PROBLEM OF A GPAMELA RUN FIXED.
394       The non-reproducibility of a GPAMELA run was due to the random number
395       initialization in the GARFIELD code. In GARFIELD by default, the initial
396       seeds of the random number generators are always the same while the random
397       number generators are called a given number of times (determined by the
398       hour of the day) during the initialization phase (see init.f subroutine in
399       the GARFIELD code for details). Follows that different runs produce
400       different results without changing the initial seeds. To have identical
401       results in different runs, the GARFIELD program has to start typing the
402       noRNDM_initialisation switch. To avoid of specifying this switch
403       by the user,
404       the GARFIELD package has been upgraded with a patch. In this way the problem
405       is partially solved because, now, the initial seeds of the random generators
406       in GARFIELD will be always the same even if the RNDM GEANT data card is
407       activated by the user for changing the initial seeds in the GPAMELA program.
408       Work is in progress for a more general correction of this problem.
409       Please, use the updated GARFIELD code released with the CVS version v4r1
410       to fix this problem.  
411    
412    
413    RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS.
414       The obsolete RNDM random number generator has been replaced by the GEANT
415       GRNDN routine in the gpxtr.F subroutine and in the npoiss.F function.
416    
417    BUG FOUND AND FIXED: the set and detector calorimeter addresses (ISCAL
418       and IDCASI variables) used in GUTREV were respectively set to a fixed
419       values of 12 and 1. The correct values of these variables are stored in
420       the GPSED common when the set and the detector ZEBRA banks are filled
421       during a run. In general the values of the set and detector addresses
422       depend on the number of active detectors in a given run. ISCAL=12 and
423       IDCASI=1 are only right when all the detectors of GPAMELA are active.
424    
425  9 December 2003, Bari  9 December 2003, Bari
426    
427  CALORIMETER SIMULATION completed! The update of the geometry and of the     CALORIMETER SIMULATION completed! The update of the geometry and of the
428  special tracking parameters and the tuning of the calorimeter have been     special tracking parameters and the tuning of the calorimeter have been
429  successfully done. A great quantity of simulated data have been produced     successfully done. A great quantity of simulated data have been produced
430  in the calorimeter for different particles (muons, electrons and pions)     in the calorimeter for different particles (muons, electrons and pions)
431  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
432  distributions of the total energy deposited in the calorimeter and the     distributions of the total energy deposited in the calorimeter and the
433  total number of strips hit have been compared with the respective     total number of strips hit have been compared with the respective
434  distributions produced by the Trieste's tuned standalone Monte Carlo     distributions produced by the Trieste's tuned standalone Monte Carlo
435  simulation program of the PAMELA calorimeter. The accord between the     simulation program of the PAMELA calorimeter. The accord between the
436  two simulations is excellent. Many thanks to Mirko for his collaboration.     two simulations is excellent. Many thanks to Mirko for his collaboration.
437    
438  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
439  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
440  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.
441    
442  18 April 2003, Bari  18 April 2003, Bari
443    
444  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
445  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
446  routine.     routine.
447  Also the length of the common /PAWC/ (parameter NWPAW) has been increased     Also the length of the common /PAWC/ (parameter NWPAW) has been increased
448  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
449  of columns times the buffer size.     of columns times the buffer size.
450    
451  10 April 2003, Bari  10 April 2003, Bari
452    
453  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
454  filled have been changed according to the electronics system of the real     filled have been changed according to the electronics system of the real
455  detector. In fact, because each silicon detector (module) consists of     detector. In fact, because each silicon detector (module) consists of
456  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
457  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
458  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
459  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
460  end of the event (gutrev.F subroutine).     end of the event (gutrev.F subroutine).
461  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
462  ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL),     ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL),
463  by a call to the GPDCAL subroutine:     by a call to the GPDCAL subroutine:
464  -ICAPLANE(i) contains the number of hit plane;     -ICAPLANE(i) contains the number of hit plane;
465  -ICASTRIP(i) contains the number of hit strip;     -ICASTRIP(i) contains the number of hit strip;
466  -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
467             positions in which a `long strip' has been hit.                positions in which a `long strip' has been hit.
468  -ENESTRIP(i) contains the deposited energy in the hit strip;     -ENESTRIP(i) contains the deposited energy in the hit strip;
469  where i is the number of hit (1<i<4224).     where i is the number of hit (1<i<4224).
470  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
471  that there is a hit for each `long strip' hit from     that there is a hit for each `long strip' hit from
472  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
473  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
474  calorimeter due to the showers. Follows that NTHCAL, which is the     calorimeter due to the showers. Follows that NTHCAL, which is the
475  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
476  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
477  be less or equal to 4224.     be less or equal to 4224.
478  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
479  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
480  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.
481  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
482  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.
483  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:
484  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.
485  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
486  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
487  `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
488  other in the third one.     other in the third one.
489    
490  The energy deposited in the calorimeter is calculated in GeV.     The energy deposited in the calorimeter is calculated in GeV.
491    
492  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
493  GSCHIT.     GSCHIT.
494    
495  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
496  of a call to the GFHITS subroutine.     of a call to the GFHITS subroutine.
497    
498  25 February 2003, Bari  25 February 2003, Bari
499    
500  BUG found:  BUG found:
501  DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP,     DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP,
502  DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2,     DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2,
503  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,
504  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,
505  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,
506  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,
507  gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc,     gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc,
508  gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc.     gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc.
509  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
510  GPSTM subroutine.     GPSTM subroutine.
511    
512  Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C     Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C
513  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
514  suggested by Mirko Boezio.     suggested by Mirko Boezio.
515    
516  Updated the value of the absorption length for silicon in the calorimeter     Updated the value of the absorption length for silicon in the calorimeter
517  and tracker although this parameter is ignored by GEANT. For this reason     and tracker although this parameter is ignored by GEANT. For this reason
518  it was equal to the radiation length.     it was equal to the radiation length.
519    
520  Updated the relative positions of the calorimeter planes. The corrected     Updated the relative positions of the calorimeter planes. The corrected
521  shifting are:     shifting are:
522    
523  first view: (Dxo,Dyo)=(0.10,0.05) cm     first view: (Dxo,Dyo)=(0.10,0.05) cm
524  second view: (Dxo,Dyo)=(-0.05,0.10) cm     second view: (Dxo,Dyo)=(-0.05,0.10) cm
525  third view: (Dxo,Dyo)=(-0.10,-0.05) cm     third view: (Dxo,Dyo)=(-0.10,-0.05) cm
526  fourth view: (Dxo,Dyo)=(0.05,-0.10) cm     fourth view: (Dxo,Dyo)=(0.05,-0.10) cm
527    
528  4 November 2002, Bari  4 November 2002, Bari
529    
530  CAS detectors distances modified  CAS detectors distances modified
531    
532  The distances between the CAS detectors have been modified based on the     The distances between the CAS detectors have been modified based on the
533  latest CAD drawings.     latest CAD drawings.
534    
535  2 November 2002, Bari  2 November 2002, Bari
536    
537  CALORIMETER geometry upgrade  CALORIMETER geometry upgrade
538    
539  The volumes CAPD and CAAD have been taken off from the calorimeter.     The volumes CAPD and CAAD have been taken off from the calorimeter.
540  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
541  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,
542  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
543  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,
544  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
545  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
546  dimensions of some volumes have been upgraded, including the external ones:     dimensions of some volumes have been upgraded, including the external ones:
547  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,
548  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
549  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
550  a reliable description of the event in the S4 scintillator.     a reliable description of the event in the S4 scintillator.
551    
552  22 October 2002, Stockholm  22 October 2002, Stockholm
553    
Legend:
Removed from v.3.4  
changed lines
  Added in v.3.23
  ViewVC Help
Powered by ViewVC 1.1.23