PamG4RunConfiguration_exp/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()));

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

return builder;  
}  
1	pamelats	1.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			/*R G4cout << "Adding Range Cuts PhysicsList " << G4endl;
58			builder->AddPhysicsList(new PamG4RangeCutsPhysicsList());*/
59
60			return builder;
61			}