ERDecayEXP1811.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 ERDecayEXP1811_H
#define ERDecayEXP1811_H

#include <vector>
#include <fstream>

#include "TGraph.h"
#include "TF1.h"
#include "TRandom3.h"
#include "TRandom2.h"
#include "TGenPhaseSpace.h"
#include "TLorentzVector.h"

#include "FairIon.h"

#include "ERDecay.h"

class ERDecayEXP1811 : public ERDecay {

public:
  ERDecayEXP1811();
  ~ERDecayEXP1811();

  /*Modifiers*/
  void SetMinStep(Double_t minStep) {fMinStep = minStep;}
  void SetTargetThickness(Double_t targetThickness) {fTargetThickness = targetThickness;}
  void SetH7Mass(Double_t mass) {f7HMass = mass; fIs7HUserMassSet = true;}
  void SetH7Exitation(Double_t excMean, Double_t fwhm, Double_t distibWeight);
  void SetDecayFile(const TString& filePath, Double_t excitationEnergyInFile /*[GeV]*/){ fDecayFilePath = filePath; }

  void SetAngularDistribution(TString ADfile);   

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

  void BeginEvent();
  void FinishEvent();

private:
  void ReactionPhaseGenerator(Double_t Ecm, Double_t h7Mass);

  Bool_t DecayPhaseGenerator(Double_t excitation);

  std::vector<TLorentzVector> ReadDecayEvent();

private:
  TRandom3       *fRnd;
  TRandom3       *fRnd2;
  
  TParticlePDG   *f8He;
  TParticlePDG   *f2H;
  TParticlePDG   *f3He;
  TParticlePDG   *f7H;
  TParticlePDG   *f3H;
  TParticlePDG   *fn;

  TLorentzVector *fLv3He; 
  TLorentzVector *fLv7H; 
  TLorentzVector *fLv3H; 
  TLorentzVector *fLvn1; 
  TLorentzVector *fLvn2; 
  TLorentzVector *fLvn3; 
  TLorentzVector *fLvn4; 

  FairIon        *fIon3He;
  FairIon        *fUnstableIon7H;

  TGenPhaseSpace  *fReactionPhaseSpace;
  TGenPhaseSpace  *fDecayPhaseSpace;
  Double_t         fTargetReactZ;
  Double_t         fMinStep;
  Double_t         fTargetThickness;
  Bool_t           fDecayFinish;

  std::vector<Double_t> f7HExcitationMean;
  std::vector<Double_t> f7HExcitationSigma; 
  std::vector<Double_t> f7HExcitationWeight;

  Double_t        f7HMass;
  Bool_t          fIs7HUserMassSet;
  Bool_t          fIs7HExcitationSet;

  TString         fDecayFilePath;
  Double_t        fDecayFileExcitation = 1. /*[GeV]*/;
  Bool_t          fDecayFileFinished;
  Int_t           fDecayFileCurrentEvent;
  std::ifstream   fDecayFile;

  TGraph *fADInput = nullptr;    
  TF1    *fADFunction = nullptr; 
  Double_t fThetaMin = 0.;
  Double_t fThetaMax = 0.;

  //ADEvaluate function is necessary for TF1 constructor
  Double_t ADEvaluate(Double_t *x, Double_t *p);

  ClassDef(ERDecayEXP1811,1)
};

#endif