Difference between revisions of "DV Creating LUND Files"

From New IAC Wiki
Jump to navigation Jump to search
Line 96: Line 96:
 
where E<sub>i</sub> and E<sub>f</sup> are the initial and final electron energies.
 
where E<sub>i</sub> and E<sub>f</sup> are the initial and final electron energies.
  
=Writing LUND files=
+
=Writing LUND files{| border=1 align=center
 +
  |+ LUND Header
 +
|-
 +
  ! Column
 +
  ! Quantity
 +
|-
 +
  | <center>1</center>
 +
  | <center>Index</center>
 +
|-
 +
  |<center>2</center>
 +
  | <center>Charge</center>
 +
|-
 +
  | <center>3</center>
 +
  | <center>'''Type (1 is active)'''</center>
 +
|-
 +
  | <center>4</center>
 +
  | <center>'''Particle ID'''</center>
 +
|-
 +
  | <center>5</center>
 +
  | <center>Parent Index</center>
 +
|-
 +
  | <center>6</center>
 +
  | <center>Daughter Index</center>
 +
|-
 +
  | <center>7</center>
 +
  | <center>'''Momentum x (GeV)'''</center>
 +
|-
 +
  | <center>8</center>
 +
  | <center>'''Momentum y (GeV)'''</center>
 +
|-
 +
  | <center>9</center>
 +
  | <center>'''Momentum z (GeV)'''</center>
 +
|-
 +
  | <center>10</center>
 +
  | <center><math>E</math></center>
 +
|-
 +
  | <center>11</center>
 +
  | <center>Mass</center>
 +
|-
 +
  | <center>12</center>
 +
  | <center>'''Vertex x (cm)'''</center>
 +
|-
 +
  | <center>13</center>
 +
  | <center>'''Vertex y (cm)'''</center>
 +
|-
 +
  | <center>14</center>
 +
  | <center>'''Vertex z (cm)'''</center>
 +
|}
  
 
=Correct for vertex location=
 
=Correct for vertex location=

Revision as of 19:12, 29 March 2016

The LUND format

LUND Header
Column Quantity
1
Number of Particles
2
Number of Target Nucleons
3
Number of Target Protons
4
Target Polarization
5
Beam Polarization
6
x
7
y
8
W
9
[math]Q^2[/math]
10
[math]\nu[/math]

Where

1:Number of Particles

This line tells how many particles follow the header line. For Moller Scattering, this number should always be 2.

2:Number of Target Nucleons

For this simulation, only an electron-electron collision is considered. This quantity is always set to 1 for the one stationary electron we consider the incident electron scattering from.

3:Number of Target Protons

For Moller Scattering, there are no target protons. This number is set to 1, but does not have any effect within the GEMC simulations.

4:Target Polarization

This represents the polarization of the target material, either positive or negative 1. This value is always set to 1 and has no effect within the GEMC simulations.

5:Beam Polarization

This represents the polarization of the electron beam, either negative or positive 1. This value is always set to positive 1.

6:x

This represents the Bjorken x scaling variable.

[math]x=\frac{-q^2}{2p\cdot q}=\frac{-q^2}{2M\nu}=\frac{Q^2}{2M\nu}[/math]


Where M is the rest mass-energy of the proton[math] M\approx 938MeV[/math]

7:y

[math]y=\frac{p\cdot q}{p\cdot k}=\frac{\nu}{E_i}=\frac{E_i-E_f}{E_i}[/math]

8:W

The invariant mass of the final hadronic system


[math]W^2\equiv (p+q)^2=M^2+2M\nu+q^2=M^2+2M\nu-Q^2[/math]

9:[math]Q^2[/math]

This represents the squared 4-momentum-transfer vector q of the exchanged virtual photon.


[math]Q^2\equiv -q^2[/math]


Where q is the momentum transfer betwwn the incident electron and target via the virtual photon.


[math]q\equiv k_i-k_f[/math]

10:[math]\nu[/math]

This represents the energy loss between scattering electrons.

[math]\nu = \frac{p\cdot q}{M}[/math]

This can be written in the Lab frame as:


[math]\nu\equiv E_i-E_f[/math]


where Ei and Ef are the initial and final electron energies.

=Writing LUND files{| border=1 align=center

 |+ LUND Header

|-

 ! Column
 ! Quantity

|-

|

1

|

Index

|-

|

2

|

Charge

|-

|

3

|

Type (1 is active)

|-

|

4

|

Particle ID

|-

|

5

|

Parent Index

|-

|

6

|

Daughter Index

|-

|

7

|

Momentum x (GeV)

|-

|

8

|

Momentum y (GeV)

|-

|

9

|

Momentum z (GeV)

|-

|

10

|

[math]E[/math]

|-

|

11

|

Mass

|-

|

12

|

Vertex x (cm)

|-

|

13

|

Vertex y (cm)

|-

|

14

|

Vertex z (cm)

|}

Correct for vertex location