#!/bin/bash

# init
function pause(){
   read -p "$*"
}

# running

./glissando -h

pause "Press ENTER to run the fastest fun example: a single event in the hot-spot model"
./glissando input_snap.dat
root -b -l -q -x "density.C(\"glissando.root\")"
gv 3D-density.eps &
gv 2D-density.eps &

pause "Press ENTER to run proton-Pb collisions (1 min.)"
./glissando input_p_Pb.dat pPb.root
root -b -l -q -x "centrality2.C(\"pPb.root\")" 
gv centrality_b.eps &

pause "Press ENTER to run deuteron-Pb collisions (1 min.)"
./glissando input_d_Pb.dat dPb.root
root -b -l -q -x "density.C(\"dPb.root\")" 
gv 3D-density.eps &
gv 2D-density.eps &

pause "Press ENTER to run a minimum-bias Pb-Pb run for LHC (8 min.)"
./glissando input_minbias_LHC.dat lhc.root 

root -b -l -q -x "centrality2.C(\"lhc.root\")"
gv centrality_b.eps &
gv centrality_nw.eps &
gv centrality_RDS.eps &

root -b -l -q -x "epsilon.C(\"lhc.root\")" 
gv epsilon.eps &
gv sigma_epsilon.eps &

root -b -l -q -x "mult.C(\"lhc.root\")" 
gv om_RDS.eps &

root -b -l -q -x "fourier.C(\"lhc.root\")" 
gv epsn.eps &

root -b -l -q -x "size.C(\"lhc.root\")" 
gv sigma_r.eps &

pause "Press ENTER to run the centrality 30-40% run for RHIC at 200 GeV (11 min.)"
./glissando input_30_40.dat glissando_30_40.root        
root -b -l -q -x "density.C(\"glissando_30_40.root\")" 
gv 3D-density.eps &
gv 2D-density.eps &
root -b -l -q -x "core_mantle.C(\"glissando_30_40.root\")" 
gv 3D-cm-density.eps &
gv 2D-cm-density.eps &
root -b -l -q -x "profile2.C(\"glissando_30_40.root\")" 
gv f0_fig.eps &
gv f2_fig.eps &
gv f4_fig.eps &
root -b -l -q -x "dxdy.C(\"glissando_30_40.root\")" 
gv dxdy.eps &


pause "Press ENTER to see other functions (4 min.)"
pause "(make sure the prepocessor parameters in Makefile_aux were at compilation as originally set)"
./glissando_aux input_aux.dat glissando_aux.root        
root -b -l -q -x "fitr.C(\"glissando_aux.root\")" 
gv fitr.eps &
root -b -l -q -x "corr.C(\"glissando_aux.root\")" 
gv corr.eps &
root -b -l -q -x "overlay.C(\"glissando_aux.root\")" 
gv overlay_pr.eps &
root -b -l -q -x "tilted.C(\"glissando_aux.root\")" 
gv tilted_pr.eps &

pause "Press ENTER for a sample run with nuclear correlations (1 min.)"
pause "(first you must get and save some files from http://www.phys.psu.edu/∼malvioli/eventgenerator/ and create one big file, for instance running cat o16-1.dat o16-2.dat o16-3.dat [more files] > o16.dat )"

./glissando_files input_O16.dat o16.root o16.dat o16.dat
root -b -l -q -x "corr.C(\"o16.root\")" 
gv corr.eps &
root -b -l -q -x "wounding_profile.C(\"o16.root\")" 
gv wounding_pr.eps &

pause "Press ENTER to use the Gaussian wounding profile  (4 min.)"
./glissando_gauss input_aux.dat glissando_gauss.root        
root -b -l -q -x "wounding_profile.C(\"glissando_gauss.root\")" 
gv wounding_pr.eps &
root -b -l -q -x "epsilon.C(\"glissando_gauss.root\")" 
gv epsilon.eps &
gv sigma_epsilon.eps &

pause "Press ENTER to run an asymmetric collision (1 min.)"
./glissando input_S_Pb.dat SPb.root
root -b -l -q -x "epsilon.C(\"SPb.root\")" 
gv epsilon.eps &
root -b -l -q -x "mult.C(\"SPb.root\")" 
gv om_RDS.eps &

pause "Press ENTER to run the intepolation code"
./interpolation SPb.root

pause "Press ENTER to generate and retrieve the full event info (LARGE OUTPUT DATA!)"
pause "(FULL=1 flag must be set in the GLISSANDO input file)"
./glissando input_full.dat full.root
./retrieve full.root


