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	cafagna	1.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