Difference between revisions of "Running Expanded LUND files in GEMC"
Jump to navigation
Jump to search
| Line 1: | Line 1: | ||
| − | Since the LUND file is limited to 75000 lines, the gemc will have to be run in batch mode; | + | Since the LUND file is limited to 75000 particles (225000 lines), the gemc will have to be run in batch mode; |
Creating a batch directory, with two subdirectories; 1)Phi_Parts, 2)submit. | Creating a batch directory, with two subdirectories; 1)Phi_Parts, 2)submit. | ||
| Line 88: | Line 88: | ||
gemc -USE_GUI=0 -Hall_Material="Vacuum" -INPUT_GEN_FILE="LUND, Phi_Parts_0000.LUND" -N=75000 eg12.gcard | gemc -USE_GUI=0 -Hall_Material="Vacuum" -INPUT_GEN_FILE="LUND, Phi_Parts_0000.LUND" -N=75000 eg12.gcard | ||
~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -i eg12.ev -config DCHB::torus=1.0 -config DCHB::solenoid=0.0 -config DCTB::kalman=true -o eg12_rec.ev -s DCHB:DCTB:EC:FTOF:EB | ~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -i eg12.ev -config DCHB::torus=1.0 -config DCHB::solenoid=0.0 -config DCTB::kalman=true -o eg12_rec.ev -s DCHB:DCTB:EC:FTOF:EB | ||
| − | ~/src/CLAS/coatjava-1.0/bin/ | + | ~/src/CLAS/coatjava-1.0/bin/run-groovy Analysis.groovy eg12_rec.0.evio |
</pre> | </pre> | ||
Latest revision as of 16:28, 6 June 2016
Since the LUND file is limited to 75000 particles (225000 lines), the gemc will have to be run in batch mode;
Creating a batch directory, with two subdirectories; 1)Phi_Parts, 2)submit.
1)Once the LUND file is broken into 178 parts, they can have the LUND extension added by:
prename 's/(Phi_Parts_\d{4})/$1.LUND/' Phi_Parts_*
Placing each of these files into its own directory, within a directory named Phi_Parts
find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;
2)Creating the submit directory, and using a c++ program, creating the needed 178
#include <iomanip>
#include <sstream>
#include <iostream>
#include <fstream>
using namespace std;
void submit() {
for(int a=0;a<2;a++)
{
for(int b=0;b<10;b++)
{
for(int c=0;c<10;c++)
{
string filename="submit0";
stringstream hundreds;
hundreds << a;
stringstream tens;
tens << b;
stringstream ones;
ones << c;
string fullname="";
fullname=filename + hundreds.str() + tens.str() + ones.str();
// cout << fullname << "\n";
ofstream myfile;
myfile.open(fullname.c_str());
myfile << "#!/bin/sh\n";
myfile << "#PBS -l nodes=1\n";
myfile << "#PBS -A FIAC\n";
myfile << "#PBS -M vanwdani@isu.edu\n";
myfile << "#PBS -m abe\n";
myfile << "#\n";
myfile << "cd /home/lds/src/CLAS/GEMC\n";
myfile << "tcsh\n";
myfile << "source setup\n";
myfile << "cd /home/lds/src/GEANT/geant4.9.6/geant4.9.6-install/bin/geant4.sh\n";
myfile << "cd /home/vanwdani/src/GEANT4/geant4.9.6/Simulations/Research/Moller/batch/Phi_Parts/Phi_Parts_0";
myfile <<a<<b<<c<<"\n";
myfile << "gemc -USE_GUI=0 -Hall_Material=\"Vacuum\" -INPUT_GEN_FILE=\"LUND, Phi_Parts_0";
myfile <<a<<b<<c;
myfile << ".LUND\" -N=75000 eg12.gcard\n";
myfile << "~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -i eg12.ev -config DCHB::torus=1.0 ";
myfile << "-config DCHB::solenoid=0.0 -config DCTB::kalman=true -o eg12_rec.ev -s DCHB:DCTB:EC:FTOF:EB\n";
myfile << "~/src/CLAS/coatjava-1.0/bin/run-groovy Analysis.groovy eg12_rec.0.evio\n";
myfile.close();
}
}
}
}
This creates the submitXXXX file
#!/bin/sh #PBS -l nodes=1 #PBS -A FIAC #PBS -M vanwdani@isu.edu #PBS -m abe # cd /home/lds/src/CLAS/GEMC tcsh source setup cd /home/lds/src/GEANT/geant4.9.6/geant4.9.6-install/bin/geant4.sh cd /home/vanwdani/src/GEANT4/geant4.9.6/Simulations/Research/Moller/batch/Phi_Parts/Phi_Parts_0000 gemc -USE_GUI=0 -Hall_Material="Vacuum" -INPUT_GEN_FILE="LUND, Phi_Parts_0000.LUND" -N=75000 eg12.gcard ~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -i eg12.ev -config DCHB::torus=1.0 -config DCHB::solenoid=0.0 -config DCTB::kalman=true -o eg12_rec.ev -s DCHB:DCTB:EC:FTOF:EB ~/src/CLAS/coatjava-1.0/bin/run-groovy Analysis.groovy eg12_rec.0.evio
Creating a file named lds-submit
