// PSaid Class Header


#ifndef _PSAIDLOWENERGY_H_
#define _PSAIDLOWENERGY_H_

#include <TGlobal.h>
#include <TROOT.h>
#include <TRandom3.h>
#include <TArrayD.h>
#include "PAngularDistribution.h"

class PSaidLowEnergy : public PAngularDistribution {

 public:
    
    PSaidLowEnergy();
    PSaidLowEnergy(const Char_t *id, const Char_t *de);
    
    PDistribution *Clone(const char *delme=NULL) const;

    Bool_t Init(void);
    Bool_t Prepare(void);
    Bool_t IsNotRejected(void);
        
    double SamplePolarAngle(double);

 private:

    double told;
    int dim;
    TArrayD y, aa;                       // private arrays for the sampling algorithm

    double dsdw(double, double);   // diff. cross section (mb/sr) by sc. angle (deg), Tlab (GeV)
    double sample(double, double); // sampling algorithm (rejection method)

    ClassDef(PSaidLowEnergy, 0) //Pluto SAID Class (low energy)

};
#endif // _PSAILOWENERGYD_H_