Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/util/src/G4HadFinalState.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 /processes/hadronic/util/src/G4HadFinalState.cc (Version 11.3.0) and /processes/hadronic/util/src/G4HadFinalState.cc (Version 10.0.p4)


  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 // Modifications:                                  26 // Modifications:
 27 // 20110810  M. Kelsey -- Store secondaries by     27 // 20110810  M. Kelsey -- Store secondaries by value, not by pointer.
 28 //    Improve constness of argument passing.       28 //    Improve constness of argument passing.  Fix up some
 29 //    functions to avoid creating temporaries.     29 //    functions to avoid creating temporaries.
 30                                                    30 
 31 #include "G4HadFinalState.hh"                      31 #include "G4HadFinalState.hh"
 32 #include "G4HadronicException.hh"                  32 #include "G4HadronicException.hh"
 33                                                    33 
 34                                                    34 
 35 G4HadFinalState::G4HadFinalState()                 35 G4HadFinalState::G4HadFinalState()
 36   : theDirection(0,0,1), theEnergy(-1), theSta     36   : theDirection(0,0,1), theEnergy(-1), theStat(isAlive), 
 37     theW(1.), theEDep(0.) {}                       37     theW(1.), theEDep(0.) {}
 38                                                    38 
                                                   >>  39 
                                                   >>  40 G4int G4HadFinalState::GetNumberOfSecondaries() const {return theSecs.size();}
                                                   >>  41 
 39 void G4HadFinalState::SetEnergyChange(G4double     42 void G4HadFinalState::SetEnergyChange(G4double anEnergy) 
 40 {                                                  43 {
 41   theEnergy=anEnergy;                              44   theEnergy=anEnergy;
 42   if(theEnergy<0)                                  45   if(theEnergy<0) 
 43     {                                              46     {
 44       G4cout << "Final state energy was: E = " <<  47       std::cout << "Final state energy was: E = "<<theEnergy<<G4endl;
 45       throw G4HadronicException(__FILE__, __LI <<  48       throw G4HadronicException(__FILE__, __LINE__, "G4HadFinalState: fatal - negative energy");
 46       "G4HadFinalState: fatal - negative energ << 
 47     }                                              49     }
 48 }                                                  50 }
 49                                                    51 
                                                   >>  52 G4double G4HadFinalState::GetEnergyChange() const {return theEnergy;}
                                                   >>  53 
                                                   >>  54 void G4HadFinalState::SetMomentumChange(const G4ThreeVector& aV) {
                                                   >>  55   theDirection=aV;
                                                   >>  56 }
                                                   >>  57 
 50 void G4HadFinalState::SetMomentumChange(G4doub     58 void G4HadFinalState::SetMomentumChange(G4double x, G4double y, G4double z) 
 51 {                                                  59 {
 52   theDirection.set(x,y,z);                         60   theDirection.set(x,y,z);
 53   if(std::abs(x*x + y*y + z*z - 1.0)>0.001) {  <<  61   if(std::fabs(theDirection.mag()-1)>0.001) 
 54     G4cout <<"We have negative theDirection.ma <<  62     {
 55      <<G4endl;                                 <<  63       G4cout <<"We have negative theDirection.mag() = "<<theDirection.mag()<<G4endl;
 56     throw G4HadronicException(__FILE__, __LINE <<  64       throw G4HadronicException(__FILE__, __LINE__, "G4HadFinalState: fatal - negative direction.mag().");
 57     "G4HadFinalState: fatal - negative directi <<  65     }
 58   }                                            <<  66 }
                                                   >>  67 
                                                   >>  68 const G4ThreeVector& G4HadFinalState::GetMomentumChange() const {return theDirection;}
                                                   >>  69 
                                                   >>  70 void G4HadFinalState::AddSecondary(G4DynamicParticle *aP) {
                                                   >>  71   // NOTE:  In-situ constructor will be optimized away (no copying)
                                                   >>  72   if (aP) theSecs.push_back(G4HadSecondary(aP));
 59 }                                                  73 }
 60                                                    74 
 61 // Concatenate lists efficiently                   75 // Concatenate lists efficiently
 62 void G4HadFinalState::AddSecondaries(const std     76 void G4HadFinalState::AddSecondaries(const std::vector<G4HadSecondary>& addSecs)
 63 {                                                  77 {
 64   theSecs.insert(theSecs.end(),addSecs.begin()     78   theSecs.insert(theSecs.end(),addSecs.begin(),addSecs.end());
 65 }                                                  79 }
 66                                                    80 
                                                   >>  81 void G4HadFinalState::SetStatusChange(G4HadFinalStateStatus aS){theStat=aS;}
                                                   >>  82 
                                                   >>  83 G4HadFinalStateStatus G4HadFinalState::GetStatusChange() const {return theStat;}
                                                   >>  84 
                                                   >>  85 void G4HadFinalState::ClearSecondaries() { 
                                                   >>  86   theSecs.clear();
                                                   >>  87 }
                                                   >>  88 
 67 void G4HadFinalState::Clear()                      89 void G4HadFinalState::Clear()
 68 {                                                  90 {
 69   theDirection.set(0,0,1);                         91   theDirection.set(0,0,1);
 70   theEnergy = -1;                                  92   theEnergy = -1;
 71   theStat = isAlive;                               93   theStat = isAlive;
 72   theW = 1.;                                       94   theW = 1.;
 73   theEDep = 0.;                                    95   theEDep = 0.;
 74   ClearSecondaries();                              96   ClearSecondaries();
 75 }                                                  97 }
 76                                                    98 
 77 //void G4HadFinalState::SecondariesAreStale()  <<  99 void G4HadFinalState::SecondariesAreStale() { /*DEPRECATED*/ }
                                                   >> 100 
                                                   >> 101 const G4LorentzRotation& G4HadFinalState::GetTrafoToLab() const {return theT;}
                                                   >> 102 
                                                   >> 103 void G4HadFinalState::SetTrafoToLab(const G4LorentzRotation & aT) {theT = aT;}
                                                   >> 104 
                                                   >> 105 void G4HadFinalState::SetWeightChange(G4double aW){ theW=aW;}
                                                   >> 106 
                                                   >> 107 G4double G4HadFinalState::GetWeightChange() const {return theW;}
 78                                                   108 
 79 G4HadSecondary * G4HadFinalState::GetSecondary    109 G4HadSecondary * G4HadFinalState::GetSecondary(size_t i) 
 80 {                                                 110 {
 81   if(i>theSecs.size()) {                       << 111   if(i>theSecs.size())
 82     throw G4HadronicException(__FILE__, __LINE << 112     {
 83     "Trying direct access to secondary beyond  << 113       throw G4HadronicException(__FILE__, __LINE__, 
 84   }                                            << 114         "Trying direct access to secondary beyond end of list");
                                                   >> 115     }
 85   return &theSecs[i];                             116   return &theSecs[i];
 86 }                                                 117 }
 87                                                   118 
 88 const G4HadSecondary* G4HadFinalState::GetSeco    119 const G4HadSecondary* G4HadFinalState::GetSecondary(size_t i) const
 89 {                                                 120 {
 90   if(i>theSecs.size()) {                       << 121   if(i>theSecs.size())
 91     throw G4HadronicException(__FILE__, __LINE << 122     {
 92     "Trying direct access to secondary beyond  << 123       throw G4HadronicException(__FILE__, __LINE__, 
 93   }                                            << 124         "Trying direct access to secondary beyond end of list");
                                                   >> 125     }
 94   return &theSecs[i];                             126   return &theSecs[i];
 95 }                                                 127 }
                                                   >> 128 
                                                   >> 129 void G4HadFinalState::SetLocalEnergyDeposit(G4double aE) {theEDep=aE;}
                                                   >> 130 
                                                   >> 131 G4double G4HadFinalState::GetLocalEnergyDeposit() const {return theEDep;}
 96                                                   132