]>
Commit | Line | Data |
---|---|---|
2817d3e2 | 1 | // $Id$ |
73dd85b2 | 2 | // Category: global |
2817d3e2 | 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 | |
c63f260d | 14 | const G4double TG4G3Defaults::fgkCutTolerance = 1. * keV; |
2817d3e2 | 15 | |
16 | // kinetic energy cuts | |
c63f260d | 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; | |
2817d3e2 | 27 | |
28 | // physics processes | |
9187be37 | 29 | const TG4G3ControlValue TG4G3Defaults::fgkPAIR = kActivate; // 1 |
30 | const TG4G3ControlValue TG4G3Defaults::fgkCOMP = kActivate; // 1 | |
31 | const TG4G3ControlValue TG4G3Defaults::fgkPHOT = kActivate; // 1 | |
32 | const TG4G3ControlValue TG4G3Defaults::fgkPFIS = kInActivate; // 0 | |
33 | const TG4G3ControlValue TG4G3Defaults::fgkDRAY = kActivate2; // 2 | |
34 | const TG4G3ControlValue TG4G3Defaults::fgkANNI = kActivate; // 1 | |
35 | const TG4G3ControlValue TG4G3Defaults::fgkBREM = kActivate; // 1 | |
36 | const TG4G3ControlValue TG4G3Defaults::fgkHADR = kActivate; // 1 | |
37 | const TG4G3ControlValue TG4G3Defaults::fgkMUNU = kInActivate; // 0 | |
38 | const TG4G3ControlValue TG4G3Defaults::fgkDCAY = kActivate; // 1 | |
39 | const TG4G3ControlValue TG4G3Defaults::fgkLOSS = kActivate2; // 2 | |
40 | const TG4G3ControlValue TG4G3Defaults::fgkMULS = kActivate; // 1 | |
2817d3e2 | 41 | |
73dd85b2 | 42 | //_____________________________________________________________________________ |
2817d3e2 | 43 | TG4G3Defaults::TG4G3Defaults() { |
44 | // | |
45 | } | |
46 | ||
73dd85b2 | 47 | //_____________________________________________________________________________ |
2817d3e2 | 48 | TG4G3Defaults::~TG4G3Defaults() { |
49 | // | |
50 | } | |
51 | ||
73dd85b2 | 52 | //_____________________________________________________________________________ |
2817d3e2 | 53 | G4double TG4G3Defaults::CutValue(G4int g3Cut) |
54 | { | |
55 | // Returns the G3 default value for the specified cut. | |
56 | // --- | |
57 | ||
58 | switch (g3Cut) { | |
59 | case kCUTGAM: | |
4e994a92 | 60 | return fgkCUTGAM; |
2817d3e2 | 61 | case kCUTELE: |
4e994a92 | 62 | return fgkCUTELE; |
2817d3e2 | 63 | case kCUTNEU: |
4e994a92 | 64 | return fgkCUTNEU; |
2817d3e2 | 65 | case kCUTHAD: |
4e994a92 | 66 | return fgkCUTHAD; |
2817d3e2 | 67 | case kCUTMUO: |
4e994a92 | 68 | return fgkCUTMUO; |
2817d3e2 | 69 | case kBCUTE: |
4e994a92 | 70 | return fgkBCUTE; |
2817d3e2 | 71 | case kBCUTM: |
4e994a92 | 72 | return fgkBCUTM; |
2817d3e2 | 73 | case kDCUTE: |
4e994a92 | 74 | return fgkDCUTE; |
2817d3e2 | 75 | case kDCUTM: |
4e994a92 | 76 | return fgkDCUTM; |
2817d3e2 | 77 | case kPPCUTM: |
4e994a92 | 78 | return fgkPPCUTM; |
2817d3e2 | 79 | default: |
80 | TG4Globals::Warning("TG4G3Defaults::CutValue: Inconsistent cut."); | |
81 | return 0.; | |
82 | } | |
83 | } | |
84 | ||
73dd85b2 | 85 | //_____________________________________________________________________________ |
9187be37 | 86 | TG4G3ControlValue TG4G3Defaults::ControlValue(G4int control) |
2817d3e2 | 87 | { |
9187be37 | 88 | // Returns the G3 default value for the specified control. |
2817d3e2 | 89 | // --- |
90 | ||
9187be37 | 91 | switch (control) { |
2817d3e2 | 92 | case kPAIR: |
4e994a92 | 93 | return fgkPAIR; |
2817d3e2 | 94 | case kCOMP: |
4e994a92 | 95 | return fgkCOMP; |
2817d3e2 | 96 | case kPHOT: |
4e994a92 | 97 | return fgkPHOT; |
2817d3e2 | 98 | case kPFIS: |
4e994a92 | 99 | return fgkPFIS; |
2817d3e2 | 100 | case kDRAY: |
4e994a92 | 101 | return fgkDRAY; |
2817d3e2 | 102 | case kANNI: |
4e994a92 | 103 | return fgkANNI; |
2817d3e2 | 104 | case kBREM: |
4e994a92 | 105 | return fgkBREM; |
2817d3e2 | 106 | case kHADR: |
4e994a92 | 107 | return fgkHADR; |
2817d3e2 | 108 | case kMUNU: |
4e994a92 | 109 | return fgkMUNU; |
2817d3e2 | 110 | case kDCAY: |
4e994a92 | 111 | return fgkDCAY; |
2817d3e2 | 112 | case kLOSS: |
4e994a92 | 113 | return fgkLOSS; |
2817d3e2 | 114 | case kMULS: |
4e994a92 | 115 | return fgkMULS; |
2817d3e2 | 116 | default: |
9187be37 | 117 | TG4Globals::Warning( |
118 | "TG4G3Defaults::ControlValue: Inconsistent control."); | |
2817d3e2 | 119 | return kUnset; |
120 | } | |
121 | } | |
122 | ||
73dd85b2 | 123 | //_____________________________________________________________________________ |
9187be37 | 124 | G4bool TG4G3Defaults::IsDefaultCut(TG4G3Cut cut, G4double value) |
2817d3e2 | 125 | { |
126 | // Tests if the parameter value is equal to the G3 default value. | |
127 | // --- | |
128 | ||
9187be37 | 129 | if (abs(value*GeV - CutValue(cut)) > fgkCutTolerance) |
2817d3e2 | 130 | return false; |
131 | else | |
132 | return true; | |
133 | } | |
134 | ||
73dd85b2 | 135 | //_____________________________________________________________________________ |
9187be37 | 136 | G4bool TG4G3Defaults::IsDefaultControl(TG4G3Control control, |
137 | G4double value) | |
2817d3e2 | 138 | { |
139 | // Tests if the parameter value is equal to the G3 default value. | |
140 | // --- | |
141 | ||
9187be37 | 142 | if (value == ControlValue(control)) |
2817d3e2 | 143 | return true; |
144 | else | |
145 | return false; | |
146 | } |