PHadronDecay.h

Go to the documentation of this file.

// Author: I. Froehlich
// Written: 27.5.2007
// Revised:
 
#ifndef _PHADRONDECAY_H_
#define _PHADRONDECAY_H_
 
#include "TF1.h"
#include "TF2.h"
#include "PChannelModel.h"
#include "PDynamicData.h"
#include "PKinematics.h"
 
class PHadronDecay : public PChannelModel  {
 
 public:
    PHadronDecay();
    PHadronDecay(const Char_t *id, const Char_t *de, Int_t key);
    PDistribution *Clone(const char*delme=NULL) const;
 
    Bool_t   Init(void);
    Bool_t   SampleMass(void);
    Bool_t   SampleMass(Double_t *mass, Int_t *didx=NULL);
 
    Bool_t   GetWidth(Double_t mass, Double_t *width, Int_t didx=-1);
 
    int      GetDepth(int i=0);
//    Bool_t GetBR(Double_t mass, Double_t *br, Double_t totalwidth=-1);
 
    Double_t Eval(Double_t x, Double_t y = 0, Double_t z = 0, Double_t t = 0) const;
    Double_t EvalPar(const Double_t *x, const Double_t *params);
    //TF1 wrapper
 
    void Use_m0_over_m(int i)    {use_m0_over_m = i;};
    void SetCutoffVersion(int i) {cutoff_version = i;};
 
 protected:
 
    int    parent_id, id1, id2;          //PIDs
    double parent_mass, mass1, mass2;  //Static Masses
    double parent_g0; //Static Width of parent
 
    //Width Configuration
    int use_fixed_delta;
    double fixed_delta;
    int angular_l;
    int cutoff_l;
    int use_m0_over_m, cutoff_version; //0=theis,1=ernst,2=none
    double w0;
 
    double HadronWidth(const double &m, const double &ma, const double &mb);
 
    ClassDef(PHadronDecay, 0)  // Decay of Hadron -> Hadron(stable) + Hadron(stable)
};
 
#endif