GlobalPhysicsList.cc

Go to the documentation of this file.

#include "GlobalPhysicsList.h"

#include <utl/ErrorLogger.h>


#include "DefaultPhysicsList.h"

using namespace tls;
using namespace std;


GlobalPhysicsList::GlobalPhysicsList()
{
  SetVerboseLevel(-10);
}


GlobalPhysicsList::~GlobalPhysicsList()
{
  INFO("destructing global physics list");
  for (vector<G4VUserPhysicsList*>::iterator it = fPhysicsVector.begin();
       it!= fPhysicsVector.end(); ++it) {
    delete (*it);
  }
  fPhysicsVector.clear();
}


void
GlobalPhysicsList::ConstructParticle()
{
  for (vector<G4VUserPhysicsList*>::iterator it = fPhysicsVector.begin();
       it!= fPhysicsVector.end(); ++it) {
    (*it)->ConstructParticle();
  }
}


void
GlobalPhysicsList::ConstructProcess()
{
  AddTransportation();

  for (vector<G4VUserPhysicsList*>::iterator it = fPhysicsVector.begin();
       it!= fPhysicsVector.end(); ++it) {
    (*it)->ConstructProcess();
  }
}


void
GlobalPhysicsList::SetVerboseLevel(bool verbosity)
{
  for (vector<G4VUserPhysicsList*>::iterator it = fPhysicsVector.begin();
       it!= fPhysicsVector.end(); ++it) {
    (*it)->SetVerboseLevel(verbosity);
  }
}


void
GlobalPhysicsList::SetCuts()
{
  this->SetCutsWithDefault();
}


const
G4VUserPhysicsList*
GlobalPhysicsList::GetPhysics(const G4String& name)
  const
{
  const unsigned int n = fPhysicsVectorNames.size();
  for (unsigned int i = 0; i != n; ++i) {
    if (fPhysicsVectorNames[i] == name)
      return fPhysicsVector[i];
  }
  return 0;
}


void
GlobalPhysicsList::RegisterPhysics(const G4String& name, G4VUserPhysicsList* userList)
{
  vector<G4String>::iterator it =
    find(fPhysicsVectorNames.begin(), fPhysicsVectorNames.end(), name);
  if (it == fPhysicsVectorNames.end()) {
    fPhysicsVectorNames.push_back(name);
    fPhysicsVector.push_back(userList);
  }
  else {
    const unsigned int i = distance(fPhysicsVectorNames.begin(), it);
    fPhysicsVector[i] = userList;
  }
}