Back to home page

EIC code displayed by LXR

 
 

    


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

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:     G4PAIPhotData
0033 //
0034 // Author:        V. Grichine based on  G4PAIModelData code for MT
0035 //
0036 // Creation date: 07.10.2013
0037 //
0038 // Modifications:
0039 //
0040 //
0041 // Class Description:
0042 //
0043 // Implementation of PAIPhotModel internal data class.
0044 // This class is extracted from G4PAIPhot in order to provide sharing
0045 // of these data between threads.
0046 //
0047 // Internal data tables are computed for proton. 
0048 //
0049 // -------------------------------------------------------------------
0050 //
0051 
0052 #ifndef G4PAIPhotData_h
0053 #define G4PAIPhotData_h 1
0054 
0055 #include <vector>
0056 #include "globals.hh"
0057 #include "G4PAIxSection.hh"
0058 #include "G4SandiaTable.hh"
0059 
0060 class G4PhysicsLogVector;
0061 class G4PhysicsTable;
0062 class G4MaterialCutsCouple;
0063 class G4PAIPhotModel;
0064 
0065 class G4PAIPhotData 
0066 {
0067 
0068 public:
0069 
0070   explicit G4PAIPhotData(G4double tmin, G4double tmax, G4int verbose);
0071 
0072   ~G4PAIPhotData();
0073 
0074   void Initialise(const G4MaterialCutsCouple*, G4double cut, G4PAIPhotModel*);
0075 
0076   G4double DEDXPerVolume(G4int coupleIndex, G4double scaledTkin,
0077              G4double cut) const;
0078 
0079   G4double CrossSectionPerVolume(G4int coupleIndex, G4double scaledTkin,
0080                  G4double tcut, G4double tmax) const;
0081 
0082   G4double GetPlasmonRatio( G4int coupleIndex, G4double scaledTkin ) const;
0083 
0084   G4double SampleAlongStepTransfer(G4int coupleIndex, G4double kinEnergy,
0085                    G4double scaledTkin,
0086                    G4double stepFactor) const;
0087   G4double SampleAlongStepPhotonTransfer(G4int coupleIndex, G4double kinEnergy,
0088                    G4double scaledTkin,
0089                    G4double stepFactor) const;
0090   G4double SampleAlongStepPlasmonTransfer(G4int coupleIndex, G4double kinEnergy,
0091                    G4double scaledTkin,
0092                    G4double stepFactor) const;
0093 
0094   G4double SamplePostStepTransfer(G4int coupleIndex, 
0095                   G4double scaledTkin) const;
0096   G4double SamplePostStepPhotonTransfer(G4int coupleIndex, 
0097                   G4double scaledTkin) const;
0098   G4double SamplePostStepPlasmonTransfer(G4int coupleIndex, 
0099                   G4double scaledTkin) const;
0100 
0101   // hide assignment operator 
0102   G4PAIPhotData & operator=(const  G4PAIPhotData &right) = delete;
0103   G4PAIPhotData(const  G4PAIPhotData&) = delete;
0104 
0105 private:
0106 
0107   G4double GetEnergyTransfer(G4int coupleIndex, size_t iPlace, 
0108                  G4double position) const;
0109   G4double GetEnergyPhotonTransfer(G4int coupleIndex, size_t iPlace, 
0110                  G4double position) const;
0111   G4double GetEnergyPlasmonTransfer(G4int coupleIndex, size_t iPlace, 
0112                  G4double position) const;
0113 
0114   G4int                fTotBin;
0115   G4double             fLowestKineticEnergy;
0116   G4double             fHighestKineticEnergy;
0117 
0118   G4PhysicsLogVector*  fParticleEnergyVector;
0119 
0120   G4PAIxSection        fPAIxSection;
0121   G4SandiaTable        fSandia;
0122 
0123   std::vector<G4PhysicsTable*>      fPAIxscBank;
0124   std::vector<G4PhysicsTable*>      fPAIphotonBank;
0125   std::vector<G4PhysicsTable*>      fPAIplasmonBank;
0126 
0127   std::vector<G4PhysicsTable*>      fPAIdEdxBank;
0128   std::vector<G4PhysicsLogVector*>  fdEdxTable;
0129 
0130   std::vector<G4PhysicsLogVector*>  fdNdxCutTable;
0131   std::vector<G4PhysicsLogVector*>  fdNdxCutPhotonTable;
0132   std::vector<G4PhysicsLogVector*>  fdNdxCutPlasmonTable;
0133 
0134   std::vector<G4PhysicsLogVector*>  fdEdxCutTable;
0135 };
0136 
0137 #endif
0138 
0139 
0140 
0141 
0142 
0143 
0144