G4ShellVacancy.cc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 //
00027 // $Id: G4ShellVacancy.cc
00028 // GEANT4 tag $Name:
00029 //
00030 // Author: Elena Guardincerri (Elena.Guardincerri@ge.infn.it)
00031 //
00032 // History:
00033 // -----------
00034 // 21 Sept 2001 Elena Guardincerri     Created
00035 // 25 Mar  2002  V.Ivanchenko          Change AverageNOfIonisations int->double
00036 //
00037 // -------------------------------------------------------------------
00038 
00039 #include "G4ShellVacancy.hh"
00040 #include "G4MaterialCutsCouple.hh"
00041 #include "G4Material.hh"
00042 #include "G4Poisson.hh"
00043 #include "G4VEMDataSet.hh"
00044 
00045 G4ShellVacancy::G4ShellVacancy()
00046 
00047 { }
00048 
00049 G4ShellVacancy::~G4ShellVacancy()
00050 {
00051   G4int size = xsis.size();
00052   for (G4int k =0; k<size; k++)
00053     {
00054       delete xsis[k];
00055       xsis[k] = 0;
00056     }
00057 }
00058 
00059 void G4ShellVacancy::AddXsiTable(G4VEMDataSet* p)
00060 {
00061   xsis.push_back(p);
00062 }
00063 
00064 std::vector<G4int> G4ShellVacancy::GenerateNumberOfIonisations(const G4MaterialCutsCouple*
00065                                                                  couple,
00066                                                                  G4double
00067                                                                  incidentEnergy,
00068                                                                  G4double eLoss) const
00069 
00070 {
00071   std::vector<G4int> numberOfIonisations;
00072   const G4Material* material = couple->GetMaterial();
00073   G4int numberOfElements = material->GetNumberOfElements();
00074 
00075   for (G4int i = 0; i<numberOfElements; i++)
00076     {
00077       G4double averageNumberOfIonisations = AverageNOfIonisations(couple,
00078                                                                   i,
00079                                                                   incidentEnergy,
00080                                                                   eLoss);
00081       G4int ionisations = 0;
00082       if(averageNumberOfIonisations > 0.0) {
00083         ionisations = (G4int) G4Poisson(averageNumberOfIonisations);
00084       }
00085 
00086       numberOfIonisations.push_back(ionisations);
00087 
00088     }
00089   return numberOfIonisations;
00090 
00091 }
00092 
00093 G4double G4ShellVacancy::AverageNOfIonisations(const G4MaterialCutsCouple* couple,
00094                                                      G4int index,
00095                                                      G4double energy,
00096                                                      G4double eLoss) const
00097 
00098 {
00099   //  G4int indexOfElementInMaterial= -1;
00100 
00101   G4double averageEnergy = energy - eLoss/2.;
00102 
00103   size_t indexInMaterialTable = couple->GetIndex();
00104 
00105   G4VEMDataSet* aSetOfXsi = xsis[indexInMaterialTable];
00106 
00107   G4double aXsi = aSetOfXsi->FindValue(averageEnergy,index);
00108 
00109   return aXsi * eLoss;
00110 }

Generated on Mon May 27 17:49:49 2013 for Geant4 by  doxygen 1.4.7