52 for (
int i = 0; i < 40; i++) {
53 h_qcss_Evtof1[i] =
new TH2F(TString::Format(
"qcss_Evtof1_%d", i),
"Energy deposited [MeV] vs TOF [ns] - all", 500, 0., 1000.,
55 h_qcss_Evtof2[i] =
new TH2F(TString::Format(
"qcss_Evtof2_%d", i),
"Energy deposited [MeV] vs TOF [ns] - only photons", 500, 0.,
57 h_qcss_Evtof3[i] =
new TH2F(TString::Format(
"qcss_Evtof3_%d", i),
"Energy deposited [MeV] vs TOF [ns] - only e+/e-", 500, 0.,
59 h_qcss_Evtof4[i] =
new TH2F(TString::Format(
"qcss_Evtof4_%d", i),
"Energy deposited [MeV] vs TOF [ns] - only e+/e-", 500, 0.,
61 h_qcss_edep[i] =
new TH1F(TString::Format(
"qcss_edep_%d", i),
"Energy deposited [MeV]", 5000, 0., 10.);
62 h_Wqcss_edep[i] =
new TH1F(TString::Format(
"Wqcss_edep_%d", i),
"Energy deposited [MeV]", 5000, 0., 10.);
65 h_qcss_hitrate0 =
new TH1F(
"qcss_hitrate0",
"Hit distributions", 100, 0., 100.);
66 h_qcss_hitrate1 =
new TH1F(
"qcss_hitrate1",
"Hit distributions", 100, 0., 100.);
67 h_qcss_hitrate2 =
new TH1F(
"qcss_hitrate2",
"Hit distributions", 100, 0., 100.);
68 h_qcss_hitrate1W =
new TH1F(
"qcss_hitrate1W",
"Hit distributions", 100, 0., 100.);
69 h_qcss_hitrate2W =
new TH1F(
"qcss_hitrate2W",
"Hit distributions", 100, 0., 100.);
77 h_qcss_rs_hitrate1 =
new TH2F(
"qcss_rs_hitrate1",
"Hit distributions vs rs", 100, 0., 100., 12, 0., 12.);
78 h_qcss_rs_hitrate2 =
new TH2F(
"qcss_rs_hitrate2",
"Hit distributions vs rs", 100, 0., 100., 12, 0., 12.);
79 h_qcss_rs_hitrate1W =
new TH2F(
"qcss_rs_hitrate1W",
"Hit distributions vs rs", 100, 0., 100., 12, 0., 12.);
80 h_qcss_rs_hitrate2W =
new TH2F(
"qcss_rs_hitrate2W",
"Hit distributions vs rs", 100, 0., 100., 12, 0., 12.);
87 for (
int i = 0; i < 40; i++) {
88 h_qcss_rate1[i] =
new TH1F(TString::Format(
"qcss_rate1_%d", i),
"PE distributions", 500, 0., 500.);
89 h_qcss_rate2[i] =
new TH1F(TString::Format(
"qcss_rate2_%d", i),
"PE distributions", 500, 0., 500.);
90 h_qcss_rate1W[i] =
new TH1F(TString::Format(
"qcss_rate1W_%d", i),
"PE distributions", 500, 0., 500.);
91 h_qcss_rate2W[i] =
new TH1F(TString::Format(
"qcss_rate2W_%d", i),
"PE distributions", 500, 0., 500.);
92 h_qcss_pe1[i] =
new TH2F(TString::Format(
"qcss_pe1_%d", i),
"PE distributions", 500, 0., 500., 100, 0., 1000.);
93 h_qcss_pe2[i] =
new TH2F(TString::Format(
"qcss_pe2_%d", i),
"PE distributions", 500, 0., 500., 100, 0., 1000.);
94 h_qcss_pe1W[i] =
new TH2F(TString::Format(
"qcss_pe1W_%d", i),
"PE distributions", 500, 0., 500., 100, 0., 1000.);
95 h_qcss_pe2W[i] =
new TH2F(TString::Format(
"qcss_pe2W_%d", i),
"PE distributions", 500, 0., 500., 100, 0., 1000.);
97 h_qcss_rs_rate1[i] =
new TH2F(TString::Format(
"qcss_rs_rate1_%d", i),
"PE distributions", 500, 0., 500., 12, 0., 12.);
98 h_qcss_rs_rate2[i] =
new TH2F(TString::Format(
"qcss_rs_rate2_%d", i),
"PE distributions", 500, 0., 500., 12, 0., 12.);
99 h_qcss_rs_rate1W[i] =
new TH2F(TString::Format(
"qcss_rs_rate1W_%d", i),
"PE distributions", 500, 0., 500., 12, 0., 12.);
100 h_qcss_rs_rate2W[i] =
new TH2F(TString::Format(
"qcss_rs_rate2W_%d", i),
"PE distributions", 500, 0., 500., 12, 0., 12.);
140 int ring_section = -1;
141 for (
const auto& MetaHit : MetaHits) {
142 rate = MetaHit.getrate();
143 ring_section = MetaHit.getring_section() - 1;
150 for (
int i = 0; i < nSimHits; i++) {
171 for (
const auto&
Hit : Hits) {
172 const int detNb =
Hit.getdetNb();
174 const int timebin =
Hit.gettime();
175 const float edep =
Hit.getedep();
176 const float pe =
Hit.getPE();
void addParam(const std::string &name, T ¶mVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module.