Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/hadronic/Hadr02/src/IonCRMCPhysics.cc

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 /// \file hadronic/Hadr02/src/IonCRMCPhysics.cc
 27 /// \brief Implementation of the IonCRMCPhysics class
 28 //
 29 //
 30 //---------------------------------------------------------------------------
 31 //
 32 // Class:     IonCRMCPhysics
 33 //
 34 // Author:    2018 Alberto Ribon
 35 //
 36 // Modified:
 37 // -  18-May-2021 Alberto Ribon : Used the latest Geant4-CRMC interface.
 38 //
 39 //---------------------------------------------------------------------------
 40 //
 41 #ifdef G4_USE_CRMC
 42 
 43 #  include "IonCRMCPhysics.hh"
 44 
 45 #  include "HadronicInelasticModelCRMC.hh"
 46 
 47 #  include "G4Alpha.hh"
 48 #  include "G4BinaryLightIonReaction.hh"
 49 #  include "G4BuilderType.hh"
 50 #  include "G4ComponentGGNuclNuclXsc.hh"
 51 #  include "G4CrossSectionInelastic.hh"
 52 #  include "G4Deuteron.hh"
 53 #  include "G4ExcitationHandler.hh"
 54 #  include "G4FTFBuilder.hh"
 55 #  include "G4GenericIon.hh"
 56 #  include "G4HadronInelasticProcess.hh"
 57 #  include "G4HadronicInteraction.hh"
 58 #  include "G4HadronicInteractionRegistry.hh"
 59 #  include "G4HadronicParameters.hh"
 60 #  include "G4He3.hh"
 61 #  include "G4IonConstructor.hh"
 62 #  include "G4IonPhysics.hh"
 63 #  include "G4ParticleDefinition.hh"
 64 #  include "G4PreCompoundModel.hh"
 65 #  include "G4ProcessManager.hh"
 66 #  include "G4SystemOfUnits.hh"
 67 #  include "G4Triton.hh"
 68 
 69 using namespace std;
 70 
 71 // factory
 72 #  include "G4PhysicsConstructorFactory.hh"
 73 //
 74 G4_DECLARE_PHYSCONSTR_FACTORY(IonCRMCPhysics);
 75 
 76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 77 
 78 const std::array<std::string, 13> IonCRMCPhysics::fModelNames = {
 79   "EPOS-LHC", "EPOS-1.99", "QGSJET-01",   "",           "", "", "SIBYLL-2.3", "QGSJETII-04", "",
 80   "",         "",          "QGSJETII-03", "DPMJET-3.06"};
 81 
 82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 83 
 84 IonCRMCPhysics::IonCRMCPhysics(G4int ver) : G4VPhysicsConstructor("ionInelasticCRMC")
 85 {
 86   fModel = 0;  //***LOOKHERE*** CRMC model: 0:EPOS-LHC, 1:EPOS-1.99, 2:QGSJET:01, 6:SIBYLL-2.3,
 87                //                           7:QGSJETII-04, 11:QGSJETII-03, 12:DPMJET-3.06
 88   fVerbose = ver;
 89   if (fVerbose > 1) G4cout << "### IonCRMCPhysics" << G4endl;
 90   SetPhysicsType(bIons);
 91 }
 92 
 93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 94 
 95 IonCRMCPhysics::~IonCRMCPhysics() {}
 96 
 97 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 98 
 99 void IonCRMCPhysics::ConstructParticle()
100 {
101   // Construct ions
102   G4IonConstructor pConstructor;
103   pConstructor.ConstructParticle();
104 }
105 
106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
107 
108 void IonCRMCPhysics::ConstructProcess()
109 {
110   fModel = 0;  //***LOOKHERE*** 0:EPOS-LHC, 1:EPOS-1.99, 2:QGSJET:01, 6:SIBYLL-2.3,
111                //               7:QGSJETII-04, 11:QGSJETII-03, 12:DPMJET-3.06
112   const G4double minCRMC =
113     100.0
114     * GeV;  //***LOOKHERE*** CRMC model is applied only above this projectile lab energy per nucleon
115   const G4double maxFTFP =
116     110.0
117     * GeV;  //***LOOKHERE*** FTFP model is applied only below this projectile lab energy per nucleon
118   G4HadronicInteraction* p = G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
119   G4PreCompoundModel* thePreCompound = static_cast<G4PreCompoundModel*>(p);
120   if (!thePreCompound) thePreCompound = new G4PreCompoundModel;
121   // Binary Cascade
122   G4HadronicInteraction* theIonBC = new G4BinaryLightIonReaction(thePreCompound);
123   theIonBC->SetMinEnergy(0.0);
124   theIonBC->SetMaxEnergy(G4HadronicParameters::Instance()->GetMaxEnergyTransitionFTF_Cascade());
125   // FTFP
126   G4FTFBuilder theBuilder("FTFP", thePreCompound);
127   G4HadronicInteraction* theFTFP = theBuilder.GetModel();
128   theFTFP->SetMinEnergy(G4HadronicParameters::Instance()->GetMinEnergyTransitionFTF_Cascade());
129   theFTFP->SetMaxEnergy(maxFTFP);
130   // CRMC
131   G4HadronicInteraction* theCRMC = new HadronicInelasticModelCRMC(fModel, fModelNames[fModel]);
132   theCRMC->SetMinEnergy(minCRMC);
133   theCRMC->SetMaxEnergy(G4HadronicParameters::Instance()->GetMaxEnergy());
134   // Cross section
135   G4CrossSectionInelastic* theXS = new G4CrossSectionInelastic(new G4ComponentGGNuclNuclXsc);
136   // Processes
137   AddProcess("dInelastic", G4Deuteron::Deuteron(), theIonBC, theFTFP, theCRMC, theXS);
138   AddProcess("tInelastic", G4Triton::Triton(), theIonBC, theFTFP, theCRMC, theXS);
139   AddProcess("He3Inelastic", G4He3::He3(), theIonBC, theFTFP, theCRMC, theXS);
140   AddProcess("alphaInelastic", G4Alpha::Alpha(), theIonBC, theFTFP, theCRMC, theXS);
141   AddProcess("ionInelastic", G4GenericIon::GenericIon(), theIonBC, theFTFP, theCRMC, theXS);
142   if (fVerbose > 1) G4cout << "IonCRMCPhysics::ConstructProcess done! " << G4endl;
143 }
144 
145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
146 
147 void IonCRMCPhysics::AddProcess(const G4String& name, G4ParticleDefinition* part,
148                                 G4HadronicInteraction* theIonBC, G4HadronicInteraction* theFTFP,
149                                 G4HadronicInteraction* theCRMC, G4VCrossSectionDataSet* xs)
150 {
151   G4HadronInelasticProcess* hadi = new G4HadronInelasticProcess(name, part);
152   G4ProcessManager* pManager = part->GetProcessManager();
153   pManager->AddDiscreteProcess(hadi);
154   if (xs) hadi->AddDataSet(xs);
155   if (theIonBC) hadi->RegisterMe(theIonBC);
156   if (theFTFP) hadi->RegisterMe(theFTFP);
157   if (theCRMC) hadi->RegisterMe(theCRMC);
158 }
159 
160 #endif  // G4_USE_CRMC
161