BDSComptonScatteringEngine.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) BDSIM Collaboration, 2001 - 2026.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef BDSCOMPTONSCATTERINGENGINE_H
#define BDSCOMPTONSCATTERINGENGINE_H
 
#include "globals.hh" // geant4 types / globals
#include "G4ThreeVector.hh"
#include "G4DynamicParticle.hh"
 
#include "CLHEP/Units/PhysicalConstants.h"
#include "G4INCLRandom.hh"
#include "G4RotationMatrix.hh"
#include "G4StokesVector.hh"
#include "CLHEP/Units/PhysicalConstants.h"
#include "G4INCLRandom.hh"
 
class BDSComptonScatteringEngine
{
public:
  BDSComptonScatteringEngine();
  ~BDSComptonScatteringEngine();
 
  G4double CrossSection(G4double photonEnergy, G4int partIn);
  void PerformCompton(const G4ThreeVector& boost, G4int partIn);
  void SetParticle(G4int partIDIn);
  
  inline void setIncomingGamma(G4LorentzVector incomingGammaIn) {incomingGamma=incomingGammaIn;}
  inline void setIncomingElectron(G4LorentzVector incomingElectronIn) {incomingElectron=incomingElectronIn;}
  inline G4LorentzVector GetScatteredGamma()     {return scatteredGamma;};
  inline G4LorentzVector GetScatteredElectron()  {return scatteredElectron;};
  G4double KleinNishinaDifferential(G4double theta);
  G4double MCMCTheta();
  G4RotationMatrix* CalculateRotation();
  inline void SetIncomingGammaPolarization(G4StokesVector incomingGammaPolarizationIn) {incomingGammaPolarization=incomingGammaPolarizationIn;}
  inline void SetIncomingElectronPolarization(G4StokesVector incomingElectronPolarizationIn) {incomingElectronPolarization=incomingElectronPolarizationIn;}
  G4double MCPhi(G4double theta, G4double scatteredEnergy);
  G4double PolarizationCrossSectionPhi(G4double theta, G4double phi, G4double scatteredEnergy);
  G4double ScatteredEnergy(G4double theta);
  G4double PolarizationCrossSectionMaxPhi(G4double theta,G4double EP);
private:
  G4LorentzVector incomingGamma;
  G4LorentzVector scatteredGamma;
  G4LorentzVector incomingElectron;
  G4LorentzVector scatteredElectron;
  G4StokesVector incomingGammaPolarization;
  G4StokesVector incomingElectronPolarization;
  G4StokesVector outgoingGammaPolarization;
  G4StokesVector outgoingElectronPolarization;
  G4double particleMass;
  G4double particleRadius;
 
  G4int partID;
};
 
#endif