73 buildAsymmetryTable(true),
74 useAsymmetryTable(true),
78 theAsymmetryTable(NULL)
91 if (theAsymmetryTable) {
92 delete theAsymmetryTable;
101 isInitialised =
true;
108 else if(10 == mType) selectedModel = emModel;
119 G4cout <<
" Total cross sections has a good parametrisation"
120 <<
" from 10 KeV to (100/Z) GeV"
121 <<
"\n Sampling according " << selectedModel->
GetName() <<
" model"
129 if(ss ==
"Klein-Nishina") mType = 0;
130 if(ss ==
"Polarized-Compton") mType = 10;
148 if (theAsymmetryTable && useAsymmetryTable) {
165 if (!VolumeIsPolarized || mfp ==
DBL_MAX)
return mfp;
170 G4cout <<
" Polarization " << GammaPolarization <<
G4endl;
171 G4cout <<
" MaterialPol. " << ElectronPolarization <<
G4endl;
183 asymmetry = aVector->
GetValue(GammaEnergy, isOutRange);
185 G4cout <<
" MaterialIndex " << midx <<
" is out of range \n";
194 G4double pol=ElectronPolarization*GammaDirection0;
196 G4double polProduct = GammaPolarization.
p3() * pol;
197 mfp *= 1. / ( 1. + polProduct * asymmetry );
218 if (theAsymmetryTable && useAsymmetryTable) {
235 if (!VolumeIsPolarized || mfp ==
DBL_MAX)
return mfp;
240 G4cout <<
" Polarization " << GammaPolarization <<
G4endl;
241 G4cout <<
" MaterialPol. " << ElectronPolarization <<
G4endl;
253 asymmetry = aVector->
GetValue(GammaEnergy, isOutRange);
255 G4cout <<
" MaterialIndex " << midx <<
" is out of range \n";
264 G4double pol=ElectronPolarization*GammaDirection0;
266 G4double polProduct = GammaPolarization.
p3() * pol;
267 mfp *= 1. / ( 1. + polProduct * asymmetry );
284 if(buildAsymmetryTable)
295 if(buildAsymmetryTable)
308 for(
size_t i=0; i<numOfCouples; ++i) {
309 if (!theAsymmetryTable)
break;
310 if (theAsymmetryTable->
GetFlag(i)) {
361 lAsymmetry=sigma2/sigma0-1.;
G4double condition(const G4ErrorSymMatrix &m)
const G4ThreeVector & GetPolarization() const
static void SetPhysicsVector(G4PhysicsTable *physTable, size_t idx, G4PhysicsVector *vec)
G4double GetValue(G4double theEnergy, G4bool &isOutRange) const
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
const G4DynamicParticle * GetDynamicParticle() const
void SetBuildTableFlag(G4bool val)
void BuildAsymmetryTable(const G4ParticleDefinition &part)
static G4PolarizationManager * GetInstance()
virtual void PreparePhysicsTable(const G4ParticleDefinition &)
G4double MaxKinEnergy() const
void BuildPhysicsTable(const G4ParticleDefinition &)
G4double PostStepGetPhysicalInteractionLength(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
G4double GetLowEdgeEnergy(size_t binNumber) const
static G4PhysicsTable * PreparePhysicsTable(G4PhysicsTable *physTable)
G4int LambdaBinning() const
virtual void BuildPhysicsTable(const G4ParticleDefinition &)
void SetHighEnergyLimit(G4double)
virtual G4double GetMeanFreePath(const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)
void SetLambdaBinning(G4int nbins)
double precision function energy(A, Z)
size_t CurrentMaterialCutsCoupleIndex() const
G4GLOB_DLL std::ostream G4cout
size_t GetTableSize() const
const G4String & GetName() const
virtual G4double PostStepGetPhysicalInteractionLength(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
const G4ThreeVector & GetMomentumDirection() const
G4double CrossSection(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
void PutValue(size_t index, G4double theValue)
G4double GetMeanFreePath(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
void SetProcessSubType(G4int)
void PreparePhysicsTable(const G4ParticleDefinition &)
G4Material * GetMaterial() const
static G4ProductionCutsTable * GetProductionCutsTable()
G4PolarizedCompton(const G4String &processName="pol-compt", G4ProcessType type=fElectromagnetic)
G4PhysicsVector * LambdaPhysicsVector(const G4MaterialCutsCouple *)
G4double ComputeAsymmetry(G4double energy, const G4MaterialCutsCouple *couple, const G4ParticleDefinition &particle, G4double cut, G4double &tAsymmetry)
void SetMaxKinEnergy(G4double e)
G4LogicalVolume * GetLogicalVolume() const
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
void SetBeamPolarization(const G4ThreeVector &pBeam)
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=0)
void SetSecondaryParticle(const G4ParticleDefinition *p)
bool IsPolarized(G4LogicalVolume *lVol) const
G4VPhysicalVolume * GetVolume() const
virtual void InitialiseProcess(const G4ParticleDefinition *)
static G4Electron * Electron()
virtual ~G4PolarizedCompton()
const G4String & GetName() const
void SetMinKinEnergy(G4double e)
G4bool GetFlag(size_t i) const
void SetLowEnergyLimit(G4double)
G4double MinKinEnergy() const
const G4ThreeVector & GetVolumePolarization(G4LogicalVolume *lVol) const
void SetTargetPolarization(const G4ThreeVector &pTarget)
void SetModel(const G4String &name)