/[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.1.1 by cafagna, Thu Jul 11 16:01:59 2002 UTC revision 3.28 by pamela, Tue Jan 29 18:25:18 2008 UTC
# Line 1  Line 1 
1    # Revision 3.21  2006/12/05 18:45:29  bruno
2    # Update of the CARD geometry
3    # Simulated also legs, table and magnetic screen
4  #  #
5  # $Id$  # $Id: v_100.txt,v 3.27 2007/10/02 17:58:32 cafagna Exp $
6    #
7    # $Log: v_100.txt,v $
8    # Revision 3.27  2007/10/02 17:58:32  cafagna
9    # Power law spectra introduced
10    #
11    # Revision 3.26  2007/03/30 15:24:54  cafagna
12    # Calo geometry modified. X-view plane shifts in X coordinates have been swapped, i.e. they do start with X=0.05 instead of X=-0.05
13    #
14    # Revision 3.25  2006/12/17 14:53:18  cafagna
15    # Full review of the CARD geometry. S1 table added
16    #
17    # Revision 3.20  2006/11/16 18:45:29  pam-ba
18    # Simulated an aluminum container for S4
19    #
20    # Revision 3.19  2006/11/16 12:01:13  pam-ba
21    # S4 dimensions corrected.
22    #
23    # Revision 3.18  2006/11/10 11:39:35  pam-ba
24    # S2 and S1 z-positions corrected, He3 and plystyrene mixture added, Top Plate geometry simulated and titanium mixture added.
25    #
26    # Revision 3.17  2006/10/13 16:36:59  pam-ba
27    # Added a new material, the cadmium, for ND
28    #
29    # Revision 3.16  2006/10/12 11:11:21  pam-ba
30    # ND geometry updated.
31    #
32    # Revision 3.15  2006/10/02 11:17:30  pam-ba
33    # NDET 'SPHE' data card meaning changed. Now it eliminates the whole PAMELA container.
34    #
35    # Revision 3.14  2006/06/30 15:38:16  pam-ba
36    # 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)
37    #
38    # Revision 3.13  2006/06/05 13:56:17  pamela
39    # Gigantic resonance added for gamma enetering in the calorimeter absorber
40    #
41    # Revision 3.12  2006/05/18 10:52:32  pam-ba
42    # TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
43    #
44    # Revision 3.11  2006/05/11 23:53:15  cafagna
45    # More bugs fixed in the CALO ntple structure filling
46    #
47    # Revision 3.10  2006/04/10 11:07:43  cafagna
48    # GEN data card updated, ZDGEN added
49    #
50    # Revision 3.9  2005/12/14 03:34:40  cafagna
51    # An update of the history and inform readme files.
52    #
53    # Revision 3.8  2005/12/14 03:16:08  cafagna
54    # Neutron detector added. Geometry and GPCALOR package
55    #
56    # Revision 3.7  2005/10/18 08:24:35  cafagna
57    # History updated
58    #
59    # Revision 3.6  2005/07/25 11:53:21  cafagna
60    # Several updates. See history for details
61    #
62    # Revision 3.5  2004/04/06 10:33:46  pamela
63    # NON-REPRODUCIBILITY problem of a GPAMELA RUN fixed; bug found and fixed filling in the hit structure of the calorimeter
64    #
65    # Revision 3.4  2003/12/17 11:32:50  pamela
66    # CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation
67    #
68    # Revision 3.3  2002/12/05 17:27:59  pamela
69    # New GARFIELD.GAR file added and GPAMELA.FFR cleaned and updated
70    #
71    # Revision 3.2  2002/12/05 10:17:42  pamela
72    # Update CAS and CALO geometries and positions. Makefile updated as well
73    #
74    # Revision 3.1.1.1  2002/07/11 16:01:59  cafagna
75    # First GPAMELA release on CVS
76  #  #
 # $Log$  
77  #  #
78  #CMZ :  3.00/00 11/02/2002  20.05.23  by  Unknown  #CMZ :  3.00/00 11/02/2002  20.05.23  by  Unknown
79  #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 15  Line 87 
87  #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
88  #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
89  #-- Author :    Francesco Cafagna   28/11/95  #-- Author :    Francesco Cafagna   28/11/95
90    
91    Jan 2008, Bari
92    
93    PRESSURIZED CONTAINER GEOMETRY:
94       Pressurized container geometry has been updated: top container TSPH
95       now has been placed at the maximum distance, quoted in the ICD,
96       from S1: 5cm. This distance corresponds to the maximum one in
97       flight due to the cover deformation due to the internal
98       pressure. Three thermal cover layers (nylon6, mylar and aluminium,
99       respectively in the new volumes COV1, COV2 and COV3) have been
100       added to the top container.
101    
102    NEW MATERIALS AND TRACKING MEDIA:
103    
104       Two new materials have been implemented: AMg6M mixture for volume
105       TSPH, and Nylon6 mixture for volume COV1.  New tracking medium have
106       been added for Nylon6 Mixture (TNYL), AMg6M Mixture (TAMG) and Iron
107       (TIRO).
108    
109    TOF GEOMETRY:
110       TOF geometry has been updated: Scintillators positions in the
111       GPAMELA reference system have been corrected shifting all
112       scintillator paddle by 23 microns.
113    
114    TRACK NUMBER RECORDED FOR TOF:
115       The TOF structure now register the stack number into the HIT
116       structure. In GUSTEP the stack number is tested, if is greater than
117       0, i.e. the particle is not the primary, the particle number is set
118       to 1000+IPART and passed to GPUTOF. Accordingly, now IPARTOF
119       variable can assume 2000 values instead of 100.
120    
121    POWER LAW SPECTRA GENERATION:
122       PLAW datacard parameters have been increased to 4. The meaning of
123       the 4th parameter is: extract random kinetic energy if PLAW(4)<0.,
124       or extract random momentum if PLAW(4)>=0..
125    
126    Sept 2007, Bari
127    
128    POWER LAW SPECTRA GENERATION:
129       It is now possible to select a power law spectra. Using the
130       datacard: PLAW, is possible to pass the spectral index and the
131       kinetic energy range. In PAWL(1) the user can pass the signed power
132       law index, i.e. PAWL -2.5, will generate a power law spectrum
133       according to exp(-2.5). Please note that the range is in kinetic
134       energy. In gukine the generated particle charateristics will be
135       retrived and the mass calculated from this infos. The KINE data
136       card assigned momentum will be overriden as well.
137    
138    PATHLENGHT CALCULATION MODIFIED IN TOF:
139       The pathlength calulation saved in the TOF hit structure has now
140       been modified. Insted of the pathlenght inside each paddle the
141       average pathlenght from the generation vertex to the hit in the
142       paddle, is saved. In this way it will be possible to calculate the
143       beta using this variable.
144    
145    ERROR IN MATERIAL ASSIGNMENT FOR NEUTRON DETECTOR:
146       The material assigned to the NDTU volume was wrong. It was
147       scintillator, IRON has been asigned now.
148    
149    March 2007, Bari
150    
151    CALO GEOMETRY:
152       An error was found in the original document describing PAMELA
153       geometry.  The shift along X axis of the X view planes are swapped,
154       i.e. they don't start with X=+0.05 but with X=-0.05.  The sequence
155       is: 0.05m -0.05, 0.05, -0.05 and so on.
156    
157    MAKEFILE FIXED:
158       Option -fdebug_kludge has been discharged
159    
160    GASATT WRNING FIXED:
161       There was a warning message from a graphic routine. The GSATT calls has been commented out.
162    
163    
164    December 2006, Bari
165    
166    CARD GEOMETRY:
167    Positions, dimensions and quotes have been updated. Indeed, the aluminium table
168    under S1, legs included, and the magnetic screen between the CARD and S2 have
169    been added.
170    
171    #-------------------------------------------------------------------------
172    
173    November 2006, Bari
174    
175    
176    S4 GEOMETRY UPDATED:
177    The x and y dimensions have been corrected because they were wrong.
178    A container of aluminum has been added having walls of thickness equal
179    to 0.1 cm.
180    
181    TOF POSITIONS DEFINITELY UPDATED:
182    S2 and S1 have been positioned again, after the simulation of the top plate.
183    Before now, the positions were put ad hoc based on the positions given by
184    Sergio Ricciarini in the document 'Main geometrical parameters of the PAMELA
185    sub-detectors' released by O. Adriani, L. Bonechi, E. Mocchiutti and
186    S. Ricciarini on the 20th of December 2005.. These positions were
187    lightly approximated.
188    
189    ND GEOMETRY COMPLETED.
190    The He3 and the polyethylene (CH2) have been added to fill
191    the volumes NDTI and NDPB in the neutron detector. Leonov gave us the
192    details of these materials.
193    
194    TOP PLATE GEOMETRY ADDED.
195    The top plate geometry has been simulated as a N2 box that has the same
196    dimensions of CATA. Inside it, at the right positions, there are S2 and CATA.
197    The real top plate is a rectangular window having a thickness of 0.5 cm
198    (volume TPTL emptied by the volume TPCV). This volume is filled with a mixture
199    of titanium (90%), aluminum (6%) and vanadium (4%). Other parts of titanium
200    are the four volumes TPTU, which are positioned at the corners of TPTL at the
201    same height of CATA, and the four volumes TPTM, which are positioned at the
202    corners of TPTL, at the same height of S2.
203    The geometry  of the top plate is based on a simplified version of the CAD
204    drawings.
205    
206    October 2006, Bari
207    
208    The geometry of the neutron detector has been updated. Some dimensions of some
209    volumes have been corrected and the aluminum cover and the aluminum boxes to
210    put the cables have been added. A new material, the cadmium, has been also
211    defined.
212    
213    September 2006, Bari
214    
215    SPHE and ND data card bugs fixed: the definition of the ND data card,
216    missing in the subroutine gpgeo.F, has been added; the meaning of the SPHE
217    data card has been changed. Before the correction the data card:
218    NDET 'SPHE' was used to delete the spherical top shell to substitute it with
219    a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA.
220    
221    
222    June 2006, Bari
223    
224    The center of the scintillator planes S22Y (variable ZPAMS22Y in gpdgeo.inc)
225    and S12X (variable ZPAMS12X in gpdgeo.inc) has been positioned at the
226    nominal height as measured in PAMELA (See the document: "Main geometrical
227    parameters of the PAMELA sub-detectors" by O. Adriani, L. Bonechi,
228    E. Mocchiutti, S. Ricciarini, 20 December 2005). Follows that the positions
229    of S21Y and S12X are higher than those in the cited document due to the fact
230    that in GPAMELA the thickness of the mylar has been considered while in the
231    document it has been neglected.
232    
233    
234    May 2006, Bari & Tor Vergata
235    
236    GIGANTIC RESONANCE FOR NEUTRON DETECTOR ADDED
237    
238       Routines to simulate the gigantic resonance of gammas in Tungsten
239       have been added.  The GPGIG routine is called in GUSTEP if a gamma
240       enter the calorimeter absorber.  This is the steering routine to
241       simulate the production of neutrons from gigantic resonance.  It
242       does checks on STEP lenght. If the range is smaller than the other
243       selected for that step, it does generate the neutron and stops the
244       gamma. Please note that the neutron has now a new particle
245       number. This is to tag the gigantic resonance neutrons.
246    
247    
248    May 2006, Bari & Florence
249    
250    CAL HIT STRUCTURE BUGS FIXED
251    
252       The maximum number of hit is now different for the two hit
253       structures: CALST and CALI. Vectors inizialization and HBOOK
254       ntple booking have been updated. The GPDCAL routine has been fixed
255       so to handle the case in wich hits stored are more than the maximum
256       number of hit.
257       In this case in the ntple up to the maximum number of hits will be stored.
258    
259    April 2006, Bari
260    
261    TOF GEOMETRY AND POSITIONS UPDATED AND NEW MIXTURES ADDED
262    
263       The TOF geometry has been modified. The following boxes have been
264       added: POL1, POL2 and POLY made of polystyrene, S11M, S12M, S21M,
265       S22M, S31M and S32M made of mylar, S1A, S2A and S3 made of air and
266       S1 and S2 made of aluminum. Each scintillator paddle has been put
267       in his mylar box and the other materials: air, polystyrene, and
268       aluminum have been added at their nominal positions.  According to
269       Naples people the araldite glue has been simulated has an air
270       gap. For this work two new materials: the Mylar (MYLAR) and the
271       polystyrene (POLYSTYRENE) with a density of 35 g/l have been
272       defined as a mixture.  The positions of the three bottom
273       scintillator planes that contain respectively the S12X, S22Y and
274       S32X paddles have been regulated according on their official
275       positions in PAMELA.
276    
277    Mar 2006, Bari
278    
279    GEN DATA CARD UPDATED
280    
281       To enable generation on a surface perpendicular to the XY plane,
282       GEN gata card has been updated addingh a new parameter: ZDGEN. This is
283       the dimension, along Z axis , of the generation surface. The Z
284       position will be randomply chosen according to: Z= ZDGEN*RNDM_NUMBER +
285       ZGEN, i.e. Z= GEN(6)*RNDM_NOMBER + GEN(3).
286    
287    Nov 2005, Bari
288    
289    GUHADR AND GUPHAD UPDATED
290    
291       To use GCALOR package the hadronic routines have been updated. The
292       inizialization routine call CALSIG, while the other calls GCALOR.
293    
294    NEW GPKEY ADDED: GPCALOR
295    
296       This logical has been added to enable the GCALOR package. This flag
297       is set to true in GPDAT if the data card: HPAK, is set to
298       'GCAL'. The gpkey.inc has been update accordingly.
299    
300    
301    NEUTRON DETECTOR ADDED. NEW DIR: GPND
302    
303       The neutron detector has been added. At the moment it is just the
304       geometry. The directory structure of the repository has been
305       updated as well. Dimensions has been taken from picture and
306       literature. A full upgrade to the drawing is needed.
307    
308    GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR, GPCALORDES
309    
310       GCALOR package contins the CALOR simulation code and an interface
311       to use it in GEANT. The important feature for us is the usage of
312       the MICAP code. This is facused on the low energy neutron
313       simulation. for details see:
314       http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html
315       This package should be distributed with the GEANT library but is
316       not up to date. I did download the latest release and stored into
317       gpcalor directory of the gpamela tree.
318       Then I did clean up the code substituting the explicit inclusion of
319       the commons with a #include cpp directive. In parallel I did
320       extract the commons to include files having the same common name. I
321       did store the include files into a newly created directory:
322       gpcalordes.
323       The Makefile has been updated accordingly.
324       Please note that to avoid conflict with CRENLIB distribution the gcalor source file has been named gpcalor.F
325       NOTE: There are still problem due to different common sizes. In
326       particular the common MICFIL is maller in the geant library
327       libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are
328       present and linked using a wrong version of the common. This still needs to be debuged.
329       NOTE2: The auxiliary files with the cross sections: chetc.dat.gz
330       and xsneut.dat.gz, have been added to the aux directory and moved
331       to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will
332       look for CERN_ROOT environment variable. If found files are
333       searched there at first, then in the working directory. A fool
334       proof policy has to be implemented to avoid problem with
335       synchronization fo these files.
336    
337    
338    The GCALOR package
339    
340    June 2005, Bari
341    
342    TOF SCINTILLATOR PADDLES UPDATED
343    
344       The dimensions and the number of the scintillator paddles for each
345       TOF planes have been updated.
346    
347    May 2005, Bari
348    
349    Some updates on the latest modification done in the past year.
350    
351    NEW DATA CARD ADDED: HFSF
352    
353       To define a policy for the random number initial seeds
354       definition. Using this card is possible to override GEANT seeds
355       defined via NRDM card. The policy is selected according to the
356       values:
357    
358       - 1: The seeds are initialized to the initial values found in a user
359            defined file or the default file: INPUTSEED.DAT
360    
361       - 2: The seeds are initialized to the final values found in a user defined
362            file or the default file: INPUTSEED.DAT
363    
364       The case 1 must be used in case the user needs to reproduce the
365       random chain of a previous run. In this case the user can save the
366       initial seeds, used in the run he would like to reproduce, in a
367       binary file and pass the filename to the program using the *FLSF
368       data card. In case the user file is not specified the default
369       INPUTSEED.DAT will be used.
370    
371       The case 2 must be used in case the user needs to chain several
372       GPAMELA run and likes to be sure he is starting the random
373       generator using the right sequence. In this case the user must
374       specify an input binary file using the *FLSF data card, otherwise
375       the INPUTSEED.DAT file will be used.
376    
377    NEW DATA CARD ADDED: *FSFI
378    
379       Using this card the user can specify the logical unit and name of
380       the file storing the initial seeds to be used to initialize the
381       random number generator. This file must be a FORTRAN binary one
382       storing four integer numbers. The first two are the number to be
383       used in the case: HFSF=1, the other two will be used in the case:
384       HFSF=2. This file can be one created by GPAMELA or by the user
385       filled with his own seeds. For this purpose an utility program:
386       writeseeds.f, has been added in the aux directory.  In case the
387       *FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will
388       be used as LUN and file name respectively.
389    
390    NEW DATA CARD ADDED: *LSFI
391    
392       Using this card the user can specify the logical unit and name of
393       the file storing the first and last seeds used in the GPAMELA
394       run. This file is a FORTRAN binary one. This file can be used as
395       input one specifying it in the *FSFI data card of the next GPAMELA
396       run.  In case the *LSFI card is not specified the default values: 26
397       and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN
398       and file name respectively.
399    
400    NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f
401    
402       These new programs have been added in the aux directory. Using these a
403       user defined seed file can be created and re-read.
404    
405    NEW VOLUMES ADDED: MSHE, BSPH; PRESSURIZED CONTAINER ADDED
406    
407       Alexey Bakaldin, in MEPHI, did add the PAMELA pressurized container to
408       the simulation. He did defined new volumes filled with aluminum and
409       placed inside the mother volume. Positions have been fine tuned by
410       Marialuigia Ambriola and compared to the CAD drawings.
411       Two new volumes have been added to simulate the container:
412       - MSHE, a tube simulating the middle part of the container
413       - BSPH, the spherical bottom part of the container
414    
415       To better simulate the upper part the SHEL volume has been modified
416       into a cone. Dimentions of the top cover: TSPH, have been modified
417       accordingly.
418    
419    DETECTOR POSITIONS REVIEWED
420    
421       All detector Z positions have been reviewd to fit into the
422       simulated pressurized container.
423    
424    TRD GEOMETRY AND CALIBRATION REVIEWD
425    
426       The TRD geometry has been deeply reviewed. Using the CAD drawings
427       the carbon fiber frames have been simulated and radiator dimentions
428       corrected. For this reason the calibration done on the beam tests
429       has been revied and new sets of calibration constants calculated
430       comparing the beam test data with the GPAMELA results. The new
431       constants are about 3% larger than the previous ones.
432    
433    TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI
434    
435       Thanks to Lorenzo Bonechi for the drawings and explanations. Now the
436       hybrd cards have been put into the simulation and the geometry updated
437       considering the dead zones in the silicon detectors. The hybrid zone
438       has been simulated as well. At the moment the hybrid is simulated as
439       a G10 plates. The full height of the tracker magnet has been
440       reviewed as well.
441    
442       The tracker ladder is now simulated inside a nitrogen box: TPAS,
443       placed inside an aluminum frame: TRPB. Each silicon ladder has been
444       simulated using two silicon blocks: TRSL, into each of this block a
445       smaller silicon detector: TPAI, has been placed inside the larger
446       silicon block TRSL. In this way the subdivided silicon ladder can
447       be upgraded with an indipendend roto-translation for each sensor.
448    
449       The TRPB aluminum frame has been enlarged to fit the external
450       magnet canister frame.
451    
452       The last plane has been flipped with a 180 degree rotation around
453       the X axis.
454    
455    TRACKER HIT STRUCTURE REVIEWED
456    
457       Taking into account the new version of the tracker geometry, the hit
458       structure for this detector has been revied.
459    
460    CALORIMETER GEOMETRY REVIEWED
461    
462       Marco Albi reviewed the calorimeter dimentions and positioning.
463    
464    
465    29 March 2004, Bari
466    
467    NON-REPRODUCIBILITY PROBLEM OF A GPAMELA RUN FIXED.
468       The non-reproducibility of a GPAMELA run was due to the random number
469       initialization in the GARFIELD code. In GARFIELD by default, the initial
470       seeds of the random number generators are always the same while the random
471       number generators are called a given number of times (determined by the
472       hour of the day) during the initialization phase (see init.f subroutine in
473       the GARFIELD code for details). Follows that different runs produce
474       different results without changing the initial seeds. To have identical
475       results in different runs, the GARFIELD program has to start typing the
476       noRNDM_initialisation switch. To avoid of specifying this switch
477       by the user,
478       the GARFIELD package has been upgraded with a patch. In this way the problem
479       is partially solved because, now, the initial seeds of the random generators
480       in GARFIELD will be always the same even if the RNDM GEANT data card is
481       activated by the user for changing the initial seeds in the GPAMELA program.
482       Work is in progress for a more general correction of this problem.
483       Please, use the updated GARFIELD code released with the CVS version v4r1
484       to fix this problem.
485    
486    
487    RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS.
488       The obsolete RNDM random number generator has been replaced by the GEANT
489       GRNDN routine in the gpxtr.F subroutine and in the npoiss.F function.
490    
491    BUG FOUND AND FIXED: the set and detector calorimeter addresses (ISCAL
492       and IDCASI variables) used in GUTREV were respectively set to a fixed
493       values of 12 and 1. The correct values of these variables are stored in
494       the GPSED common when the set and the detector ZEBRA banks are filled
495       during a run. In general the values of the set and detector addresses
496       depend on the number of active detectors in a given run. ISCAL=12 and
497       IDCASI=1 are only right when all the detectors of GPAMELA are active.
498    
499    9 December 2003, Bari
500    
501       CALORIMETER SIMULATION completed! The update of the geometry and of the
502       special tracking parameters and the tuning of the calorimeter have been
503       successfully done. A great quantity of simulated data have been produced
504       in the calorimeter for different particles (muons, electrons and pions)
505       and momenta (5 and 40 GeV/c) and the output data have been analyzed. The
506       distributions of the total energy deposited in the calorimeter and the
507       total number of strips hit have been compared with the respective
508       distributions produced by the Trieste's tuned standalone Monte Carlo
509       simulation program of the PAMELA calorimeter. The accord between the
510       two simulations is excellent. Many thanks to Mirko for his collaboration.
511    
512       Working in progress on TRD. The GARFIELD interface to the HEED program is not
513       optimized to track particle with a charge greater than one and photons. The
514       program print a warning message to advise the user when it is the case.
515    
516    18 April 2003, Bari
517    
518       The buffer size of each column of the GPAMELA Ntuple has been increased to
519       4096 and set equal to the record length, defined by a call to the HROPEN
520       routine.
521       Also the length of the common /PAWC/ (parameter NWPAW) has been increased
522       to 1.34E8, according to the rule that it has to be larger than the number
523       of columns times the buffer size.
524    
525    10 April 2003, Bari
526    
527       The variables in the HIT STRUCTURE of the CALORIMETER and their way to be
528       filled have been changed according to the electronics system of the real
529       detector. In fact, because each silicon detector (module) consists of
530       32 strips and each strip is connected to those belonging to the two detectors
531       of the same row (or column) for forming 24 cm long strips, the sum of the
532       deposited energies in the strips forming a `long strip' is now calculated for
533       each event (gpucal.F subroutine) and it is stored in a hit only at the
534       end of the event (gutrev.F subroutine).
535       The output variables of the GPAMELA en-tuple are then filled in the vectors
536       ICAPLANE(NTHCAL), ICASTRIP(NTHCAL), ENESTRIP(NTHCAL) and ICAMOD(NTHCAL),
537       by a call to the GPDCAL subroutine:
538       -ICAPLANE(i) contains the number of hit plane;
539       -ICASTRIP(i) contains the number of hit strip;
540       -ICAMOD(i) can assume different values based on the number of times and
541              positions in which a `long strip' has been hit.
542       -ENESTRIP(i) contains the deposited energy in the hit strip;
543       where i is the number of hit (1<i<4224).
544       Note that in the calorimeter each hit is filled at the end of the event and
545       that there is a hit for each `long strip' hit from
546       the particle. This use of the hit structure is different for the other
547       detectors and it has been considered to avoid a too big number of hit in the
548       calorimeter due to the showers. Follows that NTHCAL, which is the
549       max number of hit in the calorimeter, is equal to 4224, the total
550       number of `long strips'. So, for each event, the real number of hit will
551       be less or equal to 4224.
552       ICAMOD(i) is an additional information that does not exist in the real
553       detector: if the strip i (i=1,32) of the module 1 or 2 or 3
554       is hit, the value of ICAMOD(i) is respectively incremented of 1, 100, 10000.
555       Analogously it is done, if it is the strip j (j=33,64) of the modules 4, 5
556       and 6 or if it is the strip k (k=65,96) of the modules 7, 8 and 9.
557       For example if we consider the hit 1 of an event, we could read:
558       ICASTRIP(1)=30, ICAPLANE(1)=21, ENESTRIP(1)=0.5E-03 and ICAMOD(1)=10001.
559       It means that the hit 1 contains the information that in the strip 30 of the
560       plane 21 has been deposited a total energy of 0.5E-03 GeV. In addition the
561       `long strip 30' has been hit two times, one in the first module and the
562       other in the third one.
563    
564       The energy deposited in the calorimeter is calculated in GeV.
565    
566       To store the hits in the calorimeter the subroutine GSAHIT is used instead of
567       GSCHIT.
568    
569       To retrieve the hit structure the call to the routine GPRHIT is done instead
570       of a call to the GFHITS subroutine.
571    
572    25 February 2003, Bari
573    
574    BUG found:
575       DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP,
576       DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2,
577       DCUTEW, DCUTEXE variables missed in the commons: gpaer.inc, gpal.inc, gpce.inc,
578       gpcp.inc, gpfe.inc, gpg10c.inc, gpg10.inc, gpkap.inc, gpn2g.inc, gproa.inc,
579       gpscin.inc (obsolete), gpscint.inc, gpsica.inc, gptrad.inc, gpw2.inc, gpw.inc,
580       gpxe.inc, gpdaer.inc, gpdal.inc, gpdce.inc, gpdcp.inc, gpdfe.inc, gpdg10c.inc,
581       gpdg10.inc, gpdkap.inc, gpdn2g.inc, gpdroa.inc, gpdscin.inc, gpdsica.inc,
582       gpdtrad.inc, gpdw2.inc, gpdw.inc, gpdxe.inc.
583       They have been added in these commons and they have been initialized in the
584       GPSTM subroutine.
585    
586       Updated the special tracking parameters SICALO, TUNGA, KAOLINITE and G10C
587       in the subroutines gpsica.F, gpw2.F, gpw.F, gpce.F and gpg10c.F. They were
588       suggested by Mirko Boezio.
589    
590       Updated the value of the absorption length for silicon in the calorimeter
591       and tracker although this parameter is ignored by GEANT. For this reason
592       it was equal to the radiation length.
593    
594       Updated the relative positions of the calorimeter planes. The corrected
595       shifting are:
596    
597       first view: (Dxo,Dyo)=(0.10,0.05) cm
598       second view: (Dxo,Dyo)=(-0.05,0.10) cm
599       third view: (Dxo,Dyo)=(-0.10,-0.05) cm
600       fourth view: (Dxo,Dyo)=(0.05,-0.10) cm
601    
602    4 November 2002, Bari
603    
604    CAS detectors distances modified
605    
606       The distances between the CAS detectors have been modified based on the
607       latest CAD drawings.
608    
609    2 November 2002, Bari
610    
611    CALORIMETER geometry upgrade
612    
613       The volumes CAPD and CAAD have been taken off from the calorimeter.
614       In addition the logical tree has been slightly changed to make the shifts of
615       the silicon planes into the calorimeter box easier, i.e. the CAPL volume,
616       which was made of the CASI, CAKP, CAGL, C10C and CAKA volumes, has
617       been split up in the volumes CANS and CAPL. Now CANS is made of the CAKP,
618       CAGL, C10C and CAKA volumes while CAPL contains the CASI volume, that has to
619       be shifted as a function of the vertical position in the calorimeter. Also the
620       dimensions of some volumes have been upgraded, including the external ones:
621       CALB and CALS. CALS is an aluminum box of dimensions: 48.4*48.4*21.278 cm^3,
622       having side-walls 1 cm thick and a bottom of 1 mm. The real box is more
623       complicated and the configuration of the bottom should be upgraded if we want
624       a reliable description of the event in the S4 scintillator.
625    
626    22 October 2002, Stockholm
627    
628    ANTICOINC. GEOMETRY UPGRADE
629    
630       The AC geometry has been updated. The top AC scintillator (CAT) now
631       consists of 1 single sheet of scintillator with a hole in the middle
632       and the correct geometry(*). The side AC scintillators (CAS) also
633       have the correct shape. The AC scintillators are placed in aluminum
634       boxes with plastic rims inside. For these rims a 'new' material, PLAS,
635       was defined. PLAS has all the characteristics of SCIN but is
636       non-sensitive. No PMTs or PMT holders have been modelled.
637       (*)-The interfaces on CAT where the PMTs should be located are
638           slightly different from the real case.
639    
640  11 February 2002, Bari  11 February 2002, Bari
641    
642  MACRO CLEAN-UP  MACRO CLEAN-UP
# Line 116  TRACK COMMAND CALLED BY GPGARIN Line 738  TRACK COMMAND CALLED BY GPGARIN
738    
739  TRD IONIZATION ENERGY LOSS GENERATED NOW BY GARFIELD  TRD IONIZATION ENERGY LOSS GENERATED NOW BY GARFIELD
740     To generate the ionization in the TRD straw tubes the HEED program     To generate the ionization in the TRD straw tubes the HEED program
741     interfaced by GARFIELD is used (GEANT does not simulate the ionization     interfaced by GARFIELD is used (GEANT does not correctly simulate
742     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
743     the particle in the gas and then passes the coordinates, translated in     tracks the particle in the gas and then passes the coordinates,
744     the DRS, to GARFIELD. The GARFIELD subroutines are called by GPUTRD.     translated in the DRS, to GARFIELD. The GARFIELD subroutines are
745     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
746     variables EGARTRD and NGARTRD of the CWN-tplu.     are stored in the variables EGARTRD and NGARTRD of the CWN-tplu.
747    
748   1 May 2001, Bari   1 May 2001, Bari
749    
# Line 200  NEW SEQUENCES ADDED: $XPRINTPLOT,$PRINTP Line 822  NEW SEQUENCES ADDED: $XPRINTPLOT,$PRINTP
822     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:
823     NVTRD has been forced to 2 for compatibility with GPDTRD.     NVTRD has been forced to 2 for compatibility with GPDTRD.
824    
 3 april 2001, Bari  
   
   
825  28 march 2001, Bari  28 march 2001, Bari
826    
827     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
# Line 666  CAS dimensions and positions modified Line 1285  CAS dimensions and positions modified
1285  CAT dimension modified  CAT dimension modified
1286     CAT dimensions modified after PAOLO simulation and Tom design. Here     CAT dimensions modified after PAOLO simulation and Tom design. Here
1287     are the new dimensions:     are the new dimensions:
1288          DATA CATB/22.05,21.05,0.55/      DATA CATB/22.05,21.05,0.55/
1289          DATA CATL/15.75,6.675,0.55/          DATA CATL/15.75,6.675,0.55/
1290          DATA CATT/6.3,14.375,0.55/          DATA CATT/6.3,14.375,0.55/
1291          DATA ZCAT/75.75/          DATA ZCAT/75.75/

Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.28

  ViewVC Help
Powered by ViewVC 1.1.23