include/Geant4/G4EquilibriumEvaporator.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 // 20100413  M. Kelsey -- Pass G4CollisionOutput by ref to ::collide()
0028 // 20100517  M. Kelsey -- Inherit from common base class, make other colliders
0029 //      simple data members.  Rename timeToBigBang() to override
0030 //      base explosion().
0031 // 20100714  M. Kelsey -- Switch to new G4CascadeColliderBase class
0032 // 20100923  M. Kelsey -- Migrate to integer A and Z
0033 // 20100925  M. Kelsey -- Remove no longer necessary explosion() interface
0034 // 20110801  M. Kelsey -- Move "parms" buffer to data member, allocate in
0035 //      constructor.
0036 // 20110809  M. Kelsey -- Move "foutput" buffer to data member
0037 // 20130129  M. Kelsey -- Move QF interpolation to global statics for
0038 //      multi-threaded shared memory.
0039 // 20130620  Address Coverity complaint about missing copy actions
0040 // 20130621  Follow base class change to explosion() reference interface
0041 // 20130622  Inherit from G4CascadeDeexciteBase, move to deExcite() interface
0042 //      with G4Fragment
0043 // 20130808  Use new object-version of paraMaker, for thread safety
0044 // 20131001  Move QFinterpolation to data member to be thread local (not shared)
0045
0046 #ifndef G4EQUILIBRIUM_EVAPORATOR_HH
0047 #define G4EQUILIBRIUM_EVAPORATOR_HH
0048
0049 #include "G4CascadeDeexciteBase.hh"
0050 #include "G4CascadeInterpolator.hh"
0051 #include "G4CollisionOutput.hh"
0052 #include "G4Fissioner.hh"
0053 #include "G4BigBanger.hh"
0054 #include "G4InuclSpecialFunctions.hh"
0055
0056
0057 class G4EquilibriumEvaporator : public G4CascadeDeexciteBase {
0058 public:
0059   G4EquilibriumEvaporator();
0060   virtual ~G4EquilibriumEvaporator();
0061
0062   virtual void setVerboseLevel(G4int verbose);
0063
0064   virtual void deExcite(const G4Fragment& target, G4CollisionOutput& output);
0065
0066 private:
0067   // Replace base class verision with more complex conditions
0068   virtual G4bool explosion(G4int a, G4int z, G4double e) const;
0069
0070   // NOTE:  Must redeclare base-class polymorphisms
0071   virtual G4bool explosion(const G4Fragment& target) const {
0072     return G4CascadeDeexciteBase::explosion(target);
0073   }
0074
0075   G4bool goodRemnant(G4int a, G4int z) const;
0076
0077   G4InuclSpecialFunctions::paraMaker theParaMaker;
0078   G4double getE0(G4int A) const;
0079   G4double getPARLEVDEN(G4int A, G4int Z) const;
0080   G4double getQF(G4double x, G4double x2, G4int a, G4int z, G4double e) const;
0081   G4double getAF(G4double x, G4int a, G4int z, G4double e) const;
0082
0083   // Buffer for parameter sets
0084   std::pair<std::vector<G4double>, std::vector<G4double> > parms;
0085   G4CollisionOutput fission_output;
0086
0087   // Interpolation object for QF
0088   G4CascadeInterpolator<72> QFinterp;
0089
0090   G4Fissioner theFissioner;
0091   G4BigBanger theBigBanger;
0092
0093 private:
0094   // Copying of modules is forbidden
0095   G4EquilibriumEvaporator(const G4EquilibriumEvaporator&);
0096   G4EquilibriumEvaporator& operator=(const G4EquilibriumEvaporator&);
0097 };
0098
0099 #endif /* G4EQUILIBRIUM_EVAPORATOR_HH */
0100
0101