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

#include <G4OpMieHG.hh>

Inheritance diagram for G4OpMieHG:
G4VDiscreteProcess G4VProcess

Public Member Functions

 G4OpMieHG (const G4String &processName="OpMieHG", G4ProcessType type=fOptical)
 
 ~G4OpMieHG ()
 
G4bool IsApplicable (const G4ParticleDefinition &aParticleType)
 
G4double GetMeanFreePath (const G4Track &aTrack, G4double, G4ForceCondition *)
 
G4VParticleChange * PostStepDoIt (const G4Track &aTrack, const G4Step &aStep)
 
- 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 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
 
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 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
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcess * GetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 
- 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 48 of file G4OpMieHG.hh.

Constructor & Destructor Documentation

G4OpMieHG::G4OpMieHG ( const G4String &  processName = "OpMieHG",
G4ProcessType  type = fOptical 
)

Definition at line 47 of file G4OpMieHG.cc.

References fOpMieHG, G4cout, G4endl, G4VProcess::GetProcessName(), G4VProcess::SetProcessSubType(), and G4VProcess::verboseLevel.

48  : G4VDiscreteProcess(processName, type)
49 {
50  if (verboseLevel>0) {
51  G4cout << GetProcessName() << " is created " << G4endl;
52  }
53 
54  SetProcessSubType(fOpMieHG);
55 }
G4VProcess::verboseLevel
G4int verboseLevel
Definition: G4VProcess.hh:368
G4cout
G4GLOB_DLL std::ostream G4cout
G4VProcess::SetProcessSubType
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
G4VProcess::GetProcessName
const G4String & GetProcessName() const
Definition: G4VProcess.hh:408
G4endl
#define G4endl
Definition: G4ios.hh:61
fOpMieHG
Definition: G4OpProcessSubType.hh:50
G4OpMieHG::~G4OpMieHG ( )

Definition at line 57 of file G4OpMieHG.cc.

57 {}

Member Function Documentation

G4double G4OpMieHG::GetMeanFreePath ( const G4Track &  aTrack,
G4double  ,
G4ForceCondition *   
)
virtual

Implements G4VDiscreteProcess.

Definition at line 159 of file G4OpMieHG.cc.

References DBL_MAX, G4Track::GetDynamicParticle(), G4Track::GetMaterial(), G4Material::GetMaterialPropertiesTable(), G4MaterialPropertiesTable::GetProperty(), and G4DynamicParticle::GetTotalEnergy().

162 {
163  const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
164  const G4Material* aMaterial = aTrack.GetMaterial();
165 
166  G4double thePhotonEnergy = aParticle->GetTotalEnergy();
167 
168  G4double AttenuationLength = DBL_MAX;
169 
170  G4MaterialPropertiesTable* aMaterialPropertyTable =
171  aMaterial->GetMaterialPropertiesTable();
172 
173  if (aMaterialPropertyTable) {
174  G4MaterialPropertyVector* AttenuationLengthVector =
175  aMaterialPropertyTable->GetProperty("MIEHG");
176  if (AttenuationLengthVector) {
177  AttenuationLength = AttenuationLengthVector ->
178  Value(thePhotonEnergy);
179  } else {
180 // G4cout << "No Mie scattering length specified" << G4endl;
181  }
182  } else {
183 // G4cout << "No Mie scattering length specified" << G4endl;
184  }
185 
186 // G4cout << thePhotonEnergy/GeV << " \t" << AttenuationLength/m << G4endl;
187 
188  return AttenuationLength;
189 }
G4MaterialPropertiesTable::GetProperty
G4MaterialPropertyVector * GetProperty(const char *key)
G4DynamicParticle::GetTotalEnergy
G4double GetTotalEnergy() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4Track::GetMaterial
G4Material * GetMaterial() const
G4Material::GetMaterialPropertiesTable
G4MaterialPropertiesTable * GetMaterialPropertiesTable() const
Definition: G4Material.hh:250
G4double
double G4double
Definition: G4Types.hh:76
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4bool G4OpMieHG::IsApplicable ( const G4ParticleDefinition &  aParticleType)
inlinevirtual

Reimplemented from G4VProcess.

Definition at line 92 of file G4OpMieHG.hh.

References G4OpticalPhoton::OpticalPhoton().

93 {
94  return ( &aParticleType == G4OpticalPhoton::OpticalPhoton() );
95 }
G4OpticalPhoton::OpticalPhoton
static G4OpticalPhoton * OpticalPhoton()
Definition: G4OpticalPhoton.cc:87
G4VParticleChange * G4OpMieHG::PostStepDoIt ( const G4Track &  aTrack,
const G4Step &  aStep 
)
virtual

Reimplemented from G4VDiscreteProcess.

Definition at line 67 of file G4OpMieHG.cc.

References G4VProcess::aParticleChange, CLHEP::Hep3Vector::dot(), G4cout, G4endl, G4UniformRand, G4MaterialPropertiesTable::GetConstProperty(), G4Track::GetDynamicParticle(), G4Track::GetMaterial(), G4Material::GetMaterialPropertiesTable(), G4DynamicParticle::GetMomentumDirection(), G4ParticleChange::GetMomentumDirection(), G4ParticleChange::GetPolarization(), G4DynamicParticle::GetPolarization(), G4ParticleChange::Initialize(), CLHEP::Hep3Vector::mag(), python.hepunit::pi, G4VDiscreteProcess::PostStepDoIt(), G4ParticleChange::ProposeMomentumDirection(), G4ParticleChange::ProposePolarization(), CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), python.hepunit::twopi, CLHEP::Hep3Vector::unit(), and G4VProcess::verboseLevel.

68 {
69  aParticleChange.Initialize(aTrack);
70 
71  const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
72  const G4Material* aMaterial = aTrack.GetMaterial();
73  G4MaterialPropertiesTable* aMaterialPropertyTable =
74  aMaterial->GetMaterialPropertiesTable();
75 
76  G4double forward_g =
77  aMaterialPropertyTable->GetConstProperty("MIEHG_FORWARD");
78  G4double backward_g =
79  aMaterialPropertyTable->GetConstProperty("MIEHG_BACKWARD");
80  G4double ForwardRatio =
81  aMaterialPropertyTable->GetConstProperty("MIEHG_FORWARD_RATIO");
82 
83  if (verboseLevel>0) {
84  G4cout << "MIE Scattering Photon!" << G4endl;
85  G4cout << "MIE Old Momentum Direction: "
86  << aParticle->GetMomentumDirection() << G4endl;
87  G4cout << "MIE Old Polarization: "
88  << aParticle->GetPolarization() << G4endl;
89  }
90 
91  G4double gg;
92  G4int direction;
93  if (G4UniformRand()<=ForwardRatio){
94  gg = forward_g;
95  direction = 1;
96  } else {
97  gg = backward_g;
98  direction = -1;
99  }
100 
101  G4double r = G4UniformRand();
102 
103  G4double Theta;
104  //sample the direction
105  if (gg!=0) {
106  Theta = std::acos(2*r*(1+gg)*(1+gg)*(1-gg+gg*r)/((1-gg+2*gg*r)*(1-gg+2*gg*r)) -1);
107  } else {
108  Theta = std::acos(2*r-1.);
109  }
110  G4double Phi = G4UniformRand()*2*pi;
111 
112  if (direction==-1) Theta = pi - Theta; //backward scattering
113 
114  G4ThreeVector NewMomentumDirection, OldMomentumDirection;
115  G4ThreeVector OldPolarization, NewPolarization;
116 
117  NewMomentumDirection.set
118  (std::sin(Theta)*std::cos(Phi), std::sin(Theta)*std::sin(Phi), std::cos(Theta));
119  OldMomentumDirection = aParticle->GetMomentumDirection();
120  NewMomentumDirection.rotateUz(OldMomentumDirection);
121  NewMomentumDirection = NewMomentumDirection.unit();
122 
123  OldPolarization = aParticle->GetPolarization();
124  G4double constant = -1./NewMomentumDirection.dot(OldPolarization);
125 
126  NewPolarization = NewMomentumDirection + constant*OldPolarization;
127  NewPolarization = NewPolarization.unit();
128 
129  if (NewPolarization.mag()==0) {
130  r = G4UniformRand()*twopi;
131  NewPolarization.set(std::cos(r),std::sin(r),0.);
132  NewPolarization.rotateUz(NewMomentumDirection);
133  } else {
134  // There are two directions which perpendicular
135  // new momentum direction
136  if (G4UniformRand() < 0.5) NewPolarization = -NewPolarization;
137  }
138 
139  aParticleChange.ProposePolarization(NewPolarization);
140  aParticleChange.ProposeMomentumDirection(NewMomentumDirection);
141 
142  if (verboseLevel>0) {
143  G4cout << "MIE New Polarization: "
144  << NewPolarization << G4endl;
145  G4cout << "MIE Polarization Change: "
146  << *(aParticleChange.GetPolarization()) << G4endl;
147  G4cout << "MIE New Momentum Direction: "
148  << NewMomentumDirection << G4endl;
149  G4cout << "MIE Momentum Change: "
150  << *(aParticleChange.GetMomentumDirection()) << G4endl;
151  }
152 
153  return G4VDiscreteProcess::PostStepDoIt(aTrack, aStep);
154 }
CLHEP::Hep3Vector::set
void set(double x, double y, double z)
G4VProcess::verboseLevel
G4int verboseLevel
Definition: G4VProcess.hh:368
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
CLHEP::Hep3Vector::dot
double dot(const Hep3Vector &) const
G4ParticleChange::ProposePolarization
void ProposePolarization(G4double Px, G4double Py, G4double Pz)
G4ParticleChange::GetMomentumDirection
const G4ThreeVector * GetMomentumDirection() const
G4int
int G4int
Definition: G4Types.hh:78
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:87
G4cout
G4GLOB_DLL std::ostream G4cout
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
CLHEP::Hep3Vector::rotateUz
Hep3Vector & rotateUz(const Hep3Vector &)
Definition: ThreeVector.cc:72
G4Track::GetMaterial
G4Material * GetMaterial() const
G4ParticleChange::Initialize
virtual void Initialize(const G4Track &)
Definition: G4ParticleChange.cc:228
G4Material::GetMaterialPropertiesTable
G4MaterialPropertiesTable * GetMaterialPropertiesTable() const
Definition: G4Material.hh:250
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
G4DynamicParticle::GetPolarization
const G4ThreeVector & GetPolarization() const
G4VProcess::aParticleChange
G4ParticleChange aParticleChange
Definition: G4VProcess.hh:289
G4endl
#define G4endl
Definition: G4ios.hh:61
G4ParticleChange::ProposeMomentumDirection
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4double
double G4double
Definition: G4Types.hh:76
G4ParticleChange::GetPolarization
const G4ThreeVector * GetPolarization() const
G4MaterialPropertiesTable::GetConstProperty
G4double GetConstProperty(const char *key)
CLHEP::Hep3Vector::mag
double mag() const
G4VDiscreteProcess::PostStepDoIt
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
Definition: G4VDiscreteProcess.cc:112
The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:56:51 for Geant4.10 by   doxygen 1.8.7