Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/particles/leptons/src/G4NeutrinoE.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/G4NeutrinoE.cc (Version 11.3.0) and /particles/leptons/src/G4NeutrinoE.cc (Version 1.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: G4NeutrinoE.cc,v 1.1.10.1 1999/12/07 20:49:48 gunter Exp $
                                                   >>   9 // GEANT4 tag $Name: geant4-01-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 <fstream.h>
                                                   >>  23 #include <iomanip.h>
 35                                                    24 
 36 #include "G4NeutrinoE.hh"                          25 #include "G4NeutrinoE.hh"
 37                                                    26 
 38 #include "G4ParticleTable.hh"                  <<  27 // ######################################################################
 39 #include "G4String.hh"                         <<  28 // ###                         NEUTRINO                               ###
 40 #include "G4SystemOfUnits.hh"                  <<  29 // ######################################################################
                                                   >>  30 
                                                   >>  31 G4NeutrinoE::G4NeutrinoE(
                                                   >>  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 }        
 41                                                    46 
 42 G4NeutrinoE* G4NeutrinoE::theInstance = nullpt <<  47 // ......................................................................
                                                   >>  48 // ...                 static member definitions                      ...
                                                   >>  49 // ......................................................................
                                                   >>  50 //     
                                                   >>  51 //    Arguments for constructor are as follows
                                                   >>  52 //               name             mass          width         charge
                                                   >>  53 //             2*spin           parity  C-conjugation
                                                   >>  54 //          2*Isospin       2*Isospin3       G-parity
                                                   >>  55 //               type    lepton number  baryon number   PDG encoding
                                                   >>  56 //             stable         lifetime    decay table 
 43                                                    57 
 44 G4NeutrinoE* G4NeutrinoE::Definition()         <<  58 G4NeutrinoE G4NeutrinoE::theNeutrinoE(
 45 {                                              <<  59          "nu_e",          0.0*MeV,       0.0*MeV,         0.0, 
 46   if (theInstance != nullptr) return theInstan << 
 47   const G4String name = "nu_e";                << 
 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,         60         1,               0,             0,          
 66         0,               0,             0,         61         0,               0,             0,             
 67        "lepton",               1,                  62        "lepton",               1,             0,           12,
 68      true,            -1.0,          nullptr,  <<  63      true,             0.0,          NULL
 69              false,           "e"              <<  64 );
 70               );                               <<  65 
 71     // clang-format on                         <<  66 G4NeutrinoE* G4NeutrinoE::NeutrinoEDefinition() {return &theNeutrinoE;}
 72   }                                            <<  67 // initialization for static cut values
 73   theInstance = static_cast<G4NeutrinoE*>(anIn <<  68 G4double   G4NeutrinoE::theNeutrinoELengthCut = -1.0;
 74   return theInstance;                          <<  69 G4double*  G4NeutrinoE::theNeutrinoEKineticEnergyCuts = NULL;
 75 }                                              << 
 76                                                    70 
 77 G4NeutrinoE* G4NeutrinoE::NeutrinoEDefinition( << 
 78 {                                              << 
 79   return Definition();                         << 
 80 }                                              << 
 81                                                    71 
 82 G4NeutrinoE* G4NeutrinoE::NeutrinoE()          <<  72 // **********************************************************************
                                                   >>  73 // **************************** SetCuts *********************************
                                                   >>  74 // **********************************************************************
                                                   >>  75 
                                                   >>  76 void G4NeutrinoE::SetCuts(G4double aCut)
 83 {                                                  77 {
 84   return Definition();                         <<  78   theCutInMaxInteractionLength = aCut;
                                                   >>  79 
                                                   >>  80   const G4MaterialTable* materialTable = G4Material::GetMaterialTable();
                                                   >>  81   // Create the vector of cuts in energy
                                                   >>  82   // corresponding to the stopping range cut
                                                   >>  83   if(theKineticEnergyCuts) delete [] theKineticEnergyCuts;
                                                   >>  84   theKineticEnergyCuts = new G4double [materialTable->length()];
                                                   >>  85 
                                                   >>  86   // Build range vector for every material, convert cut into energy-cut,
                                                   >>  87   // fill theKineticEnergyCuts and delete the range vector
                                                   >>  88   for (G4int J=0; J<materialTable->length(); J++)
                                                   >>  89   {
                                                   >>  90     G4Material* aMaterial = (*materialTable)[J];
                                                   >>  91     theKineticEnergyCuts[J] = 0.0*keV;
                                                   >>  92   }
                                                   >>  93   theNeutrinoELengthCut = theCutInMaxInteractionLength;  
                                                   >>  94   theNeutrinoEKineticEnergyCuts = theKineticEnergyCuts;
                                                   >>  95 
                                                   >>  96   // Rebuild the physics tables for every process for this particle type
                                                   >>  97   
 85 }                                                  98 }
                                                   >>  99 
                                                   >> 100 
                                                   >> 101 
 86                                                   102