[u/mrichter/AliRoot.git] / TGeant4 / TG4PhysicsConstructorMuon.cxx

// $Id$
// Category: physics
//
// Author: I. Hrivnacova
//
// Class TG4PhysicsConstructorMuon
// -----------------------------
// See the class description in the header file.
// According to ExN04MuonPhysics.cc,v 1.2.2.1 2001/06/28 19:07:37 gunter Exp 
// GEANT4 tag Name: geant4-03-02

#include "TG4PhysicsConstructorMuon.h"
#include "TG4ProcessControlMap.h"
#include "TG4ProcessMCMap.h"
#include "TG4G3Control.h"

#include <G4ParticleDefinition.hh>
#include <G4ProcessManager.hh>
#include <G4MuonPlus.hh>
#include <G4MuonMinus.hh>
#include <G4TauMinus.hh>
#include <G4TauPlus.hh>
#include <G4NeutrinoTau.hh>
#include <G4AntiNeutrinoTau.hh>
#include <G4NeutrinoMu.hh>
#include <G4AntiNeutrinoMu.hh>

//TBR
#include <G4MuIonisation.hh>
#include <G4MuBremsstrahlung.hh>
#include <G4MuPairProduction.hh>
#include <G4hIonisation.hh>

//_____________________________________________________________________________
TG4PhysicsConstructorMuon::TG4PhysicsConstructorMuon(const G4String& name)
  : G4VPhysicsConstructor(name)
{
//
  SetVerboseLevel(1);
}

//_____________________________________________________________________________
TG4PhysicsConstructorMuon::~TG4PhysicsConstructorMuon() {
//
}

// private methods

//_____________________________________________________________________________
void TG4PhysicsConstructorMuon::ConstructProcessForMuonPlus()
{
// Constructs electromagnetic processes for mu+.
// ---
  
  // add processes
  G4ProcessManager* pManager = G4MuonPlus::MuonPlus()->GetProcessManager();
  pManager->AddProcess(&fMuPlusIonisation, ordInActive,2, 2);
  pManager->AddDiscreteProcess(&fMuPlusBremsstrahlung);
  pManager->AddDiscreteProcess(&fMuPlusPairProduction);
  pManager->AddProcess(&fMuPlusMultipleScattering);

  // set ordering
  pManager->SetProcessOrdering(&fMuPlusMultipleScattering, idxAlongStep,  1);
  pManager->SetProcessOrdering(&fMuPlusMultipleScattering, idxPostStep,  1);

  // map to G3 controls
  TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
  controlMap->Add(&fMuPlusIonisation, kLOSS); 
  controlMap->Add(&fMuPlusBremsstrahlung, kBREM); 
  controlMap->Add(&fMuPlusPairProduction, kPAIR); 
  controlMap->Add(&fMuPlusMultipleScattering, kMULS); 


  // map to AliMCProcess codes
  TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
  mcMap->Add(&fMuPlusIonisation, kPEnergyLoss); 
  mcMap->Add(&fMuPlusBremsstrahlung, kPBrem); 
  mcMap->Add(&fMuPlusPairProduction, kPPair); 
  mcMap->Add(&fMuPlusMultipleScattering, kPMultipleScattering); 
}  

//_____________________________________________________________________________
void TG4PhysicsConstructorMuon::ConstructProcessForMuonMinus()
{
// Constructs electromagnetic processes for mu-.
// ---
  
  // add processes & set ordering
  G4ProcessManager* pManager = G4MuonMinus::MuonMinus()->GetProcessManager();
  pManager->AddProcess(&fMuMinusIonisation, ordInActive,2, 2);
  pManager->AddDiscreteProcess(&fMuMinusBremsstrahlung);
  pManager->AddDiscreteProcess(&fMuMinusPairProduction);
  pManager->AddProcess(&fMuMinusMultipleScattering);

  pManager->SetProcessOrdering(&fMuMinusMultipleScattering, idxAlongStep,  1);
  pManager->SetProcessOrdering(&fMuMinusMultipleScattering, idxPostStep,  1);

  pManager->AddRestProcess(&fMuMinusCaptureAtRest);

  // map to G3 controls
  TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
  controlMap->Add(&fMuMinusIonisation, kLOSS); 
  controlMap->Add(&fMuMinusBremsstrahlung, kBREM); 
  controlMap->Add(&fMuMinusPairProduction, kPAIR); 
  controlMap->Add(&fMuMinusMultipleScattering, kMULS); 
  controlMap->Add(&fMuMinusCaptureAtRest, kMUNU); 

  // map to AliMCProcess codes
  TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
  mcMap->Add(&fMuMinusIonisation, kPEnergyLoss); 
  mcMap->Add(&fMuMinusBremsstrahlung, kPBrem); 
  mcMap->Add(&fMuMinusPairProduction, kPPair); 
  mcMap->Add(&fMuMinusMultipleScattering, kPMultipleScattering); 
  mcMap->Add(&fMuMinusCaptureAtRest, kPMuonNuclear); 
}  

//_____________________________________________________________________________
void TG4PhysicsConstructorMuon::ConstructProcessForTauPlus()
{
// Constructs electromagnetic processes for tau+.
// ---
  
  // add processes
  G4ProcessManager* pManager = G4TauPlus::TauPlus()->GetProcessManager();
  pManager->AddProcess(&fTauPlusIonisation, ordInActive,2, 2);
  pManager->AddProcess(&fTauPlusMultipleScattering);

  // set ordering
  pManager->SetProcessOrdering(&fTauPlusMultipleScattering, idxAlongStep,  1);
  pManager->SetProcessOrdering(&fTauPlusMultipleScattering, idxPostStep,  1);

  // map to G3 controls
  TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
  controlMap->Add(&fTauPlusIonisation, kLOSS); 
  controlMap->Add(&fTauPlusMultipleScattering, kMULS); 


  // map to AliMCProcess codes
  TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
  mcMap->Add(&fTauPlusIonisation, kPEnergyLoss); 
  mcMap->Add(&fTauPlusMultipleScattering, kPMultipleScattering); 
}  

//_____________________________________________________________________________
void TG4PhysicsConstructorMuon::ConstructProcessForTauMinus()
{
// Constructs electromagnetic processes for tau-.
// ---
  
  // add processes
  G4ProcessManager* pManager = G4TauMinus::TauMinus()->GetProcessManager();
  pManager->AddProcess(&fTauMinusIonisation, ordInActive,2, 2);
  pManager->AddProcess(&fTauMinusMultipleScattering);

  // set ordering
  pManager->SetProcessOrdering(&fTauMinusMultipleScattering, idxAlongStep,  1);
  pManager->SetProcessOrdering(&fTauMinusMultipleScattering, idxPostStep,  1);

  // map to G3 controls
  TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
  controlMap->Add(&fTauMinusIonisation, kLOSS); 
  controlMap->Add(&fTauMinusMultipleScattering, kMULS); 


  // map to AliMCProcess codes
  TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
  mcMap->Add(&fTauMinusIonisation, kPEnergyLoss); 
  mcMap->Add(&fTauMinusMultipleScattering, kPMultipleScattering); 
}  


// protected methods

//_____________________________________________________________________________
void TG4PhysicsConstructorMuon::ConstructParticle()
{
// Instantiates particles.
// ---

  // Mu
  G4MuonPlus::MuonPlusDefinition();
  G4MuonMinus::MuonMinusDefinition();
  G4NeutrinoMu::NeutrinoMuDefinition();
  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();

  // Tau
  G4TauMinus::TauMinusDefinition();
  G4TauPlus::TauPlusDefinition();
  G4NeutrinoTau::NeutrinoTauDefinition();
  G4AntiNeutrinoTau::AntiNeutrinoTauDefinition();
}

//_____________________________________________________________________________
void TG4PhysicsConstructorMuon::ConstructProcess()
{
// Constructs electromagnetic processes for muons.
// ---

  ConstructProcessForMuonPlus();
  ConstructProcessForMuonMinus();
  ConstructProcessForTauPlus();
  ConstructProcessForTauMinus();

  if (verboseLevel>0)
    G4cout << "### Muon physics constructed." << G4endl;
}
Commit	Line	Data
5d321328	1	// $Id$
	2	// Category: physics
	3	//
	4	// Author: I. Hrivnacova
	5	//
	6	// Class TG4PhysicsConstructorMuon
	7	// -----------------------------
	8	// See the class description in the header file.
	9	// According to ExN04MuonPhysics.cc,v 1.2.2.1 2001/06/28 19:07:37 gunter Exp
	10	// GEANT4 tag Name: geant4-03-02
	11
	12	#include "TG4PhysicsConstructorMuon.h"
	13	#include "TG4ProcessControlMap.h"
	14	#include "TG4ProcessMCMap.h"
	15	#include "TG4G3Control.h"
	16
	17	#include <G4ParticleDefinition.hh>
	18	#include <G4ProcessManager.hh>
	19	#include <G4MuonPlus.hh>
	20	#include <G4MuonMinus.hh>
	21	#include <G4TauMinus.hh>
	22	#include <G4TauPlus.hh>
	23	#include <G4NeutrinoTau.hh>
	24	#include <G4AntiNeutrinoTau.hh>
	25	#include <G4NeutrinoMu.hh>
	26	#include <G4AntiNeutrinoMu.hh>
	27
	28	//TBR
	29	#include <G4MuIonisation.hh>
	30	#include <G4MuBremsstrahlung.hh>
	31	#include <G4MuPairProduction.hh>
	32	#include <G4hIonisation.hh>
	33
	34	//_____________________________________________________________________________
	35	TG4PhysicsConstructorMuon::TG4PhysicsConstructorMuon(const G4String& name)
	36	: G4VPhysicsConstructor(name)
	37	{
	38	//
	39	SetVerboseLevel(1);
	40	}
	41
	42	//_____________________________________________________________________________
	43	TG4PhysicsConstructorMuon::~TG4PhysicsConstructorMuon() {
	44	//
	45	}
	46
	47	// private methods
	48
	49	//_____________________________________________________________________________
	50	void TG4PhysicsConstructorMuon::ConstructProcessForMuonPlus()
	51	{
	52	// Constructs electromagnetic processes for mu+.
	53	// ---
	54
	55	// add processes
	56	G4ProcessManager* pManager = G4MuonPlus::MuonPlus()->GetProcessManager();
	57	pManager->AddProcess(&fMuPlusIonisation, ordInActive,2, 2);
	58	pManager->AddDiscreteProcess(&fMuPlusBremsstrahlung);
	59	pManager->AddDiscreteProcess(&fMuPlusPairProduction);
	60	pManager->AddProcess(&fMuPlusMultipleScattering);
	61
	62	// set ordering
	63	pManager->SetProcessOrdering(&fMuPlusMultipleScattering, idxAlongStep, 1);
	64	pManager->SetProcessOrdering(&fMuPlusMultipleScattering, idxPostStep, 1);
65
66	// map to G3 controls
67	TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
68	controlMap->Add(&fMuPlusIonisation, kLOSS);
69	controlMap->Add(&fMuPlusBremsstrahlung, kBREM);
70	controlMap->Add(&fMuPlusPairProduction, kPAIR);
71	controlMap->Add(&fMuPlusMultipleScattering, kMULS);
72
73
74	// map to AliMCProcess codes
75	TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
76	mcMap->Add(&fMuPlusIonisation, kPEnergyLoss);
77	mcMap->Add(&fMuPlusBremsstrahlung, kPBrem);
78	mcMap->Add(&fMuPlusPairProduction, kPPair);
79	mcMap->Add(&fMuPlusMultipleScattering, kPMultipleScattering);
80	}
81
82	//_____________________________________________________________________________
83	void TG4PhysicsConstructorMuon::ConstructProcessForMuonMinus()
84	{
85	// Constructs electromagnetic processes for mu-.
86	// ---
87
88	// add processes & set ordering
89	G4ProcessManager* pManager = G4MuonMinus::MuonMinus()->GetProcessManager();
90	pManager->AddProcess(&fMuMinusIonisation, ordInActive,2, 2);
91	pManager->AddDiscreteProcess(&fMuMinusBremsstrahlung);
92	pManager->AddDiscreteProcess(&fMuMinusPairProduction);
93	pManager->AddProcess(&fMuMinusMultipleScattering);
94
95	pManager->SetProcessOrdering(&fMuMinusMultipleScattering, idxAlongStep, 1);
96	pManager->SetProcessOrdering(&fMuMinusMultipleScattering, idxPostStep, 1);
97
98	pManager->AddRestProcess(&fMuMinusCaptureAtRest);
99
100	// map to G3 controls
101	TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
102	controlMap->Add(&fMuMinusIonisation, kLOSS);
103	controlMap->Add(&fMuMinusBremsstrahlung, kBREM);
104	controlMap->Add(&fMuMinusPairProduction, kPAIR);
105	controlMap->Add(&fMuMinusMultipleScattering, kMULS);
106	controlMap->Add(&fMuMinusCaptureAtRest, kMUNU);
107
108	// map to AliMCProcess codes
109	TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
110	mcMap->Add(&fMuMinusIonisation, kPEnergyLoss);
111	mcMap->Add(&fMuMinusBremsstrahlung, kPBrem);
112	mcMap->Add(&fMuMinusPairProduction, kPPair);
113	mcMap->Add(&fMuMinusMultipleScattering, kPMultipleScattering);
114	mcMap->Add(&fMuMinusCaptureAtRest, kPMuonNuclear);
115	}
116
117	//_____________________________________________________________________________
118	void TG4PhysicsConstructorMuon::ConstructProcessForTauPlus()
119	{
120	// Constructs electromagnetic processes for tau+.
121	// ---
122
123	// add processes
124	G4ProcessManager* pManager = G4TauPlus::TauPlus()->GetProcessManager();
125	pManager->AddProcess(&fTauPlusIonisation, ordInActive,2, 2);
126	pManager->AddProcess(&fTauPlusMultipleScattering);
127
128	// set ordering
129	pManager->SetProcessOrdering(&fTauPlusMultipleScattering, idxAlongStep, 1);
130	pManager->SetProcessOrdering(&fTauPlusMultipleScattering, idxPostStep, 1);
131
132	// map to G3 controls
133	TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
134	controlMap->Add(&fTauPlusIonisation, kLOSS);
135	controlMap->Add(&fTauPlusMultipleScattering, kMULS);
136
137
138	// map to AliMCProcess codes
139	TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
140	mcMap->Add(&fTauPlusIonisation, kPEnergyLoss);
141	mcMap->Add(&fTauPlusMultipleScattering, kPMultipleScattering);
142	}
143
144	//_____________________________________________________________________________
145	void TG4PhysicsConstructorMuon::ConstructProcessForTauMinus()
146	{
147	// Constructs electromagnetic processes for tau-.
148	// ---
149
150	// add processes
151	G4ProcessManager* pManager = G4TauMinus::TauMinus()->GetProcessManager();
152	pManager->AddProcess(&fTauMinusIonisation, ordInActive,2, 2);
153	pManager->AddProcess(&fTauMinusMultipleScattering);
154
155	// set ordering
156	pManager->SetProcessOrdering(&fTauMinusMultipleScattering, idxAlongStep, 1);
157	pManager->SetProcessOrdering(&fTauMinusMultipleScattering, idxPostStep, 1);
158
159	// map to G3 controls
160	TG4ProcessControlMap* controlMap = TG4ProcessControlMap::Instance();
161	controlMap->Add(&fTauMinusIonisation, kLOSS);
162	controlMap->Add(&fTauMinusMultipleScattering, kMULS);
163
164
165	// map to AliMCProcess codes
166	TG4ProcessMCMap* mcMap = TG4ProcessMCMap::Instance();
167	mcMap->Add(&fTauMinusIonisation, kPEnergyLoss);
168	mcMap->Add(&fTauMinusMultipleScattering, kPMultipleScattering);
169	}
170
171
172	// protected methods
173
174	//_____________________________________________________________________________
175	void TG4PhysicsConstructorMuon::ConstructParticle()
176	{
177	// Instantiates particles.
178	// ---
179
180	// Mu
181	G4MuonPlus::MuonPlusDefinition();
182	G4MuonMinus::MuonMinusDefinition();
183	G4NeutrinoMu::NeutrinoMuDefinition();
184	G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
185
186	// Tau
187	G4TauMinus::TauMinusDefinition();
188	G4TauPlus::TauPlusDefinition();
189	G4NeutrinoTau::NeutrinoTauDefinition();
190	G4AntiNeutrinoTau::AntiNeutrinoTauDefinition();
191	}
192
193	//_____________________________________________________________________________
194	void TG4PhysicsConstructorMuon::ConstructProcess()
195	{
196	// Constructs electromagnetic processes for muons.
197	// ---
198
199	ConstructProcessForMuonPlus();
200	ConstructProcessForMuonMinus();
201	ConstructProcessForTauPlus();
202	ConstructProcessForTauMinus();
203
204	if (verboseLevel>0)
205	G4cout << "### Muon physics constructed." << G4endl;
206	}