Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 /*/ 27 Authors: 28 M. Omer and R. Hajima on 15 November 2019 29 contact: 30 omer.mohamed@jaea.go.jp and hajima.ryoichi@qst.go.jp 31 Publication Information: 32 1- M. Omer, R. Hajima, Validating polarization effects in gamma-rays elastic scattering by Monte 33 Carlo simulation, New J. Phys., vol. 21, 2019, pp. 113006 (1-10), 34 https://doi.org/10.1088/1367-2630/ab4d8a 35 */ 36 37 #ifndef G4JAEAPolarizedElasticScatteringModel_h 38 #define G4JAEAPolarizedElasticScatteringModel_h 1 39 40 #include "G4VEmModel.hh" 41 #include "G4ParticleChangeForGamma.hh" 42 #include "G4PhysicsFreeVector.hh" 43 #include "G4ProductionCutsTable.hh" 44 #include "G4DataVector.hh" 45 46 class G4JAEAPolarizedElasticScatteringModel : public G4VEmModel 47 { 48 public: 49 explicit G4JAEAPolarizedElasticScatteringModel(); 50 virtual ~G4JAEAPolarizedElasticScatteringModel(); 51 52 void Initialise(const G4ParticleDefinition*, const G4DataVector&) override; 53 54 void InitialiseLocal(const G4ParticleDefinition*, 55 G4VEmModel* masterModel) override; 56 57 void InitialiseForElement(const G4ParticleDefinition*, G4int Z) override; 58 59 G4double ComputeCrossSectionPerAtom( 60 const G4ParticleDefinition*, 61 G4double kinEnergy, 62 G4double Z, 63 G4double A=0, 64 G4double cut=0, 65 G4double emax=DBL_MAX) override; 66 67 void SampleSecondaries(std::vector<G4DynamicParticle*>*, 68 const G4MaterialCutsCouple*, 69 const G4DynamicParticle*, 70 G4double tmin, 71 G4double maxEnergy) override; 72 73 inline void SetLowEnergyThreshold(G4double val){lowEnergyLimit = val;}; 74 inline void SetPolarizationSensitvity(G4bool,G4bool,G4bool); 75 76 void SetDebugVerbosity(G4int val){verboseLevel = val;}; 77 78 G4JAEAPolarizedElasticScatteringModel & operator=(const G4JAEAPolarizedElasticScatteringModel &right) = delete; 79 G4JAEAPolarizedElasticScatteringModel(const G4JAEAPolarizedElasticScatteringModel&) = delete; 80 81 private: 82 void ReadData(std::size_t Z, const char* path = 0); 83 G4double GeneratePolarizedPhi(G4double Sigma_para,G4double Sigma_perp, G4double initial_Pol_Plane); 84 85 static const G4int maxZ = 99; 86 static G4PhysicsFreeVector* dataCS[maxZ+1]; 87 static G4DataVector* Polarized_ES_Data[maxZ+1]; 88 G4double distribution[181]; 89 G4double cdistribution[181]; 90 91 G4ParticleChangeForGamma* fParticleChange; 92 G4double lowEnergyLimit; 93 94 G4int verboseLevel; 95 96 G4bool fLinearPolarizationSensitvity1; 97 G4bool fLinearPolarizationSensitvity2; 98 G4bool fCircularPolarizationSensitvity; 99 G4bool isInitialised; 100 }; 101 102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 103 104 inline void G4JAEAPolarizedElasticScatteringModel::SetPolarizationSensitvity(G4bool linear1, 105 G4bool linear2, 106 G4bool circular) 107 { 108 fLinearPolarizationSensitvity1=linear1; 109 fLinearPolarizationSensitvity2=linear2; 110 fCircularPolarizationSensitvity=circular; 111 } 112 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 114 115 #endif 116