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 ]

Diff markup

Differences between /particles/management/include/G4MuonDecayChannelWithSpin.hh (Version 11.3.0) and /particles/management/include/G4MuonDecayChannelWithSpin.hh (Version 10.2.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 // G4MuonDecayChannelWithSpin                  <<  26 // ------------------------------------------------------------
                                                   >>  27 //      GEANT 4 class header file
 27 //                                                 28 //
 28 // Class decription:                           <<  29 //      History:
                                                   >>  30 //               17 August 2004 P.Gumplinger and T.MacPhail
                                                   >>  31 //               samples Michel spectrum including 1st order
                                                   >>  32 //               radiative corrections
                                                   >>  33 //               Reference: Florian Scheck "Muon Physics", in Physics Reports
                                                   >>  34 //                          (Review Section of Physics Letters) 44, No. 4 (1978)
                                                   >>  35 //                          187-248. North-Holland Publishing Company, Amsterdam
                                                   >>  36 //                          at page 210 cc.
 29 //                                                 37 //
 30 // This class describes muon decay kinematics. <<  38 //                          W.E. Fisher and F. Scheck, Nucl. Phys. B83 (1974) 25.
 31 // It assumes V-A coupling with 1st order radi <<  39 //
 32 // standard model parameter values, but gives  <<  40 // ------------------------------------------------------------
 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              41 #ifndef G4MuonDecayChannelWithSpin_hh
 43 #define G4MuonDecayChannelWithSpin_hh 1            42 #define G4MuonDecayChannelWithSpin_hh 1
 44                                                    43 
 45 #include "G4MuonDecayChannel.hh"               << 
 46 #include "G4ThreeVector.hh"                    << 
 47 #include "globals.hh"                          << 
 48                                                << 
 49 #include <CLHEP/Units/PhysicalConstants.h>         44 #include <CLHEP/Units/PhysicalConstants.h>
 50                                                    45 
                                                   >>  46 #include "globals.hh"
                                                   >>  47 #include "G4ThreeVector.hh"
                                                   >>  48 #include "G4MuonDecayChannel.hh"
                                                   >>  49 
 51 class G4MuonDecayChannelWithSpin : public G4Mu     50 class G4MuonDecayChannelWithSpin : public G4MuonDecayChannel
 52 {                                                  51 {
 53   public:                                      <<  52   // Class Decription
 54     G4MuonDecayChannelWithSpin(const G4String& <<  53   // This class describes muon decay kinemtics.
 55     ~G4MuonDecayChannelWithSpin() override = d <<  54   // This version assumes V-A coupling with 1st order radiative correctons,
 56                                                <<  55   //              the standard model Michel parameter values, but 
 57     G4DecayProducts* DecayIt(G4double) overrid <<  56   //              gives incorrect energy spectrum for neutrinos
 58                                                <<  57 
 59   protected:                                   <<  58 public:  // With Description
 60     // Copy constructor and assignment operato <<  59 
 61     G4MuonDecayChannelWithSpin(const G4MuonDec <<  60   //Constructors 
 62     G4MuonDecayChannelWithSpin& operator=(cons <<  61   G4MuonDecayChannelWithSpin(const G4String& theParentName,
 63                                                <<  62              G4double        theBR);
 64   private:                                     <<  63   //  Destructor
 65     G4MuonDecayChannelWithSpin() = default;    <<  64   virtual ~G4MuonDecayChannelWithSpin();
 66                                                <<  65 
 67     // Radiative Correction Factors            <<  66 protected:
 68     inline G4double F_c(G4double x, G4double x <<  67   // Copy constructor and assignment operator
 69     inline G4double F_theta(G4double x, G4doub <<  68   G4MuonDecayChannelWithSpin(const G4MuonDecayChannelWithSpin &);
 70     G4double R_c(G4double x, G4double omega);  <<  69   G4MuonDecayChannelWithSpin & operator=(const G4MuonDecayChannelWithSpin &);
 71 };                                             <<  70   
                                                   >>  71 private:
                                                   >>  72   G4MuonDecayChannelWithSpin();
                                                   >>  73 
                                                   >>  74 public:  // With Description
                                                   >>  75 
                                                   >>  76   virtual G4DecayProducts *DecayIt(G4double);
                                                   >>  77 
                                                   >>  78 private:
                                                   >>  79 // Radiative Correction Factors
                                                   >>  80 
                                                   >>  81   G4double F_c(G4double x, G4double x0);
                                                   >>  82   G4double F_theta(G4double x, G4double x0);
                                                   >>  83   G4double R_c(G4double x);
 72                                                    84 
 73 // ------------------------                    <<  85   G4double EMMU;
 74 // Inline methods                              <<  86   G4double EMASS;
 75 // ------------------------                    << 
 76                                                    87 
 77 inline G4double G4MuonDecayChannelWithSpin::F_ <<  88 };
                                                   >>  89 
                                                   >>  90 inline G4double G4MuonDecayChannelWithSpin::F_c(G4double x, G4double x0)
 78 {                                                  91 {
                                                   >>  92   G4double omega = std::log(EMMU/EMASS);
                                                   >>  93 
 79   G4double f_c;                                    94   G4double f_c;
 80                                                    95 
 81   f_c = (5. + 17. * x - 34. * x * x) * (omega  <<  96   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;         <<  97   f_c = (1.-x)/(3.*x*x)*f_c;
 83   f_c = (6. - 4. * x) * R_c(x, omega) + (6. -  <<  98   f_c = (6.-4.*x)*R_c(x)+(6.-6.*x)*std::log(x) + f_c;
 84   f_c = (CLHEP::fine_structure_const / CLHEP:: <<  99   f_c = (CLHEP::fine_structure_const/CLHEP::twopi) * (x*x-x0*x0) * f_c;
 85                                                   100 
 86   return f_c;                                     101   return f_c;
 87 }                                                 102 }
 88                                                   103 
 89 inline G4double G4MuonDecayChannelWithSpin::F_ << 104 inline G4double G4MuonDecayChannelWithSpin::F_theta(G4double x, G4double x0)
 90 {                                                 105 {
                                                   >> 106   G4double omega = std::log(EMMU/EMASS);
                                                   >> 107 
 91   G4double f_theta;                               108   G4double f_theta;
 92                                                   109 
 93   f_theta = (1. + x + 34 * x * x) * (omega + s << 110   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. -  << 111   f_theta = f_theta + ((4.*(1.-x)*(1.-x))/x)*std::log(1.-x);
 95   f_theta = (1. - x) / (3. * x * x) * f_theta; << 112   f_theta = (1.-x)/(3.*x*x) * f_theta;
 96   f_theta = (2. - 4. * x) * R_c(x, omega) + (2 << 113   f_theta = (2.-4.*x)*R_c(x)+(2.-6.*x)*std::log(x)-f_theta;
 97   f_theta = (CLHEP::fine_structure_const / CLH << 114   f_theta = (CLHEP::fine_structure_const/CLHEP::twopi) * (x*x-x0*x0) * f_theta;
 98                                                   115 
 99   return f_theta;                                 116   return f_theta;
100 }                                                 117 }
101                                                   118 
102 #endif                                            119 #endif
103                                                   120