PHadronDecayM2Mesh.h

Go to the documentation of this file.

// Author: I. Froehlich
// Written: 19.08.2026
// Revised:
 
#ifndef _PHADRONDECAYM2MESH_H_
#define _PHADRONDECAYM2MESH_H_
 
#include "TF1.h"
#include "TF2.h"
#include "PHadronDecayM2.h"
#include "PDynamicData.h"
#include "PKinematics.h"
 
class PHadronDecayM2Mesh : public PHadronDecayM2  {
 
 public:
    PHadronDecayM2Mesh();
    PHadronDecayM2Mesh(const Char_t *id, const Char_t *de, Int_t key);
    PDistribution *Clone(const char *delme=NULL) const;
 
    Bool_t Init(void);
    Bool_t Prepare(void);
    void   SubPrint(Int_t opt) const;
    int    GetDepth(int i=0);
 
 protected:
 
    //Generated by Gemini:
    static const int N_ECM = 100; // Number of ECM energy slices
    static const int N1 = 20;     // m1 resolution per slice
    static const int N2 = 20;     // m2 resolution per slice
    static const int MAX_CELLS = N1 * N2;
 
    Double_t ecm_grid_min;
    Double_t ecm_grid_max;
    Double_t ecm_step;
 
    // Cache arrays for each ECM slice
    Bool_t is_calculated[N_ECM];
    Double_t cell_ceiling[N_ECM][N1][N2];
    Double_t total_weight_sum[N_ECM];
    Double_t cell_cum_weight[N_ECM][MAX_CELLS];
    Int_t cell_i[N_ECM][MAX_CELLS];
    Int_t cell_j[N_ECM][MAX_CELLS];
    Int_t valid_cells_cnt[N_ECM];
 
    void initGridBounds();
    void calculateSlice(int slice_idx, double ecm_slice);
    //End generated by Gemini
    
    
    bool sampleM2(const double &ecm);
 
    ClassDef(PHadronDecayM2Mesh, 0)  // Hadron decay in 2 unstable products using meshes
};
 
#endif