EXP1803

/********************************************************************************

 *              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 ERDecayEXP1803_H

#define ERDecayEXP1803_H

#include "TGraph.h"

#include "TF1.h"

#include "TRandom3.h"

#include "TGenPhaseSpace.h"

#include "FairIon.h"

#include "ERDecay.h" // mother class

class ERDecayEXP1803 : public ERDecay {

public:

  ERDecayEXP1803();

  ~ERDecayEXP1803();

  /*Modifiers*/

  void SetMinStep(Double_t minStep) {fMinStep = minStep;}

  void SetTargetVolumeName(TString volumeName) {fVolumeName = volumeName;}

  void SetTargetThickness(Double_t targetThickness) {fTargetThickness = targetThickness;}

  void SetH5Mass(Double_t mass) {f5HMass = mass; fIs5HUserMassSet = true;}

  void SetH5Exitation(Double_t excMean, Double_t fwhm, Double_t distibWeight);

  void phasegen2(Double_t Ecm, Double_t h5Mass);

  void ReadADInput(TString ADfile);

public:

  Bool_t Init();

  Bool_t Stepping();

  void BeginEvent();

  void FinishEvent();

private:

  TRandom3       *fRnd;

  TParticlePDG   *f6He;

  TParticlePDG   *f2H;

  TParticlePDG   *f3He;

  TParticlePDG   *f5H;

  TParticlePDG   *f3H;

  TParticlePDG   *fn;

  TLorentzVector *flv3He; //!

  TLorentzVector *flv5H; //!

  FairIon        *fIon3He;

  FairIon        *fUnstableIon5H;

  TGenPhaseSpace  *fReactionPhaseSpace;

  TGenPhaseSpace  *fDecayPhaseSpace;

  TString          fVolumeName;

  Double_t         fTargetReactZ;

  Double_t         fMinStep;

  Double_t         fTargetThickness;

  Bool_t           fDecayFinish;

  std::vector<Double_t> f5HExcitationMean;

  std::vector<Double_t> f5HExcitationSigma; 

  std::vector<Double_t> f5HExcitationWeight;

  Double_t        f5HMass;

  Bool_t          fIs5HUserMassSet;

  Bool_t          fIs5HExcitationSet;

  TString fADFile;     //!   distribution is taken from file containing angular

  TGraph *fADInput;    //!   distribution (AD) graph containing AD input

  TF1    *fADFunction; //!   function describing AD of binary reaction

  Double_t ADEvaluate(Double_t *x, Double_t *p);

  ClassDef(ERDecayEXP1803,1)

};

#endif