Commit a677380c authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

modif calcul puiss@1420 pour avoir la ligne de base de P, adaptation plac8.pic, Reza 12/03/2015

parent 53856325
......@@ -50,30 +50,15 @@ end
###### Trace de la puissance a +/-1 MHz @ 1420 MHz
newwin 1 1 800 600
plot2d $NT time/3600 p1HI*$FNORM[0] 1 " nsta cpts notit $COLS[0] "
plot2d $NT time/3600 (p1HI-p1HI5)*$FNORM[0] 1 " nsta cpts notit $COLS[0] "
### trace de p2 ... p8
for i 1:8
j = i+1
plot2d $NT time/3600 p${j}HI*$FNORM[i] 1 "nsta cpts same notit $COLS[i] "
plot2d $NT time/3600 (p${j}HI-p${j}HI5)*$FNORM[i] 1 "nsta cpts same notit $COLS[i] "
end
### Axes et titres
setaxelabels 'Time (hours)' 'P1...P8 @1420+/-1 MHz ' 'font=helvetica,bolditalic,16 black'
settitle "$titre 1420+/-1 MHz" ' ' 'font=helvetica,bold,18 black'
###### Trace de la puissance a +/-5 MHz @ 1420 MHz
newwin 1 1 800 600
plot2d $NT time/3600 p1HI5*$FNORM[0] 1 " nsta cpts notit $COLS[0] "
### trace de p2 ... p8
for i 1:8
j = i+1
plot2d $NT time/3600 p${j}HI5*$FNORM[i] 1 "nsta cpts same notit $COLS[i] "
end
### Axes et titres
setaxelabels 'Time (hours)' 'P1...P8 @1420+/-5 MHz ' 'font=helvetica,bolditalic,16 black'
settitle "$titre 1420+/-5 MHz" ' ' 'font=helvetica,bold,18 black'
setaxelabels 'Time (hours)' 'P1...P8 excess @1420' 'font=helvetica,bolditalic,16 black'
settitle "$titre 1420+/-1 MHz - baseline (pjHI-pjHI5)" ' ' 'font=helvetica,bold,18 black'
###########################################################
......@@ -151,12 +151,15 @@ int main(int narg, char* arg[])
}
I++;
if (I==deltaIavg) {
// puissance moyenne dans la bande large
double fnorm=(1./(double)deltaIavg)/(double)(JFmax-JFmin+1);
for(int k=0; k<8; k++) xnt[2+k]*=fnorm;
// puissance moyenne ds +/5 MHz @1420 , exclu +/- 1 MHz centre a 1420
fnorm=(1./(double)deltaIavg)/(double)((JFmax21_5-JFmin21_5)-(JFmax21-JFmin21));
for(int k=0; k<8; k++) xnt[2+16+k] = (xnt[2+16+k]-xnt[2+8+k])*fnorm;
// puissance moyenne ds +/1 MHz @1420
fnorm=(1./(double)deltaIavg)/(double)(JFmax21-JFmin21+1);
for(int k=0; k<8; k++) xnt[2+8+k]*=fnorm;
fnorm=(1./(double)deltaIavg)/(double)(JFmax21_5-JFmin21_5+1);
for(int k=0; k<8; k++) xnt[2+16+k]*=fnorm;
nt.Fill(xnt);
I=0; cntnt++;
}
......
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