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

#include <G4UnknownDecay.hh>

+ Inheritance diagram for G4UnknownDecay:

Public Member Functions

 G4UnknownDecay (const G4String &processName="UnknownDecay")
 
virtual ~G4UnknownDecay ()
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &aTrack, const G4Step &aStep)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual G4bool IsApplicable (const G4ParticleDefinition &)
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
- Public Member Functions inherited from G4VDiscreteProcess
 G4VDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)
 
 G4VDiscreteProcess (G4VDiscreteProcess &)
 
virtual ~G4VDiscreteProcess ()
 
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 &)
 
- Public Member Functions inherited from G4VProcess
 G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
 
 G4VProcess (const G4VProcess &right)
 
virtual ~G4VProcess ()
 
G4int operator== (const G4VProcess &right) const
 
G4int 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 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
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 

Protected Member Functions

virtual G4VParticleChange * DecayIt (const G4Track &aTrack, const G4Step &aStep)
 
virtual G4double GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4double GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManager * aProcessManager
 
G4VParticleChange * pParticleChange
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft
 
G4double currentInteractionLength
 
G4double theInitialNumberOfInteractionLength
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType
 
G4int theProcessSubType
 
G4double thePILfactor
 
G4bool enableAtRestDoIt
 
G4bool enableAlongStepDoIt
 
G4bool enablePostStepDoIt
 
G4int verboseLevel
 

Detailed Description

Definition at line 44 of file G4UnknownDecay.hh.

Constructor & Destructor Documentation

◆ G4UnknownDecay()

G4UnknownDecay::G4UnknownDecay ( const G4String &  processName = "UnknownDecay")

Definition at line 47 of file G4UnknownDecay.cc.

48 :G4VDiscreteProcess(processName, fDecay),
49 verboseLevel(1),
50 HighestValue(20.0)
51{
52 // set Process Sub Type
53 SetProcessSubType(static_cast<int>(DECAY_Unknown));
54
55#ifdef G4VERBOSE
56 if (GetVerboseLevel()>1) {
57 G4cout << "G4UnknownDecay constructor " << " Name:" << processName << G4endl;
58 }
59#endif
60 pParticleChange = &fParticleChangeForDecay;
61}
DECAY_Unknown
@ DECAY_Unknown
Definition: G4DecayProcessType.hh:49
fDecay
@ fDecay
Definition: G4ProcessType.hh:51
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
G4UnknownDecay::GetVerboseLevel
G4int GetVerboseLevel() const
Definition: G4UnknownDecay.hh:154
G4VProcess::SetProcessSubType
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:403
G4VProcess::pParticleChange
G4VParticleChange * pParticleChange
Definition: G4VProcess.hh:283

◆ ~G4UnknownDecay()

G4UnknownDecay::~G4UnknownDecay ( )
virtual

Definition at line 63 of file G4UnknownDecay.cc.

64{
65}

Member Function Documentation

◆ BuildPhysicsTable()

void G4UnknownDecay::BuildPhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 78 of file G4UnknownDecay.cc.

79{
80 return;
81}

◆ DecayIt()

G4VParticleChange * G4UnknownDecay::DecayIt ( const G4Track &  aTrack,
const G4Step &  aStep 
)
protectedvirtual

Definition at line 83 of file G4UnknownDecay.cc.

84{
85 // The DecayIt() method returns by pointer a particle-change object.
86 // Units are expressed in GEANT4 internal units.
87
88 // Initialize ParticleChange
89 // all members of G4VParticleChange are set to equal to
90 // corresponding member in G4Track
91 fParticleChangeForDecay.Initialize(aTrack);
92
93 // get particle
94 const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
95
96 //check if thePreAssignedDecayProducts exists
97 const G4DecayProducts* o_products = (aParticle->GetPreAssignedDecayProducts());
98 G4bool isPreAssigned = (o_products != 0);
99 G4DecayProducts* products = 0;
100
101 if (!isPreAssigned ){
102 fParticleChangeForDecay.SetNumberOfSecondaries(0);
103 // Kill the parent particle
104 fParticleChangeForDecay.ProposeTrackStatus( fStopAndKill ) ;
105 fParticleChangeForDecay.ProposeLocalEnergyDeposit(0.0);
106
107 ClearNumberOfInteractionLengthLeft();
108 return &fParticleChangeForDecay ;
109 }
110
111 // copy decay products
112 products = new G4DecayProducts(*o_products);
113
114 // get parent particle information ...................................
115 G4double ParentEnergy = aParticle->GetTotalEnergy();
116 G4double ParentMass = aParticle->GetMass();
117 if (ParentEnergy < ParentMass) {
118 ParentEnergy = ParentMass;
119#ifdef G4VERBOSE
120 if (GetVerboseLevel()>1) {
121 G4cout << "G4UnknownDecay::DoIt : Total Energy is less than its mass" << G4endl;
122 G4cout << " Particle: " << aParticle->GetDefinition()->GetParticleName();
123 G4cout << " Energy:" << ParentEnergy/MeV << "[MeV]";
124 G4cout << " Mass:" << ParentMass/MeV << "[MeV]";
125 G4cout << G4endl;
126 }
127#endif
128 }
129 G4ThreeVector ParentDirection(aParticle->GetMomentumDirection());
130
131 G4double energyDeposit = 0.0;
132 G4double finalGlobalTime = aTrack.GetGlobalTime();
133 //boost all decay products to laboratory frame
134 //if the particle has traveled
135 if(aParticle->GetPreAssignedDecayProperTime()>0.) {
136 products->Boost( ParentEnergy, ParentDirection);
137 }
138
139 //add products in fParticleChangeForDecay
140 G4int numberOfSecondaries = products->entries();
141 fParticleChangeForDecay.SetNumberOfSecondaries(numberOfSecondaries);
142#ifdef G4VERBOSE
143 if (GetVerboseLevel()>1) {
144 G4cout << "G4UnknownDecay::DoIt : Decay vertex :";
145 G4cout << " Time: " << finalGlobalTime/ns << "[ns]";
146 G4cout << " X:" << (aTrack.GetPosition()).x() /cm << "[cm]";
147 G4cout << " Y:" << (aTrack.GetPosition()).y() /cm << "[cm]";
148 G4cout << " Z:" << (aTrack.GetPosition()).z() /cm << "[cm]";
149 G4cout << G4endl;
150 G4cout << "G4UnknownDecay::DoIt : decay products in Lab. Frame" << G4endl;
151 products->DumpInfo();
152 }
153#endif
154 G4int index;
155 G4ThreeVector currentPosition;
156 const G4TouchableHandle thand = aTrack.GetTouchableHandle();
157 for (index=0; index < numberOfSecondaries; index++)
158 {
159 // get current position of the track
160 currentPosition = aTrack.GetPosition();
161 // create a new track object
162 G4Track* secondary = new G4Track( products->PopProducts(),
163 finalGlobalTime ,
164 currentPosition );
165 // switch on good for tracking flag
166 secondary->SetGoodForTrackingFlag();
167 secondary->SetTouchableHandle(thand);
168 // add the secondary track in the List
169 fParticleChangeForDecay.AddSecondary(secondary);
170 }
171 delete products;
172
173 // Kill the parent particle
174 fParticleChangeForDecay.ProposeTrackStatus( fStopAndKill ) ;
175 fParticleChangeForDecay.ProposeLocalEnergyDeposit(energyDeposit);
176 fParticleChangeForDecay.ProposeGlobalTime( finalGlobalTime );
177 // reset NumberOfInteractionLengthLeft
178 ClearNumberOfInteractionLengthLeft();
179
180 return &fParticleChangeForDecay ;
181}
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:56
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4bool
bool G4bool
Definition: G4Types.hh:67
G4DecayProducts::DumpInfo
void DumpInfo() const
Definition: G4DecayProducts.cc:311
G4DecayProducts::entries
G4int entries() const
Definition: G4DecayProducts.hh:86
G4DecayProducts::PopProducts
G4DynamicParticle * PopProducts()
Definition: G4DecayProducts.cc:141
G4DecayProducts::Boost
void Boost(G4double totalEnergy, const G4ThreeVector &momentumDirection)
Definition: G4DecayProducts.cc:176
G4DynamicParticle::GetMass
G4double GetMass() const
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetPreAssignedDecayProducts
const G4DecayProducts * GetPreAssignedDecayProducts() const
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4DynamicParticle::GetTotalEnergy
G4double GetTotalEnergy() const
G4DynamicParticle::GetPreAssignedDecayProperTime
G4double GetPreAssignedDecayProperTime() const
G4ParticleChangeForDecay::Initialize
virtual void Initialize(const G4Track &)
Definition: G4ParticleChangeForDecay.cc:124
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:115
G4Track::GetPosition
const G4ThreeVector & GetPosition() const
G4Track::SetTouchableHandle
void SetTouchableHandle(const G4TouchableHandle &apValue)
G4Track::GetGlobalTime
G4double GetGlobalTime() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4Track::GetTouchableHandle
const G4TouchableHandle & GetTouchableHandle() const
G4Track::SetGoodForTrackingFlag
void SetGoodForTrackingFlag(G4bool value=true)
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4VParticleChange::AddSecondary
void AddSecondary(G4Track *aSecondary)
Definition: G4VParticleChange.cc:171
G4VParticleChange::ProposeLocalEnergyDeposit
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4VParticleChange::SetNumberOfSecondaries
void SetNumberOfSecondaries(G4int totSecondaries)
G4VProcess::ClearNumberOfInteractionLengthLeft
void ClearNumberOfInteractionLengthLeft()
Definition: G4VProcess.hh:418
ns
#define ns
Definition: xmlparse.cc:597

Referenced by PostStepDoIt().

◆ GetMeanFreePath()

G4double G4UnknownDecay::GetMeanFreePath ( const G4Track &  aTrack,
G4double  previousStepSize,
G4ForceCondition *  condition 
)
protectedvirtual

Implements G4VDiscreteProcess.

Definition at line 73 of file G4UnknownDecay.cc.

74{
75 return DBL_MIN;
76}
DBL_MIN
#define DBL_MIN
Definition: templates.hh:75

◆ GetVerboseLevel()

G4int G4UnknownDecay::GetVerboseLevel ( ) const
inline

Definition at line 154 of file G4UnknownDecay.hh.

154{ return verboseLevel; }

Referenced by DecayIt(), and G4UnknownDecay().

◆ IsApplicable()

G4bool G4UnknownDecay::IsApplicable ( const G4ParticleDefinition &  aParticleType)
virtual

Reimplemented from G4VProcess.

Definition at line 67 of file G4UnknownDecay.cc.

68{
69 if(aParticleType.GetParticleName()=="unknown") return true;
70 return false;
71}

◆ PostStepDoIt()

G4VParticleChange * G4UnknownDecay::PostStepDoIt ( const G4Track &  aTrack,
const G4Step &  aStep 
)
inlinevirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 157 of file G4UnknownDecay.hh.

161{
162 return DecayIt(aTrack, aStep);
163}
G4UnknownDecay::DecayIt
virtual G4VParticleChange * DecayIt(const G4Track &aTrack, const G4Step &aStep)
Definition: G4UnknownDecay.cc:83

◆ PostStepGetPhysicalInteractionLength()

G4double G4UnknownDecay::PostStepGetPhysicalInteractionLength ( const G4Track &  track,
G4double  previousStepSize,
G4ForceCondition *  condition 
)
inlinevirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 126 of file G4UnknownDecay.hh.

131{
132 // pre-assigned UnknownDecay time
133 G4double pTime = track.GetDynamicParticle()->GetPreAssignedDecayProperTime();
134
135 if (pTime < 0.) pTime = DBL_MIN;
136
137 // condition is set to "Not Forced"
138 *condition = NotForced;
139
140 // reminder proper time
141 G4double remainder = pTime - track.GetProperTime();
142 if (remainder <= 0.0) remainder = DBL_MIN;
143
144 // use pre-assigned Decay time to determine PIL
145 //return GetMeanFreePath(track, previousStepSize, condition);
146 return remainder*CLHEP::c_light;
147
148}
condition
G4double condition(const G4ErrorSymMatrix &m)
NotForced
@ NotForced
Definition: G4ForceCondition.hh:56
G4Track::GetProperTime
G4double GetProperTime() const

◆ SetVerboseLevel()

void G4UnknownDecay::SetVerboseLevel ( G4int  value)
inline

Definition at line 151 of file G4UnknownDecay.hh.

151{ verboseLevel = value; }
The documentation for this class was generated from the following files: