Skip to content
Snippets Groups Projects
Commit 98ffc9c8 authored by Jihee Kim's avatar Jihee Kim
Browse files

Merge branch 'master' into 'master'

Electron data set fro crystal calorimeter

See merge request EIC/datasets!12
parents baf9c2ce 4b901085
No related branches found
No related tags found
No related merge requests found
This diff is collapsed.
Source diff could not be displayed: it is too large. Options to address this: view the blob.
...@@ -11,11 +11,13 @@ ...@@ -11,11 +11,13 @@
#include <iostream> #include <iostream>
#include<random> #include<random>
#include<cmath> #include<cmath>
#include <math.h>
#include <TMath.h>
using namespace HepMC3; using namespace HepMC3;
void emcal_electrons(){ void emcal_electrons(){
WriterAscii hepmc_output("./data/emcal_electrons.hepmc"); WriterAscii hepmc_output("./data/emcal_electrons_upto30GeV_10kevents.hepmc");
int events_parsed = 0; int events_parsed = 0;
GenEvent evt(Units::GEV, Units::MM); GenEvent evt(Units::GEV, Units::MM);
...@@ -24,7 +26,7 @@ void emcal_electrons(){ ...@@ -24,7 +26,7 @@ void emcal_electrons(){
std::uniform_real_distribution<double> uniform_theta(135.0*TMath::DegToRad(),178.0*TMath::DegToRad()); // 135-178[degree] std::uniform_real_distribution<double> uniform_theta(135.0*TMath::DegToRad(),178.0*TMath::DegToRad()); // 135-178[degree]
std::uniform_real_distribution<double> uniform_phi(0.0,2*TMath::Pi()); // 360[degree] std::uniform_real_distribution<double> uniform_phi(0.0,2*TMath::Pi()); // 360[degree]
for (events_parsed = 0; events_parsed < 100; events_parsed++) { for (events_parsed = 0; events_parsed < 10000; events_parsed++) {
// FourVector(px,py,pz,e,pdgid,status) // FourVector(px,py,pz,e,pdgid,status)
// type 4 is beam // type 4 is beam
// pdgid 11 - electron // pdgid 11 - electron
...@@ -34,7 +36,7 @@ void emcal_electrons(){ ...@@ -34,7 +36,7 @@ void emcal_electrons(){
GenParticlePtr p2 = std::make_shared<GenParticle>(FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4); GenParticlePtr p2 = std::make_shared<GenParticle>(FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
// Define variables - energy, theta, phi, momentum vectors // Define variables - energy, theta, phi, momentum vectors
double p = 1.0 + events_parsed*1.0; // temp. energy range 1 GeV to 100 GeV double p = 1.0 + events_parsed*0.003; // temp. energy range 1 GeV to 20 GeV
double theta = uniform_theta(gen); double theta = uniform_theta(gen);
double phi = uniform_phi(gen); double phi = uniform_phi(gen);
double px = p*sin(theta)*cos(phi); double px = p*sin(theta)*cos(phi);
......
...@@ -26,7 +26,7 @@ void emcal_electrons_reader(){ ...@@ -26,7 +26,7 @@ void emcal_electrons_reader(){
gStyle->SetPadGridY(1); gStyle->SetPadGridY(1);
gStyle->SetPadLeftMargin(0.14); gStyle->SetPadLeftMargin(0.14);
ReaderAscii hepmc_input("./data/emcal_electrons.hepmc"); ReaderAscii hepmc_input("./data/emcal_electrons_upto30GeV_10kevents.hepmc");
int events_parsed = 0; int events_parsed = 0;
GenEvent evt(Units::GEV, Units::MM); GenEvent evt(Units::GEV, Units::MM);
......
No preview for this file type
results/emcal_electrons_energy_reader.png

12.1 KiB | W: | H:

results/emcal_electrons_energy_reader.png

12 KiB | W: | H:

results/emcal_electrons_energy_reader.png
results/emcal_electrons_energy_reader.png
results/emcal_electrons_energy_reader.png
results/emcal_electrons_energy_reader.png
  • 2-up
  • Swipe
  • Onion skin
No preview for this file type
results/emcal_electrons_eta_reader.png

12.1 KiB | W: | H:

results/emcal_electrons_eta_reader.png

12.2 KiB | W: | H:

results/emcal_electrons_eta_reader.png
results/emcal_electrons_eta_reader.png
results/emcal_electrons_eta_reader.png
results/emcal_electrons_eta_reader.png
  • 2-up
  • Swipe
  • Onion skin
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment