ERDecay.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 ERDECAYER_H
#define ERDECAYER_H

#include "TRandom3.h"
#include "TString.h"
#include "TGenPhaseSpace.h"
#include "TDatabasePDG.h"               //for TDatabasePDG
#include "TVector3.h"

class ERDecay{
public:
  ERDecay(TString name);
  ~ERDecay();
  /*Accessors*/
  TString GetName() const {return fName;}
  
  /*Modifiers*/
  
  void SetInteractionVolumeName(TString name) {fInteractionVolumeName = name;}
  
  void SetNuclearInteractionLength(Double_t lambda) {fNuclearInteractionLength = lambda;}
  
  void SetMaxPathLength(Double_t pathLength);

  TString GetInteractionVolumeName() {return fInteractionVolumeName;}

  virtual Bool_t Stepping() = 0;
  virtual void   BeginEvent();
  virtual void   FinishEvent();
  virtual Bool_t Init() = 0;

protected:
  Bool_t FindInteractionPoint();
  void   CalculateTargetParameters();
  
  Bool_t   fIsInterationPointFound;
  Double_t fDistanceToInteractPoint;
  Double_t fDistanceFromEntrance;

private:
  TString  fInteractionVolumeName;
  Double_t fNuclearInteractionLength;
  Double_t fInteractionProbability;
  // Double_t fTargetBoundBoxDiagonal;
  Double_t fNormalizingProbability;
  TRandom3 *fRnd1;
  TRandom3 *fRnd2;
  // ----Old variables and methods ------------------------------------------------
protected:
  TString fName;

  TRandom3 *fRnd;

  Bool_t fUniform;
  Float_t fUniformA;
  Float_t fUniformB;

  Bool_t fExponential;
  Float_t fExponentialStart;
  Float_t fExponentialTau;

  TString fVolumeName;
  Double_t fStep;

  Double_t fTargetMass;

  TString fInputIonName;
  TParticlePDG*   fInputIonPDG;

  Double_t fDecayPosZ;
  TLorentzVector fDecayPos;
  Bool_t fDecayFinish;

  void AddParticleToStack(Int_t pdg, TLorentzVector pos, TLorentzVector state);

public:

  void SetDecayPosZ(Double_t pos) {fDecayPosZ = pos;}
  void SetInputIon(Int_t A, Int_t Z, Int_t Q);
  void SetUniformPos(Double_t a, Double_t b);
  void SetExponentialPos(Double_t start, Double_t tau);
  void SetDecayVolume(TString name){fVolumeName = name;}
  void SetStep(Double_t step){fStep = step;}
  void SetTargetMass(Double_t targetMass){fTargetMass = targetMass;}
  TString GetName() {return fName;}
  // ------------------------------------------------------------------------------

  ClassDef(ERDecay,1)
};

#endif