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

update processes

Location:
trunk/source/processes/hadronic/models/de_excitation/evaporation/include
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4AlphaEvaporationChannel.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4AlphaEvaporationChannel.hh,v 1.3 2006/06/29 20:09:43 gunter Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4AlphaEvaporationChannel.hh,v 1.8 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    3131// by V. Lara (Nov. 1999)
    3232//
    33 
    3433
    3534#ifndef G4AlphaEvaporationChannel_h
    … …  
    6261private:
    6362
    64   G4AlphaCoulombBarrier theCoulombBarrier;
    65        
     63  G4AlphaCoulombBarrier theCoulombBarrier;     
     64
    6665  G4AlphaEvaporationProbability theEvaporationProbability;
    6766

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4AlphaEvaporationProbability.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4AlphaEvaporationProbability.hh,v 1.3 2006/06/29 20:09:45 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    31 // by V. Lara (Nov 1999)
     29// by V. Lara (Oct 1998)
    3230//
    33 
    34 
    35 
    3631#ifndef G4AlphaEvaporationProbability_h
    3732#define G4AlphaEvaporationProbability_h 1
    … …  
    3934
    4035#include "G4EvaporationProbability.hh"
    41 
     36#include "G4AlphaCoulombBarrier.hh"
    4237
    4338class G4AlphaEvaporationProbability : public G4EvaporationProbability
    … …  
    5550  G4bool operator==(const G4AlphaEvaporationProbability &right) const;
    5651  G4bool operator!=(const G4AlphaEvaporationProbability &right) const;
    57  
     52
    5853
    5954private:
    6055
    61   virtual G4double CalcAlphaParam(const G4Fragment & fragment) const
    62   { return 1.0 + CCoeficient(static_cast<G4double>(fragment.GetZ()-GetZ()));}
    63        
    64   virtual G4double CalcBetaParam(const G4Fragment & ) const
    65   { return 0.0; }
     56  virtual G4double CrossSection(const  G4Fragment & fragment, const  G4double K);
    6657
    67        
    68   G4double CCoeficient(const G4double ) const;
     58  G4double GetOpt0(const G4double K);
     59  G4double GetOpt12(const G4double K);
     60  G4double GetOpt34(const G4double K);
    6961
    70   // Excitation energy levels
    71   std::vector<G4double> ExcitEnergies;
    72   // Spin of excitation energy levels
    73   std::vector<G4int> ExcitSpins;
     62 
     63 virtual G4double CalcAlphaParam(const G4Fragment & fragment) ;
     64 
     65 virtual G4double CalcBetaParam(const G4Fragment & fragment)  ;
     66 
     67  G4double CCoeficient(const G4double aZ) ;
     68 
     69//data members
     70   
     71      G4AlphaCoulombBarrier theCoulombBarrier;
     72
     73      G4double ResidualA;
     74      G4double ResidualZ;
     75      G4double theA;
     76      G4double theZ;
     77      G4double ResidualAthrd;
     78      G4double FragmentA;
     79      G4double FragmentAthrd;
     80
     81
    7482};
    7583
    7684
    7785#endif
     86

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4DeuteronEvaporationChannel.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4DeuteronEvaporationChannel.hh,v 1.3 2006/06/29 20:09:47 gunter Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4DeuteronEvaporationChannel.hh,v 1.8 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    3131// by V. Lara (Nov. 1999)
    3232//
    33 
    3433
    3534#ifndef G4DeuteronEvaporationChannel_h

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4DeuteronEvaporationProbability.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4DeuteronEvaporationProbability.hh,v 1.3 2006/06/29 20:09:49 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    31 // by V. Lara (Nov 1999)
     29// by V. Lara (Oct 1998)
    3230//
    33 
    34 
    3531
    3632#ifndef G4DeuteronEvaporationProbability_h
    … …  
    3935
    4036#include "G4EvaporationProbability.hh"
    41 
     37#include "G4DeuteronCoulombBarrier.hh"
    4238
    4339class G4DeuteronEvaporationProbability : public G4EvaporationProbability
    … …  
    5551  G4bool operator==(const G4DeuteronEvaporationProbability &right) const;
    5652  G4bool operator!=(const G4DeuteronEvaporationProbability &right) const;
    57  
     53
    5854
    5955private:
    6056
    61   virtual G4double CalcAlphaParam(const G4Fragment & fragment) const
    62   { return 1.0 + CCoeficient(static_cast<G4double>(fragment.GetZ()-GetZ()));}
    63        
    64   virtual G4double CalcBetaParam(const G4Fragment & ) const
    65   { return 0.0; }
     57  virtual G4double CrossSection(const  G4Fragment & fragment, const  G4double K);
    6658
    67        
    68   G4double CCoeficient(const G4double aZ) const;
     59  G4double GetOpt0(const G4double K);
     60  G4double GetOpt12(const G4double K);
     61  G4double GetOpt34(const G4double K);
    6962
    70   // Excitation energy levels
    71   std::vector<G4double> ExcitEnergies;
    72   // Spin of excitation energy levels
    73   std::vector<G4int> ExcitSpins;
     63 
     64 virtual G4double CalcAlphaParam(const G4Fragment & fragment) ;
     65 
     66 virtual G4double CalcBetaParam(const G4Fragment & fragment) ;
     67 
     68  G4double CCoeficient(const G4double aZ) ;
     69 
     70//data members
     71
     72      G4DeuteronCoulombBarrier theCoulombBarrier;
     73
     74      G4double ResidualA;
     75      G4double ResidualZ;
     76      G4double theA;
     77      G4double theZ;
     78      G4double ResidualAthrd;
     79      G4double FragmentA;
     80      G4double FragmentAthrd;
     81
    7482
    7583};
     84
     85
    7686#endif
     87
     88
     89

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4Evaporation.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4Evaporation.hh,v 1.4 2007/02/14 13:37:49 ahoward Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4Evaporation.hh,v 1.5 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    … …  
    6969  void SetDefaultChannel();
    7070  void SetGEMChannel();
    71  
    72 private:
    7371
    7472#ifdef debug

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4EvaporationChannel.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4EvaporationChannel.hh,v 1.3 2006/06/29 20:09:53 gunter Exp $
    28 // GEANT4 tag $Name:  $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    … …  
    3937#include "G4VEmissionProbability.hh"
    4038#include "G4EvaporationProbability.hh"
     39#include "G4NeutronEvaporationProbability.hh"
     40#include "G4ProtonEvaporationProbability.hh"
     41#include "G4DeuteronEvaporationProbability.hh"
     42#include "G4TritonEvaporationProbability.hh"
     43#include "G4He3EvaporationProbability.hh"
     44#include "G4AlphaEvaporationProbability.hh"
    4145#include "G4VLevelDensityParameter.hh"
    4246#include "G4VCoulombBarrier.hh"
    … …  
    5155{
    5256public:
    53   // Available constructors
    54   G4EvaporationChannel(const G4int theA, const G4int theZ,
    55                        G4VEmissionProbability * aEmissionStrategy,
    56                        G4VCoulombBarrier * aCoulombBarrier);
     57  // constructor
     58
    5759
    5860  G4EvaporationChannel(const G4int theA, const G4int theZ, const G4String & aName,
    5961                       G4VEmissionProbability * aEmissionStrategy,
    60                        G4VCoulombBarrier * aCoulombBarrier);
    61 
    62   G4EvaporationChannel(const G4int theA, const G4int theZ, const G4String * aName,
    63                        G4VEmissionProbability * aEmissionStrategy,
    64                        G4VCoulombBarrier * aCoulombBarrier);
    65 
     62                       G4VCoulombBarrier * aCoulombBarrier);
    6663public:
    6764  // destructor
    … …  
    7067  void SetEmissionStrategy(G4VEmissionProbability * aEmissionStrategy)
    7168  {theEvaporationProbabilityPtr = aEmissionStrategy;}
    72  
     69
    7370  void SetCoulombBarrierStrategy(G4VCoulombBarrier * aCoulombBarrier)
    7471  {theCoulombBarrierPtr = aCoulombBarrier;}
     72 
     73
    7574 
    7675protected:
    … …  
    9392
    9493  G4FragmentVector * BreakUp(const G4Fragment & theNucleus);
    95  
    96 //   inline void SetEmissionStrategy(G4VEmissionProbability * aStrategy)
    97 //     {
    98 //       if (MyOwnEvaporationProbability) delete theEvaporationProbabilityPtr;
    99 //       theEvaporationProbabilityPtr = aStrategy;
    100 //       MyOwnEvaporationProbability = false;
    101 //     }
    10294
    103 
    104   inline void SetLevelDensityParameter(G4VLevelDensityParameter * aLevelDensity)
    105   {
    106     if (MyOwnLevelDensity) delete theLevelDensityPtr;
    107     theLevelDensityPtr = aLevelDensity;
    108     MyOwnLevelDensity = false;
    109   }
    110  
     95  // void SetLevelDensityParameter(G4VLevelDensityParameter * aLevelDensity);
     96 
    11197public:
    11298
    … …  
    119105  { return MaximalKineticEnergy; }
    120106 
    121   // ----------------------
    122  
    123107private:
    124108 
    … …  
    130114
    131115  // Samples fragment kinetic energy.
    132   G4double CalcKineticEnergy(void);
     116    G4double  GetKineticEnergy(const G4Fragment & aFragment);
    133117
    134118  // This has to be removed and put in Random Generator
    … …  
    142126  // They are intializated at object creation (constructor) time.
    143127
    144   // Atomic Number
    145   G4int A;
     128  // Atomic Number of ejectile
     129  G4int theA;
    146130
    147   // Charge
    148   G4int Z;
     131  // Charge of ejectile
     132  G4int theZ;
    149133
    150134
    151   // For evaporation probability calcualtion
     135  // For evaporation probability calcualation
    152136  G4VEmissionProbability * theEvaporationProbabilityPtr;
    153137
    154138  // For Level Density calculation
    155   G4bool MyOwnLevelDensity;
     139 // G4bool MyOwnLevelDensity;
    156140  G4VLevelDensityParameter * theLevelDensityPtr;
    157  
     141
     142
    158143  // For Coulomb Barrier calculation
    159144  G4VCoulombBarrier * theCoulombBarrierPtr;
    160145  G4double CoulombBarrier;
    161 
     146 
     147 
    162148  //---------------------------------------------------
    163149
    … …  
    166152  // the atomic number, charge and excitation energy of nucleus.
    167153
    168   // Residual Atomic Number
    169   G4int AResidual;
     154  // Residual Mass Number
     155  G4int ResidualA;
    170156
    171157  // Residual Charge
    172   G4int ZResidual;
    173 
    174 //   // Binding Energy
    175 //   G4double BindingEnergy;
    176 
    177 //      // Level Density Parameter
    178 //      G4double LevelDensityParameter;
     158  G4int ResidualZ;
    179159       
    180160  // Emission Probability
    … …  
    184164  // Maximal Kinetic Energy that can be carried by fragment
    185165  G4double MaximalKineticEnergy;
     166
     167
    186168};
    187169

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4EvaporationFactory.hh

    r819 r962  
    2626//
    2727// $Id: G4EvaporationFactory.hh,v 1.3 2006/06/29 20:09:55 gunter Exp $
    28 // GEANT4 tag $Name: $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4EvaporationProbability.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4EvaporationProbability.hh,v 1.3 2006/06/29 20:09:57 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    31 // by V. Lara (Oct 1998) 
     29// by V. Lara (Oct 1998)
    3230//
    33 
    34 
    35 
    3631#ifndef G4EvaporationProbability_h
    3732#define G4EvaporationProbability_h 1
    … …  
    4136#include "G4VLevelDensityParameter.hh"
    4237#include "G4EvaporationLevelDensityParameter.hh"
     38#include "G4VCoulombBarrier.hh"
     39#include "G4CoulombBarrier.hh"
    4340
    4441
    … …  
    4744public:
    4845  // Only available constructor
    49   G4EvaporationProbability(const G4int anA, const G4int aZ, const G4double aGamma) :
     46  G4EvaporationProbability(const G4int anA, const G4int aZ, const G4double aGamma,G4VCoulombBarrier * aCoulombBarrier) :
    5047    theA(anA),
    5148    theZ(aZ),
    52     Gamma(aGamma)
     49    Gamma(aGamma)
     50,  theCoulombBarrierptr(aCoulombBarrier)
    5351  {
    5452    theEvapLDPptr = new G4EvaporationLevelDensityParameter;
     53
     54   
    5555  }
    5656
    … …  
    5858  {
    5959    if (theEvapLDPptr != 0) delete theEvapLDPptr;
     60
    6061  }
    6162
    … …  
    6768
    6869protected:
    69 
    70   void SetExcitationEnergiesPtr(std::vector<G4double> * anExcitationEnergiesPtr)
    71   {ExcitationEnergies = anExcitationEnergiesPtr;}
    72 
    73   void SetExcitationSpinsPtr(std::vector<G4int> * anExcitationSpinsPtr)
    74   {ExcitationSpins = anExcitationSpinsPtr;}
    75 
    7670 
    7771  // Default constructor
    7872  G4EvaporationProbability() {}
     73
    7974private:
    8075  // Copy constructor
    … …  
    8681 
    8782public:
    88   G4double EmissionProbability(const G4Fragment & fragment, const G4double anEnergy);
     83
     84 G4double ProbabilityDistributionFunction( const G4Fragment & aFragment, const G4double K);
     85
     86 G4double EmissionProbability(const G4Fragment & fragment, const G4double anEnergy);
    8987
    9088private:
    9189
    92   G4double CalcProbability(const G4Fragment & fragment, const G4double MaximalKineticEnergy);
    93   virtual G4double CCoeficient(const G4double ) const {return 0.0;};
     90  G4double CalculateProbability(const G4Fragment & fragment, const G4double MaximalKineticEnergy );
    9491
    95   virtual G4double CalcAlphaParam(const G4Fragment & ) const {return 1.0;}
    96   virtual G4double CalcBetaParam(const G4Fragment & ) const {return 1.0;}
     92  G4double IntegrateEmissionProbability(const G4Fragment & aFragment, const G4double & Low, const G4double & Up );
     93
     94protected:
     95
     96 virtual G4double CrossSection( const  G4Fragment & fragment, const G4double K )= 0; 
     97
     98 virtual G4double CalcAlphaParam(const G4Fragment & fragment)=0 ;
     99 
     100 virtual G4double CalcBetaParam(const G4Fragment & fragment)=0 ;
     101
     102private:
    97103
    98104  // Data Members
    … …  
    107113  G4double Gamma;
    108114
    109   // Discrete Excitation Energies
    110   std::vector<G4double> * ExcitationEnergies;
     115//The Coulomb Barrier
     116         G4VCoulombBarrier * theCoulombBarrierptr;
    111117
    112   //
    113   std::vector<G4int> * ExcitationSpins;
    114118
    115119};
    116120
    117121
     122
     123
    118124#endif

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4He3EvaporationChannel.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4He3EvaporationChannel.hh,v 1.3 2006/06/29 20:09:59 gunter Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4He3EvaporationChannel.hh,v 1.8 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    3131// by V. Lara (Nov. 1999)
    3232//
    33 
    3433
    3534#ifndef G4He3EvaporationChannel_h
    … …  
    6160private:
    6261
    63   G4He3CoulombBarrier theCoulombBarrier;
     62    G4He3CoulombBarrier theCoulombBarrier;
    6463       
    6564  G4He3EvaporationProbability theEvaporationProbability;

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4He3EvaporationProbability.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4He3EvaporationProbability.hh,v 1.3 2006/06/29 20:10:01 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    31 // by V. Lara (Nov 1999)
     29// by V. Lara (Oct 1998)
    3230//
    33 
    34 
    35 
    3631#ifndef G4He3EvaporationProbability_h
    3732#define G4He3EvaporationProbability_h 1
    … …  
    3934
    4035#include "G4EvaporationProbability.hh"
    41 
     36#include "G4He3CoulombBarrier.hh"
    4237
    4338class G4He3EvaporationProbability : public G4EvaporationProbability
    … …  
    5550  G4bool operator==(const G4He3EvaporationProbability &right) const;
    5651  G4bool operator!=(const G4He3EvaporationProbability &right) const;
    57  
     52
    5853
    5954private:
    6055
    61   virtual G4double CalcAlphaParam(const G4Fragment & fragment) const
    62   { return 1.0 + CCoeficient(static_cast<G4double>(fragment.GetZ()-GetZ()));}
    63        
    64   virtual G4double CalcBetaParam(const G4Fragment & ) const
    65   { return 0.0; }
     56  virtual G4double CrossSection(const  G4Fragment & fragment, const  G4double K);
    6657
    67        
    68   G4double CCoeficient(const G4double aZ) const;
     58  G4double GetOpt0(const G4double K);
     59  G4double GetOpt12(const G4double K);
     60  G4double GetOpt34(const G4double K);
    6961
    70   // Excitation energy levels
    71   std::vector<G4double> ExcitEnergies;
    72   // Spin of excitation energy levels
    73   std::vector<G4int> ExcitSpins;
     62 
     63 virtual G4double CalcAlphaParam(const G4Fragment & fragment) ;
     64 
     65 virtual G4double CalcBetaParam(const G4Fragment & fragment) ;
     66 
     67  G4double CCoeficient(const G4double aZ) ;
     68 
     69//data members
     70   
     71      G4He3CoulombBarrier theCoulombBarrier;
     72
     73      G4double ResidualA;
     74      G4double ResidualZ;
     75      G4double theA;
     76      G4double theZ;
     77      G4double ResidualAthrd;
     78      G4double FragmentA;
     79      G4double FragmentAthrd;
     80
     81
    7482};
    7583
    7684
    7785#endif
     86
     87
     88
     89

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4NeutronEvaporationChannel.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4NeutronEvaporationChannel.hh,v 1.3 2006/06/29 20:10:03 gunter Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4NeutronEvaporationChannel.hh,v 1.8 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    3131// by V. Lara (Nov. 1999)
    3232//
    33 
    3433
    3534#ifndef G4NeutronEvaporationChannel_h
    … …  
    6160private:
    6261
    63   G4NeutronCoulombBarrier theCoulombBarrier;
     62 G4NeutronCoulombBarrier theCoulombBarrier;
    6463       
    6564  G4NeutronEvaporationProbability theEvaporationProbability;

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4NeutronEvaporationProbability.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4NeutronEvaporationProbability.hh,v 1.4 2006/06/29 20:10:05 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    31 // by V. Lara (Nov 1999)
     29// by V. Lara (Oct 1998)
    3230//
    33 
    34 
    3531
    3632#ifndef G4NeutronEvaporationProbability_h
    … …  
    3935
    4036#include "G4EvaporationProbability.hh"
    41 
     37#include "G4NeutronCoulombBarrier.hh"
    4238
    4339class G4NeutronEvaporationProbability : public G4EvaporationProbability
    4440{
    4541public:
    46   // Only available constructor
     42 
    4743  G4NeutronEvaporationProbability();
    4844               
    4945  ~G4NeutronEvaporationProbability() {}
    5046private: 
    51   // Copy constructor
     47 
    5248  G4NeutronEvaporationProbability(const G4NeutronEvaporationProbability &right);
    5349
    … …  
    5551  G4bool operator==(const G4NeutronEvaporationProbability &right) const;
    5652  G4bool operator!=(const G4NeutronEvaporationProbability &right) const;
    57  
    5853
    5954private:
    6055
    61   virtual G4double CalcAlphaParam(const G4Fragment & fragment) const
    62   { return 0.76+2.2/std::pow(static_cast<G4double>(fragment.GetA()-GetA()),1.0/3.0);}
    63        
    64   virtual G4double CalcBetaParam(const G4Fragment & fragment) const
    65   { return (2.12/std::pow(static_cast<G4double>(fragment.GetA()-GetA()),2.0/3.0) - 0.05)*MeV/
    66       CalcAlphaParam(fragment); }
     56  virtual G4double CrossSection(const  G4Fragment & fragment, const  G4double K);
    6757
    68   // Excitation energy levels
    69   std::vector<G4double> ExcitEnergies;
    70   // Spin of excitation energy levels
    71   std::vector<G4int> ExcitSpins;
     58  G4double GetOpt12(const G4double K);
     59  G4double GetOpt34(const G4double K);
     60
     61 virtual G4double CalcAlphaParam(const G4Fragment & fragment);
     62 
     63 virtual G4double CalcBetaParam(const G4Fragment & fragment);
     64 
     65 
     66//data members
     67
     68      G4NeutronCoulombBarrier theCoulombBarrier;
     69
     70      G4double ResidualA;
     71      G4double ResidualZ;
     72      G4double theA;
     73      G4double theZ;
     74      G4double ResidualAthrd;
     75      G4double FragmentA;
     76      G4double FragmentAthrd;
     77
     78
     79
    7280
    7381};

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4ProtonEvaporationChannel.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4ProtonEvaporationChannel.hh,v 1.3 2006/06/29 20:10:07 gunter Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4ProtonEvaporationChannel.hh,v 1.8 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    3131// by V. Lara (Nov. 1999)
    3232//
    33 
    3433
    3534#ifndef G4ProtonEvaporationChannel_h
    … …  
    6261private:
    6362
     63   G4ProtonCoulombBarrier  theCoulombBarrier;
     64 
    6465  G4ProtonEvaporationProbability  theEvaporationProbability;
    6566       
    66   G4ProtonCoulombBarrier  theCoulombBarrier;
     67
    6768};
    6869#endif

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4ProtonEvaporationProbability.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4ProtonEvaporationProbability.hh,v 1.3 2006/06/29 20:10:09 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    31 // by V. Lara (Nov 1999)
     29// by V. Lara (Oct 1998)
    3230//
    33 
    34 
    3531
    3632#ifndef G4ProtonEvaporationProbability_h
    … …  
    3935
    4036#include "G4EvaporationProbability.hh"
    41 
     37#include "G4ProtonCoulombBarrier.hh"
    4238
    4339class G4ProtonEvaporationProbability : public G4EvaporationProbability
    … …  
    5551  G4bool operator==(const G4ProtonEvaporationProbability &right) const;
    5652  G4bool operator!=(const G4ProtonEvaporationProbability &right) const;
    57  
     53
    5854
    5955private:
    6056
    61   virtual G4double CalcAlphaParam(const G4Fragment & fragment) const
    62   { return 1.0 + CCoeficient(static_cast<G4double>(fragment.GetZ()-GetZ()));}
    63        
    64   virtual G4double CalcBetaParam(const G4Fragment & ) const
    65   { return 0.0; }
     57  virtual G4double CrossSection(const  G4Fragment & fragment, const  G4double K);
    6658
    67        
    68   G4double CCoeficient(const G4double aZ) const;
     59  G4double GetOpt0(const G4double K);
     60  G4double GetOpt1(const G4double K);
     61  G4double GetOpt2(const G4double K);
     62  G4double GetOpt3(const G4double K);
    6963
    70   // Excitation energy levels
    71   std::vector<G4double> ExcitEnergies;
    72   // Spin of excitation energy levels
    73   std::vector<G4int> ExcitSpins;
     64 
     65 virtual G4double CalcAlphaParam(const G4Fragment & fragment)  ;
     66 
     67 virtual G4double CalcBetaParam(const G4Fragment & fragment) ;
     68 
     69  G4double CCoeficient(const G4double aZ) ;
     70 
     71//data members
     72
     73      G4ProtonCoulombBarrier theCoulombBarrier;
     74
     75      G4double ResidualA;
     76      G4double ResidualZ;
     77      G4double theA;
     78      G4double theZ;
     79      G4double ResidualAthrd;
     80      G4double FragmentA;
     81      G4double FragmentAthrd;
     82      G4double U;
     83
     84
     85
    7486
    7587};
     88
     89
    7690#endif

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4TritonEvaporationChannel.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4TritonEvaporationChannel.hh,v 1.3 2006/06/29 20:10:11 gunter Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4TritonEvaporationChannel.hh,v 1.8 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    3131// by V. Lara (Nov. 1999)
    3232//
    33 
    3433
    3534#ifndef G4TritonEvaporationChannel_h

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4TritonEvaporationProbability.hh

    r819 r962  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: G4TritonEvaporationProbability.hh,v 1.3 2006/06/29 20:10:13 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26//J.M. Quesada (August2008). Based on:
    2927//
    3028// Hadronic Process: Nuclear De-excitations
    31 // by V. Lara (Nov 1999)
     29// by V. Lara (Oct 1998)
    3230//
    33 
    34 
    35 
    3631#ifndef G4TritonEvaporationProbability_h
    3732#define G4TritonEvaporationProbability_h 1
    … …  
    3934
    4035#include "G4EvaporationProbability.hh"
    41 
     36#include "G4TritonCoulombBarrier.hh"
    4237
    4338class G4TritonEvaporationProbability : public G4EvaporationProbability
    … …  
    5550  G4bool operator==(const G4TritonEvaporationProbability &right) const;
    5651  G4bool operator!=(const G4TritonEvaporationProbability &right) const;
    57  
     52
    5853
    5954private:
    6055
    61   virtual G4double CalcAlphaParam(const G4Fragment & fragment) const
    62   { return 1.0 + CCoeficient(static_cast<G4double>(fragment.GetZ()-GetZ()));}
    63        
    64   virtual G4double CalcBetaParam(const G4Fragment & ) const
    65   { return 0.0; }
     56  virtual G4double CrossSection(const  G4Fragment & fragment, const  G4double K);
    6657
    67   G4double CCoeficient(const G4double aZ) const;
     58  G4double GetOpt0(const G4double K);
     59  G4double GetOpt12(const G4double K);
     60  G4double GetOpt34(const G4double K);
    6861
    69   // Excitation energy levels
    70   std::vector<G4double> ExcitEnergies;
    71   // Spin of excitation energy levels
    72   std::vector<G4int> ExcitSpins;
     62 
     63 virtual G4double CalcAlphaParam(const G4Fragment & fragment)  ;
     64 
     65 virtual G4double CalcBetaParam(const G4Fragment & fragment)  ;
     66 
     67  G4double CCoeficient(const G4double aZ) ;
     68 
     69//data members
     70
     71     
     72      G4TritonCoulombBarrier theCoulombBarrier;
     73
     74      G4double ResidualA;
     75      G4double ResidualZ;
     76      G4double theA;
     77      G4double theZ;
     78      G4double ResidualAthrd;
     79      G4double FragmentA;
     80      G4double FragmentAthrd;
     81
    7382
    7483};
    … …  
    7685
    7786#endif
     87
     88

  • trunk/source/processes/hadronic/models/de_excitation/evaporation/include/G4VEvaporation.hh

    r819 r962  
    2525//
    2626//
    27 // $Id: G4VEvaporation.hh,v 1.3 2006/06/29 20:10:15 gunter Exp $
    28 // GEANT4 tag $Name: $
     27// $Id: G4VEvaporation.hh,v 1.4 2008/09/19 13:32:54 ahoward Exp $
     28// GEANT4 tag $Name: geant4-09-02-ref-02 $
    2929//
    3030// Hadronic Process: Nuclear De-excitations
    3131// by V. Lara (Oct 1998) written from G4Evaporation.hh (May 1998)
    3232//
    33 
     33// Modif (03 September 2008) by J. M. Quesada for external choice of inverse
     34// cross section option
     35// JMQ (06 September 2008) Also external choices have been added for
     36// superimposed Coulomb barrier (if useSICBis set true, by default is false)
    3437
    3538
    … …  
    5659  virtual G4FragmentVector * BreakItUp(const G4Fragment &theNucleus) = 0;
    5760
     61  // for inverse cross section choice
     62  inline void SetOPTxs(G4int opt) { OPTxs = opt;}
     63  // for superimposed Coulomb Barrier for inverse cross sections       
     64  inline void UseSICB(G4bool use) { useSICB = use; }   
     65protected:
     66   G4int OPTxs;
     67   G4bool useSICB;
    5868
    5969};
Note: See TracChangeset for help on using the changeset viewer.