include/Geant4/G4PreCompoundModel.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 // by V. Lara
0027 //
0028 // Class Description
0029 // Model implementation for pre-equilibrium decay models in geant4.
0030 // To be used in your physics list, in case you neeed this kind of physics.
0031 // Can be used as a stand-allone model, but also in conjunction with an intra-nuclear
0032 // transport, or any of the string-parton models.
0033 // Class Description - End
0034 //
0035 // Modified:
0036 // 03.09.2008 J.M.Quesada added external choice of inverse
0037 //            cross section option.(default OPTxs=3)
0038 // 06.09.2008 J.M.Quesada external choices have been added for:
0039 //     - superimposed Coulomb barrier (if useSICB=true, default false)
0040 //     - "never go back"  hipothesis (if useNGB=true, default false)
0041 //     - soft cutoff from preeq. to equlibrium (if useSCO=true, default false)
0042 //                - CEM transition probabilities (if useCEMtr=true)
0043 // 30.10.2009 J.M.Quesada CEM transition probabilities are set as default
0044 // 20.08.2010 V.Ivanchenko Cleanup of the code
0045 // 03.01.2012 V.Ivanchenko Added pointer to G4ExcitationHandler to the
0046 //                         constructor
0047
0048 #ifndef G4PreCompoundModel_h
0049 #define G4PreCompoundModel_h 1
0050
0051 #include "G4VPreCompoundModel.hh"
0052 #include "G4Fragment.hh"
0053 #include "G4ReactionProductVector.hh"
0054 #include "G4ReactionProduct.hh"
0055 #include "G4ExcitationHandler.hh"
0056
0057 class G4PreCompoundEmission;
0058 class G4VPreCompoundTransitions;
0059 class G4NuclearLevelData;
0060 class G4ParticleDefinition;
0061
0062 class G4PreCompoundModel : public G4VPreCompoundModel
0063 {
0064 public:
0065
0066   explicit G4PreCompoundModel(G4ExcitationHandler* ptr = nullptr);
0067
0068   virtual ~G4PreCompoundModel();
0069
0070   virtual G4HadFinalState * ApplyYourself(const G4HadProjectile & thePrimary,
0071                       G4Nucleus & theNucleus) final;
0072
0073   virtual G4ReactionProductVector* DeExcite(G4Fragment& aFragment) final;
0074
0075   virtual void BuildPhysicsTable(const G4ParticleDefinition&) final;
0076
0077   virtual void InitialiseModel() final;
0078
0079   virtual void ModelDescription(std::ostream& outFile) const final;
0080   virtual void DeExciteModelDescription(std::ostream& outFile) const final;
0081
0082 private:
0083
0084   inline
0085   void PerformEquilibriumEmission(const G4Fragment & aFragment,
0086                   G4ReactionProductVector * result) const;
0087
0088   G4PreCompoundModel(const G4PreCompoundModel &) = delete;
0089   const G4PreCompoundModel& operator=(const G4PreCompoundModel &right) = delete;
0090   G4bool operator==(const G4PreCompoundModel &right) const = delete;
0091   G4bool operator!=(const G4PreCompoundModel &right) const = delete;
0092
0093   G4PreCompoundEmission* theEmission = nullptr;
0094   G4VPreCompoundTransitions* theTransition = nullptr;
0095   G4NuclearLevelData* fNuclData = nullptr;
0096
0097   const G4ParticleDefinition* proton;
0098   const G4ParticleDefinition* neutron;
0099
0100   G4double fLowLimitExc = 0.0;
0101   G4double fHighLimitExc = DBL_MAX;
0102
0103   G4bool useSCO = false;
0104   G4bool isInitialised = false;
0105   G4bool isActive = true;
0106
0107   G4int minZ = 3;
0108   G4int minA = 5;
0109   G4int modelID = -1;
0110
0111   G4HadFinalState theResult;
0112 };
0113
0114 inline void G4PreCompoundModel::PerformEquilibriumEmission(
0115             const G4Fragment & aFragment,
0116             G4ReactionProductVector * result) const
0117 {
0118   G4ReactionProductVector* deexResult =
0119     GetExcitationHandler()->BreakItUp(aFragment);
0120   result->insert(result->end(),deexResult->begin(), deexResult->end());
0121   delete deexResult;
0122 }
0123
0124 #endif
0125