TGeant4/TG4G3Defaults.cxx

   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; break;
  58     case kCUTELE:
  59       return fgkCUTELE; break;
  60     case kCUTNEU:
  61       return fgkCUTNEU; break;
  62     case kCUTHAD:
  63       return fgkCUTHAD; break;
  64     case kCUTMUO:
  65       return fgkCUTMUO; break;
  66     case kBCUTE:
  67       return fgkBCUTE;  break;
  68     case kBCUTM:
  69       return fgkBCUTM;  break;
  70     case kDCUTE:
  71       return fgkDCUTE;  break;
  72     case kDCUTM:
  73       return fgkDCUTM;  break;
  74     case kPPCUTM:
  75       return fgkPPCUTM; break;
  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; break;
  90     case kCOMP:
  91       return fgkCOMP; break;
  92     case kPHOT:
  93       return fgkPHOT; break;
  94     case kPFIS:
  95       return fgkPFIS; break;
  96     case kDRAY:
  97       return fgkDRAY; break;
  98     case kANNI:
  99       return fgkANNI; break;
 100     case kBREM:
 101       return fgkBREM; break;
 102     case kHADR:
 103       return fgkHADR; break;
 104     case kMUNU:
 105       return fgkMUNU; break;
 106     case kDCAY:
 107       return fgkDCAY; break;
 108     case kLOSS:
 109       return fgkLOSS; break;
 110     case kMULS:
 111       return fgkMULS; break;
 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 }