Back to home page

EIC code displayed by LXR

 
 

    


Warning, file /include/Geant4/G4eDPWACoulombScatteringModel.hh was not indexed or was modified since last indexation (in which case cross-reference links may be missing, inaccurate or erroneous).

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 // -------------------------------------------------------------------
0028 //
0029 // GEANT4 Class header file
0030 //
0031 //
0032 // File name:     G4eDPWACoulombScatteringModel
0033 //
0034 // Author:        Mihaly Novak
0035 //
0036 // Creation date: 02.07.2020
0037 //
0038 // Modifications:
0039 //
0040 // Class Description:
0041 //
0042 // e-/e+ Coulomb scattering model based on numerical Differential Cross Sections
0043 // (DCS) obtained by Dirac Partial Wave Analysis (DPWA) and supplied by the
0044 // G4eDPWAElasticDCS class.
0045 // The model contains the possibility to incorporate the effects of angular
0046 // deflections of sub-threshold ionisation intercations when it's described by
0047 // the condensed history model. Note, this must be inactivated (by setting the
0048 // `isscpcor` input argument of the CTR to false) when ionisation is described
0049 // with a classical, event by event based simulation model instead of usign the
0050 // condensed history approach (otherwise, the corresponding angular defelctions
0051 // will be "double counted").
0052 //
0053 // -------------------------------------------------------------------
0054 
0055 
0056 
0057 #ifndef G4eDPWACoulombScatteringModel_h
0058 #define G4eDPWACoulombScatteringModel_h 1
0059 
0060 #include "G4VEmModel.hh"
0061 #include "globals.hh"
0062 
0063 class G4eDPWAElasticDCS;
0064 class G4ParticleChangeForGamma;
0065 class G4ParticleDefinition;
0066 class G4DataVector;
0067 
0068 class G4eDPWACoulombScatteringModel : public G4VEmModel {
0069 
0070 public:
0071 
0072   /**
0073    * Constructor.
0074    *
0075    * @param[in] ismixed  Indicates if the model is for mixed or for pure single
0076    *                     Coulomb scattering. Different type of tables are pre-
0077    *                     pared for sampling polar angle of Coulomb scattering
0078    *                     for mixed and for pure single scattering models: cosine
0079    *                     of the polar scattering angle can be sampled in a
0080    *                     restriced inteval (see fMuMin parameter).
0081    * @param[in] isscpcor Indicates if scattering power correction should be used.
0082    *                     Note, scattering power correction accounts the effects
0083    *                     angular deflections due to sub-threshold ionisations
0084    *                     when ionisation is described by using condensed history
0085    *                     model (should be active only in this case).
0086    */
0087   G4eDPWACoulombScatteringModel(G4bool ismixed=false, G4bool isscpcor=true);
0088 
0089   ~G4eDPWACoulombScatteringModel() override;
0090 
0091   //
0092   // Interface methods:
0093 
0094   void     Initialise(const G4ParticleDefinition*, const G4DataVector&) override;
0095 
0096   void     InitialiseLocal(const G4ParticleDefinition*, G4VEmModel*) override;
0097 
0098   G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*, G4double ekin,
0099                                       G4double Z, G4double A, G4double prodcut,
0100                                       G4double emax) override;
0101 
0102   void     SampleSecondaries(std::vector<G4DynamicParticle*>*,
0103                              const G4MaterialCutsCouple*,
0104                              const G4DynamicParticle*,
0105                              G4double tmin,
0106                              G4double maxEnergy) override;
0107 
0108   G4double MinPrimaryEnergy(const G4Material*, const G4ParticleDefinition*,
0109                             G4double) override;
0110 
0111   void     SetTheDCS(G4eDPWAElasticDCS* theDCS) { fTheDCS = theDCS; }
0112 
0113   G4eDPWAElasticDCS* GetTheDCS() { return fTheDCS; }
0114 
0115 
0116 private:
0117 
0118   // Indicates if the model is mixed: MSC for soft (theta<theta_c), Singe
0119   // Scattering(SS) for hard scatterings(theta>theta_c). SS otherwise.
0120   // Note, that while the model provides restricted (elastic and transport)
0121   // cross sections, it's responsible to handle, i.e. provide final state,
0122   // only for the Singe Scattering part in case of a mixed model.
0123   G4bool                     fIsMixedModel;
0124   // indicates if scattering power correction should be applied: correction due
0125   // to deflection in case of sub-threshold, inelastic interactions -> only in
0126   // case of condensed history simulation of inonisation!
0127   G4bool                     fIsScpCorrection;
0128   // mu(theta)=0.5[1-cos(theta)]: the model porvides final states \in [fMuMin,1]
0129   // NOTE: `theta` for this limit can be set in `G4EmParameters::SetMscThetaLimit`
0130   G4double                   fMuMin;
0131   // the object that provides cross sections and polar angle of scattering
0132   G4eDPWAElasticDCS*         fTheDCS;
0133   // particle change
0134   G4ParticleChangeForGamma*  fParticleChange;
0135 
0136 };
0137 
0138 #endif