// GLISSANDO ver. 1.6
// typical running sequences

// installation
sh install

// sample run
./glissando.exe

// generate nuclear profile of nucleus A, no expulsion
./glissando_profile.exe input_profile_0.dat             
root 
.x fitr.C("")

// generate nuclear profile of nucleus A, expulsion distance 0.4 fm, fix-last algorithm
./glissando_profile.exe input_profile_04.dat             
root 
.x fitr.C("")

// generate nuclear profile of nucleus A, expulsion distance 0.4 fm, fix-last algorithm
./glissando_profile.exe input_profile_08.dat             
root 
.x fitr.C("")

// generate nuclear profile of nucleus A, expulsion distance 0.4 fm, return-to-beginning algorithm (SLOW!)
./glissando_profile.exe input_profile_04_ret.dat         
root 
.x fitr.C("")

./glissando.exe input_snap.dat
root
.x density.C("")

// generate centrality classes for min-bias Au+Au collisions
./glissando.exe input_minbias.dat           
root
.x centrality.C("")
.x epsilon.C("")
.x epsilon_b.C("")
.x dxdy.C("")
...

// generate most central Au+Au collisions
./glissando.exe input_0_5.dat glissando_0_5.root        
// check b distribution
root 
TBrowser a, open "glissando_0_5.root", go to "events" tree, click branch "b"
// generate 2D density plots and Fourier profiles
root 
.x density.C("glissando_0_5.root")
.x profile2.C("glissando_0_5.root")
...

// run the interpolation code
interpolation.exe glissando_0_5.root

// repeat the above and for other centrality classes  
./glissando.exe input_30_40.dat glissando_30_40.root    
...
./glissando.exe input_80_95.dat glissando_80_95.root   
...

// analysis of the hot-spot model
./glissando.exe input_hotspot.dat hs.root
root 
.x centrality.C("hs.root")
.x epsilon.C("hs.root")
.x epsilon_b.C("hs.root")
...

// analysis for asymmetric collisions
./glissando.exe input_S_Pb.dat
root 
.x centrality.C("")
.x epsilon.C("")
...

// p+Au collisions
./glissando_prot.exe
root 
.x centrality.C("")
...

// d+Au collisions
./glissando_deut.exe
root 
.x centrality.C("")
...

// LHC predictions
./glissando.exe input_minbias_LHC.dat
root
.x centrality.C("")
.x epsilon.C("")
...

// retrieving the full event info (LARGE OUTPUT DATA!)
// (FULL=1 flag must be set in GLISSANDO input file)
/.glissando.exe input_full.dat full.root
./retrieve.exe full.root
...

// Gaussian wounding profile (SLOW!)
./glissando_gauss.exe
...
