Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/particles/management/include/G4MuonDecayChannelWithSpin.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 // G4MuonDecayChannelWithSpin
 27 //
 28 // Class decription:
 29 //
 30 // This class describes muon decay kinematics.
 31 // It assumes V-A coupling with 1st order radiative corrections, the
 32 // standard model parameter values, but gives incorrect energy spectrum
 33 // for neutrinos.
 34 // References:
 35 // - Florian Scheck "Muon Physics", in Physics Reports
 36 //   (Review Section of Physics Letters) 44, No. 4 (1978)
 37 //   187-248. North-Holland Publishing Company, Amsterdam at page 210 cc.
 38 // - W.E. Fisher and F. Scheck, Nucl. Phys. B83 (1974) 25.
 39 
 40 // Authors: P.Gumplinger and T.MacPhail, 17 August 2004
 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 G4MuonDecayChannel
 52 {
 53   public:
 54     G4MuonDecayChannelWithSpin(const G4String& theParentName, G4double theBR);
 55     ~G4MuonDecayChannelWithSpin() override = default;
 56 
 57     G4DecayProducts* DecayIt(G4double) override;
 58 
 59   protected:
 60     // Copy constructor and assignment operator
 61     G4MuonDecayChannelWithSpin(const G4MuonDecayChannelWithSpin&) = default;
 62     G4MuonDecayChannelWithSpin& operator=(const G4MuonDecayChannelWithSpin&);
 63 
 64   private:
 65     G4MuonDecayChannelWithSpin() = default;
 66 
 67     // Radiative Correction Factors
 68     inline G4double F_c(G4double x, G4double x0, G4double omega);
 69     inline G4double F_theta(G4double x, G4double x0, G4double omega);
 70     G4double R_c(G4double x, G4double omega);
 71 };
 72 
 73 // ------------------------
 74 // Inline methods
 75 // ------------------------
 76 
 77 inline G4double G4MuonDecayChannelWithSpin::F_c(G4double x, G4double x0, G4double omega)
 78 {
 79   G4double f_c;
 80 
 81   f_c = (5. + 17. * x - 34. * x * x) * (omega + std::log(x)) - 22. * x + 34. * x * x;
 82   f_c = (1. - x) / (3. * x * x) * f_c;
 83   f_c = (6. - 4. * x) * R_c(x, omega) + (6. - 6. * x) * std::log(x) + f_c;
 84   f_c = (CLHEP::fine_structure_const / CLHEP::twopi) * (x * x - x0 * x0) * f_c;
 85 
 86   return f_c;
 87 }
 88 
 89 inline G4double G4MuonDecayChannelWithSpin::F_theta(G4double x, G4double x0, G4double omega)
 90 {
 91   G4double f_theta;
 92 
 93   f_theta = (1. + x + 34 * x * x) * (omega + std::log(x)) + 3. - 7. * x - 32. * x * x;
 94   f_theta = f_theta + ((4. * (1. - x) * (1. - x)) / x) * std::log(1. - x);
 95   f_theta = (1. - x) / (3. * x * x) * f_theta;
 96   f_theta = (2. - 4. * x) * R_c(x, omega) + (2. - 6. * x) * std::log(x) - f_theta;
 97   f_theta = (CLHEP::fine_structure_const / CLHEP::twopi) * (x * x - x0 * x0) * f_theta;
 98 
 99   return f_theta;
100 }
101 
102 #endif
103