Belle II Software  release-08-01-10
PhysicsObjectsMiraBelleDstModule.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 #include <dqm/modules/PhysicsObjectsMiraBelle/PhysicsObjectsMiraBelleDstModule.h>
10 #include <analysis/dataobjects/ParticleList.h>
11 #include <analysis/variables/ContinuumSuppressionVariables.h>
12 #include <analysis/variables/TrackVariables.h>
13 #include <analysis/utility/PCmsLabTransform.h>
14 #include <analysis/variables/Variables.h>
15 #include <framework/datastore/StoreObjPtr.h>
16 #include <framework/datastore/StoreArray.h>
17 #include <mdst/dataobjects/Track.h>
18 #include <mdst/dataobjects/TrackFitResult.h>
19 #include <mdst/dataobjects/KLMCluster.h>
20 #include <mdst/dataobjects/HitPatternCDC.h>
21 #include <mdst/dataobjects/HitPatternVXD.h>
22 #include <mdst/dataobjects/EventLevelTrackingInfo.h>
23 #include <mdst/dataobjects/PIDLikelihood.h>
24 #include <top/variables/TOPDigitVariables.h>
25 #include <arich/modules/arichDQM/ARICHDQMModule.h>
26 #include <arich/dataobjects/ARICHLikelihood.h>
27 #include <mdst/dataobjects/SoftwareTriggerResult.h>
28 #include <TDirectory.h>
29 #include <map>
30 
31 using namespace Belle2;
32 
33 REG_MODULE(PhysicsObjectsMiraBelleDst);
34 
35 PhysicsObjectsMiraBelleDstModule::PhysicsObjectsMiraBelleDstModule() : HistoModule()
36 {
37  setDescription("Monitor Physics Objects Quality");
38  setPropertyFlags(c_ParallelProcessingCertified);
39 
40  addParam("TriggerIdentifier", m_triggerIdentifier,
41  "Trigger identifier string used to select events for the histograms", std::string("software_trigger_cut&skim&accept_hadron"));
42  addParam("DstListName", m_dstListName, "Name of the D*+ particle list", std::string("D*+:physMiraBelle"));
43 }
44 
45 void PhysicsObjectsMiraBelleDstModule::defineHisto()
46 {
47  TDirectory* oldDir = gDirectory;
48  oldDir->mkdir("PhysicsObjectsMiraBelleDst");
49  oldDir->cd("PhysicsObjectsMiraBelleDst");
50 
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 (pid_K != NULL && pid_Pi != NULL && pid_softPi != NULL) {
216  if (isSignal_D0_InvM && isSignal_delta_m) {
217  // Signal region
218  // PID of K
219  m_h_D0_K_PID_ALL_kaon->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion, 1.));
220  if (pid_K->isAvailable(Const::SVD)) m_h_D0_K_PID_SVD_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
221  1., Const::SVD));
222  if (pid_K->isAvailable(Const::CDC)) m_h_D0_K_PID_CDC_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
223  1., Const::CDC));
224  if (pid_K->isAvailable(Const::TOP)) m_h_D0_K_PID_TOP_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
225  1., Const::TOP));
226  if (pid_K->isAvailable(Const::ARICH)) m_h_D0_K_PID_ARICH_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
227  1., Const::ARICH));
228  if (pid_K->isAvailable(Const::ECL)) m_h_D0_K_PID_ECL_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
229  1., Const::ECL));
230  if (pid_K->isAvailable(Const::KLM)) m_h_D0_K_PID_KLM_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion,
231  1., Const::KLM));
232  // PID of pi
233  m_h_D0_pi_PID_ALL_pion->Fill(pid_Pi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
234  if (pid_Pi->isAvailable(Const::SVD)) m_h_D0_pi_PID_SVD_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
235  Belle2::Const::kaon, 1., Const::SVD));
236  if (pid_Pi->isAvailable(Const::CDC)) m_h_D0_pi_PID_CDC_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
237  Belle2::Const::kaon, 1., Const::CDC));
238  if (pid_Pi->isAvailable(Const::TOP)) m_h_D0_pi_PID_TOP_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
239  Belle2::Const::kaon, 1., Const::TOP));
240  if (pid_Pi->isAvailable(Const::ARICH)) m_h_D0_pi_PID_ARICH_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
241  Belle2::Const::kaon, 1., Const::ARICH));
242  if (pid_Pi->isAvailable(Const::ECL)) m_h_D0_pi_PID_ECL_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
243  Belle2::Const::kaon, 1., Const::ECL));
244  if (pid_Pi->isAvailable(Const::KLM)) m_h_D0_pi_PID_KLM_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
245  Belle2::Const::kaon, 1., Const::KLM));
246  // PID of soft pi
247  m_h_D0_softpi_PID_ALL_pion->Fill(pid_softPi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
248  if (pid_softPi->isAvailable(Const::SVD)) m_h_D0_softpi_PID_SVD_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
249  Belle2::Const::kaon, 1., Const::SVD));
250  if (pid_softPi->isAvailable(Const::CDC)) m_h_D0_softpi_PID_CDC_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
251  Belle2::Const::kaon, 1., Const::CDC));
252  if (pid_softPi->isAvailable(Const::TOP)) m_h_D0_softpi_PID_TOP_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
253  Belle2::Const::kaon, 1., Const::TOP));
254  if (pid_softPi->isAvailable(Const::ARICH)) m_h_D0_softpi_PID_ARICH_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
255  Belle2::Const::kaon, 1., Const::ARICH));
256  if (pid_softPi->isAvailable(Const::ECL)) m_h_D0_softpi_PID_ECL_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
257  Belle2::Const::kaon, 1., Const::ECL));
258  if (pid_softPi->isAvailable(Const::KLM)) m_h_D0_softpi_PID_KLM_pion ->Fill(pid_softPi->getProbability(Belle2::Const::pion,
259  Belle2::Const::kaon, 1., Const::KLM));
260  } else {
261  // Sideband region for BG subtraction
262  // PID of K
263  m_h_sideband_D0_K_PID_ALL_kaon->Fill(pid_K->getProbability(Belle2::Const::kaon, Belle2::Const::pion, 1.));
264  if (pid_K->isAvailable(Const::SVD)) m_h_sideband_D0_K_PID_SVD_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
265  Belle2::Const::pion, 1., Const::SVD));
266  if (pid_K->isAvailable(Const::CDC)) m_h_sideband_D0_K_PID_CDC_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
267  Belle2::Const::pion, 1., Const::CDC));
268  if (pid_K->isAvailable(Const::TOP)) m_h_sideband_D0_K_PID_TOP_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
269  Belle2::Const::pion, 1., Const::TOP));
270  if (pid_K->isAvailable(Const::ARICH)) m_h_sideband_D0_K_PID_ARICH_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
271  Belle2::Const::pion, 1., Const::ARICH));
272  if (pid_K->isAvailable(Const::ECL)) m_h_sideband_D0_K_PID_ECL_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
273  Belle2::Const::pion, 1., Const::ECL));
274  if (pid_K->isAvailable(Const::KLM)) m_h_sideband_D0_K_PID_KLM_kaon ->Fill(pid_K->getProbability(Belle2::Const::kaon,
275  Belle2::Const::pion, 1., Const::KLM));
276  // PID of pi
277  m_h_sideband_D0_pi_PID_ALL_pion->Fill(pid_Pi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
278  if (pid_Pi->isAvailable(Const::SVD)) m_h_sideband_D0_pi_PID_SVD_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
279  Belle2::Const::kaon, 1., Const::SVD));
280  if (pid_Pi->isAvailable(Const::CDC)) m_h_sideband_D0_pi_PID_CDC_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
281  Belle2::Const::kaon, 1., Const::CDC));
282  if (pid_Pi->isAvailable(Const::TOP)) m_h_sideband_D0_pi_PID_TOP_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
283  Belle2::Const::kaon, 1., Const::TOP));
284  if (pid_Pi->isAvailable(Const::ARICH)) m_h_sideband_D0_pi_PID_ARICH_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
285  Belle2::Const::kaon, 1., Const::ARICH));
286  if (pid_Pi->isAvailable(Const::ECL)) m_h_sideband_D0_pi_PID_ECL_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
287  Belle2::Const::kaon, 1., Const::ECL));
288  if (pid_Pi->isAvailable(Const::KLM)) m_h_sideband_D0_pi_PID_KLM_pion ->Fill(pid_Pi->getProbability(Belle2::Const::pion,
289  Belle2::Const::kaon, 1., Const::KLM));
290  // PID of soft pi
291  m_h_sideband_D0_softpi_PID_ALL_pion->Fill(pid_softPi->getProbability(Belle2::Const::pion, Belle2::Const::kaon, 1.));
292  if (pid_softPi->isAvailable(Const::SVD)) m_h_sideband_D0_softpi_PID_SVD_pion ->Fill(pid_softPi->getProbability(
293  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::SVD));
294  if (pid_softPi->isAvailable(Const::CDC)) m_h_sideband_D0_softpi_PID_CDC_pion ->Fill(pid_softPi->getProbability(
295  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::CDC));
296  if (pid_softPi->isAvailable(Const::TOP)) m_h_sideband_D0_softpi_PID_TOP_pion ->Fill(pid_softPi->getProbability(
297  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::TOP));
298  if (pid_softPi->isAvailable(Const::ARICH)) m_h_sideband_D0_softpi_PID_ARICH_pion ->Fill(pid_softPi->getProbability(
299  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::ARICH));
300  if (pid_softPi->isAvailable(Const::ECL)) m_h_sideband_D0_softpi_PID_ECL_pion ->Fill(pid_softPi->getProbability(
301  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::ECL));
302  if (pid_softPi->isAvailable(Const::KLM)) m_h_sideband_D0_softpi_PID_KLM_pion ->Fill(pid_softPi->getProbability(
303  Belle2::Const::pion, Belle2::Const::kaon, 1., Const::KLM));
304  }
305  }
306  }
307 }
308 
309 // void PhysicsObjectsMiraBelleDstModule::getPIDInformationMiraBelle(Particle* part, float pid[3][7])
310 // {
311 // PIDLikelihood* I_like_coffee = part->getPIDLikelihood();
312 // pid[0][0] = I_like_coffee->getProbability();
313 // }
314 
315 void PhysicsObjectsMiraBelleDstModule::endRun()
316 {
317 }
318 
319 void PhysicsObjectsMiraBelleDstModule::terminate()
320 {
321 }
322 
Belle2::Const::pion
static const ChargedStable pion
charged pion particle
Definition: Const.h:652
Belle2::Const::kaon
static const ChargedStable kaon
charged kaon particle
Definition: Const.h:653
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::PIDLikelihood
Class to collect log likelihoods from TOP, ARICH, dEdx, ECL and KLM aimed for output to mdst includes...
Definition: PIDLikelihood.h:26
Belle2::PIDLikelihood::isAvailable
bool isAvailable(Const::PIDDetectorSet set) const
Check whether PID information from a given set of detectors is available.
Definition: PIDLikelihood.h:50
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:83
Belle2::Particle
Class to store reconstructed particles.
Definition: Particle.h:75
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:875
Belle2::Particle::getDaughter
const Particle * getDaughter(unsigned i) const
Returns a pointer to the i-th daughter particle.
Definition: Particle.cc:635
Belle2::Particle::getMass
double getMass() const
Returns invariant mass (= nominal for FS particles)
Definition: Particle.h:479
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:165
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:87
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:177
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:192
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:156
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:93
Belle2::PhysicsObjectsMiraBelleDstModule::m_dstListName
std::string m_dstListName
Name of the mu+ particle list.
Definition: PhysicsObjectsMiraBelleDstModule.h:204
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:162
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:198
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:150
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:84
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:183
Belle2::PhysicsObjectsMiraBelleDstModule::event
void event() override
This method is called for each event.
Definition: PhysicsObjectsMiraBelleDstModule.cc:171
Belle2::PhysicsObjectsMiraBelleDstModule::endRun
void endRun() override
This method is called if the current run ends.
Definition: PhysicsObjectsMiraBelleDstModule.cc:315
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:141
Belle2::PhysicsObjectsMiraBelleDstModule::terminate
void terminate() override
This method is called at the end of the event processing.
Definition: PhysicsObjectsMiraBelleDstModule.cc:319
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:189
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:168
Belle2::PhysicsObjectsMiraBelleDstModule::m_triggerIdentifier
std::string m_triggerIdentifier
Trigger identifier string used to select events for the histograms.
Definition: PhysicsObjectsMiraBelleDstModule.h:201
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:153
Belle2::PhysicsObjectsMiraBelleDstModule::m_h_delta_m
TH1F * m_h_delta_m
Mass difference b/w D*+ and D0.
Definition: PhysicsObjectsMiraBelleDstModule.h:72
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:90
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:138
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:159
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:75
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:78
Belle2::PhysicsObjectsMiraBelleDstModule::beginRun
void beginRun() override
Called when entering a new run.
Definition: PhysicsObjectsMiraBelleDstModule.cc:123
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:195
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:171
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:81
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:180
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:147
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:174
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:186
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:144
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:117
Belle2::PhysicsObjectsMiraBelleDstModule::defineHisto
void defineHisto() override
Definition of the histograms.
Definition: PhysicsObjectsMiraBelleDstModule.cc:45
Belle2::StoreObjPtr
Type-safe access to single objects in the data store.
Definition: StoreObjPtr.h:96
Belle2::REG_MODULE
REG_MODULE(arichBtest)
Register the Module.
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:560
Belle2::SoftwareTriggerCutResult::c_accept
@ c_accept
Accept this event.
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17