Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/particles/leptons/src/G4TauMinus.cc

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/leptons/src/G4TauMinus.cc (Version 11.3.0) and /particles/leptons/src/G4TauMinus.cc (Version 9.6.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 //
                                                   >>  27 // $Id$
                                                   >>  28 //
                                                   >>  29 // 
 26 // -------------------------------------------     30 // ----------------------------------------------------------------------
 27 //      GEANT 4 class implementation file          31 //      GEANT 4 class implementation file
 28 //                                                 32 //
 29 //      History: first implementation, based o     33 //      History: first implementation, based on object model of
 30 //      4th April 1996, G.Cosmo                    34 //      4th April 1996, G.Cosmo
 31 //      7 July 1996                   H.Kurash <<  35 //      7 July 1996                   H.Kurashige 
 32 // *******************************************     36 // **********************************************************************
 33 //  New impelemenataion as an utility class  M     37 //  New impelemenataion as an utility class  M.Asai, 26 July 2004
 34 // -------------------------------------------     38 // ----------------------------------------------------------------------
 35                                                    39 
 36 #include "G4TauMinus.hh"                           40 #include "G4TauMinus.hh"
 37                                                << 
 38 #include "G4DecayTable.hh"                     << 
 39 #include "G4ParticleTable.hh"                  << 
 40 #include "G4PhaseSpaceDecayChannel.hh"         << 
 41 #include "G4PhysicalConstants.hh"                  41 #include "G4PhysicalConstants.hh"
 42 #include "G4String.hh"                         << 
 43 #include "G4SystemOfUnits.hh"                      42 #include "G4SystemOfUnits.hh"
                                                   >>  43 #include "G4ParticleTable.hh"
                                                   >>  44 
                                                   >>  45 #include "G4PhaseSpaceDecayChannel.hh"
 44 #include "G4TauLeptonicDecayChannel.hh"            46 #include "G4TauLeptonicDecayChannel.hh"
 45 #include "G4Types.hh"                          <<  47 #include "G4DecayTable.hh"
 46 #include "G4VDecayChannel.hh"                  << 
 47                                                    48 
 48 G4TauMinus* G4TauMinus::theInstance = nullptr; <<  49 // ######################################################################
                                                   >>  50 // ###                          TAUMINUS                             ###
                                                   >>  51 // ######################################################################
                                                   >>  52 G4TauMinus* G4TauMinus::theInstance = 0;
 49                                                    53 
 50 G4TauMinus* G4TauMinus::Definition()               54 G4TauMinus* G4TauMinus::Definition()
 51 {                                                  55 {
 52   if (theInstance != nullptr) return theInstan <<  56   if (theInstance !=0) return theInstance;
 53   const G4String name = "tau-";                    57   const G4String name = "tau-";
 54   // search in particle table]                     58   // search in particle table]
 55   G4ParticleTable* pTable = G4ParticleTable::G     59   G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
 56   G4ParticleDefinition* anInstance = pTable->F     60   G4ParticleDefinition* anInstance = pTable->FindParticle(name);
 57   if (anInstance == nullptr) {                 <<  61   if (anInstance ==0)
 58     // create particle                         <<  62   {
 59     //                                         <<  63   // create particle
 60     //    Arguments for constructor are as fol <<  64   //
 61     //               name             mass     <<  65   //    Arguments for constructor are as follows
 62     //             2*spin           parity  C- <<  66   //               name             mass          width         charge
 63     //          2*Isospin       2*Isospin3     <<  67   //             2*spin           parity  C-conjugation
 64     //               type    lepton number  ba <<  68   //          2*Isospin       2*Isospin3       G-parity
 65     //             stable         lifetime     <<  69   //               type    lepton number  baryon number   PDG encoding
 66     //             shortlived      subType     <<  70   //             stable         lifetime    decay table
 67                                                <<  71   //             shortlived      subType    anti_encoding
 68     // clang-format off                        << 
 69    anInstance = new G4ParticleDefinition(          72    anInstance = new G4ParticleDefinition(
 70                  name,     1.77686*GeV,  2.267 <<  73                  name,     1.77682*GeV,  2.265e-9*MeV,    -1.*eplus, 
 71         1,               0,             0,         74         1,               0,             0,          
 72         0,               0,             0,         75         0,               0,             0,             
 73        "lepton",               1,                  76        "lepton",               1,             0,          15,
 74     false,     290.3e-6*ns,          nullptr,  <<  77     false,     290.6e-6*ns,          NULL,
 75              false,           "tau"                78              false,           "tau"
 76               );                                   79               );
 77     // clang-format on                         << 
 78                                                << 
 79     // Bohr Magnetron                          << 
 80     G4double muB = -0.5 * eplus * hbar_Planck  << 
 81                                                    80 
 82     anInstance->SetPDGMagneticMoment(muB * 1.0 <<  81    // Bohr Magnetron
                                                   >>  82    G4double muB =  -0.5*eplus*hbar_Planck/(1.77699*GeV/c_squared) ;
                                                   >>  83    
                                                   >>  84    anInstance->SetPDGMagneticMoment( muB * 2. * 1.00116);
                                                   >>  85 
                                                   >>  86   //create Decay Table 
                                                   >>  87   G4DecayTable* table = new G4DecayTable();
                                                   >>  88 
                                                   >>  89   // create decay channels
                                                   >>  90   G4VDecayChannel* mode;
                                                   >>  91   // tau- -> mu- + anti_nu_mu + nu_tau
                                                   >>  92   mode = new G4TauLeptonicDecayChannel("tau-",0.1736,"mu-");
                                                   >>  93   table->Insert(mode);
                                                   >>  94   // tau- -> e- + anti_nu_e + nu_tau
                                                   >>  95   mode = new G4TauLeptonicDecayChannel("tau-",0.1784,"e-");
                                                   >>  96   table->Insert(mode);
                                                   >>  97   // tau- -> pi- + nu_tau
                                                   >>  98   mode = new G4PhaseSpaceDecayChannel("tau-",0.1106,2,"pi-","nu_tau");
                                                   >>  99   table->Insert(mode);
                                                   >> 100   // tau- -> pi0 + pi- + nu_tau
                                                   >> 101   mode = new G4PhaseSpaceDecayChannel("tau-",0.2541,3,"pi0","pi-","nu_tau");
                                                   >> 102   table->Insert(mode);
                                                   >> 103   // tau- -> pi0 + pi0 + pi- + nu_tau
                                                   >> 104   mode = new G4PhaseSpaceDecayChannel();
                                                   >> 105   mode->SetParent("tau-");
                                                   >> 106   mode->SetBR(0.0917);
                                                   >> 107   mode->SetNumberOfDaughters(4);
                                                   >> 108   mode->SetDaughter(0,"pi0");
                                                   >> 109   mode->SetDaughter(1,"pi0");
                                                   >> 110   mode->SetDaughter(2,"pi-");
                                                   >> 111   mode->SetDaughter(3,"nu_tau");
                                                   >> 112   table->Insert(mode);
                                                   >> 113   // tau- -> pi- + pi- + pi+ + nu_tau
                                                   >> 114   mode = new G4PhaseSpaceDecayChannel();
                                                   >> 115   mode->SetParent("tau-");
                                                   >> 116   mode->SetBR(0.0946);
                                                   >> 117   mode->SetNumberOfDaughters(4);
                                                   >> 118   mode->SetDaughter(0,"pi-");
                                                   >> 119   mode->SetDaughter(1,"pi-");
                                                   >> 120   mode->SetDaughter(2,"pi+");
                                                   >> 121   mode->SetDaughter(3,"nu_tau");
                                                   >> 122   table->Insert(mode);
 83                                                   123 
 84     // create Decay Table                      << 124   anInstance->SetDecayTable(table);
 85     auto table = new G4DecayTable();           << 
 86                                                << 
 87     // create decay channels                   << 
 88     G4VDecayChannel* mode;                     << 
 89     // tau- -> mu- + anti_nu_mu + nu_tau       << 
 90     mode = new G4TauLeptonicDecayChannel("tau- << 
 91     table->Insert(mode);                       << 
 92     // tau- -> e- + anti_nu_e + nu_tau         << 
 93     mode = new G4TauLeptonicDecayChannel("tau- << 
 94     table->Insert(mode);                       << 
 95     // tau- -> pi- + nu_tau                    << 
 96     mode = new G4PhaseSpaceDecayChannel("tau-" << 
 97     table->Insert(mode);                       << 
 98     // tau- -> pi0 + pi- + nu_tau              << 
 99     mode = new G4PhaseSpaceDecayChannel("tau-" << 
100     table->Insert(mode);                       << 
101     // tau- -> pi0 + pi0 + pi- + nu_tau        << 
102     mode = new G4PhaseSpaceDecayChannel();     << 
103     mode->SetParent("tau-");                   << 
104     mode->SetBR(0.0917);                       << 
105     mode->SetNumberOfDaughters(4);             << 
106     mode->SetDaughter(0, "pi0");               << 
107     mode->SetDaughter(1, "pi0");               << 
108     mode->SetDaughter(2, "pi-");               << 
109     mode->SetDaughter(3, "nu_tau");            << 
110     table->Insert(mode);                       << 
111     // tau- -> pi- + pi- + pi+ + nu_tau        << 
112     mode = new G4PhaseSpaceDecayChannel();     << 
113     mode->SetParent("tau-");                   << 
114     mode->SetBR(0.0946);                       << 
115     mode->SetNumberOfDaughters(4);             << 
116     mode->SetDaughter(0, "pi-");               << 
117     mode->SetDaughter(1, "pi-");               << 
118     mode->SetDaughter(2, "pi+");               << 
119     mode->SetDaughter(3, "nu_tau");            << 
120     table->Insert(mode);                       << 
121                                                << 
122     anInstance->SetDecayTable(table);          << 
123   }                                               125   }
124   theInstance = static_cast<G4TauMinus*>(anIns << 126   theInstance = reinterpret_cast<G4TauMinus*>(anInstance);
125   return theInstance;                             127   return theInstance;
126 }                                                 128 }
127                                                   129 
128 G4TauMinus* G4TauMinus::TauMinusDefinition()   << 130 G4TauMinus*  G4TauMinus::TauMinusDefinition()
129 {                                                 131 {
130   return Definition();                            132   return Definition();
131 }                                                 133 }
132                                                   134 
133 G4TauMinus* G4TauMinus::TauMinus()             << 135 G4TauMinus*  G4TauMinus::TauMinus()
134 {                                                 136 {
135   return Definition();                            137   return Definition();
136 }                                                 138 }
                                                   >> 139 
                                                   >> 140 
137                                                   141