Created new electron dataset for crystal calorimeter

emcal_electrons.cxx

@@ -11,11 +11,13 @@
 #include <iostream>
 #include<random>
 #include<cmath>
+#include <math.h>
+#include <TMath.h>
 
 using namespace HepMC3;
 
 void emcal_electrons(){
-        WriterAscii hepmc_output("./data/emcal_electrons.hepmc");
+        WriterAscii hepmc_output("./data/emcal_electrons_upto30GeV_10kevents.hepmc");
         int        events_parsed = 0;
         GenEvent   evt(Units::GEV, Units::MM);
 
@@ -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_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)
         // type 4 is beam
         // pdgid 11 - electron
@@ -34,7 +36,7 @@ void emcal_electrons(){
         GenParticlePtr p2 = std::make_shared<GenParticle>(FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
 
 	// 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 phi   = uniform_phi(gen); 
 	double px    = p*sin(theta)*cos(phi);

emcal_electrons_reader.cxx

@@ -26,7 +26,7 @@ void emcal_electrons_reader(){
   	gStyle->SetPadGridY(1);
   	gStyle->SetPadLeftMargin(0.14);
 
-	ReaderAscii hepmc_input("./data/emcal_electrons.hepmc");
+	ReaderAscii hepmc_input("./data/emcal_electrons_upto30GeV_10kevents.hepmc");
 	int        events_parsed = 0;
 	GenEvent   evt(Units::GEV, Units::MM);