PHadronDecayM1N.h

Go to the documentation of this file.

// Author: I. Froehlich
// Written: 22.01.2010
// Revised:
 
#ifndef _PHADRONDECAYM1N_H_
#define _PHADRONDECAYM1N_H_
 
#include "TF1.h"
#include "TF2.h"
#include "TGenPhaseSpace.h"
#include "PChannelModel.h"
#include "PDynamicData.h"
#include "PKinematics.h"
#include "PMesh.h"
 
 
#define M1N_MAX_DAUGHTERS 10
#define M1N_PARENT_GRID   0.1
 
class PHadronDecayM1N : public PChannelModel  {
 
 public:
    PHadronDecayM1N();
    PHadronDecayM1N(const Char_t *id, const Char_t *de, Int_t key);
    PDistribution *Clone(const char *delme=NULL) const;
 
 
    Bool_t Prepare(void);
    Bool_t Init(void);
    Bool_t SampleMass(void);
    using  PChannelModel::SampleMass;
 
    int      GetDepth(int i=0);
 
    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 SubPrint(Int_t opt) const ;  //Print sub-models
 
 protected:
 
    //int didx1,didx2,didx3;
    PParticle * parent, *daughters[M1N_MAX_DAUGHTERS], *unstable_particle;
    Int_t daughter_pos, unstable_daughter;
    PChannelModel *unstable_model;
 
    Int_t unstable_pos,
    daughter_id[M1N_MAX_DAUGHTERS],
    daughter_didx[M1N_MAX_DAUGHTERS],
    unstable_id, parent_id, unstable_didx;
    Double_t unstable_width, daughter_masses[M1N_MAX_DAUGHTERS], parent_mass, old_parent_mass;
 
    PMesh *mesh;
    TGenPhaseSpace *event;
 
    void UpdateMesh(void);
    Int_t maxmesh;
 
    ClassDef(PHadronDecayM1N, 0)  // Decay of Hadron -> 1 unstable + N stable
};
 
#endif