Ignore:
Timestamp:
Apr 6, 2009, 12:21:12 PM (15 years ago)
Author:
garnier
Message:

update processes

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/source/processes/electromagnetic/lowenergy/src/G4FinalStateExcitationMillerGreen.cc

    r819 r961  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4FinalStateExcitationMillerGreen.cc,v 1.2 2007/11/09 20:11:04 pia Exp $
    28 // GEANT4 tag $Name:  $
    29 //
    30 // Contact Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch)
    31 //
    32 // Reference: TNS Geant4-DNA paper
    33 // Reference for implementation model: NIM. 155, pp. 145-156, 1978
    34 
    35 // History:
    36 // -----------
    37 // Date         Name              Modification
    38 // 28 Apr 2007  M.G. Pia          Created in compliance with design described in TNS paper
    39 //
    40 // -------------------------------------------------------------------
    41 
    42 // Class description:
    43 // Reference: TNS Geant4-DNA paper
    44 // S. Chauvie et al., Geant4 physics processes for microdosimetry simulation:
    45 // design foundation and implementation of the first set of models,
    46 // IEEE Trans. Nucl. Sci., vol. 54, no. 6, Dec. 2007.
    47 // Further documentation available from http://www.ge.infn.it/geant4/dna
    48 
    49 // -------------------------------------------------------------------
    50 
     26// $Id: G4FinalStateExcitationMillerGreen.cc,v 1.3 2008/07/14 20:47:34 sincerti Exp $
     27// GEANT4 tag $Name: geant4-09-02-ref-02 $
    5128
    5229#include "G4FinalStateExcitationMillerGreen.hh"
    53 #include "G4Track.hh"
    54 #include "G4Step.hh"
    55 #include "G4DynamicParticle.hh"
    56 #include "Randomize.hh"
    5730
    58 #include "G4ParticleTypes.hh"
    59 #include "G4ParticleDefinition.hh"
    60 #include "G4Electron.hh"
    61 #include "G4SystemOfUnits.hh"
    62 #include "G4ParticleMomentum.hh"
     31//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    6332
    6433G4FinalStateExcitationMillerGreen::G4FinalStateExcitationMillerGreen()
    6534{
    66   name = "ExcitationMillerGreen";
    6735  lowEnergyLimit = 10 * eV;
    6836  highEnergyLimit = 10 * MeV;
    6937}
    7038
     39//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    7140
    7241G4FinalStateExcitationMillerGreen::~G4FinalStateExcitationMillerGreen()
    7342{}
    7443 
     44//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    7545
    7646const G4FinalStateProduct& G4FinalStateExcitationMillerGreen::GenerateFinalState(const G4Track& track, const G4Step& /* step */)
    7747{
    78   // Clear previous secondaries, energy deposit and particle kill status
    7948  product.Clear();
    8049
    8150  const G4DynamicParticle* particle = track.GetDynamicParticle();
    8251
    83   // Kinetic energy of primary particle
    8452  G4double k = particle->GetKineticEnergy();
    8553
    86   // Select excitation level on the basis of partial excitation cross section
    8754  G4int level = cross.RandomSelect(k,track.GetDefinition());
    88   // Excitation energy corresponding to the selected level
    8955  G4double excitationEnergy = waterStructure.ExcitationEnergy(level);
    9056  G4double newEnergy = k - excitationEnergy;
    9157 
    92   // ---- SI ---- Test on newEnergy
    9358  if (newEnergy > 0)
    94     {
    95       // Deposit excitation energy locally, modify primary energy accordingly
    96       // Particle direction is unchanged
    97       product.ModifyPrimaryParticle(particle->GetMomentumDirection(),newEnergy);
    98       product.AddEnergyDeposit(excitationEnergy);
    99     }
    100 
    101   // ---- SI ---- Particle is not modified by default otherwise
    102 /*
    103   else
    104     {
    105       // Primary particle is killed
    106       product.KillPrimaryParticle();
    107     }
    108 */
     59  {
     60    product.ModifyPrimaryParticle(particle->GetMomentumDirection(),newEnergy);
     61    product.AddEnergyDeposit(excitationEnergy);
     62  }
    10963
    11064  return product;
Note: See TracChangeset for help on using the changeset viewer.