Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2025-01-18 09:58:20

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 file
0030 //
0031 //
0032 // File name:  G4Generator2BN
0033 //
0034 // Author:     Andreia Trindade (andreia@lip.pt)
0035 //             Pedro Rodrigues  (psilva@lip.pt)
0036 //             Luis Peralta      (luis@lip.pt)
0037 // 
0038 // Creation date: 2 June 2003
0039 //
0040 // Modifications: 
0041 // 02 Jun 2003                           First implementation acording with new design
0042 // 07 Nov 2003                           Lockup tables for fast initialization
0043 //               
0044 //
0045 // Class Description: 
0046 //
0047 // Concrete class for Bremsstrahlung Angular Distribution Generation - 2BN Distribution
0048 // Further documentation available from http://www.ge.infn.it/geant4/lowE
0049 
0050 // -------------------------------------------------------------------
0051 //
0052 
0053 #ifndef G4Generator2BN_h
0054 #define G4Generator2BN_h 1
0055 
0056 #include "G4ios.hh"
0057 #include "globals.hh"
0058 #include "G4VEmAngularDistribution.hh"
0059 #include "G4Generator2BS.hh"
0060 
0061 class G4Generator2BN : public G4VEmAngularDistribution
0062 {
0063 
0064 public:
0065   explicit G4Generator2BN(const G4String& name = "");
0066 
0067   virtual ~G4Generator2BN(){;};
0068 
0069   G4ThreeVector& SampleDirection(const G4DynamicParticle* dp,
0070                  G4double out_energy,
0071                  G4int Z,
0072                  const G4Material* mat = nullptr) override;
0073   
0074   void PrintGeneratorInformation() const override;
0075 
0076   void SetInterpolationThetaIncrement(G4double increment) {dtheta = increment;};
0077   G4double GetInterpolationThetaIncrement() {return dtheta;};
0078 
0079   void SetGammaCutValue(G4double cutValue) {kcut = cutValue;};
0080   G4double GetGammaCutValue() {return kcut;};
0081 
0082   void ConstructMajorantSurface();
0083 
0084   G4Generator2BN & operator=(const  G4Generator2BN &right) = delete;
0085   G4Generator2BN(const  G4Generator2BN&) = delete;
0086 
0087 protected:
0088 
0089   G4double CalculateFkt(G4double k, G4double theta, G4double A, G4double c) const;
0090   G4double Calculatedsdkdt(G4double kout, G4double theta, G4double Eel) const;
0091 
0092 private:
0093   G4Generator2BS fGenerator2BS; 
0094   static G4double Atab[320];
0095   static G4double ctab[320];
0096   G4double b;
0097  
0098   G4double kmin, Ekmin;
0099   G4double dtheta;
0100   G4double kcut;
0101  
0102   G4int index_min, index_max;
0103   G4int nwarn;
0104 
0105 };
0106 
0107 #endif
0108