Vanwdani at 02:38, 23 January 2018

2018-01-23T02:38:56Z

Vanwdani: Created page with "

  #include   #include   #include   #include    void LUND2ParticleGun() {         int nlines,file_parts,step;         char LUND_file..."

2018-01-17T15:56:13Z

Created page with "<pre> #include <math.h> #include <TRandom3.h> #include <stdio.h> #include <stdlib.h> void LUND2ParticleGun() { int nlines,file_parts,step; char LUND_file..."

New page

<pre>
#include <math.h>
#include <TRandom3.h>
#include <stdio.h>
#include <stdlib.h>

void LUND2ParticleGun()
{
int nlines,file_parts,step;
char LUND_file[200];
char path[200];
char Run_file[20];
char mac[20]="Run";
double Px,Py,Pz,E;
double px,py,pz,e;
double weight;

double M,X,Y,Z;
double m,x,y,z;

double Number_of_particles,Number_of_target_nucleons,Number_of_target_protons,Target_polarization,Beam_polarization,W,Q2,v;
double Index1,Charge1,Type1,PID1,Parent_Index1,Daughter_Index1;
double Index2,Charge2,Type2,PID2,Parent_Index2,Daughter_Index2;

double extra;

cout << "Number of LUND files to convert?" << endl;

cin >> file_parts ;

cout << "Path?" << endl;

cin >> path;
step=0;
for(int i=0;i<file_parts;i++)
{
ifstream in;
if(step<10)
{
sprintf(LUND_file,"%s_00%d.LUND",path,step);
in.open(LUND_file);
}
if(step>9 && step<100)
{
sprintf(LUND_file,"%s_0%d.LUND",path,step);
in.open(LUND_file);
}
if(step>99 && step<1000)
{
sprintf(LUND_file,"%s_%d.LUND",path,step);
in.open(LUND_file);
}
ofstream out;
if(step<10)
{
sprintf(Run_file,"%s00%d.mac",mac,step);
out.open(Run_file);
}
if(step>9 && step<100)
{
sprintf(Run_file,"%s0%d.mac",mac,step);
out.open(Run_file);
}
if(step>99 && step<1000)
{
sprintf(Run_file,"%s%d.mac",mac,step);
out.open(Run_file);
}

out << "random/setSeeds 12345679" << endl;
for(nlines=0;nlines<1000;nlines++)
{
//Read in LUND variables
in >> Number_of_particles >> Number_of_target_nucleons >> Number_of_target_protons >> extra >> Target_polarization >> Beam_polarization >> weight >>
y >> W >> Q2 >> v ;

in >> Index1 >> Charge1 >> Type1 >> PID1 >> Parent_Index1 >> Daughter_Index1 >> Px >> Py >> Pz >> E >> M >> X >> Y >> Z;

in >> Index2 >> Charge2 >> Type2 >> PID2 >> Parent_Index2 >> Daughter_Index2 >> px >> py >> pz >> e >> m >> x >> y >> z;

//Write particle gun macro
out << "/event/verbose 1" << endl;

out << "/tracking/verbose 1" << endl;

out << "/gun/particle e-" << endl;

out << "/gun/energy "<< e << endl;

out << "/gun/position " << x << "\t" << y << "\t" << z << "\t" << "cm" << endl;

out << "/gun/momentum " << px << "\t" << py << "\t" << pz << "\t" << "MeV" << endl;

out << "/run/beamOn 1" << endl;

}//End individual LUND loop
step=step+1;
in.close();
out.close();
}//End loop over all LUNDS
}//End main

LUND2ParticleGun.C - Revision history

Vanwdani at 02:38, 23 January 2018

Vanwdani: Created page with "
#include #include #include #include void LUND2ParticleGun() { int nlines,file_parts,step; char LUND_file..."

@@ Line 70: / Line 70: @@
+                 }
-                 out << "random/setSeeds 12345679" << endl;
+                 out << "/random/setSeeds 12345679 0" << //1st numberis the initial random seed and 2nd is the number of times a random seed is sampled before starting the simulation  endl;
                  for(nlines=0;nlines<1000;nlines++)
+                 {
@@ Line 89: / Line 89: @@
                          out << "/gun/particle e-" << endl;
                          out << "/gun/position " << x << "\t" << y << "\t" << z << "\t" << "cm" << endl;