KVGenPhaseSpace.h

 
#ifndef __KVGENPHASESPACE_H
#define __KVGENPHASESPACE_H
 
#include "KVBase.h"
#include "KVNucleus.h"
#include "KVEvent.h"
 
 
 
class KVGenPhaseSpace : public KVBase {
private:
   void init();
 
protected:
    TObject* fMCSampler;
    KVNucleus fCompound;
    KVEvent* fEvent;
    Int_t fMult;
    std::vector<Double_t> fMasses;
    Bool_t fOK;
    Double_t fEtot;
 
   Bool_t CheckBreakUpChannel();
   virtual void InitialiseMCSampler();
 
public:
   KVGenPhaseSpace();
   KVGenPhaseSpace(const Char_t* name, const Char_t* title = "");
   virtual ~KVGenPhaseSpace();
   void Copy(TObject& obj) const override;
 
   Bool_t SetBreakUpChannel(const KVNucleus& CN, KVEvent* e);
 
   virtual Double_t Generate();
 
   Double_t GetAvailableEnergy() const
   {
      return fEtot;
   }
 
   ClassDefOverride(KVGenPhaseSpace, 1) //Generate momenta for an event using microcanonical phase space sampling
};
 
#endif