[E16-e1f] lund format output for Gagik's simulater

[Harut Avakian]

Dear All,
Here are the steps on using Gagik's framework for fast simulation:

1) We provide  a file with candidates of ep->e'p'\pi+\pi- with all possible combinations
ep->e'p'\pi+X,ep->e'p'\pi-X,ep->e'\pi+\pi-X with mass of the X in the right place for exclusivity (this is needed 1 for e16 and 1 for e1f).

2) We are providing N MC generated   files  with our physics (ep->e'X,ep->\ep+X,ep\gamma,... any other events also in LUND format). This can be for every ongoing process.

->From "fast simulator" we expect output of reconstructed N files with given acceptance in the form of bos files (like data).


Please use the gsim lund format for each event (example below):

  1-line N_tracks, target_Z,Q^2,nu,heliciti-prot,helicity-lepton (N_track important other not critical, can put 0 if not available )
  
  For next lines we have the track infos for those N_tracks, so for i-th track:
>                   write(41,110) istat, lund_id, parent_lund_id, line_number_for_part, line_number_for parent
>                   write(41,112) (p(i,j),j=1,5)
>                   write(41,112) (V(i,j),j=1,5)
>
The format is general, but we will not need last 2 values (show the number of lines from original LUND generated). If not available you can put
the parent_id also 0
The istat is 1 for physical particles and other number for decaying. For the sample ep\pi+\pi- all will have 1, the lund_id is the lund_id of the
track (11,211,-211,...),  parent_lund_id, is the lund_id of the parent (ex. Delta's if pi+ is coming from \Delta decay) if not defined put 0.
basically for the 1st line

1, lund_id, parent_lund_id, 0 0

second line are momenta
> -------
> c P(I,1) $p_x$, momentum in the $x$ direction, in GeV/$c$.
> c P(I,2) $p_y$, momentum in the $y$ direction, in GeV/$c$.
> c P(I,3) $p_z$, momentum in the $z$ direction, in GeV/$c$.
> c P(I,4) $E$, energy, in GeV.
> c P(I,5) $m$, mass, in GeV/$c^2$.

third line are vertexes (you can put 0 for   V(I,4) and V(I,5) for now.
> c V(I,1) $x$ position of production vertex, in cm.
> c V(I,2) $y$ position of production vertex, in cm.
> c V(I,3) $z$ position of production vertex, in cm.
> c V(I,4) time of production, in mm/$c$ ($\approx 3.33 \times 10^{-12}$ s).
> c V(I,5) proper lifetime of particle, in mm/$c$
> The formats are for first line
> 110           FORMAT(2x,10I10)
> 112           FORMAT(2x,10F10.4)
>
> Example of gsim LUND format entry
> -----
>      4    1.0000    1.3165    2.9650    0.0000    1.0000
>           1        11        11         4         1
>     -0.1123   -0.7612    2.4124    2.5322    0.0005
>     -0.3618    0.4774  -26.8623    0.0000    0.0000
>           1       321        92        13        12
>      0.3841    0.0127    2.6674    2.7398    0.4936
>     -0.3618    0.4774  -26.8623    0.0000    0.0000
>           1      2212      3122        15        14
>     -0.1295    0.6146    0.6244    1.2903    0.9383
>     -0.5216    1.1491  -26.1124    0.0000    0.0000
>           1      -211      3122        16        14
>     -0.0593    0.1789    0.2615    0.3513    0.1396
>     -0.5216    1.1491  -26.1124    0.0000    0.0000
> =======================================================


Let me know if you have any problems.
Best,
Harut