GeantChecklist

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 talk
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