forked from sPHENIX-Collaboration/coresoftware
/
PHG4ParticleGeneratorD0.h
49 lines (40 loc) · 1.26 KB
/
PHG4ParticleGeneratorD0.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
// Tell emacs that this is a C++ source
// -*- C++ -*-.
#ifndef G4MAIN_PHG4PARTICLEGENERATORD0_H
#define G4MAIN_PHG4PARTICLEGENERATORD0_H
#include "PHG4ParticleGeneratorBase.h"
#include <string> // for string
class PHCompositeNode;
class TF1;
class PHG4ParticleGeneratorD0 : public PHG4ParticleGeneratorBase
{
public:
PHG4ParticleGeneratorD0(const std::string &name = "D0GEN");
~PHG4ParticleGeneratorD0() override {}
int InitRun(PHCompositeNode *topNode) override;
int process_event(PHCompositeNode *topNode) override;
void set_eta_range(const double eta_min, const double eta_max);
void set_rapidity_range(const double y_min, const double y_max);
void set_mom_range(const double mom_min, const double mom_max);
void set_pt_range(const double pt_min, const double pt_max);
void set_vtx_zrange(const double zmin, const double zmax);
void set_mass(const double mass);
protected:
double vtx_zmin = -10.;
double vtx_zmax = 10.;
double y_min = 0.;
double y_max = 0.;
double eta_min = -1.;
double eta_max = 1.;
double mom_min = 0.;
double mom_max = 10.;
double pt_min = 4.;
double pt_max = 4.;
double mass = 1.86486;
double m1 = 0.493677;
double m2 = 0.13957018;
TF1 *fsin = nullptr;
TF1 *frap = nullptr;
TF1 *fpt = nullptr;
};
#endif