physics/include/G4RDGenerator2BN.hh

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:  G4RDGenerator2BN
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 G4RDGenerator2BN_h
0054 #define G4RDGenerator2BN_h 1
0055
0056 #include "G4ios.hh"
0057 #include "globals.hh"
0058 #include "G4RDVBremAngularDistribution.hh"
0059
0060 class G4RDGenerator2BN : public G4RDVBremAngularDistribution
0061 {
0062
0063 public:
0064
0065   G4RDGenerator2BN(const G4String& name);
0066
0067   ~G4RDGenerator2BN();
0068
0069   G4double PolarAngle(const G4double initial_energy,
0070               const G4double final_energy,
0071               const G4int Z);
0072
0073   void PrintGeneratorInformation() const;
0074
0075 public:
0076
0077   void SetInterpolationThetaIncrement(G4double increment) {dtheta = increment;};
0078   G4double GetInterpolationThetaIncrement() {return dtheta;};
0079
0080   void SetGammaCutValue(G4double cutValue) {kcut = cutValue;};
0081   G4double GetGammaCutValue() {return kcut;};
0082
0083   void ConstructMajorantSurface();
0084
0085 protected:
0086
0087   G4double CalculateFkt(G4double k, G4double theta, G4double A, G4double c) const;
0088   G4double Calculatedsdkdt(G4double kout, G4double theta, G4double Eel) const;
0089   G4double Generate2BN(G4double Ek, G4double k) const;
0090
0091
0092 private:
0093
0094   G4double b;
0095   G4int index_min, index_max;
0096   G4double kmin, Ekmin;
0097   G4double dtheta;
0098   G4double kcut;
0099   static G4double Atab[320];
0100   static G4double ctab[320];
0101
0102   // hide assignment operator
0103   G4RDGenerator2BN & operator=(const  G4RDGenerator2BN &right);
0104   G4RDGenerator2BN(const  G4RDGenerator2BN&);
0105
0106 };
0107
0108 #endif
0109