Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2025-01-18 09:58:15

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 //J.M. Quesada (August2008). Based on:
0027 //
0028 // Hadronic Process: Nuclear De-excitations
0029 // by V. Lara (Oct 1998)
0030 //
0031 // 17-11-2010 V.Ivanchenko in constructor replace G4VEmissionProbability by 
0032 //            G4EvaporationProbability and do not new and delete probability
0033 //            object at each call; use G4Pow
0034 // 16-07-2019 V.Ivanchenko use C++11 
0035 
0036 #ifndef G4EvaporationChannel_h
0037 #define G4EvaporationChannel_h 1
0038 
0039 #include "G4VEvaporationChannel.hh"
0040 
0041 class G4EvaporationProbability;
0042 class G4CoulombBarrier;
0043 class G4NuclearLevelData;
0044 
0045 class G4EvaporationChannel : public G4VEvaporationChannel
0046 {
0047 public:
0048 
0049   explicit G4EvaporationChannel(G4int A, G4int Z, 
0050                         G4EvaporationProbability*);
0051 
0052   ~G4EvaporationChannel() override;
0053 
0054   void Initialise() override;
0055 
0056   G4double GetEmissionProbability(G4Fragment* fragment) override; 
0057   
0058   G4Fragment* EmittedFragment(G4Fragment* theNucleus) override;
0059 
0060   G4double ComputeInverseXSection(G4Fragment*, G4double kinEnergy) override;
0061 
0062   G4double ComputeProbability(G4Fragment*, G4double kinEnergy) override;
0063 
0064   inline G4int GetZ() const { return theZ; };
0065 
0066   inline G4int GetA() const { return theA; };
0067 
0068   inline G4EvaporationProbability* GetEvaporationProbability()
0069   { return theProbability; }
0070 
0071   G4EvaporationChannel(const G4EvaporationChannel & right) = delete;
0072   const G4EvaporationChannel & operator=
0073   (const G4EvaporationChannel & right) = delete;
0074   G4bool operator==(const G4EvaporationChannel & right) const = delete;
0075   G4bool operator!=(const G4EvaporationChannel & right) const = delete;
0076 
0077 private: 
0078 
0079   G4NuclearLevelData* theLevelData;
0080 
0081   // For evaporation probability calcualation
0082   G4EvaporationProbability* theProbability;
0083 
0084   // For Coulomb Barrier calculation
0085   G4CoulombBarrier* theCoulombBarrier;
0086 
0087   // This data member define the channel. 
0088   // They are initialised at object creation (constructor) time.
0089   G4int theA;
0090   G4int theZ;
0091   G4int resA = 0;
0092   G4int resZ = 0;
0093 
0094   G4int secID;  // Creator model ID for this model
0095   
0096   G4double mass = 0.0;
0097   G4double resMass = 0.0;
0098   G4double ekinmax = 0.0;
0099   G4double bCoulomb = 0.0;
0100   G4double evapMass;
0101   G4double evapMass2;
0102 };
0103 
0104 #endif