Skip to content
GitLab
Projects
Groups
Snippets
Help
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Open sidebar
BAORadio
AnaPAON4
Commits
eb539467
Commit
eb539467
authored
Feb 27, 2019
by
perdereau
Browse files
Merge branch 'master' of
https://gitlab.in2p3.fr/baoradio/AnaPAON4
parents
47e0916b
9aa276e6
Changes
24
Expand all
Hide whitespace changes
Inline
Side-by-side
Showing
20 changed files
with
1956 additions
and
165 deletions
+1956
-165
AnaFringes/fringeV6.cc
AnaFringes/fringeV6.cc
+138
-69
AnaFringes/makefile
AnaFringes/makefile
+2
-2
Scripts/ckfitres.pic
Scripts/ckfitres.pic
+79
-31
Scripts/p4sat.cmd
Scripts/p4sat.cmd
+135
-1
TrackingSrc/altaz/.idea/altaz.iml
TrackingSrc/altaz/.idea/altaz.iml
+11
-0
TrackingSrc/altaz/.idea/encodings.xml
TrackingSrc/altaz/.idea/encodings.xml
+4
-0
TrackingSrc/altaz/.idea/misc.xml
TrackingSrc/altaz/.idea/misc.xml
+4
-0
TrackingSrc/altaz/.idea/modules.xml
TrackingSrc/altaz/.idea/modules.xml
+8
-0
TrackingSrc/altaz/.idea/workspace.xml
TrackingSrc/altaz/.idea/workspace.xml
+178
-0
TrackingSrc/altaz/CasA_survey_Exemple.txt
TrackingSrc/altaz/CasA_survey_Exemple.txt
+245
-0
TrackingSrc/altaz/altaz.py
TrackingSrc/altaz/altaz.py
+238
-0
TrackingSrc/altaz/tracksrc.sh
TrackingSrc/altaz/tracksrc.sh
+20
-0
acbeam.h
acbeam.h
+22
-0
cxbeam.h
cxbeam.h
+41
-5
gacfit.h
gacfit.h
+105
-18
gcxfit.h
gcxfit.h
+98
-26
gcxfitbaseline.h
gcxfitbaseline.h
+354
-0
makefile
makefile
+39
-11
p4autils.cc
p4autils.cc
+2
-2
tkfit_minuit.h
tkfit_minuit.h
+233
-0
No files found.
AnaFringes/fringeV6.cc
View file @
eb539467
...
...
@@ -72,7 +72,7 @@ int nloop = 0;
#define POSXYZ_FIXED 1
//JEC 20/7/18 auto-corr azimut contraint >0 (en general non => 0)
#define AZIMUT_FIXED
1
#define AZIMUT_FIXED
0
...
...
@@ -295,6 +295,7 @@ struct PARAM {
r_8
decSrcSec
;
TimeStamp
MeridianDate
;
//passage au meridien de la source
TimeStamp
ReferenceDate
;
//JEC 18/2/19 reference des temps pour TimeVec
Now
LocalGeoTime
;
int
year0
;
...
...
@@ -1068,6 +1069,9 @@ void fitAutoCrossCorr(){
TimeStamp
meridianT
=
param
.
MeridianDate
;
//JEC 18/2/19 new input the reference date HH::MM::SS
TimeStamp
referenceDate
=
param
.
ReferenceDate
;
r_4
hourAngleMax
=
param
.
hourAngleMax
;
r_4
hourAngleMin
=
param
.
hourAngleMin
;
...
...
@@ -1107,14 +1111,14 @@ void fitAutoCrossCorr(){
// ATTENTION !!!!!!!!! Repere Oriente vers North positif X, vers West positif Y, Z vertical
vector
<
r_8
>
posXCh
(
Ndish
),
posYCh
(
Ndish
),
posZCh
(
Ndish
);
//Geometrical initial values
//Geometrical initial values
(units = meters)
posXCh
[
0
]
=
0.
;
posYCh
[
0
]
=
0.
;
posZCh
[
0
]
=
0.
;
posXCh
[
1
]
=
0.
;
posYCh
[
1
]
=
5.993
;
posZCh
[
1
]
=
0.
;
posZCh
[
1
]
=
0.
;
//was 0.100 pour tester => phase se translate
posXCh
[
2
]
=
-
5.996
;
posYCh
[
2
]
=
-
4.380
;
...
...
@@ -1272,18 +1276,53 @@ void fitAutoCrossCorr(){
//RAVec
//JEC 18/2/19 : attention we have during long data acquisition many [0h,24h] a cyclic RAVec
//2: select indices according to TimeVec: #seconds since a 0h0 of a given day
TArray
<
r_8
>
RAVec
;
fIn
>>
PPFNameTag
(
"RAVec"
)
>>
RAVec
;
TArray
<
r_8
>
TimeVec
;
fIn
>>
PPFNameTag
(
"TimeVec"
)
>>
TimeVec
;
//number of second since reference date of source meridian transit
r_8
tOfMeridianSinceRef
=
TimeStamp
::
TimeDifferenceSeconds
(
meridianT
,
referenceDate
);
//select +/- 5h arround
r_8
rangeOfTime
=
5.0
*
3600.
;
//sec
sa_size_t
ifirstTime
=
0
;
sa_size_t
ilastTime
=
0
;
{
bool
first
=
true
;
bool
last
=
true
;
for
(
sa_size_t
i
=
0
;
i
<
TimeVec
.
Size
();
i
++
){
if
(
fabs
(
TimeVec
(
i
)
-
tOfMeridianSinceRef
)
<=
rangeOfTime
){
if
(
first
){
first
=
false
;
ifirstTime
=
i
;
last
=
true
;
}
if
(
!
first
&&
last
)
ilastTime
=
i
;
}
else
{
last
=
false
;
}
}
}
//find ifirstTime/ilastTime
if
(
ifirstTime
==
ilastTime
)
throw
RangeCheckError
(
"Time selection not valid"
);
cout
<<
"Select Time/RA range in index ["
<<
ifirstTime
<<
", "
<<
ilastTime
<<
"]"
<<
endl
;
TVector
<
r_4
>
HAVec
(
RAVec
.
Size
());
//JEC 9/6/18 seems that xephem routine has bugs
for
(
sa_size_t
i
=
0
;
i
<
RAVec
.
Size
();
i
++
){
HAVec
(
i
)
=
hrrad
(
RAVec
(
i
)
-
raSrc
);
//hour angle = RA-raSrc hour -> radian
}
//Selection des Angles Horaires
r_8
haMin
=
Angle
(
hourAngleMin
,
Angle
::
Degree
).
ToRadian
();
r_8
haMAX
=
Angle
(
hourAngleMax
,
Angle
::
Degree
).
ToRadian
();
sa_size_t
ifirstHA
=
0
;
...
...
@@ -1292,7 +1331,10 @@ void fitAutoCrossCorr(){
//Select the first transit of the source during observation
bool
first
=
true
;
bool
last
=
true
;
for
(
sa_size_t
i
=
0
;
i
<
HAVec
.
Size
();
i
++
){
//JEC 18/2/19 restriction of the good HA range according to the ifirstTime/ilastTime range
//was for(sa_size_t i=0; i<HAVec.Size();i++){
for
(
sa_size_t
i
=
ifirstTime
;
i
<
ilastTime
;
i
++
){
if
(
HAVec
(
i
)
<
haMAX
&&
HAVec
(
i
)
>
haMin
){
if
(
first
){
first
=
false
;
...
...
@@ -2889,25 +2931,38 @@ int main(int narg, char* arg[]) {
string
source
=
"CasA31oct17"
;
string
meridianDate
=
"2017-10-31T20:00:00"
;
//passage au meridien format YYYY-MM-DDTHH:MM:SS
//JEC 18/2/19 reference date
string
referenceDate
=
"2019-01-14T00:00:00"
;
//YYYY-MM-DDTHH:MM:SS
string
inputFile
=
""
;
string
inputParamFile
=
""
;
string
outputFile
=
""
;
//CasA RA-DEC
r_8
raSrcHour
=
23.
;
//RA H:M:S.S
r_8
raSrcMin
=
23.
;
r_8
raSrcSec
=
27.9
;
r_8
decSrcDeg
=
58
;
//DEC D:M:S.S
r_8
decSrcMin
=
48
;
r_8
decSrcSec
=
42
;
r_4
raSrcHour
=
23.
;
//RA H:M:S.S
r_4
raSrcMin
=
23.
;
r_4
raSrcSec
=
27.9
;
r_4
decSrcDeg
=
58
;
//DEC D:M:S.S
r_4
decSrcMin
=
48
;
r_4
decSrcSec
=
42
;
//Crabe RA-DEC
// r_4 raSrcHour = 5.; //RA H:M:S.S
// r_4 raSrcMin = 34.;
// r_4 raSrcSec = 31.94;
// r_4 decSrcDeg = 22; //DEC D:M:S.S
// r_4 decSrcMin = 0;
// r_4 decSrcSec = 52.2;
//CygA RA-DEC
// r_
8
raSrcHour = 19.; //RA H:M:S.S
// r_
8
raSrcMin = 59.;
// r_
8
raSrcSec = 28.3;
// r_
8
decSrcDeg = 40.; //DEC D:M:S.S
// r_
8
decSrcMin = 44.;
// r_
8
decSrcSec = 2.0;
// r_
4
raSrcHour = 19.; //RA H:M:S.S
// r_
4
raSrcMin = 59.;
// r_
4
raSrcSec = 28.3;
// r_
4
decSrcDeg = 40.; //DEC D:M:S.S
// r_
4
decSrcMin = 44.;
// r_
4
decSrcSec = 2.0;
r_4
FreqMin
=
1250.0
;
//MHz //Freq minimal, maxi, step OR Freq. centrale + bandwidth
r_4
FreqMax
=
1500.0
;
//MHz
...
...
@@ -2929,7 +2984,9 @@ int main(int narg, char* arg[]) {
int
ka
=
1
;
while
(
ka
<
narg
)
{
cout
<<
"option "
<<
arg
[
ka
]
<<
" process...."
<<
endl
;
if
(
strcmp
(
arg
[
ka
],
"-h"
)
==
0
)
{
cout
<<
"usage ex: ./bin/fringeV2 -hasel -10,10 >& tmp
\n
"
;
cout
<<
"options:
\n
"
;
...
...
@@ -2944,6 +3001,7 @@ int main(int narg, char* arg[]) {
cout
<<
"-ra: Right Ascension (J2000 H:M:S.S) (ex. 23:23:27.9)
\n
"
;
cout
<<
"-dec: Declinaison (J2000 D:M:S.S) (ex. 58:48:42.0)
\n
"
;
cout
<<
"-meridian: approximate meridian date (ex. 2017-10-18T21:00:00)
\n
"
;
cout
<<
"-refdate: date de reference (ex. 2019-01-14T00:00:00)
\n
"
;
cout
<<
"-f: freq. range MHz freqMin,freqMax,freqBand (ex. 1300,1450,1): notice freqBand not used
\n
"
;
cout
<<
"-hasel: hour angle selected in [hamin,hamax] (Degree): (ex. -10,10)
\n
"
;
cout
<<
"-type: type de cross REal or IMag [RE] (WARNING PARTIE_REELLE define sould be tuned accordingly for the moment)
\n
"
;
...
...
@@ -2986,11 +3044,14 @@ int main(int narg, char* arg[]) {
ka
+=
2
;
}
else
if
(
strcmp
(
arg
[
ka
],
"-ra"
)
==
0
)
{
scanf
(
arg
[
ka
+
1
],
"%f:%f:%f"
,
&
raSrcHour
,
&
raSrcMin
,
&
raSrcSec
);
cout
<<
"option -ra <"
<<
arg
[
ka
+
1
]
<<
">"
<<
endl
;
sscanf
(
arg
[
ka
+
1
],
"%f:%f:%f"
,
&
raSrcHour
,
&
raSrcMin
,
&
raSrcSec
);
cout
<<
raSrcHour
<<
","
<<
raSrcMin
<<
","
<<
raSrcSec
<<
endl
;
ka
+=
2
;
}
else
if
(
strcmp
(
arg
[
ka
],
"-dec"
)
==
0
)
{
scanf
(
arg
[
ka
+
1
],
"%f:%f:%f"
,
&
decSrcDeg
,
&
decSrcMin
,
&
decSrcSec
);
cout
<<
"option -dec <"
<<
arg
[
ka
+
1
]
<<
">"
<<
endl
;
sscanf
(
arg
[
ka
+
1
],
"%f:%f:%f"
,
&
decSrcDeg
,
&
decSrcMin
,
&
decSrcSec
);
ka
+=
2
;
}
else
if
(
strcmp
(
arg
[
ka
],
"-f"
)
==
0
){
...
...
@@ -3002,6 +3063,10 @@ int main(int narg, char* arg[]) {
meridianDate
=
arg
[
ka
+
1
];
ka
+=
2
;
}
else
if
(
strcmp
(
arg
[
ka
],
"-refdate"
)
==
0
){
referenceDate
=
arg
[
ka
+
1
];
ka
+=
2
;
}
else
if
(
strcmp
(
arg
[
ka
],
"-hasel"
)
==
0
){
sscanf
(
arg
[
ka
+
1
],
"%f,%f"
,
&
hourAngleMin
,
&
hourAngleMax
);
ka
+=
2
;
...
...
@@ -3046,7 +3111,7 @@ int main(int narg, char* arg[]) {
}
//Computation of the exact median passage
r_8
raSrc
=
HdecfrHMS
(
raSrcHour
,
raSrcMin
,
raSrcSec
);
//J2000
r_8
raSrc
=
HdecfrHMS
(
(
int
)
raSrcHour
,
(
int
)
raSrcMin
,
(
r_8
)
raSrcSec
);
//J2000
Now
lGeoTime
;
//Old value
...
...
@@ -3093,76 +3158,80 @@ int main(int narg, char* arg[]) {
//calcul du passage au meridien
r_8
day0
=
iday0
+
HdecfrHMS
(
0
,
0
,
0
)
/
24.
;
//noon
r_8
day0
=
iday0
+
HdecfrHMS
(
hr
0
,
0
,
0
)
/
24.
;
//noon
//
{ //OLD
//
double dt = 0;
//
double lst;
//
#define MAXLOOPS 10
//
#define MAXERR (1./3600.) /* hours */
//
int i =0;
//
do {
//
day0 += dt/24.;
//
lGeoTime.n_mjd = MJDfrDate(day0,month0,year0);
//
// calcul passage meridien
//
now_lst(&lGeoTime,&lst);
//
dt = raSrc-lst;
//
cout << "METHOD OLD: meridien dt = " << dt << " " << ToStringHdec(dt) << endl;
//
} while (++i < MAXLOOPS && fabs(dt) > MAXERR);
//
if (i == MAXLOOPS) throw NotFoundExc("Meridian date not found");
{
//OLD
double
dt
=
0
;
double
lst
;
#define MAXLOOPS 10
#define MAXERR (1./3600.)
/* hours */
int
i
=
0
;
do
{
day0
+=
dt
/
24.
;
lGeoTime
.
n_mjd
=
MJDfrDate
(
day0
,
month0
,
year0
);
// calcul passage meridien
now_lst
(
&
lGeoTime
,
&
lst
);
dt
=
raSrc
-
lst
;
cout
<<
"METHOD OLD: meridien dt = "
<<
dt
<<
" "
<<
ToStringHdec
(
dt
)
<<
endl
;
}
while
(
++
i
<
MAXLOOPS
&&
fabs
(
dt
)
>
MAXERR
);
if
(
i
==
MAXLOOPS
)
throw
NotFoundExc
(
"Meridian date not found"
);
//
#undef MAXLOOPS
//
#undef MAXERR
#undef MAXLOOPS
#undef MAXERR
//
meridianT.SetHour(ToStringHdec((day0-iday0)*24.));
meridianT
.
SetHour
(
ToStringHdec
((
day0
-
iday0
)
*
24.
));
//
}//OLD
}
//OLD
// cout << "OLD passage ay meridien " << meridianT << endl;
{
//NEW
//
{//NEW
double
dt
=
0
;
#define MAXLOOPS 10
#define MAXERR (1./3600.)
/* hours */
int
i
=
0
;
//
double dt = 0;
//
#define MAXLOOPS 10
//
#define MAXERR (1./3600.) /* hours */
//
int i =0;
r_8
hr
=
0
;
//
r_8 hr = 0;
JMA
day
;
day
.
Annee
=
year0
;
day
.
Mois
=
month0
;
day
.
Jour
=
iday0
;
do
{
hr
+=
dt
;
HMS
hour
=
HtoHMS
(
hr
);
HMS
tsid
;
tsid
=
TUtoTSid
(
day
,
hour
);
r_8
ra
=
HMStoH
(
tsid
);
// le ra si lobservatoire etait a Greenwich (longitude 0 deg)
ra
+=
lGeoTime
.
n_lng
/
degrad
(
15.
);
//
JMA day;
//
day.Annee = year0;
//
day.Mois = month0;
//
day.Jour = iday0;
//
do {
//
hr += dt;
//
HMS hour = HtoHMS(hr);
//
HMS tsid;
//
tsid = TUtoTSid(day, hour);
//
r_8 ra = HMStoH(tsid); // le ra si lobservatoire etait a Greenwich (longitude 0 deg)
//
ra += lGeoTime.n_lng/degrad(15.);
dt
=
raSrc
-
ra
;
cout
<<
"METHOD NEW: meridien dt = "
<<
dt
<<
" "
<<
ToStringHdec
(
dt
)
<<
endl
;
//
dt = raSrc-ra;
//
cout << "METHOD NEW: meridien dt = " << dt << " " << ToStringHdec(dt) << endl;
}
while
(
++
i
<
MAXLOOPS
&&
fabs
(
dt
)
>
MAXERR
);
if
(
i
==
MAXLOOPS
)
throw
NotFoundExc
(
"Meridian date not found"
);
//
} while (++i < MAXLOOPS && fabs(dt) > MAXERR);
//
if (i == MAXLOOPS) throw NotFoundExc("Meridian date not found");
#undef MAXLOOPS
#undef MAXERR
//
#undef MAXLOOPS
//
#undef MAXERR
meridianT
.
SetHour
(
ToStringHdec
(
hr
));
//
meridianT.SetHour(ToStringHdec(hr));
}
//NEW
//
}//NEW
param
.
MeridianDate
=
meridianT
;
//JEC 18/2/19
TimeStamp
referenceT
(
referenceDate
);
param
.
ReferenceDate
=
referenceT
;
cout
<<
"APana initial parameter:
\n
"
<<
"Input dir: "
<<
param
.
inputDir
<<
"
\n
"
...
...
@@ -3184,7 +3253,7 @@ int main(int narg, char* arg[]) {
<<
param
.
raSrcHour
<<
":"
<<
param
.
raSrcMin
<<
":"
<<
param
.
raSrcSec
<<
"), DEC d:m:s ("
<<
param
.
decSrcDeg
<<
":"
<<
param
.
decSrcMin
<<
":"
<<
param
.
decSrcSec
<<
")
\n
"
<<
"Passage au meridien Nancay a: "
<<
param
.
MeridianDate
<<
endl
;
cout
<<
"Reference date a: "
<<
param
.
ReferenceDate
<<
endl
;
//JEC 18/2/19
cout
<<
" Action = "
<<
action
<<
" Type of cross REal or IMag: "
<<
param
.
Type
<<
endl
;
...
...
AnaFringes/makefile
View file @
eb539467
...
...
@@ -3,8 +3,8 @@
include
$(SOPHYABASE)/include/sophyamake.inc
### CCIN2P3
MINUITDIR = /sps/baoradio/JEC/PAON4/Minuit2-5.34.14/
MINUITDIR
=
/Users/campagne/Travail/kits/Minuit2-5.34.14/
MINUITDIR
=
/sps/baoradio/JEC/PAON4/Minuit2-5.34.14/
### MacJEC
MINUITDIR = /Users/campagne/Travail/kits/Minuit2-5.34.14/
MINUITLIB
=
-L
$(MINUITDIR)
/lib
MINUITINC
=
-I
$(MINUITDIR)
/include
...
...
Scripts/ckfitres.pic
View file @
eb539467
#########################################################################
### Script piapp pour verifier la sortie du programme trkacfit.cc
### Script piapp pour verifier la sortie du programme trkac
x
fit.cc
### Fit de parametres de D_dish, pointage des antennes, phases
### sur auto-cor et cross-cor
### Usage : exec ckfitres BaseFileName NTrk
### BaseFileName --> fichiers BaseFileName.ppf et BaseFileName_cx.ppf
#########################################################################
if ( $# < 2 ) then
echo 'ckfitres / missing arguments '
echo ' Usage ckfitres BaseFileName NTrk (or NSat) '
echo '
BaseFileName without .ppf, will open BaseFileName.ppf and BaseFileName_cx.ppf
'
return
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
delobjs *
if ( $# < 1 ) then
echo '
ckfitres.pic / missing arguments
'
echo ' Usage: exec ckfitres.pic PPFName_fit6cx '
endif
set flnm $1
set flnmA ${flnm}.ppf
set flnmB ${flnm}_cx.ppf
set NTrk $2
JJMAX = $NTrk+1
echo ' BaseFileName=' $flnm ' NTrk=' $NTrk ' --> JJMAX=' $JJMAX
delobjs *
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
echo '--- Opening ' $flnmA
openppf $flnmA
for j 1:$JJMAX
newwin 2 2 600 500
for i 1:5
disp ac_${i}_${j} 'black'
disp simac_${i}_${j} 'blue same'
defscript ckacx ' Check trkacxfit.cc result : Auto-Correlation and Cross-Cor fit for individual baselines ckacx BaseFileName NumSeq NTrk'
if ( $# < 3 ) then
echo 'ckacx / missing arguments '
echo ' Usage: ckacx BaseFileName NumSeq NbTrk (or NbSource) '
echo ' BaseFileName without .ppf, will open ac_NumSeq_BaseFileName.ppf and cx_NumSeq_BaseFileName.ppf'
return
endif
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
set flnm $1
set seq $2
set NTrk $3
JJMAX = $NTrk+1
set flnmA ac_${seq}_${flnm}.ppf
set flnmB cx_${seq}_${flnm}.ppf
echo ' BaseFileName=' $flnm ' NumSeq=' $seq ' NTrk=' $NTrk ' --> JJMAX=' $JJMAX
echo '---ckacx--- Opening ' $flnmA
openppf $flnmA
for j 1:$JJMAX
newwin 2 2 600 500
for i 1:5
disp ac_${i}_${j} 'black'
disp simac_${i}_${j} 'blue same'
end
end
echo '---ckacx--- Opening ' $flnmB
openppf $flnmB
for j 1:$JJMAX
newwin 3 2 800 500
for i 1:7
disp cx_${i}_${j} 'cdreal black'
disp simcx_${i}_${j} 'cdreal same blue'
disp cx_${i}_${j} 'cdimag same red'
disp simcx_${i}_${j} 'cdimag same orange'
end
end
end
echo '--- Opening ' $flnmB
openppf $flnmB
for j 1:$JJMAX
endscript
defscript ck6cx ' Check trkacxfit.cc result : Simultaneous 6 Cross-Cor fit : ck6cx PPFFileName NumSeq NTrk'
if ( $# < 3 ) then
echo 'ckacx / missing arguments '
echo ' Usage ckacx BaseFileName NumSeq NTrk (or NSource) '
echo ' BaseFileName without .ppf, will open ac_NumSeq_BaseFileName.ppf and cx_NumSeq_BaseFileName.ppf'
return
endif
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
defscript ck6 ' ck6 numseq (1...n) numtrak (1...n) '
set trk $2
set seq $1
echo '--- ck6 script: NumSeq=' $seq ' NumTrack=' $trk
newwin 3 2 800 500
for i 1:7
disp cx_${i}_${
j
} 'cdreal
black
'
disp simcx_${i}_${
j
} 'cdreal same blue'
disp cx_${i}_${
j
} 'cdimag same
red
'
disp simcx_${i}_${
j
} 'cdimag same orange'
for i 1:7
disp cx_${i}_${
trk}_${seq
} 'cdreal '
disp simcx_${i}_${
trk}_${seq
} 'cdreal same blue'
disp cx_${i}_${
trk}_${seq
} 'cdimag
red
same'
disp simcx_${i}_${
trk}_${seq
} 'cdimag same orange'
end
end
endscript
set flnm6cx $1
echo '----ckfitres.pic---- Opening file:' $flnm6cx
openppf $flnm6cx
echo '======= Use scripts ckacx ck6 to plot obs.visi / expected after fit :'
echo '---- Cx6Fit : ck6 numseq (1...n) numtrak (1...n)'
echo '---- AC,Cx Fits: ckacx BaseFileName NumSeq NbTrk (or NbSource) '
Scripts/p4sat.cmd
View file @
eb539467
######################################################################
#######
Exemples
de
commandes
pour
analyses
des
satellites
ds
PAON4
#######
Analyse
PAON4
,
########
1
ere
version
,
R
.
Ansari
,
Novembre
2018
########
1
ere
version
,
R
.
Ansari
,
Novembre
2018
###
MAJ
Decembre
2018
,
Fev
2018
###
Voir
les
commandes
pour
l
'analyse la plus recente en fin de fichier
######################################################################
#####################
...
...
@@ -125,3 +127,135 @@ ${PRGBASE}/Objs/tfm2dt ChezOlivier/Outputs/cor_B12_6dec18.ppf 1275,1282 dt_B12.p
$
{
PRGBASE
}
/Objs/trkacfit -D
4
.5
-sdec
10
-ckf
tstB12
.ppf
-out
B12_fac
.txt
dt_B12
,
trk_B12_43055
,
1178
,
1215
,
1278
.5
dt_B12
,
trk_B12_37847
,
1665
,
1695
,
1278
.5
##################################################################################################
#####
Fevrier
2019
################
###
Analyse
acquisition
longue
duree
Janv
2019
,
a
partir
du
fichier
LD_goodsat
.txt
d
'Olivier
##################################################################################################
### AU CC-IN2P3 : se mettre ds /sps/baoradio/Reza/AnaP4/SatLD14jan19
source /pbs/throng/baoradio/Library/Sophya/env.csh
cd /sps/baoradio/Reza/AnaP4/SatLD14jan19
set PRGBASE = /sps/baoradio/Reza/AnaP4/AnaPAON4
### Pour compiler le code au CC
cd $PRGBASE
git pull
make
## Pour compiler avec MINUIT
setenv MINUITDIR /pbs/throng/baoradio/Library/ExtLibs/pro/
make mn_trkacxfit
########################
cd /sps/baoradio/Reza/AnaP4/SatLD14jan19
## Extraction du signal pour les satellites entre 1275-1282 MHz
${PRGBASE}/Objs/tfm2dt ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1275,1282 dt_LD14jan19.ppf
## Temps en minutes 1115,1160
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2019/01/14 18:56:00" -H 0.,79.08 -K 40889 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_40889_20190114.txt trk_40889_20190114.ppf 2019/01/14
## Temps en minutes 1460,1500
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2019/01/15 0:43:00" -H 0.,79.08 -K 43566 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_43566_20190115.txt trk_43566_20190115.ppf 2019/01/14
## Temps en minutes 2455,2495
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2019/01/15 17:16:00" -H 0.,79.08 -K 43055 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_43055_20190115.txt trk_43055_20190115.ppf 2019/01/14
## Temps en minutes 4130,4175
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2019/01/16 21:15:00" -H 180.,83.35 -K 43057 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_43057_20190116.txt trk_43057_20190116.ppf 2019/01/14
## Temps en minutes 6120,6160
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2019/01/18 6:20:00" -H 180.,83.35 -K 41175 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_41175_20190118.txt trk_41175_20190118.ppf 2019/01/14
## Temps en minutes 11520+720,+760 : 12240,12280
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2019/01/22 12:20:00" -H 180.,55 -K 40890 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_40890_20190122.txt trk_40890_20190122.ppf 2019/01/14
## Temps en minutes 11520+850,+890 : 12370,12410
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2019/01/22 14:35:00" -H 180.,55 -K 41859 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_41859_20190122.txt trk_41859_20190122.ppf 2019/01/14
## Extraction du signal pour les sources du ciel entre 1342-1350 MHz -> FreqCentral=1346 MHz
${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1342,1350 dtsky1346_LD14jan19.ppf
## Extraction du signal pour les sources du ciel entre 1392-1400 MHz -> FreqCentral=1396 MHz
${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
#### 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 )
${PRGBASE}/Satellites/Objs/trk2dt $f.txt $f.ppf 2019/01/14
end
foreach f ( M1_survey_2019-1-19 M1_survey_2019-1-20 M87_survey_2019-1-21 M87_survey_2019-1-22 M87_survey_2019-1-23 )
${PRGBASE}/Satellites/Objs/trk2dt $f.txt $f.ppf 2019/01/14
end
##### Intervalle de temps de transit (en minutes, since 14/01/2019 0h00)
CasA_survey_2019-1-14 : 920,960
CasA_survey_2019-1-15 : 2355,2395
CasA_survey_2019-1-16 : 3790,3830
CygA_survey_2019-1-17 : 5025,5065
CygA_survey_2019-1-18 : 6460,6500
#####
#########################
####### Fichiers .d : Tracks-list
#### ldtrk_casA.d
--- Observation Longue duree Janvier 2019
Sources avec observation vers CasA (z=+11.42)
@zenang 10
@trk dt_LD14jan19 1115,1160 1278.5 trk_40889_20190114
@trk dt_LD14jan19 1460,1500 1278.5 trk_43566_20190115
@trk dt_LD14jan19 2455,2495 1278.5 trk_43055_20190115
# Sky at 1346 MHz
@trk dtsky1346_LD14jan19 920,960 1346 CasA_survey_2019-1-14
@trk dtsky1346_LD14jan19 2355,2395 1346 CasA_survey_2019-1-15
@trk dtsky1346_LD14jan19 3790,3830 1346 CasA_survey_2019-1-16
# Sky at 1396 MHz
@trk dtsky1396_LD14jan19 920,960 1396 CasA_survey_2019-1-14
@trk dtsky1396_LD14jan19 2355,2395 1396 CasA_survey_2019-1-15
@trk dtsky1396_LD14jan19 3790,3830 1396 CasA_survey_2019-1-16
# Sky at 1416 MHz
@trk dtsky1416_LD14jan19 920,960 1416 CasA_survey_2019-1-14
@trk dtsky1416_LD14jan19 2355,2395 1416 CasA_survey_2019-1-15
@trk dtsky1416_LD14jan19 3790,3830 1416 CasA_survey_2019-1-16
#### ldtrk_cygA.d
--- Observation Longue duree Janvier 2019
Sources avec observation vers CygA (z=-6.5)
@zenang -6.
@trk dt_LD14jan19 4130,4175 1278.5 trk_43057_20190116
@trk dt_LD14jan19 6120,6160 1278.5 trk_41175_20190118
@trk dtsky1346_LD14jan19 5025,5065 1346 CygA_survey_2019-1-17
@trk dtsky1346_LD14jan19 6460,6500 1346 CygA_survey_2019-1-18
@trk dtsky1396_LD14jan19 5025,5065 1396 CygA_survey_2019-1-17
@trk dtsky1396_LD14jan19 6460,6500 1396 CygA_survey_2019-1-18
@trk dtsky1416_LD14jan19 5025,5065 1416 CygA_survey_2019-1-17
@trk dtsky1416_LD14jan19 6460,6500 1416 CygA_survey_2019-1-18
#### ldtrk_virA.d
--- Observation Longue duree Janvier 2019
Sources avec observation vers VirgoA (z=-35)
@zenang -35
@trk dt_LD14jan19 12240,12280 1278.5 trk_40890_20190122
@trk dt_LD14jan19 12370,12410 1278.5 trk_41859_20190122
##### Exemple d'
execution
de
commandes
$
{
PRGBASE
}
/Objs/trkacxfit -docx
6
s
-docx
6
f
-D
4
.5
-out
ldj
.txt
-ckf
ldj
.ppf
ldtrk_casA
.d
ldtrk_cygA
.d
ldtrk_virA
.d
#
Ou
$
{
PRGBASE
}
/Objs/trkacxfit -docx
6
s
-docx
6
f
-D
4
.5
-out
casA
.txt
-ckf
casA
.ppf
ldtrk_casA
.d
$
{
PRGBASE
}
/Objs/trkacxfit -nocxf -D
4
.5
-out
casA
.txt
-ckf
casA
.ppf
ldtrk_casA
.d
####
Utiliser
le
script
piapp
:
ckfitres
.pic
et
les
scripts
definis
dedans
pour
verifier
les
resultats