PSaid.h

Go to the documentation of this file.

// Author: M.A. Kagarlis
// Written:  02.02.00
// Revised:  19.06.03  R.H.
// Revised:  26.06.07  IF
// PSaid Class Header
 
// R. Arndt's p+p elastic scattering parametrization
#ifndef _PSAID_H_
#define _PSAID_H_
 
#include <TGlobal.h>
#include <TROOT.h>
#include <TRandom3.h>
#include <TArrayD.h>
#include "PAngularDistribution.h"
 
class PSaid : public PAngularDistribution {
 
 public:
 
    PSaid();
    PSaid(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);
 
    ~PSaid();
 
 private:
 
    double told;
 
    int dim;
    //TRandom3 *REngine;                   // private M-Twistor random number engine
    TArrayD y, aa;                       // private arrays for the sampling algorithm
    void defaults();
 
    double dsdw(double, double);   // diff. cross section (mb/sr) by sc. angle (deg), Tlab (GeV)
    double sample(double, double); // sampling algorithm (rejection method)
 
    ClassDef(PSaid, 0) //Pluto SAID Class
 
};
#endif // _PSAID_H_