]>
Commit | Line | Data |
---|---|---|
1 | // $Id$ | |
2 | // Category: physics | |
3 | // | |
4 | // See the class description in the header file. | |
5 | ||
6 | #include "TG4G3Defaults.h" | |
7 | #include "TG4Globals.h" | |
8 | ||
9 | #include <math.h> | |
10 | ||
11 | // static const data members | |
12 | ||
13 | // precision tolerance | |
14 | const G4double TG4G3Defaults::fgkCutTolerance = 1. * keV; | |
15 | ||
16 | // kinetic energy cuts | |
17 | const G4double TG4G3Defaults::fgkCUTGAM = 0.001 * GeV; | |
18 | const G4double TG4G3Defaults::fgkCUTELE = 0.001 * GeV; | |
19 | const G4double TG4G3Defaults::fgkCUTNEU = 0.01 * GeV; | |
20 | const G4double TG4G3Defaults::fgkCUTHAD = 0.01 * GeV; | |
21 | const G4double TG4G3Defaults::fgkCUTMUO = 0.01 * GeV; | |
22 | const G4double TG4G3Defaults::fgkBCUTE = fgkCUTGAM; | |
23 | const G4double TG4G3Defaults::fgkBCUTM = fgkCUTGAM; | |
24 | const G4double TG4G3Defaults::fgkDCUTE = 10. * TeV; | |
25 | const G4double TG4G3Defaults::fgkDCUTM = 10. * TeV; | |
26 | const G4double TG4G3Defaults::fgkPPCUTM = 0.01 * GeV; | |
27 | ||
28 | // physics processes | |
29 | const TG3FlagValue TG4G3Defaults::fgkPAIR = kActivate; // 1 | |
30 | const TG3FlagValue TG4G3Defaults::fgkCOMP = kActivate; // 1 | |
31 | const TG3FlagValue TG4G3Defaults::fgkPHOT = kActivate; // 1 | |
32 | const TG3FlagValue TG4G3Defaults::fgkPFIS = kInActivate; // 0 | |
33 | const TG3FlagValue TG4G3Defaults::fgkDRAY = kActivate2; // 2 | |
34 | const TG3FlagValue TG4G3Defaults::fgkANNI = kActivate; // 1 | |
35 | const TG3FlagValue TG4G3Defaults::fgkBREM = kActivate; // 1 | |
36 | const TG3FlagValue TG4G3Defaults::fgkHADR = kActivate; // 1 | |
37 | const TG3FlagValue TG4G3Defaults::fgkMUNU = kInActivate; // 0 | |
38 | const TG3FlagValue TG4G3Defaults::fgkDCAY = kActivate; // 1 | |
39 | const TG3FlagValue TG4G3Defaults::fgkLOSS = kActivate2; // 2 | |
40 | const TG3FlagValue TG4G3Defaults::fgkMULS = kActivate; // 1 | |
41 | ||
42 | TG4G3Defaults::TG4G3Defaults() { | |
43 | // | |
44 | } | |
45 | ||
46 | TG4G3Defaults::~TG4G3Defaults() { | |
47 | // | |
48 | } | |
49 | ||
50 | G4double TG4G3Defaults::CutValue(G4int g3Cut) | |
51 | { | |
52 | // Returns the G3 default value for the specified cut. | |
53 | // --- | |
54 | ||
55 | switch (g3Cut) { | |
56 | case kCUTGAM: | |
57 | return fgkCUTGAM; | |
58 | case kCUTELE: | |
59 | return fgkCUTELE; | |
60 | case kCUTNEU: | |
61 | return fgkCUTNEU; | |
62 | case kCUTHAD: | |
63 | return fgkCUTHAD; | |
64 | case kCUTMUO: | |
65 | return fgkCUTMUO; | |
66 | case kBCUTE: | |
67 | return fgkBCUTE; | |
68 | case kBCUTM: | |
69 | return fgkBCUTM; | |
70 | case kDCUTE: | |
71 | return fgkDCUTE; | |
72 | case kDCUTM: | |
73 | return fgkDCUTM; | |
74 | case kPPCUTM: | |
75 | return fgkPPCUTM; | |
76 | default: | |
77 | TG4Globals::Warning("TG4G3Defaults::CutValue: Inconsistent cut."); | |
78 | return 0.; | |
79 | } | |
80 | } | |
81 | ||
82 | TG3FlagValue TG4G3Defaults::FlagValue(G4int g3Flag) | |
83 | { | |
84 | // Returns the G3 default value for the specified flag. | |
85 | // --- | |
86 | ||
87 | switch (g3Flag) { | |
88 | case kPAIR: | |
89 | return fgkPAIR; | |
90 | case kCOMP: | |
91 | return fgkCOMP; | |
92 | case kPHOT: | |
93 | return fgkPHOT; | |
94 | case kPFIS: | |
95 | return fgkPFIS; | |
96 | case kDRAY: | |
97 | return fgkDRAY; | |
98 | case kANNI: | |
99 | return fgkANNI; | |
100 | case kBREM: | |
101 | return fgkBREM; | |
102 | case kHADR: | |
103 | return fgkHADR; | |
104 | case kMUNU: | |
105 | return fgkMUNU; | |
106 | case kDCAY: | |
107 | return fgkDCAY; | |
108 | case kLOSS: | |
109 | return fgkLOSS; | |
110 | case kMULS: | |
111 | return fgkMULS; | |
112 | default: | |
113 | TG4Globals::Warning("TG4G3Defaults::FlagValue: Inconsistent flag."); | |
114 | return kUnset; | |
115 | } | |
116 | } | |
117 | ||
118 | G4bool TG4G3Defaults::IsDefaultCut(TG3Cut g3Cut, G4double value) | |
119 | { | |
120 | // Tests if the parameter value is equal to the G3 default value. | |
121 | // --- | |
122 | ||
123 | if (abs(value*GeV - CutValue(g3Cut)) > fgkCutTolerance) | |
124 | return false; | |
125 | else | |
126 | return true; | |
127 | } | |
128 | ||
129 | G4bool TG4G3Defaults::IsDefaultFlag(TG3Flag g3Flag, G4double value) | |
130 | { | |
131 | // Tests if the parameter value is equal to the G3 default value. | |
132 | // --- | |
133 | ||
134 | if (value == FlagValue(g3Flag)) | |
135 | return true; | |
136 | else | |
137 | return false; | |
138 | } |