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 // G4MuonDecayChannelWithSpin 27 // 28 // Class decription: 29 // 30 // This class describes muon decay kinematics. 31 // It assumes V-A coupling with 1st order radi 32 // standard model parameter values, but gives 33 // for neutrinos. 34 // References: 35 // - Florian Scheck "Muon Physics", in Physics 36 // (Review Section of Physics Letters) 44, N 37 // 187-248. North-Holland Publishing Company 38 // - W.E. Fisher and F. Scheck, Nucl. Phys. B8 39 40 // Authors: P.Gumplinger and T.MacPhail, 17 Au 41 // ------------------------------------------- 42 #ifndef G4MuonDecayChannelWithSpin_hh 43 #define G4MuonDecayChannelWithSpin_hh 1 44 45 #include "G4MuonDecayChannel.hh" 46 #include "G4ThreeVector.hh" 47 #include "globals.hh" 48 49 #include <CLHEP/Units/PhysicalConstants.h> 50 51 class G4MuonDecayChannelWithSpin : public G4Mu 52 { 53 public: 54 G4MuonDecayChannelWithSpin(const G4String& 55 ~G4MuonDecayChannelWithSpin() override = d 56 57 G4DecayProducts* DecayIt(G4double) overrid 58 59 protected: 60 // Copy constructor and assignment operato 61 G4MuonDecayChannelWithSpin(const G4MuonDec 62 G4MuonDecayChannelWithSpin& operator=(cons 63 64 private: 65 G4MuonDecayChannelWithSpin() = default; 66 67 // Radiative Correction Factors 68 inline G4double F_c(G4double x, G4double x 69 inline G4double F_theta(G4double x, G4doub 70 G4double R_c(G4double x, G4double omega); 71 }; 72 73 // ------------------------ 74 // Inline methods 75 // ------------------------ 76 77 inline G4double G4MuonDecayChannelWithSpin::F_ 78 { 79 G4double f_c; 80 81 f_c = (5. + 17. * x - 34. * x * x) * (omega 82 f_c = (1. - x) / (3. * x * x) * f_c; 83 f_c = (6. - 4. * x) * R_c(x, omega) + (6. - 84 f_c = (CLHEP::fine_structure_const / CLHEP:: 85 86 return f_c; 87 } 88 89 inline G4double G4MuonDecayChannelWithSpin::F_ 90 { 91 G4double f_theta; 92 93 f_theta = (1. + x + 34 * x * x) * (omega + s 94 f_theta = f_theta + ((4. * (1. - x) * (1. - 95 f_theta = (1. - x) / (3. * x * x) * f_theta; 96 f_theta = (2. - 4. * x) * R_c(x, omega) + (2 97 f_theta = (CLHEP::fine_structure_const / CLH 98 99 return f_theta; 100 } 101 102 #endif 103