Belle II Software  release-05-01-25
PhysicsObjectsMiraBelleDstModule.cc
1 /**************************************************************************
2  * BASF2 (Belle Analysis Framework 2) *
3  * Copyright(C) 2020 Belle II Collaboration *
4  * *
5  * Author: The Belle II Collaboration *
6  * Contributors: Shun Watanuki *
7  * *
8  * This software is provided "as is" without any warranty. *
9  **************************************************************************/
10 
11 #include <dqm/modules/PhysicsObjectsMiraBelle/PhysicsObjectsMiraBelleDstModule.h>
12 #include <analysis/dataobjects/ParticleList.h>
13 #include <analysis/variables/ContinuumSuppressionVariables.h>
14 #include <analysis/variables/TrackVariables.h>
15 #include <analysis/utility/PCmsLabTransform.h>
16 #include <analysis/variables/Variables.h>
17 #include <framework/datastore/StoreObjPtr.h>
18 #include <framework/datastore/StoreArray.h>
19 #include <mdst/dataobjects/Track.h>
20 #include <mdst/dataobjects/TrackFitResult.h>
21 #include <mdst/dataobjects/KLMCluster.h>
22 #include <mdst/dataobjects/HitPatternCDC.h>
23 #include <mdst/dataobjects/HitPatternVXD.h>
24 #include <mdst/dataobjects/EventLevelTrackingInfo.h>
25 #include <mdst/dataobjects/PIDLikelihood.h>
26 #include <top/variables/TOPDigitVariables.h>
27 #include <arich/modules/arichDQM/ARICHDQMModule.h>
28 #include <arich/dataobjects/ARICHLikelihood.h>
29 #include <mdst/dataobjects/SoftwareTriggerResult.h>
30 #include <TDirectory.h>
31 #include <map>
32 
33 using namespace Belle2;
34 
35 REG_MODULE(PhysicsObjectsMiraBelleDst)
36 
37 PhysicsObjectsMiraBelleDstModule::PhysicsObjectsMiraBelleDstModule() : HistoModule()
38 {
39  setDescription("Monitor Physics Objects Quality");
40  setPropertyFlags(c_ParallelProcessingCertified);
41 
42  addParam("TriggerIdentifier", m_triggerIdentifier,
43  "Trigger identifier string used to select events for the histograms", std::string("software_trigger_cut&skim&accept_hadron"));
44  addParam("DstListName", m_dstListName, "Name of the D*+ particle list", std::string("D*+:physMiraBelle"));
45 }
46 
47 void PhysicsObjectsMiraBelleDstModule::defineHisto()
48 {
49  TDirectory* oldDir = gDirectory;
50  oldDir->mkdir("PhysicsObjectsMiraBelleDst")->cd();
51  // Mass distributions
52  m_h_D0_InvM = new TH1F("hist_D0_InvM", "Signal enhanced;m_h_D0_InvM;", 50, 1.81, 1.95);
53  m_h_delta_m = new TH1F("hist_delta_m", "Signal enhanced;delta_m;", 50, 0.14, 0.16);
54  // Soft pion
55  m_h_D0_softpi_PID_ALL_pion = new TH1F("hist_D0_softpi_PID_ALL_pion", "D0_softpi_PID_ALL_pion;D0_softpi_PID_ALL_pion;", 50, 0, 1);
56  m_h_D0_softpi_PID_SVD_pion = new TH1F("hist_D0_softpi_PID_SVD_pion", "D0_softpi_PID_SVD_pion;D0_softpi_PID_SVD_pion;", 50, 0, 1);
57  m_h_D0_softpi_PID_CDC_pion = new TH1F("hist_D0_softpi_PID_CDC_pion", "D0_softpi_PID_CDC_pion;D0_softpi_PID_CDC_pion;", 50, 0, 1);
58  m_h_D0_softpi_PID_TOP_pion = new TH1F("hist_D0_softpi_PID_TOP_pion", "D0_softpi_PID_TOP_pion;D0_softpi_PID_TOP_pion;", 50, 0, 1);
59  m_h_D0_softpi_PID_ARICH_pion = new TH1F("hist_D0_softpi_PID_ARICH_pion", "D0_softpi_PID_ARICH_pion;D0_softpi_PID_ARICH_pion;", 50,
60  0,
61  1);
62  m_h_D0_softpi_PID_ECL_pion = new TH1F("hist_D0_softpi_PID_ECL_pion", "D0_softpi_PID_ECL_pion;D0_softpi_PID_ECL_pion;", 50, 0, 1);
63  m_h_D0_softpi_PID_KLM_pion = new TH1F("hist_D0_softpi_PID_KLM_pion", "D0_softpi_PID_KLM_pion;D0_softpi_PID_KLM_pion;", 50, 0, 1);
64  m_h_sideband_D0_softpi_PID_ALL_pion = new TH1F("hist_sideband_D0_softpi_PID_ALL_pion",
65  "D0_softpi_PID_ALL_pion;D0_softpi_PID_ALL_pion;", 50, 0, 1);
66  m_h_sideband_D0_softpi_PID_SVD_pion = new TH1F("hist_sideband_D0_softpi_PID_SVD_pion",
67  "D0_softpi_PID_SVD_pion;D0_softpi_PID_SVD_pion;", 50, 0, 1);
68  m_h_sideband_D0_softpi_PID_CDC_pion = new TH1F("hist_sideband_D0_softpi_PID_CDC_pion",
69  "D0_softpi_PID_CDC_pion;D0_softpi_PID_CDC_pion;", 50, 0, 1);
70  m_h_sideband_D0_softpi_PID_TOP_pion = new TH1F("hist_sideband_D0_softpi_PID_TOP_pion",
71  "D0_softpi_PID_TOP_pion;D0_softpi_PID_TOP_pion;", 50, 0, 1);
72  m_h_sideband_D0_softpi_PID_ARICH_pion = new TH1F("hist_sideband_D0_softpi_PID_ARICH_pion",
73  "D0_softpi_PID_ARICH_pion;D0_softpi_PID_ARICH_pion;", 50, 0, 1);
74  m_h_sideband_D0_softpi_PID_ECL_pion = new TH1F("hist_sideband_D0_softpi_PID_ECL_pion",
75  "D0_softpi_PID_ECL_pion;D0_softpi_PID_ECL_pion;", 50, 0, 1);
76  m_h_sideband_D0_softpi_PID_KLM_pion = new TH1F("hist_sideband_D0_softpi_PID_KLM_pion",
77  "D0_softpi_PID_KLM_pion;D0_softpi_PID_KLM_pion;", 50, 0, 1);
78  // pi from D0
79  m_h_D0_pi_PID_ALL_pion = new TH1F("hist_D0_pi_PID_ALL_pion", "D0_pi_PID_ALL_pion;D0_pi_PID_ALL_pion;", 50, 0, 1);
80  m_h_D0_pi_PID_SVD_pion = new TH1F("hist_D0_pi_PID_SVD_pion", "D0_pi_PID_SVD_pion;D0_pi_PID_SVD_pion;", 50, 0, 1);
81  m_h_D0_pi_PID_CDC_pion = new TH1F("hist_D0_pi_PID_CDC_pion", "D0_pi_PID_CDC_pion;D0_pi_PID_CDC_pion;", 50, 0, 1);
82  m_h_D0_pi_PID_TOP_pion = new TH1F("hist_D0_pi_PID_TOP_pion", "D0_pi_PID_TOP_pion;D0_pi_PID_TOP_pion;", 50, 0, 1);
83  m_h_D0_pi_PID_ARICH_pion = new TH1F("hist_D0_pi_PID_ARICH_pion", "D0_pi_PID_ARICH_pion;D0_pi_PID_ARICH_pion;", 50, 0, 1);
84  m_h_D0_pi_PID_ECL_pion = new TH1F("hist_D0_pi_PID_ECL_pion", "D0_pi_PID_ECL_pion;D0_pi_PID_ECL_pion;", 50, 0, 1);
85  m_h_D0_pi_PID_KLM_pion = new TH1F("hist_D0_pi_PID_KLM_pion", "D0_pi_PID_KLM_pion;D0_pi_PID_KLM_pion;", 50, 0, 1);
86  m_h_sideband_D0_pi_PID_ALL_pion = new TH1F("hist_sideband_D0_pi_PID_ALL_pion", "D0_pi_PID_ALL_pion;D0_pi_PID_ALL_pion;", 50, 0, 1);
87  m_h_sideband_D0_pi_PID_SVD_pion = new TH1F("hist_sideband_D0_pi_PID_SVD_pion", "D0_pi_PID_SVD_pion;D0_pi_PID_SVD_pion;", 50, 0, 1);
88  m_h_sideband_D0_pi_PID_CDC_pion = new TH1F("hist_sideband_D0_pi_PID_CDC_pion", "D0_pi_PID_CDC_pion;D0_pi_PID_CDC_pion;", 50, 0, 1);
89  m_h_sideband_D0_pi_PID_TOP_pion = new TH1F("hist_sideband_D0_pi_PID_TOP_pion", "D0_pi_PID_TOP_pion;D0_pi_PID_TOP_pion;", 50, 0, 1);
90  m_h_sideband_D0_pi_PID_ARICH_pion = new TH1F("hist_sideband_D0_pi_PID_ARICH_pion", "D0_pi_PID_ARICH_pion;D0_pi_PID_ARICH_pion;", 50,
91  0, 1);
92  m_h_sideband_D0_pi_PID_ECL_pion = new TH1F("hist_sideband_D0_pi_PID_ECL_pion", "D0_pi_PID_ECL_pion;D0_pi_PID_ECL_pion;", 50, 0, 1);
93  m_h_sideband_D0_pi_PID_KLM_pion = new TH1F("hist_sideband_D0_pi_PID_KLM_pion", "D0_pi_PID_KLM_pion;D0_pi_PID_KLM_pion;", 50, 0, 1);
94  // K from D0
95  m_h_D0_K_PID_ALL_kaon = new TH1F("hist_D0_K_PID_ALL_kaon", "D0_K_PID_ALL_kaon;D0_K_PID_ALL_kaon;", 50, 0, 1);
96  m_h_D0_K_PID_SVD_kaon = new TH1F("hist_D0_K_PID_SVD_kaon", "D0_K_PID_SVD_kaon;D0_K_PID_SVD_kaon;", 50, 0, 1);
97  m_h_D0_K_PID_CDC_kaon = new TH1F("hist_D0_K_PID_CDC_kaon", "D0_K_PID_CDC_kaon;D0_K_PID_CDC_kaon;", 50, 0, 1);
98  m_h_D0_K_PID_TOP_kaon = new TH1F("hist_D0_K_PID_TOP_kaon", "D0_K_PID_TOP_kaon;D0_K_PID_TOP_kaon;", 50, 0, 1);
99  m_h_D0_K_PID_ARICH_kaon = new TH1F("hist_D0_K_PID_ARICH_kaon", "D0_K_PID_ARICH_kaon;D0_K_PID_ARICH_kaon;", 50, 0, 1);
100  m_h_D0_K_PID_ECL_kaon = new TH1F("hist_D0_K_PID_ECL_kaon", "D0_K_PID_ECL_kaon;D0_K_PID_ECL_kaon;", 50, 0, 1);
101  m_h_D0_K_PID_KLM_kaon = new TH1F("hist_D0_K_PID_KLM_kaon", "D0_K_PID_KLM_kaon;D0_K_PID_KLM_kaon;", 50, 0, 1);
102  m_h_sideband_D0_K_PID_ALL_kaon = new TH1F("hist_sideband_D0_K_PID_ALL_kaon", "D0_K_PID_ALL_kaon;D0_K_PID_ALL_kaon;", 50, 0, 1);
103  m_h_sideband_D0_K_PID_SVD_kaon = new TH1F("hist_sideband_D0_K_PID_SVD_kaon", "D0_K_PID_SVD_kaon;D0_K_PID_SVD_kaon;", 50, 0, 1);
104  m_h_sideband_D0_K_PID_CDC_kaon = new TH1F("hist_sideband_D0_K_PID_CDC_kaon", "D0_K_PID_CDC_kaon;D0_K_PID_CDC_kaon;", 50, 0, 1);
105  m_h_sideband_D0_K_PID_TOP_kaon = new TH1F("hist_sideband_D0_K_PID_TOP_kaon", "D0_K_PID_TOP_kaon;D0_K_PID_TOP_kaon;", 50, 0, 1);
106  m_h_sideband_D0_K_PID_ARICH_kaon = new TH1F("hist_sideband_D0_K_PID_ARICH_kaon", "D0_K_PID_ARICH_kaon;D0_K_PID_ARICH_kaon;", 50, 0,
107  1);
108  m_h_sideband_D0_K_PID_ECL_kaon = new TH1F("hist_sideband_D0_K_PID_ECL_kaon", "D0_K_PID_ECL_kaon;D0_K_PID_ECL_kaon;", 50, 0, 1);
109  m_h_sideband_D0_K_PID_KLM_kaon = new TH1F("hist_sideband_D0_K_PID_KLM_kaon", "D0_K_PID_KLM_kaon;D0_K_PID_KLM_kaon;", 50, 0, 1);
110 
111  oldDir->cd();
112 }
113 
114 
115 void PhysicsObjectsMiraBelleDstModule::initialize()
116 {
117  REG_HISTOGRAM
118 
119  StoreObjPtr<SoftwareTriggerResult> result;
120  result.isOptional();
121 }
122 
123 void PhysicsObjectsMiraBelleDstModule::beginRun()
124 {
125  m_h_D0_InvM->Reset();
126  m_h_delta_m->Reset();
127  m_h_D0_softpi_PID_ALL_pion->Reset();
128  m_h_D0_softpi_PID_SVD_pion->Reset();
129  m_h_D0_softpi_PID_CDC_pion->Reset();
130  m_h_D0_softpi_PID_TOP_pion->Reset();
131  m_h_D0_softpi_PID_ARICH_pion->Reset();
132  m_h_D0_softpi_PID_ECL_pion->Reset();
133  m_h_D0_softpi_PID_KLM_pion->Reset();
134  m_h_D0_pi_PID_ALL_pion->Reset();
135  m_h_D0_pi_PID_SVD_pion->Reset();
136  m_h_D0_pi_PID_CDC_pion->Reset();
137  m_h_D0_pi_PID_TOP_pion->Reset();
138  m_h_D0_pi_PID_ARICH_pion->Reset();
139  m_h_D0_pi_PID_ECL_pion->Reset();
140  m_h_D0_pi_PID_KLM_pion->Reset();
141  m_h_D0_K_PID_ALL_kaon->Reset();
142  m_h_D0_K_PID_SVD_kaon->Reset();
143  m_h_D0_K_PID_CDC_kaon->Reset();
144  m_h_D0_K_PID_TOP_kaon->Reset();
145  m_h_D0_K_PID_ARICH_kaon->Reset();
146  m_h_D0_K_PID_ECL_kaon->Reset();
147  m_h_D0_K_PID_KLM_kaon->Reset();
148  m_h_sideband_D0_softpi_PID_ALL_pion->Reset();
149  m_h_sideband_D0_softpi_PID_SVD_pion->Reset();
150  m_h_sideband_D0_softpi_PID_CDC_pion->Reset();
151  m_h_sideband_D0_softpi_PID_TOP_pion->Reset();
152  m_h_sideband_D0_softpi_PID_ARICH_pion->Reset();
153  m_h_sideband_D0_softpi_PID_ECL_pion->Reset();
154  m_h_sideband_D0_softpi_PID_KLM_pion->Reset();
155  m_h_sideband_D0_pi_PID_ALL_pion->Reset();
156  m_h_sideband_D0_pi_PID_SVD_pion->Reset();
157  m_h_sideband_D0_pi_PID_CDC_pion->Reset();
158  m_h_sideband_D0_pi_PID_TOP_pion->Reset();
159  m_h_sideband_D0_pi_PID_ARICH_pion->Reset();
160  m_h_sideband_D0_pi_PID_ECL_pion->Reset();
161  m_h_sideband_D0_pi_PID_KLM_pion->Reset();
162  m_h_sideband_D0_K_PID_ALL_kaon->Reset();
163  m_h_sideband_D0_K_PID_SVD_kaon->Reset();
164  m_h_sideband_D0_K_PID_CDC_kaon->Reset();
165  m_h_sideband_D0_K_PID_TOP_kaon->Reset();
166  m_h_sideband_D0_K_PID_ARICH_kaon->Reset();
167  m_h_sideband_D0_K_PID_ECL_kaon->Reset();
168  m_h_sideband_D0_K_PID_KLM_kaon->Reset();
169 }
170 
171 void PhysicsObjectsMiraBelleDstModule::event()
172 {
173 
174  StoreObjPtr<SoftwareTriggerResult> result;
175  if (!result.isValid()) {
176  B2WARNING("SoftwareTriggerResult object not available but needed to select events for the histograms.");
177  return;
178  }
179 
180  const std::map<std::string, int>& results = result->getResults();
181  if (results.find(m_triggerIdentifier) == results.end()) {
182  B2WARNING("PhysicsObjectsMiraBelleDst: Can't find trigger identifier: " << m_triggerIdentifier);
183  return;
184  }
185 
186  // apply software trigger
187  const bool accepted = (result->getResult(m_triggerIdentifier) == SoftwareTriggerCutResult::c_accept);
188  if (accepted == false) return;
189 
190  // get D* candidates
191  StoreObjPtr<ParticleList> dstParticles(m_dstListName);
192 
193  for (unsigned int i = 0; i < dstParticles->getListSize(); i++) {
194  const Particle* dst = dstParticles->getParticle(i);
195  const Particle* d0 = dst->getDaughter(0);
196  float dst_mass = dst->getMass();
197  float d0_mass = d0->getMass();
198  float delta_m = dst_mass - d0_mass;
199  if (delta_m < 0.14 || 0.16 < delta_m || d0_mass < 1.81 || 1.95 < d0_mass) continue;
200  // True if the event is in signal region
201  bool isSignal_D0_InvM(false), isSignal_delta_m(false);
202  // D0 and D*+ mass
203  if (0.143 < delta_m && delta_m < 0.147) {
204  m_h_D0_InvM->Fill(d0_mass);
205  isSignal_D0_InvM = true;
206  }
207  if (1.83 < d0_mass && d0_mass < 1.89) {
208  m_h_delta_m->Fill(delta_m);
209  isSignal_delta_m = true;
210  }
211  const PIDLikelihood* pid_K = d0->getDaughter(0)->getPIDLikelihood();
212  const PIDLikelihood* pid_Pi = d0->getDaughter(1)->getPIDLikelihood();
213  const PIDLikelihood* pid_softPi = dst->getDaughter(1)->getPIDLikelihood();
214 
215  if (isSignal_D0_InvM && isSignal_delta_m) {
216  // Signal region
217  // PID of K
218  m_h_D0_K_PID_ALL_kaon->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion, 1.));
219  if (pid_K->isAvailable(Const::SVD)) m_h_D0_K_PID_SVD_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
220  1., Const::SVD));
221  if (pid_K->isAvailable(Const::CDC)) m_h_D0_K_PID_CDC_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
222  1., Const::CDC));
223  if (pid_K->isAvailable(Const::TOP)) m_h_D0_K_PID_TOP_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
224  1., Const::TOP));
225  if (pid_K->isAvailable(Const::ARICH)) m_h_D0_K_PID_ARICH_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
226  1., Const::ARICH));
227  if (pid_K->isAvailable(Const::ECL)) m_h_D0_K_PID_ECL_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
228  1., Const::ECL));
229  if (pid_K->isAvailable(Const::KLM)) m_h_D0_K_PID_KLM_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
230  1., Const::KLM));
231  // PID of pi
232  m_h_D0_pi_PID_ALL_pion->Fill(pid_Pi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
233  if (pid_Pi->isAvailable(Const::SVD)) m_h_D0_pi_PID_SVD_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
234  Belle2::Const::kaon, 1., Const::SVD));
235  if (pid_Pi->isAvailable(Const::CDC)) m_h_D0_pi_PID_CDC_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
236  Belle2::Const::kaon, 1., Const::CDC));
237  if (pid_Pi->isAvailable(Const::TOP)) m_h_D0_pi_PID_TOP_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
238  Belle2::Const::kaon, 1., Const::TOP));
239  if (pid_Pi->isAvailable(Const::ARICH)) m_h_D0_pi_PID_ARICH_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
240  Belle2::Const::kaon, 1., Const::ARICH));
241  if (pid_Pi->isAvailable(Const::ECL)) m_h_D0_pi_PID_ECL_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
242  Belle2::Const::kaon, 1., Const::ECL));
243  if (pid_Pi->isAvailable(Const::KLM)) m_h_D0_pi_PID_KLM_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
244  Belle2::Const::kaon, 1., Const::KLM));
245  // PID of soft pi
246  m_h_D0_softpi_PID_ALL_pion->Fill(pid_softPi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
247  if (pid_softPi->isAvailable(Const::SVD)) m_h_D0_softpi_PID_SVD_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
248  Belle2::Const::kaon, 1., Const::SVD));
249  if (pid_softPi->isAvailable(Const::CDC)) m_h_D0_softpi_PID_CDC_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
250  Belle2::Const::kaon, 1., Const::CDC));
251  if (pid_softPi->isAvailable(Const::TOP)) m_h_D0_softpi_PID_TOP_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
252  Belle2::Const::kaon, 1., Const::TOP));
253  if (pid_softPi->isAvailable(Const::ARICH)) m_h_D0_softpi_PID_ARICH_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
254  Belle2::Const::kaon, 1., Const::ARICH));
255  if (pid_softPi->isAvailable(Const::ECL)) m_h_D0_softpi_PID_ECL_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
256  Belle2::Const::kaon, 1., Const::ECL));
257  if (pid_softPi->isAvailable(Const::KLM)) m_h_D0_softpi_PID_KLM_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
258  Belle2::Const::kaon, 1., Const::KLM));
259  } else {
260  // Sideband region for BG subtraction
261  // PID of K
262  m_h_sideband_D0_K_PID_ALL_kaon->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion, 1.));
263  if (pid_K->isAvailable(Const::SVD)) m_h_sideband_D0_K_PID_SVD_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
264  Belle2::Const::pion, 1., Const::SVD));
265  if (pid_K->isAvailable(Const::CDC)) m_h_sideband_D0_K_PID_CDC_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
266  Belle2::Const::pion, 1., Const::CDC));
267  if (pid_K->isAvailable(Const::TOP)) m_h_sideband_D0_K_PID_TOP_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
268  Belle2::Const::pion, 1., Const::TOP));
269  if (pid_K->isAvailable(Const::ARICH)) m_h_sideband_D0_K_PID_ARICH_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
270  Belle2::Const::pion, 1., Const::ARICH));
271  if (pid_K->isAvailable(Const::ECL)) m_h_sideband_D0_K_PID_ECL_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
272  Belle2::Const::pion, 1., Const::ECL));
273  if (pid_K->isAvailable(Const::KLM)) m_h_sideband_D0_K_PID_KLM_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
274  Belle2::Const::pion, 1., Const::KLM));
275  // PID of pi
276  m_h_sideband_D0_pi_PID_ALL_pion->Fill(pid_Pi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
277  if (pid_Pi->isAvailable(Const::SVD)) m_h_sideband_D0_pi_PID_SVD_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
278  Belle2::Const::kaon, 1., Const::SVD));
279  if (pid_Pi->isAvailable(Const::CDC)) m_h_sideband_D0_pi_PID_CDC_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
280  Belle2::Const::kaon, 1., Const::CDC));
281  if (pid_Pi->isAvailable(Const::TOP)) m_h_sideband_D0_pi_PID_TOP_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
282  Belle2::Const::kaon, 1., Const::TOP));
283  if (pid_Pi->isAvailable(Const::ARICH)) m_h_sideband_D0_pi_PID_ARICH_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
284  Belle2::Const::kaon, 1., Const::ARICH));
285  if (pid_Pi->isAvailable(Const::ECL)) m_h_sideband_D0_pi_PID_ECL_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
286  Belle2::Const::kaon, 1., Const::ECL));
287  if (pid_Pi->isAvailable(Const::KLM)) m_h_sideband_D0_pi_PID_KLM_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
288  Belle2::Const::kaon, 1., Const::KLM));
289  // PID of soft pi
290  m_h_sideband_D0_softpi_PID_ALL_pion->Fill(pid_softPi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
291  if (pid_softPi->isAvailable(Const::SVD)) m_h_sideband_D0_softpi_PID_SVD_pion ->Fill(pid_softPi->getProbability(
292  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::SVD));
293  if (pid_softPi->isAvailable(Const::CDC)) m_h_sideband_D0_softpi_PID_CDC_pion ->Fill(pid_softPi->getProbability(
294  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::CDC));
295  if (pid_softPi->isAvailable(Const::TOP)) m_h_sideband_D0_softpi_PID_TOP_pion ->Fill(pid_softPi->getProbability(
296  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::TOP));
297  if (pid_softPi->isAvailable(Const::ARICH)) m_h_sideband_D0_softpi_PID_ARICH_pion ->Fill(pid_softPi->getProbability(
298  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::ARICH));
299  if (pid_softPi->isAvailable(Const::ECL)) m_h_sideband_D0_softpi_PID_ECL_pion ->Fill(pid_softPi->getProbability(
300  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::ECL));
301  if (pid_softPi->isAvailable(Const::KLM)) m_h_sideband_D0_softpi_PID_KLM_pion ->Fill(pid_softPi->getProbability(
302  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::KLM));
303  }
304  }
305 }
306 
307 // void PhysicsObjectsMiraBelleDstModule::getPIDInformationMiraBelle(Particle* part, float pid[3][7])
308 // {
309 // PIDLikelihood* I_like_coffee = part->getPIDLikelihood();
310 // pid[0][0] = I_like_coffee->getProbability();
311 // }
312 
313 void PhysicsObjectsMiraBelleDstModule::endRun()
314 {
315 }
316 
317 void PhysicsObjectsMiraBelleDstModule::terminate()
318 {
319 }
320 
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_pi_PID_ECL_pion
TH1F * m_h_sideband_D0_pi_PID_ECL_pion
PID of ECL for pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:159
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_softpi_PID_ARICH_pion
TH1F * m_h_D0_softpi_PID_ARICH_pion
PID of ARICH for soft pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:72
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_K_PID_ALL_kaon
TH1F * m_h_D0_K_PID_ALL_kaon
PID of all detectors for kaon.
Definition: PhysicsObjectsMiraBelleDstModule.h:102
Belle2::PhysicsObjectsMiraBelleDst2Module::m_triggerIdentifier
std::string m_triggerIdentifier
Trigger identifier string used to select events for the histograms.
Definition: PhysicsObjectsMiraBelleDst2Module.h:57
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_softpi_PID_ALL_pion
TH1F * m_h_D0_softpi_PID_ALL_pion
PID of all detectors for soft pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:60
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_softpi_PID_TOP_pion
TH1F * m_h_D0_softpi_PID_TOP_pion
PID of TOP for soft pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:69
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_pi_PID_SVD_pion
TH1F * m_h_D0_pi_PID_SVD_pion
PID of SVD for pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:84
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_pi_PID_CDC_pion
TH1F * m_h_sideband_D0_pi_PID_CDC_pion
PID of CDC for pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:150
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_K_PID_ECL_kaon
TH1F * m_h_D0_K_PID_ECL_kaon
PID of ECL for kaon.
Definition: PhysicsObjectsMiraBelleDstModule.h:117
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_pi_PID_ECL_pion
TH1F * m_h_D0_pi_PID_ECL_pion
PID of ECL for pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:96
Belle2::Module::setDescription
void setDescription(const std::string &description)
Sets the description of the module.
Definition: Module.cc:216
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_K_PID_ARICH_kaon
TH1F * m_h_sideband_D0_K_PID_ARICH_kaon
PID of ARICH for kaon in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:177
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_K_PID_ECL_kaon
TH1F * m_h_sideband_D0_K_PID_ECL_kaon
PID of ECL for kaon in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:180
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_pi_PID_CDC_pion
TH1F * m_h_D0_pi_PID_CDC_pion
PID of CDC for pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:87
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_K_PID_TOP_kaon
TH1F * m_h_D0_K_PID_TOP_kaon
PID of TOP for kaon.
Definition: PhysicsObjectsMiraBelleDstModule.h:111
REG_MODULE
#define REG_MODULE(moduleName)
Register the given module (without 'Module' suffix) with the framework.
Definition: Module.h:652
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:82
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_pi_PID_KLM_pion
TH1F * m_h_sideband_D0_pi_PID_KLM_pion
PID of KLM for pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:162
Belle2::PhysicsObjectsMiraBelleDstModule::event
void event() override
Function to process event record.
Definition: PhysicsObjectsMiraBelleDstModule.cc:171
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_softpi_PID_KLM_pion
TH1F * m_h_D0_softpi_PID_KLM_pion
PID of KLM for soft pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:78
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_K_PID_CDC_kaon
TH1F * m_h_sideband_D0_K_PID_CDC_kaon
PID of CDC for kaon in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:171
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_K_PID_SVD_kaon
TH1F * m_h_D0_K_PID_SVD_kaon
PID of SVD for kaon.
Definition: PhysicsObjectsMiraBelleDstModule.h:105
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_K_PID_CDC_kaon
TH1F * m_h_D0_K_PID_CDC_kaon
PID of CDC for kaon.
Definition: PhysicsObjectsMiraBelleDstModule.h:108
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_softpi_PID_SVD_pion
TH1F * m_h_D0_softpi_PID_SVD_pion
PID of SVD for soft pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:63
Belle2::PIDLikelihood
Class to collect log likelihoods from TOP, ARICH, dEdx, ECL and KLM aimed for output to mdst includes...
Definition: PIDLikelihood.h:37
Belle2::Particle::getDaughter
const Particle * getDaughter(unsigned i) const
Returns a pointer to the i-th daughter particle.
Definition: Particle.cc:595
Belle2::PhysicsObjectsMiraBelleDstModule::m_dstListName
std::string m_dstListName
Name of the mu+ particle list.
Definition: PhysicsObjectsMiraBelleDstModule.h:189
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_softpi_PID_TOP_pion
TH1F * m_h_sideband_D0_softpi_PID_TOP_pion
PID of TOP for soft pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:132
Belle2::PIDLikelihood::isAvailable
bool isAvailable(Const::PIDDetectorSet set) const
Check whether PID information from a given set of detectors is available.
Definition: PIDLikelihood.h:61
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_softpi_PID_CDC_pion
TH1F * m_h_sideband_D0_softpi_PID_CDC_pion
PID of CDC for soft pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:129
Belle2::SoftwareTriggerCutResult::c_accept
@ c_accept
Accept this event.
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_softpi_PID_SVD_pion
TH1F * m_h_sideband_D0_softpi_PID_SVD_pion
PID of SVD for soft pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:126
Belle2::Const::kaon
static const ChargedStable kaon
charged kaon particle
Definition: Const.h:536
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_InvM
TH1F * m_h_D0_InvM
D0 invariant mass.
Definition: PhysicsObjectsMiraBelleDstModule.h:54
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_softpi_PID_ALL_pion
TH1F * m_h_sideband_D0_softpi_PID_ALL_pion
PID of all detectors for soft pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:123
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_pi_PID_ARICH_pion
TH1F * m_h_sideband_D0_pi_PID_ARICH_pion
PID of ARICH for pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:156
Belle2::Module::setPropertyFlags
void setPropertyFlags(unsigned int propertyFlags)
Sets the flags for the module properties.
Definition: Module.cc:210
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_softpi_PID_ECL_pion
TH1F * m_h_D0_softpi_PID_ECL_pion
PID of ECL for soft pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:75
Belle2::Const::pion
static const ChargedStable pion
charged pion particle
Definition: Const.h:535
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_pi_PID_ALL_pion
TH1F * m_h_D0_pi_PID_ALL_pion
PID of all for pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:81
Belle2::Particle::getPIDLikelihood
const PIDLikelihood * getPIDLikelihood() const
Returns the pointer to the PIDLikelihood object that is related to the Track, which was used to creat...
Definition: Particle.cc:796
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_pi_PID_ARICH_pion
TH1F * m_h_D0_pi_PID_ARICH_pion
PID of ARICH for pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:93
Belle2::PhysicsObjectsMiraBelleDstModule::beginRun
void beginRun() override
Function to process begin_run record.
Definition: PhysicsObjectsMiraBelleDstModule.cc:123
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_pi_PID_ALL_pion
TH1F * m_h_sideband_D0_pi_PID_ALL_pion
PID of all detectors for pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:144
Belle2::Particle::getMass
float getMass() const
Returns invariant mass (= nominal for FS particles)
Definition: Particle.h:433
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:19
Belle2::StoreObjPtr
Type-safe access to single objects in the data store.
Definition: ParticleList.h:33
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_softpi_PID_ECL_pion
TH1F * m_h_sideband_D0_softpi_PID_ECL_pion
PID of ECL for soft pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:138
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_pi_PID_SVD_pion
TH1F * m_h_sideband_D0_pi_PID_SVD_pion
PID of SVD for pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:147
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_pi_PID_KLM_pion
TH1F * m_h_D0_pi_PID_KLM_pion
PID of KLM for pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:99
Belle2::PhysicsObjectsMiraBelleDstModule::terminate
void terminate() override
Function to terminate module.
Definition: PhysicsObjectsMiraBelleDstModule.cc:317
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_K_PID_TOP_kaon
TH1F * m_h_sideband_D0_K_PID_TOP_kaon
PID of TOP for kaon in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:174
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_K_PID_ARICH_kaon
TH1F * m_h_D0_K_PID_ARICH_kaon
PID of ARICH for kaon.
Definition: PhysicsObjectsMiraBelleDstModule.h:114
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_K_PID_SVD_kaon
TH1F * m_h_sideband_D0_K_PID_SVD_kaon
PID of SVD for kaon in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:168
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_K_PID_KLM_kaon
TH1F * m_h_D0_K_PID_KLM_kaon
PID of KLM for kaon.
Definition: PhysicsObjectsMiraBelleDstModule.h:120
Belle2::PhysicsObjectsMiraBelleDst2Module::m_dstListName
std::string m_dstListName
Name of the mu+ particle list.
Definition: PhysicsObjectsMiraBelleDst2Module.h:60
Belle2::PhysicsObjectsMiraBelleDstModule::initialize
void initialize() override
Function for dynamic initialization of module.
Definition: PhysicsObjectsMiraBelleDstModule.cc:115
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_softpi_PID_CDC_pion
TH1F * m_h_D0_softpi_PID_CDC_pion
PID of CDC for soft pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:66
Belle2::PhysicsObjectsMiraBelleDstModule::defineHisto
void defineHisto() override
Function to define histograms.
Definition: PhysicsObjectsMiraBelleDstModule.cc:47
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_delta_m
TH1F * m_h_delta_m
Mass difference b/w D*+ and D0.
Definition: PhysicsObjectsMiraBelleDstModule.h:57
Belle2::Particle
Class to store reconstructed particles.
Definition: Particle.h:77
Belle2::Module::addParam
void addParam(const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module.
Definition: Module.h:562
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_pi_PID_TOP_pion
TH1F * m_h_sideband_D0_pi_PID_TOP_pion
PID of TOP for pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:153
Belle2::PhysicsObjectsMiraBelleDstModule::m_triggerIdentifier
std::string m_triggerIdentifier
Trigger identifier string used to select events for the histograms.
Definition: PhysicsObjectsMiraBelleDstModule.h:186
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_softpi_PID_ARICH_pion
TH1F * m_h_sideband_D0_softpi_PID_ARICH_pion
PID of ARICH for soft pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:135
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_D0_pi_PID_TOP_pion
TH1F * m_h_D0_pi_PID_TOP_pion
PID of TOP for pion.
Definition: PhysicsObjectsMiraBelleDstModule.h:90
Belle2::PhysicsObjectsMiraBelleDstModule::endRun
void endRun() override
Function to process end_run record.
Definition: PhysicsObjectsMiraBelleDstModule.cc:313
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_softpi_PID_KLM_pion
TH1F * m_h_sideband_D0_softpi_PID_KLM_pion
PID of KLM for soft pion in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:141
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::PIDLikelihood::getProbability
double getProbability(const Const::ChargedStable &p1, const Const::ChargedStable &p2, Const::PIDDetectorSet set=Const::PIDDetectorSet::set()) const
Return combined likelihood probability for a particle being p1 and not p2, assuming equal prior proba...
Definition: PIDLikelihood.h:94
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_K_PID_KLM_kaon
TH1F * m_h_sideband_D0_K_PID_KLM_kaon
PID of KLM for kaon in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:183
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_sideband_D0_K_PID_ALL_kaon
TH1F * m_h_sideband_D0_K_PID_ALL_kaon
PID of all detectors for kaon in sideband.
Definition: PhysicsObjectsMiraBelleDstModule.h:165