|
| 1 | +#!/usr/bin/env python3 |
| 2 | +"""Test to exercise DataSource thoroughly""" |
| 3 | + |
| 4 | +import sys |
| 5 | +import ROOT |
| 6 | +from podio.data_source import CreateDataFrame # pylint: disable=import-error, no-name-in-module |
| 7 | + |
| 8 | +input_file = sys.argv[1] if len(sys.argv) > 1 else "example_frame.root" |
| 9 | +snapshot_file = sys.argv[2] if len(sys.argv) > 2 else "datasource_snapshot.root" |
| 10 | + |
| 11 | +if ROOT.gSystem.Load("libTestDataModelDict") < 0: |
| 12 | + raise RuntimeError("Could not load TestDataModel dictionary") |
| 13 | + |
| 14 | +ROOT.gInterpreter.ProcessLine("using namespace ROOT::VecOps;") |
| 15 | + |
| 16 | +ROOT.gInterpreter.Declare( |
| 17 | + """ |
| 18 | +#include <datamodel/ExampleHitCollection.h> |
| 19 | +#include <datamodel/ExampleClusterCollection.h> |
| 20 | +#include <datamodel/EventInfoCollection.h> |
| 21 | +#include <datamodel/TestLinkCollection.h> |
| 22 | +""" |
| 23 | +) |
| 24 | + |
| 25 | +# Declare helpers that extract quantities from the test collections |
| 26 | +ROOT.gInterpreter.Declare( |
| 27 | + """ |
| 28 | +RVec<double> getHitEnergies(const ExampleHitCollection& hits) { |
| 29 | + RVec<double> v; |
| 30 | + v.reserve(hits.size()); |
| 31 | + for (const auto& h : hits) { |
| 32 | + v.push_back(h.energy()); |
| 33 | + } |
| 34 | + return v; |
| 35 | +} |
| 36 | +
|
| 37 | +RVec<double> getClusterEnergies(const ExampleClusterCollection& clusters) { |
| 38 | + RVec<double> v; |
| 39 | + v.reserve(clusters.size()); |
| 40 | + for (const auto& c : clusters) { |
| 41 | + v.push_back(c.energy()); |
| 42 | + } |
| 43 | + return v; |
| 44 | +} |
| 45 | +
|
| 46 | +int getEventNumber(const EventInfoCollection& info) { |
| 47 | + return info[0].Number(); |
| 48 | +} |
| 49 | +
|
| 50 | +RVec<double> getCompositeClusterHitEnergies(const ExampleClusterCollection& clusters) { |
| 51 | + RVec<double> v; |
| 52 | + for (const auto& hit : clusters[2].Hits()) { |
| 53 | + v.push_back(hit.energy()); |
| 54 | + } |
| 55 | + return v; |
| 56 | +} |
| 57 | +
|
| 58 | +RVec<double> getSubClusterEnergies(const ExampleClusterCollection& clusters) { |
| 59 | + RVec<double> v; |
| 60 | + for (const auto& sub : clusters[2].Clusters()) { |
| 61 | + v.push_back(sub.energy()); |
| 62 | + } |
| 63 | + return v; |
| 64 | +} |
| 65 | +
|
| 66 | +// Follow each link From (hit) and To (cluster) relations and return: |
| 67 | +// [weight0, hit_energy0, cluster_energy0, weight1, hit_energy1, cluster_energy1] |
| 68 | +RVec<double> getLinkInfo(const TestLinkCollection& links) { |
| 69 | + RVec<double> v; |
| 70 | + for (const auto& link : links) { |
| 71 | + v.push_back(link.getWeight()); |
| 72 | + v.push_back(link.getFrom().energy()); |
| 73 | + v.push_back(link.getTo().energy()); |
| 74 | + } |
| 75 | + return v; |
| 76 | +} |
| 77 | +""" |
| 78 | +) |
| 79 | + |
| 80 | +rdf = CreateDataFrame(input_file) |
| 81 | + |
| 82 | +assert rdf.Count().GetValue() == 10, f"Expected 10 events in {input_file}" |
| 83 | + |
| 84 | +rdf = ( |
| 85 | + rdf.Define("hit_energies", "getHitEnergies(hits)") |
| 86 | + .Define("hit_energy_0", "hit_energies[0]") |
| 87 | + .Define("hit_energy_1", "hit_energies[1]") |
| 88 | + .Define("hit_energy_sum", "Sum(hit_energies)") |
| 89 | + .Define("cluster_energies", "getClusterEnergies(clusters)") |
| 90 | + .Define("cluster_energy_0", "cluster_energies[0]") |
| 91 | + .Define("cluster_energy_1", "cluster_energies[1]") |
| 92 | + .Define("cluster_energy_2", "cluster_energies[2]") |
| 93 | + .Define("event_number", "getEventNumber(info)") |
| 94 | + .Define("composite_hit_energies", "getCompositeClusterHitEnergies(clusters)") |
| 95 | + .Define("composite_hit_energy_0", "composite_hit_energies[0]") |
| 96 | + .Define("composite_hit_energy_1", "composite_hit_energies[1]") |
| 97 | + .Define("sub_cluster_energies", "getSubClusterEnergies(clusters)") |
| 98 | + .Define("sub_cluster_energy_0", "sub_cluster_energies[0]") |
| 99 | + .Define("sub_cluster_energy_1", "sub_cluster_energies[1]") |
| 100 | + .Define("link_info", "getLinkInfo(links)") |
| 101 | + .Define("link0_weight", "link_info[0]") |
| 102 | + .Define("link0_from_energy", "link_info[1]") |
| 103 | + .Define("link0_to_energy", "link_info[2]") |
| 104 | + .Define("link1_weight", "link_info[3]") |
| 105 | + .Define("link1_from_energy", "link_info[4]") |
| 106 | + .Define("link1_to_energy", "link_info[5]") |
| 107 | +) |
| 108 | + |
| 109 | +# print(rdf.Describe()) |
| 110 | +# print(rdf.Display(["event_number", "hit_energy_0", "cluster_energy_2", "composite_hit_energy_0", "sub_cluster_energy_0"]).AsString()) |
| 111 | + |
| 112 | +# 10 events, i = 0..9 |
| 113 | +# hit_energy_0 = 23+i |
| 114 | +mean_hit0 = rdf.Mean("hit_energy_0").GetValue() |
| 115 | +assert abs(mean_hit0 - 27.5) < 1e-9, f"Mean of hit_energy_0 should be 27.5, got {mean_hit0}" |
| 116 | + |
| 117 | +# cluster_energy_2 = 35+2*i |
| 118 | +mean_clu2 = rdf.Mean("cluster_energy_2").GetValue() |
| 119 | +assert abs(mean_clu2 - 44.0) < 1e-9, f"Mean of cluster_energy_2 should be 44.0, got {mean_clu2}" |
| 120 | + |
| 121 | +# event_number = i |
| 122 | +mean_evtnum = rdf.Mean("event_number").GetValue() |
| 123 | +assert abs(mean_evtnum - 4.5) < 1e-9, f"Mean of event_number should be 4.5, got {mean_evtnum}" |
| 124 | + |
| 125 | +# composite_hit_energy_0 follows clusters[2].Hits()[0] == hits[0], energy = 23+i -> Mean = 27.5 |
| 126 | +mean_chit0 = rdf.Mean("composite_hit_energy_0").GetValue() |
| 127 | +assert ( |
| 128 | + abs(mean_chit0 - 27.5) < 1e-9 |
| 129 | +), f"Mean of composite_hit_energy_0 (cluster->hit relation) should be 27.5, got {mean_chit0}" |
| 130 | + |
| 131 | +# sub_cluster_energy_0 follows clusters[2].Clusters()[0] == clusters[0], energy = 23+i -> Mean = 27.5 |
| 132 | +mean_sub0 = rdf.Mean("sub_cluster_energy_0").GetValue() |
| 133 | +assert ( |
| 134 | + abs(mean_sub0 - 27.5) < 1e-9 |
| 135 | +), f"Mean of sub_cluster_energy_0 (cluster->cluster relation) should be 27.5, got {mean_sub0}" |
| 136 | + |
| 137 | +# link0_weight is always 0.0, link1_weight is always 0.5 |
| 138 | +mean_w0 = rdf.Mean("link0_weight").GetValue() |
| 139 | +assert abs(mean_w0 - 0.0) < 1e-9, f"Mean of link0_weight should be 0.0, got {mean_w0}" |
| 140 | +mean_w1 = rdf.Mean("link1_weight").GetValue() |
| 141 | +assert abs(mean_w1 - 0.5) < 1e-9, f"Mean of link1_weight should be 0.5, got {mean_w1}" |
| 142 | + |
| 143 | +# link0 goes hits[0]->clusters[1]: from-energy = 23+i -> Mean = 27.5, |
| 144 | +# to-energy = 12+i -> Mean = 16.5 |
| 145 | +mean_l0from = rdf.Mean("link0_from_energy").GetValue() |
| 146 | +assert ( |
| 147 | + abs(mean_l0from - 27.5) < 1e-9 |
| 148 | +), f"Mean of link0_from_energy (link->hit relation) should be 27.5, got {mean_l0from}" |
| 149 | +mean_l0to = rdf.Mean("link0_to_energy").GetValue() |
| 150 | +assert ( |
| 151 | + abs(mean_l0to - 16.5) < 1e-9 |
| 152 | +), f"Mean of link0_to_energy (link->cluster relation) should be 16.5, got {mean_l0to}" |
| 153 | + |
| 154 | +rdf_even = rdf.Filter("event_number % 2 == 0") |
| 155 | +assert rdf_even.Count().GetValue() == 5, "Expected 5 even-numbered events after filter" |
| 156 | + |
| 157 | +columns = [ |
| 158 | + "event_number", |
| 159 | + "hit_energy_0", |
| 160 | + "hit_energy_1", |
| 161 | + "hit_energy_sum", |
| 162 | + "cluster_energy_0", |
| 163 | + "cluster_energy_1", |
| 164 | + "cluster_energy_2", |
| 165 | + "composite_hit_energy_0", |
| 166 | + "composite_hit_energy_1", |
| 167 | + "sub_cluster_energy_0", |
| 168 | + "sub_cluster_energy_1", |
| 169 | + "link0_weight", |
| 170 | + "link0_from_energy", |
| 171 | + "link0_to_energy", |
| 172 | + "link1_weight", |
| 173 | + "link1_from_energy", |
| 174 | + "link1_to_energy", |
| 175 | +] |
| 176 | + |
| 177 | +rdf.Snapshot("events", snapshot_file, columns) |
| 178 | + |
| 179 | +print(f"All assertions passed, snapshot written to {snapshot_file}") |
0 commit comments