TG4HadronPhysicsList.cxx

Go to the documentation of this file.

//------------------------------------------------
// The Geant4 Virtual Monte Carlo package
// Copyright (C) 2007 - 2014 Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see geant4_vmc/LICENSE.
// Contact: root-vmc@cern.ch
//-------------------------------------------------

 
#include "TG4HadronPhysicsList.h"
 
#include <G4SystemOfUnits.hh>
 
#include "G4PhysListFactory.hh"
 
const G4double TG4HadronPhysicsList::fgkDefaultCutValue = 1.0 * mm;
 
//
// static methods
//
 
//_____________________________________________________________________________
G4String TG4HadronPhysicsList::AvailableHadronSelections()
{
 
  G4PhysListFactory phyListFactory;
  const std::vector<G4String>& availablePhysLists =
    phyListFactory.AvailablePhysLists();
 
  G4String selections;
  for (G4int i = 0; i < G4int(availablePhysLists.size()); ++i) {
    selections += availablePhysLists[i];
    selections += " ";
  }
  // selections += "ShieldingLEND ";
 
  return selections;
}
 
//_____________________________________________________________________________
G4String TG4HadronPhysicsList::AvailableEMSelections()
{
 
  G4PhysListFactory phyListFactory;
  const std::vector<G4String>& availablePhysListsEM =
    phyListFactory.AvailablePhysListsEM();
 
  G4String selections;
  for (G4int i = 0; i < G4int(availablePhysListsEM.size()); ++i) {
    selections += availablePhysListsEM[i];
    selections += " ";
  }
 
  return selections;
}
 
//_____________________________________________________________________________
G4bool TG4HadronPhysicsList::IsAvailableSelection(const G4String& selection)
{
 
  G4PhysListFactory physListFactory;
  return physListFactory.IsReferencePhysList(selection);
}
 
//
// ctors, dtor
//
 
//_____________________________________________________________________________
TG4HadronPhysicsList::TG4HadronPhysicsList(const G4String& selection)
  : G4VUserPhysicsList(), TG4Verbose("hadronPhysicsList")
{
 
  Configure(selection);
 
  defaultCutValue = fgkDefaultCutValue;
 
  SetVerboseLevel(TG4Verbose::VerboseLevel());
}
 
//_____________________________________________________________________________
TG4HadronPhysicsList::~TG4HadronPhysicsList()
{
 
  // delete fExtDecayer;
  // fExtDecayer is deleted in G4Decay destructor
}
 
//
// private methods
//
 
//_____________________________________________________________________________
void TG4HadronPhysicsList::Configure(const G4String& sel)
{
 
  G4PhysListFactory phyListFactory;
 
  fPhysicsList = phyListFactory.GetReferencePhysList(sel);
}
 
//
// public methods
//
 
//_____________________________________________________________________________
void TG4HadronPhysicsList::ConstructParticle()
{
 
  // create processes for registered physics
  fPhysicsList->ConstructParticle();
}
 
//_____________________________________________________________________________
void TG4HadronPhysicsList::ConstructProcess()
{
 
  // create processes for registered physics
  fPhysicsList->ConstructProcess();
 
  if (VerboseLevel() > 0)
    G4cout << "### Hadron physics constructed. " << G4endl;
}
 
//_____________________________________________________________________________
G4int TG4HadronPhysicsList::VerboseLevel() const
{
 
  return TG4Verbose::VerboseLevel();
}
 
//_____________________________________________________________________________
void TG4HadronPhysicsList::VerboseLevel(G4int level)
{
 
  TG4Verbose::VerboseLevel(level);
  fPhysicsList->SetVerboseLevel(level);
}
 
//_____________________________________________________________________________
void TG4HadronPhysicsList::SetRangeCut(G4double value)
{
 
  defaultCutValue = value;
}