Commit 7557c090 authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

Correction index freq pour 1420MHz, Reza 12/03/2015

parent a677380c
......@@ -91,9 +91,9 @@ int main(int narg, char* arg[])
<< " Imin,max,step="<<Imin<<","<<Imax<<","<<Istep<<" OutFile:"<<outfile<<"\n"
<< JFmin<<" <=NumFreq<= "<<JFmax<<" DeltaIAvg="<<deltaIavg<<" PrtLev="<<prtlev<<endl;
// Bande de frequence de 2 MHz autour du 21 cm (1420 MHz -> JF=2758)
sa_size_t JFmin21=2748,JFmax21=2768; // +/- 1 MHz autour de 1420 MHz
sa_size_t JFmin21_5=2676,JFmax21_5=2840; // +/- 5 MHz autour de 1420 MHz
// Bande de frequence de 2 MHz autour du 21 cm (1420 MHz -> JF=2792)
sa_size_t JFmin21=2776,JFmax21=2808; // +/- 1 MHz autour de 1420 MHz
sa_size_t JFmin21_5=2710,JFmax21_5=2874; // +/- 5 MHz autour de 1420 MHz
cout << " pjHI +/-1 MHz @1420 MHz " << JFmin21<<" <=NumFreq<= "<<JFmax21
<< " pjHI +/-5 MHz @1420 MHz " << JFmin21_5<<" <=NumFreq<= "<<JFmax21_5 << endl;
......@@ -137,17 +137,11 @@ int main(int narg, char* arg[])
cout << "visi2ntac/Info: dateobs="<<dateobs<<" SecondsPart()="<<dateobs.SecondsPart()<<endl;
for(int k=0; k<24; k++) xnt[2+k]=0.;
}
for(int k=0; k<8; k++) { // integration en bande large
for(int k=0; k<8; k++) {
TVector< complex<r_4> > vac = vismtx.Row(KVAC[k]);
xnt[2+k] += vac(freqs).Sum().real();
}
for(int k=0; k<8; k++) { // integration a +/- 1 MHz @ 1420 MHz
TVector< complex<r_4> > vac = vismtx.Row(KVAC[k]);
xnt[2+8+k] += vac(freqs21).Sum().real();
}
for(int k=0; k<8; k++) { // integration a +/- 5 MHz @ 1420 MHz
TVector< complex<r_4> > vac = vismtx.Row(KVAC[k]);
xnt[2+16+k] += vac(freqs21_5).Sum().real();
xnt[2+k] += vac(freqs).Sum().real(); // integration en bande large
xnt[2+8+k] += vac(freqs21).Sum().real(); // integration a +/- 1 MHz @ 1420 MHz
xnt[2+16+k] += vac(freqs21_5).Sum().real(); // integration a +/- 5 MHz @ 1420 MHz
}
I++;
if (I==deltaIavg) {
......
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