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

Merge branch 'master' into 'master'

Change in the number of events

See merge request !8
parents d171b9cf 3ff541ea
No related branches found
No related tags found
1 merge request!8Change in the number of events
Source diff could not be displayed: it is too large. Options to address this: view the blob.
...@@ -15,17 +15,17 @@ ...@@ -15,17 +15,17 @@
using namespace HepMC3; using namespace HepMC3;
void zdc_photons(){ void zdc_photons(){
WriterAscii hepmc_output("./data/zdc_photons.hepmc"); WriterAscii hepmc_output("./data/zdc_photons_evt500.hepmc");
int events_parsed = 0; int events_parsed = 0;
GenEvent evt(Units::GEV, Units::MM); GenEvent evt(Units::GEV, Units::MM);
// Random number pulled from the env variable // Random number pulled from the env variable
std::mt19937 gen( std::stoi(std::getenv("SEED")) ); std::mt19937 gen( std::stoi(std::getenv("SEED")) );
std::uniform_real_distribution<double> uniform_theta(0.0,2.0*TMath::DegToRad()); // 2[degree] std::uniform_real_distribution<double> uniform_theta(0.0,0.1*TMath::DegToRad()); // 0.1[degree]
std::uniform_real_distribution<double> uniform_phi(0.0,2.0*TMath::Pi()); // 360[degree] std::uniform_real_distribution<double> uniform_phi(0.0,2.0*TMath::Pi()); // 360[degree]
for (events_parsed = 0; events_parsed < 100; events_parsed++) { for (events_parsed = 0; events_parsed < 500; 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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment