Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
G4HadronicProcess Class Reference

#include <G4HadronicProcess.hh>

+ Inheritance diagram for G4HadronicProcess:

Public Member Functions

 G4HadronicProcess (const G4String &processName="Hadronic", G4ProcessType procType=fHadronic)
 
 G4HadronicProcess (const G4String &processName, G4HadronicProcessType subType)
 
 ~G4HadronicProcess () override
 
void RegisterMe (G4HadronicInteraction *a)
 
G4double GetElementCrossSection (const G4DynamicParticle *part, const G4Element *elm, const G4Material *mat=nullptr)
 
G4double GetMicroscopicCrossSection (const G4DynamicParticle *part, const G4Element *elm, const G4Material *mat=nullptr)
 
void StartTracking (G4Track *track) override
 
G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double, G4ForceCondition *) override
 
G4VParticleChange * PostStepDoIt (const G4Track &aTrack, const G4Step &aStep) override
 
void PreparePhysicsTable (const G4ParticleDefinition &) override
 
void BuildPhysicsTable (const G4ParticleDefinition &) override
 
void DumpPhysicsTable (const G4ParticleDefinition &p)
 
void AddDataSet (G4VCrossSectionDataSet *aDataSet)
 
std::vector< G4HadronicInteraction * > & GetHadronicInteractionList ()
 
G4HadronicInteraction * GetHadronicModel (const G4String &)
 
G4HadronicInteraction * GetHadronicInteraction () const
 
G4double GetMeanFreePath (const G4Track &aTrack, G4double, G4ForceCondition *) override
 
const G4Nucleus * GetTargetNucleus () const
 
G4Nucleus * GetTargetNucleusPointer ()
 
const G4Isotope * GetTargetIsotope ()
 
G4double ComputeCrossSection (const G4ParticleDefinition *, const G4Material *, const G4double kinEnergy)
 
G4HadXSType CrossSectionType () const
 
void SetCrossSectionType (G4HadXSType val)
 
void ProcessDescription (std::ostream &outFile) const override
 
void BiasCrossSectionByFactor (G4double aScale)
 
void MultiplyCrossSectionBy (G4double factor)
 
G4double CrossSectionFactor () const
 
void SetIntegral (G4bool val)
 
void SetEpReportLevel (G4int level)
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
std::pair< G4double, G4double > GetEnergyMomentumCheckLevels () const
 
G4CrossSectionDataStore * GetCrossSectionDataStore ()
 
std::vector< G4TwoPeaksHadXS * > * TwoPeaksXS () const
 
std::vector< G4double > * EnergyOfCrossSectionMax () const
 
G4HadronicProcess & operator= (const G4HadronicProcess &right)=delete
 
 G4HadronicProcess (const G4HadronicProcess &)=delete
 
- Public Member Functions inherited from G4VDiscreteProcess
 G4VDiscreteProcess (const G4String &aName, G4ProcessType aType=fNotDefined)
 
 G4VDiscreteProcess (G4VDiscreteProcess &)
 
virtual ~G4VDiscreteProcess ()
 
G4VDiscreteProcess & operator= (const G4VDiscreteProcess &)=delete
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &, const G4Step &)
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4VParticleChange * AlongStepDoIt (const G4Track &, const G4Step &)
 
virtual G4double GetCrossSection (const G4double, const G4MaterialCutsCouple *)
 
virtual G4double MinPrimaryEnergy (const G4ParticleDefinition *, const G4Material *)
 
- Public Member Functions inherited from G4VProcess
 G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
 
 G4VProcess (const G4VProcess &right)
 
virtual ~G4VProcess ()
 
G4VProcess & operator= (const G4VProcess &)=delete
 
G4bool operator== (const G4VProcess &right) const
 
G4bool operator!= (const G4VProcess &right) const
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChange * AlongStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)=0
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)=0
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)=0
 
G4double GetCurrentInteractionLength () const
 
void SetPILfactor (G4double value)
 
G4double GetPILfactor () const
 
G4double AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
 
G4double AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
 
G4double PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4bool IsApplicable (const G4ParticleDefinition &)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void PreparePhysicsTable (const G4ParticleDefinition &)
 
virtual G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
virtual G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
const G4String & GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
 
const G4String & GetProcessName () const
 
G4ProcessType GetProcessType () const
 
void SetProcessType (G4ProcessType)
 
G4int GetProcessSubType () const
 
void SetProcessSubType (G4int)
 
virtual const G4VProcess * GetCreatorProcess () const
 
virtual void StartTracking (G4Track *)
 
virtual void EndTracking ()
 
virtual void SetProcessManager (const G4ProcessManager *)
 
virtual const G4ProcessManager * GetProcessManager ()
 
virtual void ResetNumberOfInteractionLengthLeft ()
 
G4double GetNumberOfInteractionLengthLeft () const
 
G4double GetTotalNumberOfInteractionLengthTraversed () const
 
G4bool isAtRestDoItIsEnabled () const
 
G4bool isAlongStepDoItIsEnabled () const
 
G4bool isPostStepDoItIsEnabled () const
 
virtual void DumpInfo () const
 
virtual void ProcessDescription (std::ostream &outfile) const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcess * GetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Protected Member Functions

G4HadronicInteraction * ChooseHadronicInteraction (const G4HadProjectile &aHadProjectile, G4Nucleus &aTargetNucleus, const G4Material *aMaterial, const G4Element *anElement)
 
G4double GetLastCrossSection ()
 
void FillResult (G4HadFinalState *aR, const G4Track &aT)
 
void DumpState (const G4Track &, const G4String &, G4ExceptionDescription &)
 
G4HadFinalState * CheckResult (const G4HadProjectile &thePro, const G4Nucleus &targetNucleus, G4HadFinalState *result)
 
void CheckEnergyMomentumConservation (const G4Track &, const G4Nucleus &)
 
virtual G4double GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Protected Attributes

G4HadProjectile thePro
 
G4ParticleChange * theTotalResult
 
G4CrossSectionDataStore * theCrossSectionDataStore
 
G4double fWeight = 1.0
 
G4double aScaleFactor = 1.0
 
G4double theLastCrossSection = 0.0
 
G4double mfpKinEnergy = DBL_MAX
 
G4long epReportLevel = 0
 
G4HadXSType fXSType = fHadNoIntegral
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManager * aProcessManager = nullptr
 
G4VParticleChange * pParticleChange = nullptr
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft = -1.0
 
G4double currentInteractionLength = -1.0
 
G4double theInitialNumberOfInteractionLength = -1.0
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType = fNotDefined
 
G4int theProcessSubType = -1
 
G4double thePILfactor = 1.0
 
G4int verboseLevel = 0
 
G4bool enableAtRestDoIt = true
 
G4bool enableAlongStepDoIt = true
 
G4bool enablePostStepDoIt = true
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 

Detailed Description

Definition at line 73 of file G4HadronicProcess.hh.

Constructor & Destructor Documentation

◆ G4HadronicProcess() [1/3]

G4HadronicProcess::G4HadronicProcess ( const G4String &  processName = "Hadronic",
G4ProcessType  procType = fHadronic 
)

Definition at line 87 of file G4HadronicProcess.cc.

89 : G4VDiscreteProcess(processName, procType)
90{
91 SetProcessSubType(fHadronInelastic); // Default unless subclass changes
92 InitialiseLocal();
93}
fHadronInelastic
@ fHadronInelastic
Definition: G4HadronicProcessType.hh:45
G4VProcess::SetProcessSubType
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:410

◆ G4HadronicProcess() [2/3]

G4HadronicProcess::G4HadronicProcess ( const G4String &  processName,
G4HadronicProcessType  subType 
)

Definition at line 95 of file G4HadronicProcess.cc.

97 : G4VDiscreteProcess(processName, fHadronic)
98{
99 SetProcessSubType(aHadSubType);
100 InitialiseLocal();
101}
fHadronic
@ fHadronic
Definition: G4ProcessType.hh:43

◆ ~G4HadronicProcess()

G4HadronicProcess::~G4HadronicProcess ( )
override

Definition at line 103 of file G4HadronicProcess.cc.

104{
105 theProcessStore->DeRegister(this);
106 delete theTotalResult;
107 delete theCrossSectionDataStore;
108 if(isMaster) {
109 delete fXSpeaks;
110 delete theEnergyOfCrossSectionMax;
111 }
112}
G4HadronicProcessStore::DeRegister
void DeRegister(G4HadronicProcess *)
Definition: G4HadronicProcessStore.cc:467
G4HadronicProcess::theTotalResult
G4ParticleChange * theTotalResult
Definition: G4HadronicProcess.hh:220
G4HadronicProcess::theCrossSectionDataStore
G4CrossSectionDataStore * theCrossSectionDataStore
Definition: G4HadronicProcess.hh:221

◆ G4HadronicProcess() [3/3]

G4HadronicProcess::G4HadronicProcess ( const G4HadronicProcess &  )
delete

Member Function Documentation

◆ AddDataSet()

void G4HadronicProcess::AddDataSet ( G4VCrossSectionDataSet *  aDataSet)

Definition at line 890 of file G4HadronicProcess.cc.

891{
892 theCrossSectionDataStore->AddDataSet(aDataSet);
893}
G4CrossSectionDataStore::AddDataSet
void AddDataSet(G4VCrossSectionDataSet *)
Definition: G4CrossSectionDataStore.cc:398

Referenced by G4HadProcesses::AddCaptureCrossSection(), G4HadProcesses::AddElasticCrossSection(), G4HadProcesses::AddFissionCrossSection(), G4HadProcesses::AddInelasticCrossSection(), G4HadronElasticPhysics::AddXSection(), G4BertiniElectroNuclearBuilder::Build(), G4LENDBertiniGammaElectroNuclearBuilder::Build(), G4NeutronLENDBuilder::Build(), G4NeutronPHPBuilder::Build(), G4AlphaPHPBuilder::Build(), G4BertiniKaonBuilder::Build(), G4BertiniPiKBuilder::Build(), G4BertiniPionBuilder::Build(), G4BinaryPiKBuilder::Build(), G4BinaryPionBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4FTFBinaryNeutronBuilder::Build(), G4FTFBinaryPiKBuilder::Build(), G4FTFBinaryPionBuilder::Build(), G4FTFBinaryProtonBuilder::Build(), G4FTFPAntiBarionBuilder::Build(), G4FTFPNeutronBuilder::Build(), G4FTFPPiKBuilder::Build(), G4FTFPPionBuilder::Build(), G4FTFPProtonBuilder::Build(), G4He3PHPBuilder::Build(), G4HyperonFTFPBuilder::Build(), G4HyperonQGSPBuilder::Build(), G4INCLXXNeutronBuilder::Build(), G4INCLXXPionBuilder::Build(), G4INCLXXProtonBuilder::Build(), G4PrecoNeutronBuilder::Build(), G4PrecoProtonBuilder::Build(), G4ProtonPHPBuilder::Build(), G4QGSBinaryNeutronBuilder::Build(), G4QGSBinaryPiKBuilder::Build(), G4QGSBinaryPionBuilder::Build(), G4QGSBinaryProtonBuilder::Build(), G4QGSPAntiBarionBuilder::Build(), G4QGSPLundStrFragmProtonBuilder::Build(), G4QGSPNeutronBuilder::Build(), G4QGSPPiKBuilder::Build(), G4QGSPPionBuilder::Build(), G4TritonPHPBuilder::Build(), G4QGSPProtonBuilder::Build(), G4HadProcesses::BuildNeutronElastic(), G4HadProcesses::BuildNeutronInelasticAndCapture(), LBE::ConstructHad(), G4NeutronCrossSectionXS::ConstructProcess(), G4HadronElasticPhysicsHP::ConstructProcess(), G4HadronElasticPhysicsLEND::ConstructProcess(), G4HadronElasticPhysicsPHP::ConstructProcess(), G4EmExtraPhysics::ConstructProcess(), G4ChargeExchangePhysics::ConstructProcess(), G4HadronDElasticPhysics::ConstructProcess(), G4HadronElasticPhysics::ConstructProcess(), G4HadronHElasticPhysics::ConstructProcess(), G4IonElasticPhysics::ConstructProcess(), G4ThermalNeutrons::ConstructProcess(), G4HadronInelasticQBBC::ConstructProcess(), G4HadronPhysicsFTFQGSP_BERT::ConstructProcess(), G4ChargeExchangeProcess::G4ChargeExchangeProcess(), G4MuonNuclearProcess::G4MuonNuclearProcess(), G4NeutronCaptureProcess::G4NeutronCaptureProcess(), G4NeutronFissionProcess::G4NeutronFissionProcess(), G4HadronPhysicsFTFP_BERT::Neutron(), G4HadronPhysicsQGSP_BERT::Neutron(), G4HadronPhysicsQGSP_BIC::Neutron(), G4HadronPhysicsFTF_BIC::Neutron(), G4HadronPhysicsQGS_BIC::Neutron(), G4NeutronGeneralProcess::SetCaptureProcess(), G4NeutronGeneralProcess::SetElasticProcess(), and G4NeutronGeneralProcess::SetInelasticProcess().

◆ BiasCrossSectionByFactor()

void G4HadronicProcess::BiasCrossSectionByFactor ( G4double  aScale)

Definition at line 612 of file G4HadronicProcess.cc.

613{
614 if (aScale <= 0.0) {
615 G4ExceptionDescription ed;
616 ed << " Wrong biasing factor " << aScale << " for " << GetProcessName();
617 G4Exception("G4HadronicProcess::BiasCrossSectionByFactor", "had010",
618 JustWarning, ed, "Cross-section bias is ignored");
619 } else {
620 aScaleFactor = aScale;
621 }
622}
JustWarning
@ JustWarning
Definition: G4ExceptionSeverity.hh:65
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:59
G4ExceptionDescription
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
G4VProcess::GetProcessName
const G4String & GetProcessName() const
Definition: G4VProcess.hh:386

Referenced by MultiplyCrossSectionBy().

◆ BuildPhysicsTable()

void G4HadronicProcess::BuildPhysicsTable ( const G4ParticleDefinition &  p)
overridevirtual

Reimplemented from G4VProcess.

Reimplemented in G4HadronStoppingProcess, G4NeutronGeneralProcess, and G4ChargeExchangeProcess.

Definition at line 178 of file G4HadronicProcess.cc.

179{
180 if(firstParticle != &p) { return; }
181
182 theCrossSectionDataStore->BuildPhysicsTable(p);
183 theEnergyRangeManager.BuildPhysicsTable(p);
184 G4HadronicParameters* param = G4HadronicParameters::Instance();
185
186 G4int subtype = GetProcessSubType();
187 if(useIntegralXS) {
188 if(subtype == fHadronInelastic) {
189 useIntegralXS = param->EnableIntegralInelasticXS();
190 } else if(subtype == fHadronElastic) {
191 useIntegralXS = param->EnableIntegralElasticXS();
192 }
193 }
194 fXSType = fHadNoIntegral;
195
196 // check particle for integral method
197 if(isMaster) {
198 G4double charge = p.GetPDGCharge()/eplus;
199 G4bool isLepton = (p.GetLeptonNumber() != 0);
200 G4bool ok = (p.GetAtomicNumber() != 0 || p.GetPDGMass() < GeV);
201
202 // select cross section shape
203 if(charge != 0.0 && useIntegralXS && !isLepton && ok) {
204 G4double tmax = param->GetMaxEnergy();
205 fXSType = (charge > 0.0) ? fHadIncreasing : fHadDecreasing;
206 currentParticle = firstParticle;
207 // initialisation in the master thread
208 G4int pdg = p.GetPDGEncoding();
209 if(std::abs(pdg) == 211) {
210 fXSType = fHadTwoPeaks;
211 } else if(pdg == 321) {
212 fXSType = fHadOnePeak;
213 } else if(pdg == 2212) {
214 fXSType = fHadTwoPeaks;
215 }
216 delete theEnergyOfCrossSectionMax;
217 theEnergyOfCrossSectionMax = nullptr;
218 if(fXSType == fHadTwoPeaks) {
219 delete fXSpeaks;
220 fXSpeaks =
221 G4HadXSHelper::FillPeaksStructure(this, &p, minKinEnergy, tmax);
222 if(nullptr == fXSpeaks) {
223 fXSType = fHadOnePeak;
224 }
225 }
226 if(fXSType == fHadOnePeak) {
227 theEnergyOfCrossSectionMax =
228 G4HadXSHelper::FindCrossSectionMax(this, &p, minKinEnergy, tmax);
229 if(nullptr == theEnergyOfCrossSectionMax) {
230 fXSType = fHadIncreasing;
231 }
232 }
233 }
234 } else {
235 if(nullptr == masterProcess) {
236 masterProcess =
237 dynamic_cast<const G4HadronicProcess*>(GetMasterProcess());
238 }
239 if(nullptr == masterProcess) {
240 G4cout << "G4HadronicProcess::BuildPhysicsTable: for "
241 << GetProcessName() << " and " << p.GetParticleName()
242 << " fail due to undefined pointer to the master process"
243 << G4endl;
244 } else {
245 // initialisation in worker threads
246 fXSType = masterProcess->CrossSectionType();
247 fXSpeaks = masterProcess->TwoPeaksXS();
248 theEnergyOfCrossSectionMax =
249 masterProcess->EnergyOfCrossSectionMax();
250 }
251 }
252 if(isMaster && 1 < param->GetVerboseLevel()) {
253 G4cout << "G4HadronicProcess::BuildPhysicsTable: for "
254 << GetProcessName() << " and " << p.GetParticleName()
255 << " typeXS=" << fXSType << G4endl;
256 }
257 G4HadronicProcessStore::Instance()->PrintInfo(&p);
258}
fHadTwoPeaks
@ fHadTwoPeaks
Definition: G4HadXSTypes.hh:49
fHadIncreasing
@ fHadIncreasing
Definition: G4HadXSTypes.hh:46
fHadDecreasing
@ fHadDecreasing
Definition: G4HadXSTypes.hh:47
fHadNoIntegral
@ fHadNoIntegral
Definition: G4HadXSTypes.hh:45
fHadOnePeak
@ fHadOnePeak
Definition: G4HadXSTypes.hh:48
fHadronElastic
@ fHadronElastic
Definition: G4HadronicProcessType.hh:43
G4double
double G4double
Definition: G4Types.hh:83
G4bool
bool G4bool
Definition: G4Types.hh:86
G4int
int G4int
Definition: G4Types.hh:85
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4CrossSectionDataStore::BuildPhysicsTable
void BuildPhysicsTable(const G4ParticleDefinition &)
Definition: G4CrossSectionDataStore.cc:272
G4EnergyRangeManager::BuildPhysicsTable
void BuildPhysicsTable(const G4ParticleDefinition &)
Definition: G4EnergyRangeManager.cc:178
G4HadXSHelper::FillPeaksStructure
static std::vector< G4TwoPeaksHadXS * > * FillPeaksStructure(G4HadronicProcess *, const G4ParticleDefinition *, const G4double tmin, const G4double tmax)
Definition: G4HadXSHelper.cc:97
G4HadXSHelper::FindCrossSectionMax
static std::vector< G4double > * FindCrossSectionMax(G4HadronicProcess *, const G4ParticleDefinition *, const G4double tmin, const G4double tmax)
Definition: G4HadXSHelper.cc:43
G4HadronicParameters::EnableIntegralInelasticXS
G4bool EnableIntegralInelasticXS() const
Definition: G4HadronicParameters.hh:249
G4HadronicParameters::Instance
static G4HadronicParameters * Instance()
Definition: G4HadronicParameters.cc:54
G4HadronicParameters::EnableIntegralElasticXS
G4bool EnableIntegralElasticXS() const
Definition: G4HadronicParameters.hh:253
G4HadronicParameters::GetMaxEnergy
G4double GetMaxEnergy() const
Definition: G4HadronicParameters.hh:178
G4HadronicProcessStore::Instance
static G4HadronicProcessStore * Instance()
Definition: G4HadronicProcessStore.cc:68
G4HadronicProcessStore::PrintInfo
void PrintInfo(const G4ParticleDefinition *)
Definition: G4HadronicProcessStore.cc:562
G4HadronicProcess::TwoPeaksXS
std::vector< G4TwoPeaksHadXS * > * TwoPeaksXS() const
Definition: G4HadronicProcess.hh:343
G4HadronicProcess::CrossSectionType
G4HadXSType CrossSectionType() const
Definition: G4HadronicProcess.hh:295
G4HadronicProcess::EnergyOfCrossSectionMax
std::vector< G4double > * EnergyOfCrossSectionMax() const
Definition: G4HadronicProcess.hh:349
G4ParticleDefinition::GetAtomicNumber
G4int GetAtomicNumber() const
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:113
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:109
G4VProcess::GetMasterProcess
const G4VProcess * GetMasterProcess() const
Definition: G4VProcess.hh:522
G4VProcess::GetProcessSubType
G4int GetProcessSubType() const
Definition: G4VProcess.hh:404

Referenced by G4GammaGeneralProcess::BuildPhysicsTable(), G4NeutronGeneralProcess::BuildPhysicsTable(), and G4ChargeExchangeProcess::BuildPhysicsTable().

◆ CheckEnergyMomentumConservation()

void G4HadronicProcess::CheckEnergyMomentumConservation ( const G4Track &  aTrack,
const G4Nucleus &  aNucleus 
)
protected

Definition at line 715 of file G4HadronicProcess.cc.

717{
718 G4int target_A=aNucleus.GetA_asInt();
719 G4int target_Z=aNucleus.GetZ_asInt();
720 G4double targetMass = G4NucleiProperties::GetNuclearMass(target_A,target_Z);
721 G4LorentzVector target4mom(0, 0, 0, targetMass
722 + nICelectrons*CLHEP::electron_mass_c2);
723
724 G4LorentzVector projectile4mom = aTrack.GetDynamicParticle()->Get4Momentum();
725 G4int track_A = aTrack.GetDefinition()->GetBaryonNumber();
726 G4int track_Z = G4lrint(aTrack.GetDefinition()->GetPDGCharge());
727
728 G4int initial_A = target_A + track_A;
729 G4int initial_Z = target_Z + track_Z - nICelectrons;
730
731 G4LorentzVector initial4mom = projectile4mom + target4mom;
732
733 // Compute final-state momentum for scattering and "do nothing" results
734 G4LorentzVector final4mom;
735 G4int final_A(0), final_Z(0);
736
737 G4int nSec = theTotalResult->GetNumberOfSecondaries();
738 if (theTotalResult->GetTrackStatus() != fStopAndKill) { // If it is Alive
739 // Either interaction didn't complete, returned "do nothing" state
740 // or the primary survived the interaction (e.g. electro-nucleus )
741
742 // Interaction didn't complete, returned "do nothing" state
743 // - or suppressed recoil (e.g. Neutron elastic )
744 final4mom = initial4mom;
745 final_A = initial_A;
746 final_Z = initial_Z;
747 if (nSec > 0) {
748 // The primary remains in final state (e.g. electro-nucleus )
749 // Use the final energy / momentum
750 const G4ThreeVector& v = *theTotalResult->GetMomentumDirection();
751 G4double ekin = theTotalResult->GetEnergy();
752 G4double mass = aTrack.GetDefinition()->GetPDGMass();
753 G4double ptot = std::sqrt(ekin*(ekin + 2*mass));
754 final4mom.set(ptot*v.x(), ptot*v.y(), ptot*v.z(), mass + ekin);
755 final_A = track_A;
756 final_Z = track_Z;
757 // Expect that the target nucleus will have interacted,
758 // and its products, including recoil, will be included in secondaries.
759 }
760 }
761 if( nSec > 0 ) {
762 G4Track* sec;
763
764 for (G4int i = 0; i < nSec; i++) {
765 sec = theTotalResult->GetSecondary(i);
766 final4mom += sec->GetDynamicParticle()->Get4Momentum();
767 final_A += sec->GetDefinition()->GetBaryonNumber();
768 final_Z += G4lrint(sec->GetDefinition()->GetPDGCharge());
769 }
770 }
771
772 // Get level-checking information (used to cut-off relative checks)
773 G4String processName = GetProcessName();
774 G4HadronicInteraction* theModel = GetHadronicInteraction();
775 G4String modelName("none");
776 if (theModel) modelName = theModel->GetModelName();
777 std::pair<G4double, G4double> checkLevels = epCheckLevels;
778 if (!levelsSetByProcess) {
779 if (theModel) checkLevels = theModel->GetEnergyMomentumCheckLevels();
780 checkLevels.first= std::min(checkLevels.first, epCheckLevels.first);
781 checkLevels.second=std::min(checkLevels.second, epCheckLevels.second);
782 }
783
784 // Compute absolute total-energy difference, and relative kinetic-energy
785 G4bool checkRelative = (aTrack.GetKineticEnergy() > checkLevels.second);
786
787 G4LorentzVector diff = initial4mom - final4mom;
788 G4double absolute = diff.e();
789 G4double relative = checkRelative ? absolute/aTrack.GetKineticEnergy() : 0.;
790
791 G4double absolute_mom = diff.vect().mag();
792 G4double relative_mom = checkRelative ? absolute_mom/aTrack.GetMomentum().mag() : 0.;
793
794 // Evaluate relative and absolute conservation
795 G4bool relPass = true;
796 G4String relResult = "pass";
797 if ( std::abs(relative) > checkLevels.first
798 || std::abs(relative_mom) > checkLevels.first) {
799 relPass = false;
800 relResult = checkRelative ? "fail" : "N/A";
801 }
802
803 G4bool absPass = true;
804 G4String absResult = "pass";
805 if ( std::abs(absolute) > checkLevels.second
806 || std::abs(absolute_mom) > checkLevels.second ) {
807 absPass = false ;
808 absResult = "fail";
809 }
810
811 G4bool chargePass = true;
812 G4String chargeResult = "pass";
813 if ( (initial_A-final_A)!=0
814 || (initial_Z-final_Z)!=0 ) {
815 chargePass = checkLevels.second < DBL_MAX ? false : true;
816 chargeResult = "fail";
817 }
818
819 G4bool conservationPass = (relPass || absPass) && chargePass;
820
821 std::stringstream Myout;
822 G4bool Myout_notempty(false);
823 // Options for level of reporting detail:
824 // 0. off
825 // 1. report only when E/p not conserved
826 // 2. report regardless of E/p conservation
827 // 3. report only when E/p not conserved, with model names, process names, and limits
828 // 4. report regardless of E/p conservation, with model names, process names, and limits
829 // negative -1.., as above, but send output to stderr
830
831 if( std::abs(epReportLevel) == 4
832 || ( std::abs(epReportLevel) == 3 && ! conservationPass ) ){
833 Myout << " Process: " << processName << " , Model: " << modelName << G4endl;
834 Myout << " Primary: " << aTrack.GetParticleDefinition()->GetParticleName()
835 << " (" << aTrack.GetParticleDefinition()->GetPDGEncoding() << "),"
836 << " E= " << aTrack.GetDynamicParticle()->Get4Momentum().e()
837 << ", target nucleus (" << aNucleus.GetZ_asInt() << ","
838 << aNucleus.GetA_asInt() << ")" << G4endl;
839 Myout_notempty=true;
840 }
841 if ( std::abs(epReportLevel) == 4
842 || std::abs(epReportLevel) == 2
843 || ! conservationPass ){
844
845 Myout << " "<< relResult <<" relative, limit " << checkLevels.first << ", values E/T(0) = "
846 << relative << " p/p(0)= " << relative_mom << G4endl;
847 Myout << " "<< absResult << " absolute, limit (MeV) " << checkLevels.second/MeV << ", values E / p (MeV) = "
848 << absolute/MeV << " / " << absolute_mom/MeV << " 3mom: " << (diff.vect())*1./MeV << G4endl;
849 Myout << " "<< chargeResult << " charge/baryon number balance " << (initial_Z-final_Z) << " / " << (initial_A-final_A) << " "<< G4endl;
850 Myout_notempty=true;
851
852 }
853 Myout.flush();
854 if ( Myout_notempty ) {
855 if (epReportLevel > 0) G4cout << Myout.str()<< G4endl;
856 else if (epReportLevel < 0) G4cerr << Myout.str()<< G4endl;
857 }
858}
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:46
G4cerr
G4GLOB_DLL std::ostream G4cerr
CLHEP::Hep3Vector::z
double z() const
CLHEP::Hep3Vector::x
double x() const
CLHEP::Hep3Vector::y
double y() const
CLHEP::Hep3Vector::mag
double mag() const
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
CLHEP::HepLorentzVector::set
void set(double x, double y, double z, double t)
G4DynamicParticle::Get4Momentum
G4LorentzVector Get4Momentum() const
G4HadronicInteraction::GetEnergyMomentumCheckLevels
virtual std::pair< G4double, G4double > GetEnergyMomentumCheckLevels() const
Definition: G4HadronicInteraction.cc:217
G4HadronicInteraction::GetModelName
const G4String & GetModelName() const
Definition: G4HadronicInteraction.hh:115
G4HadronicProcess::GetHadronicInteraction
G4HadronicInteraction * GetHadronicInteraction() const
Definition: G4HadronicProcess.hh:278
G4NucleiProperties::GetNuclearMass
static G4double GetNuclearMass(const G4double A, const G4double Z)
Definition: G4NucleiProperties.cc:50
G4Nucleus::GetA_asInt
G4int GetA_asInt() const
Definition: G4Nucleus.hh:99
G4Nucleus::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Nucleus.hh:105
G4ParticleChange::GetEnergy
G4double GetEnergy() const
G4ParticleChange::GetMomentumDirection
const G4ThreeVector * GetMomentumDirection() const
G4Track::GetParticleDefinition
const G4ParticleDefinition * GetParticleDefinition() const
G4Track::GetMomentum
G4ThreeVector GetMomentum() const
G4Track::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4Track::GetKineticEnergy
G4double GetKineticEnergy() const
G4VParticleChange::GetNumberOfSecondaries
G4int GetNumberOfSecondaries() const
G4VParticleChange::GetSecondary
G4Track * GetSecondary(G4int anIndex) const
G4VParticleChange::GetTrackStatus
G4TrackStatus GetTrackStatus() const
G4lrint
int G4lrint(double ad)
Definition: templates.hh:134
DBL_MAX
#define DBL_MAX
Definition: templates.hh:62

Referenced by G4HadronStoppingProcess::AtRestDoIt(), PostStepDoIt(), and G4HadronElasticProcess::PostStepDoIt().

◆ CheckResult()

G4HadFinalState * G4HadronicProcess::CheckResult ( const G4HadProjectile &  thePro,
const G4Nucleus &  targetNucleus,
G4HadFinalState *  result 
)
protected

Definition at line 624 of file G4HadronicProcess.cc.

627{
628 // check for catastrophic energy non-conservation
629 // to re-sample the interaction
630 G4HadronicInteraction * theModel = GetHadronicInteraction();
631 G4double nuclearMass(0);
632 if (nullptr != theModel) {
633
634 // Compute final-state total energy
635 G4double finalE(0.);
636 G4int nSec = (G4int)result->GetNumberOfSecondaries();
637
638 nuclearMass = G4NucleiProperties::GetNuclearMass(aNucleus.GetA_asInt(),
639 aNucleus.GetZ_asInt());
640 if (result->GetStatusChange() != stopAndKill) {
641 // Interaction didn't complete, returned "do nothing" state
642 // and reset nucleus or the primary survived the interaction
643 // (e.g. electro-nuclear ) => keep nucleus
644 finalE=result->GetLocalEnergyDeposit() +
645 aPro.GetDefinition()->GetPDGMass() + result->GetEnergyChange();
646 if( nSec == 0 ){
647 // Since there are no secondaries, there is no recoil nucleus.
648 // To check energy balance we must neglect the initial nucleus too.
649 nuclearMass=0.0;
650 }
651 }
652 for (G4int i = 0; i < nSec; ++i) {
653 G4DynamicParticle *pdyn=result->GetSecondary(i)->GetParticle();
654 finalE += pdyn->GetTotalEnergy();
655 G4double mass_pdg=pdyn->GetDefinition()->GetPDGMass();
656 G4double mass_dyn=pdyn->GetMass();
657 if ( std::abs(mass_pdg - mass_dyn) > 0.1*mass_pdg + 1.*MeV ) {
658 // If it is shortlived, then a difference less than 3 times the width is acceptable
659 if ( pdyn->GetDefinition()->IsShortLived() &&
660 std::abs(mass_pdg - mass_dyn) < 3.0*pdyn->GetDefinition()->GetPDGWidth() ) {
661 continue;
662 }
663 result->Clear();
664 result = nullptr;
665 G4ExceptionDescription desc;
666 desc << "Warning: Secondary with off-shell dynamic mass detected: "
667 << G4endl
668 << " " << pdyn->GetDefinition()->GetParticleName()
669 << ", PDG mass: " << mass_pdg << ", dynamic mass: "
670 << mass_dyn << G4endl
671 << (epReportLevel<0 ? "abort the event"
672 : "re-sample the interaction") << G4endl
673 << " Process / Model: " << GetProcessName()<< " / "
674 << theModel->GetModelName() << G4endl
675 << " Primary: " << aPro.GetDefinition()->GetParticleName()
676 << " (" << aPro.GetDefinition()->GetPDGEncoding() << "), "
677 << " E= " << aPro.Get4Momentum().e()
678 << ", target nucleus (" << aNucleus.GetZ_asInt() << ", "
679 << aNucleus.GetA_asInt() << ")" << G4endl;
680 G4Exception("G4HadronicProcess:CheckResult()", "had012",
681 epReportLevel<0 ? EventMustBeAborted : JustWarning,desc);
682 // must return here.....
683 return result;
684 }
685 }
686 G4double deltaE= nuclearMass + aPro.GetTotalEnergy() - finalE;
687
688 std::pair<G4double, G4double> checkLevels =
689 theModel->GetFatalEnergyCheckLevels(); // (relative, absolute)
690 if (std::abs(deltaE) > checkLevels.second &&
691 std::abs(deltaE) > checkLevels.first*aPro.GetKineticEnergy()){
692 // do not delete result, this is a pointer to a data member;
693 result->Clear();
694 result = nullptr;
695 G4ExceptionDescription desc;
696 desc << "Warning: Bad energy non-conservation detected, will "
697 << (epReportLevel<0 ? "abort the event"
698 : "re-sample the interaction") << G4endl
699 << " Process / Model: " << GetProcessName()<< " / "
700 << theModel->GetModelName() << G4endl
701 << " Primary: " << aPro.GetDefinition()->GetParticleName()
702 << " (" << aPro.GetDefinition()->GetPDGEncoding() << "), "
703 << " E= " << aPro.Get4Momentum().e()
704 << ", target nucleus (" << aNucleus.GetZ_asInt() << ", "
705 << aNucleus.GetA_asInt() << ")" << G4endl
706 << " E(initial - final) = " << deltaE << " MeV." << G4endl;
707 G4Exception("G4HadronicProcess:CheckResult()", "had012",
708 epReportLevel<0 ? EventMustBeAborted : JustWarning,desc);
709 }
710 }
711 return result;
712}
EventMustBeAborted
@ EventMustBeAborted
Definition: G4ExceptionSeverity.hh:64
stopAndKill
@ stopAndKill
Definition: G4HadFinalState.hh:42
G4DynamicParticle::GetMass
G4double GetMass() const
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4DynamicParticle::GetTotalEnergy
G4double GetTotalEnergy() const
G4HadFinalState::GetEnergyChange
G4double GetEnergyChange() const
Definition: G4HadFinalState.hh:54
G4HadFinalState::GetStatusChange
G4HadFinalStateStatus GetStatusChange() const
Definition: G4HadFinalState.hh:69
G4HadFinalState::GetLocalEnergyDeposit
G4double GetLocalEnergyDeposit() const
Definition: G4HadFinalState.hh:84
G4HadFinalState::GetNumberOfSecondaries
std::size_t GetNumberOfSecondaries() const
Definition: G4HadFinalState.hh:51
G4HadFinalState::GetSecondary
G4HadSecondary * GetSecondary(size_t i)
Definition: G4HadFinalState.cc:79
G4HadSecondary::GetParticle
G4DynamicParticle * GetParticle()
Definition: G4HadSecondary.hh:43
G4HadronicInteraction::GetFatalEnergyCheckLevels
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels() const
Definition: G4HadronicInteraction.cc:210
G4ParticleDefinition::GetPDGWidth
G4double GetPDGWidth() const
Definition: G4ParticleDefinition.hh:112

Referenced by G4HadronStoppingProcess::AtRestDoIt(), PostStepDoIt(), and G4HadronElasticProcess::PostStepDoIt().

◆ ChooseHadronicInteraction()

G4HadronicInteraction * G4HadronicProcess::ChooseHadronicInteraction ( const G4HadProjectile &  aHadProjectile,
G4Nucleus &  aTargetNucleus,
const G4Material *  aMaterial,
const G4Element *  anElement 
)
inlineprotected

Definition at line 354 of file G4HadronicProcess.hh.

359{
360 return theEnergyRangeManager.GetHadronicInteraction(aHadProjectile,
361 aTargetNucleus,
362 aMaterial,anElement);
363}
G4EnergyRangeManager::GetHadronicInteraction
G4HadronicInteraction * GetHadronicInteraction(const G4HadProjectile &aHadProjectile, G4Nucleus &aTargetNucleus, const G4Material *aMaterial, const G4Element *anElement) const
Definition: G4EnergyRangeManager.cc:59

Referenced by G4HadronStoppingProcess::AtRestDoIt(), PostStepDoIt(), and G4HadronElasticProcess::PostStepDoIt().

◆ ComputeCrossSection()

G4double G4HadronicProcess::ComputeCrossSection ( const G4ParticleDefinition *  part,
const G4Material *  mat,
const G4double  kinEnergy 
)

Definition at line 913 of file G4HadronicProcess.cc.

916{
917 auto dp = new G4DynamicParticle(part, unitVector, kinEnergy);
918 G4double xs = theCrossSectionDataStore->ComputeCrossSection(dp, mat);
919 delete dp;
920 return xs;
921}
G4CrossSectionDataStore::ComputeCrossSection
G4double ComputeCrossSection(const G4DynamicParticle *, const G4Material *)
Definition: G4CrossSectionDataStore.cc:68

Referenced by G4HadXSHelper::FillPeaksStructure(), and G4HadXSHelper::FindCrossSectionMax().

◆ CrossSectionFactor()

G4double G4HadronicProcess::CrossSectionFactor ( ) const
inline

Definition at line 306 of file G4HadronicProcess.hh.

307{
308 return aScaleFactor;
309}

◆ CrossSectionType()

G4HadXSType G4HadronicProcess::CrossSectionType ( ) const
inline

Definition at line 295 of file G4HadronicProcess.hh.

296{
297 return fXSType;
298}

Referenced by BuildPhysicsTable().

◆ DumpPhysicsTable()

void G4HadronicProcess::DumpPhysicsTable ( const G4ParticleDefinition &  p)

Definition at line 885 of file G4HadronicProcess.cc.

886{
887 theCrossSectionDataStore->DumpPhysicsTable(p);
888}
G4CrossSectionDataStore::DumpPhysicsTable
void DumpPhysicsTable(const G4ParticleDefinition &)
Definition: G4CrossSectionDataStore.cc:304

◆ DumpState()

void G4HadronicProcess::DumpState ( const G4Track &  aTrack,
const G4String &  method,
G4ExceptionDescription &  ed 
)
protected

Definition at line 860 of file G4HadronicProcess.cc.

863{
864 ed << "Unrecoverable error in the method " << method << " of "
865 << GetProcessName() << G4endl;
866 ed << "TrackID= "<< aTrack.GetTrackID() << " ParentID= "
867 << aTrack.GetParentID()
868 << " " << aTrack.GetParticleDefinition()->GetParticleName()
869 << G4endl;
870 ed << "Ekin(GeV)= " << aTrack.GetKineticEnergy()/CLHEP::GeV
871 << "; direction= " << aTrack.GetMomentumDirection() << G4endl;
872 ed << "Position(mm)= " << aTrack.GetPosition()/CLHEP::mm << ";";
873
874 if (aTrack.GetMaterial()) {
875 ed << " material " << aTrack.GetMaterial()->GetName();
876 }
877 ed << G4endl;
878
879 if (aTrack.GetVolume()) {
880 ed << "PhysicalVolume <" << aTrack.GetVolume()->GetName()
881 << ">" << G4endl;
882 }
883}
G4Material::GetName
const G4String & GetName() const
Definition: G4Material.hh:172
G4Track::GetTrackID
G4int GetTrackID() const
G4Track::GetVolume
G4VPhysicalVolume * GetVolume() const
G4Track::GetPosition
const G4ThreeVector & GetPosition() const
G4Track::GetMaterial
G4Material * GetMaterial() const
G4Track::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4Track::GetParentID
G4int GetParentID() const
G4VPhysicalVolume::GetName
const G4String & GetName() const

Referenced by G4HadronStoppingProcess::AtRestDoIt(), FillResult(), G4NeutrinoElectronProcess::GetMeanFreePath(), G4NeutrinoElectronProcess::PostStepDoIt(), PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), G4MuNeutrinoNucleusProcess::PostStepDoIt(), and G4TauNeutrinoNucleusProcess::PostStepDoIt().

◆ EnergyOfCrossSectionMax()

std::vector< G4double > * G4HadronicProcess::EnergyOfCrossSectionMax ( ) const
inline

Definition at line 349 of file G4HadronicProcess.hh.

350{
351 return theEnergyOfCrossSectionMax;
352}

Referenced by BuildPhysicsTable().

◆ FillResult()

void G4HadronicProcess::FillResult ( G4HadFinalState *  aR,
const G4Track &  aT 
)
protected

Definition at line 507 of file G4HadronicProcess.cc.

508{
509 theTotalResult->ProposeLocalEnergyDeposit(aR->GetLocalEnergyDeposit());
510 const G4ThreeVector& dir = aT.GetMomentumDirection();
511
512 G4double efinal = std::max(aR->GetEnergyChange(), 0.0);
513
514 // check status of primary
515 if(aR->GetStatusChange() == stopAndKill) {
516 theTotalResult->ProposeTrackStatus(fStopAndKill);
517 theTotalResult->ProposeEnergy( 0.0 );
518
519 // check its final energy
520 } else if(0.0 == efinal) {
521 theTotalResult->ProposeEnergy( 0.0 );
522 if(aT.GetParticleDefinition()->GetProcessManager()
523 ->GetAtRestProcessVector()->size() > 0)
524 { theTotalResult->ProposeTrackStatus(fStopButAlive); }
525 else { theTotalResult->ProposeTrackStatus(fStopAndKill); }
526
527 // primary is not killed apply rotation and Lorentz transformation
528 } else {
529 theTotalResult->ProposeTrackStatus(fAlive);
530 G4ThreeVector newDir = aR->GetMomentumChange();
531 newDir.rotateUz(dir);
532 theTotalResult->ProposeMomentumDirection(newDir);
533 theTotalResult->ProposeEnergy(efinal);
534 }
535 //G4cout << "FillResult: Efinal= " << efinal << " status= "
536 // << theTotalResult->GetTrackStatus()
537 // << " fKill= " << fStopAndKill << G4endl;
538
539 // check secondaries
540 nICelectrons = 0;
541 G4int nSec = (G4int)aR->GetNumberOfSecondaries();
542 theTotalResult->SetNumberOfSecondaries(nSec);
543 G4double time0 = aT.GetGlobalTime();
544
545 for (G4int i = 0; i < nSec; ++i) {
546 G4DynamicParticle* dynParticle = aR->GetSecondary(i)->GetParticle();
547
548 // apply rotation
549 G4ThreeVector newDir = dynParticle->GetMomentumDirection();
550 newDir.rotateUz(dir);
551 dynParticle->SetMomentumDirection(newDir);
552
553 // check if secondary is on the mass shell
554 const G4ParticleDefinition* part = dynParticle->GetDefinition();
555 G4double mass = part->GetPDGMass();
556 G4double dmass= dynParticle->GetMass();
557 const G4double delta_mass_lim = 1.0*CLHEP::keV;
558 const G4double delta_ekin = 0.001*CLHEP::eV;
559 if(std::abs(dmass - mass) > delta_mass_lim) {
560 G4double e =
561 std::max(dynParticle->GetKineticEnergy() + dmass - mass, delta_ekin);
562 if(G4HadronicProcess_debug_flag) {
563 G4ExceptionDescription ed;
564 ed << "TrackID= "<< aT.GetTrackID()
565 << " " << aT.GetParticleDefinition()->GetParticleName()
566 << " Target Z= " << targetNucleus.GetZ_asInt() << " A= "
567 << targetNucleus.GetA_asInt()
568 << " Ekin(GeV)= " << aT.GetKineticEnergy()/CLHEP::GeV
569 << "\n Secondary is out of mass shell: " << part->GetParticleName()
570 << " EkinNew(MeV)= " << e
571 << " DeltaMass(MeV)= " << dmass - mass << G4endl;
572 G4Exception("G4HadronicProcess::FillResults", "had012", JustWarning, ed);
573 }
574 dynParticle->SetKineticEnergy(e);
575 dynParticle->SetMass(mass);
576 }
577 G4int idModel = aR->GetSecondary(i)->GetCreatorModelID();
578 if(part->GetPDGEncoding() == 11) { ++nICelectrons; }
579
580 // time of interaction starts from zero + global time
581 G4double time = std::max(aR->GetSecondary(i)->GetTime(), 0.0) + time0;
582
583 G4Track* track = new G4Track(dynParticle, time, aT.GetPosition());
584 track->SetCreatorModelID(idModel);
585 track->SetParentResonanceDef(aR->GetSecondary(i)->GetParentResonanceDef());
586 track->SetParentResonanceID(aR->GetSecondary(i)->GetParentResonanceID());
587 G4double newWeight = fWeight*aR->GetSecondary(i)->GetWeight();
588 track->SetWeight(newWeight);
589 track->SetTouchableHandle(aT.GetTouchableHandle());
590 theTotalResult->AddSecondary(track);
591 if (G4HadronicProcess_debug_flag) {
592 G4double e = dynParticle->GetKineticEnergy();
593 if (e == 0.0) {
594 G4ExceptionDescription ed;
595 DumpState(aT,"Secondary has zero energy",ed);
596 ed << "Secondary " << part->GetParticleName()
597 << G4endl;
598 G4Exception("G4HadronicProcess::FillResults", "had011",
599 JustWarning,ed);
600 }
601 }
602 }
603 aR->Clear();
604 // G4cout << "FillResults done nICe= " << nICelectrons << G4endl;
605}
fAlive
@ fAlive
Definition: G4TrackStatus.hh:42
fStopButAlive
@ fStopButAlive
Definition: G4TrackStatus.hh:44
CLHEP::Hep3Vector::rotateUz
Hep3Vector & rotateUz(const Hep3Vector &)
Definition: ThreeVector.cc:33
G4DynamicParticle::SetMomentumDirection
void SetMomentumDirection(const G4ThreeVector &aDirection)
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::SetMass
void SetMass(G4double mass)
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4DynamicParticle::SetKineticEnergy
void SetKineticEnergy(G4double aEnergy)
G4HadFinalState::GetMomentumChange
const G4ThreeVector & GetMomentumChange() const
Definition: G4HadFinalState.hh:59
G4HadSecondary::GetParentResonanceID
G4int GetParentResonanceID() const
Definition: G4HadSecondary.hh:53
G4HadSecondary::GetWeight
G4double GetWeight() const
Definition: G4HadSecondary.hh:45
G4HadSecondary::GetParentResonanceDef
const G4ParticleDefinition * GetParentResonanceDef() const
Definition: G4HadSecondary.hh:51
G4HadSecondary::GetTime
G4double GetTime() const
Definition: G4HadSecondary.hh:48
G4HadSecondary::GetCreatorModelID
G4int GetCreatorModelID() const
Definition: G4HadSecondary.hh:50
G4HadronicProcess::DumpState
void DumpState(const G4Track &, const G4String &, G4ExceptionDescription &)
Definition: G4HadronicProcess.cc:860
G4ParticleChange::AddSecondary
void AddSecondary(G4Track *aSecondary)
Definition: G4ParticleChange.cc:100
G4ParticleChange::ProposeEnergy
void ProposeEnergy(G4double finalEnergy)
G4ParticleChange::ProposeMomentumDirection
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
Definition: G4ParticleDefinition.cc:238
G4ProcessManager::GetAtRestProcessVector
G4ProcessVector * GetAtRestProcessVector(G4ProcessVectorTypeIndex typ=typeGPIL) const
G4ProcessVector::size
std::size_t size() const
G4Track::SetWeight
void SetWeight(G4double aValue)
G4Track::SetParentResonanceID
void SetParentResonanceID(const G4int parentID)
G4Track::SetTouchableHandle
void SetTouchableHandle(const G4TouchableHandle &apValue)
G4Track::GetGlobalTime
G4double GetGlobalTime() const
G4Track::GetTouchableHandle
const G4TouchableHandle & GetTouchableHandle() const
G4Track::SetCreatorModelID
void SetCreatorModelID(const G4int id)
G4Track::SetParentResonanceDef
void SetParentResonanceDef(const G4ParticleDefinition *parent)
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4VParticleChange::ProposeLocalEnergyDeposit
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4VParticleChange::SetNumberOfSecondaries
void SetNumberOfSecondaries(G4int totSecondaries)

Referenced by G4NeutrinoElectronProcess::PostStepDoIt(), PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), G4MuNeutrinoNucleusProcess::PostStepDoIt(), and G4TauNeutrinoNucleusProcess::PostStepDoIt().

◆ GetCrossSectionDataStore()

G4CrossSectionDataStore * G4HadronicProcess::GetCrossSectionDataStore ( )
inline

Definition at line 337 of file G4HadronicProcess.hh.

338{
339 return theCrossSectionDataStore;
340}

Referenced by G4GammaGeneralProcess::BuildPhysicsTable(), G4ChargeExchangeProcess::BuildPhysicsTable(), G4ElectronNuclearProcess::G4ElectronNuclearProcess(), G4PositronNuclearProcess::G4PositronNuclearProcess(), G4NeutrinoElectronProcess::GetMeanFreePath(), G4ElNeutrinoNucleusProcess::GetMeanFreePath(), G4MuNeutrinoNucleusProcess::GetMeanFreePath(), G4TauNeutrinoNucleusProcess::GetMeanFreePath(), G4HadronPhysicsShielding::Neutron(), G4NeutrinoElectronProcess::PostStepDoIt(), G4HadronicProcessStore::PrintHtml(), G4GammaGeneralProcess::SelectHadProcess(), G4NeutronGeneralProcess::SetCaptureProcess(), G4NeutronGeneralProcess::SetElasticProcess(), and G4NeutronGeneralProcess::SetInelasticProcess().

◆ GetElementCrossSection()

G4double G4HadronicProcess::GetElementCrossSection ( const G4DynamicParticle *  part,
const G4Element *  elm,
const G4Material *  mat = nullptr 
)

Definition at line 148 of file G4HadronicProcess.cc.

151{
152 if(nullptr == mat)
153 {
154 static const G4int nmax = 5;
155 if(nMatWarn < nmax) {
156 ++nMatWarn;
157 G4ExceptionDescription ed;
158 ed << "Cannot compute Element x-section for " << GetProcessName()
159 << " because no material defined \n"
160 << " Please, specify material pointer or define simple material"
161 << " for Z= " << elm->GetZasInt();
162 G4Exception("G4HadronicProcess::GetElementCrossSection", "had066",
163 JustWarning, ed);
164 }
165 }
166 return theCrossSectionDataStore->GetCrossSection(dp, elm, mat);
167}
G4CrossSectionDataStore::GetCrossSection
G4double GetCrossSection(const G4DynamicParticle *, const G4Material *)
Definition: G4CrossSectionDataStore.hh:149
G4Element::GetZasInt
G4int GetZasInt() const
Definition: G4Element.hh:132

Referenced by G4HadronicProcessStore::GetCaptureCrossSectionPerAtom(), G4HadronicProcessStore::GetChargeExchangeCrossSectionPerAtom(), G4HadronicProcessStore::GetElasticCrossSectionPerAtom(), G4HadronicProcessStore::GetFissionCrossSectionPerAtom(), G4HadronicProcessStore::GetInelasticCrossSectionPerAtom(), and GetMicroscopicCrossSection().

◆ GetEnergyMomentumCheckLevels()

std::pair< G4double, G4double > G4HadronicProcess::GetEnergyMomentumCheckLevels ( ) const
inline

Definition at line 331 of file G4HadronicProcess.hh.

332{
333 return epCheckLevels;
334}

Referenced by G4HadronicProcessStore::SetProcessAbsLevel(), and G4HadronicProcessStore::SetProcessRelLevel().

◆ GetHadronicInteraction()

G4HadronicInteraction * G4HadronicProcess::GetHadronicInteraction ( ) const
inline

Definition at line 278 of file G4HadronicProcess.hh.

279{
280 return theInteraction;
281}

Referenced by CheckEnergyMomentumConservation(), and CheckResult().

◆ GetHadronicInteractionList()

std::vector< G4HadronicInteraction * > & G4HadronicProcess::GetHadronicInteractionList ( )

Definition at line 896 of file G4HadronicProcess.cc.

897{
898 return theEnergyRangeManager.GetHadronicInteractionList();
899}
G4EnergyRangeManager::GetHadronicInteractionList
std::vector< G4HadronicInteraction * > & GetHadronicInteractionList()
Definition: G4EnergyRangeManager.cc:156

Referenced by G4ThermalNeutrons::ConstructProcess(), G4HadronElasticPhysics::GetElasticModel(), G4NeutrinoElectronProcess::PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4MuNeutrinoNucleusProcess::PostStepDoIt(), and G4TauNeutrinoNucleusProcess::PostStepDoIt().

◆ GetHadronicModel()

G4HadronicInteraction * G4HadronicProcess::GetHadronicModel ( const G4String &  modelName)

Definition at line 902 of file G4HadronicProcess.cc.

903{
904 std::vector<G4HadronicInteraction*>& list
905 = theEnergyRangeManager.GetHadronicInteractionList();
906 for (auto & mod : list) {
907 if (mod->GetModelName() == modelName) return mod;
908 }
909 return nullptr;
910}

◆ GetLastCrossSection()

G4double G4HadronicProcess::GetLastCrossSection ( )
inlineprotected

Definition at line 370 of file G4HadronicProcess.hh.

371{
372 return theLastCrossSection;
373}

◆ GetMeanFreePath()

G4double G4HadronicProcess::GetMeanFreePath ( const G4Track &  aTrack,
G4double  ,
G4ForceCondition *   
)
overridevirtual

Implements G4VDiscreteProcess.

Reimplemented in G4NeutrinoElectronProcess, G4ElNeutrinoNucleusProcess, G4MuNeutrinoNucleusProcess, G4TauNeutrinoNucleusProcess, and G4NeutronGeneralProcess.

Definition at line 314 of file G4HadronicProcess.cc.

317{
318 G4double xs = aScaleFactor*theCrossSectionDataStore
319 ->ComputeCrossSection(aTrack.GetDynamicParticle(),aTrack.GetMaterial());
320 return (xs > 0.0) ? 1.0/xs : DBL_MAX;
321}

◆ GetMicroscopicCrossSection()

G4double G4HadronicProcess::GetMicroscopicCrossSection ( const G4DynamicParticle *  part,
const G4Element *  elm,
const G4Material *  mat = nullptr 
)
inline

Definition at line 269 of file G4HadronicProcess.hh.

273{
274 return GetElementCrossSection(part, elm, mat);
275}
G4HadronicProcess::GetElementCrossSection
G4double GetElementCrossSection(const G4DynamicParticle *part, const G4Element *elm, const G4Material *mat=nullptr)
Definition: G4HadronicProcess.cc:148

◆ GetTargetIsotope()

const G4Isotope * G4HadronicProcess::GetTargetIsotope ( )
inline

Definition at line 289 of file G4HadronicProcess.hh.

290{
291 return targetNucleus.GetIsotope();
292}
G4Nucleus::GetIsotope
const G4Isotope * GetIsotope()
Definition: G4Nucleus.hh:111

◆ GetTargetNucleus()

const G4Nucleus * G4HadronicProcess::GetTargetNucleus ( ) const
inline

Definition at line 284 of file G4HadronicProcess.hh.

285{
286 return &targetNucleus;
287}

◆ GetTargetNucleusPointer()

G4Nucleus * G4HadronicProcess::GetTargetNucleusPointer ( )
inline

Definition at line 365 of file G4HadronicProcess.hh.

366{
367 return &targetNucleus;
368}

Referenced by G4HadronStoppingProcess::AtRestDoIt(), G4NeutrinoElectronProcess::PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), G4MuNeutrinoNucleusProcess::PostStepDoIt(), and G4TauNeutrinoNucleusProcess::PostStepDoIt().

◆ MultiplyCrossSectionBy()

void G4HadronicProcess::MultiplyCrossSectionBy ( G4double  factor)

Definition at line 607 of file G4HadronicProcess.cc.

608{
609 BiasCrossSectionByFactor(factor);
610}
G4HadronicProcess::BiasCrossSectionByFactor
void BiasCrossSectionByFactor(G4double aScale)
Definition: G4HadronicProcess.cc:612

Referenced by G4HadProcesses::BuildNeutronElastic(), G4HadProcesses::BuildNeutronInelasticAndCapture(), G4HadronDElasticPhysics::ConstructProcess(), G4HadronElasticPhysics::ConstructProcess(), G4HadronHElasticPhysics::ConstructProcess(), G4HadronInelasticQBBC::ConstructProcess(), G4HadronPhysicsFTFQGSP_BERT::ConstructProcess(), G4HadronPhysicsFTFP_BERT::Kaon(), G4HadronPhysicsFTF_BIC::Kaon(), G4HadronPhysicsINCLXX::Kaon(), G4HadronPhysicsFTFP_BERT::Neutron(), G4HadronPhysicsQGSP_BERT::Neutron(), G4HadronPhysicsQGSP_BIC::Neutron(), G4HadronPhysicsFTF_BIC::Neutron(), G4HadronPhysicsFTFP_BERT_HP::Neutron(), G4HadronPhysicsINCLXX::Neutron(), G4HadronPhysicsQGS_BIC::Neutron(), G4HadronPhysicsQGSP_BERT_HP::Neutron(), G4HadronPhysicsQGSP_BIC_HP::Neutron(), G4HadronPhysicsShielding::Neutron(), G4HadronPhysicsFTFP_BERT::Pion(), G4HadronPhysicsQGSP_BERT::Pion(), G4HadronPhysicsQGSP_BIC::Pion(), G4HadronPhysicsFTF_BIC::Pion(), G4HadronPhysicsINCLXX::Pion(), G4HadronPhysicsQGS_BIC::Pion(), G4HadronPhysicsFTFP_BERT::Proton(), G4HadronPhysicsQGSP_BERT::Proton(), G4HadronPhysicsQGSP_BIC::Proton(), G4HadronPhysicsFTF_BIC::Proton(), G4HadronPhysicsINCLXX::Proton(), G4HadronPhysicsNuBeam::Proton(), G4HadronPhysicsQGS_BIC::Proton(), and G4HadronPhysicsQGSP_BIC_AllHP::Proton().

◆ operator=()

G4HadronicProcess & G4HadronicProcess::operator= ( const G4HadronicProcess &  right)
delete

◆ PostStepDoIt()

G4VParticleChange * G4HadronicProcess::PostStepDoIt ( const G4Track &  aTrack,
const G4Step &  aStep 
)
overridevirtual

Reimplemented from G4VDiscreteProcess.

Reimplemented in G4NeutronGeneralProcess, G4NeutrinoElectronProcess, G4ElNeutrinoNucleusProcess, G4HadronElasticProcess, G4MuNeutrinoNucleusProcess, and G4TauNeutrinoNucleusProcess.

Definition at line 324 of file G4HadronicProcess.cc.

325{
326 theNumberOfInteractionLengthLeft = -1.0;
327
328 //G4cout << "PostStepDoIt " << aTrack.GetDefinition()->GetParticleName()
329 // << " Ekin= " << aTrack.GetKineticEnergy() << G4endl;
330 // if primary is not Alive then do nothing
331 theTotalResult->Clear();
332 theTotalResult->Initialize(aTrack);
333 fWeight = aTrack.GetWeight();
334 theTotalResult->ProposeWeight(fWeight);
335 if(aTrack.GetTrackStatus() != fAlive) { return theTotalResult; }
336
337 // Find cross section at end of step and check if <= 0
338 //
339 const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
340 const G4Material* aMaterial = aTrack.GetMaterial();
341
342 // check only for charged particles
343 if(fXSType != fHadNoIntegral) {
344 mfpKinEnergy = DBL_MAX;
345 G4double xs = aScaleFactor*
346 theCrossSectionDataStore->ComputeCrossSection(aParticle,aMaterial);
347 //G4cout << "xs=" << xs << " xs0=" << theLastCrossSection
348 // << " " << aMaterial->GetName() << G4endl;
349 if(xs < theLastCrossSection*G4UniformRand()) {
350 // No interaction
351 return theTotalResult;
352 }
353 }
354
355 const G4Element* anElement =
356 theCrossSectionDataStore->SampleZandA(aParticle,aMaterial,targetNucleus);
357
358 // Next check for illegal track status
359 //
360 if (aTrack.GetTrackStatus() != fAlive &&
361 aTrack.GetTrackStatus() != fSuspend) {
362 if (aTrack.GetTrackStatus() == fStopAndKill ||
363 aTrack.GetTrackStatus() == fKillTrackAndSecondaries ||
364 aTrack.GetTrackStatus() == fPostponeToNextEvent) {
365 G4ExceptionDescription ed;
366 ed << "G4HadronicProcess: track in unusable state - "
367 << aTrack.GetTrackStatus() << G4endl;
368 ed << "G4HadronicProcess: returning unchanged track " << G4endl;
369 DumpState(aTrack,"PostStepDoIt",ed);
370 G4Exception("G4HadronicProcess::PostStepDoIt", "had004", JustWarning, ed);
371 }
372 // No warning for fStopButAlive which is a legal status here
373 return theTotalResult;
374 }
375
376 // Initialize the hadronic projectile from the track
377 thePro.Initialise(aTrack);
378
379 theInteraction = ChooseHadronicInteraction(thePro, targetNucleus,
380 aMaterial, anElement);
381 if(nullptr == theInteraction) {
382 G4ExceptionDescription ed;
383 ed << "Target element "<<anElement->GetName()<<" Z= "
384 << targetNucleus.GetZ_asInt() << " A= "
385 << targetNucleus.GetA_asInt() << G4endl;
386 DumpState(aTrack,"ChooseHadronicInteraction",ed);
387 ed << " No HadronicInteraction found out" << G4endl;
388 G4Exception("G4HadronicProcess::PostStepDoIt", "had005",
389 FatalException, ed);
390 return theTotalResult;
391 }
392
393 G4HadFinalState* result = nullptr;
394 G4int reentryCount = 0;
395 /*
396 G4cout << "### " << aParticle->GetDefinition()->GetParticleName()
397 << " Ekin(MeV)= " << aParticle->GetKineticEnergy()
398 << " Z= " << targetNucleus.GetZ_asInt()
399 << " A= " << targetNucleus.GetA_asInt()
400 << " by " << theInteraction->GetModelName()
401 << G4endl;
402 */
403 do
404 {
405 try
406 {
407 // Save random engine if requested for debugging
408 if (G4Hadronic_Random_File) {
409 CLHEP::HepRandom::saveEngineStatus(G4Hadronic_Random_File);
410 }
411 // Call the interaction
412 result = theInteraction->ApplyYourself( thePro, targetNucleus);
413 ++reentryCount;
414 }
415 catch(G4HadronicException & aR)
416 {
417 G4ExceptionDescription ed;
418 aR.Report(ed);
419 ed << "Call for " << theInteraction->GetModelName() << G4endl;
420 ed << "Target element "<<anElement->GetName()<<" Z= "
421 << targetNucleus.GetZ_asInt()
422 << " A= " << targetNucleus.GetA_asInt() << G4endl;
423 DumpState(aTrack,"ApplyYourself",ed);
424 ed << " ApplyYourself failed" << G4endl;
425 G4Exception("G4HadronicProcess::PostStepDoIt", "had006", FatalException,
426 ed);
427 }
428
429 // Check the result for catastrophic energy non-conservation
430 result = CheckResult(thePro, targetNucleus, result);
431
432 if(reentryCount>100) {
433 G4ExceptionDescription ed;
434 ed << "Call for " << theInteraction->GetModelName() << G4endl;
435 ed << "Target element "<<anElement->GetName()<<" Z= "
436 << targetNucleus.GetZ_asInt()
437 << " A= " << targetNucleus.GetA_asInt() << G4endl;
438 DumpState(aTrack,"ApplyYourself",ed);
439 ed << " ApplyYourself does not completed after 100 attempts" << G4endl;
440 G4Exception("G4HadronicProcess::PostStepDoIt", "had006", FatalException,
441 ed);
442 }
443 }
444 while(!result); /* Loop checking, 30-Oct-2015, G.Folger */
445
446 // Check whether kaon0 or anti_kaon0 are present between the secondaries:
447 // if this is the case, transform them into either kaon0S or kaon0L,
448 // with equal, 50% probability, keeping their dynamical masses (and
449 // the other kinematical properties).
450 // When this happens - very rarely - a "JustWarning" exception is thrown.
451 G4int nSec = (G4int)result->GetNumberOfSecondaries();
452 if ( nSec > 0 ) {
453 for ( G4int i = 0; i < nSec; ++i ) {
454 auto dynamicParticle = result->GetSecondary(i)->GetParticle();
455 auto part = dynamicParticle->GetParticleDefinition();
456 if ( part == G4KaonZero::Definition() ||
457 part == G4AntiKaonZero::Definition() ) {
458 G4ParticleDefinition* newPart;
459 if( G4UniformRand() > 0.5 ) { newPart = G4KaonZeroShort::Definition(); }
460 else { newPart = G4KaonZeroLong::Definition(); }
461 dynamicParticle->SetDefinition( newPart );
462 if(nKaonWarn < 5) {
463 ++nKaonWarn;
464 G4ExceptionDescription ed;
465 ed << " Hadronic model " << theInteraction->GetModelName() << G4endl;
466 ed << " created " << part->GetParticleName() << G4endl;
467 ed << " -> forced to be " << newPart->GetParticleName() << G4endl;
468 G4Exception( "G4HadronicProcess::PostStepDoIt", "had007", JustWarning, ed );
469 }
470 }
471 }
472 }
473
474 result->SetTrafoToLab(thePro.GetTrafoToLab());
475 FillResult(result, aTrack);
476
477 if (epReportLevel != 0) {
478 CheckEnergyMomentumConservation(aTrack, targetNucleus);
479 }
480 //G4cout << "PostStepDoIt done nICelectrons= " << nICelectrons << G4endl;
481 return theTotalResult;
482}
FatalException
@ FatalException
Definition: G4ExceptionSeverity.hh:61
fKillTrackAndSecondaries
@ fKillTrackAndSecondaries
Definition: G4TrackStatus.hh:48
fSuspend
@ fSuspend
Definition: G4TrackStatus.hh:50
fPostponeToNextEvent
@ fPostponeToNextEvent
Definition: G4TrackStatus.hh:52
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
CLHEP::HepRandom::saveEngineStatus
static void saveEngineStatus(const char filename[]="Config.conf")
Definition: Random.cc:278
G4AntiKaonZero::Definition
static G4AntiKaonZero * Definition()
Definition: G4AntiKaonZero.cc:52
G4CrossSectionDataStore::SampleZandA
const G4Element * SampleZandA(const G4DynamicParticle *, const G4Material *, G4Nucleus &target)
Definition: G4CrossSectionDataStore.cc:192
G4DynamicParticle::GetParticleDefinition
const G4ParticleDefinition * GetParticleDefinition() const
G4Element::GetName
const G4String & GetName() const
Definition: G4Element.hh:127
G4HadFinalState::SetTrafoToLab
void SetTrafoToLab(const G4LorentzRotation &aT)
Definition: G4HadFinalState.hh:74
G4HadProjectile::Initialise
void Initialise(const G4Track &aT)
Definition: G4HadProjectile.cc:54
G4HadProjectile::GetTrafoToLab
G4LorentzRotation & GetTrafoToLab()
Definition: G4HadProjectile.hh:96
G4HadronicException::Report
void Report(std::ostream &aS) const
Definition: G4HadronicException.cc:52
G4HadronicInteraction::ApplyYourself
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
Definition: G4HadronicInteraction.cc:63
G4HadronicProcess::FillResult
void FillResult(G4HadFinalState *aR, const G4Track &aT)
Definition: G4HadronicProcess.cc:507
G4HadronicProcess::thePro
G4HadProjectile thePro
Definition: G4HadronicProcess.hh:218
G4HadronicProcess::CheckResult
G4HadFinalState * CheckResult(const G4HadProjectile &thePro, const G4Nucleus &targetNucleus, G4HadFinalState *result)
Definition: G4HadronicProcess.cc:624
G4HadronicProcess::ChooseHadronicInteraction
G4HadronicInteraction * ChooseHadronicInteraction(const G4HadProjectile &aHadProjectile, G4Nucleus &aTargetNucleus, const G4Material *aMaterial, const G4Element *anElement)
Definition: G4HadronicProcess.hh:355
G4HadronicProcess::CheckEnergyMomentumConservation
void CheckEnergyMomentumConservation(const G4Track &, const G4Nucleus &)
Definition: G4HadronicProcess.cc:715
G4KaonZeroLong::Definition
static G4KaonZeroLong * Definition()
Definition: G4KaonZeroLong.cc:55
G4KaonZeroShort::Definition
static G4KaonZeroShort * Definition()
Definition: G4KaonZeroShort.cc:52
G4KaonZero::Definition
static G4KaonZero * Definition()
Definition: G4KaonZero.cc:52
G4ParticleChange::Initialize
void Initialize(const G4Track &) override
Definition: G4ParticleChange.cc:107
G4Track::GetTrackStatus
G4TrackStatus GetTrackStatus() const
G4Track::GetWeight
G4double GetWeight() const
G4VParticleChange::ProposeWeight
void ProposeWeight(G4double finalWeight)
G4VProcess::theNumberOfInteractionLengthLeft
G4double theNumberOfInteractionLengthLeft
Definition: G4VProcess.hh:335

Referenced by G4NeutronGeneralProcess::PostStepDoIt().

◆ PostStepGetPhysicalInteractionLength()

G4double G4HadronicProcess::PostStepGetPhysicalInteractionLength ( const G4Track &  track,
G4double  previousStepSize,
G4ForceCondition *  condition 
)
overridevirtual

Reimplemented from G4VDiscreteProcess.

Reimplemented in G4HadronStoppingProcess, and G4NeutronGeneralProcess.

Definition at line 268 of file G4HadronicProcess.cc.

272{
273 *condition = NotForced;
274
275 const G4Material* mat = track.GetMaterial();
276 if(mat != currentMat) {
277 currentMat = mat;
278 mfpKinEnergy = DBL_MAX;
279 matIdx = (G4int)track.GetMaterial()->GetIndex();
280 }
281 /*
282 G4cout << GetProcessName() << " E=" << track.GetKineticEnergy()
283 << " " << currentParticle->GetParticleName()
284 << " lastxs=" << theLastCrossSection
285 << " lastmfp=" << theMFP << G4endl;
286 */
287 UpdateCrossSectionAndMFP(track.GetKineticEnergy());
288 /*
289 G4cout << " xs=" << theLastCrossSection
290 << " mfp=" << theMFP << " nleft=" << theNumberOfInteractionLengthLeft
291 << G4endl;
292 */
293 // zero cross section
294 if(theLastCrossSection <= 0.0) {
295 theNumberOfInteractionLengthLeft = -1.0;
296 currentInteractionLength = DBL_MAX;
297 return DBL_MAX;
298 }
299
300 // non-zero cross section
301 if (theNumberOfInteractionLengthLeft < 0.0) {
302 theNumberOfInteractionLengthLeft = -G4Log( G4UniformRand() );
303 theInitialNumberOfInteractionLength = theNumberOfInteractionLengthLeft;
304 } else {
305 theNumberOfInteractionLengthLeft -=
306 previousStepSize/currentInteractionLength;
307 theNumberOfInteractionLengthLeft =
308 std::max(theNumberOfInteractionLengthLeft, 0.0);
309 }
310 currentInteractionLength = theMFP;
311 return theNumberOfInteractionLengthLeft*theMFP;
312}
condition
G4double condition(const G4ErrorSymMatrix &m)
NotForced
@ NotForced
Definition: G4ForceCondition.hh:47
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:227
G4Material::GetIndex
size_t GetIndex() const
Definition: G4Material.hh:255
G4VProcess::currentInteractionLength
G4double currentInteractionLength
Definition: G4VProcess.hh:339
G4VProcess::theInitialNumberOfInteractionLength
G4double theInitialNumberOfInteractionLength
Definition: G4VProcess.hh:342

◆ PreparePhysicsTable()

void G4HadronicProcess::PreparePhysicsTable ( const G4ParticleDefinition &  p)
overridevirtual

Reimplemented from G4VProcess.

Reimplemented in G4NeutrinoElectronProcess, G4HadronStoppingProcess, G4ElNeutrinoNucleusProcess, G4MuNeutrinoNucleusProcess, G4NeutronGeneralProcess, and G4TauNeutrinoNucleusProcess.

Definition at line 169 of file G4HadronicProcess.cc.

170{
171 if(std::getenv("G4HadronicProcess_debug")) {
172 G4HadronicProcess_debug_flag = true;
173 }
174 if(nullptr == firstParticle) { firstParticle = &p; }
175 theProcessStore->RegisterParticle(this, &p);
176}
G4HadronicProcessStore::RegisterParticle
void RegisterParticle(G4HadronicProcess *, const G4ParticleDefinition *)
Definition: G4HadronicProcessStore.cc:413

Referenced by G4NeutrinoElectronProcess::PreparePhysicsTable(), G4GammaGeneralProcess::PreparePhysicsTable(), G4ElNeutrinoNucleusProcess::PreparePhysicsTable(), G4MuNeutrinoNucleusProcess::PreparePhysicsTable(), G4NeutronGeneralProcess::PreparePhysicsTable(), and G4TauNeutrinoNucleusProcess::PreparePhysicsTable().

◆ ProcessDescription()

void G4HadronicProcess::ProcessDescription ( std::ostream &  outFile) const
overridevirtual

Reimplemented from G4VProcess.

Reimplemented in G4ElectronNuclearProcess, G4MuonNuclearProcess, G4NeutrinoElectronProcess, G4NeutronFissionProcess, G4PositronNuclearProcess, G4HadronicAbsorptionBertini, G4HadronicAbsorptionFritiof, G4HadronicAbsorptionFritiofWithBinaryCascade, G4HadronStoppingProcess, G4MuonMinusCapture, G4NeutronCaptureProcess, G4ElNeutrinoNucleusProcess, G4HadronElasticProcess, G4MuNeutrinoNucleusProcess, G4NeutronGeneralProcess, and G4TauNeutrinoNucleusProcess.

Definition at line 484 of file G4HadronicProcess.cc.

485{
486 outFile << "The description for this process has not been written yet.\n";
487}

Referenced by G4HadronicProcessStore::PrintHtml(), G4GammaGeneralProcess::ProcessDescription(), and G4NeutronGeneralProcess::ProcessDescription().

◆ RegisterMe()

void G4HadronicProcess::RegisterMe ( G4HadronicInteraction *  a)

Definition at line 140 of file G4HadronicProcess.cc.

141{
142 if(nullptr == a) { return; }
143 theEnergyRangeManager.RegisterMe( a );
144 G4HadronicProcessStore::Instance()->RegisterInteraction(this, a);
145}
G4EnergyRangeManager::RegisterMe
void RegisterMe(G4HadronicInteraction *a)
Definition: G4EnergyRangeManager.cc:46
G4HadronicProcessStore::RegisterInteraction
void RegisterInteraction(G4HadronicProcess *, G4HadronicInteraction *)
Definition: G4HadronicProcessStore.cc:448

Referenced by G4BertiniElectroNuclearBuilder::Build(), G4LENDBertiniGammaElectroNuclearBuilder::Build(), G4NeutronLENDBuilder::Build(), G4NeutronPHPBuilder::Build(), G4AlphaPHPBuilder::Build(), G4BertiniKaonBuilder::Build(), G4BertiniNeutronBuilder::Build(), G4BertiniPiKBuilder::Build(), G4BertiniPionBuilder::Build(), G4BertiniProtonBuilder::Build(), G4BinaryAlphaBuilder::Build(), G4BinaryDeuteronBuilder::Build(), G4BinaryHe3Builder::Build(), G4BinaryNeutronBuilder::Build(), G4BinaryPiKBuilder::Build(), G4BinaryPionBuilder::Build(), G4BinaryProtonBuilder::Build(), G4BinaryTritonBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4FTFBinaryKaonBuilder::Build(), G4FTFBinaryNeutronBuilder::Build(), G4FTFBinaryPiKBuilder::Build(), G4FTFBinaryPionBuilder::Build(), G4FTFBinaryProtonBuilder::Build(), G4FTFPAntiBarionBuilder::Build(), G4FTFPKaonBuilder::Build(), G4FTFPNeutronBuilder::Build(), G4FTFPPiKBuilder::Build(), G4FTFPPionBuilder::Build(), G4FTFPProtonBuilder::Build(), G4He3PHPBuilder::Build(), G4HyperonFTFPBuilder::Build(), G4HyperonQGSPBuilder::Build(), G4INCLXXNeutronBuilder::Build(), G4INCLXXPionBuilder::Build(), G4INCLXXProtonBuilder::Build(), G4PrecoNeutronBuilder::Build(), G4PrecoProtonBuilder::Build(), G4ProtonPHPBuilder::Build(), G4QGSBinaryKaonBuilder::Build(), G4QGSBinaryNeutronBuilder::Build(), G4QGSBinaryPiKBuilder::Build(), G4QGSBinaryPionBuilder::Build(), G4QGSBinaryProtonBuilder::Build(), G4QGSPAntiBarionBuilder::Build(), G4QGSPKaonBuilder::Build(), G4QGSPLundStrFragmProtonBuilder::Build(), G4QGSPNeutronBuilder::Build(), G4QGSPPiKBuilder::Build(), G4QGSPPionBuilder::Build(), G4TritonPHPBuilder::Build(), G4QGSPProtonBuilder::Build(), G4HadProcesses::BuildNeutronInelasticAndCapture(), LBE::ConstructHad(), G4HadronElasticPhysicsHP::ConstructProcess(), G4HadronElasticPhysicsLEND::ConstructProcess(), G4HadronElasticPhysicsPHP::ConstructProcess(), G4EmExtraPhysics::ConstructProcess(), G4ChargeExchangePhysics::ConstructProcess(), G4HadronDElasticPhysics::ConstructProcess(), G4HadronElasticPhysics::ConstructProcess(), G4HadronHElasticPhysics::ConstructProcess(), G4IonElasticPhysics::ConstructProcess(), G4ThermalNeutrons::ConstructProcess(), G4HadronInelasticQBBC::ConstructProcess(), G4HadronPhysicsFTFQGSP_BERT::ConstructProcess(), G4HadronicAbsorptionBertini::G4HadronicAbsorptionBertini(), G4HadronicAbsorptionFritiof::G4HadronicAbsorptionFritiof(), G4HadronicAbsorptionFritiofWithBinaryCascade::G4HadronicAbsorptionFritiofWithBinaryCascade(), G4MuonMinusCapture::G4MuonMinusCapture(), G4HadronPhysicsFTFP_BERT::Neutron(), G4HadronPhysicsQGSP_BERT::Neutron(), G4HadronPhysicsQGSP_BIC::Neutron(), G4HadronPhysicsFTF_BIC::Neutron(), G4HadronPhysicsFTFP_BERT_HP::Neutron(), G4HadronPhysicsINCLXX::Neutron(), G4HadronPhysicsQGS_BIC::Neutron(), G4HadronPhysicsQGSP_BERT_HP::Neutron(), G4HadronPhysicsQGSP_BIC_HP::Neutron(), and G4HadronPhysicsShielding::Neutron().

◆ SetCrossSectionType()

void G4HadronicProcess::SetCrossSectionType ( G4HadXSType  val)
inline

Definition at line 301 of file G4HadronicProcess.hh.

302{
303 fXSType = val;
304}

◆ SetEnergyMomentumCheckLevels()

void G4HadronicProcess::SetEnergyMomentumCheckLevels ( G4double  relativeLevel,
G4double  absoluteLevel 
)
inline

Definition at line 322 of file G4HadronicProcess.hh.

324{
325 epCheckLevels.first = relativeLevel;
326 epCheckLevels.second = absoluteLevel;
327 levelsSetByProcess = true;
328}

Referenced by G4HadronicProcessStore::SetProcessAbsLevel(), and G4HadronicProcessStore::SetProcessRelLevel().

◆ SetEpReportLevel()

void G4HadronicProcess::SetEpReportLevel ( G4int  level)
inline

Definition at line 316 of file G4HadronicProcess.hh.

317{
318 epReportLevel = level;
319}

◆ SetIntegral()

void G4HadronicProcess::SetIntegral ( G4bool  val)
inline

Definition at line 311 of file G4HadronicProcess.hh.

312{
313 useIntegralXS = val;
314}

◆ StartTracking()

void G4HadronicProcess::StartTracking ( G4Track *  track)
overridevirtual

Reimplemented from G4VProcess.

Reimplemented in G4NeutronGeneralProcess.

Definition at line 260 of file G4HadronicProcess.cc.

261{
262 currentMat = nullptr;
263 currentParticle = track->GetDefinition();
264 fDynParticle = track->GetDynamicParticle();
265 theNumberOfInteractionLengthLeft = -1.0;
266}

◆ TwoPeaksXS()

std::vector< G4TwoPeaksHadXS * > * G4HadronicProcess::TwoPeaksXS ( ) const
inline

Definition at line 343 of file G4HadronicProcess.hh.

344{
345 return fXSpeaks;
346}

Referenced by BuildPhysicsTable().

Member Data Documentation

◆ aScaleFactor

G4double G4HadronicProcess::aScaleFactor = 1.0
protected

Definition at line 224 of file G4HadronicProcess.hh.

Referenced by BiasCrossSectionByFactor(), CrossSectionFactor(), GetMeanFreePath(), PostStepDoIt(), and G4HadronElasticProcess::PostStepDoIt().

◆ epReportLevel

G4long G4HadronicProcess::epReportLevel = 0
protected

Definition at line 227 of file G4HadronicProcess.hh.

Referenced by G4HadronStoppingProcess::AtRestDoIt(), CheckEnergyMomentumConservation(), CheckResult(), PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), and SetEpReportLevel().

◆ fWeight

G4double G4HadronicProcess::fWeight = 1.0
protected

Definition at line 223 of file G4HadronicProcess.hh.

Referenced by FillResult(), and PostStepDoIt().

◆ fXSType

G4HadXSType G4HadronicProcess::fXSType = fHadNoIntegral
protected

Definition at line 229 of file G4HadronicProcess.hh.

Referenced by BuildPhysicsTable(), CrossSectionType(), PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), and SetCrossSectionType().

◆ mfpKinEnergy

G4double G4HadronicProcess::mfpKinEnergy = DBL_MAX
protected

Definition at line 226 of file G4HadronicProcess.hh.

Referenced by PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), and PostStepGetPhysicalInteractionLength().

◆ theCrossSectionDataStore

G4CrossSectionDataStore* G4HadronicProcess::theCrossSectionDataStore
protected

Definition at line 221 of file G4HadronicProcess.hh.

Referenced by AddDataSet(), BuildPhysicsTable(), ComputeCrossSection(), DumpPhysicsTable(), GetCrossSectionDataStore(), GetElementCrossSection(), GetMeanFreePath(), PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), and ~G4HadronicProcess().

◆ theLastCrossSection

G4double G4HadronicProcess::theLastCrossSection = 0.0
protected

Definition at line 225 of file G4HadronicProcess.hh.

Referenced by GetLastCrossSection(), PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), and PostStepGetPhysicalInteractionLength().

◆ thePro

G4HadProjectile G4HadronicProcess::thePro
protected

Definition at line 218 of file G4HadronicProcess.hh.

Referenced by G4HadronStoppingProcess::AtRestDoIt(), G4NeutrinoElectronProcess::PostStepDoIt(), PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4MuNeutrinoNucleusProcess::PostStepDoIt(), and G4TauNeutrinoNucleusProcess::PostStepDoIt().

◆ theTotalResult

G4ParticleChange* G4HadronicProcess::theTotalResult
protected

Definition at line 220 of file G4HadronicProcess.hh.

Referenced by G4HadronStoppingProcess::AtRestDoIt(), CheckEnergyMomentumConservation(), FillResult(), G4NeutronGeneralProcess::PostStepDoIt(), G4NeutrinoElectronProcess::PostStepDoIt(), PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), G4MuNeutrinoNucleusProcess::PostStepDoIt(), G4TauNeutrinoNucleusProcess::PostStepDoIt(), and ~G4HadronicProcess().

The documentation for this class was generated from the following files: