Fix hard coded beam energies in various places

cbravo135 · May 2, 2018 · cc6066e · cc6066e
1 parent 8e3eaf0
commit cc6066e
Show file tree

Hide file tree

Showing 3 changed files with 8 additions and 8 deletions.
diff --git a/include/LHEWriter.h b/include/LHEWriter.h
@@ -13,7 +13,7 @@ class LHEWriter
         ofstream oF;
 
     public:
-        LHEWriter(string fName);
+        LHEWriter(string fName, double sqrtS);
         ~LHEWriter();
 
         int writeEvent(vector<particle> outParts, double Q);

diff --git a/src/BaryoGEN.cc b/src/BaryoGEN.cc
@@ -94,7 +94,7 @@ int main(int argc, char* argv[])
 
     //Init output files
     TFile *myF = new TFile((ofName+".root").c_str(),"RECREATE","Holds daughters from sphaleron decay");
-    LHEWriter lheF(ofName);
+    LHEWriter lheF(ofName, SQRTS);
 
     double minx = thr*thr/(SQRTS*SQRTS);
 
@@ -213,12 +213,12 @@ int main(int argc, char* argv[])
         particle partBuf = partBase->getParticle(iq1);
         if(iq1 != partBuf.pid) cout << "iq1 = " << iq1 << " != partBuf.pid = " << partBuf.pid << endl;
         partBuf.color = 501;
-        partBuf.p4v.SetXYZM(0.0,0.0,x1*6500,partBuf.mass);
+        partBuf.p4v.SetXYZM(0.0,0.0,x1*SQRTS/2.0,partBuf.mass);
         inParts.push_back(partBuf); //Push first incoming quark
 
         partBuf = partBase->getParticle(iq2);
         partBuf.color = 502;
-        partBuf.p4v.SetXYZM(0.0,0.0,x2*(-6500),partBuf.mass);
+        partBuf.p4v.SetXYZM(0.0,0.0,x2*SQRTS/-2.0,partBuf.mass);
         inParts.push_back(partBuf); //Push second incoming quark
 
         int colNow = 503;//This is used to keep track of the color line numbers already used
@@ -237,8 +237,8 @@ int main(int argc, char* argv[])
             masses[i] = confBuf[i].mass;
         }
 
-        u1.SetXYZM(0.0,0.0,x1*6500,UQ_MASS);
-        u2.SetXYZM(0.0,0.0,x2*(-6500),UQ_MASS);
+        u1.SetXYZM(0.0,0.0,x1*SQRTS/2.0,inParts[0].mass);
+        u2.SetXYZM(0.0,0.0,x2*SQRTS/-2.0,inParts[1].mass);
         mom = u1 + u2;
         momM = mom.M();
         pz = mom.Pz();

diff --git a/src/LHEWriter.cc b/src/LHEWriter.cc
@@ -3,15 +3,15 @@
 #include <math.h>
 using namespace std;
 
-LHEWriter::LHEWriter(string fName)
+LHEWriter::LHEWriter(string fName, double sqrtS)
 {
     oF.open((fName+".lhe").c_str());
     oF << std::scientific;
     oF << "<LesHouchesEvents version=\"1.0\">" << endl;
     oF << "<header>" << endl;
     oF << "</header>" << endl;
     oF << "<init>" << endl;
-    oF << "\t2212 2212 0.65000000000e+04 0.65000000000e+04 292200 292200 292200 292200 3 1" << endl;//beam conditions information
+    oF << "\t2212 2212 " << sqrtS/2.0 << " " << sqrtS/2.0 << " 292200 292200 292200 292200 3 1" << endl;//beam conditions information
     oF << "\t7.3 1.0 1.0 7000" << endl;//sphaleron process information
     oF << "</init>" << endl;
 }