Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/track/src/G4ParticleChangeForLoss.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 /track/src/G4ParticleChangeForLoss.cc (Version 11.3.0) and /track/src/G4ParticleChangeForLoss.cc (Version 9.5.p2)


  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 // G4ParticleChangeForLoss class implementatio << 
 27 //                                                 26 //
 28 // Author: Hisaya Kurashige, 23 March 1998     <<  27 // $Id: G4ParticleChangeForLoss.cc,v 1.18 2010-07-21 09:30:15 gcosmo Exp $
 29 // Revision: Vladimir Ivantchenko, 16 January  <<  28 // GEANT4 tag $Name: not supported by cvs2svn $
 30 // ------------------------------------------- <<  29 //
 31                                                <<  30 //
                                                   >>  31 // --------------------------------------------------------------
                                                   >>  32 //  GEANT 4 class implementation file
                                                   >>  33 //
                                                   >>  34 // ------------------------------------------------------------
                                                   >>  35 //   Implemented for the new scheme                 23 Mar. 1998  H.Kurahige
                                                   >>  36 // --------------------------------------------------------------
                                                   >>  37 //
                                                   >>  38 //   Modified:
                                                   >>  39 //   16.01.04 V.Ivanchenko update for model variant of energy loss
                                                   >>  40 //   15.04.05 V.Ivanchenko inline update methods
                                                   >>  41 //   28.08.06 V.Ivanchenko Add access to current track and polarizaion
                                                   >>  42 //
                                                   >>  43 // ------------------------------------------------------------
                                                   >>  44 //
 32 #include "G4ParticleChangeForLoss.hh"              45 #include "G4ParticleChangeForLoss.hh"
 33 #include "G4SystemOfUnits.hh"                  <<  46 #include "G4Track.hh"
                                                   >>  47 #include "G4Step.hh"
                                                   >>  48 #include "G4DynamicParticle.hh"
 34 #include "G4ExceptionSeverity.hh"                  49 #include "G4ExceptionSeverity.hh"
                                                   >>  50 #include "G4VelocityTable.hh"
 35                                                    51 
 36 // ------------------------------------------- << 
 37 G4ParticleChangeForLoss::G4ParticleChangeForLo     52 G4ParticleChangeForLoss::G4ParticleChangeForLoss()
                                                   >>  53   : G4VParticleChange(), currentTrack(0), proposedKinEnergy(0.),
                                                   >>  54     lowEnergyLimit(1.0*eV), currentCharge(0.)
 38 {                                                  55 {
 39   // Disable flag that is enabled in G4VPartic <<  56   theSteppingControlFlag = NormalCondition;
 40   debugFlag = false;                               57   debugFlag = false;
 41   SetNumberOfSecondaries(1);                   <<  58 #ifdef G4VERBOSE
                                                   >>  59   if (verboseLevel>2) {
                                                   >>  60     G4cout << "G4ParticleChangeForLoss::G4ParticleChangeForLoss() " << G4endl;
                                                   >>  61   }
                                                   >>  62 #endif
                                                   >>  63 }
                                                   >>  64 
                                                   >>  65 G4ParticleChangeForLoss::~G4ParticleChangeForLoss()
                                                   >>  66 {
                                                   >>  67 #ifdef G4VERBOSE
                                                   >>  68   if (verboseLevel>2) {
                                                   >>  69     G4cout << "G4ParticleChangeForLoss::~G4ParticleChangeForLoss() " << G4endl;
                                                   >>  70   }
                                                   >>  71 #endif
 42 }                                                  72 }
 43                                                    73 
 44 // ------------------------------------------- <<  74 G4ParticleChangeForLoss::
                                                   >>  75 G4ParticleChangeForLoss(const G4ParticleChangeForLoss &right)
                                                   >>  76   : G4VParticleChange(right)
                                                   >>  77 {
                                                   >>  78   if (verboseLevel>1) {
                                                   >>  79     G4cout << "G4ParticleChangeForLoss::  copy constructor is called " << G4endl;
                                                   >>  80   }
                                                   >>  81   currentTrack = right.currentTrack;
                                                   >>  82   proposedKinEnergy = right.proposedKinEnergy;
                                                   >>  83   lowEnergyLimit = right.lowEnergyLimit;
                                                   >>  84   currentCharge = right.currentCharge;
                                                   >>  85   proposedMomentumDirection = right.proposedMomentumDirection;
                                                   >>  86 }
                                                   >>  87 
                                                   >>  88 // assignment operator
                                                   >>  89 G4ParticleChangeForLoss & G4ParticleChangeForLoss::operator=(
                                                   >>  90                                    const G4ParticleChangeForLoss &right)
                                                   >>  91 {
                                                   >>  92 #ifdef G4VERBOSE
                                                   >>  93   if (verboseLevel>1) {
                                                   >>  94     G4cout << "G4ParticleChangeForLoss:: assignment operator is called " << G4endl;
                                                   >>  95   }
                                                   >>  96 #endif
                                                   >>  97 
                                                   >>  98   if (this != &right) {
                                                   >>  99     if (theNumberOfSecondaries>0) {
                                                   >> 100 #ifdef G4VERBOSE
                                                   >> 101       if (verboseLevel>0) {
                                                   >> 102   G4cout << "G4ParticleChangeForLoss: assignment operator Warning  ";
                                                   >> 103   G4cout << "theListOfSecondaries is not empty ";
                                                   >> 104       }
                                                   >> 105 #endif
                                                   >> 106        for (G4int index= 0; index<theNumberOfSecondaries; index++){
                                                   >> 107    if ( (*theListOfSecondaries)[index] ) delete (*theListOfSecondaries)[index] ;
                                                   >> 108        }
                                                   >> 109     }
                                                   >> 110     delete theListOfSecondaries; 
                                                   >> 111     theListOfSecondaries =  new G4TrackFastVector();
                                                   >> 112     theNumberOfSecondaries = right.theNumberOfSecondaries;
                                                   >> 113     for (G4int index = 0; index<theNumberOfSecondaries; index++){
                                                   >> 114       G4Track* newTrack =  new G4Track(*((*right.theListOfSecondaries)[index] ));
                                                   >> 115       theListOfSecondaries->SetElement(index, newTrack);          }
                                                   >> 116 
                                                   >> 117     theStatusChange = right.theStatusChange;
                                                   >> 118     theLocalEnergyDeposit = right.theLocalEnergyDeposit;
                                                   >> 119     theSteppingControlFlag = right.theSteppingControlFlag;
                                                   >> 120     theParentWeight = right.theParentWeight;
                                                   >> 121     isParentWeightProposed = right.isParentWeightProposed;
                                                   >> 122     isParentWeightSetByProcess = right.isParentWeightSetByProcess;
                                                   >> 123     fSetSecondaryWeightByProcess = right.fSetSecondaryWeightByProcess;
                                                   >> 124 
                                                   >> 125     currentTrack = right.currentTrack;
                                                   >> 126     proposedKinEnergy = right.proposedKinEnergy;
                                                   >> 127     currentCharge = right.currentCharge;
                                                   >> 128     proposedMomentumDirection = right.proposedMomentumDirection;
                                                   >> 129   }
                                                   >> 130   return *this;
                                                   >> 131 }
                                                   >> 132 
                                                   >> 133 //----------------------------------------------------------------
                                                   >> 134 // methods for printing messages
                                                   >> 135 //
                                                   >> 136 
 45 void G4ParticleChangeForLoss::DumpInfo() const    137 void G4ParticleChangeForLoss::DumpInfo() const
 46 {                                                 138 {
 47   // use base-class DumpInfo                   << 139 // use base-class DumpInfo
 48   G4VParticleChange::DumpInfo();                  140   G4VParticleChange::DumpInfo();
 49                                                   141 
 50   G4long oldprc = G4cout.precision(8);         << 142   G4int oldprc = G4cout.precision(3);
 51   G4cout << "      --------------------------- << 143   G4cout << "        Charge (eplus)   : "
 52   G4cout << "        G4ParticleChangeForLoss p << 144        << std::setw(20) << currentCharge/eplus
 53   G4cout << "        Charge (eplus)   : " << s << 145        << G4endl;
 54          << currentCharge / eplus << G4endl;   << 146   G4cout << "        Kinetic Energy (MeV): "
 55   G4cout << "        Kinetic Energy (MeV): " < << 147        << std::setw(20) << proposedKinEnergy/MeV
 56          << proposedKinEnergy / MeV << G4endl; << 148        << G4endl;
 57   G4cout << "        Momentum Direct - x : " < << 149   G4cout << "        Momentum Direct - x : "
 58          << proposedMomentumDirection.x() << G << 150        << std::setw(20) << proposedMomentumDirection.x()
 59   G4cout << "        Momentum Direct - y : " < << 151        << G4endl;
 60          << proposedMomentumDirection.y() << G << 152   G4cout << "        Momentum Direct - y : "
 61   G4cout << "        Momentum Direct - z : " < << 153        << std::setw(20) << proposedMomentumDirection.y()
 62          << proposedMomentumDirection.z() << G << 154        << G4endl;
                                                   >> 155   G4cout << "        Momentum Direct - z : "
                                                   >> 156        << std::setw(20) << proposedMomentumDirection.z()
                                                   >> 157        << G4endl;
 63   G4cout.precision(oldprc);                       158   G4cout.precision(oldprc);
 64 }                                                 159 }
 65                                                   160 
 66 // ------------------------------------------- << 161 G4bool G4ParticleChangeForLoss::CheckIt(const G4Track& aTrack)
 67 G4Step* G4ParticleChangeForLoss::UpdateStepFor << 
 68 {                                                 162 {
 69   const G4StepPoint* pPreStepPoint = pStep->Ge << 163   G4bool    itsOK = true;
 70   G4StepPoint* pPostStepPoint = pStep->GetPost << 164   G4bool    exitWithError = false;
 71                                                   165 
 72   // accumulate change of the kinetic energy   << 166   G4double  accuracy;
 73   G4double preKinEnergy = pPreStepPoint->GetKi << 
 74   G4double kinEnergy =                         << 
 75     pPostStepPoint->GetKineticEnergy() + (prop << 
 76                                                   167 
 77   pPostStepPoint->SetCharge(currentCharge);    << 168   // Energy should not be lager than initial value
                                                   >> 169   accuracy = ( proposedKinEnergy - aTrack.GetKineticEnergy())/MeV;
                                                   >> 170   if (accuracy > accuracyForWarning) {
                                                   >> 171 #ifdef G4VERBOSE
                                                   >> 172     G4cout << "G4ParticleChangeForLoss::CheckIt: ";
                                                   >> 173     G4cout << "KinEnergy become larger than the initial value!" << G4endl;
                                                   >> 174     G4cout << "  Difference:  " << accuracy  << "[MeV] " <<G4endl;
                                                   >> 175 #endif
                                                   >> 176     itsOK = false;
                                                   >> 177     if (accuracy > accuracyForException) exitWithError = true;
                                                   >> 178   }
 78                                                   179 
 79   // calculate velocity                        << 180   // dump out information of this particle change
 80   if(kinEnergy > 0.0)                          << 181 #ifdef G4VERBOSE
 81   {                                            << 182   if (!itsOK) {
 82     pPostStepPoint->SetKineticEnergy(kinEnergy << 183     G4cout << "G4ParticleChangeForLoss::CheckIt " << G4endl;
 83                                                << 184     DumpInfo();
 84     // assuming that mass>0, zero mass particl << 
 85     pPostStepPoint->SetVelocity(CLHEP::c_light << 
 86   }                                            << 
 87   else                                         << 
 88   {                                            << 
 89     pPostStepPoint->SetKineticEnergy(0.0);     << 
 90     pPostStepPoint->SetVelocity(0.0);          << 
 91   }                                            << 
 92                                                << 
 93   if(isParentWeightProposed)                   << 
 94   {                                            << 
 95     pPostStepPoint->SetWeight(theParentWeight) << 
 96   }                                               185   }
                                                   >> 186 #endif
 97                                                   187 
 98   pStep->AddTotalEnergyDeposit(theLocalEnergyD << 188   // Exit with error
 99   pStep->AddNonIonizingEnergyDeposit(theNonIon << 189   if (exitWithError) {
100   return pStep;                                << 190     G4Exception("G4ParticleChangeForLoss::CheckIt",
                                                   >> 191     "TRACK004", EventMustBeAborted,
                                                   >> 192     "energy was  illegal");
                                                   >> 193   }
                                                   >> 194 
                                                   >> 195   //correction
                                                   >> 196   if (!itsOK) {
                                                   >> 197     proposedKinEnergy = aTrack.GetKineticEnergy();
                                                   >> 198   }
                                                   >> 199 
                                                   >> 200   itsOK = (itsOK) && G4VParticleChange::CheckIt(aTrack);
                                                   >> 201   return itsOK;
101 }                                                 202 }
102                                                   203 
103 // ------------------------------------------- << 204 //----------------------------------------------------------------
104 G4Step* G4ParticleChangeForLoss::UpdateStepFor << 205 // methods for updating G4Step
                                                   >> 206 //
                                                   >> 207 
                                                   >> 208 G4Step* G4ParticleChangeForLoss::UpdateStepForAlongStep(G4Step* pStep)
105 {                                                 209 {
106   G4StepPoint* pPostStepPoint = pStep->GetPost    210   G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
107                                                   211 
108   pPostStepPoint->SetCharge(currentCharge);    << 212   // accumulate change of the kinetic energy
109   pPostStepPoint->SetMomentumDirection(propose << 213   G4double kinEnergy = pPostStepPoint->GetKineticEnergy() +
110   if(proposedKinEnergy > 0.0)                  << 214     (proposedKinEnergy - pStep->GetPreStepPoint()->GetKineticEnergy());
111   {                                            << 
112     pPostStepPoint->SetKineticEnergy(proposedK << 
113                                                   215 
114     // assuming that mass>0, zero mass particl << 216   // update kinetic energy and charge
115     pPostStepPoint->SetVelocity(CLHEP::c_light << 217   if (kinEnergy < lowEnergyLimit) {
                                                   >> 218     theLocalEnergyDeposit += kinEnergy;
                                                   >> 219     kinEnergy = 0.0;
                                                   >> 220   } else {
                                                   >> 221     pPostStepPoint->SetCharge( currentCharge );
116   }                                               222   }
117   else                                         << 223   pPostStepPoint->SetKineticEnergy( kinEnergy );
118   {                                            << 224   // calculate velocity
119     pPostStepPoint->SetKineticEnergy(0.0);     << 225   pStep->GetTrack()->SetKineticEnergy(kinEnergy); 
120     pPostStepPoint->SetVelocity(0.0);          << 226   pPostStepPoint->SetVelocity(pStep->GetTrack()->CalculateVelocity());
                                                   >> 227   pStep->GetTrack()->SetKineticEnergy(pStep->GetPreStepPoint()->GetKineticEnergy()); 
                                                   >> 228 
                                                   >> 229   if (isParentWeightProposed) {
                                                   >> 230     // update weight 
                                                   >> 231     G4StepPoint* pPreStepPoint = pStep->GetPreStepPoint();
                                                   >> 232     G4double newWeight= theParentWeight/(pPreStepPoint->GetWeight())
                                                   >> 233       * (pPostStepPoint->GetWeight());
                                                   >> 234     if (isParentWeightSetByProcess) pPostStepPoint->SetWeight( newWeight );
                                                   >> 235     if (!fSetSecondaryWeightByProcess) {    
                                                   >> 236       //Set weight of secondary tracks
                                                   >> 237       for (G4int index= 0; index<theNumberOfSecondaries; index++){
                                                   >> 238   if ( (*theListOfSecondaries)[index] ) {
                                                   >> 239     ((*theListOfSecondaries)[index])->SetWeight(newWeight); ;
                                                   >> 240   }
                                                   >> 241       }
                                                   >> 242     }
121   }                                               243   }
122   pPostStepPoint->SetPolarization(proposedPola << 
123                                                   244 
124   if(isParentWeightProposed)                   << 245   pStep->AddTotalEnergyDeposit( theLocalEnergyDeposit );
125   {                                            << 246   pStep->AddNonIonizingEnergyDeposit( theNonIonizingEnergyDeposit );
126     pPostStepPoint->SetWeight(theParentWeight) << 247   return pStep;
                                                   >> 248 }
                                                   >> 249 
                                                   >> 250 G4Step* G4ParticleChangeForLoss::UpdateStepForPostStep(G4Step* pStep)
                                                   >> 251 {
                                                   >> 252   G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
                                                   >> 253   pPostStepPoint->SetCharge( currentCharge );
                                                   >> 254   pPostStepPoint->SetMomentumDirection( proposedMomentumDirection );
                                                   >> 255   pPostStepPoint->SetKineticEnergy( proposedKinEnergy );
                                                   >> 256   pStep->GetTrack()->SetKineticEnergy( proposedKinEnergy );
                                                   >> 257   pPostStepPoint->SetVelocity(pStep->GetTrack()->CalculateVelocity());
                                                   >> 258   pPostStepPoint->SetPolarization( proposedPolarization );
                                                   >> 259 
                                                   >> 260   if (isParentWeightProposed) {
                                                   >> 261     if (isParentWeightSetByProcess) pPostStepPoint->SetWeight( theParentWeight );
                                                   >> 262     if (!fSetSecondaryWeightByProcess) {    
                                                   >> 263       // Set weight of secondary tracks
                                                   >> 264       for (G4int index= 0; index<theNumberOfSecondaries; index++){
                                                   >> 265   if ( (*theListOfSecondaries)[index] ) {
                                                   >> 266     ((*theListOfSecondaries)[index])->SetWeight( theParentWeight );
                                                   >> 267   }
                                                   >> 268       }
                                                   >> 269     }
127   }                                               270   }
128                                                   271 
129   pStep->AddTotalEnergyDeposit(theLocalEnergyD << 272   pStep->AddTotalEnergyDeposit( theLocalEnergyDeposit );
130   pStep->AddNonIonizingEnergyDeposit(theNonIon << 273   pStep->AddNonIonizingEnergyDeposit( theNonIonizingEnergyDeposit );
131   return pStep;                                   274   return pStep;
132 }                                                 275 }
                                                   >> 276 
133                                                   277