ERElasticScattering.h

/********************************************************************************
 *              Copyright (C) Joint Institute for Nuclear Research              *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *         GNU Lesser General Public Licence version 3 (LGPL) version 3,        *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/

#ifndef ERElasticScattering_H
#define ERElasticScattering_H

#include <ERDecay.h> // mother class

#include <TString.h>
#include <TLorentzVector.h>

enum ERInteractionParticipant{ kPROJECTILE, kTARGET, kEJECTILE}; // TODO It has to be in the header of the future class "Interaction.h"

class TF1;
class TParticlePDG;

class ERElasticScattering : public ERDecay
{
public:
  ERElasticScattering(TString name);

  virtual ~ERElasticScattering();

  void SetTargetIon(Int_t A, Int_t Z, Int_t Q);

  void SetThetaCDF(TString fileName) { fThetaFileName = fileName; }

  void SetThetaRange(Double_t th1, Double_t th2, ERInteractionParticipant DetIonType = kPROJECTILE);

  void SetLabThetaRange(Double_t thetaCenter, Double_t dTheta, ERInteractionParticipant DetIonType = kPROJECTILE, 
                        Bool_t relMod = kTRUE, Double_t BeamAvE = 0.);

  void SetPhiRange(Double_t phi1, Double_t phi2) { fPhiMin = phi1; fPhiMax = phi2; }

  void SetProjectileIonMass(Double_t mass) { fProjectileIonMass = mass; }

  void SetTargetIonMass(Double_t mass) { fTargetIonMass = mass; }


  Double_t GetProjectileIonMass() const { return fProjectileIonMass; }

  Double_t GetTargetIonMass() const { return fTargetIonMass; }

  Double_t GetdPhi() const { return (fPhiMax - fPhiMin); }

  Double_t GetdThetaCDF() const { return fCDFmax - fCDFmin; }

  Int_t GetInteractNumInTarget() const { return fInteractNumInTarget; }

  Double_t GetThetaCMMean() const;

public:
  Bool_t Init();
  Bool_t Stepping();

private:
  void ThetaRangesLab2CM();

  void ThetaRangesLab2CMRelativistic();

  Bool_t ThetaCDFRead();

  Double_t ThetaGen();

  Bool_t DefineOfIonsMasses();

private:
  ERInteractionParticipant  fDetectionIonType=kPROJECTILE;  

  TString         fThetaFileName="";           
  TString         fTargetIonName="";           

  TParticlePDG*   fTargetIonPDG=NULL;          

  TF1*            fThetaCDF=NULL;              
  TF1*            fThetaInvCDF=NULL;           

  Double_t        fThetaMinCM=0.;              
  Double_t        fThetaMaxCM=180.;            
  Double_t        fThetaRangeCenter=0.;        
  Double_t        fThetaRangedTheta=0.;        
  Double_t        fPhiMin=0.;                  
  Double_t        fPhiMax=360.;                
  Double_t        fCDFmin=0.;                  
  Double_t        fCDFmax=1.;                  
  Double_t        fProjectileIonMass=-1.;      
  Double_t        fTargetIonMass=0.;           
  Double_t        fBeamAverageEnergy=0.;       
  Bool_t          fRelativisticMode=kTRUE;     

  Double_t        fThetaCMSum=0.;              
  Int_t           fInteractNumInTarget=0;      
 
  ClassDef(ERElasticScattering, 1);
};

#endif // ERElasticScattering_H