Skip to content
Snippets Groups Projects

feat: transition from eicd to edm4eic

Merged Wouter Deconinck requested to merge eicd-to-edm4eic into master
13 files
+ 140
140
Compare changes
  • Side-by-side
  • Inline
Files
13
@@ -10,8 +10,8 @@
@@ -10,8 +10,8 @@
#include "edm4hep/MCParticleCollection.h"
#include "edm4hep/MCParticleCollection.h"
#include "edm4hep/SimCalorimeterHitCollection.h"
#include "edm4hep/SimCalorimeterHitCollection.h"
#include "eicd/ClusterCollection.h"
#include "edm4eic/ClusterCollection.h"
#include "eicd/ClusterData.h"
#include "edm4eic/ClusterData.h"
#include "TCanvas.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TStyle.h"
@@ -23,7 +23,7 @@
@@ -23,7 +23,7 @@
#include "TH1D.h"
#include "TH1D.h"
#include "TVector3.h"
#include "TVector3.h"
R__LOAD_LIBRARY(libeicd.so)
R__LOAD_LIBRARY(libedm4eic.so)
#include "fmt/core.h"
#include "fmt/core.h"
#include "DD4hep/Detector.h"
#include "DD4hep/Detector.h"
#include "DDG4/Geant4Data.h"
#include "DDG4/Geant4Data.h"
@@ -82,7 +82,7 @@ void analysis_zdc_particles(
@@ -82,7 +82,7 @@ void analysis_zdc_particles(
auto theta = std::atan2(std::hypot(p.momentum.x, p.momentum.y), p.momentum.z);
auto theta = std::atan2(std::hypot(p.momentum.x, p.momentum.y), p.momentum.z);
return theta;
return theta;
};
};
auto Thetarec = [](const std::vector<eicd::ClusterData> & input) {
auto Thetarec = [](const std::vector<edm4eic::ClusterData> & input) {
std::vector<float> output;
std::vector<float> output;
for (const auto &c: input) {
for (const auto &c: input) {
auto r = c.position;
auto r = c.position;
@@ -99,7 +99,7 @@ void analysis_zdc_particles(
@@ -99,7 +99,7 @@ void analysis_zdc_particles(
if (phi < 0) phi += 2.*M_PI;
if (phi < 0) phi += 2.*M_PI;
return phi;
return phi;
};
};
auto Phirec = [](const std::vector<eicd::ClusterData> & input) {
auto Phirec = [](const std::vector<edm4eic::ClusterData> & input) {
std::vector<float> output;
std::vector<float> output;
for (const auto &c: input) {
for (const auto &c: input) {
auto r = c.position;
auto r = c.position;
@@ -117,7 +117,7 @@ void analysis_zdc_particles(
@@ -117,7 +117,7 @@ void analysis_zdc_particles(
auto eta = -std::log(std::tan(theta/2));
auto eta = -std::log(std::tan(theta/2));
return eta;
return eta;
};
};
auto Etarec = [](const std::vector<eicd::ClusterData> & input) {
auto Etarec = [](const std::vector<edm4eic::ClusterData> & input) {
std::vector<float> output;
std::vector<float> output;
for (const auto &c: input) {
for (const auto &c: input) {
auto r = c.position;
auto r = c.position;
Loading