Geant4 Cross Reference

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


                                                   >>   1 // This code implementation is the intellectual property of
                                                   >>   2 // the GEANT4 collaboration.
  1 //                                                  3 //
  2 // ******************************************* <<   4 // By copying, distributing or modifying the Program (or any work
  3 // * License and Disclaimer                    <<   5 // based on the Program) you indicate your acceptance of this statement,
  4 // *                                           <<   6 // and all its terms.
  5 // * The  Geant4 software  is  copyright of th << 
  6 // * the Geant4 Collaboration.  It is provided << 
  7 // * conditions of the Geant4 Software License << 
  8 // * LICENSE and available at  http://cern.ch/ << 
  9 // * include a list of copyright holders.      << 
 10 // *                                           << 
 11 // * Neither the authors of this software syst << 
 12 // * institutes,nor the agencies providing fin << 
 13 // * work  make  any representation or  warran << 
 14 // * regarding  this  software system or assum << 
 15 // * use.  Please see the license in the file  << 
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                           << 
 18 // * This  code  implementation is the result  << 
 19 // * technical work of the GEANT4 collaboratio << 
 20 // * By using,  copying,  modifying or  distri << 
 21 // * any work based  on the software)  you  ag << 
 22 // * use  in  resulting  scientific  publicati << 
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // ******************************************* << 
 25 //                                                  7 //
                                                   >>   8 // $Id: G4NeutrinoMu.cc,v 1.3 2000/02/27 06:23:41 kurasige Exp $
                                                   >>   9 // GEANT4 tag $Name: geant4-02-00 $
                                                   >>  10 //
                                                   >>  11 // 
 26 // -------------------------------------------     12 // ----------------------------------------------------------------------
 27 //      GEANT 4 class implementation file          13 //      GEANT 4 class implementation file
 28 //                                                 14 //
                                                   >>  15 //      For information related to this code contact:
                                                   >>  16 //      CERN, CN Division, ASD Group
 29 //      History: first implementation, based o     17 //      History: first implementation, based on object model of
 30 //      4th April 1996, G.Cosmo                <<  18 //      4th April 1996, G.Cosmo 
 31 //      by             H.Kurashige,7  July 199     19 //      by             H.Kurashige,7  July 1996
 32 // *******************************************     20 // **********************************************************************
 33 //  New impelemenataion as an utility class  M <<  21 
 34 // ------------------------------------------- <<  22 #include "g4std/fstream"
                                                   >>  23 #include "g4std/iomanip"
 35                                                    24 
 36 #include "G4NeutrinoMu.hh"                         25 #include "G4NeutrinoMu.hh"
 37                                                    26 
 38 #include "G4ParticleTable.hh"                  <<  27 // ######################################################################
 39 #include "G4String.hh"                         <<  28 // ###                     MU NEUTRINO                               ###
 40 #include "G4SystemOfUnits.hh"                  <<  29 // ######################################################################
                                                   >>  30 
                                                   >>  31 G4NeutrinoMu::G4NeutrinoMu(
                                                   >>  32        const G4String&     aName,        G4double            mass,
                                                   >>  33        G4double            width,        G4double            charge,   
                                                   >>  34        G4int               iSpin,        G4int               iParity,    
                                                   >>  35        G4int               iConjugation, G4int               iIsospin,   
                                                   >>  36        G4int               iIsospin3,    G4int               gParity,
                                                   >>  37        const G4String&     pType,        G4int               lepton,      
                                                   >>  38        G4int               baryon,       G4int               encoding,
                                                   >>  39        G4bool              stable,       G4double            lifetime,
                                                   >>  40        G4DecayTable        *decaytable )
                                                   >>  41  : G4VLepton( aName,mass,width,charge,iSpin,iParity,
                                                   >>  42               iConjugation,iIsospin,iIsospin3,gParity,pType,
                                                   >>  43               lepton,baryon,encoding,stable,lifetime,decaytable )
                                                   >>  44 {
                                                   >>  45   SetParticleSubType("mu");
                                                   >>  46 }        
 41                                                    47 
 42 G4NeutrinoMu* G4NeutrinoMu::theInstance = null <<  48 // ......................................................................
                                                   >>  49 // ...                 static member definitions                      ...
                                                   >>  50 // ......................................................................
                                                   >>  51 //     
                                                   >>  52 //    Arguments for constructor are as follows
                                                   >>  53 //               name             mass          width         charge
                                                   >>  54 //             2*spin           parity  C-conjugation
                                                   >>  55 //          2*Isospin       2*Isospin3       G-parity
                                                   >>  56 //               type    lepton number  baryon number   PDG encoding
                                                   >>  57 //             stable         lifetime    decay table 
 43                                                    58 
 44 G4NeutrinoMu* G4NeutrinoMu::Definition()       <<  59 G4NeutrinoMu G4NeutrinoMu::theNeutrinoMu(
 45 {                                              <<  60               "nu_mu",         0.0*MeV,       0.0*MeV,         0.0, 
 46   if (theInstance != nullptr) return theInstan << 
 47   const G4String name = "nu_mu";               << 
 48   // search in particle table]                 << 
 49   G4ParticleTable* pTable = G4ParticleTable::G << 
 50   G4ParticleDefinition* anInstance = pTable->F << 
 51   if (anInstance == nullptr) {                 << 
 52     // create particle                         << 
 53     //                                         << 
 54     //    Arguments for constructor are as fol << 
 55     //               name             mass     << 
 56     //             2*spin           parity  C- << 
 57     //          2*Isospin       2*Isospin3     << 
 58     //               type    lepton number  ba << 
 59     //             stable         lifetime     << 
 60     //             shortlived      subType     << 
 61                                                << 
 62     // clang-format off                        << 
 63    anInstance = new G4ParticleDefinition(      << 
 64                  name,         0.0*MeV,        << 
 65         1,               0,             0,         61         1,               0,             0,          
 66         0,               0,             0,         62         0,               0,             0,             
 67        "lepton",               1,                  63        "lepton",               1,             0,           14,
 68      true,            -1.0,          nullptr,  <<  64      true,             0.0,          NULL
 69              false,           "mu"             <<  65 );
 70               );                               << 
 71     // clang-format on                         << 
 72   }                                            << 
 73   theInstance = static_cast<G4NeutrinoMu*>(anI << 
 74   return theInstance;                          << 
 75 }                                              << 
 76                                                    66 
 77 G4NeutrinoMu* G4NeutrinoMu::NeutrinoMuDefiniti     67 G4NeutrinoMu* G4NeutrinoMu::NeutrinoMuDefinition()
 78 {                                                  68 {
 79   return Definition();                         <<  69   return &theNeutrinoMu;
 80 }                                                  70 }
                                                   >>  71 G4double   G4NeutrinoMu::theNeutrinoMuLengthCut = -1.0;
                                                   >>  72 G4double*  G4NeutrinoMu::theNeutrinoMuKineticEnergyCuts = NULL;
 81                                                    73 
 82 G4NeutrinoMu* G4NeutrinoMu::NeutrinoMu()       <<  74 // **********************************************************************
                                                   >>  75 // **************************** SetCuts *********************************
                                                   >>  76 // **********************************************************************
                                                   >>  77 void G4NeutrinoMu::SetCuts(G4double aCut)
 83 {                                                  78 {
 84   return Definition();                         <<  79   theCutInMaxInteractionLength = aCut;
                                                   >>  80 
                                                   >>  81   const G4MaterialTable* materialTable = G4Material::GetMaterialTable();
                                                   >>  82   // Create the vector of cuts in energy
                                                   >>  83   // corresponding to the stopping range cut
                                                   >>  84   if(theKineticEnergyCuts) delete [] theKineticEnergyCuts;
                                                   >>  85   theKineticEnergyCuts = new G4double [materialTable->length()];
                                                   >>  86 
                                                   >>  87   // Build range vector for every material, convert cut into energy-cut,
                                                   >>  88   // fill theKineticEnergyCuts and delete the range vector
                                                   >>  89   for (G4int J=0; J<materialTable->length(); J++)
                                                   >>  90   {
                                                   >>  91     G4Material* aMaterial = (*materialTable)[J];
                                                   >>  92     theKineticEnergyCuts[J] = 0.0*keV;
                                                   >>  93   }
                                                   >>  94   theNeutrinoMuLengthCut = theCutInMaxInteractionLength;  
                                                   >>  95   theNeutrinoMuKineticEnergyCuts = theKineticEnergyCuts;
                                                   >>  96   // Rebuild the physics tables for every process for this particle type
                                                   >>  97   
 85 }                                                  98 }
                                                   >>  99 
 86                                                   100