76 buildAsymmetryTable(true),
77 useAsymmetryTable(true),
80 targetPolarization(0.0,0.0,0.0)
139 G4cout <<
" Total cross sections has a good parametrisation"
140 <<
" from 10 KeV to (100/Z) GeV"
149 if(ss ==
"Klein-Nishina") {
mType = 0; }
150 if(ss ==
"Polarized-Compton") {
mType = 10; }
166 G4cout <<
"G4PolarizedCompton::MeanFreePath: " << mfp /
mm <<
" mm " <<
G4endl;
194 G4double prvLength = iLength*satFact;
197 std::max(nLength - previousStepSize/prvLength, 0.0);
202 G4cout <<
"G4PolarizedCompton::PostStepGPIL: "
203 << std::setprecision(8) << x/
mm <<
" mm;" <<
G4endl
204 <<
" unpolarized value: "
205 << std::setprecision(8) << x0/
mm <<
" mm." <<
G4endl;
232 if (VolumeIsPolarized) {
235 G4cout <<
"G4PolarizedCompton::ComputeSaturationFactor: " <<
G4endl;
237 G4cout <<
" Polarization " << GammaPolarization <<
G4endl;
238 G4cout <<
" MaterialPol. " << ElectronPolarization <<
G4endl;
246 if(midx < theAsymmetryTable->size()) {
247 aVector = (*theAsymmetryTable)(midx);
257 G4double pol = ElectronPolarization*GammaDirection0;
258 G4double polProduct = GammaPolarization.
p3() * pol;
259 factor /= (1. + polProduct * asymmetry);
267 ed <<
"Problem with asymmetry table: material index " << midx
268 <<
" is out of range or the table is not filled";
269 G4Exception(
"G4PolarizedComptonModel::ComputeSaturationFactor",
"em0048",
286 if(masterProcess && masterProcess !=
this) { isMaster =
false; }
311 for(
size_t i=0; i<numOfCouples; ++i) {
325 for (
G4int j = 0; j <= nbins; ++j ) {
365 lAsymmetry = sigma2/sigma0-1.;
G4double Energy(size_t index) const
static void SetPhysicsVector(G4PhysicsTable *physTable, size_t idx, G4PhysicsVector *vec)
void SetModel(const G4String &name)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
CLHEP::Hep3Vector G4ThreeVector
std::ostringstream G4ExceptionDescription
G4LogicalVolume * GetLogicalVolume() const
static G4PhysicsTable * PreparePhysicsTable(G4PhysicsTable *physTable)
static constexpr double MeV
void SetHighEnergyLimit(G4double)
G4PolarizedCompton(const G4String &processName="pol-compt", G4ProcessType type=fElectromagnetic)
static constexpr double mm
G4PolarizedComptonModel * emModel
const G4ThreeVector & GetMomentumDirection() const
const G4VProcess * GetMasterProcess() const
void SetBuildTableFlag(G4bool val)
G4double MaxKinEnergy() const
G4VEmModel * EmModel(size_t index=0) const
const G4ThreeVector & GetPolarization() const
static const G4double emax
G4double Value(G4double theEnergy, size_t &lastidx) const
G4double condition(const G4ErrorSymMatrix &m)
virtual G4double PostStepGetPhysicalInteractionLength(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
void BuildAsymmetryTable(const G4ParticleDefinition &part)
G4double ComputeAsymmetry(G4double energy, const G4MaterialCutsCouple *couple, const G4ParticleDefinition &particle, G4double cut, G4double &tAsymmetry)
G4double CrossSection(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4double currentInteractionLength
size_t CurrentMaterialCutsCoupleIndex() const
G4bool buildAsymmetryTable
G4double ComputeSaturationFactor(const G4Track &aTrack)
virtual G4bool IsApplicable(const G4ParticleDefinition &) override
virtual ~G4PolarizedCompton()
void SetEmModel(G4VEmModel *, G4int index=0)
G4double MinKinEnergy() const
G4int LambdaBinning() const
virtual G4double PostStepGetPhysicalInteractionLength(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
void SetSplineFlag(G4bool val)
void SetSecondaryParticle(const G4ParticleDefinition *p)
G4Material * GetMaterial() const
static G4Electron * Electron()
size_t GetTableSize() const
G4double MinKinEnergy() const
void SetMinKinEnergyPrim(G4double e)
virtual void PrintInfo() override
G4double MaxKinEnergy() const
void SetTargetPolarization(const G4ThreeVector &pTarget)
static G4ProductionCutsTable * GetProductionCutsTable()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
virtual G4double GetMeanFreePath(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
void SetBeamPolarization(const G4ThreeVector &pBeam)
const G4String & GetName() const
void SetLowEnergyLimit(G4double)
virtual void BuildPhysicsTable(const G4ParticleDefinition &) override
void SetProcessSubType(G4int)
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
const G4ThreeVector & GetVolumePolarization(G4LogicalVolume *lVol) const
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=nullptr)
virtual void BuildPhysicsTable(const G4ParticleDefinition &) override
G4bool GetFlag(size_t i) const
G4double theNumberOfInteractionLengthLeft
void SetStartFromNullFlag(G4bool val)
static G4PolarizationManager * GetInstance()
virtual void InitialiseProcess(const G4ParticleDefinition *) override
static G4PhysicsTable * theAsymmetryTable
virtual G4double GetMeanFreePath(const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition) override
const G4DynamicParticle * GetDynamicParticle() const
static G4EmParameters * Instance()
void PutValue(size_t index, G4double theValue)
const G4String & GetName() const
const G4String & GetName() const
G4VPhysicalVolume * GetVolume() const
bool IsPolarized(G4LogicalVolume *lVol) const