### Script piapp pour faire tourner le programme de simulation 
### de reconstruction de beam par la technique de double FFT
###  Reza , Sep 2006
#  Le code (2 classes) se trouve dans mbeamcyl.h mbeamcyl.cc
#  et une classe auxiliaire brsrc.h brsrc.cc 
#  Classe  MultiBeamCyl et la classe utilitaire BRSourceGen
# --------------------
# 0/ Verification, modif des parametres
#  Modifier si besoin le code dans trecbeam.cc qui instancie
#  les classes (MultiBeamCyl en particulier) avec le parametres
#  necessaires (nb de recepteurs, nb d'echantillons ...) 
#  Avant de faire la simulation et la reconstruction de beam

# 1/ Faire tourner le programme - reconstruction plan X-freq (en dehors de piapp)
Objs/treccyl X onecyl.ppf


# 2/ Affichage resultat
#  Afficher le plan source reconstruit (les frequences le 
# long de l'axe Y (les lignes de la matrice), et les directions
# des sources (angles) suivant l'axe Ox (colonnes)
# Le zero des angles se trouve au milieu
delobjs *
openppf onecyl.ppf
disp recsrcplane 'lut=sqrt,0.,0.1 invy grey32 win'

# 2.b/ le ntuple ntsrc contient la liste des sources de depart
print ntsrc
# Attention, ce qui est reconstruit n'est pas directement l'angle, mais
# sin(angle)*frequence - pour une matrice 512lignesx256colonnes
plot2dw ntsrc sin(angrad)*freq*0.25*256+128 freqred*512*2 amp 1 'nsta same red marker=circle,7'
addtitle 'PlanSource (Y=freq,X=AngNS)' ' ' 'font=times,bold,18 red'
w2eps recbeam.eps


# 3/ Faire tourner le programme - reconstruction cube X-Y-freq (en dehors de piapp)
Objs/treccyl XY fourcyl.ppf

# 4/ Affichage resultats
delobjs *
openppf fourcyl.ppf
listobjs

disp recsrcplane0 'lut=sqrt,0.,0.1 invy grey32 win'
plot2dw ntsrc sin(angrad)*freq*0.25*256+128 freqred*512*2 amp 1 'nsta same red marker=circle,7'
disp recsrcplane3 'lut=sqrt,0.,0.1 invy grey32 win'
plot2dw ntsrc sin(angrad)*freq*0.25*256+128 freqred*512*2 amp 1 'nsta same red marker=circle,7'