| Geant4 Cross Reference |
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 // G4ParticleChangeForMSC class implementation <<
27 // 26 //
28 // Author: Hisaya Kurashige, 23 March 1998 << 27 // $Id$
29 // Revision: Vladimir Ivantchenko, 16 January << 28 //
30 // 23 August 2 << 29 //
31 // ------------------------------------------- << 30 // --------------------------------------------------------------
>> 31 // GEANT 4 class implementation file
>> 32 //
>> 33 //
>> 34 //
>> 35 // ------------------------------------------------------------
>> 36 // Implemented for the new scheme 23 Mar. 1998 H.Kurahige
>> 37 // Update for model variant of msc 16 Jan 2004 V.Ivanchenko
>> 38 // --------------------------------------------------------------
32 39
33 #include "G4ParticleChangeForMSC.hh" 40 #include "G4ParticleChangeForMSC.hh"
34 #include "G4SystemOfUnits.hh" 41 #include "G4SystemOfUnits.hh"
35 #include "G4Track.hh" 42 #include "G4Track.hh"
36 #include "G4Step.hh" 43 #include "G4Step.hh"
>> 44 #include "G4DynamicParticle.hh"
37 #include "G4ExceptionSeverity.hh" 45 #include "G4ExceptionSeverity.hh"
38 46
39 // ------------------------------------------- << 47 G4ParticleChangeForMSC::G4ParticleChangeForMSC()
40 G4ParticleChangeForMSC::G4ParticleChangeForMSC << 48 : G4VParticleChange()
>> 49 {
>> 50 #ifdef G4VERBOSE
>> 51 if (verboseLevel>2) {
>> 52 G4cout << "G4ParticleChangeForMSC::G4ParticleChangeForMSC() " << G4endl;
>> 53 }
>> 54 #endif
>> 55 }
>> 56
>> 57 G4ParticleChangeForMSC::~G4ParticleChangeForMSC()
>> 58 {
>> 59 #ifdef G4VERBOSE
>> 60 if (verboseLevel>2) {
>> 61 G4cout << "G4ParticleChangeForMSC::~G4ParticleChangeForMSC() " << G4endl;
>> 62 }
>> 63 #endif
>> 64 }
>> 65
>> 66 G4ParticleChangeForMSC::
>> 67 G4ParticleChangeForMSC(const G4ParticleChangeForMSC &right)
>> 68 : G4VParticleChange(right)
>> 69 {
>> 70 if (verboseLevel>1) {
>> 71 G4cout << "G4ParticleChangeForMSC:: copy constructor is called " << G4endl;
>> 72 }
>> 73 theMomentumDirection = right.theMomentumDirection;
>> 74 thePosition = right.thePosition;
>> 75 }
>> 76
>> 77 // assignment operator
>> 78 G4ParticleChangeForMSC & G4ParticleChangeForMSC::operator=(
>> 79 const G4ParticleChangeForMSC &right)
41 { 80 {
42 // Disable flag that is enabled in G4VPartic << 81 if (verboseLevel>1) {
43 debugFlag = false; << 82 G4cout << "G4ParticleChangeForMSC:: assignment operator is called " << G4endl;
>> 83 }
>> 84 if (this != &right)
>> 85 {
>> 86 theListOfSecondaries = right.theListOfSecondaries;
>> 87 theSizeOftheListOfSecondaries = right.theSizeOftheListOfSecondaries;
>> 88 theNumberOfSecondaries = right.theNumberOfSecondaries;
>> 89 theStatusChange = right.theStatusChange;
>> 90 theLocalEnergyDeposit = right.theLocalEnergyDeposit;
>> 91 theSteppingControlFlag = right.theSteppingControlFlag;
>> 92 theTrueStepLength = right.theTrueStepLength;
>> 93
>> 94 theMomentumDirection = right.theMomentumDirection;
>> 95 thePosition = right.thePosition;
>> 96 }
>> 97 return *this;
44 } 98 }
45 99
46 // ------------------------------------------- << 100 //----------------------------------------------------------------
>> 101 // methods for updating G4Step
>> 102 //
>> 103
47 G4Step* G4ParticleChangeForMSC::UpdateStepForA 104 G4Step* G4ParticleChangeForMSC::UpdateStepForAlongStep(G4Step* pStep)
48 { 105 {
49 // update the G4Step specific attributes << 106 // Update the G4Step specific attributes
50 pStep->SetStepLength(theTrueStepLength); 107 pStep->SetStepLength(theTrueStepLength);
>> 108 theStatusChange = pStep->GetTrack()->GetTrackStatus();
>> 109
>> 110 // Multiple scattering calculates the final state of the particle
>> 111 G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
>> 112
>> 113 // update momentum direction
>> 114 pPostStepPoint->SetMomentumDirection(theMomentumDirection);
>> 115
>> 116 // update position
>> 117 pPostStepPoint->SetPosition( thePosition );
>> 118 return pStep;
>> 119 }
51 120
52 // multiple scattering calculates the final << 121 G4Step* G4ParticleChangeForMSC::UpdateStepForPostStep(G4Step* pStep)
>> 122 {
>> 123 // Multiple scattering calculates the final state of the particle
53 G4StepPoint* pPostStepPoint = pStep->GetPost 124 G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
54 125
55 // update momentum direction 126 // update momentum direction
56 pPostStepPoint->SetMomentumDirection(theMome 127 pPostStepPoint->SetMomentumDirection(theMomentumDirection);
57 128
58 // update position 129 // update position
59 pPostStepPoint->SetPosition(thePosition); << 130 pPostStepPoint->SetPosition( thePosition );
>> 131
>> 132 // Update the G4Step specific attributes
60 return pStep; 133 return pStep;
61 } 134 }
>> 135
>> 136
>> 137 //----------------------------------------------------------------
>> 138 // methods for printing messages
>> 139 //
>> 140
>> 141 void G4ParticleChangeForMSC::DumpInfo() const
>> 142 {
>> 143 // use base-class DumpInfo
>> 144 G4VParticleChange::DumpInfo();
>> 145
>> 146 G4int oldprc = G4cout.precision(3);
>> 147 G4cout << " Position - x (mm) : "
>> 148 << std::setw(20) << thePosition.x()/mm
>> 149 << G4endl;
>> 150 G4cout << " Position - y (mm) : "
>> 151 << std::setw(20) << thePosition.y()/mm
>> 152 << G4endl;
>> 153 G4cout << " Position - z (mm) : "
>> 154 << std::setw(20) << thePosition.z()/mm
>> 155 << G4endl;
>> 156 G4cout << " Momentum Direct - x : "
>> 157 << std::setw(20) << theMomentumDirection.x()
>> 158 << G4endl;
>> 159 G4cout << " Momentum Direct - y : "
>> 160 << std::setw(20) << theMomentumDirection.y()
>> 161 << G4endl;
>> 162 G4cout << " Momentum Direct - z : "
>> 163 << std::setw(20) << theMomentumDirection.z()
>> 164 << G4endl;
>> 165 G4cout.precision(oldprc);
>> 166 }
>> 167
>> 168
>> 169 G4bool G4ParticleChangeForMSC::CheckIt(const G4Track& aTrack)
>> 170 {
>> 171 G4bool itsOK = true;
>> 172 G4bool exitWithError = false;
>> 173
>> 174 G4double accuracy;
>> 175
>> 176 // check
>> 177
>> 178 // MomentumDirection should be unit vector
>> 179 accuracy = std::fabs(theMomentumDirection.mag2()-1.0);
>> 180 if (accuracy > accuracyForWarning) {
>> 181 itsOK = false;
>> 182 exitWithError = (accuracy > accuracyForException);
>> 183 #ifdef G4VERBOSE
>> 184 G4cout << " G4ParticleChangeForMSC::CheckIt : ";
>> 185 G4cout << "the Momentum Change is not unit vector !!"
>> 186 << " Difference: " << accuracy << G4endl;
>> 187 G4cout << aTrack.GetDefinition()->GetParticleName()
>> 188 << " E=" << aTrack.GetKineticEnergy()/MeV
>> 189 << " pos=" << aTrack.GetPosition().x()/m
>> 190 << ", " << aTrack.GetPosition().y()/m
>> 191 << ", " << aTrack.GetPosition().z()/m
>> 192 <<G4endl;
>> 193 #endif
>> 194 }
>> 195
>> 196 // dump out information of this particle change
>> 197 #ifdef G4VERBOSE
>> 198 if (!itsOK) DumpInfo();
>> 199 #endif
>> 200
>> 201 // Exit with error
>> 202 if (exitWithError) {
>> 203 G4Exception("G4ParticleChangeForMSC::CheckIt",
>> 204 "300",
>> 205 EventMustBeAborted,
>> 206 "momentum direction was illegal");
>> 207 }
>> 208 //correction
>> 209 if (!itsOK) {
>> 210 G4double vmag = theMomentumDirection.mag();
>> 211 theMomentumDirection = (1./vmag)*theMomentumDirection;
>> 212 }
>> 213
>> 214 itsOK = (itsOK) && G4VParticleChange::CheckIt(aTrack);
>> 215 return itsOK;
>> 216 }
>> 217
>> 218
>> 219
62 220
63 221