PamG4/src/PhysicsList.cc



#include "PhysicsList.hh"

#include "globals.hh"
#include "G4ParticleDefinition.hh"
#include "G4ParticleWithCuts.hh"
#include "G4ProcessManager.hh"
#include "G4ProcessVector.hh"
#include "G4ParticleTypes.hh"
#include "G4ParticleTable.hh"

#include "G4Material.hh"
#include "G4MaterialTable.hh"
#include "G4ios.hh"

#include "GeneralPhysics.hh"
#include "EMPhysics.hh"
#include "MuonPhysics.hh"
#include "HadronPhysics.hh"
#include "IonPhysics.hh"

PhysicsList::PhysicsList():  G4VModularPhysicsList()
{
  // default cut value  (1.0mm) 
  defaultCutValue = 0.1*mm;
  SetVerboseLevel(1);

  // General Physics
  RegisterPhysics( new GeneralPhysics("general") );

  // EM Physics
  RegisterPhysics( new EMPhysics("standard EM"));

  // Muon Physics
  RegisterPhysics(  new MuonPhysics("muon"));

   // Hadron Physics
  RegisterPhysics(  new HadronPhysics("hadron"));

  // Ion Physics
  RegisterPhysics( new IonPhysics("ion"));


}

PhysicsList::~PhysicsList()
{
}

void PhysicsList::SetCuts()
{
  //  " G4VUserPhysicsList::SetCutsWithDefault" method sets 
  //   the default cut value for all particle types 
  SetCutsWithDefault();   
  G4cout<<"CutLength : "<<G4BestUnit(defaultCutValue,"Length")<<G4endl;
  if(verboseLevel>0) DumpCutValuesTable();
}


1
2
3	#include "PhysicsList.hh"
4
5	#include "globals.hh"
6	#include "G4ParticleDefinition.hh"
7	#include "G4ParticleWithCuts.hh"
8	#include "G4ProcessManager.hh"
9	#include "G4ProcessVector.hh"
10	#include "G4ParticleTypes.hh"
11	#include "G4ParticleTable.hh"
12
13	#include "G4Material.hh"
14	#include "G4MaterialTable.hh"
15	#include "G4ios.hh"
16
17	#include "GeneralPhysics.hh"
18	#include "EMPhysics.hh"
19	#include "MuonPhysics.hh"
20	#include "HadronPhysics.hh"
21	#include "IonPhysics.hh"
22
23	PhysicsList::PhysicsList(): G4VModularPhysicsList()
24	{
25	// default cut value (1.0mm)
26	defaultCutValue = 0.1*mm;
27	SetVerboseLevel(1);
28
29	// General Physics
30	RegisterPhysics( new GeneralPhysics("general") );
31
32	// EM Physics
33	RegisterPhysics( new EMPhysics("standard EM"));
34
35	// Muon Physics
36	RegisterPhysics( new MuonPhysics("muon"));
37
38	// Hadron Physics
39	RegisterPhysics( new HadronPhysics("hadron"));
40
41	// Ion Physics
42	RegisterPhysics( new IonPhysics("ion"));
43
44
45	}
46
47	PhysicsList::~PhysicsList()
48	{
49	}
50
51	void PhysicsList::SetCuts()
52	{
53	// " G4VUserPhysicsList::SetCutsWithDefault" method sets
54	// the default cut value for all particle types
55	SetCutsWithDefault();
56	G4cout<<"CutLength : "<<G4BestUnit(defaultCutValue,"Length")<<G4endl;
57	if(verboseLevel>0) DumpCutValuesTable();
58	}
59
60
61