Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // ------------------------------------------- 27 // 28 // Geant4 class file 29 // 30 // File name: G4PolarizedBremsstrahlungMod 31 // 32 // Author: Karim Laihem 33 34 #include "G4PolarizedBremsstrahlungModel.hh" 35 36 #include "G4ParticleChangeForLoss.hh" 37 #include "G4PolarizationHelper.hh" 38 #include "G4PolarizedBremsstrahlungXS.hh" 39 40 G4PolarizedBremsstrahlungModel::G4PolarizedBre 41 const G4ParticleDefinition* p, const G4Strin 42 : G4SeltzerBergerModel(p, nam) 43 , fCrossSectionCalculator(nullptr) 44 {} 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 G4PolarizedBremsstrahlungModel::~G4PolarizedBr 48 { 49 if(fCrossSectionCalculator) 50 delete fCrossSectionCalculator; 51 } 52 53 //....oooOO0OOooo........oooOO0OOooo........oo 54 void G4PolarizedBremsstrahlungModel::Initialis 55 56 { 57 G4SeltzerBergerModel::Initialise(p, d); 58 if(!fCrossSectionCalculator) 59 fCrossSectionCalculator = new G4PolarizedB 60 } 61 62 //....oooOO0OOooo........oooOO0OOooo........oo 63 // The emitted gamma energy is sampled using a 64 void G4PolarizedBremsstrahlungModel::SampleSec 65 std::vector<G4DynamicParticle*>* vdp, const 66 const G4DynamicParticle* dp, G4double tmin, 67 { 68 G4SeltzerBergerModel::SampleSecondaries(vdp, 69 std::size_t num = vdp->size(); 70 71 if(num > 0) 72 { 73 G4double lepEnergy0 = dp->GetKineticEnergy 74 G4double gamEnergy1 = (*vdp)[0]->GetKineti 75 G4double sintheta = 76 dp->GetMomentumDirection().cross((*vdp)[ 77 if(sintheta > 1.) 78 sintheta = 1.; 79 80 G4StokesVector beamPol = G4StokesVector(dp 81 82 // determine interaction plane 83 G4ThreeVector nInteractionFrame = G4Polari 84 dp->GetMomentumDirection(), 85 fParticleChange->GetProposedMomentumDire 86 87 // transform polarization into interaction 88 beamPol.InvRotateAz(nInteractionFrame, dp- 89 90 // calculate polarization transfer 91 fCrossSectionCalculator->SetMaterial( 92 GetCurrentElement()->GetN(), // number 93 GetCurrentElement()->GetZ(), GetCurrentE 94 fCrossSectionCalculator->Initialize(lepEne 95 beamPo 96 97 // determine final state polarization 98 G4StokesVector newBeamPol = fCrossSectionC 99 newBeamPol.RotateAz(nInteractionFrame, 100 fParticleChange->GetPr 101 fParticleChange->ProposePolarization(newBe 102 103 if(num != 1) 104 { 105 G4ExceptionDescription ed; 106 ed << num << " secondaries in polarized 107 G4Exception("G4PolarizedBremsstrahlungMo 108 JustWarning, ed); 109 } 110 for(std::size_t i = 0; i < num; ++i) 111 { 112 G4StokesVector photonPol = fCrossSection 113 photonPol.SetPhoton(); 114 photonPol.RotateAz(nInteractionFrame, (* 115 (*vdp)[i]->SetPolarization(photonPol.p1( 116 photonPol.p3( 117 } 118 } 119 return; 120 } 121