Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/inclxx/utils/src/G4INCLConfig.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 // INCL++ intra-nuclear cascade model
 27 // Alain Boudard, CEA-Saclay, France
 28 // Joseph Cugnon, University of Liege, Belgium
 29 // Jean-Christophe David, CEA-Saclay, France
 30 // Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
 31 // Sylvie Leray, CEA-Saclay, France
 32 // Davide Mancusi, CEA-Saclay, France
 33 //
 34 #define INCLXX_IN_GEANT4_MODE 1
 35 
 36 #include "globals.hh"
 37 
 38 #include "G4INCLParticleType.hh"
 39 #include "G4INCLConfig.hh"
 40 #include "G4INCLParticleSpecies.hh"
 41 #include "G4INCLParticleTable.hh"
 42 
 43 namespace G4INCL {
 44   
 45   Config::Config() {
 46     init();
 47   }
 48 
 49   Config::~Config() {}
 50 
 51   void Config::init() {
 52     verbosity = 1;
 53     logFileName = "-";
 54     inputFileName = "";
 55     title = "INCL default run title";
 56     nShots = 1000;
 57     naturalTarget = false;
 58     projectileString = "proton";
 59     projectileSpecies = G4INCL::Proton;
 60     projectileKineticEnergy = 1000.0;
 61     verboseEvent = -1;
 62     randomSeeds = "";
 63     randomSeedVector.push_back(666);
 64     randomSeedVector.push_back(777);
 65     randomSeedVector.push_back(1234);
 66     pauliString = "strict-statistical";
 67     pauliType = StrictStatisticalPauli;
 68     CDPP = true;
 69     coulombString = "non-relativistic";
 70     coulombType = NonRelativisticCoulomb;
 71     potentialString = "isospin-energy";
 72     potentialType = IsospinEnergyPotential;
 73     pionPotential = true;
 74     localEnergyBBString = "first-collision";
 75     localEnergyBBType = FirstCollisionLocalEnergy;
 76     localEnergyPiString = "first-collision";
 77     localEnergyPiType = FirstCollisionLocalEnergy;
 78     deExcitationString = "none";
 79     deExcitationType = DeExcitationNone;
 80     clusterAlgorithmString = "intercomparison";
 81     clusterAlgorithmType = IntercomparisonClusterAlgorithm;
 82     clusterMaxMass = 8;
 83     backToSpectator = true;
 84     useRealMasses = true;
 85     impactParameter = -1.;
 86     separationEnergyString = "INCL";
 87     separationEnergyType = INCLSeparationEnergy;
 88     fermiMomentumString = "constant";
 89     fermiMomentumType = ConstantFermiMomentum;
 90     fermiMomentum = -1.;
 91     cutNN = 1910.;
 92 #ifdef INCL_DEEXCITATION_FERMI_BREAKUP
 93     maxMassFermiBreakUp = 16;
 94     maxChargeFermiBreakUp = 8;
 95 #endif
 96     rpCorrelationCoefficient = 0.98;
 97     rpCorrelationCoefficientProton = 0.5;
 98     rpCorrelationCoefficientNeutron = 0.73;
 99     neutronSkin = 0.;
100     neutronHalo = 0.;
101     refraction=false;
102     phaseSpaceGenerator = "Raubold-Lynch";
103     phaseSpaceGeneratorType = RauboldLynchType;
104     cascadeAction = "default";
105     cascadeActionType = DefaultActionType;
106     randomNumberGenerator = "Ranecu";
107     rngType = RanecuType;
108     autosaveFrequency = 10000;
109     maxNumberMultipions = -1;
110     crossSectionsString = "antiparticles";
111     crossSectionsType = AntiparticlesCrossSections;
112     hadronizationTime = 0.;
113 #ifdef INCL_ROOT_USE
114     conciseROOTTree = false;
115 #endif
116     inverseKinematics = false;
117     decayTimeThreshold = 1.e-20;
118     bias = 1.;
119     atrestThreshold = 200.;
120   }
121 
122   std::string Config::summary() {
123     std::stringstream message;
124     message << "INCL++ version " << getVersionString() << '\n';
125     if(projectileSpecies.theType != Composite)
126       message << "Projectile: " << ParticleTable::getName(projectileSpecies) << '\n';
127     else
128       message << "Projectile: composite, A=" << projectileSpecies.theA << ", Z=" << projectileSpecies.theZ << ", S=" << projectileSpecies.theS << '\n';
129     message << "  energy = " << projectileKineticEnergy << '\n';
130     if(targetSpecies.theA>0)
131       message << "Target: A = " << targetSpecies.theA << " Z = " << targetSpecies.theZ << " S = " << targetSpecies.theS << '\n';
132     else
133       message << "Target: natural isotopic composition, Z = " << targetSpecies.theZ << '\n';
134     message << "Number of requested shots = " << nShots << '\n';
135     return message.str();
136   }
137 
138 }
139