PhysicsListFactory.cc

Go to the documentation of this file.

// $Header: /cvs/lcd/slic/src/PhysicsListFactory.cc,v 1.22 2007/10/22 13:26:07 wenzel Exp $
#include "PhysicsListFactory.hh"

// geant4 physics lists
#ifdef HAVE_G4PHYSICSLISTS
#include "FTFC.hh"
#include "FTFP.hh"
#ifdef HAVE_FTFP_EMV
#include "FTFP_EMV.hh"
#endif
#include "LBE.hh"
#include "LHEP.hh"
#include "LHEP_BERT.hh"
#include "LHEP_BERT_HP.hh"
#include "LHEP_EMV.hh"
#include "LHEP_PRECO_HP.hh"
#include "QBBC.hh"
#include "QGSC.hh"
#include "QGSC_EFLOW.hh"
#include "QGSC_EMV.hh"
#include "QGSP.hh"
#include "QGSP_BERT.hh"
#ifdef HAVE_QGSP_BERT_EMV
#include "QGSP_BERT_EMV.hh"
#endif
#include "QGSP_BERT_HP.hh"
#ifdef HAVE_QGSP_BERT_NQE
#include "QGSP_BERT_NQE.hh"
#endif
#ifdef HAVE_QGSP_BERT_TRV
#include "QGSP_BERT_TRV.hh"
#endif
#include "QGSP_BIC.hh"
#include "QGSP_BIC_HP.hh"
#include "QGSP_EMV.hh"
#ifdef HAVE_QGSP_EMV_NQE
#include "QGSP_EMV_NQE.hh"
#endif
#include "QGSP_EMX.hh"
#ifdef HAVE_QGSP_NQE
#include "QGSP_NQE.hh"
#endif
#include "QGSP_QEL.hh"
#endif

// lcphys 
#ifdef SLIC_USE_LCPHYS
#include "LCPhysicsList.hh"
#endif

// slic
#include "PhysicsListBuilder.hh"

using namespace std;

#define BUILDER_LCPHYS(list,name,description) this->addBuilder( new TPhysicsListBuilder<list>( name, description ) );

#define BUILDER(list,description) this->addBuilder( new TPhysicsListBuilder<list>( #list, description ) );

namespace slic
{

  PhysicsListFactory::PhysicsListFactory()
    : Module("PhysicsListFactory")
  {
    initializeBuilders();
  }
  
  PhysicsListFactory::~PhysicsListFactory()
  {
    for ( BuilderMap::const_iterator it = m_builders.begin();
          it != m_builders.end();
          it++ ) {
      delete it->second;
    }      
    m_builders.clear();
  }
  
  G4VUserPhysicsList* PhysicsListFactory::create( const std::string& name )
  {    
    log().okay("Creating physics list <" + name + ">.");

    G4VUserPhysicsList* physicslist = 0;
    
    VPhysicsListBuilder* builder = m_builders[ name ];

    if ( builder != 0 ) {
      physicslist = builder->create();
    }     
    else {
      G4Exception("ERROR: No builder found for physics list <" + name + "> .");
    }

    return physicslist;
  }
  
  void PhysicsListFactory::printAvailableLists()
  {
    log() << LOG::okay << LOG::done;
    log() << LOG::okay << "----GEANT4 PHYSICS LISTS----" << LOG::endl << LOG::done;
    log() << LOG::okay << "NAME           DESCRIPTION" << LOG::done;
    log() << LOG::okay << "----------------------------------------------------------------------" << LOG::done;
    for ( BuilderMap::const_iterator it = m_builders.begin();
          it != m_builders.end();
          it++ ) {
      log().width(15);
      log() << LOG::okay << left << it->first;
      log() << LOG::okay << it->second->getDescription() << LOG::done;
    }
  }

  void PhysicsListFactory::initializeBuilders()
  {

#ifdef SLIC_USE_LCPHYS
    BUILDER_LCPHYS(LCPhysicsList, "LCPhys", "ILC Physics List by Dennis Wright.")
#endif

#ifdef HAVE_G4PHYSICSLISTS
    BUILDER(FTFC,          "Fritiof with CHIPS")
    BUILDER(FTFP,          "Fritiof with precompound")
#ifdef HAVE_FTFP_EMV
    BUILDER(FTFP_EMV,      "FTFP with v 7.1 EM processes")
#endif
    BUILDER(LBE,           "Low backgrounds experiment with low energy EM")
    BUILDER(LHEP,          "Low and High Energy Parameterised")
    BUILDER(LHEP_BERT,     "LHEP with Bertini Cascade")
    BUILDER(LHEP_BERT_HP,  "LHEP with Bertini Cascade and high precision neutrons")
    BUILDER(LHEP_EMV,      "LHEP with v 7.1 EM processes")
    BUILDER(LHEP_PRECO_HP, "LHEP_PRECO with HP")
    BUILDER(QBBC,          "Uses best cross section data within QGSP")
    BUILDER(QGSC,          "Quark Gluon-String (QGS) with CHIPS")
    BUILDER(QGSC_EFLOW,    "QGSC using energy flow methods from CHIPS model")
    BUILDER(QGSC_EMV,      "QGSC with v 7.1 EM processes")
    BUILDER(QGSP,          "QGS with precompound")
    BUILDER(QGSP_BERT,     "QGSP with Bertini Cascade")
#ifdef HAVE_QGSP_BERT_EMV
    BUILDER(QGSP_BERT_EMV, "QGSP_BERT with v 7.1 EM processes")
#endif
    BUILDER(QGSP_BERT_HP,  "QGSP with Bertini Cascade and high precision neutrons")
#ifdef HAVE_QGSP_BERT_NQE
    BUILDER(QGSP_BERT_NQE, "QGSP_BERT with quasi-elastic processes removed")
#endif
#ifdef HAVE_QGSP_BERT_TRV
    BUILDER(QGSP_BERT_TRV, "QGSP_BERT with transition to Bertini at 5 GeV")
#endif
    BUILDER(QGSP_BIC,      "QGSP with Binary Cascade")
    BUILDER(QGSP_BIC_HP,   "QGSP with Binary Cascade")
    BUILDER(QGSP_EMV,      "LHEP with v 7.1 EM processes")
#ifdef HAVE_QGSP_EMV_NQE
    BUILDER(QGSP_EMV_NQE,  "QGSP with v 7.1 EM processes and quasi-elastic processes removed")
#endif
    BUILDER(QGSP_EMX,      "QGSP with v 7.2 EM processes")
#ifdef HAVE_QGSP_NQE
    BUILDER(QGSP_NQE,      "QGSP with quasi-elastic processes removed")
#endif
    BUILDER(QGSP_QEL,      "QGSP with elastic cross sections and model from CHIPS")
#endif
  }

  void PhysicsListFactory::addBuilder( VPhysicsListBuilder* builder)
  {
    const G4String& name = builder->getName();

    if ( m_builders.find( name ) != m_builders.end() ) {
      log() << LOG::warning << "Ignoring duplicate physics list builder for <" << name << "> ." << LOG::done;
    }
    else {
      m_builders[ name ] = builder;
    }
  }
}

/*

physics lists in Geant4 9.0

*/

---