Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/nudex/include/G4NuDEXNeutronCaptureModel.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 //
 27 // -------------------------------------------------------------------
 28 //
 29 //      GEANT4 header file 
 30 //
 31 //      File name:     G4NuDEXNeutronCaptureModel
 32 //
 33 //      Author:        E.Mendoza & A.Ribon
 34 // 
 35 //      Creation date: 29 May 2024
 36 //
 37 //      Description:   This class (a proxy of the class G4NuDEX) uses
 38 //                     the NuDEX model to produce gammas and internal
 39 //                     conversion electrons from neutron capture.
 40 //                     Whenever NuDEX is not applicable, G4PhotonEvaporation
 41 //                     is used.
 42 //                     The implementation of this class follows the code
 43 //                     of the class G4NeutronRadCapture.
 44 //
 45 //      Modifications:
 46 //      
 47 // -------------------------------------------------------------------
 48 //
 49 // Class to use NuDEX model inside Geant4
 50 // 
 51 
 52 #ifndef G4NUDEXNEUTRONCAPTUREMODEL_HH
 53 #define G4NUDEXNEUTRONCAPTUREMODEL_HH 1
 54 
 55 #include "globals.hh"
 56 #include "G4HadronicInteraction.hh"
 57 #include "G4HadProjectile.hh"
 58 #include "G4Nucleus.hh"
 59 
 60 class G4NuDEXStatisticalNucleus;
 61 class G4VEvaporationChannel;
 62 
 63 
 64 #define G4NUDEX_MAXZA 120000
 65 
 66 
 67 class G4NuDEXNeutronCaptureModel : public G4HadronicInteraction {
 68   public:
 69     explicit G4NuDEXNeutronCaptureModel();
 70     virtual ~G4NuDEXNeutronCaptureModel();
 71   
 72     virtual G4HadFinalState* ApplyYourself( const G4HadProjectile &aTrack, G4Nucleus &targetNucleus ) final;
 73     virtual void InitialiseModel() final;
 74 
 75   private:
 76     G4NuDEXNeutronCaptureModel & operator=( const G4NuDEXNeutronCaptureModel &right ) = delete;
 77     G4NuDEXNeutronCaptureModel( const G4NuDEXNeutronCaptureModel& ) = delete;
 78 
 79     G4int GenerateNeutronCaptureCascade( G4int theZ, G4int theA, G4double NeutronEnergy, G4int InitialLevel,
 80                    std::vector< char >& pType, std::vector< G4double >& pEnergy, std::vector< G4double >& pTime );
 81 
 82     // Initial level for neutron capture. If jspinx2v < 0 it is sampled according to the 2J+1 rule
 83     // l-spin = 0, 1, 2 --> s-wave, p-wave, d-wave ...
 84     G4int SelectInitialLevel( G4int theCompoundZ, G4int theCompoundA, G4double NeutronEnergy, G4int lspin, G4int jspinx2 );
 85     G4int SampleJ( G4int theCompoundZ, G4int theCompoundA, G4int lspin );
 86     G4int GetAllowedJx2values( G4int theCompoundZ, G4int theCompoundA, G4int lspin, G4int* jx2vals );
 87 
 88     const G4NuDEXStatisticalNucleus* GetStatisticalNucleus( G4int za ) { return theStatisticalNucleus[za]; }
 89     G4int Init( G4int theZA, unsigned int seed1 = 0, unsigned int seed2 = 0, unsigned int seed3 = 0 );
 90     void SetBandWidth( G4double bandWidth ) { BandWidth = bandWidth; }
 91     void SetBrOption( G4int brOption ) { BrOption = brOption; }  
 92 
 93     G4NuDEXStatisticalNucleus* theStatisticalNucleus[G4NUDEX_MAXZA];
 94     G4int HasData[G4NUDEX_MAXZA];  // -1:no; 0:don't know; 1:yes
 95     G4String NuDEXLibDirectory;
 96     G4int BrOption;
 97     G4double BandWidth;
 98 
 99     G4int secID;  // creator model ID for the other secondaries produced by this model
100     G4double lowestEnergyLimit;
101     G4double minExcitation; 
102     G4VEvaporationChannel* photonEvaporation;  // Needed when NuDEX is not applicable
103 };
104 
105 #endif
106