/[PAMELA software]/gpamela/history/v_100.txt

Diff of /gpamela/history/v_100.txt

Parent Directory | Revision Log | View Patch Patch

-revision 3.10 by cafagna,
Mon Apr 10 11:07:43 2006 UTC
+revision 3.28 by pamela,
Tue Jan 29 18:25:18 2008 UTC
 Line 1
+ # Revision 3.21  2006/12/05 18:45:29  bruno
+ # Update of the CARD geometry
+ # Simulated also legs, table and magnetic screen
  #
- # $Id: v_100.txt,v 3.9 2005/12/14 03:34:40 cafagna Exp $
+ # $Id: v_100.txt,v 3.27 2007/10/02 17:58:32 cafagna Exp $
  #
  # $Log: v_100.txt,v $
+ # Revision 3.27  2007/10/02 17:58:32  cafagna
+ # Power law spectra introduced
+ #
+ # Revision 3.26  2007/03/30 15:24:54  cafagna
+ # 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
+ #
+ # Revision 3.25  2006/12/17 14:53:18  cafagna
+ # Full review of the CARD geometry. S1 table added
+ #
+ # Revision 3.20  2006/11/16 18:45:29  pam-ba
+ # Simulated an aluminum container for S4
+ #
+ # Revision 3.19  2006/11/16 12:01:13  pam-ba
+ # S4 dimensions corrected.
+ #
+ # Revision 3.18  2006/11/10 11:39:35  pam-ba
+ # S2 and S1 z-positions corrected, He3 and plystyrene mixture added, Top Plate geometry simulated and titanium mixture added.
+ #
+ # Revision 3.17  2006/10/13 16:36:59  pam-ba
+ # Added a new material, the cadmium, for ND
+ #
+ # Revision 3.16  2006/10/12 11:11:21  pam-ba
+ # ND geometry updated.
+ #
+ # Revision 3.15  2006/10/02 11:17:30  pam-ba
+ # NDET 'SPHE' data card meaning changed. Now it eliminates the whole PAMELA container.
+ #
+ # Revision 3.14  2006/06/30 15:38:16  pam-ba
+ # 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)
+ #
+ # Revision 3.13  2006/06/05 13:56:17  pamela
+ # Gigantic resonance added for gamma enetering in the calorimeter absorber
+ #
+ # Revision 3.12  2006/05/18 10:52:32  pam-ba
+ # TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
+ #
+ # Revision 3.11  2006/05/11 23:53:15  cafagna
+ # More bugs fixed in the CALO ntple structure filling
+ #
+ # Revision 3.10  2006/04/10 11:07:43  cafagna
+ # GEN data card updated, ZDGEN added
+ #
  # Revision 3.9  2005/12/14 03:34:40  cafagna
  # An update of the history and inform readme files.
  #
-Line 42
+Line 86
  #CMZ :  1.00/03 30/04/96  12.23.59  by  Francesco Cafagna
  #CMZ :  1.00/02 05/04/96  15.31.25  by  Francesco Cafagna
  #CMZ :  1.00/01 28/11/95  18.51.23  by  Francesco Cafagna
  #-- Author :    Francesco Cafagna   28/11/95
+ Jan 2008, Bari
+ PRESSURIZED CONTAINER GEOMETRY:
+    Pressurized container geometry has been updated: top container TSPH
+    now has been placed at the maximum distance, quoted in the ICD,
+    from S1: 5cm. This distance corresponds to the maximum one in
+    flight due to the cover deformation due to the internal
+    pressure. Three thermal cover layers (nylon6, mylar and aluminium,
+    respectively in the new volumes COV1, COV2 and COV3) have been
+    added to the top container.
+ NEW MATERIALS AND TRACKING MEDIA:
+    Two new materials have been implemented: AMg6M mixture for volume
+    TSPH, and Nylon6 mixture for volume COV1.  New tracking medium have
+    been added for Nylon6 Mixture (TNYL), AMg6M Mixture (TAMG) and Iron
+    (TIRO).
+ TOF GEOMETRY:
+    TOF geometry has been updated: Scintillators positions in the
+    GPAMELA reference system have been corrected shifting all
+    scintillator paddle by 23 microns.
+ TRACK NUMBER RECORDED FOR TOF:
+    The TOF structure now register the stack number into the HIT
+    structure. In GUSTEP the stack number is tested, if is greater than
+, i.e. the particle is not the primary, the particle number is set
+    to 1000+IPART and passed to GPUTOF. Accordingly, now IPARTOF
+    variable can assume 2000 values instead of 100.
+ POWER LAW SPECTRA GENERATION:
+    PLAW datacard parameters have been increased to 4. The meaning of
+    the 4th parameter is: extract random kinetic energy if PLAW(4)<0.,
+    or extract random momentum if PLAW(4)>=0..
+ Sept 2007, Bari
+ POWER LAW SPECTRA GENERATION:
+    It is now possible to select a power law spectra. Using the
+    datacard: PLAW, is possible to pass the spectral index and the
+    kinetic energy range. In PAWL(1) the user can pass the signed power
+    law index, i.e. PAWL -2.5, will generate a power law spectrum
+    according to exp(-2.5). Please note that the range is in kinetic
+    energy. In gukine the generated particle charateristics will be
+    retrived and the mass calculated from this infos. The KINE data
+    card assigned momentum will be overriden as well.
+ PATHLENGHT CALCULATION MODIFIED IN TOF:
+    The pathlength calulation saved in the TOF hit structure has now
+    been modified. Insted of the pathlenght inside each paddle the
+    average pathlenght from the generation vertex to the hit in the
+    paddle, is saved. In this way it will be possible to calculate the
+    beta using this variable.
+ ERROR IN MATERIAL ASSIGNMENT FOR NEUTRON DETECTOR:
+    The material assigned to the NDTU volume was wrong. It was
+    scintillator, IRON has been asigned now.
+ March 2007, Bari
+ CALO GEOMETRY:
+    An error was found in the original document describing PAMELA
+    geometry.  The shift along X axis of the X view planes are swapped,
+    i.e. they don't start with X=+0.05 but with X=-0.05.  The sequence
+    is: 0.05m -0.05, 0.05, -0.05 and so on.
+ MAKEFILE FIXED:
+    Option -fdebug_kludge has been discharged
+ GASATT WRNING FIXED:
+    There was a warning message from a graphic routine. The GSATT calls has been commented out.
+ December 2006, Bari
+ CARD GEOMETRY:
+ Positions, dimensions and quotes have been updated. Indeed, the aluminium table
+ under S1, legs included, and the magnetic screen between the CARD and S2 have
+ been added.
+ #-------------------------------------------------------------------------
+ November 2006, Bari
+ S4 GEOMETRY UPDATED:
+ The x and y dimensions have been corrected because they were wrong.
+ A container of aluminum has been added having walls of thickness equal
+ to 0.1 cm.
+ TOF POSITIONS DEFINITELY UPDATED:
+ S2 and S1 have been positioned again, after the simulation of the top plate.
+ Before now, the positions were put ad hoc based on the positions given by
+ Sergio Ricciarini in the document 'Main geometrical parameters of the PAMELA
+ sub-detectors' released by O. Adriani, L. Bonechi, E. Mocchiutti and
+ S. Ricciarini on the 20th of December 2005.. These positions were
+ lightly approximated.
+ ND GEOMETRY COMPLETED.
+ The He3 and the polyethylene (CH2) have been added to fill
+ the volumes NDTI and NDPB in the neutron detector. Leonov gave us the
+ details of these materials.
+ TOP PLATE GEOMETRY ADDED.
+ The top plate geometry has been simulated as a N2 box that has the same
+ dimensions of CATA. Inside it, at the right positions, there are S2 and CATA.
+ The real top plate is a rectangular window having a thickness of 0.5 cm
+ (volume TPTL emptied by the volume TPCV). This volume is filled with a mixture
+ of titanium (90%), aluminum (6%) and vanadium (4%). Other parts of titanium
+ are the four volumes TPTU, which are positioned at the corners of TPTL at the
+ same height of CATA, and the four volumes TPTM, which are positioned at the
+ corners of TPTL, at the same height of S2.
+ The geometry  of the top plate is based on a simplified version of the CAD
+ drawings.
+ October 2006, Bari
+ The geometry of the neutron detector has been updated. Some dimensions of some
+ volumes have been corrected and the aluminum cover and the aluminum boxes to
+ put the cables have been added. A new material, the cadmium, has been also
+ defined.
+ September 2006, Bari
+ SPHE and ND data card bugs fixed: the definition of the ND data card,
+ missing in the subroutine gpgeo.F, has been added; the meaning of the SPHE
+ data card has been changed. Before the correction the data card:
+ NDET 'SPHE' was used to delete the spherical top shell to substitute it with
+ a flat one.Now NDET 'SPHE' eliminates the whole container of PAMELA.
+ June 2006, Bari
+ The center of the scintillator planes S22Y (variable ZPAMS22Y in gpdgeo.inc)
+ and S12X (variable ZPAMS12X in gpdgeo.inc) has been positioned at the
+ nominal height as measured in PAMELA (See the document: "Main geometrical
+ parameters of the PAMELA sub-detectors" by O. Adriani, L. Bonechi,
+ E. Mocchiutti, S. Ricciarini, 20 December 2005). Follows that the positions
+ of S21Y and S12X are higher than those in the cited document due to the fact
+ that in GPAMELA the thickness of the mylar has been considered while in the
+ document it has been neglected.
+ May 2006, Bari & Tor Vergata
+ GIGANTIC RESONANCE FOR NEUTRON DETECTOR ADDED
+    Routines to simulate the gigantic resonance of gammas in Tungsten
+    have been added.  The GPGIG routine is called in GUSTEP if a gamma
+    enter the calorimeter absorber.  This is the steering routine to
+    simulate the production of neutrons from gigantic resonance.  It
+    does checks on STEP lenght. If the range is smaller than the other
+    selected for that step, it does generate the neutron and stops the
+    gamma. Please note that the neutron has now a new particle
+    number. This is to tag the gigantic resonance neutrons.
+ May 2006, Bari & Florence
+ CAL HIT STRUCTURE BUGS FIXED
+    The maximum number of hit is now different for the two hit
+    structures: CALST and CALI. Vectors inizialization and HBOOK
+    ntple booking have been updated. The GPDCAL routine has been fixed
+    so to handle the case in wich hits stored are more than the maximum
+    number of hit.
+    In this case in the ntple up to the maximum number of hits will be stored.
+ April 2006, Bari
+ TOF GEOMETRY AND POSITIONS UPDATED AND NEW MIXTURES ADDED
+    The TOF geometry has been modified. The following boxes have been
+    added: POL1, POL2 and POLY made of polystyrene, S11M, S12M, S21M,
+    S22M, S31M and S32M made of mylar, S1A, S2A and S3 made of air and
+    S1 and S2 made of aluminum. Each scintillator paddle has been put
+    in his mylar box and the other materials: air, polystyrene, and
+    aluminum have been added at their nominal positions.  According to
+    Naples people the araldite glue has been simulated has an air
+    gap. For this work two new materials: the Mylar (MYLAR) and the
+    polystyrene (POLYSTYRENE) with a density of 35 g/l have been
+    defined as a mixture.  The positions of the three bottom
+    scintillator planes that contain respectively the S12X, S22Y and
+    S32X paddles have been regulated according on their official
+    positions in PAMELA.
  Mar 2006, Bari
-Line 67 
 NEW GPKEY ADDED: GPCALOR
+Line 297 
 NEW GPKEY ADDED: GPCALOR
     is set to true in GPDAT if the data card: HPAK, is set to
     'GCAL'. The gpkey.inc has been update accordingly.
  NEUTRON DETECTOR ADDED. NEW DIR: GPND
     The neutron detector has been added. At the moment it is just the
-Line 84 
 GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR,
+Line 314 
 GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR,
     http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html
     This package should be distributed with the GEANT library but is
     not up to date. I did download the latest release and stored into
     gpcalor directory of the gpamela tree.
     Then I did clean up the code substituting the explicit inclusion of
     the commons with a #include cpp directive. In parallel I did
     extract the commons to include files having the same common name. I
-Line 95 
 GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR,
+Line 325 
 GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR,
     NOTE: There are still problem due to different common sizes. In
     particular the common MICFIL is maller in the geant library
     libgeant.a . There the subroutines: gmorin, gmxsec, gmplxs, are
     present and linked using a wrong version of the common. This still needs to be debuged.
     NOTE2: The auxiliary files with the cross sections: chetc.dat.gz
     and xsneut.dat.gz, have been added to the aux directory and moved
     to the working directory, i.e. GPAMELA_BIN. The GCALOR routine will
-Line 103 
 GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR,
+Line 333 
 GCALOR PACKAGE ADDED. NEW DIRs: GPCALOR,
     searched there at first, then in the working directory. A fool
     proof policy has to be implemented to avoid problem with
     synchronization fo these files.
- The GCALOR package
+ The GCALOR package
+ June 2005, Bari
+ TOF SCINTILLATOR PADDLES UPDATED
+    The dimensions and the number of the scintillator paddles for each
+    TOF planes have been updated.
  May 2005, Bari
  Some updates on the latest modification done in the past year.
  NEW DATA CARD ADDED: HFSF
     To define a policy for the random number initial seeds
     definition. Using this card is possible to override GEANT seeds
     defined via NRDM card. The policy is selected according to the
     values:
     - 1: The seeds are initialized to the initial values found in a user
          defined file or the default file: INPUTSEED.DAT
     - 2: The seeds are initialized to the final values found in a user defined
          file or the default file: INPUTSEED.DAT
     The case 1 must be used in case the user needs to reproduce the
-Line 129 
 NEW DATA CARD ADDED: HFSF
+Line 367 
 NEW DATA CARD ADDED: HFSF
     binary file and pass the filename to the program using the *FLSF
     data card. In case the user file is not specified the default
     INPUTSEED.DAT will be used.
     The case 2 must be used in case the user needs to chain several
     GPAMELA run and likes to be sure he is starting the random
     generator using the right sequence. In this case the user must
-Line 147 
 NEW DATA CARD ADDED: *FSFI
+Line 385 
 NEW DATA CARD ADDED: *FSFI
     filled with his own seeds. For this purpose an utility program:
     writeseeds.f, has been added in the aux directory.  In case the
     *FSFI card is not specified the default values: 24 and INPUTSEEDS.DAT, will
     be used as LUN and file name respectively.
  NEW DATA CARD ADDED: *LSFI
     Using this card the user can specify the logical unit and name of
-Line 156 
 NEW DATA CARD ADDED: *LSFI
+Line 394 
 NEW DATA CARD ADDED: *LSFI
     run. This file is a FORTRAN binary one. This file can be used as
     input one specifying it in the *FSFI data card of the next GPAMELA
     run.  In case the *LSFI card is not specified the default values: 26
     and HBOOKFILENAME.DAT (as sepified in *HFI), will be used as LUN
     and file name respectively.
  NEW UTILITY PROGRAMS ADDED: writeseeds.f, readseeds.f
     These new programs have been added in the aux directory. Using these a
-Line 191 
 TRD GEOMETRY AND CALIBRATION REVIEWD
+Line 429 
 TRD GEOMETRY AND CALIBRATION REVIEWD
     has been revied and new sets of calibration constants calculated
     comparing the beam test data with the GPAMELA results. The new
     constants are about 3% larger than the previous ones.
  TRACKER GEOMETRY REVIEWED. NEW VOLUME DEFINED: THBP, TPAS, TPAI
     Thanks to Lorenzo Bonechi for the drawings and explanations. Now the
     hybrd cards have been put into the simulation and the geometry updated
     considering the dead zones in the silicon detectors. The hybrid zone
     has been simulated as well. At the moment the hybrid is simulated as
     a G10 plates. The full height of the tracker magnet has been
     reviewed as well.
     The tracker ladder is now simulated inside a nitrogen box: TPAS,
     placed inside an aluminum frame: TRPB. Each silicon ladder has been
-Line 207 
 TRACKER GEOMETRY REVIEWED. NEW VOLUME DE
+Line 445 
 TRACKER GEOMETRY REVIEWED. NEW VOLUME DE
     smaller silicon detector: TPAI, has been placed inside the larger
     silicon block TRSL. In this way the subdivided silicon ladder can
     be upgraded with an indipendend roto-translation for each sensor.
     The TRPB aluminum frame has been enlarged to fit the external
     magnet canister frame.
     The last plane has been flipped with a 180 degree rotation around
     the X axis.
  TRACKER HIT STRUCTURE REVIEWED
     Taking into account the new version of the tracker geometry, the hit
     structure for this detector has been revied.
  CALORIMETER GEOMETRY REVIEWED
     Marco Albi reviewed the calorimeter dimentions and positioning.
 March 2004, Bari
-Line 235 
 NON-REPRODUCIBILITY PROBLEM OF A GPAMELA
+Line 473 
 NON-REPRODUCIBILITY PROBLEM OF A GPAMELA
     the GARFIELD code for details). Follows that different runs produce
     different results without changing the initial seeds. To have identical
     results in different runs, the GARFIELD program has to start typing the
     noRNDM_initialisation switch. To avoid of specifying this switch
     by the user,
     the GARFIELD package has been upgraded with a patch. In this way the problem
     is partially solved because, now, the initial seeds of the random generators
-Line 243 
 NON-REPRODUCIBILITY PROBLEM OF A GPAMELA
+Line 481 
 NON-REPRODUCIBILITY PROBLEM OF A GPAMELA
     activated by the user for changing the initial seeds in the GPAMELA program.
     Work is in progress for a more general correction of this problem.
     Please, use the updated GARFIELD code released with the CVS version v4r1
     to fix this problem.
  RNDM ROUTINE REPLACED BY THE GRNDM ROUTINE IN GPXTR AND NPOISS.
-Line 256 
 BUG FOUND AND FIXED: the set and detecto
+Line 494 
 BUG FOUND AND FIXED: the set and detecto
     the GPSED common when the set and the detector ZEBRA banks are filled
     during a run. In general the values of the set and detector addresses
     depend on the number of active detectors in a given run. ISCAL=12 and
     IDCASI=1 are only right when all the detectors of GPAMELA are active.
 December 2003, Bari
-Line 300 
 BUG FOUND AND FIXED: the set and detecto
+Line 538 
 BUG FOUND AND FIXED: the set and detecto
     -ICAPLANE(i) contains the number of hit plane;
     -ICASTRIP(i) contains the number of hit strip;
     -ICAMOD(i) can assume different values based on the number of times and
-               positions in which a `long strip' has been hit.
+           positions in which a `long strip' has been hit.
     -ENESTRIP(i) contains the deposited energy in the hit strip;
     where i is the number of hit (1<i<4224).
     Note that in the calorimeter each hit is filled at the end of the event and
-Line 326 
 BUG FOUND AND FIXED: the set and detecto
+Line 564 
 BUG FOUND AND FIXED: the set and detecto
     The energy deposited in the calorimeter is calculated in GeV.
     To store the hits in the calorimeter the subroutine GSAHIT is used instead of
     GSCHIT.
     To retrieve the hit structure the call to the routine GPRHIT is done instead
     of a call to the GFHITS subroutine.
 February 2003, Bari
  BUG found:
     DCUTEAER, DCUTEAL, DCUTECE, DCUTECP, DCUTEFE, DCUTEG10C, DCUTEG10, DCUTEKAP,
     DCUTEN2G, DCUTEROA, DCUTESCIN, DCUTESICA, DCUTETRAD, DCUTEW2,
-Line 356 
 BUG found:
+Line 594 
 BUG found:
     Updated the relative positions of the calorimeter planes. The corrected
     shifting are:
     first view: (Dxo,Dyo)=(0.10,0.05) cm
     second view: (Dxo,Dyo)=(-0.05,0.10) cm
     third view: (Dxo,Dyo)=(-0.10,-0.05) cm
     fourth view: (Dxo,Dyo)=(0.05,-0.10) cm
 November 2002, Bari
  CAS detectors distances modified
     The distances between the CAS detectors have been modified based on the
     latest CAD drawings.
 November 2002, Bari
-Line 383 
 CALORIMETER geometry upgrade
+Line 621 
 CALORIMETER geometry upgrade
     CALB and CALS. CALS is an aluminum box of dimensions: 48.4*48.4*21.278 cm^3,
     having side-walls 1 cm thick and a bottom of 1 mm. The real box is more
     complicated and the configuration of the bottom should be upgraded if we want
     a reliable description of the event in the S4 scintillator.
 October 2002, Stockholm
-Line 1047 
 CAS dimensions and positions modified
+Line 1285 
 CAS dimensions and positions modified
  CAT dimension modified
     CAT dimensions modified after PAOLO simulation and Tom design. Here
     are the new dimensions:
-         DATA CATB/22.05,21.05,0.55/
+     DATA CATB/22.05,21.05,0.55/
          DATA CATL/15.75,6.675,0.55/
          DATA CATT/6.3,14.375,0.55/
          DATA ZCAT/75.75/

 Legend:



Removed from v.3.10
 


changed lines


 
Added in v.3.28
 Legend:



Removed from v.3.10
 


changed lines


 
Added in v.3.28
-Removed from v.3.10
+Added in v.3.28

	ViewVC Help
Powered by ViewVC 1.1.23