Update of SSD simulation and reconstruction code by Boris and Enrico.

[u/mrichter/AliRoot.git] / TGeant4 / TG4G3Defaults.cxx

// $Id$
// Category: physics
//
// See the class description in the header file.

#include "TG4G3Defaults.h"
#include "TG4Globals.h"

#include <math.h>

// static const data members

// precision tolerance
const G4double TG4G3Defaults::fgkCutTolerance = 1. * keV;

// kinetic energy cuts
const G4double TG4G3Defaults::fgkCUTGAM = 0.001 * GeV;
const G4double TG4G3Defaults::fgkCUTELE = 0.001 * GeV;
const G4double TG4G3Defaults::fgkCUTNEU = 0.01 * GeV;
const G4double TG4G3Defaults::fgkCUTHAD = 0.01 * GeV;
const G4double TG4G3Defaults::fgkCUTMUO = 0.01 * GeV;
const G4double TG4G3Defaults::fgkBCUTE  = fgkCUTGAM;
const G4double TG4G3Defaults::fgkBCUTM  = fgkCUTGAM;
const G4double TG4G3Defaults::fgkDCUTE  = 10. * TeV;
const G4double TG4G3Defaults::fgkDCUTM  = 10. * TeV;
const G4double TG4G3Defaults::fgkPPCUTM = 0.01 * GeV;

// physics processes
const TG4G3ControlValue TG4G3Defaults::fgkPAIR = kActivate;   // 1
const TG4G3ControlValue TG4G3Defaults::fgkCOMP = kActivate;   // 1
const TG4G3ControlValue TG4G3Defaults::fgkPHOT = kActivate;   // 1
const TG4G3ControlValue TG4G3Defaults::fgkPFIS = kInActivate; // 0
const TG4G3ControlValue TG4G3Defaults::fgkDRAY = kActivate2;  // 2
const TG4G3ControlValue TG4G3Defaults::fgkANNI = kActivate;   // 1
const TG4G3ControlValue TG4G3Defaults::fgkBREM = kActivate;   // 1
const TG4G3ControlValue TG4G3Defaults::fgkHADR = kActivate;   // 1
const TG4G3ControlValue TG4G3Defaults::fgkMUNU = kInActivate; // 0
const TG4G3ControlValue TG4G3Defaults::fgkDCAY = kActivate;   // 1
const TG4G3ControlValue TG4G3Defaults::fgkLOSS = kActivate2;  // 2
const TG4G3ControlValue TG4G3Defaults::fgkMULS = kActivate;   // 1

TG4G3Defaults::TG4G3Defaults() {
//
}
  
TG4G3Defaults::~TG4G3Defaults() {
//
}

G4double TG4G3Defaults::CutValue(G4int g3Cut)
{
// Returns the G3 default value for the specified cut.
// ---

  switch (g3Cut) {
    case kCUTGAM:  
      return fgkCUTGAM;
    case kCUTELE:  
      return fgkCUTELE;
    case kCUTNEU:  
      return fgkCUTNEU;
    case kCUTHAD:  
      return fgkCUTHAD;
    case kCUTMUO:  
      return fgkCUTMUO;
    case kBCUTE:   
      return fgkBCUTE; 
    case kBCUTM:   
      return fgkBCUTM; 
    case kDCUTE:   
      return fgkDCUTE; 
    case kDCUTM:   
      return fgkDCUTM; 
    case kPPCUTM:  
      return fgkPPCUTM;
    default:
      TG4Globals::Warning("TG4G3Defaults::CutValue: Inconsistent cut.");
      return 0.;      
  }
}          

TG4G3ControlValue TG4G3Defaults::ControlValue(G4int control)
{
// Returns the G3 default value for the specified control.
// ---

  switch (control) {
    case kPAIR:
      return fgkPAIR;
    case kCOMP:
      return fgkCOMP;
    case kPHOT:
      return fgkPHOT;
    case kPFIS:
      return fgkPFIS;
    case kDRAY:
      return fgkDRAY;
    case kANNI:
      return fgkANNI;
    case kBREM:
      return fgkBREM;
    case kHADR:
      return fgkHADR;
    case kMUNU:
      return fgkMUNU;
    case kDCAY:
      return fgkDCAY;
    case kLOSS:
      return fgkLOSS;
    case kMULS:
      return fgkMULS;
    default:
      TG4Globals::Warning(
        "TG4G3Defaults::ControlValue: Inconsistent control.");
      return kUnset;      
  }
}          

G4bool TG4G3Defaults::IsDefaultCut(TG4G3Cut cut, G4double value)
{
// Tests if the parameter value is equal to the G3 default value.
// ---

  if (abs(value*GeV - CutValue(cut)) > fgkCutTolerance) 
    return false;
  else  
    return true;
}

G4bool TG4G3Defaults::IsDefaultControl(TG4G3Control control,
                                              G4double value)
{
// Tests if the parameter value is equal to the G3 default value.
// ---

  if (value == ControlValue(control)) 
    return true;
  else  
    return false;
}