Geant4  10.02.p03
Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | List of all members
G4AdjointAlongStepWeightCorrection Class Reference

#include <G4AdjointAlongStepWeightCorrection.hh>

Inheritance diagram for G4AdjointAlongStepWeightCorrection:
Collaboration diagram for G4AdjointAlongStepWeightCorrection:

Public Member Functions

 G4AdjointAlongStepWeightCorrection (const G4String &name="ContinuousWeightCorrection", G4ProcessType type=fElectromagnetic)
 
virtual ~G4AdjointAlongStepWeightCorrection ()
 
void PreparePhysicsTable (const G4ParticleDefinition &)
 
void BuildPhysicsTable (const G4ParticleDefinition &)
 
G4VParticleChange * AlongStepDoIt (const G4Track &, const G4Step &)
 
- Public Member Functions inherited from G4VContinuousProcess
 G4VContinuousProcess (const G4String &, G4ProcessType aType=fNotDefined)
 
 G4VContinuousProcess (G4VContinuousProcess &)
 
virtual ~G4VContinuousProcess ()
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &, G4double, G4ForceCondition *)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4VParticleChange * PostStepDoIt (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
 
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 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
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcess * GetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Protected Member Functions

virtual G4double GetContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)
 
- Protected Member Functions inherited from G4VContinuousProcess
void SetGPILSelection (G4GPILSelection selection)
 
G4GPILSelection GetGPILSelection () const
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Protected Attributes

G4ParticleChange * fParticleChange
 
- 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
 

Private Member Functions

void DefineMaterial (const G4MaterialCutsCouple *couple)
 
 G4AdjointAlongStepWeightCorrection (G4AdjointAlongStepWeightCorrection &)
 
G4AdjointAlongStepWeightCorrection & operator= (const G4AdjointAlongStepWeightCorrection &right)
 

Private Attributes

const G4Material * currentMaterial
 
const G4MaterialCutsCouple * currentCouple
 
size_t currentMaterialIndex
 
G4double preStepKinEnergy
 

Additional Inherited Members

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

Detailed Description

Definition at line 70 of file G4AdjointAlongStepWeightCorrection.hh.

Constructor & Destructor Documentation

◆ G4AdjointAlongStepWeightCorrection() [1/2]

G4AdjointAlongStepWeightCorrection::G4AdjointAlongStepWeightCorrection ( const G4String &  name = "ContinuousWeightCorrection",
G4ProcessType  type = fElectromagnetic 
)

Definition at line 37 of file G4AdjointAlongStepWeightCorrection.cc.

38  : G4VContinuousProcess(name, type)
39 {fParticleChange = new G4ParticleChange();
40  currentMaterialIndex=0;
41  preStepKinEnergy=1.;
42  currentCouple=0;
43 }

◆ ~G4AdjointAlongStepWeightCorrection()

G4AdjointAlongStepWeightCorrection::~G4AdjointAlongStepWeightCorrection ( )
virtual

Definition at line 48 of file G4AdjointAlongStepWeightCorrection.cc.

49 {delete fParticleChange;
50 }

◆ G4AdjointAlongStepWeightCorrection() [2/2]

G4AdjointAlongStepWeightCorrection::G4AdjointAlongStepWeightCorrection ( G4AdjointAlongStepWeightCorrection &  )
private

Member Function Documentation

◆ AlongStepDoIt()

G4VParticleChange * G4AdjointAlongStepWeightCorrection::AlongStepDoIt ( const G4Track &  track,
const G4Step &  step 
)
virtual

Reimplemented from G4VContinuousProcess.

Definition at line 66 of file G4AdjointAlongStepWeightCorrection.cc.

68 {
69 
70  fParticleChange->Initialize(track);
71 
72  // Get the actual (true) Step length
73  //----------------------------------
74  G4double length = step.GetStepLength();
75 
76 
77  G4double Tkin = step.GetPostStepPoint()->GetKineticEnergy();
78  G4ParticleDefinition* thePartDef= const_cast<G4ParticleDefinition*> (track.GetDynamicParticle()->GetDefinition());
79  G4double weight_correction=G4AdjointCSManager::GetAdjointCSManager()->GetContinuousWeightCorrection(thePartDef,
80  preStepKinEnergy,Tkin, currentCouple,length);
81 
82 
83 
84 
85  //Caution!!!
86  // It is important to select the weight of the post_step_point
87  // as the current weight and not the weight of the track, as t
88  // the weight of the track is changed after having applied all
89  // the along_step_do_it.
90 
91  // G4double new_weight=weight_correction*track.GetWeight(); //old
92  G4double new_weight=weight_correction*step.GetPostStepPoint()->GetWeight();
93 
94 
95 
96  //if (weight_correction >2.) new_weight=1.e-300;
97 
98 
99  //The following test check for zero weight.
100  //This happens after weight correction of gamma for photo electric effect.
101  //When the new weight is 0 it will be later on consider as nan by G4.
102  //Therefore we do put a lower limit of 1.e-300. for new_weight
103  //Correction by L.Desorgher on 15 July 2009
104  if (new_weight==0 || (new_weight<=0 && new_weight>0)){
105  //G4cout<<new_weight<<'\t'<<weight_correction<<'\t'<<track.GetWeight()<<G4endl;
106  new_weight=1.e-300;
107  }
108 
109  //G4cout<<new_weight<<'\t'<<weight_correction<<'\t'<<track.GetWeight()<<G4endl;
110  fParticleChange->SetParentWeightByProcess(false);
111  fParticleChange->SetSecondaryWeightByProcess(false);
112  fParticleChange->ProposeParentWeight(new_weight);
113 
114 
115  return fParticleChange;
116 
117 }
G4AdjointCSManager::GetContinuousWeightCorrection
G4double GetContinuousWeightCorrection(G4ParticleDefinition *aPartDef, G4double PreStepEkin, G4double AfterStepEkin, const G4MaterialCutsCouple *aCouple, G4double step_length)
Definition: G4AdjointCSManager.cc:504
G4double
double G4double
Definition: G4Types.hh:76
G4AdjointCSManager::GetAdjointCSManager
static G4AdjointCSManager * GetAdjointCSManager()
Definition: G4AdjointCSManager.cc:61
Here is the call graph for this function:

◆ BuildPhysicsTable()

void G4AdjointAlongStepWeightCorrection::BuildPhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 61 of file G4AdjointAlongStepWeightCorrection.cc.

62 {;
63 }

◆ DefineMaterial()

void G4AdjointAlongStepWeightCorrection::DefineMaterial ( const G4MaterialCutsCouple *  couple)
inlineprivate

Definition at line 124 of file G4AdjointAlongStepWeightCorrection.hh.

126 {
127  if(couple != currentCouple) {
128  currentCouple = couple;
129  currentMaterial = couple->GetMaterial();
130  currentMaterialIndex = couple->GetIndex();
131 
132  }
133 }
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:150
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetContinuousStepLimit()

G4double G4AdjointAlongStepWeightCorrection::GetContinuousStepLimit ( const G4Track &  track,
G4double  previousStepSize,
G4double  currentMinimumStep,
G4double &  currentSafety 
)
protectedvirtual

Implements G4VContinuousProcess.

Definition at line 120 of file G4AdjointAlongStepWeightCorrection.cc.

122 {
123  G4double x = DBL_MAX;
124  DefineMaterial(track.GetMaterialCutsCouple());
125  preStepKinEnergy = track.GetKineticEnergy();
126  return x;
127 }
G4AdjointAlongStepWeightCorrection::DefineMaterial
void DefineMaterial(const G4MaterialCutsCouple *couple)
Definition: G4AdjointAlongStepWeightCorrection.hh:124
G4double
double G4double
Definition: G4Types.hh:76
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
Here is the call graph for this function:

◆ operator=()

G4AdjointAlongStepWeightCorrection& G4AdjointAlongStepWeightCorrection::operator= ( const G4AdjointAlongStepWeightCorrection &  right)
private

◆ PreparePhysicsTable()

void G4AdjointAlongStepWeightCorrection::PreparePhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 53 of file G4AdjointAlongStepWeightCorrection.cc.

55 {
56 ;
57 }

Member Data Documentation

◆ currentCouple

const G4MaterialCutsCouple* G4AdjointAlongStepWeightCorrection::currentCouple
private

Definition at line 118 of file G4AdjointAlongStepWeightCorrection.hh.

◆ currentMaterial

const G4Material* G4AdjointAlongStepWeightCorrection::currentMaterial
private

Definition at line 117 of file G4AdjointAlongStepWeightCorrection.hh.

◆ currentMaterialIndex

size_t G4AdjointAlongStepWeightCorrection::currentMaterialIndex
private

Definition at line 119 of file G4AdjointAlongStepWeightCorrection.hh.

◆ fParticleChange

G4ParticleChange* G4AdjointAlongStepWeightCorrection::fParticleChange
protected

Definition at line 112 of file G4AdjointAlongStepWeightCorrection.hh.

◆ preStepKinEnergy

G4double G4AdjointAlongStepWeightCorrection::preStepKinEnergy
private

Definition at line 120 of file G4AdjointAlongStepWeightCorrection.hh.

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