The results shown on this page use Jas4pp data-analysis package.

wget http://atlaswww.hep.anl.gov/asc/jas4pp/download/current.php -O jas4pp.tgz
tar -zvxf jas4pp.tgz
cd jas4pp
source ./setup.sh # takes 5 sec for first-time optimization

You can study this geometry using sifcch3 link. This geometry corresponds to rfull005 tag SLCIO files.

We will download samples in directories dataNgev using reg. expressions. We will use 2 threads and 100 files (if available) from each sample.

hs-get pgun_pi_1to10000gev%rfull005 data1gev 2 100 pi1gev         # 1 GeV sample
hs-get pgun_pi_1to10000gev%rfull005 data10gev 2 100 pi10gev       # 10 GeV sample
hs-get pgun_pi_1to10000gev%rfull005 data100gev 2 100 pi100gev     # 100 GeV sample
hs-get pgun_pi_1to10000gev%rfull005 data1000gev 2 100 pi1000gev   # 1000 GeV sample
hs-get pgun_pi_1to10000gev%rfull005 data10000gev 2 100 pi10000gev # 10000 GeV sample

Run a Python script given here: pgun1.py. To process these samples, run this sequence. Each scripts creates a PDF file with the image.

fpad pgun1.py data1gev
fpad pgun1.py data10gev
fpad pgun1.py data100gev
fpad pgun1.py data1000gev
fpad pgun1.py data10000gev

A typical plot for 100 GeV single pion is shown here