DV Creating LUND Files - Revision history

Vanwdani: /* The LUND format */

2017-09-08T18:47:42Z

The LUND format

Vanwdani: /* The LUND format */

2017-09-08T17:30:40Z

The LUND format

Vanwdani: /* The LUND format */

2017-09-08T16:24:47Z

The LUND format

Vanwdani: /* The LUND format */

2017-09-08T16:04:55Z

The LUND format

Vanwdani at 15:35, 8 September 2017

2017-09-08T15:35:03Z

Vanwdani at 19:54, 5 April 2016

2016-04-05T19:54:00Z

Vanwdani: /* Corrections */

2016-03-30T15:56:00Z

Corrections

Vanwdani: /* Corrections */

2016-03-30T15:54:24Z

Corrections

Vanwdani: /* Writing LUND files */

2016-03-30T15:48:37Z

Writing LUND files

Vanwdani: /* References */

2016-03-29T23:39:58Z

References

@@ Line 23: / Line 23: @@
 GEMC 2.4 had been taking the 5th element on the first line of each LUND event, the beam polization for the 5th entry in the header array (10,5).  Newer GEMC is grabbing the 6th element of the first line of each LUND event.
-This can be fixed by ensuring that each LUND file starts with a blank line before the first LUND event.
+This can be fixed by adding an additional column with value 1 every 3rd row to the LUND file.
 <pre>

← Older revision		Revision as of 17:30, 8 September 2017
Line 24:		Line 24:

	This can be fixed by ensuring that each LUND file starts with a blank line before the first LUND event.		This can be fixed by ensuring that each LUND file starts with a blank line before the first LUND event.
		+
		+	<pre>
		+	sed -i -e '1i\ \' LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut/LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut_.LUND
		+	</pre>

	[[File:HeaderGEMC2_6fix.png]]		[[File:HeaderGEMC2_6fix.png]]

@@ Line 22: / Line 22: @@
 GEMC 2.4 had been taking the 5th element on the first line of each LUND event, the beam polization for the 5th entry in the header array (10,5).  Newer GEMC is grabbing the 6th element of the first line of each LUND event.
 ==The Header==
 {| border=1 align=center

@@ Line 4: / Line 4: @@
 The LUND format has extra variables that are not utilized within GEMC.  Only the '''BOLD''' variables are necessary within GEMC simulations.[1]
-NOTE: Newer versions of GEMC use different format header
+NOTE: Newer versions of GEMC use different format header.  GEMC 2.4 and 2.6 are run for the same 10 events under the same circumstances with identical gcard files.  The ouput is displayed in text format instead of evio.
 [[File:HeaderGEMC2_4.png]]
 [[File:HeaderLUND.png]]
 [[File:HeaderGEMC2_6.png]]
 ==The Header==
 {| border=1 align=center

← Older revision		Revision as of 15:35, 8 September 2017
Line 3:		Line 3:

	The LUND format has extra variables that are not utilized within GEMC. Only the '''BOLD''' variables are necessary within GEMC simulations.[1]		The LUND format has extra variables that are not utilized within GEMC. Only the '''BOLD''' variables are necessary within GEMC simulations.[1]
		+
		+	NOTE: Newer versions of GEMC use different format header
		+
		+	[[File:HeaderGEMC2_4.png]]
		+
		+
		+	[[File:HeaderLUND.png]]
		+
		+
		+	[[File:HeaderGEMC2_6.png]]

	==The Header==		==The Header==

← Older revision		Revision as of 19:54, 5 April 2016
Line 203:		Line 203:
	===14:Vertex z===		===14:Vertex z===
	This represents the z position within the eg12 detector at time t=0.		This represents the z position within the eg12 detector at time t=0.
−
−	~~=Writing LUND files=~~
−	~~Reading the GEANT4 simulation output in,~~
−	~~<pre>~~
−	~~void MollerLUND()~~
−	{
−	~~struct evt_t~~
−	{
−	~~Long_t event;~~
−	~~Float_t IntKE, IntMom[3],IntPos[3],FnlKE,FnlMom[3],FnlPos[3],MolKE,MolMom[3],MolPos[3];~~
−	};
−	~~ifstream in;~~
−	~~in.open("MollerScattering_NH3_4e8incident.dat");~~
−	~~evt_t evt;~~
−	~~Int_t nlines=0;~~
−	~~TFile *f = new TFile("MollerScattering_NH3_4e7incident.root","RECREATE");~~
−	~~TTree *tree0 = new TTree("MollerScattering_NH3_4e7","Moller data from ascii file");~~
−	~~FILE *f0;~~
−	~~f0=fopen("MollerScattering_NH3_4e7.LUND","w");~~
−
−	~~tree0->Branch("evt",&evt.event,"event/I:IntKE/F:IntPx:IntPy:IntPz:IntPosx:IntPosy:IntPosz:FnlKE:FnlPx:FnlPy:FnlPz:FnlPosx:FnlPosy:FnlPosz:~~
−	~~MolKE:MolPx:MolPy:MolPz:MolPosx:MolPosy:MolPosz");~~
−
−	~~while(in.good())~~
−	{
−	~~evt.event=nlines;~~
−	~~in >> evt.IntKE >> evt.IntMom[0] >> evt.IntMom[1] >> evt.IntMom[2] >> evt.IntPos[0]~~
−	~~>> evt.IntPos[1] >> evt.IntPos[2] >> evt.FnlKE >> evt.FnlMom[0] >> evt.FnlMom[1] >> evt.FnlMom[2]~~
−	~~>> evt.FnlPos[0] >> evt.FnlPos[1] >> evt.FnlPos[2] >> evt.MolKE >> evt.MolMom[0] >> evt.MolMom[1]~~
−	~~>> evt.MolMom[2] >> evt.MolPos[0] >> evt.MolPos[1] >> evt.MolPos[2]; double Nu;~~
−	~~double x_bj; // Bjorken-x~~
−	~~double y;~~
−	~~double W ; // missing mass~~
−
−	~~// Momentum and Energy are measured in MeV!~~
−
−	~~evt.FnlKE=sqrt(evt.FnlMom[0]evt.FnlMom[0]+evt.FnlMom[1]evt.FnlMom[1]+evt.FnlMom[2]evt.FnlMom[2]+0.5110.511);~~
−	~~evt.MolKE=sqrt(evt.MolMom[0]evt.MolMom[0]+evt.MolMom[1]evt.MolMom[1]+evt.MolMom[2]evt.MolMom[2]+0.5110.511);~~
−
−
−	~~Nu=evt.IntKE-evt.FnlKE;~~
−	~~Qsqrd=4evt.IntKEevt.FnlKE(1-evt.FnlMom[2]/evt.FnlKE)/2; / should be final momentum and not final energy*/~~
−
−	~~x_bj=Qsqrd/(20.938Nu);~~
−	~~y=Nu/evt.IntKE;~~
−	~~W=0.9380.938+20.938*Nu-Qsqrd;~~
−	~~if(W>0)~~
−	~~W=sqrt(W);~~
−
−	//
−	~~// MeV MUST be converted to GeV for LUND format!~~
−
−	~~double px,py,pz;~~
−	~~double Px,Py,Pz;~~
−	~~double KE,ke;~~
−
−	~~Px=evt.FnlMom[0]/10;~~
−	~~Py=evt.FnlMom[1]/10;~~
−	~~Pz=evt.FnlMom[2]/10;~~
−	~~px=evt.MolMom[0]/100;~~
−	~~py=evt.MolMom[1]/100;~~
−	~~pz=evt.MolMom[2]/100;~~
−
−	~~KE=evt.FnlKE/10;~~
−	~~ke=evt.MolKE/100;~~
−
−
−	~~fprintf(f0,"%d\t%g\t%g\t%g\t%g\t%g\t%g\t%g\t%g\t%g\n",2,1.,1.,1.,1.,x_bj,y,W,Qsqrd,Nu);~~
−	~~fprintf(f0, "%d\t%g\t%g\t%g\t%g\t%g\t%f\t%f\t%f\t%f\t%g\t%g\t%g\t%g\n",~~
−	~~1,-1.,1.,11.,0.,0.,Px,Py,Pz,KE,0.000511, 0.,0.,0.);~~
−	~~fprintf(f0, "%d\t%g\t%g\t%g\t%g\t%g\t%f\t%f\t%f\t%f\t%g\t%g\t%g\t%g\n",~~
−	~~2,-1.,1.,11.,0.,0.,px,py,pz,ke,0.000511, 0.,0.,0.);~~
−
−
−	~~nlines++;~~
−	~~tree0->Fill();~~
−	}
−	~~printf("NumberLines=%d\n",nlines);~~
−	~~tree0->Print();~~
−	~~tree0->Write();~~
−	~~delete tree0;~~
−
−	}
−	~~</pre>~~
−
−
−	~~The file is written as seen below.~~
−
−	~~[[File:LUND_view.png]]~~
−
−	~~=Corrections=~~
−	The prefered coordinates for the reconstruction within GEMC would have the vertex location at (0, 0, 0). The Moller event data from GEANT4 contains the location within that simulation where the scattering occured. These variables will need to be renormalized within the LUND file to zero. Only the 14th column within the 2nd and 3rd rows out of every 3 rows needs to be changed. Using a script:
−
−
−	~~<pre>awk '{if('NR%3==2'\|\|'NR%3==0') {$14="0"}; print $0}' MollerScattering_NH3_4e7.LUND>MollerScattering.LUND</pre>~~
−
−	~~The LUND format has a limit of 75000 lines, so to account for this by creating 21 new files with the suffix of aa-au are created:~~
−
−	~~<pre>split -l 75000 MollerScattering.LUND MollerScattering_NH3_4e7_ </pre>~~
−
−	~~We will use 5 of the 21 sets of 75000 entries in GEMC trials involving different solenoid field strengths. We can use the simple command:~~
−
−	~~<pre> mv MollerScattering_NH3_4e7_aa MollerScattering_NH3_4e7_aa.LUND</pre>~~
−
−	~~Repeating this for the other 4 files.~~
−
−	~~=References=~~
−	~~[1] M. Ungaro, ''et al.,'' Jefferson Lab:The LUND format, https://gemc.jlab.org/gemc/html/documentation/generator/lund.html (2016)~~
−
−	~~[2] F. Halzen, and A. Martin, ''Quarks and Leptons:An Introductory Course in Modern Particle Physics,'' John Wiley & Sons (1984)~~
−
−	~~[3] W. Tung, ''Bjorken scaling,'' http://www.scholarpedia.org/article/Bjorken_scaling (2009)~~
−
−	~~[4] L. Garren, ''et al.,'' MONTE CARLO PARTICLE NUMBERING SCHEME, http://pdg.lbl.gov/2007/reviews/montecarlorpp.pdf (2007)~~

← Older revision		Revision as of 15:56, 30 March 2016
Line 306:		Line 306:

	<pre> mv MollerScattering_NH3_4e7_aa MollerScattering_NH3_4e7_aa.LUND</pre>		<pre> mv MollerScattering_NH3_4e7_aa MollerScattering_NH3_4e7_aa.LUND</pre>
		+
		+	Repeating this for the other 4 files.

	=References=		=References=

@@ Line 297: / Line 297: @@
-<pre>awk '{if('NR%3==2'||'NR%3==0') {$14="0"}; print $0}' MollerScattering_NH3_4e8.LUND>MollerScattering.LUND</pre>
+<pre>awk '{if('NR%3==2'||'NR%3==0') {$14="0"}; print $0}' MollerScattering_NH3_4e7.LUND>MollerScattering.LUND</pre>
-The LUND format has a limit of 75000 lines, so to account for this
+The LUND format has a limit of 75000 lines, so to account for this by creating 21 new files with the suffix of aa-au are created:
-<pre>split -l 75000 MollerScattering.LUND MollerScattering_NH3_4e9 </pre>
+<pre>split -l 75000 MollerScattering.LUND MollerScattering_NH3_4e7_ </pre>
 =References=

@@ Line 218: / Line 218: @@
          evt_t evt;
          Int_t nlines=0;
-         TFile *f = new TFile("MollerScattering_NH3_4e8incident.root","RECREATE");
+         TFile *f = new TFile("MollerScattering_NH3_4e7incident.root","RECREATE");
-         TTree *tree0 = new TTree("MollerScattering_NH3_4e8","Moller data from ascii file");
+         TTree *tree0 = new TTree("MollerScattering_NH3_4e7","Moller data from ascii file");
          FILE *f0;
-         f0=fopen("MollerScattering_NH3_4e8.LUND","w");
+         f0=fopen("MollerScattering_NH3_4e7.LUND","w");
          tree0->Branch("evt",&evt.event,"event/I:IntKE/F:IntPx:IntPy:IntPz:IntPosx:IntPosy:IntPosz:FnlKE:FnlPx:FnlPy:FnlPz:FnlPosx:FnlPosy:FnlPosz:

@@ Line 315: / Line 315: @@
 ----
-Links
+=Links=
 [[DV_RunGroupC_Moller#Creating_LUND_Files|Back]]