Import Geant4 4.0.2 source tree

Author: gcosmo

ReleaseNotes/Patch4.4.0-2.txt

@@ -0,0 +1,157 @@
+
+                 Geant4 4.0 - patch-02 Release Notes
+                 -----------------------------------
+
+                                                      28th February 2002
+
+List of fixes included in this public patch to release 4.0, to be added on
+top of patch-01:
+
+  o Config:
+    ------
+    - liblist.c: added flags __CYGWIN__ and __CYGWIN32__ for triggering proper
+      parsing algorithm on Windows systems, also on most recent versions of
+      Cygwin.
+    - WIN32-VC.gmk: use g77 instead of f77 as default fortran compiler.
+    - sys/HP-aCC.gmk: removed include path to
+      global/management/include/std/wrappers, which was unused and not
+      necessary.
+
+  o Decay processes:
+    ---------------
+    - Removed a physics table from G4Decay.
+
+  o Electromagnetic processes:
+    -------------------------
+    - Fixed bug in G4Ve(h)EnergyLoss, G4VMuEnergyLoss:
+      o Adopt MinDeltaCut in-range instead of in-energy !
+    - G4VEnergyLoss:
+      o Set subSecFlag to false. NO default subcutoff generation.
+      o Set finalRange to 1mm as default (the actual value is computed in
+        G4VeEnergyLoss, G4VhEnergyLoss and G4VMuEnergyLoss).
+
+  o Event:
+    -----
+  - Optimised implementation of SetNewValue() in G4ParticleGunMessenger and 
+    G4GeneralParticleSourceMessenger to allow porting on WIN32-VC for optimised
+    compilation. 
+
+  o Geometry:
+    --------
+  - Management:
+    o G4GeometryManager.cc: fixed handling of ios manipulators and precision
+      to reset to their original state in G4GeometryManager::ReportVoxelStats().
+      Fixes problem report #343.
+  - Volumes:
+    o G4GeometryMessenger.cc: fixed local variable name not to clash with
+      class data member. 
+  - solids/Boolean:
+    o G4ReflectionFactory.cc: fixed bug causing repetition of daughter volumes.
+    o G4ReflectedSolid[.hh.cc]: 
+      - Defined new methods Set/GetTransform3D() and Set/GetDirectTransform3D() 
+        for defining/accessing the reflected transformation.
+    o Moved utility methods for G4AffineTransform to protected.
+  - solids/BREPS:
+    o Bug fixes to G4BREPSolidPolyhedra: 
+      - Added enhancement allowing to define RMIN and RMAX to be zero with the 
+        exception of forbidden two consecutive zero RMAX values. 
+      - Refactoring in constructor gained two new private methods for creation
+        of trapezoidal and triangular side planes of G4BREPSolidPolyhedra. 
+      - Rewritten code for delta phi >= 2*PI case with full classification of
+        possible cases. 
+      - Fixed sense of end planes. 
+      - Fixed delta phi < 2*PI case where some old code remained after bug
+        fix #266. 
+      - Added new private method for creation of planar surfaces gained 
+        after code refactoring in constructor. 
+      - Fixed bug for problem report #354 by adding the code handling the
+        cases where user supplies for the same z-section different RMIN or
+        RMAX values.
+      - Added more preconditions checks 
+      - Added code forbidding "teeth" in polyhedra.
+    o G4ConvexHull.hh: moved initialisation of 'next' to constructor's body. 
+  - solids/CSG:
+    o G4Sphere.cc: bug-fix in Inside() function. Fixes problem report #216. 
+    o G4Para.cc: fixed bug in G4Para::Inside() in calculation of xt,
+      according to problem report #357.
+  - solids/STEPinterface:
+    o G4RepresentationRelationshipCreator.cc: fixed unsafe usage of
+      '<' operator. 
+
+  o Hadronic processes:
+    ------------------
+    o Fixes in pre-equilibrium and de-excitation generator models.
+    o Fixed bug in doppler broadening in neutron-transport model:
+      - Fixed sign error in calculation of effective luminosity.
+    o Stopping: trivial bug-fix in pi-absorption at rest. 
+    o Util:
+      - Trivial bug-fix, affecting K-absorption at rest.
+      - G4Nucleus: moved definition of constructors/destructor from header to
+        source implementation file. Unresolved symbols were detected at link
+        time on SUN-CC for static constants used to initialise members in
+        inline constructors. 
+    o CHIPS model: fixed accuracy problem.
+    o Optimised implementation in G4Fancy3DNucleus.cc to allow porting on
+      WIN32-VC for optimised compilation. 
+
+  o Interfaces:
+    ----------
+    - Fixed for odd behavior of G4UIcsh in Windows systems.
+
+  o Materials:
+    ---------
+    - Optimised implementation of GetIndex() for fIndexTable in G4Material,
+      G4Element and G4Isotope.
+
+  o Particles:
+    ---------
+    - Fixed bug in the decay table for G4KaonMinus.
+
+  o Track & Tracking:
+    ----------------
+    - track:
+      o Added G4Track::CopyTrackInfo().
+      o Optimised G4ParticleChange::AddSecondary() to properly use touchable
+        handles.
+      o Fixed a bug of changing parent weight in
+        G4ParticleChange::UpdateStepForAlongStep().
+      o Optimised implementation in G4Track::GetVelocity().
+    - tracking:
+      o G4SteppingManager2.cc:
+        - fix for extra steps when the ApplyCuts flag is set 
+        - G4SteppingManager.hh: changed scope for 'SizeOfSelectedDoItVector'
+          to allow porting on WIN-VC compiler.
+      o Fix in G4SteppingManager::SetInitialStep() for applyCuts. 
+
+  o Visualization:
+    -------------
+    - HepRep: added check for missing hierarchy layers in HepRepXMLWriter. 
+    - Bug fix in modeling: removed Begin/EndPrimitives from G4HitsModel and
+      G4TrajectoriesModel. This was causing nested calls to Begin/EndPrimitives
+      (and thus to glNewList/EndList in OpenGL stored mode), since the
+      trajectories and hits Draw methods themselves have Begin/EndPrimitives.
+
+  o Examples:
+    --------
+    - extended/electromagnetic/testem1,3,4,5:
+      o Improvements in PrimaryGenerator and macros
+    - novice/N03:
+      o New macros for tutorials/exercises.
+
+ ----------------------------------------------------------------------------
+
+                             Technical Notes
+                             --------------- 
+
+  o  This patch should be applied on top of "patch-01" for release 4.0.
+  o  Technical notes distributed for release 4.0 are also applicable and
+     valid for this patch.
+
+The code and rebuilt binary libraries for release 4.0 are available through
+our "Source Code" Web page:
+- see our <A href="http://cern.ch/geant4">Geant4 Home Page</A>.
+
+Please refer to <A
+href="http://cern.ch/geant4/G4UsersDocuments/Overview/html/index.html">Geant4
+User Documentation</A> for further information about using Geant4.
+

config/History

@@ -1,5 +1,5 @@
-$Id: History,v 1.177 2001/12/14 08:37:39 gcosmo Exp $
-$Name: geant4-04-00 $
+$Id: History,v 1.180 2002/02/07 10:18:32 gcosmo Exp $
+$Name: geant4-04-00-patch-02 $
 
 -------------------------------------------------------------------
 
@@ -18,6 +18,15 @@ committal in the CVS repository !
      * Reverse chronological order (last date on top), please *
      ----------------------------------------------------------
 
+07th February 2002 Gabriele Cosmo (config-V04-00-01)
+- liblist.c: added flags __CYGWIN__ and __CYGWIN32__ for triggering proper
+  parsing algorithm on Windows systems, also on most recent versions of Cygwin.
+- WIN32-VC.gmk: use g77 instead of f77 as default fortran compiler.
+
+25th January 2002 Gabriele Cosmo (config-V04-00-00)
+- sys/HP-aCC.gmk: removed include path to global/management/include/std/wrappers,
+  which was unused and not neccessary.
+
 10th December 2001 Gabriele Cosmo (config-V03-02-15)
 - architecture.gmk:
   o updated list of supported compilers

config/liblist.c

@@ -1,4 +1,4 @@
-/* $Id: liblist.c,v 1.9 1999/11/19 18:01:40 morita Exp $ */
+/* $Id: liblist.c,v 1.11 2002/02/06 18:44:36 radoone Exp $ */
 
 /*
 Given a "libname.map" file on standard input and a list or directory
@@ -121,7 +121,7 @@ int main (int argc, char** argv) {
   int i,optl=0,swapping,c,rargc;
   FILE *fp;
 
-#ifdef _WIN32
+#if defined ( _WIN32 ) || defined ( __CYGWIN__ ) || defined ( __CYGWIN32__ )
   char *ntg4tmp=0,*ntg4tmp1=0;
   int nti; 
 #endif
@@ -263,7 +263,7 @@ int main (int argc, char** argv) {
 
   if(optl)fprintf(stderr,"  Reading dependency files...\n");
 
-#ifdef _WIN32
+#if defined ( _WIN32 ) || defined ( __CYGWIN__ ) || defined ( __CYGWIN32__ )
       ntg4tmp=getenv("G4TMP");
       if ( ! ntg4tmp ) 
         {
@@ -278,7 +278,7 @@ int main (int argc, char** argv) {
       fp=fopen(rargv[i],"r");
       fgets(buffer,BUFSIZE,fp);
 
-#ifdef _WIN32
+#if defined ( _WIN32 ) || defined ( __CYGWIN__ ) || defined ( __CYGWIN32__ )
       ptr=strchr(ntg4tmp1,':');
       if ( ptr ) *(ptr+1)='\0';
 
@@ -362,7 +362,7 @@ int main (int argc, char** argv) {
 	    fgets(buffer,BUFSIZE,fp);
 	    bufferPtr=buffer;
 
-#ifdef _WIN32
+#if defined ( _WIN32 ) || defined ( __CYGWIN__ ) || defined ( __CYGWIN32__ )
 	    while ( ptr=strchr(buffer,'\\') ) *ptr='/';
 
 	    while (ntg4tmp1 &&  (ptr=strstr(buffer,ntg4tmp1)) )
@@ -376,7 +376,7 @@ int main (int argc, char** argv) {
       }
     }
 
-#ifdef _WIN32
+#if defined ( _WIN32 ) || defined ( __CYGWIN__ ) || defined ( __CYGWIN32__ )
       free(ntg4tmp1);
 #endif
 

config/sys/HP-aCC.gmk

@@ -39,7 +39,6 @@ ifeq ($(G4SYSTEM),HP-aCC)
     OGLLIBS   += -L/opt/graphics/common/lib -lXwindow -lhpgfx -lXhp11 -ldld
   endif
   CPPFLAGS += -DG4NOT_ISO_DELETES
-  CPPFLAGS += -I$(G4BASE)/global/management/include/g4std/wrappers/$(G4SYSTEM)
   ifdef G4USE_OSPACE
     CPPFLAGS += -DOS_HPUX_10_20 -DOS_NEW_CHECK -DOS_STL_ASSERT
     CPPFLAGS += -DOS_NO_WSTRING -DOS_NO_ALLOCATORS -D__HPACC_NOEH

config/sys/WIN32-VC.gmk

@@ -40,7 +40,7 @@ ifeq ($(G4SYSTEM),WIN32-VC)
 
   AR := LIB
   ECHO := echo -e
-  FC := f77
+  FC := g77
 #  FCLIBS := -lf2c
   ifndef OGLFLAGS
     OGLFLAGS  :=

examples/extended/electromagnetic/TestEm1/History

@@ -1,4 +1,4 @@
-$Id: History,v 1.15 2001/11/29 11:28:06 maire Exp $
+$Id: History,v 1.17 2002/02/06 14:06:05 maire Exp $
 -------------------------------------------------------------------
 
      =========================================================
@@ -14,6 +14,14 @@ track of all tags.
      ----------------------------------------------------------
      * Reverse chronological order (last date on top), please *
      ----------------------------------------------------------
+06-02-02 mma (testem1-V04-00-00)
+     
+06-12-01  mma
+- In PrimaryGeneratorAction, possibility to fix the vertex at front face, 
+  and to randomize the beam spot.
+- command /gun/rndm in PrimaryGeneratorMessenger
+- similar functionality in geant3: gukine.F
+- change output precision in RunAction (+uglast.F)
 
 29-11-01  mma (testem1-V03-02-02) 
 - EventAction: DrawTrack "none"

examples/extended/electromagnetic/TestEm1/README

@@ -1,4 +1,4 @@
-$Id: README,v 1.4 2001/10/26 12:51:21 maire Exp $
+$Id: README,v 1.5 2001/12/07 11:49:08 maire Exp $
 -------------------------------------------------------------------
 
      =========================================================
@@ -42,7 +42,11 @@ $Id: README,v 1.4 2001/10/26 12:51:21 maire Exp $
  	in the PrimaryGeneratorAction class, and can be changed via the G4 
  	build-in commands of ParticleGun class (see the macros provided with 
  	this example).
- 	
+	
+ 	In addition one can choose randomly the impact point of the incident
+ 	particle. The corresponding interactive command is built in
+ 	PrimaryGeneratorMessenger class.
+	 	
  	A RUN is a set of events.
 
 

examples/extended/electromagnetic/TestEm1/TestEm1.cc

@@ -21,8 +21,8 @@
 // ********************************************************************
 //
 //
-// $Id: TestEm1.cc,v 1.7 2001/11/29 11:28:06 maire Exp $
-// GEANT4 tag $Name: geant4-04-00 $
+// $Id: TestEm1.cc,v 1.8 2001/12/07 11:49:08 maire Exp $
+// GEANT4 tag $Name: geant4-04-00-patch-02 $
 // 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 
@@ -68,7 +68,7 @@ int main(int argc,char** argv) {
   Em1DetectorConstruction* det;
   runManager->SetUserInitialization(det = new Em1DetectorConstruction);
   runManager->SetUserInitialization(new Em1PhysicsList(det));
-  runManager->SetUserAction(new Em1PrimaryGeneratorAction);
+  runManager->SetUserAction(new Em1PrimaryGeneratorAction(det));
 
 #ifdef Em1NoOptimize  
   #ifdef G4VIS_USE

examples/extended/electromagnetic/TestEm1/TestEm1.out

@@ -1,6 +1,6 @@
 **********************************************
- Geant4 version $Name: geant4-04-00-patch-01 $
-                                (31-Jan-2002)
+ Geant4 version $Name: geant4-04-00-patch-02 $
+                                (28-Feb-2002)
              Copyright : Geant4 Collaboration
 **********************************************
 Visualization Manager instantiating...
@@ -22,7 +22,7 @@ Current available graphics systems are:
 # 
 /run/initialize
 
-***** Table : Nb of materials = 9 *****
+***** Table : Nb of materials = 10 *****
 
  Material:      Air     density:  1.290 mg/cm3  temperature: 273.15 K  pressure:   1.00 atm  RadLength: 285.161 m  
    --->  Element: Nitrogen  N   Z =  7.0   N =  14.0   A =  14.01 g/mole  fractionMass:  70.00 %  Abundance  72.71 %
@@ -44,6 +44,9 @@ Current available graphics systems are:
  Material:  Silicon     density:  2.330 g/cm3   temperature: 273.15 K  pressure:   1.00 atm  RadLength:   9.368 cm 
    --->  Element:  Silicon      Z = 14.0   N =  28.1   A =  28.09 g/mole  fractionMass: 100.00 %  Abundance 100.00 %
 
+ Material: Germanium     density:  5.323 g/cm3   temperature: 273.15 K  pressure:   1.00 atm  RadLength:   2.301 cm 
+   --->  Element: Germanium      Z = 32.0   N =  72.6   A =  72.61 g/mole  fractionMass: 100.00 %  Abundance 100.00 %
+
  Material:     Iron     density:  7.870 g/cm3   temperature: 273.15 K  pressure:   1.00 atm  RadLength:   1.759 cm 
    --->  Element:     Iron      Z = 26.0   N =  55.9   A =  55.85 g/mole  fractionMass: 100.00 %  Abundance 100.00 %
 
@@ -125,6 +128,7 @@ Cut in range
         liquidArgon         2.2 mm         6.05 mum
           Aluminium         2.2 mm         6.05 mum
             Silicon         2.2 mm         6.05 mum
+          Germanium         2.2 mm         6.05 mum
                Iron         2.2 mm         6.05 mum
                Lead         2.2 mm         6.05 mum
             Uranium         2.2 mm         6.05 mum
@@ -135,6 +139,7 @@ Cut in energy
         liquidArgon        9.08 keV        3.98 keV
           Aluminium        9.99 keV        9.97 keV
             Silicon          10 keV        8.45 keV
+          Germanium        31.6 keV        32.7 keV
                Iron        31.1 keV        43.3 keV
                Lead         124 keV        42.8 keV
             Uranium         143 keV        53.4 keV
@@ -164,14 +169,14 @@ Start Run processing.
 Run terminated.
 Run Summary
   Number of events processed : 100
-  User=32s Real=62s Sys=2.3s
+  User=34s Real=37s Sys=2.3s
 
- nb tracks/event   neutral:   27.29   charged:  692.97
- nb  steps/event   neutral:  145.49   charged: 1331.63
+ nb tracks/event   neutral:    27.2900   charged:   692.9700
+ nb  steps/event   neutral:   145.4900   charged:  1331.6300
 
  nb of process calls per event: 
-       eIoni    eBrem    compt     phot     conv  annihil
-     1305.51    24.92   118.20    26.09     1.20     1.20
+          eIoni       eBrem       compt        phot        conv     annihil
+      1305.5100     24.9200    118.2000     26.0900      1.2000      1.2000
 
 --------- Ranecu engine status ---------
  Initial seed (index) = 0

examples/extended/electromagnetic/TestEm1/geant3/phot.dat

@@ -0,0 +1,19 @@
+LIST
+C
+C       Photo electric; gamma 1 MeV
+C 
+CALOR   10 (Germanium) 1. (size) 0. (field)
+TRIG    100000
+KINE    1 (gamma)   0.001 (Ekine)   1. (xBeam)  0.8 (rBeam)
+DEBUG   10  5 1000
+SWIT    0 (draw) 0 (save)
+CUTS    10.e-6 (cutgam)    10.e-6 (cutele)  3*10.e-3 (cutneu/had/muo)
+      2*10.e+3 (bcute/m) 2*10.e+3 (dcute/m)   10.e-3 (ppcutm)
+ABAN    0
+PHOT    1
+COMP    0
+PAIR    0
+BREM    0
+MULS    1
+LOSS    1
+TIME    2=1.