G4AdjointhIonisationModel.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: G4AdjointhIonisationModel.hh 69844 2013-05-16 09:19:33Z gcosmo $
00027 //
00029 //      Module:         G4AdjointhIonisationModel
00030 //      Author:         L. Desorgher
00031 //      Organisation:   SpaceIT GmbH
00032 //      Contract:       ESA contract 21435/08/NL/AT
00033 //      Customer:       ESA/ESTEC
00035 //
00036 // CHANGE HISTORY
00037 // --------------
00038 //      ChangeHistory: 
00039 //              13th February 2009 creation by L. Desorgher  
00040 //              10 November 2009   Implementation of the rapid sampling.                        
00041 //
00042 //-------------------------------------------------------------
00043 //      Documentation:
00044 //              Adjoint EM model for discrete reverse hadron ionisation. Tested at the moment only for protons.
00045 //
00046 
00047 #ifndef G4AdjointhIonisationModel_h
00048 #define G4AdjointhIonisationModel_h 1
00049 
00050 #include "globals.hh"
00051 #include "G4DynamicParticle.hh"
00052 #include "G4ParticleDefinition.hh"
00053 #include "G4MaterialCutsCouple.hh"
00054 #include "G4Material.hh"
00055 #include "G4Element.hh"
00056 #include "G4ElementVector.hh"
00057 #include "Randomize.hh"
00058 #include "G4ParticleDefinition.hh"
00059 #include "G4VEmModel.hh"
00060 #include "G4Electron.hh"
00061 #include "G4Gamma.hh"
00062 #include "G4ProductionCutsTable.hh"
00063 #include "G4VEmAdjointModel.hh"
00064 class G4PhysicsTable;
00065 class G4Region;
00066 class G4VParticleChange;
00067 class G4ParticleChange;
00068 class G4Track;
00069 class G4AdjointCSMatrix;
00070 
00071 
00072 class G4AdjointhIonisationModel: public G4VEmAdjointModel
00073 {
00074 
00075 public:
00076 
00077   G4AdjointhIonisationModel(G4ParticleDefinition* projectileDefinition);
00078 
00079   virtual ~G4AdjointhIonisationModel();
00080 
00081   virtual void SampleSecondaries(const G4Track& aTrack,
00082                                 G4bool IsScatProjToProjCase,
00083                                 G4ParticleChange* fParticleChange);
00084   void RapidSampleSecondaries(const G4Track& aTrack,
00085                                 G4bool IsScatProjToProjCase,
00086                                 G4ParticleChange* fParticleChange);
00087   virtual G4double DiffCrossSectionPerAtomPrimToSecond(
00088                                       G4double kinEnergyProj,  // kinetic energy of the primary particle before the interaction 
00089                                       G4double kinEnergyProd, // kinetic energy of the secondary particle 
00090                                       G4double Z, 
00091                                       G4double A = 0.);
00092  
00093   virtual G4double AdjointCrossSection(const G4MaterialCutsCouple* aCouple,
00094                                              G4double primEnergy,
00095                                              G4bool IsScatProjToProjCase);                                    
00096  
00097   //Set/Get methods
00098   //------------------
00099   
00100   virtual G4double GetSecondAdjEnergyMaxForScatProjToProjCase(G4double PrimAdjEnergy);
00101   virtual G4double GetSecondAdjEnergyMinForScatProjToProjCase(G4double PrimAdjEnergy,G4double Tcut=0);
00102   virtual G4double GetSecondAdjEnergyMaxForProdToProjCase(G4double PrimAdjEnergy);
00103   virtual G4double GetSecondAdjEnergyMinForProdToProjCase(G4double PrimAdjEnergy);
00104  
00105 
00106 private: //Methods
00107      
00108     
00109   void DefineProjectileProperty();
00110   
00111   //projectile property
00112   G4double mass;
00113   G4double tlimit;
00114   G4double spin;
00115   G4double magMoment2;
00116   G4double chargeSquare;
00117   G4double ratio, ratio2; 
00118   G4double one_plus_ratio_2;
00119   G4double formfact;
00120   G4bool   isIon;
00121   G4double one_minus_ratio_2;
00122   
00123   
00124 
00125   G4VEmModel* theBraggDirectEMModel;    
00126   //G4double term_Cross1, term_Cross2;
00127 
00128 
00129 
00130    
00131    
00132 
00133   
00134 };
00135 
00136 
00137 #endif
00138 

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