Re: pions & protons

[Eugene Chudakov <gen@jlab.org>]

Hi,

not knowing if there is a public repository to store various 
simulation results, I made one on the work disk at JLab, in
a directory mc, made by Alex Somov:

/work/halld/mc/kin
/work/halld/mc/kin/bggen

ls -lF /work/halld/mc/kin/bggen 
-rw-r--r-- 1 gen halld 239 2009-04-30 18:08 README
drwxr-sr-x 2 gen halld  11 2009-04-30 18:04 v1

The directory v1 is a copy of a directory where I ran bggen
for 395000 events (1 sec of experiment running at high luminosity), 
for photons from 0.15 to 12. GeV.

ls -lF /work/halld/mc/kin/bggen/v1
-rw-r--r-- 1 gen a-staff 131062269 2008-01-17 20:22 bggen.hddm
-rw-r--r-- 1 gen a-staff    827392 2008-03-12 14:00 bggen.his
-rw-r--r-- 1 gen a-staff 101679104 2008-03-12 14:00 bggen.nt
lrwxrwxrwx 1 gen a-staff         7 2008-01-18 10:15 fort.15 -> run.ffr
-rw-r--r-- 1 gen a-staff     73088 2008-03-12 14:00 log
-rw-r--r-- 1 gen a-staff       758 2008-01-17 20:12 particle.dat
-rw-r--r-- 1 gen a-staff      1057 2008-01-17 20:12 pythia.dat
-rw-r--r-- 1 gen a-staff      1104 2008-01-17 20:12 pythia-geant.map
-rw-r--r-- 1 gen a-staff       676 2008-03-12 13:54 run.ffr

In order to analyze the ntuple bggen.nt one can use the directory bggen/paw
from svn:
src/programs/Simulation/bggen/paw
which contains several useful functions.

cd paw/
paw
PAW>h/fil 2 bggen.nt 0 -x
PAW>nt/lis //lun2
PAW>nt/pri //lun2/9
PAW>2dh 302 'p-th' 180 0 180 120 0 12
PAW>nt/pl //lun2/9.part_kin.f(2,8,0,302)
PAW>h/pl 302

The histogram 302 contains a theta-p plot for all the pi+ particles
in the ntuple (part_kin.f contains a brief instruction). One needs
to replace 8 to 9 to see the pi-, or 14 to see the protons. In order
to see an decaying particle, say Delta++ one can use:
PAW>nt/pl //lun2/9.part_kin.f(2,0,2224,302)


PAW>nt/loop //lun2/9 8<pin(3,1)<9.and.bgg_pri.f(0)>-1 100 1

will print a few first events in the photon beam range from 
8 to 9 GeV.

Eugene

On Thu, 30 Apr 2009, Zisis Papandreou wrote:

> Hello folks:
>
> In connection to the simulations we are doing in Regina, can you please 
> confirm/provide the following:
>
> - Eugene, the ntuple.  Can you please tell us how to fetch it?  Email is 
> small, or else please put it on one of the HallD disks and let us know where.
> - Elton/Beni:  I found GlueX-doc-999 on the DocDB portal.  There are also the 
> final talks for the PID Review (Eugene's and Curtis').  Anything else 
> pertinent to pion/proton yields?
>
> Thanks, Zisis...
>
> ---
> Dr. Zisis Papandreou		|  email: zisis@uregina.ca
> Department of Physics	|  tel. : (306) 585-5379
> University of Regina		|  tel. : (306) 585-4149
> Regina, SK  S4S 0A2		|  fax. : (306) 585-5659
> World Wide Web:    http://www.phys.uregina.ca/sparro/zisis/
>
>
>