Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 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 // 26 #include "G4ChannelingOptrMultiParticleChangeC 27 #include "G4ChannelingOptrChangeCrossSection.h 28 #include "G4BiasingProcessInterface.hh" 29 30 #include "G4ProcessManager.hh" 31 #include "G4ParticleDefinition.hh" 32 #include "G4ParticleTable.hh" 33 34 #include "G4SystemOfUnits.hh" 35 36 //....oooOO0OOooo........oooOO0OOooo........oo 37 38 G4ChannelingOptrMultiParticleChangeCrossSectio 39 G4VBiasingOperator("ChannelingChangeXS-Many"), 40 fCurrentOperator(0), 41 fnInteractions(0){ 42 AddChargedParticles(); 43 } 44 45 //....oooOO0OOooo........oooOO0OOooo........oo 46 47 void G4ChannelingOptrMultiParticleChangeCrossS 48 const G4ParticleDefinition* particle = 49 G4ParticleTable::GetParticleTable()->FindP 50 51 if ( particle == 0 ) 52 { 53 G4ExceptionDescription ed; 54 ed << "Particle `" << particleName << "' 55 G4Exception("G4ChannelingOptrMultiPartic 56 "G4Channeling", 57 JustWarning, 58 ed); 59 return; 60 } 61 62 G4ChannelingOptrChangeCrossSection* optr = 63 new G4ChannelingOptrChangeCrossSection(par 64 fParticlesToBias.push_back( particle ); 65 fBOptrForParticle[ particle ] = optr; 66 } 67 68 //....oooOO0OOooo........oooOO0OOooo........oo 69 70 void G4ChannelingOptrMultiParticleChangeCrossS 71 G4ParticleTable::G4PTblDicIterator* aParti 72 (G4ParticleTable::GetParticleTable())- 73 74 aParticleIterator->reset(); 75 76 while( (*aParticleIterator)() ){ 77 G4ParticleDefinition* particle = aPart 78 if (particle->GetPDGCharge() !=0) { 79 AddParticle(particle->GetParticleN 80 } 81 } 82 } 83 84 //....oooOO0OOooo........oooOO0OOooo........oo 85 86 G4VBiasingOperation* 87 G4ChannelingOptrMultiParticleChangeCrossSectio 88 ProposeOccurenceBiasingOperation(const G4Track 89 const G4Biasi 90 91 if ( fCurrentOperator ) return fCurrentOpera 92 GetProposedOccuren 93 else return 0; 94 } 95 96 //....oooOO0OOooo........oooOO0OOooo........oo 97 98 void G4ChannelingOptrMultiParticleChangeCrossS 99 const G4ParticleDefinition* definition = tra 100 std::map < const G4ParticleDefinition*, G4Ch 101 it = fBOptrForParticle.find( definition ); 102 fCurrentOperator = 0; 103 if ( it != fBOptrForParticle.end() ) fCurren 104 fnInteractions = 0; 105 } 106 107 //....oooOO0OOooo........oooOO0OOooo........oo 108 109 void 110 G4ChannelingOptrMultiParticleChangeCrossSectio 111 OperationApplied( const G4BiasingProcessInterf 112 G4BiasingAppliedCase 113 G4VBiasingOperation* 114 G4double 115 G4VBiasingOperation* 116 const G4VParticleChange* 117 fnInteractions++; 118 if ( fCurrentOperator ) fCurrentOperator->Re 119 120 121 122 123 124 } 125 126 //....oooOO0OOooo........oooOO0OOooo........oo 127