Geant4  10.01.p02
G4CascadeParameters.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4CascadeParameters.cc 72016 2013-07-03 16:24:15Z mkelsey $
27 // Encapsulate all user-configurable parameters with associated envvars
28 //
29 // 20120912 M. Kelsey -- Add interface to support UI commands
30 // 20130304 M. Kelsey -- Add flag to collect and display cascade structure
31 // 20130308 M. Kelsey -- Add flag to use separate 3-body momentum generators
32 // 20130421 M. Kelsey -- Add flag for CHECK_ECONS, replacing #ifdef's
33 // 20130702 M. Kelsey -- Add flag to use N-body phase-space generator
34 // 20140311 G. Cosmo -- Implement standard (non-const) singleton pattern
35 // 20140929 M. Kelsey -- Enable some parameters as default true (must be set
36 // '0' for false): PreCompound, phase-space, clustering,
37 // trailing effect
38 // 20141111 M. Kelsey -- Revert defaults for PreCompound, phase-space,
39 // and trailing effect.
40 // 20141121 Use G4AutoDelete to avoid end-of-thread memory leaks
41 
42 #include "G4CascadeParameters.hh"
44 #include "G4AutoDelete.hh"
45 #include <stdlib.h>
46 #include <iostream>
47 using std::endl;
48 
49 
50 // Singleton accessor
51 
53 
55  if (!fpInstance) {
58  }
59 
60  return fpInstance;
61 }
62 
63 
64 // Constructor initializes everything once
65 
66 #define OLD_RADIUS_UNITS (3.3836/1.2) // Used with NucModel params
67 
69  : G4CASCADE_VERBOSE(getenv("G4CASCADE_VERBOSE")),
70  G4CASCADE_CHECK_ECONS(getenv("G4CASCADE_CHECK_ECONS")),
71  G4CASCADE_USE_PRECOMPOUND(getenv("G4CASCADE_USE_PRECOMPOUND")),
72  G4CASCADE_DO_COALESCENCE(getenv("G4CASCADE_DO_COALESCENCE")),
73  G4CASCADE_PIN_ABSORPTION(getenv("G4CASCADE_PIN_ABSORPTION")),
74  G4CASCADE_SHOW_HISTORY(getenv("G4CASCADE_SHOW_HISTORY")),
75  G4CASCADE_USE_3BODYMOM(getenv("G4CASCADE_USE_3BODYMOM")),
76  G4CASCADE_USE_PHASESPACE(getenv("G4CASCADE_USE_PHASESPACE")),
77  G4CASCADE_RANDOM_FILE(getenv("G4CASCADE_RANDOM_FILE")),
78  G4NUCMODEL_USE_BEST(getenv("G4NUCMODEL_USE_BEST")),
79  G4NUCMODEL_RAD_2PAR(getenv("G4NUCMODEL_RAD_2PAR")),
80  G4NUCMODEL_RAD_SCALE(getenv("G4NUCMODEL_RAD_SCALE")),
81  G4NUCMODEL_RAD_SMALL(getenv("G4NUCMODEL_RAD_SMALL")),
82  G4NUCMODEL_RAD_ALPHA(getenv("G4NUCMODEL_RAD_ALPHA")),
83  G4NUCMODEL_RAD_TRAILING(getenv("G4NUCMODEL_RAD_TRAILING")),
84  G4NUCMODEL_FERMI_SCALE(getenv("G4NUCMODEL_FERMI_SCALE")),
85  G4NUCMODEL_XSEC_SCALE(getenv("G4NUCMODEL_XSEC_SCALE")),
86  G4NUCMODEL_GAMMAQD(getenv("G4NUCMODEL_GAMMAQD")),
87  DPMAX_2CLUSTER(getenv("DPMAX_2CLUSTER")),
88  DPMAX_3CLUSTER(getenv("DPMAX_3CLUSTER")),
89  DPMAX_4CLUSTER(getenv("DPMAX_4CLUSTER")),
90  messenger(0) {
92  Initialize();
93 }
94 
101  G4CASCADE_DO_COALESCENCE[0]!='0');
106  G4CASCADE_USE_PHASESPACE[0]!='0');
111  : (BEST_PAR?1.0:OLD_RADIUS_UNITS));
113  : (BEST_PAR?1.992:(8.0/OLD_RADIUS_UNITS))) * RADIUS_SCALE);
115  : (BEST_PAR?0.84:0.70));
117  : 0.) * RADIUS_SCALE);
119  : (BEST_PAR?0.685:(1.932/OLD_RADIUS_UNITS))) * RADIUS_SCALE);
121  : (BEST_PAR?0.1:1.0) );
123  DPMAX_DOUBLET = (DPMAX_2CLUSTER ? strtod(DPMAX_2CLUSTER,0) : 0.090);
124  DPMAX_TRIPLET = (DPMAX_3CLUSTER ? strtod(DPMAX_3CLUSTER,0) : 0.108);
125  DPMAX_ALPHA = (DPMAX_4CLUSTER ? strtod(DPMAX_4CLUSTER,0) : 0.115);
126 }
127 
129  delete messenger;
130 }
131 
132 
133 // Report any non-default parameters (used by G4CascadeInterface)
134 
135 void G4CascadeParameters::DumpConfig(std::ostream& os) const {
136  if (G4CASCADE_VERBOSE)
137  os << "G4CASCADE_VERBOSE = " << G4CASCADE_VERBOSE << endl;
139  os << "G4CASCADE_CHECK_ECONS = " << G4CASCADE_CHECK_ECONS << endl;
141  os << "G4CASCADE_USE_PRECOMPOUND = " << G4CASCADE_USE_PRECOMPOUND << endl;
143  os << "G4CASCADE_DO_COALESCENCE = " << G4CASCADE_DO_COALESCENCE << endl;
145  os << "G4CASCADE_PIN_ABSORPTION = " << G4CASCADE_PIN_ABSORPTION << endl;
147  os << "G4CASCADE_SHOW_HISTORY = " << G4CASCADE_SHOW_HISTORY << endl;
149  os << "G4CASCADE_USE_3BODYMOM = " << G4CASCADE_USE_3BODYMOM << endl;
151  os << "G4CASCADE_USE_PHASESPACE = " << G4CASCADE_USE_PHASESPACE << endl;
153  os << "G4CASCADE_RANDOM_FILE = " << G4CASCADE_RANDOM_FILE << endl;
155  os << "G4NUCMODEL_USE_BEST = " << G4NUCMODEL_USE_BEST << endl;
157  os << "G4NUCMODEL_RAD_2PAR = " << G4NUCMODEL_RAD_2PAR << endl;
159  os << "G4NUCMODEL_RAD_SCALE = " << G4NUCMODEL_RAD_SCALE << endl;
161  os << "G4NUCMODEL_RAD_SMALL = " << G4NUCMODEL_RAD_SMALL << endl;
163  os << "G4NUCMODEL_RAD_ALPHA = " << G4NUCMODEL_RAD_ALPHA << endl;
165  os << "G4NUCMODEL_RAD_TRAILING = " << G4NUCMODEL_RAD_TRAILING << endl;
167  os << "G4NUCMODEL_FERMI_SCALE = " << G4NUCMODEL_FERMI_SCALE << endl;
169  os << "G4NUCMODEL_XSEC_SCALE = " << G4NUCMODEL_XSEC_SCALE << endl;
170  if (G4NUCMODEL_GAMMAQD)
171  os << "G4NUCMODEL_GAMMAQD = " << G4NUCMODEL_GAMMAQD << endl;
172  if (DPMAX_2CLUSTER)
173  os << "DPMAX_2CLUSTER = " << DPMAX_2CLUSTER << endl;
174  if (DPMAX_3CLUSTER)
175  os << "DPMAX_3CLUSTER = " << DPMAX_3CLUSTER << endl;
176  if (DPMAX_4CLUSTER)
177  os << "DPMAX_4CLUSTER = " << DPMAX_4CLUSTER << endl;
178 }
const char * G4NUCMODEL_USE_BEST
void DumpConfig(std::ostream &os) const
const char * G4NUCMODEL_RAD_SMALL
const char * G4CASCADE_VERBOSE
#define OLD_RADIUS_UNITS
const char * G4CASCADE_SHOW_HISTORY
static G4CascadeParameters * fpInstance
const char * G4NUCMODEL_XSEC_SCALE
const char * G4CASCADE_DO_COALESCENCE
void Register(T *inst)
Definition: G4AutoDelete.hh:65
const char * G4NUCMODEL_RAD_TRAILING
const char * G4NUCMODEL_RAD_2PAR
const char * G4NUCMODEL_RAD_SCALE
const char * G4NUCMODEL_FERMI_SCALE
const char * G4CASCADE_PIN_ABSORPTION
const char * G4CASCADE_CHECK_ECONS
friend class G4CascadeParamMessenger
static const G4CascadeParameters * Instance()
const char * G4CASCADE_USE_3BODYMOM
const char * G4CASCADE_USE_PRECOMPOUND
const char * G4CASCADE_USE_PHASESPACE
const char * G4CASCADE_RANDOM_FILE
G4CascadeParamMessenger * messenger
const char * G4NUCMODEL_GAMMAQD
const char * G4NUCMODEL_RAD_ALPHA