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 // 27 /// \file GB04BOptnBremSplitting.cc 28 /// \brief Implementation of the GB04BOptnBrem 29 30 #include "GB04BOptnBremSplitting.hh" 31 32 #include "G4BiasingProcessInterface.hh" 33 #include "G4ParticleChangeForLoss.hh" 34 35 //....oooOO0OOooo........oooOO0OOooo........oo 36 37 GB04BOptnBremSplitting::GB04BOptnBremSplitting 38 : G4VBiasingOperation(name), fSplittingFacto 39 {} 40 41 //....oooOO0OOooo........oooOO0OOooo........oo 42 43 GB04BOptnBremSplitting::~GB04BOptnBremSplittin 44 45 G4VParticleChange* 46 GB04BOptnBremSplitting::ApplyFinalStateBiasing 47 48 { 49 // -- Collect brem. process (wrapped process 50 G4VParticleChange* processFinalState = 51 callingProcess->GetWrappedProcess()->PostS 52 53 // -- if no splitting requested, let the bre 54 // -- generated final state: 55 if (fSplittingFactor == 1) return processFin 56 57 // -- a special case here: the brem. process 58 // -- over the step due to energy loss by so 59 // -- returning an unchanged incoming electr 60 // -- We respect this correction, and if no 61 // -- case is happening: 62 if (processFinalState->GetNumberOfSecondarie 63 64 // -- Now start the biasing: 65 // -- - the electron state will be taken a 66 // -- process, hence the one stored in a 67 // -- This state will be stored in our f 68 // -- - the photon accompagnying the elect 69 // -- - we will then do fSplittingFactor - 70 // -- fSplittingFactor - 1 additionnal g 71 // -- fParticleChange object. 72 73 // -- We called the brem. process above. Its 74 // -- a "G4ParticleChangeForLoss" object. We 75 // -- methods of the concrete G4ParticleChan 76 G4ParticleChangeForLoss* actualParticleChang 77 78 fParticleChange.Initialize(*track); 79 80 // -- Store electron final state: 81 fParticleChange.ProposeTrackStatus(actualPar 82 fParticleChange.ProposeEnergy(actualParticle 83 fParticleChange.ProposeMomentumDirection(act 84 85 // -- Now deal with the gamma's: 86 // -- their common weight: 87 G4double gammaWeight = track->GetWeight() / 88 89 // -- inform we will have fSplittingFactor g 90 fParticleChange.SetNumberOfSecondaries(fSpli 91 92 // -- inform we take care of secondaries wei 93 // -- secondaries are by default given the p 94 fParticleChange.SetSecondaryWeightByProcess( 95 96 // -- Store first gamma: 97 G4Track* gammaTrack = actualParticleChange-> 98 gammaTrack->SetWeight(gammaWeight); 99 fParticleChange.AddSecondary(gammaTrack); 100 // -- and clean-up the brem. process particl 101 actualParticleChange->Clear(); 102 103 // -- now start the fSplittingFactor-1 calls 104 // -- related gamma: 105 G4int nCalls = 1; 106 while (nCalls < fSplittingFactor) { 107 // ( note: we don't need to cast to actual 108 // secondary particles are from base cla 109 processFinalState = callingProcess->GetWra 110 if (processFinalState->GetNumberOfSecondar 111 gammaTrack = processFinalState->GetSecon 112 gammaTrack->SetWeight(gammaWeight); 113 fParticleChange.AddSecondary(gammaTrack) 114 nCalls++; 115 } 116 // -- very rare special case: we ignore fo 117 else if (processFinalState->GetNumberOfSec 118 for (G4int i = 0; i < processFinalState- 119 delete processFinalState->GetSecondary 120 } 121 processFinalState->Clear(); 122 } 123 124 // -- we are done: 125 return &fParticleChange; 126 } 127 128 //....oooOO0OOooo........oooOO0OOooo........oo 129