CustomPhysicsTable.cxx

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \file CustomPhysicsTable.cxx
//
/// \author  bjpjones@mit.edu
////////////////////////////////////////////////////////////////////////
// See header file for full description

#include "messagefacility/MessageLogger/MessageLogger.h"

#include "larsim/LegacyLArG4/CustomPhysicsFactory.hh"
#include "larsim/LegacyLArG4/CustomPhysicsTable.hh"

namespace larg4 {

  static CustomPhysicsTable* TheCustomPhysicsTable;

  //-----------------------------------------------------------------
  CustomPhysicsTable::CustomPhysicsTable(CustomPhysicsFactoryBase* Factory)
  {
    if (!TheCustomPhysicsTable) {
      TheCustomPhysicsTable = new CustomPhysicsTable;
      TheCustomPhysicsTable->AddPhysics(Factory);
    }
    else {
      TheCustomPhysicsTable->AddPhysics(Factory);
    }
  }

  //-----------------------------------------------------------------
  std::vector<std::string>
  CustomPhysicsTable::GetAvailablePhysicsList()
  {
    std::vector<std::string> ReturnVector;
    for (std::map<std::string, CustomPhysicsFactoryBase*>::const_iterator i =
           TheCustomPhysicsTable->theTable.begin();
         i != TheCustomPhysicsTable->theTable.end();
         i++) {
      ReturnVector.push_back((*i).first);
    }
    return ReturnVector;
  }

  //-----------------------------------------------------------------
  bool
  CustomPhysicsTable::IsPhysicsAvailable(std::string const& PhysicsName)
  {
    if (!TheCustomPhysicsTable->theTable[PhysicsName])
      return false;
    else
      return true;
  }

  //-----------------------------------------------------------------
  G4VPhysicsConstructor*
  CustomPhysicsTable::GetPhysicsConstructor(std::string const& PhysicsName)
  {
    if (IsPhysicsAvailable(PhysicsName)) {
      return TheCustomPhysicsTable->theTable[PhysicsName]->Build();
    }

    return nullptr;
  }

  //-----------------------------------------------------------------
  void
  CustomPhysicsTable::AddPhysics(CustomPhysicsFactoryBase* Factory)
  {

    if (IsPhysicsAvailable(Factory->GetName()))
      mf::LogWarning("CustomPhysicsTable") << "Physics constructor being overwritten"
                                           << " in CustomPhysicsTable";
    TheCustomPhysicsTable->theTable[Factory->GetName()] = Factory;
    MF_LOG_DEBUG("CustomPhysicsTable")
      << "CustomPhysicsTable : Physics Table registering new physics " << Factory->GetName();
  }

}