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 // Geant4 Class file 27 // 28 // File name: G4PolarizedComptonXS 29 // 30 // Author: Andreas Schaelicke 31 // 32 // Class Description: 33 // determine the polarization of the final 34 // in a Compton scattering process employing 35 // cross section by F.W.Lipps & H.A.Tolhoek 36 // ( Physica 20 (1954) 395 ) 37 38 #ifndef G4PolarizedComptonXS_h 39 #define G4PolarizedComptonXS_h 1 40 41 #include "globals.hh" 42 #include "G4StokesVector.hh" 43 #include "G4VPolarizedXS.hh" 44 45 class G4PolarizedComptonXS : public G4VPolariz 46 { 47 public: 48 G4PolarizedComptonXS(); 49 ~G4PolarizedComptonXS() override; 50 51 // prepares the ingredients for the calculat 52 // dependent differential cross section 53 // the kinematics is fixed (X - incoming pho 54 // mass, eps - outgoing photon energy in uni 55 // and polarization of the incoming particle 56 // a flag specifies the extent to which pola 57 void Initialize(G4double eps, G4double X, G4 58 const G4StokesVector& p0, co 59 G4int flag = 0) override; 60 61 // returns the differential cross section fo 62 // of the final state particles to be used i 63 // polarization transfer the calculation has 64 // Initialize() prior to the first call of t 65 G4double XSection(const G4StokesVector& pol2 66 const G4StokesVector& pol3 67 // total cross section 68 G4double TotalXSection(G4double xmin, G4doub 69 const G4StokesVector& 70 const G4StokesVector& 71 72 // return expected mean polarisation 73 G4StokesVector GetPol2() override; 74 G4StokesVector GetPol3() override; 75 76 G4PolarizedComptonXS& operator=(const G4Pola 77 G4PolarizedComptonXS(const G4PolarizedCompto 78 79 private: 80 void DefineCoefficients(const G4StokesVector 81 const G4StokesVector 82 83 static constexpr G4double re2 = 84 CLHEP::classic_electr_radius * CLHEP::clas 85 (4. * CLHEP::pi / CLHEP::hbarc) * (4. * CL 86 87 // these variables store the information nec 88 // differential cross section for arbitrary 89 // polarizations (used in XSection): 90 // - part depending on the polarization of t 91 G4ThreeVector fPhi2; 92 // - part depending on the polarization of t 93 G4ThreeVector fPhi3; 94 // - polarization independent part 95 G4double fPhi0; 96 // - product of polarizations of initial par 97 G4double polxx, polyy, polzz, polxz, polzx, 98 99 // G4double diffXSFactor, totalXSFactor; 100 G4double fPolXS, fUnpXS; 101 }; 102 103 #endif 104