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 // $Id$ 00027 // 00028 // ---------------- G4QAtomicElectronScattering header ---------------- 00029 // by Mikhail Kossov, December 2003. 00030 // Header of G4QAtomicElectronScattering class of the CHIPS Simulation Branch in GEANT4 00031 // ------------------------------------------------------------------------------- 00032 // This is a unique CHIPS class for the Nuclear Interactions with Atomic electrons. 00033 // ------------------------------------------------------------------------------- 00034 // **************************************************************************************** 00035 // ********* This HEADER is temporary moved from the photolepton_hadron directory ********* 00036 // ******* DO NOT MAKE ANY CHANGE! With time it'll move back to photolepton...(M.K.) ****** 00037 // **************************************************************************************** 00038 // Short description: CHIPS is re3sponsible for photo- and lepto-nuclear 00039 // reactions. In particular for thr electron-nuclear reactions. At High 00040 // Energies the nucleons (including neutrons) and nuclear fragments can 00041 // interact with atomic electrons (reversed electro-nuclear reaction - 00042 // antilab), while they are missing the nucler-nuclear (ion-ion) reac- 00043 // tions. This nucleo-electron process comes "for-free" in CHIPS, as the 00044 // cross-sections of the interaction is known from the electro-nuclear 00045 // reactions. The only problem is to move the output from the antilab to 00046 // lab system. This is what this process is aiming to do. It can be used 00047 // for the ion transport in Geant4. 00048 // --------------------------------------------------------------------- 00049 00050 #ifndef G4QAtomicElectronScattering_hh 00051 #define G4QAtomicElectronScattering_hh 00052 00053 // GEANT4 Headers 00054 #include "globals.hh" 00055 #include "G4ios.hh" 00056 #include "Randomize.hh" 00057 #include "G4VDiscreteProcess.hh" 00058 #include "G4Track.hh" 00059 #include "G4Step.hh" 00060 #include "G4ParticleTypes.hh" 00061 #include "G4VParticleChange.hh" 00062 #include "G4ParticleDefinition.hh" 00063 #include "G4DynamicParticle.hh" 00064 #include "G4ThreeVector.hh" 00065 #include "G4LorentzVector.hh" 00066 00067 // CHIPS Headers 00068 #include "G4QEnvironment.hh" 00069 #include "G4VQCrossSection.hh" 00070 #include "G4QIsotope.hh" 00071 #include "G4QElectronNuclearCrossSection.hh" 00072 #include "G4QPhotonNuclearCrossSection.hh" 00073 #include "G4QMuonNuclearCrossSection.hh" 00074 #include "G4QTauNuclearCrossSection.hh" 00075 #include "G4QPDGToG4Particle.hh" 00076 00077 class G4QAtomicElectronScattering : public G4VDiscreteProcess 00078 { 00079 public: 00080 00081 // Constructor 00082 G4QAtomicElectronScattering(const G4String& processName ="CHIPSNuclearCollision"); 00083 00084 // Destructor 00085 ~G4QAtomicElectronScattering(); 00086 00087 G4bool IsApplicable(const G4ParticleDefinition& particle); 00088 00089 G4double GetMeanFreePath(const G4Track& aTrack, G4double previousStepSize, 00090 G4ForceCondition* condition); 00091 // It returns the MeanFreePath of the process for the current track : 00092 // (energy, material) 00093 // The previousStepSize and G4ForceCondition* are not used. 00094 // This function overloads a virtual function of the base class. 00095 // It is invoked by the ProcessManager of the Particle. 00096 00097 00098 G4VParticleChange* PostStepDoIt(const G4Track& aTrack, const G4Step& aStep); 00099 // It computes the final state of the process (at end of step), 00100 // returned as a ParticleChange object. 00101 // This function overloads a virtual function of the base class. 00102 // It is invoked by the ProcessManager of the Particle. 00103 00104 00105 G4LorentzVector GetEnegryMomentumConservation(); 00106 00107 G4int GetNumberOfNeutronsInTarget(); 00108 00109 // Static functions 00110 static void SetManual(); 00111 static void SetStandard(); 00112 static void SetParameters(G4double temper=180., G4double ssin2g=.1, G4double etaetap=.3, 00113 G4double fN=0., G4double fD=0., G4double cP=1., G4double mR=1., 00114 G4int npCHIPSWorld=234, G4double solAn=.5, G4bool efFlag=false, 00115 G4double piTh=141.4,G4double mpi2=20000.,G4double dinum=1880.); 00116 00117 private: 00118 00119 // Hide assignment operator as private 00120 G4QAtomicElectronScattering& operator=(const G4QAtomicElectronScattering &right); 00121 00122 // Copy constructor 00123 G4QAtomicElectronScattering(const G4QAtomicElectronScattering&); 00124 00125 // BODY 00126 // Static Parameters 00127 static G4bool manualFlag; // If false then standard parameters are used 00128 static G4int nPartCWorld; // The#of particles for hadronization (limit of A of fragm.) 00129 // -> Parameters of the G4Quasmon class: 00130 static G4double Temperature; // Quasmon Temperature 00131 static G4double SSin2Gluons; // Percent of ssbar sea in a constituen gluon 00132 static G4double EtaEtaprime; // Part of eta-prime in all etas 00133 // -> Parameters of the G4QNucleus class: 00134 static G4double freeNuc; // probability of the quasi-free baryon on surface 00135 static G4double freeDib; // probability of the quasi-free dibaryon on surface 00136 static G4double clustProb; // clusterization probability in dense region 00137 static G4double mediRatio; // relative vacuum hadronization probability 00138 // -> Parameters of the G4QEnvironment class: 00139 static G4bool EnergyFlux; // Flag for Energy Flux use instead of Multy Quasmon 00140 static G4double SolidAngle; // Part of Solid Angle to capture secondaries(@@A-dep) 00141 static G4double PiPrThresh; // Pion Production Threshold for gammas 00142 static G4double M2ShiftVir; // Shift for M2=-Q2=m_pi^2 of the virtual gamma 00143 static G4double DiNuclMass; // Double Nucleon Mass for virtual normalization 00144 // 00145 // Working parameters 00146 G4VQCrossSection* theCS; 00147 G4LorentzVector EnMomConservation; // Residual of Energy/Momentum Cons. 00148 G4int nOfNeutrons; // #of neutrons in the target nucleus 00149 00150 // Modifires for the reaction 00151 G4double Time; // Time shift of the capture reaction 00152 G4double EnergyDeposition; // Energy deposited in the reaction 00153 }; 00154 #endif 00155