Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 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 H 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 = FirstCollisionLocalEne 76 localEnergyPiString = "first-collision"; 77 localEnergyPiType = FirstCollisionLocalEne 78 deExcitationString = "none"; 79 deExcitationType = DeExcitationNone; 80 clusterAlgorithmString = "intercomparison" 81 clusterAlgorithmType = IntercomparisonClus 82 clusterMaxMass = 8; 83 backToSpectator = true; 84 useRealMasses = true; 85 impactParameter = -1.; 86 separationEnergyString = "INCL"; 87 separationEnergyType = INCLSeparationEnerg 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 = AntiparticlesCrossSect 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 " << getVersion 125 if(projectileSpecies.theType != Composite) 126 message << "Projectile: " << ParticleTab 127 else 128 message << "Projectile: composite, A=" < 129 message << " energy = " << projectileKine 130 if(targetSpecies.theA>0) 131 message << "Target: A = " << targetSpeci 132 else 133 message << "Target: natural isotopic com 134 message << "Number of requested shots = " 135 return message.str(); 136 } 137 138 } 139