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

Legend:
Removed from v.3.3  
changed lines
  Added in v.3.27

  ViewVC Help
Powered by ViewVC 1.1.23