Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Member Functions | Protected Attributes
G4AdjointAlongStepWeightCorrection Class Reference

#include <G4AdjointAlongStepWeightCorrection.hh>

Inheritance diagram for G4AdjointAlongStepWeightCorrection:
G4VContinuousProcess G4VProcess

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
 

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::G4AdjointAlongStepWeightCorrection ( const G4String &  name = "ContinuousWeightCorrection",
G4ProcessType  type = fElectromagnetic 
)

Definition at line 37 of file G4AdjointAlongStepWeightCorrection.cc.

References fParticleChange.

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

Definition at line 48 of file G4AdjointAlongStepWeightCorrection.cc.

References fParticleChange.

49 {delete fParticleChange;
50 }

Member Function Documentation

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

Reimplemented from G4VContinuousProcess.

Definition at line 66 of file G4AdjointAlongStepWeightCorrection.cc.

References fParticleChange, G4AdjointCSManager::GetAdjointCSManager(), G4AdjointCSManager::GetContinuousWeightCorrection(), G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4StepPoint::GetKineticEnergy(), G4Step::GetPostStepPoint(), G4Step::GetStepLength(), G4StepPoint::GetWeight(), G4ParticleChange::Initialize(), G4VParticleChange::ProposeParentWeight(), G4VParticleChange::SetParentWeightByProcess(), and G4VParticleChange::SetSecondaryWeightByProcess().

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 }
G4StepPoint::GetWeight
G4double GetWeight() const
G4Step::GetStepLength
G4double GetStepLength() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4VParticleChange::ProposeParentWeight
void ProposeParentWeight(G4double finalWeight)
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4VParticleChange::SetSecondaryWeightByProcess
void SetSecondaryWeightByProcess(G4bool)
G4VParticleChange::SetParentWeightByProcess
void SetParentWeightByProcess(G4bool)
Definition: G4VParticleChange.cc:522
G4ParticleChange::Initialize
virtual void Initialize(const G4Track &)
Definition: G4ParticleChange.cc:228
G4AdjointCSManager::GetContinuousWeightCorrection
G4double GetContinuousWeightCorrection(G4ParticleDefinition *aPartDef, G4double PreStepEkin, G4double AfterStepEkin, const G4MaterialCutsCouple *aCouple, G4double step_length)
Definition: G4AdjointCSManager.cc:502
G4Step::GetPostStepPoint
G4StepPoint * GetPostStepPoint() const
G4StepPoint::GetKineticEnergy
G4double GetKineticEnergy() const
G4double
double G4double
Definition: G4Types.hh:76
G4AdjointCSManager::GetAdjointCSManager
static G4AdjointCSManager * GetAdjointCSManager()
Definition: G4AdjointCSManager.cc:61
void G4AdjointAlongStepWeightCorrection::BuildPhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 61 of file G4AdjointAlongStepWeightCorrection.cc.

62 {;
63 }
G4double G4AdjointAlongStepWeightCorrection::GetContinuousStepLimit ( const G4Track &  track,
G4double  previousStepSize,
G4double  currentMinimumStep,
G4double &  currentSafety 
)
protectedvirtual

Implements G4VContinuousProcess.

Definition at line 120 of file G4AdjointAlongStepWeightCorrection.cc.

References DBL_MAX, G4Track::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), and test::x.

122 {
123  G4double x = DBL_MAX;
124  DefineMaterial(track.GetMaterialCutsCouple());
125  preStepKinEnergy = track.GetKineticEnergy();
126  return x;
127 }
G4Track::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
G4Track::GetKineticEnergy
G4double GetKineticEnergy() const
G4double
double G4double
Definition: G4Types.hh:76
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
void G4AdjointAlongStepWeightCorrection::PreparePhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 53 of file G4AdjointAlongStepWeightCorrection.cc.

55 {
56 ;
57 }

Field Documentation

G4ParticleChange* G4AdjointAlongStepWeightCorrection::fParticleChange
protected

Definition at line 112 of file G4AdjointAlongStepWeightCorrection.hh.

Referenced by AlongStepDoIt(), G4AdjointAlongStepWeightCorrection(), and ~G4AdjointAlongStepWeightCorrection().

The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:56:10 for Geant4.10 by   doxygen 1.8.7