Belle II Software  release-08-02-04
TRGEFFDQMModule.cc
1 /**************************************************************************
2  * basf2 (Belle II Analysis Software Framework) *
3  * Author: The Belle II Collaboration *
4  * *
5  * See git log for contributors and copyright holders. *
6  * This file is licensed under LGPL-3.0, see LICENSE.md. *
7  **************************************************************************/
8 
9 // Own include
10 #include <trg/gdl/modules/trggdlDQM/TRGEFFDQMModule.h>
11 
12 // Dataobject classes
13 #include <TF1.h>
14 #include <Math/Vector3D.h>
15 #include <TVector3.h>
16 #include <TDirectory.h>
17 #include <fstream>
18 #include <math.h>
19 #include <mdst/dataobjects/HitPatternCDC.h>
20 
21 using namespace std;
22 using namespace Belle2;
23 
24 //-----------------------------------------------------------------
25 // Register module
26 //-----------------------------------------------------------------
27 
28 REG_MODULE(TRGEFFDQM);
29 
30 TRGEFFDQMModule::TRGEFFDQMModule() : HistoModule()
31 {
32  // set module description (e.g. insert text)
33  setDescription("Make kinematics dependent efficiency plot");
34  setPropertyFlags(c_ParallelProcessingCertified);
35 }
36 
37 TRGEFFDQMModule::~TRGEFFDQMModule()
38 {
39 }
40 
41 void TRGEFFDQMModule::defineHisto()
42 {
43  TDirectory* oldDir = gDirectory;
44  oldDir->mkdir("TRGEFF");
45  oldDir->cd("TRGEFF");
46 
47  m_hPhi = new TH1F("hPhi", "", 36, -180.0, 180.0);
48  m_hPhi_psnecl = new TH1F("hPhi_psnecl", "", 36, -180.0, 180.0);
49  m_hPhi_psnecl_ftdf = new TH1F("hPhi_psnecl_ftdf", "", 36, -180.0, 180.0);
50 
51  m_hPt = new TH1F("hPt", "", 50, 0, 5);
52  m_hPt_psnecl = new TH1F("hPt_psnecl", "", 50, 0, 5);
53  m_hPt_psnecl_ftdf = new TH1F("hPt_psnecl_ftdf", "", 50, 0, 5);
54 
55  m_nobha_hPt = new TH1F("nobha_hPt", "", 50, 0, 5);
56  m_nobha_hPt_psnecl = new TH1F("nobha_hPt_psnecl", "", 50, 0, 5);
57  m_nobha_hPt_psnecl_ftdf = new TH1F("nobha_hPt_psnecl_ftdf", "", 50, 0, 5);
58 
59  m_hP3_y = new TH1F("hP3_y", "", 50, 0, 5);
60  m_hP3_y_psnecl = new TH1F("hP3_y_psnecl", "", 50, 0, 5);
61  m_hP3_y_psnecl_ftdf = new TH1F("hP3_y_psnecl_ftdf", "", 50, 0, 5);
62 
63  m_hP3_z = new TH1F("hP3_z", "", 50, 0, 5);
64  m_hP3_z_psnecl = new TH1F("hP3_z_psnecl", "", 50, 0, 5);
65  m_hP3_z_psnecl_ftdf = new TH1F("hP3_z_psnecl_ftdf", "", 50, 0, 5);
66 
67  m_nobha_hP3_y = new TH1F("nobha_hP3_y", "", 50, 0, 5);
68  m_nobha_hP3_y_psnecl = new TH1F("nobha_hP3_y_psnecl", "", 50, 0, 5);
69  m_nobha_hP3_y_psnecl_ftdf = new TH1F("nobha_hP3_y_psnecl_ftdf", "", 50, 0, 5);
70 
71  m_nobha_hP3_z = new TH1F("nobha_hP3_z", "", 50, 0, 5);
72  m_nobha_hP3_z_psnecl = new TH1F("nobha_hP3_z_psnecl", "", 50, 0, 5);
73  m_nobha_hP3_z_psnecl_ftdf = new TH1F("nobha_hP3_z_psnecl_ftdf", "", 50, 0, 5);
74 
75  m_fyo_dphi = new TH1F("fyo_dphi", "", 18, 0., 180.);
76  m_fyo_dphi_psnecl = new TH1F("fyo_dphi_psnecl", "", 18, 0., 180.);
77  m_fyo_dphi_psnecl_ftdf = new TH1F("fyo_dphi_psnecl_ftdf", "", 18, 0., 180.);
78 
79  m_nobha_fyo_dphi = new TH1F("nobha_fyo_dphi", "", 18, 0., 180.);
80  m_nobha_fyo_dphi_psnecl = new TH1F("nobha_fyo_dphi_psnecl", "", 18, 0., 180.);
81  m_nobha_fyo_dphi_psnecl_ftdf = new TH1F("nobha_fyo_dphi_psnecl_ftdf", "", 18, 0., 180.);
82 
83  m_stt_phi = new TH1F("stt_phi", "", 36, -180., 180.);
84  m_stt_phi_psnecl = new TH1F("stt_phi_psnecl", "", 36, -180., 180.);
85  m_stt_phi_psnecl_ftdf = new TH1F("stt_phi_psnecl_ftdf", "", 36, -180., 180.);
86 
87  m_stt_P3 = new TH1F("stt_P3", "", 50, 0, 5);
88  m_stt_P3_psnecl = new TH1F("stt_P3_psnecl", "", 50, 0, 5);
89  m_stt_P3_psnecl_ftdf = new TH1F("stt_P3_psnecl_ftdf", "", 50, 0, 5);
90 
91  m_stt_theta = new TH1F("stt_theta", "", 18, 0, 180);
92  m_stt_theta_psnecl = new TH1F("stt_theta_psnecl", "", 18, 0, 180);
93  m_stt_theta_psnecl_ftdf = new TH1F("stt_theta_psnecl_ftdf", "", 18, 0, 180);
94 
95  m_nobha_stt_phi = new TH1F("nobha_stt_phi", "", 36, -180., 180.);
96  m_nobha_stt_phi_psnecl = new TH1F("nobha_stt_phi_psnecl", "", 36, -180., 180.);
97  m_nobha_stt_phi_psnecl_ftdf = new TH1F("nobha_stt_phi_psnecl_ftdf", "", 36, -180., 180.);
98 
99  m_nobha_stt_P3 = new TH1F("nobha_stt_P3", "", 50, 0, 5);
100  m_nobha_stt_P3_psnecl = new TH1F("nobha_stt_P3_psnecl", "", 50, 0, 5);
101  m_nobha_stt_P3_psnecl_ftdf = new TH1F("nobha_stt_P3_psnecl_ftdf", "", 50, 0, 5);
102 
103  m_nobha_stt_theta = new TH1F("nobha_stt_theta", "", 18, 0, 180);
104  m_nobha_stt_theta_psnecl = new TH1F("nobha_stt_theta_psnecl", "", 18, 0, 180);
105  m_nobha_stt_theta_psnecl_ftdf = new TH1F("nobha_stt_theta_psnecl_ftdf", "", 18, 0, 180);
106 
107  m_hie_E = new TH1F("hie_E", "", 60, 0, 12);
108  m_hie_E_psnecl = new TH1F("hie_E_psnecl", "", 60, 0, 12);
109  m_hie_E_psnecl_ftdf = new TH1F("hie_E_psnecl_ftdf", "", 60, 0, 12);
110 
111  m_nobha_hie_E = new TH1F("nobha_hie_E", "", 60, 0, 12);
112  m_nobha_hie_E_psnecl = new TH1F("nobha_hie_E_psnecl", "", 60, 0, 12);
113  m_nobha_hie_E_psnecl_ftdf = new TH1F("nobha_hie_E_psnecl_ftdf", "", 60, 0, 12);
114 
115  m_ecltiming_E = new TH1F("ecltiming_E", "", 60, 0, 12);
116  m_ecltiming_E_psnecl = new TH1F("ecltiming_E_psnecl", "", 60, 0, 12);
117  m_ecltiming_E_psnecl_ftdf = new TH1F("ecltiming_E_psnecl_ftdf", "", 60, 0, 12);
118 
119  m_ecltiming_theta = new TH1F("ecltiming_theta", "", 18, 0, 180);
120  m_ecltiming_theta_psnecl = new TH1F("ecltiming_theta_psnecl", "", 18, 0, 180);
121  m_ecltiming_theta_psnecl_ftdf = new TH1F("ecltiming_theta_psnecl_ftdf", "", 18, 0, 180);
122 
123  m_ecltiming_phi = new TH1F("ecltiming_phi", "", 36, -180., 180.);
124  m_ecltiming_phi_psnecl = new TH1F("ecltiming_phi_psnecl", "", 36, -180., 180.);
125  m_ecltiming_phi_psnecl_ftdf = new TH1F("ecltiming_phi_psnecl_ftdf", "", 36, -180., 180.);
126 
127  m_klmhit_phi = new TH1F("klmhit_phi", "", 18, -180., 180.);
128  m_klmhit_phi_psnecl = new TH1F("klmhit_phi_psnecl", "", 18, -180., 180.);
129  m_klmhit_phi_psnecl_ftdf = new TH1F("klmhit_phi_psnecl_ftdf", "", 18, -180., 180.);
130 
131  m_klmhit_theta = new TH1F("klmhit_theta", "", 18, 0, 180);
132  m_klmhit_theta_psnecl = new TH1F("klmhit_theta_psnecl", "", 18, 0, 180);
133  m_klmhit_theta_psnecl_ftdf = new TH1F("klmhit_theta_psnecl_ftdf", "", 18, 0, 180);
134 
135  m_eklmhit_phi = new TH1F("eklmhit_phi", "", 18, -180., 180.);
136  m_eklmhit_phi_psnecl = new TH1F("eklmhit_phi_psnecl", "", 18, -180., 180.);
137  m_eklmhit_phi_psnecl_ftdf = new TH1F("eklmhit_phi_psnecl_ftdf", "", 18, -180., 180.);
138 
139  m_eklmhit_theta = new TH1F("eklmhit_theta", "", 18, 0, 180);
140  m_eklmhit_theta_psnecl = new TH1F("eklmhit_theta_psnecl", "", 18, 0, 180);
141  m_eklmhit_theta_psnecl_ftdf = new TH1F("eklmhit_theta_psnecl_ftdf", "", 18, 0, 180);
142 
143 
144 
145 
146 
147 
148 
149  oldDir->cd();
150 }
151 
152 void TRGEFFDQMModule::initialize()
153 {
154  REG_HISTOGRAM
155 
156  if (!m_Tracks.isOptional()) {
157  B2WARNING("Missing Tracks array");
158  return;
159  }
160  if (!m_ECLClusters.isOptional()) {
161  B2WARNING("Missing ECLCLuster array");
162  return;
163  }
164  if (!m_KLMClusters.isOptional()) {
165  B2WARNING("Missing KLMCLuster array");
166  return;
167  }
168  if (!m_trgSummary.isOptional()) {
169  B2WARNING("Missing TRGSummary");
170  return;
171  }
172 }
173 
174 void TRGEFFDQMModule::beginRun()
175 {
176  if (!m_RecoTracks.isOptional()) {
177  B2DEBUG(22, "Missing recoTracks array in beginRun() ");
178  return;
179  }
180 
181  m_hPhi->Reset();
182 }
183 
184 void TRGEFFDQMModule::event()
185 {
186  // bool debug = 1; // if(debug)cout<<"line "<<__LINE__<<endl;
187 
188  if (!m_trgSummary.isValid()) {
189  B2WARNING("TRGSummary object not available but require to estimate trg efficiency");
190  return;
191  }
192 
193  if (!m_TrgResult.isValid()) {
194  B2WARNING("SoftwareTriggerResult object not available but require to select bhabha/mumu/hadron events skim");
195  return;
196  }
197 
198  const std::map<std::string, int>& fresults = m_TrgResult->getResults();
199  if ((fresults.find("software_trigger_cut&skim&accept_bhabha") == fresults.end())
200  || (fresults.find("software_trigger_cut&skim&accept_hadron") == fresults.end())) {
201  B2WARNING("TRGEFFDQMModule: Can't find required bhabha or mumu or hadron trigger identifier");
202  return;
203  }
204 
205  const bool IsBhabha = (m_TrgResult->getResult("software_trigger_cut&skim&accept_bhabha") == SoftwareTriggerCutResult::c_accept);
206  const bool IsHadron = (m_TrgResult->getResult("software_trigger_cut&skim&accept_hadron") == SoftwareTriggerCutResult::c_accept);
207 
208 
209  /*///////////////////////////////////////////////////////////////////
212 
213  // calculate the total energy
214  double E_ecl_all = 0; // the ECL total energy
215  double E_ecl_hie =
216  0; // the ECL total energy in the thetaID range 2<=ThetaID<=15 (corresponds to 22.49<=theta<=126.80) for ehigh bit
217  for (const auto& test_b2eclcluster : m_ECLClusters) {
218  if (!(test_b2eclcluster.hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons))) continue;
219  double energy = test_b2eclcluster.getEnergyRaw();
220  double theta = test_b2eclcluster.getTheta() / Unit::deg;
221 
222  if (energy < 0.1) continue;
223 
224  E_ecl_all = E_ecl_all + energy;
225  if (theta >= 22.49 && theta <= 126.8) {
226  E_ecl_hie = E_ecl_hie + energy;
227  }
228  }
229 
230  bool trg_hie_psncdc = 0; // for ECL energy trigger, for hie
231  bool trg_hie_Eecl = 0; // for ECL energy trigger, for hie
232  bool trg_ecltiming_psncdc = 0; // for ECL energy trigger, for ecltiming
233  bool trg_ecltiming_Eecl = 0; // for ECL energy trigger. for ecltiming
234  bool trg_nobha_hie_Eecl = 0;
235 
236  for (const auto& b2eclcluster : m_ECLClusters) {
237  if (!(b2eclcluster.hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons))) continue;
238  double phi = b2eclcluster.getPhi() / Unit::deg;
239  double theta = b2eclcluster.getTheta() / Unit::deg;
240  double energy = b2eclcluster.getEnergyRaw();
241 
242  if (energy < 0.1) continue;
243 
244  bool trg_psncdc = m_trgSummary->testPsnm("ffy") || m_trgSummary->testPsnm("fyo") || m_trgSummary->testPsnm("stt");
245  bool trg_hie = m_trgSummary->testFtdl("hie");
246  bool trg_hie_nobha = m_trgSummary->testInput("ehigh"); // remove the bha_veto for hie bit
247  bool trg_ecltiming = m_trgSummary->testFtdl("ecltiming");
248 
249  m_ecltiming_theta->Fill(theta);
250  m_ecltiming_phi->Fill(phi);
251 
252  if (trg_psncdc) {
253  trg_hie_psncdc = 1;
254  trg_ecltiming_psncdc = 1;
255 
256  m_ecltiming_theta_psnecl->Fill(theta);
257  m_ecltiming_phi_psnecl->Fill(phi);
258  }
259 
260  if (trg_psncdc && trg_hie) {
261  trg_hie_Eecl = 1;
262  }
263  if (trg_psncdc && trg_hie_nobha) {
264  trg_nobha_hie_Eecl = 1;
265  }
266  if (trg_psncdc && trg_ecltiming) {
267  trg_ecltiming_Eecl = 1;
268 
269  m_ecltiming_theta_psnecl_ftdf->Fill(theta);
270  m_ecltiming_phi_psnecl_ftdf->Fill(phi);
271  }
272  }
273 
274  // sum of the ECL total energy
275  m_hie_E->Fill(E_ecl_hie);
276  m_nobha_hie_E->Fill(E_ecl_hie);
277  if (trg_hie_psncdc) {
278  m_hie_E_psnecl->Fill(E_ecl_hie);
279  m_nobha_hie_E_psnecl->Fill(E_ecl_hie);
280  }
281  if (trg_hie_Eecl) {
282  m_hie_E_psnecl_ftdf->Fill(E_ecl_hie);
283  }
284  if (trg_nobha_hie_Eecl) {
285  m_nobha_hie_E_psnecl_ftdf->Fill(E_ecl_hie);
286  }
287 
288  m_ecltiming_E->Fill(E_ecl_all);
289  if (trg_ecltiming_psncdc) {
290  m_ecltiming_E_psnecl->Fill(E_ecl_all);
291  }
292  if (trg_ecltiming_Eecl) {
293  m_ecltiming_E_psnecl_ftdf->Fill(E_ecl_all);
294  }
295 
296 
297 
298 
299 
300 
301  /*/
304  for (const auto& b2klmcluster : m_KLMClusters) {
305  int nlayer = b2klmcluster.getLayers();
306  // ROOT::Math::XYZVector position = b2klmcluster.getClusterPosition();
307 
308  if (nlayer <= 6)
309  continue;
310 
311  // double p3 = b2klmcluster.getMomentum().R(); // 3-momentum
312  double theta = b2klmcluster.getMomentum().Theta() / Unit::deg;
313  double phiDegree = b2klmcluster.getMomentum().Phi() / Unit::deg;
314 
315  bool trg_KLMecl = m_trgSummary->testPsnm("hie") || m_trgSummary->testPsnm("c4") || m_trgSummary->testPsnm("eclmumu") ||
316  m_trgSummary->testPsnm("lml1") || m_trgSummary->testPsnm("lml2") || m_trgSummary->testPsnm("lml6") ||
317  m_trgSummary->testPsnm("lml7") || m_trgSummary->testPsnm("lml8") || m_trgSummary->testPsnm("lml9") ||
318  m_trgSummary->testPsnm("lml10");
319  // these trigger flags are not available in older data
320  try {
321  trg_KLMecl |= m_trgSummary->testPsnm("ecltaub2b3") || m_trgSummary->testPsnm("hie4");
322  } catch (const std::exception&) {
323  }
324 
325 
326  bool trg_klmhit = m_trgSummary->testFtdl("klmhit");
327  bool trg_eklmhit = m_trgSummary->testFtdl("eklmhit");
328 
329  m_klmhit_phi->Fill(phiDegree);
330  m_klmhit_theta->Fill(theta);
331  m_eklmhit_phi->Fill(phiDegree);
332  m_eklmhit_theta->Fill(theta);
333 
334  if (trg_KLMecl) {
335  m_klmhit_phi_psnecl->Fill(phiDegree);
336  m_klmhit_theta_psnecl->Fill(theta);
337  m_eklmhit_phi_psnecl->Fill(phiDegree);
338  m_eklmhit_theta_psnecl->Fill(theta);
339  }
340  if (trg_KLMecl && trg_klmhit) {
341  m_klmhit_phi_psnecl_ftdf->Fill(phiDegree);
342  m_klmhit_theta_psnecl_ftdf->Fill(theta);
343  }
344  if (trg_KLMecl && trg_eklmhit) {
345  m_eklmhit_phi_psnecl_ftdf->Fill(phiDegree);
346  m_eklmhit_theta_psnecl_ftdf->Fill(theta);
347  }
348  }
349 
350 
351 
352  /*///////////////////////////////////////////////////////////////////
355  // int N_tracks = m_Tracks.getEntries(); // number of tracks
356 
357  vector<double> p_stt_P3_psnecl_ftdf, p_stt_P3_psnecl, p_stt_P3, phi_fyo_dphi, phi_fyo_dphi_psnecl, phi_fyo_dphi_psnecl_ftdf ;
358  p_stt_P3_psnecl_ftdf.clear();
359  p_stt_P3_psnecl.clear();
360  p_stt_P3.clear();
361 
362  phi_fyo_dphi.clear();
363  phi_fyo_dphi_psnecl.clear();
364  phi_fyo_dphi_psnecl_ftdf.clear();
365 
366  vector<double> p_nobha_stt_P3_psnecl_ftdf, p_nobha_stt_P3_psnecl, p_nobha_stt_P3, phi_nobha_fyo_dphi, phi_nobha_fyo_dphi_psnecl,
367  phi_nobha_fyo_dphi_psnecl_ftdf ;
368  p_nobha_stt_P3_psnecl_ftdf.clear();
369  p_nobha_stt_P3_psnecl.clear();
370  p_nobha_stt_P3.clear();
371 
372  phi_nobha_fyo_dphi.clear();
373  phi_nobha_fyo_dphi_psnecl.clear();
374  phi_nobha_fyo_dphi_psnecl_ftdf.clear();
375 
376  int nitrack = 0; // the i track of the m_Tracks
377  for (const auto& b2track : m_Tracks) {
378  const Belle2::TrackFitResult* fitresult = b2track.getTrackFitResultWithClosestMass(Const::pion);
379  if (!fitresult) {
380  B2WARNING("No track fit result found.");
381  nitrack++;
382  continue;
383  }
384 
385  // require high NDF track
386  int ndf = fitresult->getNDF();
387  if (ndf < 20) {
388  nitrack++;
389  continue;
390  }
391 
392  // IP tracks at barrel
393  if (fabs(fitresult->getD0()) < 1.0 && fabs(fitresult->getZ0()) < 1.0 && fitresult->getHitPatternCDC().getLastLayer() > 50
394  && fitresult->getHitPatternCDC().getFirstLayer() < 5) {
395  double phiDegree = fitresult->getPhi() / Unit::deg;
396  double pt = fitresult->getTransverseMomentum();
397  double p3 = fitresult->getMomentum().R(); // 3-momentum
398  double theta = fitresult->getMomentum().Theta() / Unit::deg;
399 
400  bool trg_psnecl = m_trgSummary->testPsnm("hie") || m_trgSummary->testPsnm("c4") || m_trgSummary->testPsnm("eclmumu") ||
401  m_trgSummary->testPsnm("lml1") || m_trgSummary->testPsnm("lml2") || m_trgSummary->testPsnm("lml6") ||
402  m_trgSummary->testPsnm("lml7") || m_trgSummary->testPsnm("lml8") || m_trgSummary->testPsnm("lml9") ||
403  m_trgSummary->testPsnm("lml10");
404  // these trigger flags are not available in older data
405  try {
406  trg_psnecl |= m_trgSummary->testPsnm("ecltaub2b3") || m_trgSummary->testPsnm("hie4");
407  } catch (const std::exception&) {
408  }
409 
410  bool trg_ftdf = m_trgSummary->testFtdl("f");
411 
412  // for f bit, reomove the Bhabha_veto
413  bool trg_itdt2 = (m_trgSummary->testInput("t2_0") || m_trgSummary->testInput("t2_1") || m_trgSummary->testInput("t2_2")
414  || m_trgSummary->testInput("t2_3"));
415 
416  // for z, reomove the Bhabha_veto
417  bool trg_itdt3 = (m_trgSummary->testInput("t3_0") || m_trgSummary->testInput("t3_1") || m_trgSummary->testInput("t3_2")
418  || m_trgSummary->testInput("t3_3"));
419 
420  // for y, reomove the Bhabha_veto
421  bool trg_itdt4 = (m_trgSummary->testInput("ty_0") || m_trgSummary->testInput("ty_1") || m_trgSummary->testInput("ty_2")
422  || m_trgSummary->testInput("ty_3"));
423 
424  // (t3>0 and !bhaveto and !veto) for z
425  bool trg_ftdz = m_trgSummary->testFtdl("z");
426 
427  // (ty>0 and !bhaveto and !veto) for y
428  bool trg_ftdy = m_trgSummary->testFtdl("y");
429 
430  // typ and !bha veto and !veto
431  bool trg_stt = m_trgSummary->testFtdl("stt");
432 
433  // remove bha_veto
434  bool trg_stt_nobha = m_trgSummary->testInput("typ") ;
435 
436 
437  // require pt > 0.3 GeV
438  if (pt > 0.3) {
439  m_hPhi->Fill(phiDegree);
440  if (trg_psnecl) {
441  m_hPhi_psnecl->Fill(phiDegree);
442  }
443  if (trg_psnecl && trg_ftdf) {
444  m_hPhi_psnecl_ftdf->Fill(phiDegree);
445  }
446  }
447 
448  m_hPt->Fill(pt);
449  m_nobha_hPt->Fill(pt);
450 
451  m_hP3_z->Fill(p3);
452  m_hP3_y->Fill(p3);
453 
454  m_nobha_hP3_z->Fill(p3);
455  m_nobha_hP3_y->Fill(p3);
456 
457  m_stt_theta->Fill(theta);
458  m_stt_phi->Fill(phiDegree);
459  p_stt_P3.push_back(p3);
460 
461  m_nobha_stt_theta->Fill(theta);
462  m_nobha_stt_phi->Fill(phiDegree);
463  p_nobha_stt_P3.push_back(p3);
464 
465 
466  if (trg_psnecl) {
467  m_hPt_psnecl->Fill(pt);
468  m_nobha_hPt_psnecl->Fill(pt);
469 
470  m_hP3_z_psnecl->Fill(p3); // for z bit
471  m_hP3_y_psnecl->Fill(p3); // for y bit
472 
473  m_nobha_hP3_z_psnecl->Fill(p3); // remove bhabha veto for z bit
474  m_nobha_hP3_y_psnecl->Fill(p3); // remove bhabha veto for y bit
475 
476  m_stt_phi_psnecl->Fill(phiDegree);
477  p_stt_P3_psnecl.push_back(p3);
478  m_stt_theta_psnecl->Fill(theta);
479 
480  m_nobha_stt_phi_psnecl->Fill(phiDegree);
481  p_nobha_stt_P3_psnecl.push_back(p3);
482  m_nobha_stt_theta_psnecl->Fill(theta);
483  }
484 
485  if (trg_psnecl && trg_ftdf) {
486  m_hPt_psnecl_ftdf->Fill(pt);
487  }
488  if (trg_psnecl && trg_itdt2) {
489  m_nobha_hPt_psnecl_ftdf->Fill(pt);
490  }
491 
492  if (trg_psnecl && trg_ftdz) {
493  m_hP3_z_psnecl_ftdf->Fill(p3);
494  }
495  if (trg_psnecl && trg_ftdy) {
496  m_hP3_y_psnecl_ftdf->Fill(p3);
497  }
498  if (trg_psnecl && trg_itdt3) {
499  m_nobha_hP3_z_psnecl_ftdf->Fill(p3);
500  }
501  if (trg_psnecl && trg_itdt4) {
502  m_nobha_hP3_y_psnecl_ftdf->Fill(p3);
503  }
504 
505  if (trg_psnecl && trg_stt) {
506  m_stt_phi_psnecl_ftdf->Fill(phiDegree);
507  p_stt_P3_psnecl_ftdf.push_back(p3);
508  m_stt_theta_psnecl_ftdf->Fill(theta);
509  }
510  if (trg_psnecl && trg_stt_nobha) {
511  m_nobha_stt_phi_psnecl_ftdf->Fill(phiDegree);
512  p_nobha_stt_P3_psnecl_ftdf.push_back(p3);
513  m_nobha_stt_theta_psnecl_ftdf->Fill(theta);
514  }
515 
516 
517  // the following is for fyo \deleta_phi
518  int njtrack = 0; // the j track of the m_Tracks
519  for (const auto& j_b2track : m_Tracks) {
520  if (nitrack >= njtrack) {
521  njtrack++;
522  continue;
523  }
524 
525  const Belle2::TrackFitResult* jfitresult = j_b2track.getTrackFitResultWithClosestMass(Const::pion);
526  if (!jfitresult) {
527  B2WARNING("No track fit result found.");
528  njtrack++;
529  continue;
530  }
531 
532  // require high NDF track
533  int jndf = jfitresult->getNDF();
534  if (jndf < 20) {
535  njtrack++;
536  continue;
537  }
538 
539  // IP tracks at barrel
540  if (fabs(jfitresult->getD0()) < 1.0 && fabs(jfitresult->getZ0()) < 1.0 && jfitresult->getHitPatternCDC().getLastLayer() > 50
541  && jfitresult->getHitPatternCDC().getFirstLayer() < 5) {
542  double jrk_phiDegree = jfitresult->getPhi() / Unit::deg;
543  double deltea_phi = fabs(phiDegree - jrk_phiDegree);
544  double dphi = deltea_phi;
545 
546  if (deltea_phi > 180)
547  dphi = 360 - deltea_phi;
548 
549  bool trg_fyo = m_trgSummary->testFtdl("fyo");
550  bool trg_fyo_nobha = (m_trgSummary->testInput("t2_1") || m_trgSummary->testInput("t2_2") || m_trgSummary->testInput("t2_3"))
551  &&
552  (m_trgSummary->testInput("ty_0") || m_trgSummary->testInput("ty_1") || m_trgSummary->testInput("ty_2")
553  || m_trgSummary->testInput("ty_3")) &&
554  m_trgSummary->testInput("cdc_open90");
555 
556  // cout<<"i = "<<nitrack <<" j= "<<njtrack <<" phiDegree = "<<phiDegree<<" jrk_phiDegree = "<<jrk_phiDegree<<" dphi "<<dphi<<endl;
557 
558  phi_fyo_dphi.push_back(dphi);
559  phi_nobha_fyo_dphi.push_back(dphi);
560 
561  if (trg_psnecl) {
562  phi_fyo_dphi_psnecl.push_back(dphi);
563  phi_nobha_fyo_dphi_psnecl.push_back(dphi);
564  }
565  if (trg_psnecl && trg_fyo) {
566  phi_fyo_dphi_psnecl_ftdf.push_back(dphi);
567  }
568  if (trg_psnecl && trg_fyo_nobha) {
569  phi_nobha_fyo_dphi_psnecl_ftdf.push_back(dphi);
570  }
571  }
572  njtrack++;
573  }
574  }
575  nitrack++;
576  }
577 
578 
579  // the largest cdc_open angle in an event for fyo bit
580  if (phi_fyo_dphi_psnecl_ftdf.size() != 0) {
581  auto max_it = std::max_element(phi_fyo_dphi_psnecl_ftdf.begin(), phi_fyo_dphi_psnecl_ftdf.end());
582  double max_value = *max_it;
583  m_fyo_dphi_psnecl_ftdf->Fill(max_value);
584  }
585  if (phi_fyo_dphi_psnecl.size() != 0) {
586  auto max_it = std::max_element(phi_fyo_dphi_psnecl.begin(), phi_fyo_dphi_psnecl.end());
587  double max_value = *max_it;
588  m_fyo_dphi_psnecl->Fill(max_value);
589  }
590  if (phi_fyo_dphi.size() != 0) {
591  auto max_it = std::max_element(phi_fyo_dphi.begin(), phi_fyo_dphi.end());
592  double max_value = *max_it;
593  m_fyo_dphi->Fill(max_value);
594  }
595 
596  //
597  if (phi_nobha_fyo_dphi_psnecl_ftdf.size() != 0) {
598  auto max_it = std::max_element(phi_nobha_fyo_dphi_psnecl_ftdf.begin(), phi_nobha_fyo_dphi_psnecl_ftdf.end());
599  double max_value = *max_it;
600  m_nobha_fyo_dphi_psnecl_ftdf->Fill(max_value);
601  }
602  if (phi_nobha_fyo_dphi_psnecl.size() != 0) {
603  auto max_it = std::max_element(phi_nobha_fyo_dphi_psnecl.begin(), phi_nobha_fyo_dphi_psnecl.end());
604  double max_value = *max_it;
605  m_nobha_fyo_dphi_psnecl->Fill(max_value);
606  }
607  if (phi_nobha_fyo_dphi.size() != 0) {
608  auto max_it = std::max_element(phi_nobha_fyo_dphi.begin(), phi_nobha_fyo_dphi.end());
609  double max_value = *max_it;
610  m_nobha_fyo_dphi->Fill(max_value);
611  }
612 
613 
614  // the largest momentum track p in an event for stt bit
615  if (p_stt_P3_psnecl_ftdf.size() != 0) {
616  auto max_it = std::max_element(p_stt_P3_psnecl_ftdf.begin(), p_stt_P3_psnecl_ftdf.end());
617  double max_value = *max_it;
618  m_stt_P3_psnecl_ftdf->Fill(max_value);
619  }
620  if (p_stt_P3_psnecl.size() != 0) {
621  auto max_it = std::max_element(p_stt_P3_psnecl.begin(), p_stt_P3_psnecl.end());
622  double max_value = *max_it;
623  m_stt_P3_psnecl->Fill(max_value);
624  }
625  if (p_stt_P3.size() != 0) {
626  auto max_it = std::max_element(p_stt_P3.begin(), p_stt_P3.end());
627  double max_value = *max_it;
628  m_stt_P3->Fill(max_value);
629  }
630 
631  //
632  if (p_nobha_stt_P3_psnecl_ftdf.size() != 0) {
633  auto max_it = std::max_element(p_nobha_stt_P3_psnecl_ftdf.begin(), p_nobha_stt_P3_psnecl_ftdf.end());
634  double max_value = *max_it;
635  m_nobha_stt_P3_psnecl_ftdf->Fill(max_value);
636  }
637  if (p_nobha_stt_P3_psnecl.size() != 0) {
638  auto max_it = std::max_element(p_nobha_stt_P3_psnecl.begin(), p_nobha_stt_P3_psnecl.end());
639  double max_value = *max_it;
640  m_nobha_stt_P3_psnecl->Fill(max_value);
641  }
642  if (p_nobha_stt_P3.size() != 0) {
643  auto max_it = std::max_element(p_nobha_stt_P3.begin(), p_nobha_stt_P3.end());
644  double max_value = *max_it;
645  m_nobha_stt_P3->Fill(max_value);
646  }
647 
648 
649 
650 
651 }
652 
653 
654 
655 void TRGEFFDQMModule::endRun()
656 {
657 }
658 
659 void TRGEFFDQMModule::terminate()
660 {
661 }
Belle2::HistoModule
HistoModule.h is supposed to be used instead of Module.h for the modules with histogram definitions t...
Definition: HistoModule.h:29
Belle2::Module::setDescription
void setDescription(const std::string &description)
Sets the description of the module.
Definition: Module.cc:214
Belle2::Module::setPropertyFlags
void setPropertyFlags(unsigned int propertyFlags)
Sets the flags for the module properties.
Definition: Module.cc:208
Belle2::Module::c_ParallelProcessingCertified
@ c_ParallelProcessingCertified
This module can be run in parallel processing mode safely (All I/O must be done through the data stor...
Definition: Module.h:80
Belle2::StoreAccessorBase::isOptional
bool isOptional(const std::string &name="")
Tell the DataStore about an optional input.
Definition: StoreAccessorBase.h:93
Belle2::TRGEFFDQMModule::initialize
void initialize() override
Initialize the Module.
Definition: TRGEFFDQMModule.cc:152
Belle2::TRGEFFDQMModule::m_hPhi_psnecl_ftdf
TH1F * m_hPhi_psnecl_ftdf
Histogram of cdc phi of IP tracks with ecl and f bit.
Definition: TRGEFFDQMModule.h:103
Belle2::TRGEFFDQMModule::~TRGEFFDQMModule
virtual ~TRGEFFDQMModule()
Destructor.
Definition: TRGEFFDQMModule.cc:37
Belle2::TRGEFFDQMModule::event
void event() override
Event processor.
Definition: TRGEFFDQMModule.cc:184
Belle2::TRGEFFDQMModule::m_trgSummary
StoreObjPtr< TRGSummary > m_trgSummary
Trigger summary.
Definition: TRGEFFDQMModule.h:93
Belle2::TRGEFFDQMModule::m_hPhi_psnecl
TH1F * m_hPhi_psnecl
Histogram of cdc phi of IP tracks with ecl bit.
Definition: TRGEFFDQMModule.h:102
Belle2::TRGEFFDQMModule::beginRun
void beginRun() override
Called when entering a new run.
Definition: TRGEFFDQMModule.cc:174
Belle2::TRGEFFDQMModule::m_ECLClusters
StoreArray< ECLCluster > m_ECLClusters
ECL Clusters.
Definition: TRGEFFDQMModule.h:95
Belle2::TRGEFFDQMModule::m_hPhi
TH1F * m_hPhi
Histogram of cdc phi of IP tracks.
Definition: TRGEFFDQMModule.h:101
Belle2::TRGEFFDQMModule::m_KLMClusters
StoreArray< KLMCluster > m_KLMClusters
KLM Clusters.
Definition: TRGEFFDQMModule.h:96
Belle2::TRGEFFDQMModule::m_RecoTracks
StoreArray< RecoTrack > m_RecoTracks
RecoTracks.
Definition: TRGEFFDQMModule.h:97
Belle2::TRGEFFDQMModule::m_Tracks
StoreArray< Track > m_Tracks
Tracks.
Definition: TRGEFFDQMModule.h:94
Belle2::TRGEFFDQMModule::m_TrgResult
StoreObjPtr< SoftwareTriggerResult > m_TrgResult
Store array for Trigger selection.
Definition: TRGEFFDQMModule.h:99
Belle2::TRGEFFDQMModule::defineHisto
void defineHisto() override
Histogram definitions.
Definition: TRGEFFDQMModule.cc:41
Belle2::Unit::deg
static const double deg
degree to radians
Definition: Unit.h:109
REG_MODULE
#define REG_MODULE(moduleName)
Register the given module (without 'Module' suffix) with the framework.
Definition: Module.h:650
Belle2::SoftwareTriggerCutResult::c_accept
@ c_accept
Accept this event.
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17