G4PenelopeAnnihilationModel.hh

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 // $Id$
00027 //
00028 // Author: Luciano Pandola
00029 //
00030 // History:
00031 // -----------
00032 // 29 Oct 2008   L. Pandola   1st implementation. Migration from EM process 
00033 //                            to EM model
00034 //
00035 // -------------------------------------------------------------------
00036 //
00037 // Class description:
00038 // Low Energy Electromagnetic Physics, Positron Annihilation 
00039 // with Penelope Model
00040 // -------------------------------------------------------------------
00041 
00042 #ifndef G4PENELOPEANNIHILATIONMODEL_HH
00043 #define G4PENELOPEANNIHILATIONMODEL_HH 1
00044 
00045 #include "globals.hh"
00046 #include "G4VEmModel.hh"
00047 #include "G4DataVector.hh"
00048 #include "G4ParticleChangeForGamma.hh"
00049 
00050 class G4ParticleDefinition;
00051 class G4DynamicParticle;
00052 class G4MaterialCutsCouple;
00053 class G4Material;
00054 
00055 class G4PenelopeAnnihilationModel : public G4VEmModel 
00056 {
00057 public:
00058   
00059   G4PenelopeAnnihilationModel(const G4ParticleDefinition* p=0,
00060                          const G4String& processName ="PenAnnih");
00061   
00062   virtual ~G4PenelopeAnnihilationModel();
00063 
00064   virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
00065 
00066   virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
00067                                               G4double kinEnergy,
00068                                               G4double Z,
00069                                               G4double A=0,
00070                                               G4double cut=0,
00071                                               G4double emax=DBL_MAX);
00072 
00073   virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00074                                  const G4MaterialCutsCouple*,
00075                                  const G4DynamicParticle*,
00076                                  G4double tmin,
00077                                  G4double maxEnergy);
00078 
00079   void SetVerbosityLevel(G4int lev){verboseLevel = lev;};
00080   G4int GetVerbosityLevel(){return verboseLevel;};
00081 
00082 protected:
00083   G4ParticleChangeForGamma* fParticleChange;
00084 
00085 private:
00086   G4double ComputeCrossSectionPerElectron(G4double energy);
00087 
00088   G4PenelopeAnnihilationModel & operator=(const G4PenelopeAnnihilationModel &right);
00089   G4PenelopeAnnihilationModel(const G4PenelopeAnnihilationModel&);
00090 
00091   G4int verboseLevel;
00092   G4bool isInitialised;
00093 
00094   G4double fIntrinsicLowEnergyLimit;
00095   G4double fIntrinsicHighEnergyLimit;
00096 
00097   //Stored here so it is not recalculated every time
00098   G4double fPielr2;
00099 };
00100 
00101 #endif
00102 

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