00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // Short description: Interface of QGSC to CHIPS (all soft hadrons) 00027 // ---------------------------------------------------------------- 00028 // 00029 #ifndef G4StringChipsParticleLevelInterface_h 00030 #define G4StringChipsParticleLevelInterface_h 00031 00032 #include "G4VIntraNuclearTransportModel.hh" 00033 #include "G4ChiralInvariantPhaseSpace.hh" 00034 00035 // Open this if you wish histogramming of the impact parameter issue 00036 //#define hdebug_SCPLI 00037 00038 class G4StringChipsParticleLevelInterface : public G4VIntraNuclearTransportModel 00039 { 00040 public: 00041 G4StringChipsParticleLevelInterface(); 00042 virtual G4HadFinalState* ApplyYourself(const G4HadProjectile& aTrack, 00043 G4Nucleus& theNucleus); 00044 00045 virtual G4ReactionProductVector* Propagate(G4KineticTrackVector* theSecondaries, 00046 G4V3DNucleus* theNucleus); 00047 private: 00048 00049 G4ChiralInvariantPhaseSpace theModel; 00050 G4double theEnergyLossPerFermi; 00051 00052 G4double theInnerCoreDensityCut; 00053 00054 G4double fractionOfSingleQuasiFreeNucleons; 00055 G4double fractionOfPairedQuasiFreeNucleons; 00056 G4double clusteringCoefficient; 00057 G4double temperature; 00058 G4double halfTheStrangenessOfSee; 00059 G4double etaToEtaPrime; 00060 G4double fusionToExchange; 00061 G4int nop; 00062 00063 #ifdef hdebug_SCPLI 00064 //Static variables for histogramming 00065 static const G4int nbh; 00066 static G4double bhmax; 00067 static G4double bhdb; 00068 static G4double ehmax; 00069 static G4double ehde; 00070 static G4double toth; 00071 static G4int bover; 00072 static G4int eover; 00073 static G4int* bhis; 00074 static G4int* ehis; 00075 public: 00076 //Static functions 00077 static void Reset() 00078 { 00079 bhdb=bhmax/nbh; 00080 ehde=bhmax/nbh; 00081 toth=0.; 00082 bover=0; 00083 eover=0; 00084 for(G4int i=0; i<nbh; i++) 00085 { 00086 bhis[i]=0; 00087 ehis[i]=0; 00088 } 00089 } 00090 static void SetMaxB(G4double mB) {bhmax=mB;} 00091 static void SetMaxE(G4double mE) {ehmax=mE;} 00092 static G4int GetB(G4int i){return bhis[i];} 00093 static G4int GetE(G4int i){return ehis[i];} 00094 static G4double GetTot() {return toth;} 00095 static G4int GetNbn() {return nbh;} 00096 static G4double GetDB() {return bhdb;} 00097 static G4double GetDE() {return ehde;} 00098 static G4int GetBov() {return bover;} 00099 static G4int GetEov() {return eover;} 00100 #endif 00101 }; 00102 #endif
1.4.7