GeantChecklist

From Mu2eWiki
Jump to navigation Jump to search

Introduction

So much of the understanding of the performance of mu2e depends on simulation studies that a change in geant version could have far-reaching effects. This page contains a checklist for tests to perform before moving to a new geant. Typically, some tests would be skipped because the expert decides the effort to validate would not be justified for the low level of risk.

Note on MT build

4/4/18

  • The geant MT branch code can be compiled against (and run using) a sequential version of G4 without the '-DG4MULTITHREADED' flag in the SConscript file.
  • The MT branch code code can be compiled against an MT version of G4 with the addition of the '-DG4MULTITHREADED' flag in the SConscript file.
  • MT version can be run in 3 modes: (1) in sequential mode, (2) in MT mode using a single thread in which there is a master thread that controls a single worker thread, and (3) in MT mode.


v4_10_4

Started Mar 1, 2018, proposed update to v4_10_4 from v4_10_2_p03e. 3/21/18 - effort paused - now working on v4_10_4 EMZ. On 4/4/18, the decision was taken to restart this effort and delay EMZ, which is more complicated.

To make a release with this change, edit setup.sh and change this line

setup -B geant4 v4_10_2_p03e -q${MU2E_UPS_QUALIFIERS}${MU2E_G4_GRAPHICS_QUALIFIER}

to

setup -B geant4 v4_10_4 -qcl23:${MU2E_UPS_QUALIFIERS}${MU2E_G4_GRAPHICS_QUALIFIER}


How to change the stepper to the new one requiring fewer field evaluations. This is off by default and is not part of this checklist.

physics.producers.g4run.physics.stepper : "G4DormandPrince745"

clhep 2.4.0.0 (not used in this checklist evaluation) will enable the switch from HepJamesRandom to MixMax random number engine (art sets the default)


Results from tests, 1M conversion electrons, 300K beam stage1

Sample CPU (s) MEM (GB) file size (MB)
old beam 9351 1.611 6.290
new beam 8483 (-9%) 1.715 (+6%) 6.397(+2%)
old CE 6795 1.742 918.396
new CE 6180 (-9%) 1.804 (+4%) 909.979(-1%)


beam comparison
C.E. comparison
stopped muons

Sign off Responsible Topic Notes
3/14/18 Krzysztof release notes presented available
3/14/18 Krzysztof supporting code (new lists, options) changes enabling new features done in February and March
3/14/18 Ray high-stats conversion electrons comparison see links above and talk
3/14/18 Ray high-stats beam stage 1 comparison see links above and talk
3/14/18 Ray/Krzysztof CPU time and product sizes see above and talk
3/14/18 Ray/Andrei stopped muon rates and shape see links above and talk
4/5/18 Rob magnetic field transport the checks done so far are adequate
4/16/18 Dave (LBL) tracker charge deposition/hit resolution roughly 10% more tracker hits than in v4_10_2
4/16/18 Dave (LBL) tracker background hit rates
4/16/18 Dave (LBL) tracker integrated charge
4/16/18 Dave (LBL) tracker electronics radiation A full radiation dose test takes weeks to perform. I validated that the tracker photon interactions are reasonable
4/16/18 Bertrand calorimeter resolution (email) I just looked at the various plots produced for v4_10_4, and they all look fine for the calorimeter.
4/16/18 Bertrand calorimeter electronics radiation
Ralf CRV light model
Yuri neutron production rates and transport
Yuri CRV deadtime
4/16/18 Dave (LBL) conversion electron efficiency Efficiency and quality are slightly degraded
4/14/18 Jim/Yaqian STM efficiency lead flourescence and brem bkg checked
4/14/18 Jim/Yaqian STM background lead flourescence and brem bkg checked
4/10/18 Andrei ExtMon efficiency Skip this time (the baseline has not been checked in years). This is feasible to do in the future.
4/10/18 Andrei ExtMon background Skip. A complete study would require very large amount of CPU and effort. Unlikely to be feasible for GEANT upgrade checks. A simpler proxy like e.g. neutron flux in the ExtMon room might be feasible.
Ralf cosmic background
Bob p-bar background
4/10/18 Andrei beam electron background Skip. May be feasible in the future.
4/19/18 Dave (Lou) Pion backgrounds It would require a dedicated run to do a "real" test. However, if I look at muon stops as a proxy for pion stops, then I see nothing here that makes me panic.

v4_10_4 EMZ

Started Mar 21, 2018, proposed update to v4_10_4 EMZ from v4_10_2_p03e. At the code/soft meeting on 4/4/18, we decided to delay moving to EMZ since the very different file sizes caused complications with no obvious benefit. This will be considered at a later date, along with regional minRangeCut's.

How to trigger Mu2e specific electromagnetic option4 (_EMZ) with modified msc model transition energy (committed 3/16/18)

physics.producers.g4run.debug.diagLevel         : 1
# the above is to see the ranges printed
physics.producers.g4run.physics.modifyEMOption4 : true
physics.producers.g4run.physics.physicsListName : "ShieldingM_EMZ"

How to change the stepper to the new one requiring fewer field evaluations (off, except where noted):

physics.producers.g4run.physics.stepper : "G4DormandPrince745"

clhep 2.4.0.0 (not used yet) will enable the switch from HepJamesRandom to MixMax random number engine (art sets the default)


Results from tests, 1M conversion electrons, 300K beam stage1

Sample CPU (s) MEM (GB) file size (MB)
old beam 7325 1.76 630
new beam 7896(+8%) 1.79(+2%) 718(+13%)
old CE 6844 1.75 92088
new CE 7099(+4%) 1.84(+5%) 187456(+100%)


beam stage 1 comparison
CE comparison
stops

Why are conversion electron files so much larger:

arProductsizes:
old:
        Size     Size/Entry   Fraction  Data Product Name
   505449709        50545.0      0.548  mu2e::SimParticlemv_g4run__ceSimReco.
   240938209        24093.8      0.261  mu2e::StepPointMCs_g4run_calorimeter_ceSimReco.
    30774582         3077.5      0.033  mu2e::CaloShowerSteps_CaloShowerStepFromStepPt_calorimeter_ceSimReco.
    28011825         2801.2      0.030  mu2e::StepPointMCs_g4run_tracker_ceSimReco.

new:
         Size     Size/Entry   Fraction  Data Product Name
   514208625        51420.9      0.277  mu2e::SimParticlemv_g4run__ceSimReco.
   745787106        74578.7      0.402  mu2e::StepPointMCs_g4run_calorimeter_ceSimReco.
   236050707        23605.1      0.127  mu2e::CaloShowerSteps_CaloShowerStepFromStepPt_calorimeter_ceSimReco.
   181788056        18178.8      0.098  mu2e::StepPointMCs_g4run_tracker_ceSimReco.



Sign off Responsible Topic Notes
3/21/18 Krzysztof release notes presented talk
3/14/18 Krzysztof supporting code (new lists, options) changes enabling new features done in February and March
Ray high-stats conversion electrons comparison see links above
Ray high-stats beam stage 1 comparison see links above
Ray/Krzysztof CPU time and product sizes see above
Ray/Andrei stopped muon rates and shape see link above
Rob magnetic field transport
Dave (LBL) tracker charge deposition/hit resolution talk
Dave (LBL) tracker background hit rates
Dave (LBL) tracker integrated charge
Dave (LBL) tracker electronics radiation
Bertrand calorimeter resolution 4/16/18 probably OK
Bertrand calorimeter electronics radiation 4/16/18 this needs more study
Ralf CRV light model
Yuri neutron production rates and transport
Yuri CRV deadtime
Dave (LBL) conversion electron efficiency
Jim STM efficiency
Jim STM background
Andrei ExtMon efficiency
Andrei ExtMon background
Ralf cosmic background
Bob p-bar background
Andrei beam electron background
Dave (Lou) Pion backgrounds