PamG4RunConfiguration/src/PamG4RunConfiguration.cxx

#include <TG4ComposedPhysicsList.h>
#include "PamG4RunConfiguration.h"
#include <TG4EmPhysicsList.h>
#include <TG4HadronPhysicsList.h>
#include <TG4OpticalPhysicsList.h>
#include <TG4SpecialPhysicsList.h>
#include <TG4Globals.h>


#include "PamG4RangeCutsPhysicsList.h"
ClassImp(PamG4RunConfiguration)


PamG4RunConfiguration::PamG4RunConfiguration(const TString& userGeometry,
                                             const TString& physicsList,
                                             const TString& specialProcess)
  : TG4RunConfiguration(userGeometry, physicsList, specialProcess) 
{
  
}


PamG4RunConfiguration::~PamG4RunConfiguration()
{
/// Destructor
}


G4VUserPhysicsList* PamG4RunConfiguration::CreatePhysicsList()
{
// Create default Geant4 VMC physics list

  TG4ComposedPhysicsList* builder = new TG4ComposedPhysicsList();
  
  TString token1 = TG4Globals::GetToken(0, fPhysicsListSelection);
  TString token2 = TG4Globals::GetToken(1, fPhysicsListSelection);
  
  if ( token1 == "emStandard" ) {
  G4cout << "Adding EMPhysicsList" << G4endl;
  builder->AddPhysicsList(new TG4EmPhysicsList());
}  
else {
  G4cout << "Adding HadronPhysicsList " << token1.Data() << G4endl;
  builder->AddPhysicsList(new TG4HadronPhysicsList(token1.Data()));
}  
  
if ( token2 == "optical" ) {
  G4cout << "Adding OpticalPhysicsList " << G4endl;
  builder->AddPhysicsList(new TG4OpticalPhysicsList());
}  
  
 // add option here
  G4cout << "Adding SpecialPhysicsList " << G4endl;
  builder->AddPhysicsList(new TG4SpecialPhysicsList(
                               fSpecialProcessSelection.Data()));

  G4cout << "Adding Range Cuts PhysicsList " << G4endl;
  builder->AddPhysicsList(new PamG4RangeCutsPhysicsList());

return builder;  
}  
1	#include <TG4ComposedPhysicsList.h>
2	#include "PamG4RunConfiguration.h"
3	#include <TG4EmPhysicsList.h>
4	#include <TG4HadronPhysicsList.h>
5	#include <TG4OpticalPhysicsList.h>
6	#include <TG4SpecialPhysicsList.h>
7	#include <TG4Globals.h>
8
9
10	#include "PamG4RangeCutsPhysicsList.h"
11	ClassImp(PamG4RunConfiguration)
12
13
14	PamG4RunConfiguration::PamG4RunConfiguration(const TString& userGeometry,
15	const TString& physicsList,
16	const TString& specialProcess)
17	: TG4RunConfiguration(userGeometry, physicsList, specialProcess)
18	{
19
20	}
21
22
23	PamG4RunConfiguration::~PamG4RunConfiguration()
24	{
25	/// Destructor
26	}
27
28
29	G4VUserPhysicsList* PamG4RunConfiguration::CreatePhysicsList()
30	{
31	// Create default Geant4 VMC physics list
32
33	TG4ComposedPhysicsList* builder = new TG4ComposedPhysicsList();
34
35	TString token1 = TG4Globals::GetToken(0, fPhysicsListSelection);
36	TString token2 = TG4Globals::GetToken(1, fPhysicsListSelection);
37
38	if ( token1 == "emStandard" ) {
39	G4cout << "Adding EMPhysicsList" << G4endl;
40	builder->AddPhysicsList(new TG4EmPhysicsList());
41	}
42	else {
43	G4cout << "Adding HadronPhysicsList " << token1.Data() << G4endl;
44	builder->AddPhysicsList(new TG4HadronPhysicsList(token1.Data()));
45	}
46
47	if ( token2 == "optical" ) {
48	G4cout << "Adding OpticalPhysicsList " << G4endl;
49	builder->AddPhysicsList(new TG4OpticalPhysicsList());
50	}
51
52	// add option here
53	G4cout << "Adding SpecialPhysicsList " << G4endl;
54	builder->AddPhysicsList(new TG4SpecialPhysicsList(
55	fSpecialProcessSelection.Data()));
56
57	G4cout << "Adding Range Cuts PhysicsList " << G4endl;
58	builder->AddPhysicsList(new PamG4RangeCutsPhysicsList());
59
60	return builder;
61	}