Skip to content
Snippets Groups Projects
Commit f29ca606 authored by Alex Jentsch's avatar Alex Jentsch
Browse files

Change the generator file to generate protons in the right momentum/angular range.

parent 25b5c216
No related branches found
No related tags found
No related merge requests found
This commit is part of merge request !195. Comments created here will be created in the context of that merge request.
......@@ -17,8 +17,12 @@ using namespace HepMC3;
void gen_far_forward_protons(int n_events = 100,
const char* out_fname = "far_forward_protons.hepmc")
{
double cos_theta_min = std::cos(0.3*(M_PI/180.0));
double cos_theta_max = std::cos(2.0*(M_PI/180.0));
double minMomentum = 270.0; //GeV
double maxMomentum = 275.0; //GeV
double cos_theta_min = std::cos(0.0*(M_PI/180.0)); //we need MRAD for the units - syntax is degrees*Pi/180
double cos_theta_max = std::cos(1.15*(M_PI/180.0)); //FF region is 0.0 to 20.0 mrad
WriterAscii hepmc_output(out_fname);
int events_parsed = 0;
......@@ -39,7 +43,7 @@ void gen_far_forward_protons(int n_events = 100,
FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
// Define momentum
Double_t p = r1->Uniform(1.0, 10.0);
Double_t p = r1->Uniform(minMomentum, maxMomentum);
Double_t phi = r1->Uniform(0.0, 2.0 * M_PI);
Double_t costh = r1->Uniform(cos_theta_min, cos_theta_max);
Double_t th = std::acos(costh);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment