Changeset 3165 in Sophya for trunk/Cosmo/RadioBeam/recbeam.pic

Timestamp:

Feb 1, 2007, 3:15:16 PM (19 years ago)

Author:

ansari

Message:

petites ameliorations + reglage parametres -> version +/- utilisable , Reza 01/02/2007

File:

• trunk/Cosmo/RadioBeam/recbeam.pic (modified) (3 diffs)

trunk/Cosmo/RadioBeam/recbeam.pic

@@ -1 @@
-### Script piapp pour faire tourner le programme de simulation 
+### Script piapp pour afficher les resultats du programme 
 ### 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
-
+#  Modifier si besoin le code dans treccyl.cc et faire ensuite make
+csh> make 
+# Le programme se trouve dans Objs/treccyl 
 # 1/ Faire tourner le programme - reconstruction plan X-freq (en dehors de piapp)
 Objs/treccyl X onecyl.ppf
@@ -35 +32 @@
 
 # 3/ Faire tourner le programme - reconstruction cube X-Y-freq (en dehors de piapp)
-Objs/treccyl XY ncyl.ppf
+Objs/treccyl XY12 n12cyl.ppf
+OU 
+Objs/treccyl XY5 n5cyl.ppf
 
 # 4/ Affichage resultats
 delobjs *
-openppf ncyl.ppf
+openppf n12cyl.ppf
+OU 
+openppf n5cyl.ppf
+
 listobjs
 
+## Pour afficher les plans reconstruits (angX-freq) sur les cylindres 0 , 3 avec les sources
+## ATTENTION, si on change les parametres, la conversion angle->num pixel doit etre adapte
 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'
@@ -47 +51 @@
 plot2dw ntsrc sin(angrad)*freq*0.25*256+128 freqred*512*2 amp 1 'nsta same red marker=circle,7'
 
+## Pour afficher trois plans (angX-angY) avec les sources, autour des trois frequences des sources 
 zone 1 3
-disp recXYf0 'zoomx2 grey32 lut=sqrt,0.,1'
-plot2dw ntsrc sin(angrad)*freq*0.25*256+128 sin(angyrad)*freq*0.25*21+10  amp fabs(freqred-0.1)<0.01 'nsta same red marker=circle,7'
-disp recXYf1 'zoomx2 grey32 lut=sqrt,0.,1'
-disp recXYf2 'zoomx2 grey32 lut=sqrt,0.,1'
+disp recXYf0 'zoomx2 invy grey32 lut=sqrt,0.,1.5'
+plot2dw ntsrc sin(angrad)*freq*0.25*256+128 600.*angyrad/M_PI+10  amp fabs(freqred-0.1)<0.01 'nsta same red marker=circle,7'
+disp recXYf1 'zoomx2 invy grey32 lut=sqrt,0.,1.5'
+plot2dw ntsrc sin(angrad)*freq*0.25*256+128 600.*angyrad/M_PI+10  amp fabs(freqred-0.27)<0.01 'nsta same red marker=circle,7'
+disp recXYf2 'zoomx2 invy grey32 lut=sqrt,0.,1.5'
+plot2dw ntsrc sin(angrad)*freq*0.25*256+128 600.*angyrad/M_PI+10  amp fabs(freqred-0.38)<0.01 'nsta same red marker=circle,7'
 
-
-
+## On peut aussi extraire une coupe angX-angY de la boite 3D reconstruit
 objaoper recsrcbox slicexy 101
 set cobj /home/slicexy_101