PhysicsLIst

Dear VMC team,

I am trying to use a custom PhysicsList of G4 and have some problems with it. As an example I took
http://www.slac.stanford.edu/comp/physics/geant4/slac_physics_lists/ilc/ physics list and created a new class NRunConfiguration:

//_____________________________________________________________________________
NRunConfiguration::NRunConfiguration()
  : TG4RunConfiguration() {
//

  CreateUserConfiguration();
}

//_____________________________________________________________________________
NRunConfiguration::~NRunConfiguration(){
//

}

//
// protected methods
//

//_____________________________________________________________________________
void NRunConfiguration::CreateUserConfiguration() {
// Creates the mandatory Geant4 classes and
// the other user action classes.
// ---

    // create mandatory Geant4 classes
    fPhysicsList = new LCPhysicsList();
    cout << "LCPhysicsList is activated" << endl; }

in the main program I do:

    TG4RunConfiguration* runConfiguration = new NRunConfiguration();

    // TGeant4
    TGeant4 *geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
and finally my TMCApplication object calls InitMC in which a gMC->BuildPhysics() call is done.
The program is compiled and runs smoothly but apparently it does not have hadronic interactions while a message ###
You are using the Geant4 VMC hadronics physics equivalent to ExN04HadronPhysics from Geant4 (version 6.0) Note that this hadronic physics list is not optimized for any particular usage.
If you wish to have a starting point tailored for a particular area of work, please use one of the available physics lists by use-case.  More information can also be found from the Geant4 HyperNews. ###
is released what I see strange since I do not call anymore VMC default PhysicsList constructor. Could you please explain what is a safe way to use: 1. custom PhysicsList
2. activate one or two new modules like QGSP for hadronic interactions. It seems I understand how to do it in a pure Geant4 environment while I am not sure I understand how it works in VMC interface.

Many thanks in advance, Dmitry Received on Mon May 08 2006 - 16:12:10 MEST