Run in GEMC - Revision history

Vanwdani at 19:06, 14 August 2017

2017-08-14T19:06:33Z

Vanwdani: /* Simulating with GEMC */

2017-08-14T19:04:01Z

Simulating with GEMC

Vanwdani: /* Simulating with GEMC */

2017-08-14T19:01:20Z

Simulating with GEMC

Vanwdani: /* Simulating with GEMC */

2017-08-14T18:59:07Z

Simulating with GEMC

Vanwdani: /* Simulating with GEMC */

2017-08-14T18:55:02Z

Simulating with GEMC

Vanwdani: /* Simulating with GEMC */

2017-08-14T18:46:53Z

Simulating with GEMC

Vanwdani at 18:39, 14 August 2017

2017-08-14T18:39:50Z

Vanwdani: /* Simulating with GEMC */

2016-08-04T22:32:39Z

Simulating with GEMC

Vanwdani: /* Simulating with GEMC */

2016-08-04T22:32:08Z

Simulating with GEMC

Vanwdani: /* Simulating with GEMC */

2016-08-04T16:46:52Z

Simulating with GEMC

← Older revision		Revision as of 19:06, 14 August 2017
Line 1:		Line 1:
	Since the LUND file is limited to 75000 particles (225000 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.~~
−
−
−	~~'''1)'''Once the LUND file is broken into 178 parts, they can have the LUND extension added by:~~
−
−	~~<pre>prename 's/(Phi_Parts_\d{4})/$1.LUND/' Phi_Parts_*</pre>~~
−
−	~~Placing each of these files into its own directory, within a directory named Phi_Parts~~
−
−	~~<pre>find . -name ".LUND" -exec sh -c 'mkdir "${1%.}" ; mv "$1" "${1%.*}" ' _ {} \;</pre>~~
−
−
−	~~'''2)'''Creating the submit directory, and using a c++ program, creating the needed 178~~
−
−	~~<pre>#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();~~
−	}
−	}
−	}
−
−	~~}</pre>~~
−
−
−	~~This creates the submitXXXX file~~
−
−	~~<pre>#!/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~~
−
−	~~</pre>~~
−
−
−	~~Creating a file named lds-submit~~
−
−	~~[[File:Screen_Shot_2016-03-15_at_2.36.26_PM.png]]~~

@@ Line 133: / Line 133: @@
          <detector name="../../../../../../clas12/dc/dc"            factory="TEXT" variation="ccdb"/>
          <detector name="../../../../../../clas12/ftof/ftof"            factory="TEXT" variation="java"/>
-         <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
+         <!--option name="SCALE_FIELD" value="clas12-torus-big, -1"/-->
 <!-- Solenoid apparatus and field -->

@@ Line 111: / Line 111: @@
 <pre>
-split -d -l 3000 -a 3 LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.LUND LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_
+split -d -l 3000 -a 3 LUND_Spread_LH2_IsotropicPhi.LUND LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_
 </pre>

@@ Line 111: / Line 111: @@
 <pre>
-split -d -l 3000 -a 2 LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.LUND LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_
+split -d -l 3000 -a 3 LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.LUND LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_
 </pre>
@@ Line 117: / Line 117: @@
 <pre>
-prename 's/(LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_\d{2})/$1.LUND/' LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_*
+prename 's/(LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_\d{3})/$1.LUND/' LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_*
 </pre>

← Older revision		Revision as of 18:55, 14 August 2017
Line 103:		Line 103:

	<center>[[File:Mol_E_Theta_Lab_Spread.png]]</center>		<center>[[File:Mol_E_Theta_Lab_Spread.png]]</center>
		+
		+	<pre>~/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut> wc -l LUND_Spread_LH2_IsotropicPhi.LUND
		+	315000 LUND_Spread_LH2_IsotropicPhi.LUND</pre>

@@ Line 105: / Line 105: @@
-Spliting the LUND file into 38 parts of 3000 events, or 9000 lines
+Spliting the LUND file into 105 parts of 1000 events, or 3000 lines
 <pre>
-split -d -l 9000 -a 2 LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut.LUND LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_
+split -d -l 3000 -a 2 LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.LUND LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_
 </pre>
@@ Line 114: / Line 114: @@
 <pre>
-prename 's/(LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_\d{2})/$1.LUND/' LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_*
+prename 's/(LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_\d{2})/$1.LUND/' LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_*
 </pre>
-Placing each of these files into its own directory, within a directory named LH2_0Sol_p100Tor_11GeV_Phi10deg_ShieldOut
+Placing each of these files into its own directory, within a directory named LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut
 <pre>find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;</pre>

@@ Line 1: / Line 1: @@
-==[[Running_Expanded_LUND_files_in_GEMC| Old runs]]==
+Since the LUND file is limited to 75000 particles (225000 lines), the gemc will have to be run in batch mode;
 =11GeV incident electron=
 ==Simulating with GEMC==

@@ Line 10: / Line 10: @@
 <pre>
-split -d -l 9000 -a 2 LH2_0Sol_p100Tor_11GeV_Phi10deg_ShieldOut.LUND LH2_0Sol_p100Tor_11GeV_Phi10deg_ShieldOut_
+split -d -l 9000 -a 2 LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut.LUND LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_
 </pre>
@@ Line 16: / Line 16: @@
 <pre>
-prename 's/(LH2_0Sol_p100Tor_11GeV_Phi10deg_ShieldOut_\d{2})/$1.LUND/' LH2_0Sol_p100Tor_11GeV_Phi10deg_ShieldOut_*
+prename 's/(LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_\d{2})/$1.LUND/' LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_*
 </pre>

@@ Line 7: / Line 7: @@
-Spliting the LUND file into 5 parts of 25000 events, or 75000 lines
+Spliting the LUND file into 38 parts of 3000 events, or 9000 lines
 <pre>
-split -d -l 75000 -a 2 LUND_Spread.LUND LUND_Spread_
+split -d -l 9000 -a 2 LH2_Sol0T_Tor5T_11GeV_Phi10deg.LUND LH2_Sol0T_Tor5T_11GeV_Phi10deg_
 </pre>
@@ Line 16: / Line 16: @@
 <pre>
-prename 's/(LUND_Spread_\d{2})/$1.LUND/' LUND_Spread_*
+prename 's/(LH2_Sol0T_Tor5T_11GeV_Phi10deg_\d{2})/$1.LUND/' LH2_Sol0T_Tor5T_11GeV_Phi10deg_*
 </pre>