Test for issue 1134 #1139

MarkusFrankATcernch · 2023-06-23T10:24:03Z

BEGINRELEASENOTES

Add test for issue Abnormal material budget of DDCAD-imported files #1134

ENDRELEASENOTES

MarkusFrankATcernch · 2023-06-23T10:30:05Z

Simulation output from ddsim as described here

gdb --args `which python` `which ddsim` --compactFile ../checkout/examples/DDCAD/compact/DD4hep_Issue_1134.xml --enableG4Gun --runType shell
GNU gdb (Ubuntu 12.1-0ubuntu1~22.04) 12.1
Copyright (C) 2022 Free Software Foundation, Inc.
License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>
This is free software: you are free to change and redistribute it.
There is NO WARRANTY, to the extent permitted by law.
Type "show copying" and "show warranty" for details.
This GDB was configured as "x86_64-linux-gnu".
Type "show configuration" for configuration details.
For bug reporting instructions, please see:
<https://www.gnu.org/software/gdb/bugs/>.
Find the GDB manual and other documentation resources online at:
    <http://www.gnu.org/software/gdb/documentation/>.

For help, type "help".
Type "apropos word" to search for commands related to "word"...
Reading symbols from /home/frankm/bin/python...
(No debugging symbols found in /home/frankm/bin/python)
(gdb) r
Starting program: /home/frankm/bin/python /home/frankm/SW/DD4hep_head_dbg.root_vmaster.g4_11.00.b01_MT.G4_units/DD4hep/bin/ddsim --compactFile ../checkout/examples/DDCAD/compact/DD4hep_Issue_1134.xml --enableG4Gun --runType shell
[Thread debugging using libthread_db enabled]
Using host libthread_db library "/lib/x86_64-linux-gnu/libthread_db.so.1".
warning: File "/data/SW/root_vmaster_dbg/lib/libCore.so-gdb.py" auto-loading has been declined by your `auto-load safe-path' set to "$debugdir:$datadir/auto-load".
To enable execution of this file add
	add-auto-load-safe-path /data/SW/root_vmaster_dbg/lib/libCore.so-gdb.py
line to your configuration file "/home/frankm/.config/gdb/gdbinit".
To completely disable this security protection add
	set auto-load safe-path /
line to your configuration file "/home/frankm/.config/gdb/gdbinit".
For more information about this security protection see the
"Auto-loading safe path" section in the GDB manual.  E.g., run from the shell:
	info "(gdb)Auto-loading safe path"
[Detaching after vfork from child process 96003]
[Detaching after vfork from child process 96006]
[Detaching after vfork from child process 96010]
[Detaching after vfork from child process 96012]
[Detaching after vfork from child process 96025]
[New Thread 0x7fffd4ae2640 (LWP 96036)]
PersistencyIO    INFO  +++ Set Streamer to dd4hep::OpaqueDataBlock
DD4hep           INFO  ++ Using globally Geant4 unit system (mm,ns,MeV)
Info in <TGeoManager>: Changing system of units to Geant4 units (mm, ns, MeV).
Info in <TGeoManager::TGeoManager>: Geometry default, Detector Geometry created
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
XMLLoader        INFO  +++ Processing XML file: file:/data/SW/DD4hep_CHECKOUT/checkout/examples/DDCAD/compact/DD4hep_Issue_1134.xml
DD4hep           WARN  ++ STD conditions NOT defined by client. NTP defaults taken.
Detector         INFO  *********** Created World volume with size: 3000 3000 3000
Info in <TGeoManager::SetTopVolume>: Top volume is world_volume. Master volume is world_volume
ASSIMPReader     INFO  +++ tessellated_0     Material: DefaultMaterial   Viualization: NONE
Facet: {{109.753, 41.2677, 1068.83}, {109.803, 40.6582, 1069.38}, {109.803, 40.6582, 1069.38}} is degenerated
Facet: {{57.5019, -6.55271, 1069.14}, {57.5019, -6.55271, 1069.14}, {61.5019, -6.5531, 1069.14}} is degenerated
Facet: {{101.713, 40.6004, 1070.51}, {109.803, 40.6582, 1069.38}, {109.803, 40.6582, 1069.38}} is degenerated
Facet: {{57.5019, -6.55271, 1069.14}, {57.5019, -1.98234, 1066.39}, {57.5019, -6.55271, 1069.14}} is degenerated
Facet: {{565.4, 40.6582, 1069.38}, {565.4, 40.6582, 1069.38}, {565.45, 41.2677, 1068.83}} is degenerated
Facet: {{617.702, -6.5531, 1069.14}, {613.702, -6.5531, 1069.14}, {617.702, -6.5531, 1069.14}} is degenerated
Facet: {{573.491, 40.6004, 1070.51}, {565.4, 40.6582, 1069.38}, {565.4, 40.6582, 1069.38}} is degenerated
Facet: {{617.702, -3.7916, 1067.82}, {617.702, -6.5531, 1069.14}, {617.702, -6.5531, 1069.14}} is degenerated
ASSIMPReader           +++ Read 1 meshes from /home/frankm/SW/DD4hep_head_dbg.root_vmaster.g4_11.00.b01_MT.G4_units/DD4hep/examples/DDCAD/models/STL/DD4hep_Issue_1134.stl
TestShape        INFO  Created successfull shape of type: CAD_MultiVolume [CAD_MultiVolume]
TestShape        ERROR BAD shape type: CAD_MultiVolume <->  Instance test: FAILED
Compact          INFO  ++ Converted subdetector:Shape_STL of type DD4hep_TestShape_Creator 
Info in <TGeoManager::CheckGeometry>: Fixing runtime shapes...
Info in <TGeoManager::CheckGeometry>: ...Nothing to fix
Info in <TGeoManager::CloseGeometry>: Counting nodes...
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CloseGeometry>: Building cache...
Info in <TGeoManager::CountLevels>: max level = 3, max placements = 1
Info in <TGeoManager::CloseGeometry>: 4 nodes/ 5 volume UID's in Detector Geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
Detector         INFO  +++ Patching names of anonymous shapes....
DDG4             INFO +++ Imported 10 global values to namespace:DDG4
Geant4Kernel           OutputLevel:  2
Geant4Kernel           UI:           UI
Geant4Kernel           NumEvents:    10
Geant4Kernel           NumThreads:   0
DDG4             INFO +++  List of sensitive detectors:
Geant4UI         INFO  +++ UI> Install Geant4 control directory:/ddg4/UI/
Geant4Kernel     INFO  ++ Registered global action UI of type dd4hep::sim::Geant4UIManager
Geant4UI         INFO  +++ MagFieldTrackingSetup> Install Geant4 control directory:/ddg4/MagFieldTrackingSetup/
Geant4UI         INFO  +++ DetectorConstructionAction> Install Geant4 control directory:/ddg4/DetectorConstructionAction/
Geant4UI         INFO  +++ ConstructSD> Install Geant4 control directory:/ddg4/ConstructSD/
Geant4Kernel     INFO  ++ Registered global action RunInit of type dd4hep::sim::Test::Geant4TestRunAction
Geant4UI         INFO  +++ RunAction> Install Geant4 control directory:/ddg4/RunAction/
DDSim.Helper.Random INFO Choosing random seed for you: 631465373
Geant4UI         INFO  +++ LcioOutput> Install Geant4 control directory:/ddg4/LcioOutput/
Geant4UI         INFO  +++ EventAction> Install Geant4 control directory:/ddg4/EventAction/
DDSim            INFO ++++ Adding Geant4 Particle Gun ++++
DDSim            INFO Disabling the PrimaryHandler
Geant4UI         INFO  +++ GeneratorAction> Install Geant4 control directory:/ddg4/GeneratorAction/
Geant4UI         INFO  +++ Gun> Install Geant4 control directory:/ddg4/Gun/
Geant4UI         INFO  +++ InteractionMerger> Install Geant4 control directory:/ddg4/InteractionMerger/
Geant4UI         INFO  +++ TrackingAction> Install Geant4 control directory:/ddg4/TrackingAction/
Geant4UI         INFO  +++ SteppingAction> Install Geant4 control directory:/ddg4/SteppingAction/
Geant4UI         INFO  +++ ParticleHandler> Install Geant4 control directory:/ddg4/ParticleHandler/
DDSim.Helper.ParticleHandler INFO  *** definition of tracker region *** 
DDSim.Helper.ParticleHandler INFO     tracker_region_zmin = -10000
DDSim.Helper.ParticleHandler INFO     tracker_region_zmax = 10000
DDSim.Helper.ParticleHandler INFO     tracker_region_rmax = 5000
DDSim.Helper.ParticleHandler INFO  ************************************ 
DDSim.Helper.Filter INFO ReqFilt {'edep0', 'edep1kev'}
Geant4UI         INFO  +++ PhysicsList> Install Geant4 control directory:/ddg4/PhysicsList/
PhysicsList            +++ Dump of physics list component(s)
PhysicsList            +++ Extension name       FTFP_BERT
PhysicsList            +++ Transportation flag: 0
PhysicsList            +++ Program decays:      0
PhysicsList            +++ RangeCut:            0.700000
Geant4UI         INFO  +++ GlobalRangeCut> Install Geant4 control directory:/ddg4/GlobalRangeCut/
FieldSetup       INFO  Geant4 magnetic field tracking configured.
FieldSetup       INFO  G4MagIntegratorStepper:ClassicalRK4 G4Mag_EqRhs:Mag_UsualEqRhs
FieldSetup       INFO  Epsilon:[min:0.000050 mm max:0.001000 mm]
FieldSetup       INFO  Delta:[chord:0.250000 1-step:0.010000 intersect:0.001000] LargestStep 10000.000000 mm


          ################################
          !!! G4Backtrace is activated !!!
          ################################


**************************************************************
 Geant4 version Name: geant4-11-00-beta-01 [MT]   (25-June-2021)
                       Copyright : Geant4 Collaboration
                      References : NIM A 506 (2003), 250-303
                                 : IEEE-TNS 53 (2006), 270-278
                                 : NIM A 835 (2016), 186-225
                             WWW : http://geant4.org/
**************************************************************

Geant4UI         INFO  +++ Geant4RunManager> Install Geant4 control directory:/ddg4/Geant4RunManager/
Geant4RunManager WARN  +++ Configured run manager of type: G4RunManager.
Geant4Kernel     WARN  +++ Multi-threaded mode requested, but not supported by this compilation of Geant4.
Geant4Kernel     WARN  +++ Falling back to single threaded mode.
Geant4Exec       WARN  +++ Only 0 subdetectors present. You sure you loaded the geometry properly?
G4PhysListFactory::GetReferencePhysList <FTFP_BERT>  EMoption= 0
<<< Geant4 Physics List simulation engine: FTFP_BERT


Myphysics        INFO  Registered Geant4 physics constructor G4StepLimiterPhysics to physics list
UserInitialization INFO  +++ Executing Geant4UserActionInitialization::Build. Context:0x555577fc5860 Kernel:0x5555778ac460 [-1]
UI               INFO  +++ Starting G4UIExecutive 'DDG4' of type tcsh....
Available UI session types: [ GAG, tcsh, csh ]

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (109.753,41.2677,1068.83)
P1 = (109.803,40.6582,1069.38)
P2 = (109.803,40.6582,1069.38)
Side1 length (P0->P1) = 0.824599
Side2 length (P1->P2) = 0
Side3 length (P2->P0) = 0.824599
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (109.753,41.2677,1068.83)
P[1]      = (109.803,40.6582,1069.38)
P[2]      = (109.803,40.6582,1069.38)
*********************************************************************

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (57.5019,-6.55271,1069.14)
P1 = (57.5019,-6.55271,1069.14)
P2 = (61.5019,-6.5531,1069.14)
Side1 length (P0->P1) = 0
Side2 length (P1->P2) = 4
Side3 length (P2->P0) = 4
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (57.5019,-6.55271,1069.14)
P[1]      = (57.5019,-6.55271,1069.14)
P[2]      = (61.5019,-6.5531,1069.14)
*********************************************************************

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (101.713,40.6004,1070.51)
P1 = (109.803,40.6582,1069.38)
P2 = (109.803,40.6582,1069.38)
Side1 length (P0->P1) = 8.16956
Side2 length (P1->P2) = 0
Side3 length (P2->P0) = 8.16956
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (101.713,40.6004,1070.51)
P[1]      = (109.803,40.6582,1069.38)
P[2]      = (109.803,40.6582,1069.38)
*********************************************************************

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (57.5019,-6.55271,1069.14)
P1 = (57.5019,-1.98234,1066.39)
P2 = (57.5019,-6.55271,1069.14)
Side1 length (P0->P1) = 5.33397
Side2 length (P1->P2) = 5.33397
Side3 length (P2->P0) = 0
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (57.5019,-6.55271,1069.14)
P[1]      = (57.5019,-1.98234,1066.39)
P[2]      = (57.5019,-6.55271,1069.14)
*********************************************************************

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (565.4,40.6582,1069.38)
P1 = (565.4,40.6582,1069.38)
P2 = (565.45,41.2677,1068.83)
Side1 length (P0->P1) = 0
Side2 length (P1->P2) = 0.824597
Side3 length (P2->P0) = 0.824597
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (565.4,40.6582,1069.38)
P[1]      = (565.4,40.6582,1069.38)
P[2]      = (565.45,41.2677,1068.83)
*********************************************************************

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (617.702,-6.5531,1069.14)
P1 = (613.702,-6.5531,1069.14)
P2 = (617.702,-6.5531,1069.14)
Side1 length (P0->P1) = 3.99998
Side2 length (P1->P2) = 3.99998
Side3 length (P2->P0) = 0
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (617.702,-6.5531,1069.14)
P[1]      = (613.702,-6.5531,1069.14)
P[2]      = (617.702,-6.5531,1069.14)
*********************************************************************

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (573.491,40.6004,1070.51)
P1 = (565.4,40.6582,1069.38)
P2 = (565.4,40.6582,1069.38)
Side1 length (P0->P1) = 8.16952
Side2 length (P1->P2) = 0
Side3 length (P2->P0) = 8.16952
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (573.491,40.6004,1070.51)
P[1]      = (565.4,40.6582,1069.38)
P[2]      = (565.4,40.6582,1069.38)
*********************************************************************

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1001
      issued by : G4TriangularFacet::G4TriangularFacet()
Facet is too small or too narrow.
Triangle area = 0
P0 = (617.702,-3.7916,1067.82)
P1 = (617.702,-6.5531,1069.14)
P2 = (617.702,-6.5531,1069.14)
Side1 length (P0->P1) = 3.05857
Side2 length (P1->P2) = 0
Side3 length (P2->P0) = 3.05857
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4TessellatedSolid::AddFacet()
Attempt to add facet not properly defined.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


*********************************************************************
FACET TYPE       = G4TriangularFacet
ABSOLUTE VECTORS = 
P[0]      = (617.702,-3.7916,1067.82)
P[1]      = (617.702,-6.5531,1069.14)
P[2]      = (617.702,-6.5531,1069.14)
*********************************************************************
Geant4Converter  INFO  +++  Successfully converted geometry to Geant4.
Geant4UI         INFO  +++ ConstructGeo> Install Geant4 control directory:/ddg4/ConstructGeo/

 hInelastic FTFP_BERT : threshold between BERT and FTFP is over the interval 
 for pions :   3 to 6 GeV
 for kaons :   3 to 6 GeV
 for proton :  3 to 6 GeV
 for neutron : 3 to 6 GeV

### Adding tracking cuts for neutron  TimeCut(ns)= 10000  KinEnergyCut(MeV)= 0
GlobalRangeCut   INFO  +++ Set PhysicsList RangeCut to 0.700000 mm
Idle> /gun/particle geantino
Idle> /tracking/verbose 1
Idle>  /gun/number 10    
Idle> /run/beamOn

phot:  for gamma SubType=12 BuildTable=0
      LambdaPrime table from 200 keV to 100 TeV in 61 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
 LivermorePhElectric : Emin=    0 meV Emax=  100 TeV  SauterGavrila Fluo

compt:  for gamma SubType=13 BuildTable=1
      Lambda table from 100 eV  to 1 MeV, 7 bins/decade, spline: 1
      LambdaPrime table from 1 MeV to 100 TeV in 56 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
       Klein-Nishina : Emin=    0 meV Emax=  100 TeV

conv:  for gamma SubType=14 BuildTable=1
      Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     BetheHeitlerLPM : Emin=    0 meV Emax=  100 TeV  ModifiedTsai

Rayl:  for gamma SubType=11 BuildTable=1
      Lambda table from 100 eV  to 100 keV, 7 bins/decade, spline: 0
      LambdaPrime table from 100 keV to 100 TeV in 63 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LivermoreRayleigh : Emin=    0 meV Emax=  100 TeV  CullenGenerator

msc:  for e-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 meV Emax=  100 MeV Nbins=42 100 eV  - 100 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
        WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=42 100 MeV - 100 TeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

eIoni:  for e-  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        MollerBhabha : Emin=    0 meV Emax=  100 TeV

eBrem:  for e-  XStype:4  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    0 meV Emax=    1 GeV  ModifiedTsai
            eBremLPM : Emin=    1 GeV Emax=  100 TeV  ModifiedTsai

CoulombScat:  for e- XStype:3 SubType=1 BuildTable=1
      Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 1
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=  100 MeV Emax=  100 TeV

msc:  for e+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 meV Emax=  100 MeV Nbins=42 100 eV  - 100 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
        WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=42 100 MeV - 100 TeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

eIoni:  for e+  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        MollerBhabha : Emin=    0 meV Emax=  100 TeV

eBrem:  for e+  XStype:4  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    0 meV Emax=    1 GeV  ModifiedTsai
            eBremLPM : Emin=    1 GeV Emax=  100 TeV  ModifiedTsai

annihil:  for e+ XStype:2 SubType=5 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            eplus2gg : Emin=    0 meV Emax=  100 TeV

CoulombScat:  for e+ XStype:3 SubType=1 BuildTable=1
      Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 1
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=  100 MeV Emax=  100 TeV

msc:  for proton  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for proton  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 meV Emax=    2 MeV
          BetheBloch : Emin=    2 MeV Emax=  100 TeV

hBrems:  for proton  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

hPairProd:  for proton  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 17x1001; from 7.50618 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for proton XStype:3 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 1
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

msc:  for GenericIon  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 meV Emax=  100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

ionIoni:  for GenericIon  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.02
      Stopping Power data for 17 ion/material pairs
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            BraggIon : Emin=    0 meV Emax=    2 MeV
          BetheBloch : Emin=    2 MeV Emax=  100 TeV

msc:  for alpha  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 meV Emax=  100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

ionIoni:  for alpha  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.02
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            BraggIon : Emin=    0 meV Emax=7.9452 MeV
          BetheBloch : Emin=7.9452 MeV Emax=  100 TeV

msc:  for anti_proton  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for anti_proton  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 meV Emax=    2 MeV
          BetheBloch : Emin=    2 MeV Emax=  100 TeV

hBrems:  for anti_proton  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

hPairProd:  for anti_proton  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 17x1001; from 7.50618 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for anti_proton XStype:3 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 1
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

msc:  for kaon+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for kaon+  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 meV Emax=1.05231 MeV
          BetheBloch : Emin=1.05231 MeV Emax=  100 TeV

hBrems:  for kaon+  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

hPairProd:  for kaon+  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 18x1001; from 3.94942 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for kaon+ XStype:3 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 1
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

msc:  for kaon-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for kaon-  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 meV Emax=1.05231 MeV
          BetheBloch : Emin=1.05231 MeV Emax=  100 TeV

hBrems:  for kaon-  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

hPairProd:  for kaon-  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 18x1001; from 3.94942 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for kaon- XStype:3 SubType=1 BuildTable=1
      Used Lambda table of kaon+
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

msc:  for mu+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

muIoni:  for mu+  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 meV Emax=  200 keV
          BetheBloch : Emin=  200 keV Emax=    1 GeV
        MuBetheBloch : Emin=    1 GeV Emax=  100 TeV

muBrems:  for mu+  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
              MuBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

muPairProd:  for mu+  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 21x1001; from 1 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          muPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for mu+ XStype:3 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 1
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

msc:  for mu-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

muIoni:  for mu-  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 meV Emax=  200 keV
          BetheBloch : Emin=  200 keV Emax=    1 GeV
        MuBetheBloch : Emin=    1 GeV Emax=  100 TeV

muBrems:  for mu-  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
              MuBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

muPairProd:  for mu-  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 21x1001; from 1 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          muPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for mu- XStype:3 SubType=1 BuildTable=1
      Used Lambda table of mu+
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

msc:  for pi+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for pi+  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 meV Emax=297.505 keV
          BetheBloch : Emin=297.505 keV Emax=  100 TeV

hBrems:  for pi+  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

hPairProd:  for pi+  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 20x1001; from 1.11656 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for pi+ XStype:3 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 1
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

msc:  for pi-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 meV Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for pi-  XStype:1  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 meV Emax=297.505 keV
          BetheBloch : Emin=297.505 keV Emax=  100 TeV

hBrems:  for pi-  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

hPairProd:  for pi-  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 20x1001; from 1.11656 GeV to 100 TeV 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 meV Emax=  100 TeV  ModifiedMephi

CoulombScat:  for pi- XStype:3 SubType=1 BuildTable=1
      Used Lambda table of pi+
      ThetaMin(p) < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 meV Emax=  100 TeV

====================================================================
                  HADRONIC PROCESSES SUMMARY (verbose level 1)

---------------------------------------------------
                           Hadronic Processes for neutron

  Process: hadElastic
        Model:             hElasticCHIPS: 0 meV ---> 100 TeV
     Cr_sctns:        G4NeutronElasticXS: 0 meV ---> 100 TeV

  Process: neutronInelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:      G4NeutronInelasticXS: 0 meV ---> 100 TeV

  Process: nCapture
        Model:               nRadCapture: 0 meV ---> 100 TeV
     Cr_sctns:        G4NeutronCaptureXS: 0 meV ---> 100 TeV

  Process: nKiller

---------------------------------------------------
                           Hadronic Processes for B-

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: B-Inelastic
        Model:                      FTFP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for D-

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: D-Inelastic
        Model:                      FTFP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for GenericIon

  Process: ionInelastic
        Model:  Binary Light Ion Cascade: 0 meV/n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

---------------------------------------------------
                           Hadronic Processes for He3

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

  Process: He3Inelastic
        Model:  Binary Light Ion Cascade: 0 meV/n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

---------------------------------------------------
                           Hadronic Processes for alpha

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

  Process: alphaInelastic
        Model:  Binary Light Ion Cascade: 0 meV/n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

---------------------------------------------------
                           Hadronic Processes for anti_He3

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: anti_He3Inelastic
        Model:                      FTFP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: hFritiofCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for anti_alpha

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: anti_alphaInelastic
        Model:                      FTFP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: hFritiofCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for anti_deuteron

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: anti_deuteronInelastic
        Model:                      FTFP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: hFritiofCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for anti_lambda

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: anti_lambdaInelastic
        Model:                      FTFP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: hFritiofCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for anti_neutron

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100.1 MeV
        Model:              AntiAElastic: 100 MeV ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: anti_neutronInelastic
        Model:                      FTFP: 0 meV ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: hFritiofCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for anti_proton

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100.1 MeV
        Model:              AntiAElastic: 100 MeV ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: anti_protonInelastic
        Model:                      FTFP: 0 meV ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: hFritiofCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for anti_triton

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: anti_tritonInelastic
        Model:                      FTFP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 meV ---> 25.6 PeV

  Process: hFritiofCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for deuteron

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

  Process: dInelastic
        Model:  Binary Light Ion Cascade: 0 meV/n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

---------------------------------------------------
                           Hadronic Processes for e+

  Process: electronNuclear
        Model:   G4ElectroVDNuclearModel: 0 meV ---> 1 PeV
     Cr_sctns:          ElectroNuclearXS: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for e-

  Process: electronNuclear
        Model:   G4ElectroVDNuclearModel: 0 meV ---> 1 PeV
     Cr_sctns:          ElectroNuclearXS: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for gamma

  Process: photonNuclear
        Model:               GammaNPreco: 0 meV ---> 200 MeV
        Model:            BertiniCascade: 199 MeV ---> 6 GeV
        Model:           TheoFSGenerator: 3 GeV ---> 100 TeV
     Cr_sctns:            PhotoNuclearXS: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for kaon+

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: kaon+Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for kaon-

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: kaon-Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: hBertiniCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for lambda

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: lambdaInelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for mu+

  Process: muonNuclear
        Model:      G4MuonVDNuclearModel: 0 meV ---> 1 PeV
     Cr_sctns:     KokoulinMuonNuclearXS: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for mu-

  Process: muonNuclear
        Model:      G4MuonVDNuclearModel: 0 meV ---> 1 PeV
     Cr_sctns:     KokoulinMuonNuclearXS: 0 meV ---> 100 TeV

  Process: muMinusCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for pi+

  Process: hadElastic
        Model:           hElasticGlauber: 0 meV ---> 100 TeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 meV ---> 100 TeV

  Process: pi+Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for pi-

  Process: hadElastic
        Model:           hElasticGlauber: 0 meV ---> 100 TeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 meV ---> 100 TeV

  Process: pi-Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 meV ---> 100 TeV

  Process: hBertiniCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for proton

  Process: hadElastic
        Model:             hElasticCHIPS: 0 meV ---> 100 TeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 meV ---> 100 TeV

  Process: protonInelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 meV ---> 100 TeV

---------------------------------------------------
                           Hadronic Processes for sigma-

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: sigma-Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 meV ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 meV ---> 100 TeV

  Process: hBertiniCaptureAtRest

---------------------------------------------------
                           Hadronic Processes for triton

  Process: hadElastic
        Model:              hElasticLHEP: 0 meV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

  Process: tInelastic
        Model:  Binary Light Ion Cascade: 0 meV/n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 meV ---> 25.6 PeV

================================================================
=======================================================================
======       Pre-compound/De-excitation Physics Parameters     ========
=======================================================================
Type of pre-compound inverse x-section              3
Pre-compound model active                           1
Pre-compound excitation low energy (MeV)            0.1
Pre-compound excitation high energy (MeV)           30
Type of de-excitation inverse x-section             3
Type of de-excitation factory                       Evaporation+GEM
Number of de-excitation channels                    68
Min excitation energy (keV)                         0.01
Min energy per nucleon for multifragmentation (MeV) 2e+05
Limit excitation energy for Fermi BreakUp (MeV)     20
Level density (1/MeV)                               0.075
Use simple level density model                      1
Use discrete excitation energy of the residual      0
Time limit for long lived isomeres (ns)             1
Isomer production flag                              1
Internal e- conversion flag                         1
Store e- internal conversion data                   0
Electron internal conversion ID                     3
Correlated gamma emission flag                      0
Max 2J for sampling of angular correlations         10
=======================================================================
GenerationInit   INFO  +++ Initializing event 1. Within run:0 event 1.
ParticleHandler  INFO  +++ Event 0 Begin event action. Access event related information.

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 10,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 9,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 8,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 7,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 6,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 5,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 4,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 3,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 2,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = geantino,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     1e+03        0        0         0 world_volume_1 initStep
    1    3e+03        0        0     1e+03        0    3e+03     3e+03  OutOfWorld Transportation
Geant4OutputAction WARN  +++ The value of G4HCofThisEvent is NULL. No collections saved!
GenerationInit   WARN  +++ Finished run 0 after 1 events (1 events in total)
Idle> 

Geant4Kernel     INFO  ++ Terminate Geant4 and delete associated actions.
DDSim            INFO DDSim            INFO  Total Time:   27.05 s (User), 9.24 s (System)
[Thread 0x7fffd4ae2640 (LWP 96036) exited]
[Inferior 1 (process 96000) exited normally]
(gdb) q
12:27:07-frankm~/SW/DD4hep_head_dbg.root_vmaster.g4_11.00.b01_MT.G4_units/build/checkout (master)$

MarkusFrankATcernch · 2023-06-23T10:40:38Z

This is how the loaded geometry looks like:

Is however slightly imperfect. If moved, it crashes ROOT.

MarkusFrankATcernch · 2023-06-23T10:48:14Z

Is this now what one expects?
Or do I do something entirely wrong?

github-actions · 2023-06-23T12:21:07Z

Test Results

      6 files       6 suites 6h 13m 54s ⏱️
  356 tests   355 ✔️ 0 💤 1 ❌
1 052 runs 1 049 ✔️ 0 💤 3 ❌

For more details on these failures, see this check.

Results for commit 04d8cdf.

♻️ This comment has been updated with latest results.

MarkusFrankATcernch · 2023-06-23T14:57:26Z

Result of the overlap check:

MarkusFrankATcernch mentioned this pull request Jun 23, 2023

Abnormal material budget of DDCAD-imported files #1134

Closed

MarkusFrankATcernch added 6 commits June 23, 2023 17:44

Add Armin's example from issue 1134

9d7d880

Add example readme

3b7d1ba

Add example readme

cd677de

Add example readme

d4cc7af

Add example readme

3aeaf10

Add examples as tests for issue 1134

bd5c461

MarkusFrankATcernch force-pushed the master branch from cdcf69c to bd5c461 Compare June 23, 2023 15:44

MarkusFrankATcernch mentioned this pull request Jun 23, 2023

How to check for overlaps when using DDCAD #1121

Closed

MarkusFrankATcernch enabled auto-merge (rebase) June 23, 2023 15:59

Add examples as tests for issue 1134

04d8cdf

MarkusFrankATcernch changed the title ~~WIP: Test for issue 1134~~ Test for issue 1134 Jun 23, 2023

Fix regex of test DDCAD_Issue1134_g4overlap_check

f54f541

MarkusFrankATcernch merged commit 21df99a into AIDASoft:master Jun 24, 2023

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Test for issue 1134 #1139

Test for issue 1134 #1139

MarkusFrankATcernch commented Jun 23, 2023 •

edited by andresailer

Loading

MarkusFrankATcernch commented Jun 23, 2023

MarkusFrankATcernch commented Jun 23, 2023

MarkusFrankATcernch commented Jun 23, 2023

github-actions bot commented Jun 23, 2023 •

edited

Loading

MarkusFrankATcernch commented Jun 23, 2023

Test for issue 1134 #1139

Test for issue 1134 #1139

Conversation

MarkusFrankATcernch commented Jun 23, 2023 • edited by andresailer Loading

MarkusFrankATcernch commented Jun 23, 2023

MarkusFrankATcernch commented Jun 23, 2023

MarkusFrankATcernch commented Jun 23, 2023

github-actions bot commented Jun 23, 2023 • edited Loading

Test Results

MarkusFrankATcernch commented Jun 23, 2023

MarkusFrankATcernch commented Jun 23, 2023 •

edited by andresailer

Loading

github-actions bot commented Jun 23, 2023 •

edited

Loading