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 5.0.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
                                                   >>  23 //
                                                   >>  24 // $Id: G4TauMinus.cc,v 1.8 2002/03/08 08:48:31 kurasige Exp $
                                                   >>  25 // GEANT4 tag $Name: geant4-05-00 $
                                                   >>  26 //
                                                   >>  27 // 
 26 // -------------------------------------------     28 // ----------------------------------------------------------------------
 27 //      GEANT 4 class implementation file          29 //      GEANT 4 class implementation file
 28 //                                                 30 //
 29 //      History: first implementation, based o     31 //      History: first implementation, based on object model of
 30 //      4th April 1996, G.Cosmo                    32 //      4th April 1996, G.Cosmo
 31 //      7 July 1996                   H.Kurash <<  33 //      7 July 1996                   H.Kurashige 
 32 // *******************************************     34 // **********************************************************************
 33 //  New impelemenataion as an utility class  M <<  35 
 34 // ------------------------------------------- <<  36 #include "g4std/fstream"
                                                   >>  37 #include "g4std/iomanip"
 35                                                    38 
 36 #include "G4TauMinus.hh"                           39 #include "G4TauMinus.hh"
 37                                                    40 
 38 #include "G4DecayTable.hh"                     << 
 39 #include "G4ParticleTable.hh"                  << 
 40 #include "G4PhaseSpaceDecayChannel.hh"             41 #include "G4PhaseSpaceDecayChannel.hh"
 41 #include "G4PhysicalConstants.hh"              << 
 42 #include "G4String.hh"                         << 
 43 #include "G4SystemOfUnits.hh"                  << 
 44 #include "G4TauLeptonicDecayChannel.hh"            42 #include "G4TauLeptonicDecayChannel.hh"
 45 #include "G4Types.hh"                          <<  43 #include "G4DecayTable.hh"
 46 #include "G4VDecayChannel.hh"                  << 
 47                                                << 
 48 G4TauMinus* G4TauMinus::theInstance = nullptr; << 
 49                                                    44 
 50 G4TauMinus* G4TauMinus::Definition()           <<  45 // ######################################################################
                                                   >>  46 // ###                          TAUMINUS                             ###
                                                   >>  47 // ######################################################################
                                                   >>  48 
                                                   >>  49 G4TauMinus::G4TauMinus(
                                                   >>  50        const G4String&     aName,        G4double            mass,
                                                   >>  51        G4double            width,        G4double            charge,   
                                                   >>  52        G4int               iSpin,        G4int               iParity,    
                                                   >>  53        G4int               iConjugation, G4int               iIsospin,   
                                                   >>  54        G4int               iIsospin3,    G4int               gParity,
                                                   >>  55        const G4String&     pType,        G4int               lepton,      
                                                   >>  56        G4int               baryon,       G4int               encoding,
                                                   >>  57        G4bool              stable,       G4double            lifetime,
                                                   >>  58        G4DecayTable        *decaytable )
                                                   >>  59  : G4VLepton( aName,mass,width,charge,iSpin,iParity,
                                                   >>  60               iConjugation,iIsospin,iIsospin3,gParity,pType,
                                                   >>  61               lepton,baryon,encoding,stable,lifetime,decaytable )
 51 {                                                  62 {
 52   if (theInstance != nullptr) return theInstan <<  63   SetParticleSubType("tau");
 53   const G4String name = "tau-";                <<  64   SetPDGStable(false);
 54   // search in particle table]                 <<  65 
 55   G4ParticleTable* pTable = G4ParticleTable::G <<  66   //create Decay Table 
 56   G4ParticleDefinition* anInstance = pTable->F <<  67   G4DecayTable*   table = GetDecayTable();
 57   if (anInstance == nullptr) {                 <<  68   if (table!=NULL) delete table;
 58     // create particle                         <<  69   table = new G4DecayTable();
 59     //                                         <<  70 
 60     //    Arguments for constructor are as fol <<  71   // create decay channels
 61     //               name             mass     <<  72   G4VDecayChannel** mode = new G4VDecayChannel*[6];
 62     //             2*spin           parity  C- <<  73   // tau- -> mu- + anti_nu_mu + nu_tau
 63     //          2*Isospin       2*Isospin3     <<  74   mode[0] = new G4TauLeptonicDecayChannel("tau-",0.174,"mu-");
 64     //               type    lepton number  ba <<  75   // tau- -> e- + anti_nu_e + nu_tau
 65     //             stable         lifetime     <<  76   mode[1] = new G4TauLeptonicDecayChannel("tau-",0.178,"e-");
 66     //             shortlived      subType     <<  77   // tau- -> pi- + nu_tau
 67                                                <<  78   mode[2] = new G4PhaseSpaceDecayChannel("tau-",0.113,2,"pi-","nu_tau");
 68     // clang-format off                        <<  79   // tau- -> pi0 + pi0 + pi- + nu_tau
 69    anInstance = new G4ParticleDefinition(      <<  80   mode[3] = new G4PhaseSpaceDecayChannel("tau-",0.252,3,"pi0","pi-","nu_tau");
 70                  name,     1.77686*GeV,  2.267 <<  81   // tau- -> pi0 + pi0 + pi- + nu_tau
                                                   >>  82   mode[4] = new G4PhaseSpaceDecayChannel();
                                                   >>  83   mode[4]->SetParent("tau-");
                                                   >>  84   mode[4]->SetBR(0.093);
                                                   >>  85   mode[4]->SetNumberOfDaughters(4);
                                                   >>  86   mode[4]->SetDaughter(0,"pi0");
                                                   >>  87   mode[4]->SetDaughter(1,"pi0");
                                                   >>  88   mode[4]->SetDaughter(2,"pi-");
                                                   >>  89   mode[4]->SetDaughter(3,"nu_tau");
                                                   >>  90   // tau- -> pi- + pi- + pi+ + nu_tau
                                                   >>  91   mode[5] = new G4PhaseSpaceDecayChannel();
                                                   >>  92   mode[5]->SetParent("tau-");
                                                   >>  93   mode[5]->SetBR(0.098);
                                                   >>  94   mode[5]->SetNumberOfDaughters(4);
                                                   >>  95   mode[5]->SetDaughter(0,"pi-");
                                                   >>  96   mode[5]->SetDaughter(1,"pi-");
                                                   >>  97   mode[5]->SetDaughter(2,"pi+");
                                                   >>  98   mode[5]->SetDaughter(3,"nu_tau");
                                                   >>  99 
                                                   >> 100   for (G4int index=0; index <6; index++ ) table->Insert(mode[index]);  
                                                   >> 101   delete [] mode;
                                                   >> 102   
                                                   >> 103   SetDecayTable(table);
                                                   >> 104 }
                                                   >> 105 
                                                   >> 106 // ......................................................................
                                                   >> 107 // ...                 static member definitions                      ...
                                                   >> 108 // ......................................................................
                                                   >> 109 //     
                                                   >> 110 //    Arguments for constructor are as follows
                                                   >> 111 //               name             mass          width         charge
                                                   >> 112 //             2*spin           parity  C-conjugation
                                                   >> 113 //          2*Isospin       2*Isospin3       G-parity
                                                   >> 114 //               type    lepton number  baryon number   PDG encoding
                                                   >> 115 //             stable         lifetime    decay table 
                                                   >> 116 
                                                   >> 117 G4TauMinus G4TauMinus::theTauMinus(
                                                   >> 118     "tau-",    1.77705*GeV,   2.27e-9*MeV,    -1.*eplus, 
 71         1,               0,             0,        119         1,               0,             0,          
 72         0,               0,             0,        120         0,               0,             0,             
 73        "lepton",               1,                 121        "lepton",               1,             0,          15,
 74     false,     290.3e-6*ns,          nullptr,  << 122      true,     295.6e-6*ns,          NULL
 75              false,           "tau"            << 123 );
 76               );                               << 
 77     // clang-format on                         << 
 78                                                << 
 79     // Bohr Magnetron                          << 
 80     G4double muB = -0.5 * eplus * hbar_Planck  << 
 81                                                << 
 82     anInstance->SetPDGMagneticMoment(muB * 1.0 << 
 83                                                << 
 84     // create Decay 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   }                                            << 
124   theInstance = static_cast<G4TauMinus*>(anIns << 
125   return theInstance;                          << 
126 }                                              << 
127                                                   124 
128 G4TauMinus* G4TauMinus::TauMinusDefinition()   << 125 G4TauMinus* G4TauMinus::TauMinusDefinition() {return &theTauMinus;}
129 {                                              << 
130   return Definition();                         << 
131 }                                              << 
132                                                   126 
133 G4TauMinus* G4TauMinus::TauMinus()             << 127 G4TauMinus*  G4TauMinus::TauMinus()
134 {                                              << 128 { 
135   return Definition();                         << 129   return &theTauMinus; 
136 }                                                 130 }
                                                   >> 131 
                                                   >> 132 
137                                                   133