Commit 9d32df75 authored by Jean-Eric Campagne's avatar Jean-Eric Campagne
Browse files

Merge branch 'master' of gitlab.in2p3.fr:baoradio/AnaPAON4

parents faf18836 b1f4993a
......@@ -2,7 +2,7 @@
####### Exemples de commandes pour analyses des satellites ds PAON4
####### Analyse PAON4 ,
######## 1ere version , R. Ansari, Novembre 2018
### MAJ Decembre 2018 , Fev 2018
### MAJ Decembre 2018 , Fev 2019 , Juillet 2019
### Voir les commandes pour l'analyse la plus recente en fin de fichier
######################################################################
......@@ -202,6 +202,8 @@ ${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1342,1350 dt
${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1392,1400 dtsky1396_LD14jan19.ppf
## Extraction du signal pour les sources du ciel entre 1410-1418 MHz -> FreqCentral=1416 MHz
${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1410,1418 dtsky1416_LD14jan19.ppf
## Extraction du signal pour les sources du ciel entre 1435-1442 MHz -> FreqCentral=1439 MHz
${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1435,1442 dtsky1439_LD14jan19.ppf
#### Fabrication des DataTables track a partir des fichiers texte fourni par JEC
foreach f ( CasA_survey_2019-1-14 CasA_survey_2019-1-15 CasA_survey_2019-1-16 CygA_survey_2019-1-17 CygA_survey_2019-1-18 )
......@@ -326,3 +328,80 @@ piapp> ckacx $nom 1 3
piapp> ckacx $nom 1 3
### Utiliser le script pltsattrack.pic pour le trace des traces des satellites
#############################################################
############## Analyse de LD14Jan19 pour la polarisation V
cd /sps/baoradio/Reza/AnaP4/SatLD14jan19
set PRGBASE = /sps/baoradio/Reza/AnaP4/AnaPAON4
set PGP4 = $PRGBASE/Satellites/Objs/predictsatsgp4
set TRK2DT = $PRGBASE/Satellites/Objs/trk2dt
set TFM2DT = $PRGBASE/Objs/tfm2dt
set TRKFIT = $PRGBASE/Objs/trkacxfit
setenv DATAS /sps/baoradio/Olivier/PAON4/Work/LDOutputs/
## Extraction du signal pour les satellites entre 1275-1282 MHz (Galileo)
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1275,1282 dt_LD14jan19_V.ppf
## Extraction du signal pour les satellites entre 1270-1274 MHz (GPS) raie fine a 1272 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1270,1274 dt1272_LD14jan19_V.ppf
## Extraction du signal pour les satellites entre 1258-1266 MHz (GPS) en evitant la raie fine a 1272 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1258,1266 dt1262_LD14jan19_V.ppf
## Extraction du signal pour les satellites entre 1318-1328 MHz (GPS et Galileo faible intensite mais au milieu de la bande de frequence)
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1318,1328 dt1323_LD14jan19_V.ppf
## Extraction du signal pour les sources du ciel entre 1342-1350 MHz -> FreqCentral=1346 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1342,1350 dtsky1346_LD14jan19_V.ppf
## Extraction du signal pour les sources du ciel entre 1392-1400 MHz -> FreqCentral=1396 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1392,1400 dtsky1396_LD14jan19_V.ppf
## Extraction du signal pour les sources du ciel entre 1410-1418 MHz -> FreqCentral=1416 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1410,1418 dtsky1416_LD14jan19_V.ppf
## Extraction du signal pour les sources du ciel entre 1435-1442 MHz -> FreqCentral=1439 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1435,1442 dtsky1439_LD14jan19_V.ppf
#### Les datacards pour les traces vers cygA et casA
# On fabrique ldtrk_casA_V.d et ldtrk_cygA_V.d a partir de ldtrk_casA.d et ldtrk_cygA.d
#### ldtrk_casA_V.d (Donnees Polar V)
--- Observation Longue duree Janvier 2019
Sources avec observation vers CasA (z=+11.42)
@zenang 11.4
@trk dt_LD14jan19_V 1115,1160 1278.5 trk_40889_20190114
@trk dt_LD14jan19_V 1460,1500 1278.5 trk_43566_20190115
@trk dt_LD14jan19_V 2455,2495 1278.5 trk_43055_20190115
# Sky at 1346 MHz
@trk dtsky1346_LD14jan19_V 920,960 1346 CasA_survey_2019-1-14 NOAC
@trk dtsky1346_LD14jan19_V 2355,2395 1346 CasA_survey_2019-1-15 NOAC
@trk dtsky1346_LD14jan19_V 3790,3830 1346 CasA_survey_2019-1-16 NOAC
# Sky at 1396 MHz
@trk dtsky1396_LD14jan19_V 920,960 1396 CasA_survey_2019-1-14 NOAC
@trk dtsky1396_LD14jan19_V 2355,2395 1396 CasA_survey_2019-1-15 NOAC
@trk dtsky1396_LD14jan19_V 3790,3830 1396 CasA_survey_2019-1-16 NOAC
# Sky at 1416 MHz
@trk dtsky1416_LD14jan19_V 920,960 1416 CasA_survey_2019-1-14 NOAC
@trk dtsky1416_LD14jan19_V 2355,2395 1416 CasA_survey_2019-1-15 NOAC
@trk dtsky1416_LD14jan19_V 3790,3830 1416 CasA_survey_2019-1-16 NOAC
#### ldtrk_cygA_V.d (Donnees Polar V)
--- Observation Longue duree Janvier 2019
Sources avec observation vers CygA (z=-6.5)
@zenang -6.6
@trk dt_LD14jan19_V 4130,4175 1278.5 trk_43057_20190116
@trk dt_LD14jan19_V 6120,6160 1278.5 trk_41175_20190118
# Sky at 1346 MHz
@trk dtsky1346_LD14jan19_V 5025,5065 1346 CygA_survey_2019-1-17 NOAC
@trk dtsky1346_LD14jan19_V 6460,6500 1346 CygA_survey_2019-1-18 NOAC
# Sky at 1396 MHz
@trk dtsky1396_LD14jan19_V 5025,5065 1396 CygA_survey_2019-1-17 NOAC
@trk dtsky1396_LD14jan19_V 6460,6500 1396 CygA_survey_2019-1-18 NOAC
# Sky at 1416 MHz
@trk dtsky1416_LD14jan19_V 5025,5065 1416 CygA_survey_2019-1-17 NOAC
@trk dtsky1416_LD14jan19_V 6460,6500 1416 CygA_survey_2019-1-18 NOAC
### Ajustement pointage et phases
$TRKFIT -phifreq -D 4.5 -out casA_V.txt -ckf casA_V.ppf ldtrk_casA_V.d > resfitcasA_V.log
$TRKFIT -phifreq -D 4.5 -out cygA_V.txt -ckf cygA_V.ppf ldtrk_cygA_V.d > resfitcygA_V.log
$TRKFIT -phifreq -D 4.5 -out casA.txt -ckf casA.ppf ldtrk_casA.d > resfitcasA.log
$TRKFIT -phifreq -D 4.5 -out cygA.txt -ckf cygA.ppf ldtrk_cygA.d > resfitcygA.log
set f = cygA
$TRKFIT -D 4.5 -out $f.txt -ckf $f.ppf ldtrk_$f.d
###################################################################
#### Analyse dephasage ajustes a partir des satellites
#### AnaPAON4 - Reza , Juillet 2019
###################################################################
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
##### Le NTuple des phases AVANT correction des baselines
newnt nt z p2 p3 p4 p23 p24 p34 cph flg
ntfrascii nt ntphase.txt
print nt
set T (M_PI*z/180.)
# zone 1 2
# plot2d nt z p2 flg<1 'marker=fbox,7 nsta'
# plot2d nt z 130+360*(7*cos($T)+1.3*sin($T))/23 1 'marker=circle,9 same nsta red'
# plot2d nt z p34 flg<1 'marker=fbox,7 nsta'
# plot2d nt z 282+360*(2*cos($T)+1.4*sin($T))/23 1 'marker=circle,9 same nsta red'
# plot2d nt z p3+cph flg<1 'marker=fbox,7 nsta'
# Define the trigger delay phase correction offset
CPHOFF = 0.
newwin 2 1 1000 500
echo '======== Fitting Phase Phi12 ===> Position 2-1 '
exptovec Vz nt z flg<2 'nodisp'
exptovec Vp2 nt p2 flg<2 'nodisp'
del C
del Vph
cp Vp2 Vph
c++execfrf myfit.cc myfit.h
vec2var C viC
set vvC ( $viC )
plot2d nt z p2 flg<2 'marker=fbox,9 nsta notit'
settitle ' Phase Phi-12 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
# vecplot Vz Vp2 'marker=fbox,7 nocpts'
func $vvC[0]+($vvC[1]*cos(x*M_PI/180.)+$vvC[2]*sin(x*M_PI/180.))*360./23. -40 15 50 'line=solid,2 same red'
set lines ( "Phase_12=C(0)=$vvC[0] deg." "DeltaZ_12=C(1)=$vvC[1] cm" "DeltaY_12=C(2)=$vvC[2] cm" )
set cols ( red red red )
textdrawer lines cols 'font=helvetica,bold,16 black'
echo '======== Fitting Phase Phi13 ===> Position 3-1 '
exptovec Vz nt z flg<1 'nodisp'
exptovec Vp3 nt p3+cph+$CPHOFF flg<1 'nodisp'
del C
del Vph
cp Vp3 Vph
c++execfrf myfit.cc myfit.h
vec2var C viC
set vvC ( $viC )
plot2d nt z p3+cph+$CPHOFF flg<1 'marker=fbox,9 nsta notit'
settitle ' Phase Phi-13 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
func $vvC[0]+($vvC[1]*cos(x*M_PI/180.)+$vvC[2]*sin(x*M_PI/180.))*360./23. -40 15 50 'line=solid,2 same red'
set lines ( "Phase_13=C(0)=$vvC[0] deg." "DeltaZ_13=C(1)=$vvC[1] cm" "DeltaY_13=C(2)=$vvC[2] cm" )
set cols ( red red red )
textdrawer lines cols 'font=helvetica,bold,16 black'
w2ps
newwin 2 1 1000 500
echo '======== Fitting Phase Phi14 ===> Position 4-1 '
exptovec Vz nt z flg<1 'nodisp'
exptovec Vp4 nt p4+cph+$CPHOFF flg<1 'nodisp'
del C
del Vph
cp Vp4 Vph
c++execfrf myfit.cc myfit.h
vec2var C viC
set vvC ( $viC )
plot2d nt z p4+cph+$CPHOFF flg<1 'marker=fbox,9 nsta notit'
settitle ' Phase Phi-14 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
func $vvC[0]+($vvC[1]*cos(x*M_PI/180.)+$vvC[2]*sin(x*M_PI/180.))*360./23. -40 15 50 'line=solid,2 same red'
set lines ( "Phase=C(0)=$vvC[0] deg" "DeltaZ_14=C(1)=$vvC[1] cm" "DeltaY_14=C(2)=$vvC[2] cm" )
set cols ( red red red )
textdrawer lines cols 'font=helvetica,bold,16 black'
echo '======== Fitting Phase Phi34 ===> Position 4-3 '
exptovec Vz nt z flg<1 'nodisp'
exptovec Vp34 nt p34 flg<1 'nodisp'
del C
del Vph
cp Vp34 Vph
c++execfrf myfit.cc myfit.h
vec2var C viC
set vvC ( $viC )
plot2d nt z p34 flg<1 'marker=fbox,9 nsta notit'
settitle ' Phase Phi-34 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
func $vvC[0]+($vvC[1]*cos(x*M_PI/180.)+$vvC[2]*sin(x*M_PI/180.))*360./23. -40 15 50 'line=solid,2 same red'
set lines ( "Phase_34=C(0)=$vvC[0] deg." "DeltaZ_34=C(1)=$vvC[1] cm" "DeltaY_34=C(2)=$vvC[2] cm" )
set cols ( red red red )
textdrawer lines cols 'font=helvetica,bold,16 black'
w2ps
newwin 2 1 1000 500
echo '======== Fitting Phase Phi23 ===> Position 3-2 '
exptovec Vz nt z flg<1 'nodisp'
exptovec Vp23 nt p23+cph+$CPHOFF flg<1 'nodisp'
del C
del Vph
cp Vp23 Vph
c++execfrf myfit.cc myfit.h
vec2var C viC
set vvC ( $viC )
plot2d nt z p23+cph+$CPHOFF flg<1 'marker=fbox,9 nsta notit'
settitle ' Phase Phi-23 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
func $vvC[0]+($vvC[1]*cos(x*M_PI/180.)+$vvC[2]*sin(x*M_PI/180.))*360./23. -40 15 50 'line=solid,2 same red'
set lines ( "Phase_23=C(0)=$vvC[0] deg." "DeltaZ_23=C(1)=$vvC[1] cm" "DeltaY_23=C(2)=$vvC[2] cm" )
set cols ( red red red )
textdrawer lines cols 'font=helvetica,bold,16 black'
echo '======== Fitting Phase Phi24 ===> Position 4-2 '
exptovec Vz nt z flg<1 'nodisp'
exptovec Vp24 nt p24+cph+$CPHOFF flg<1 'nodisp'
del C
del Vph
cp Vp24 Vph
c++execfrf myfit.cc myfit.h
vec2var C viC
set vvC ( $viC )
plot2d nt z p24+cph+$CPHOFF flg<1 'marker=fbox,9 nsta notit'
settitle ' Phase Phi-24 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
func $vvC[0]+($vvC[1]*cos(x*M_PI/180.)+$vvC[2]*sin(x*M_PI/180.))*360./23. -40 15 50 'line=solid,2 same red'
set lines ( "Phase_24=C(0)=$vvC[0] deg." "DeltaZ_24=C(1)=$vvC[1] cm" "DeltaY_24=C(2)=$vvC[2] cm" )
set cols ( red red red )
textdrawer lines cols 'font=helvetica,bold,16 black'
w2ps
#################################################################
##### Le NTuple des phases apres correction des baselines
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
newnt ntc z p2 p3 p4 p23 p24 p34 cph flg
ntfrascii ntc ntphasecor.txt
print ntc
newwin
plot2d ntc z cph+$CPHOFF flg<3.5 'marker=box,11 navyblue nsta notit'
plot2d ntc z cph+$CPHOFF flg<0.5 'marker=fcircle,7 red same nsta notit'
settitle 'Trigger delay correction (3,4)-(1,2) 2ns->20.5 deg' ' ' 'font=helvetica,bold,14 black'
setaxelabels 'Zenith Angle (deg)' 'TriggerDelay (deg)' 'font=helvetica,bolditalic,16 black'
w2ps
newwin 2 1 1000 500
plot2d ntc z p2 flg<3.5 'marker=box,11 navyblue nsta notit'
plot2d ntc z p2 flg<0.5 'marker=fcircle,7 red same nsta notit'
settitle 'BaselineCor Phase Phi-12 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,14 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
plot2d ntc z p3+cph+$CPHOFF flg<3.5 'marker=box,11 navyblue nsta notit'
plot2d ntc z p3+cph+$CPHOFF flg<0.5 'marker=fcircle,7 red same nsta notit'
settitle 'BaselineCor Phase Phi-13 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,14 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
w2ps
newwin 2 1 1000 500
plot2d ntc z p4+cph+$CPHOFF flg<3.5 'marker=box,11 navyblue nsta notit'
plot2d ntc z p4+cph+$CPHOFF flg<0.5 'marker=fcircle,7 red same nsta notit'
settitle 'BaselineCor Phase Phi-14 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,14 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
plot2d ntc z p23+cph+$CPHOFF-360 flg<3.5 'marker=box,11 navyblue nsta notit'
plot2d ntc z p23+cph+$CPHOFF-360 flg<0.5 'marker=fcircle,7 red same nsta notit'
settitle 'BaselineCor Phase Phi-23 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,14 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
w2ps
newwin 2 1 1000 500
plot2d ntc z p24+cph+$CPHOFF-360 flg<3.5 'marker=box,11 navyblue nsta notit'
plot2d ntc z p24+cph+$CPHOFF-360 flg<0.5 'marker=fcircle,7 red same nsta notit'
settitle 'BaselineCor Phase Phi-24 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,14 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
plot2d ntc z ((p34<20)?p34:p34-360.) flg<3.5 'marker=box,11 navyblue nsta notit'
plot2d ntc z ((p34<20)?p34:p34-360.) flg<0.5 'marker=fcircle,7 red same nsta notit'
settitle 'BaselineCor Phase Phi-34 = f(ZenithAngle) ' ' ' 'font=helvetica,bold,14 black'
setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 black'
w2ps
######################################################################################
### Analyse ajustement phases en fonction de la frequence sur LD14jan19
######################################################################################
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
newnt phase pol src az freq p2 p3 p4 p23 p24 p34
ntfrascii phase phase_HV.txt
print phase
set CLNCUT (fabs(freq-1323)>1.)
set PHI0_H ( 168. 65. 50. )
set APHI_H ( -95. 382. 439. )
set PHI0_V ( 95. 335 390. )
set APHI_V ( 0. -260. -180. )
newwin 3 2 1200 800
for i 1:3
if ( $i == 1 ) then
set TIT 'Polar H , fcircle->CasA, box->CygA'
set PHI0 ( $PHI0_H )
set APHI ( $APHI_H )
else
set TIT 'Polar V , fcircle->CasA, box->CygA'
set PHI0 ( $PHI0_V )
set APHI ( $APHI_V )
endif
for j 0:3
jj = $j+2
set VP p$jj
if ( $j == 2 ) then
set VP ((p4<50.)?p4+360.:p4)
endif
set phi0 $PHI0[j]
set aphi $APHI[j]
plot2d phase freq $VP (pol==$i)&&(src==1)&&$CLNCUT 'marker=fcircle,9 navyblue notit nsta'
plot2d phase freq $VP (pol==$i)&&(src==2) 'marker=box,11 same navyblue notit nsta'
func $phi0+(x-1250.)/250.*$aphi 1250. 1500. 10 'line=solid,2 same red'
setaxelabels 'Freq (MHz)' " Phase Phi_1$jj (deg)" 'font=helvetica,bolditalic,16 black'
settitle "$TIT" ' ' 'font=helvetica,bolditalic,16 black'
end
end
###################################################################
#### Analyse directions pointage ajustees a partir des satellites
#### AnaPAON4 - Reza , Juillet 2019
###################################################################
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
newnt pnt dset z num D tet phi dtet dphi flg
ntfrascii pnt pointnt.txt
print pnt
set tlines ( "Dish #1" "Dish #2" "Dish #3" "Dish #4" )
set tcols ( black blue red orange )
set cols ( black black blue red orange )
newwin 2 2 1200 800
for i 1:5
plot2dw pnt z D dset (num==$i) "marker=fcircle,11 colbr128 nosizescale colscale nsta showcmap=right"
settitle " D_Dish[$i] = f(ZenithAngle) Color->DataSet" ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' "D_Dish[$i]" 'font=helvetica,bolditalic,16 black'
end
w2ps
set PHI (dphi*M_PI/180.)
newwin 2 2 1200 800
for i 1:5
plot2dw pnt z dtet*cos($PHI) dset (num==$i) "marker=ftriangle,11 colbr128 nosizescale colscale nsta showcmap=right"
settitle 'EW Pointing-Shift Color->DataSet' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'EW-Point-shift (deg)' 'font=helvetica,bolditalic,16 black'
end
w2ps
newwin 2 2 1200 800
for i 1:5
plot2dw pnt z dtet*sin($PHI) dset (num==$i) "marker=ftriangle,11 colbr128 nosizescale colscale nsta showcmap=right"
settitle 'NS Pointing-Shift Color->DataSet' ' ' 'font=helvetica,bold,16 black'
setaxelabels 'Zenith Angle (deg)' 'NS-Point-shift (deg)' 'font=helvetica,bolditalic,16 black'
end
w2ps
for i 1:5
exptovec DPEW$i pnt dtet*cos($PHI) (num==$i) 'nodisp'
exptovec DPNS$i pnt dtet*sin($PHI) (num==$i) 'nodisp'
end
exptovec VDS pnt dset (num==1) 'nodisp'
newnt rpnt DPEW1 DPEW2 DPEW3 DPEW4 DPNS1 DPNS2 DPNS3 DPNS4 dset
c++exec \
double xnt[15]; \
for(sa_size_t k=0; k<DPEW1.Size(); k++) { \
xnt[0]=DPEW1(k); xnt[1]=DPEW2(k); xnt[2]=DPEW3(k); xnt[3]=DPEW4(k); \
xnt[4]=DPNS1(k); xnt[5]=DPNS2(k); xnt[6]=DPNS3(k); xnt[7]=DPNS4(k); \
xnt[8]=VDS(k); rpnt.Fill(xnt); }
print rpnt
newwin 3 2 1200 800
for i 1:5
js = $i+1
for j $js:5
plot2dw rpnt DPEW$i DPEW$j dset 1 "marker=ftriangle,11 colbr128 nosizescale colscale nsta showcmap=right notit"
settitle "DPEW$j vs DPEW$i Color->DataSet" ' ' 'font=helvetica,bold,16 black'
end
end
w2ps
newwin 3 2 1200 800
for i 1:5
js = $i+1
for j $js:5
plot2dw rpnt DPNS$i DPNS$j dset 1 "marker=ftriangle,11 colbr128 nosizescale colscale nsta showcmap=right notit"
settitle "DPNS$j vs DPNS$i Color->DataSet" ' ' 'font=helvetica,bold,16 black'
end
end
w2ps
######################################################################################
### Analyse Pointage H-V sur LD14jan19
######################################################################################
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
newnt pointhv src az num D_H tet_H phi_H dtet_H dphi_H D_V tet_V phi_V dtet_V dphi_V
ntfrascii pointhv pointageHV_LD14jan19.txt
print pointhv
newwin 2 2 1200 800
for i 1:5
set PHI (dphi_H*M_PI/180.)
plot2dw pointhv dtet_H*cos($PHI) dtet_H*sin($PHI) az (num==$i) "marker=fcircle,15 colbr128 cmapdrw.font=helvetica,bold,14 nosizescale colscale notit nsta showcmap=right xylimits=-1.5,1.5,-1.5,1.5"
set PHI (dphi_V*M_PI/180.)
plot2dw pointhv dtet_V*cos($PHI) dtet_V*sin($PHI) az (num==$i) "marker=ftriangle,15 colbr128 nosizescale colscale nsta same"
settitle "Ant: $i Col->za H=circle V=triangle" ' ' 'font=helvetica,bold,16 black'
setaxelabels 'EW-Point-shift (deg)' 'NS-Point-Shift (deg) 'font=helvetica,bolditalic,16 black'
end
......@@ -141,10 +141,16 @@ public:
double Value(UnitVector const& uvo)
{
// circular beam response
double alp=acos(bdir_.Psc(uvo))*DoL_;
double alp=acos(bdir_.Psc(uvo))*DoL_;
// remplacer acos(...) par sin(acos(...)) ds l'expression ci-dessus pour avoir
if (fggauss_) {
// JEC , wiki BAORadio 2.1767 = 0.22054*Pi^2
return ( exp(-2.1767*alp*alp)*NormFac_ ); // Jiao : Check the factor 2 in the exp( ) -> exp(-2.*alp*alp)
// Reza / Sep 2019 : le facteur 2.61854 assure la meme largeur FWHM entre gaussian et Bessel
// On pose : a = D/lambda * gamma ( ou D/lambda * sin(gamma)
// Beam Bessel : beamJ(a) = ( 2 J1(pi a) / (pi a) )^2 Beam Gaussian : beamG(a) = Exp[ - K a^2 ]
// BeamJ(a=1.238098) = 0.5 et beamG(a=1.238098) = 0.5 pour K = 2.61854
// K=2.61854 au lieu de K=2.1767 (valeur de JEC , mais avec SphericalBesselJ[1...]
//---- return ( exp(-2.1767*alp*alp)*NormFac_ ); // Jiao : Check the factor 2 in the exp( ) -> exp(-2.*alp*alp)
return ( exp(-2.61854*alp*alp)*NormFac_ ); // Jiao : Check the factor 2 in the exp( ) -> exp(-2.*alp*alp)
}
else {
if (alp<1.e-19) return NormFac_;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment