4 // See the class description in the header file.
6 #include "TG4G3Defaults.h"
7 #include "TG4Globals.h"
11 // static const data members
13 // precision tolerance
14 const G4double TG4G3Defaults::fgCutTolerance = 1. * keV;
16 // kinetic energy cuts
17 const G4double TG4G3Defaults::fgCUTGAM = 0.001 * GeV;
18 const G4double TG4G3Defaults::fgCUTELE = 0.001 * GeV;
19 const G4double TG4G3Defaults::fgCUTNEU = 0.01 * GeV;
20 const G4double TG4G3Defaults::fgCUTHAD = 0.01 * GeV;
21 const G4double TG4G3Defaults::fgCUTMUO = 0.01 * GeV;
22 const G4double TG4G3Defaults::fgBCUTE = fgCUTGAM;
23 const G4double TG4G3Defaults::fgBCUTM = fgCUTGAM;
24 const G4double TG4G3Defaults::fgDCUTE = 10. * TeV;
25 const G4double TG4G3Defaults::fgDCUTM = 10. * TeV;
26 const G4double TG4G3Defaults::fgPPCUTM = 0.01 * GeV;
29 const TG3FlagValue TG4G3Defaults::fgPAIR = kActivate; // 1
30 const TG3FlagValue TG4G3Defaults::fgCOMP = kActivate; // 1
31 const TG3FlagValue TG4G3Defaults::fgPHOT = kActivate; // 1
32 const TG3FlagValue TG4G3Defaults::fgPFIS = kInActivate; // 0
33 const TG3FlagValue TG4G3Defaults::fgDRAY = kActivate2; // 2
34 const TG3FlagValue TG4G3Defaults::fgANNI = kActivate; // 1
35 const TG3FlagValue TG4G3Defaults::fgBREM = kActivate; // 1
36 const TG3FlagValue TG4G3Defaults::fgHADR = kActivate; // 1
37 const TG3FlagValue TG4G3Defaults::fgMUNU = kInActivate; // 0
38 const TG3FlagValue TG4G3Defaults::fgDCAY = kActivate; // 1
39 const TG3FlagValue TG4G3Defaults::fgLOSS = kActivate2; // 2
40 const TG3FlagValue TG4G3Defaults::fgMULS = kActivate; // 1
42 TG4G3Defaults::TG4G3Defaults() {
46 TG4G3Defaults::~TG4G3Defaults() {
50 G4double TG4G3Defaults::CutValue(G4int g3Cut)
52 // Returns the G3 default value for the specified cut.
57 return fgCUTGAM; break;
59 return fgCUTELE; break;
61 return fgCUTNEU; break;
63 return fgCUTHAD; break;
65 return fgCUTMUO; break;
67 return fgBCUTE; break;
69 return fgBCUTM; break;
71 return fgDCUTE; break;
73 return fgDCUTM; break;
75 return fgPPCUTM; break;
77 TG4Globals::Warning("TG4G3Defaults::CutValue: Inconsistent cut.");
82 TG3FlagValue TG4G3Defaults::FlagValue(G4int g3Flag)
84 // Returns the G3 default value for the specified flag.
101 return fgBREM; break;
103 return fgHADR; break;
105 return fgMUNU; break;
107 return fgDCAY; break;
109 return fgLOSS; break;
111 return fgMULS; break;
113 TG4Globals::Warning("TG4G3Defaults::FlagValue: Inconsistent flag.");
118 G4bool TG4G3Defaults::IsDefaultCut(TG3Cut g3Cut, G4double value)
120 // Tests if the parameter value is equal to the G3 default value.
123 if (abs(value*GeV - CutValue(g3Cut)) > fgCutTolerance)
129 G4bool TG4G3Defaults::IsDefaultFlag(TG3Flag g3Flag, G4double value)
131 // Tests if the parameter value is equal to the G3 default value.
134 if (value == FlagValue(g3Flag))