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

Diff of /gpamela/history/v_100.txt

Parent Directory | Revision Log | View Patch Patch

-revision 3.21 by pam-ba,
Tue Nov 28 10:26:15 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.20 2006/11/16 18:45:29 pam-ba 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
  #
-Line 74
+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 by hand based on the positions given by
+ 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 because the thicknesses of the mylar and of the glue were
+ lightly approximated.
- not considered in that document. Also the height of S3 has been modified as
- suggested from the Naples PAMELA collaborators.
  ND GEOMETRY COMPLETED.
  The He3 and the polyethylene (CH2) have been added to fill
-Line 103 
 TOP PLATE GEOMETRY ADDED.
+Line 195 
 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
-Line 124 
 SPHE and ND data card bugs fixed: the de
+Line 216 
 SPHE and ND data card bugs fixed: the de
  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
-Line 136 
 parameters of the PAMELA sub-detectors"
+Line 228 
 parameters of the PAMELA sub-detectors"
  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
-Line 162 
 CAL HIT STRUCTURE BUGS FIXED
+Line 254 
 CAL HIT STRUCTURE BUGS FIXED
     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
-Line 205 
 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 222 
 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 233 
 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 241 
 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
-Line 254 
 TOF SCINTILLATOR PADDLES UPDATED
+Line 346 
 TOF SCINTILLATOR PADDLES 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 275 
 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 293 
 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 302 
 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 337 
 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 353 
 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 381 
 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 389 
 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 402 
 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 446 
 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 472 
 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 502 
 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 529 
 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 1193 
 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.21
 


changed lines


 
Added in v.3.28
 Legend:



Removed from v.3.21
 


changed lines


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

	ViewVC Help
Powered by ViewVC 1.1.23