Belle II Software development
roe_variables.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#include <gtest/gtest.h>
9#include "utilities/TestParticleFactory.h"
10
11#include <analysis/VariableManager/Manager.h>
12#include <analysis/dataobjects/Particle.h>
13#include <analysis/dataobjects/ParticleList.h>
14#include <analysis/dataobjects/RestOfEvent.h>
15#include <analysis/utility/PCmsLabTransform.h>
16
17#include <mdst/dataobjects/Track.h>
18#include <mdst/dataobjects/ECLCluster.h>
19#include <mdst/dataobjects/KLMCluster.h>
20#include <mdst/dataobjects/PIDLikelihood.h>
21
22#include <framework/datastore/StoreArray.h>
23
24using namespace std;
25using namespace Belle2;
26using namespace Belle2::Variable;
27using namespace ROOT::Math;
28
29namespace {
30 class ROEVariablesTest : public ::testing::Test {
31 protected:
33 void SetUp() override
34 {
35
36 StoreObjPtr<ParticleList> pi0ParticleList("pi0:vartest");
37 StoreObjPtr<ParticleList> b0ParticleList("B0:vartest");
38 DataStore::Instance().setInitializeActive(true);
39 pi0ParticleList.registerInDataStore(DataStore::c_DontWriteOut);
40 b0ParticleList.registerInDataStore(DataStore::c_DontWriteOut);
41 StoreArray<ECLCluster> myECLClusters;
42 StoreArray<KLMCluster> myKLMClusters;
43 StoreArray<TrackFitResult> myTFRs;
44 StoreArray<Track> myTracks;
45 StoreArray<Particle> myParticles;
46 StoreArray<RestOfEvent> myROEs;
47 StoreObjPtr<RestOfEvent> roeobjptr;
48 StoreArray<PIDLikelihood> myPIDLikelihoods;
49 myECLClusters.registerInDataStore();
50 myKLMClusters.registerInDataStore();
51 myTFRs.registerInDataStore();
52 myTracks.registerInDataStore();
53 myParticles.registerInDataStore();
54 myROEs.registerInDataStore();
55 myPIDLikelihoods.registerInDataStore();
56 myParticles.registerRelationTo(myROEs);
57 myTracks.registerRelationTo(myPIDLikelihoods);
58 roeobjptr.registerInDataStore("RestOfEvent", DataStore::c_DontWriteOut);
59 DataStore::Instance().setInitializeActive(false);
60
61 pi0ParticleList.create();
62 pi0ParticleList->initialize(111, "pi0:vartest");
63 b0ParticleList.create();
64 b0ParticleList->initialize(521, "B0:vartest");
65
66 PCmsLabTransform T;
67
68 TestUtilities::TestParticleFactory factory;
69 ROOT::Math::XYZVector ipposition(0, 0, 0);
70 double halfEcms = T.getCMSEnergy() / 2;
71
72 PxPyPzEVector e_momentum(0., 0, halfEcms / 2, halfEcms / 2);
73 e_momentum = T.rotateCmsToLab() * e_momentum;
74 PxPyPzEVector p_momentum(0., 0, -halfEcms / 2, halfEcms / 2);
75 p_momentum = T.rotateCmsToLab() * p_momentum;
76
77 PxPyPzEVector b0_momentum(0, 0, 0, halfEcms);
78 b0_momentum = T.rotateCmsToLab() * b0_momentum;
79 factory.produceParticle(string("^B0 -> e- e+"), b0_momentum, ipposition);
80
81 myParticles[0]->set4Vector(e_momentum);
82 myParticles[1]->set4Vector(p_momentum);
83
84 myParticles[0]->print(); // e-
85 PxPyPzEVector fsp1_momentum(0., 0, halfEcms / 4, halfEcms / 4);
86 fsp1_momentum = T.rotateCmsToLab() * fsp1_momentum;
87 PxPyPzEVector fsp2_momentum(0., 0, -halfEcms / 4, halfEcms / 4);
88 fsp2_momentum = T.rotateCmsToLab() * fsp2_momentum;
89 PxPyPzEVector kl_momentum(0., 0, 0.1, 0.5);
90 kl_momentum = T.rotateCmsToLab() * fsp2_momentum;
91 factory.produceParticle(string("^B0 -> [pi0 -> gamma gamma] [K_S0 -> pi+ pi-]"), b0_momentum, ipposition);
92 KLMCluster myROEKLM;
93 KLMCluster* savedROEKLM = myKLMClusters.appendNew(myROEKLM);
94 Particle* roeKLMParticle = myParticles.appendNew(savedROEKLM);
95
96 // Set momentum of daughters
97 myParticles[3]->set4Vector(fsp1_momentum); // gamma
98 myParticles[4]->set4Vector(fsp2_momentum); // gamma
99 myECLClusters[0]->setEnergy(fsp1_momentum.E()); // gamma
100 myECLClusters[1]->setEnergy(fsp2_momentum.E()); // gamma
101 myParticles[6]->set4Vector(fsp1_momentum); // pi+
102 myParticles[7]->set4Vector(fsp2_momentum); // pi-
103 myParticles[3]->print(); // gamma
104 myParticles[4]->print(); // gamma
105 myParticles[7]->set4Vector(kl_momentum); // K_L0
106 myParticles[7]->print(); // K_L0
107 // Create ROE
108 RestOfEvent roe;
109 vector<const Particle*> roeParticlesToAdd;
110 // Add particles to the ROE
111 roeParticlesToAdd.push_back(myParticles[3]);
112 roeParticlesToAdd.push_back(myParticles[4]);
113 roeParticlesToAdd.push_back(myParticles[6]);
114 roeParticlesToAdd.push_back(myParticles[7]);
115 roeParticlesToAdd.push_back(roeKLMParticle);
116 roe.addParticles(roeParticlesToAdd);
117 RestOfEvent* savedROE = myROEs.appendNew(roe);
118 // Add relation ROE - particle (first B)
119 myParticles[2]->addRelationTo(savedROE); // B0
120 savedROE->print();
121 savedROE->initializeMask("my_mask", "test");
122 std::shared_ptr<Variable::Cut> chargedSelection = std::shared_ptr<Variable::Cut>
123 (Variable::Cut::compile("charge > 0")); // - exclude pi
124 std::shared_ptr<Variable::Cut> photonSelection = std::shared_ptr<Variable::Cut>
125 (Variable::Cut::compile("p > 1.5")); // - exclude gamma
126 std::shared_ptr<Variable::Cut> klSelection = std::shared_ptr<Variable::Cut>
127 (Variable::Cut::compile("E < 0")); // - exclude Klong
128 // Add mask, which should have 1 gamma, 1 pi+, 0 K_L0
129 savedROE->updateMaskWithCuts("my_mask", chargedSelection, photonSelection, klSelection);
130 savedROE->print("my_mask");
131 // Add pi0 from ROE to particle list
132 pi0ParticleList->addParticle(5, 111, Particle::c_Unflavored);
133 b0ParticleList->addParticle(2, 521, Particle::c_Unflavored);
134 }
135
137 void TearDown() override
138 {
139 DataStore::Instance().reset();
140 }
141 };
142
143 /*
144 * Test ROE particle composition variables
145 */
146 TEST_F(ROEVariablesTest, ROEParticleCompositionVariables)
147 {
148 StoreArray<Particle> myParticles{};
149 auto part = myParticles[2]; // B0
150 auto* var = Manager::Instance().getVariable("nROE_Charged()");
151 ASSERT_NE(var, nullptr);
152 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 2);
153
154 var = Manager::Instance().getVariable("nROE_Charged(my_mask)");
155 ASSERT_NE(var, nullptr);
156 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
157
158 var = Manager::Instance().getVariable("nROE_Charged(all)");
159 ASSERT_NE(var, nullptr);
160 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 2);
161
162 var = Manager::Instance().getVariable("nROE_Charged(my_mask, 13)");
163 ASSERT_NE(var, nullptr);
164 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 0);
165
166 var = Manager::Instance().getVariable("nROE_Charged(all, 13)");
167 ASSERT_NE(var, nullptr);
168 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 0);
169
170 var = Manager::Instance().getVariable("nROE_Charged(my_mask, 211)");
171 ASSERT_NE(var, nullptr);
172 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
173
174 var = Manager::Instance().getVariable("nROE_Photons()");
175 ASSERT_NE(var, nullptr);
176 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 2);
177
178 var = Manager::Instance().getVariable("nROE_Photons(my_mask)");
179 ASSERT_NE(var, nullptr);
180 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
181
182 var = Manager::Instance().getVariable("nROE_NeutralHadrons()");
183 ASSERT_NE(var, nullptr);
184 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
185
186 var = Manager::Instance().getVariable("nROE_NeutralHadrons(my_mask)");
187 ASSERT_NE(var, nullptr);
188 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 0);
189
190 }
191 /*
192 * Test ROE track/cluster composition variables
193 */
194 TEST_F(ROEVariablesTest, ROETrackClusterCompositionVariables)
195 {
196 StoreArray<Particle> myParticles{};
197 auto part = myParticles[2]; // B0
198
199 auto* var = Manager::Instance().getVariable("nROE_Tracks()");
200 ASSERT_NE(var, nullptr);
201 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 2);
202
203 var = Manager::Instance().getVariable("nROE_Tracks(my_mask)");
204 ASSERT_NE(var, nullptr);
205 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
206
207 var = Manager::Instance().getVariable("nROE_Tracks(all)");
208 ASSERT_NE(var, nullptr);
209 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 2);
210
211 var = Manager::Instance().getVariable("nROE_ECLClusters()");
212 ASSERT_NE(var, nullptr);
213 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 2);
214
215 var = Manager::Instance().getVariable("nROE_ECLClusters(my_mask)");
216 ASSERT_NE(var, nullptr);
217 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
218
219 var = Manager::Instance().getVariable("nROE_NeutralECLClusters()");
220 ASSERT_NE(var, nullptr);
221 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 2);
222
223 var = Manager::Instance().getVariable("nROE_NeutralECLClusters(my_mask)");
224 ASSERT_NE(var, nullptr);
225 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
226
227 var = Manager::Instance().getVariable("nROE_KLMClusters");
228 ASSERT_NE(var, nullptr);
229 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
230 }
231 /*
232 * Test ROE recoil frame variable
233 */
234 TEST_F(ROEVariablesTest, ROERecoilFrameVariable)
235 {
236 StoreArray<RestOfEvent> myROEs{};
237 StoreArray<Particle> myParticles{};
238 auto part = myParticles[2]; // B0
239 auto partNotROE = myParticles[0]; // electron has no ROE
240
241 auto* var = Manager::Instance().getVariable("useROERecoilFrame(E)");
242 ASSERT_NE(var, nullptr);
243 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 5.2986979);
244
245 DataStore::StoreEntry& roeobjptr = DataStore::Instance().getStoreEntryMap(DataStore::c_Event).at("RestOfEvent");
246 roeobjptr.object = myROEs[0];
247 roeobjptr.ptr = myROEs[0];
248
249 var = Manager::Instance().getVariable("useROERecoilFrame(E)");
250 ASSERT_NE(var, nullptr);
251 EXPECT_FLOAT_EQ(std::get<double>(var->function(partNotROE)), 2.801749);
252
253 // Clear ptr at the end
254 roeobjptr.object = nullptr;
255 roeobjptr.ptr = nullptr;
256 }
257 /*
258 * Test ROE kinematics variables
259 */
260 TEST_F(ROEVariablesTest, ROEKinematicsVariables)
261 {
262 StoreArray<Particle> myParticles;
263 auto part = myParticles[2]; // B0
264
265 // Tag side 4 vector
266 PxPyPzEVector roe4Vec(0, 0, 0, 0);
267 roe4Vec += myParticles[3]->get4Vector();
268 roe4Vec += myParticles[4]->get4Vector();
269 roe4Vec += myParticles[6]->get4Vector();
270 roe4Vec += myParticles[7]->get4Vector();
271 // Tag side 4 vector in mask
272 PxPyPzEVector mask4Vec(0, 0, 0, 0);
273 mask4Vec += myParticles[3]->get4Vector();
274 mask4Vec += myParticles[6]->get4Vector();
275 // Signal side 4 vector
276 PxPyPzEVector sig4Vec = part->get4Vector();
277
278 PCmsLabTransform T;
279 double E0 = T.getCMSEnergy() / 2;
280 auto roe4VecCMS = T.rotateLabToCms() * roe4Vec;
281 auto mask4VecCMS = T.rotateLabToCms() * mask4Vec;
282 auto sig4VecCMS = T.rotateLabToCms() * sig4Vec;
283
284 auto* var = Manager::Instance().getVariable("roeCharge()");
285 ASSERT_NE(var, nullptr);
286 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 0.0);
287
288 var = Manager::Instance().getVariable("roeCharge(all)");
289 ASSERT_NE(var, nullptr);
290 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 0.0);
291
292 var = Manager::Instance().getVariable("roeCharge(my_mask)");
293 ASSERT_NE(var, nullptr);
294 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 1.0);
295
296 var = Manager::Instance().getVariable("roeEextra()");
297 ASSERT_NE(var, nullptr);
298 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), myParticles[3]->getEnergy() + myParticles[4]->getEnergy());
299
300 var = Manager::Instance().getVariable("roeE()");
301 ASSERT_NE(var, nullptr);
302 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4Vec.E());
303
304 var = Manager::Instance().getVariable("useCMSFrame(roeE())");
305 ASSERT_NE(var, nullptr);
306 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4VecCMS.E());
307
308 var = Manager::Instance().getVariable("roeM()");
309 ASSERT_NE(var, nullptr);
310 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4Vec.mag());
311
312 var = Manager::Instance().getVariable("roeP()");
313 ASSERT_NE(var, nullptr);
314 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4Vec.P());
315
316 var = Manager::Instance().getVariable("useCMSFrame(roeP())");
317 ASSERT_NE(var, nullptr);
318 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4VecCMS.P());
319
320 var = Manager::Instance().getVariable("roePTheta()");
321 ASSERT_NE(var, nullptr);
322 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4Vec.Theta());
323
324 var = Manager::Instance().getVariable("useCMSFrame(roePTheta())");
325 ASSERT_NE(var, nullptr);
326 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4VecCMS.Theta());
327
328 var = Manager::Instance().getVariable("roeDeltae()");
329 ASSERT_NE(var, nullptr);
330 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), roe4VecCMS.E() - E0);
331
332 var = Manager::Instance().getVariable("roeDeltae(my_mask)");
333 ASSERT_NE(var, nullptr);
334 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), mask4VecCMS.E() - E0);
335
336 var = Manager::Instance().getVariable("roeMbc()");
337 ASSERT_NE(var, nullptr);
338 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), TMath::Sqrt(E0 * E0 - roe4VecCMS.P2()));
339
340 var = Manager::Instance().getVariable("roeMbc(my_mask)");
341 ASSERT_NE(var, nullptr);
342 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), TMath::Sqrt(E0 * E0 - mask4VecCMS.P2()));
343
344 var = Manager::Instance().getVariable("useCMSFrame(weDeltae(my_mask))");
345 ASSERT_NE(var, nullptr);
346 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), mask4VecCMS.E() + sig4VecCMS.E() - E0);
347
348 var = Manager::Instance().getVariable("useCMSFrame(weMbc(my_mask))");
349 ASSERT_NE(var, nullptr);
350 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), TMath::Sqrt(E0 * E0 - mask4VecCMS.P2()));
351
352 var = Manager::Instance().getVariable("useCMSFrame(weMbc(all))");
353 ASSERT_NE(var, nullptr);
354 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), TMath::Sqrt(E0 * E0 - roe4VecCMS.P2()));
355
356 PxPyPzEVector miss4VecCMS = - (sig4VecCMS + mask4VecCMS);
357 miss4VecCMS.SetE(2 * E0 - (sig4VecCMS.E() + mask4VecCMS.E()));
358
359 var = Manager::Instance().getVariable("useCMSFrame(weMissM2(my_mask,0))");
360 ASSERT_NE(var, nullptr);
361 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 5.3290705e-15);
362
363 var = Manager::Instance().getVariable("useCMSFrame(weMissP(my_mask,0))");
364 ASSERT_NE(var, nullptr);
365 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), miss4VecCMS.P());
366
367 var = Manager::Instance().getVariable("useCMSFrame(weMissPx(my_mask,0))");
368 ASSERT_NE(var, nullptr);
369 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 0.);
370
371 var = Manager::Instance().getVariable("useCMSFrame(weMissPy(my_mask,0))");
372 ASSERT_NE(var, nullptr);
373 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 0.);
374
375 var = Manager::Instance().getVariable("useCMSFrame(weMissPz(my_mask,0))");
376 ASSERT_NE(var, nullptr);
377 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), - E0 / 2);
378
379 var = Manager::Instance().getVariable("useCMSFrame(weMissE(my_mask,0))");
380 ASSERT_NE(var, nullptr);
381 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), miss4VecCMS.E());
382
383 }
384 /*
385 * Test specific kinematic variables
386 */
387 TEST_F(ROEVariablesTest, ROESpecificKinematicVariables)
388 {
389 StoreArray<Particle> myParticles{};
390 auto part = myParticles[2]; // B0
391
392 auto* var = Manager::Instance().getVariable("weCosThetaEll()");
393 ASSERT_NE(var, nullptr);
394 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), -0.99899036);
395
396 var = Manager::Instance().getVariable("weXiZ()");
397 ASSERT_NE(var, nullptr);
398 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 0.31122509);
399
400 var = Manager::Instance().getVariable("useCMSFrame(weQ2lnuSimple(my_mask))");
401 ASSERT_NE(var, nullptr);
402 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 0.);
403
404 // FIXME: More complicated test setup is required to pass abs(cos_angle_nu) < 1
405 var = Manager::Instance().getVariable("useCMSFrame(weQ2lnu(my_mask))");
406 ASSERT_NE(var, nullptr);
407 EXPECT_FLOAT_EQ(std::get<double>(var->function(part)), 1.9539925e-14);
408 }
409 /*
410 * Test isInROE variables
411 */
412 TEST_F(ROEVariablesTest, IsInROEVariables)
413 {
414 StoreArray<Particle> myParticles{};
415 auto part = myParticles[2]; // B0
416 StoreArray<RestOfEvent> myROEs{};
417
418 DataStore::StoreEntry& roeobjptr = DataStore::Instance().getStoreEntryMap(DataStore::c_Event).at("RestOfEvent");
419 roeobjptr.object = myROEs[0];
420 roeobjptr.ptr = myROEs[0];
421
422 auto partROE1 = myParticles[3]; // gamma from ROE
423 auto partROE2 = myParticles[5]; // pi from ROE
424 auto partROE3 = myParticles[10]; // K_L0 from ROE
425 auto partROEnotFromMask = myParticles[7]; // pi from ROE not from mask
426 partROEnotFromMask->print(); // gamma
427
428 auto partNotROE = myParticles[0]; // e NOT from ROE
429
430 auto* var = Manager::Instance().getVariable("isInRestOfEvent");
431 ASSERT_NE(var, nullptr);
432 EXPECT_EQ(std::get<bool>(var->function(partROE1)), 1);
433
434 var = Manager::Instance().getVariable("isInRestOfEvent");
435 ASSERT_NE(var, nullptr);
436 EXPECT_EQ(std::get<bool>(var->function(partROE2)), 1);
437
438 var = Manager::Instance().getVariable("isInRestOfEvent");
439 ASSERT_NE(var, nullptr);
440 EXPECT_EQ(std::get<bool>(var->function(partROE3)), 1);
441
442 var = Manager::Instance().getVariable("isInRestOfEvent");
443 ASSERT_NE(var, nullptr);
444 EXPECT_EQ(std::get<bool>(var->function(partNotROE)), 0);
445
446 var = Manager::Instance().getVariable("isInRestOfEvent");
447 ASSERT_NE(var, nullptr);
448 EXPECT_EQ(std::get<bool>(var->function(partNotROE)), 0);
449
450 var = Manager::Instance().getVariable("passesROEMask(my_mask)");
451 ASSERT_NE(var, nullptr);
452 EXPECT_FLOAT_EQ(std::get<bool>(var->function(partROE1)), 1);
453
454 var = Manager::Instance().getVariable("passesROEMask(all)");
455 ASSERT_NE(var, nullptr);
456 EXPECT_FLOAT_EQ(std::get<bool>(var->function(partROE1)), 1);
457
458 var = Manager::Instance().getVariable("passesROEMask(all)");
459 ASSERT_NE(var, nullptr);
460 EXPECT_FLOAT_EQ(std::get<bool>(var->function(partNotROE)), 0);
461
462 var = Manager::Instance().getVariable("passesROEMask(my_mask)");
463 ASSERT_NE(var, nullptr);
464 EXPECT_FLOAT_EQ(std::get<bool>(var->function(partROEnotFromMask)), 0);
465
466 var = Manager::Instance().getVariable("nROE_ParticlesInList(pi0:vartest)");
467 ASSERT_NE(var, nullptr);
468 EXPECT_FLOAT_EQ(std::get<int>(var->function(part)), 1);
469
470 var = Manager::Instance().getVariable("currentROEIsInList(B0:vartest)");
471 ASSERT_NE(var, nullptr);
472 EXPECT_FLOAT_EQ(std::get<bool>(var->function(partROE1)), 1);
473
474 var = Manager::Instance().getVariable("currentROEIsInList(pi0:vartest)");
475 ASSERT_NE(var, nullptr);
476 EXPECT_FLOAT_EQ(std::get<bool>(var->function(partROE1)), 0);
477
478 var = Manager::Instance().getVariable("particleRelatedToCurrentROE(PDG)");
479 ASSERT_NE(var, nullptr);
480 EXPECT_FLOAT_EQ(std::get<double>(var->function(partROE1)), 511.0);
481
482 // Clear ptr at the end
483 roeobjptr.object = nullptr;
484 roeobjptr.ptr = nullptr;
485 }
486}
Belle2::DataStore::c_DontWriteOut
@ c_DontWriteOut
Object/array should be NOT saved by output modules.
Definition DataStore.h:71
Belle2::DataStore::getStoreEntryMap
StoreEntryMap & getStoreEntryMap(EDurability durability)
Get a reference to the object/array map.
Definition DataStore.h:325
Belle2::DataStore::c_Event
@ c_Event
Different object in each event, all objects/arrays are invalidated after event() function has been ca...
Definition DataStore.h:59
Belle2::DataStore::StoreEntry
Belle2::StoreEntry StoreEntry
Wraps a stored array/object, stored under unique (name, durability) key.
Definition DataStore.h:84
Belle2::DataStore::Instance
static DataStore & Instance()
Instance of singleton Store.
Definition DataStore.cc:53
Belle2::DataStore::setInitializeActive
void setInitializeActive(bool active)
Setter for m_initializeActive.
Definition DataStore.cc:93
Belle2::DataStore::reset
void reset(EDurability durability)
Frees memory occupied by data store items and removes all objects from the map.
Definition DataStore.cc:85
Belle2::GeneralCut< Manager >::compile
static std::unique_ptr< GeneralCut > compile(const std::string &cut)
Definition GeneralCut.h:84
Belle2::PCmsLabTransform
Class to hold Lorentz transformations from/to CMS and boost vector.
Definition PCmsLabTransform.h:29
Belle2::PCmsLabTransform::getCMSEnergy
double getCMSEnergy() const
Returns CMS energy of e+e- (aka.
Definition PCmsLabTransform.h:49
Belle2::PCmsLabTransform::rotateCmsToLab
const ROOT::Math::LorentzRotation rotateCmsToLab() const
Returns Lorentz transformation from CMS to Lab.
Definition PCmsLabTransform.h:75
Belle2::Particle::c_Unflavored
@ c_Unflavored
Is its own antiparticle or we don't know whether it is a particle/antiparticle.
Definition Particle.h:97
Belle2::RestOfEvent::print
void print(const std::string &maskName=c_defaultMaskName, bool unpackComposite=true) const
Prints the contents of a RestOfEvent object to screen.
Definition RestOfEvent.cc:386
Belle2::RestOfEvent::initializeMask
void initializeMask(const std::string &name, const std::string &origin="unknown")
Initialize new mask.
Definition RestOfEvent.cc:129
Belle2::RestOfEvent::updateMaskWithCuts
void updateMaskWithCuts(const std::string &name, const std::shared_ptr< Variable::Cut > &trackCut=nullptr, const std::shared_ptr< Variable::Cut > &eclCut=nullptr, const std::shared_ptr< Variable::Cut > &klmCut=nullptr, bool updateExisting=false)
Update mask with cuts.
Definition RestOfEvent.cc:186
Belle2::RestOfEvent::addParticles
void addParticles(const std::vector< const Particle * > &particle)
Add StoreArray indices of given Particles to the list of unused particles in the event.
Definition RestOfEvent.cc:19
Belle2::StoreAccessorBase::registerInDataStore
bool registerInDataStore(DataStore::EStoreFlags storeFlags=DataStore::c_WriteOut)
Register the object/array in the DataStore.
Definition StoreAccessorBase.h:52
Belle2::StoreArray
Accessor to arrays stored in the data store.
Definition StoreArray.h:113
Belle2::StoreArray::appendNew
T * appendNew()
Construct a new T object at the end of the array.
Definition StoreArray.h:246
Belle2::StoreArray::registerRelationTo
bool registerRelationTo(const StoreArray< TO > &toArray, DataStore::EDurability durability=DataStore::c_Event, DataStore::EStoreFlags storeFlags=DataStore::c_WriteOut, const std::string &namedRelation="") const
Register a relation to the given StoreArray.
Definition StoreArray.h:140
Belle2::Variable::Manager::getVariable
const Var * getVariable(std::string name)
Get the variable belonging to the given key.
Definition Manager.cc:58
Belle2::Variable::Manager::Instance
static Manager & Instance()
get singleton instance.
Definition Manager.cc:26
TestUtilities::TestParticleFactory::produceParticle
const Belle2::Particle * produceParticle(const std::string &decayString, const ROOT::Math::PxPyPzEVector &momentum, const ROOT::Math::XYZVector &vertex)
Main method to produce particles.
Definition TestParticleFactory.h:41
Belle2
Abstract base class for different kinds of events.
Definition MillepedeAlgorithm.h:17
std
STL namespace.
Belle2::StoreEntry::ptr
TObject * ptr
The pointer to the returned object, either equal to 'object' or null, depending on whether the object...
Definition StoreEntry.h:51
Belle2::StoreEntry::object
TObject * object
The pointer to the actual object.
Definition StoreEntry.h:48