LinearCompton.h

Go to the documentation of this file.

#ifndef OPAL_LINEAR_COMPTON_HH
#define OPAL_LINEAR_COMPTON_HH
 
#include "OPALTypes.h"
 
#include <cstdint>
#include <random>
 
namespace Physics {
    namespace LinearCompton {
 
        struct SamplingKernel {
            double electronTotalEnergyGeV =
                    0.0;  
            double laserPhotonEnergyGeV =
                    0.0;  
            Vector_t<double, 3> beamDirection = Vector_t<double, 3>(
                    0.0);  
            Vector_t<double, 3> laserDirection = Vector_t<double, 3>(
                    0.0);  
            double incomingPhotonEnergyERFGeV =
                    0.0;  
            double rejectionUpperBoundSolidAngleERF =
                    0.0;  
        };
 
        struct SampledEvent {
            double scatteringCosineERF =
                    0.0;              
            double azimuthERF = 0.0;  
            double scatteredPhotonEnergyERFGeV =
                    0.0;  
            double scatteredPhotonEnergyLabGeV =
                    0.0;  
            Vector_t<double, 3> scatteredPhotonDirectionLab = Vector_t<double, 3>(
                    0.0);  
        };
 
        double photonEnergyFromWavelengthGeV(double wavelength_m);
 
        double electronGamma(double electronTotalEnergyGeV);
 
        double electronBeta(double electronTotalEnergyGeV);
 
        double restFrameIncomingPhotonEnergyGeV(
                double electronTotalEnergyGeV, double laserPhotonEnergyGeV,
                const Vector_t<double, 3>& beamDirection,
                const Vector_t<double, 3>& laserDirection);
 
        double invariantKappa(double incomingPhotonEnergyERFGeV);
 
        double invariantX(double incomingPhotonEnergyERFGeV);
 
        double invariantX(
                double electronTotalEnergyGeV, double laserPhotonEnergyGeV,
                const Vector_t<double, 3>& beamDirection,
                const Vector_t<double, 3>& laserDirection);
 
        double scatteredPhotonEnergyMinERFGeV(double incomingPhotonEnergyERFGeV);
 
        double scatteredPhotonEnergyMaxERFGeV(double incomingPhotonEnergyERFGeV);
 
        double scatteredPhotonEnergyERFGeV(
                double incomingPhotonEnergyERFGeV, double scatteringCosineERF);
 
        double differentialCrossSectionSolidAngleERF(
                double incomingPhotonEnergyERFGeV, double scatteringCosineERF);
 
        Vector_t<double, 3> restFrameIncomingPhotonDirection(
                double electronTotalEnergyGeV, const Vector_t<double, 3>& beamDirection,
                const Vector_t<double, 3>& laserDirection);
 
        Vector_t<double, 3> labPhotonDirection(
                double electronTotalEnergyGeV, const Vector_t<double, 3>& beamDirection,
                const Vector_t<double, 3>& laserDirection, double scatteringCosineERF,
                double azimuthERF);
 
        double labPhotonEnergyGeV(
                double electronTotalEnergyGeV, double laserPhotonEnergyGeV,
                const Vector_t<double, 3>& beamDirection, const Vector_t<double, 3>& laserDirection,
                double scatteringCosineERF, double azimuthERF);
 
        double differentialCrossSectionOmegaERF(
                double incomingPhotonEnergyERFGeV, double scatteredPhotonEnergyERFGeV);
 
        double totalCrossSection(double incomingPhotonEnergyERFGeV);
 
        double thomsonCrossSection();
 
        double restFrameScatteringCosineForLabForwardPhoton(
                double electronTotalEnergyGeV, const Vector_t<double, 3>& beamDirection,
                const Vector_t<double, 3>& laserDirection);
 
        double labForwardPhotonEnergyGeV(
                double electronTotalEnergyGeV, double laserPhotonEnergyGeV,
                const Vector_t<double, 3>& beamDirection,
                const Vector_t<double, 3>& laserDirection);
 
        std::mt19937_64 makeHostRandomEngine(std::uint64_t streamIndex = 0);
 
        SamplingKernel makeSamplingKernel(
                double electronTotalEnergyGeV, double laserPhotonEnergyGeV,
                const Vector_t<double, 3>& beamDirection,
                const Vector_t<double, 3>& laserDirection);
 
        SampledEvent sampleEvent(const SamplingKernel& kernel, std::mt19937_64& engine);
 
    }  // namespace LinearCompton
}  // namespace Physics
 
#endif  // OPAL_LINEAR_COMPTON_HH