Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/particles/leptons/src/G4TauPlus.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/G4TauPlus.cc (Version 11.3.0) and /particles/leptons/src/G4TauPlus.cc (Version 4.0)


  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: G4TauPlus.cc,v 1.7 2001/10/16 08:16:19 kurasige Exp $
                                                   >>  25 // GEANT4 tag $Name: geant4-04-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 "G4TauPlus.hh"                            39 #include "G4TauPlus.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 "G4DecayTable.hh"
 42 #include "G4String.hh"                         << 
 43 #include "G4SystemOfUnits.hh"                  << 
 44 #include "G4TauLeptonicDecayChannel.hh"        << 
 45 #include "G4Types.hh"                          << 
 46 #include "G4VDecayChannel.hh"                  << 
 47                                                << 
 48 G4TauPlus* G4TauPlus::theInstance = nullptr;   << 
 49                                                    43 
 50 G4TauPlus* G4TauPlus::Definition()             <<  44 // ######################################################################
                                                   >>  45 // ###                          TAPLUS                                ###
                                                   >>  46 // ######################################################################
                                                   >>  47 
                                                   >>  48 G4TauPlus::G4TauPlus(
                                                   >>  49        const G4String&     aName,        G4double            mass,
                                                   >>  50        G4double            width,        G4double            charge,   
                                                   >>  51        G4int               iSpin,        G4int               iParity,    
                                                   >>  52        G4int               iConjugation, G4int               iIsospin,   
                                                   >>  53        G4int               iIsospin3,    G4int               gParity,
                                                   >>  54        const G4String&     pType,        G4int               lepton,      
                                                   >>  55        G4int               baryon,       G4int               encoding,
                                                   >>  56        G4bool              stable,       G4double            lifetime,
                                                   >>  57        G4DecayTable        *decaytable )
                                                   >>  58  : G4VLepton( aName,mass,width,charge,iSpin,iParity,
                                                   >>  59               iConjugation,iIsospin,iIsospin3,gParity,pType,
                                                   >>  60               lepton,baryon,encoding,stable,lifetime,decaytable )
 51 {                                                  61 {
 52   if (theInstance != nullptr) return theInstan <<  62   SetParticleSubType("tau");
 53   const G4String name = "tau+";                <<  63   SetPDGStable(false);
 54   // search in particle table]                 <<  64 
 55   G4ParticleTable* pTable = G4ParticleTable::G <<  65   //create Decay Table 
 56   G4ParticleDefinition* anInstance = pTable->F <<  66   G4DecayTable*   table = GetDecayTable();
 57   if (anInstance == nullptr) {                 <<  67   if (table!=NULL) delete table;
 58     // create particle                         <<  68   table = new G4DecayTable();
 59     //                                         <<  69 
 60     //    Arguments for constructor are as fol <<  70   // create decay channels
 61     //               name             mass     <<  71   G4VDecayChannel** mode = new G4VDecayChannel*[6];
 62     //             2*spin           parity  C- <<  72   // tau+ -> mu+ + nu_mu + anti_nu_tau
 63     //          2*Isospin       2*Isospin3     <<  73   mode[0] = new G4PhaseSpaceDecayChannel("tau+",0.174,3,"mu+","nu_mu","anti_nu_tau");
 64     //               type    lepton number  ba <<  74   // tau+ -> e+ + nu_e + anti_nu_tau
 65     //             stable         lifetime     <<  75   mode[1] = new G4PhaseSpaceDecayChannel("tau+",0.178,3,"e+","nu_e","anti_nu_tau");
 66     //             shortlived      subType     <<  76   // tau+ -> pi+ + anti_nu_tau
 67                                                <<  77   mode[2] = new G4PhaseSpaceDecayChannel("tau+",0.113,2,"pi+","anti_nu_tau");
 68     // clang-format off                        <<  78   // tau+ -> pi0 + pi0 + pi+ + anti_nu_tau
 69    anInstance = new G4ParticleDefinition(      <<  79   mode[3] = new G4PhaseSpaceDecayChannel("tau+",0.252,3,"pi0","pi+","anti_nu_tau");
 70                  name,     1.77686*GeV,  2.267 <<  80   // tau+ -> pi0 + pi0 + pi+ + anti_nu_tau
                                                   >>  81   mode[4] = new G4PhaseSpaceDecayChannel();
                                                   >>  82   mode[4]->SetParent("tau+");
                                                   >>  83   mode[4]->SetBR(0.093);
                                                   >>  84   mode[4]->SetNumberOfDaughters(4);
                                                   >>  85   mode[4]->SetDaughter(0,"pi0");
                                                   >>  86   mode[4]->SetDaughter(1,"pi0");
                                                   >>  87   mode[4]->SetDaughter(2,"pi+");
                                                   >>  88   mode[4]->SetDaughter(3,"anti_nu_tau");
                                                   >>  89   // tau+ -> pi+ + pi+ + pi- + anti_nu_tau
                                                   >>  90   mode[5] = new G4PhaseSpaceDecayChannel();
                                                   >>  91   mode[5]->SetParent("tau+");
                                                   >>  92   mode[5]->SetBR(0.098);
                                                   >>  93   mode[5]->SetNumberOfDaughters(4);
                                                   >>  94   mode[5]->SetDaughter(0,"pi+");
                                                   >>  95   mode[5]->SetDaughter(1,"pi+");
                                                   >>  96   mode[5]->SetDaughter(2,"pi-");
                                                   >>  97   mode[5]->SetDaughter(3,"anti_nu_tau");
                                                   >>  98 
                                                   >>  99   for (G4int index=0; index <6; index++ ) table->Insert(mode[index]);  
                                                   >> 100 
                                                   >> 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 G4TauPlus G4TauPlus::theTauPlus(
                                                   >> 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+ + nu_mu + anti_nu_tau       << 
 90     mode = new G4TauLeptonicDecayChannel("tau+ << 
 91     table->Insert(mode);                       << 
 92     // tau+ -> e+ + nu_e + anti_nu_tau         << 
 93     mode = new G4TauLeptonicDecayChannel("tau+ << 
 94     table->Insert(mode);                       << 
 95     // tau+ -> pi+ + anti_nu_tau               << 
 96     mode = new G4PhaseSpaceDecayChannel("tau+" << 
 97     table->Insert(mode);                       << 
 98     // tau+ -> pi0 + pi0 + pi+ + anti_nu_tau   << 
 99     mode = new G4PhaseSpaceDecayChannel("tau+" << 
100     table->Insert(mode);                       << 
101     // tau+ -> pi0 + pi0 + pi+ + anti_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, "anti_nu_tau");       << 
110     table->Insert(mode);                       << 
111     // tau+ -> pi+ + pi+ + pi- + anti_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, "anti_nu_tau");       << 
120     table->Insert(mode);                       << 
121                                                << 
122     anInstance->SetDecayTable(table);          << 
123   }                                            << 
124   theInstance = static_cast<G4TauPlus*>(anInst << 
125   return theInstance;                          << 
126 }                                              << 
127                                                   124 
128 G4TauPlus* G4TauPlus::TauPlusDefinition()      << 125 G4TauPlus* G4TauPlus::TauPlusDefinition(){return &theTauPlus;}
129 {                                              << 
130   return Definition();                         << 
131 }                                              << 
132                                                   126 
133 G4TauPlus* G4TauPlus::TauPlus()                << 127 G4TauPlus*  G4TauPlus::TauPlus()
134 {                                              << 128 {  
135   return Definition();                         << 129   return &theTauPlus; 
136 }                                                 130 }
                                                   >> 131 
137                                                   132