- Timestamp:
- Nov 28, 2011, 10:04:36 AM (13 years ago)
- Location:
- BAORadio/AmasNancay/trunk
- Files:
-
- 4 added
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/etude_mergeAna.pic
r607 r611 1 2 set path "/sps/baoradio/AmasNancay/AST" 1 3 2 4 set source $1 … … 9 11 10 12 set fcalib 1410 11 set f1420low 1420.2 12 set f1420high 1420.6 13 14 13 14 # Abell 85, 500 cycles 15 if ( ${source} == "abell85" ) then 16 set f1420low 1420.2 17 set f1420high 1420.6 18 set packVal "1 10 25 50 100" 19 endif 20 21 # Abell 1205, 600 cycles 22 if ( ${source} == "abell1205" ) then 23 set f1420low 1420.1 24 set f1420high 1420.6 25 set packVal "1 10 25 50 100 120" 26 endif 27 28 # Abell 2440, 210 cycles 29 if ( ${source} == "abell2440" ) then 30 set f1420low 1420.4 31 set f1420high 1420.7 32 set packVal "1 10 25 42" 33 endif 34 35 echo "f1420low=${f1420low}, f1420high=${f1420high}" 36 echo "Packing values="${packVal} 15 37 16 38 nbins1420 = floor((${f1420high}-${f1420low})/250*8192) 17 18 39 echo $nbins1420 19 20 40 21 41 clearscript diffOnOff … … 26 46 clearscript sigmaRaw1400a1420 27 47 clearscript evolAt1410a1415OffOn 48 28 49 ############################################################################ 29 50 defscript evolAtCalibFreq … … 90 111 91 112 92 graphicatt "xylimits=0, 500,-0.01,0.01"113 graphicatt "xylimits=0,${ncycles},-0.02,0.02" 93 114 94 115 newwin 1 2 … … 159 180 line2vec vecCh1All $linCh1All 160 181 161 #packing values only valid for 500 cycles...162 set packVal "1 10 25 50 100"163 182 set sigmaCh0 "" 164 183 set errsigCh0 "" 165 184 set sigmaCh1 "" 166 185 set errsigCh1 "" 167 168 186 169 187 foreach ipack ( $packVal ) … … 236 254 line2vec vecCh1All $linCh1All 237 255 238 #packing values only valid for 500 cycles...239 set packVal "1 10 25 50 100"240 256 set sigmaCh0 "" 241 257 set errsigCh0 "" … … 329 345 330 346 331 graphicatt "xylimits=0, 500,-0.01,0.01"347 graphicatt "xylimits=0,${ncycles},-0.01,0.01" 332 348 333 349 newwin 1 2 … … 394 410 395 411 396 graphicatt "xylimits=0, 500,-0.1,0.1"412 graphicatt "xylimits=0,${ncycles},-0.1,0.1" 397 413 398 414 newwin 1 2 … … 487 503 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 488 504 489 490 491 505 graphicatt "xylimits=1250,1500,-0.01,0.01" 492 506 del snocal0 snocal1 … … 499 513 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 500 514 501 del sruncal0 sruncal1502 objaoper meanPerRunCalib row 0 sruncal0503 objaoper meanPerRunCalib row 1 sruncal1504 newwin 1 1505 plot2d sruncal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"506 plot2d sruncal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"507 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt508 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt509 510 newwin 1 1511 del scyclecal0 scyclecal1512 objaoper meanPerCycleCalib row 0 scyclecal0513 objaoper meanPerCycleCalib row 1 scyclecal1514 plot2d scyclecal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"515 plot2d scyclecal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"516 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt517 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt515 #del sruncal0 sruncal1 516 #objaoper meanPerRunCalib row 0 sruncal0 517 #objaoper meanPerRunCalib row 1 sruncal1 518 #newwin 1 1 519 #plot2d sruncal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta" 520 #plot2d sruncal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta" 521 #settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 522 #setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 523 524 #newwin 1 1 525 #del scyclecal0 scyclecal1 526 #objaoper meanPerCycleCalib row 0 scyclecal0 527 #objaoper meanPerCycleCalib row 1 scyclecal1 528 #plot2d scyclecal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta" 529 #plot2d scyclecal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta" 530 #settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 531 #setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 518 532 519 533 … … 526 540 #use minuscule for source label 527 541 528 529 530 531 openppf onoffsurvey_${source}-${ncycles}Cycles.ppf 542 openppf ${path}/onoffsurvey_${source}-${ncycles}Cycles.ppf 532 543 print onoffevol 544 545 #openppf ${path}/rawOnOffDiff_${source}-${ncycles}Cycles.ppf 546 #ls 533 547 534 548 set defatt "font=helvetica,bold,20 fixedfonsize" … … 546 560 #evolAt1420SideFreq 547 561 #IntAt1420 562 #sigmaRaw1420Side 563 #sigmaRaw1400a1420 564 #evolAt1400a1420OffOn -
BAORadio/AmasNancay/trunk/mergeAnaFiles.cc
r610 r611 330 330 331 331 //mean ON-OFF over the list of cycles 332 TMatrix<r_4> meanDiffONOFFovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); 332 TMatrix<r_4> meanDiffONOFFovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); //set to 0 333 333 TMatrix<r_4> meanDiffONOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); //set to 0 334 334 TMatrix<r_4> meanDiffONOFF_perRunCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); //set to 0 … … 397 397 string srcLower = tokens2[0]; 398 398 399 400 401 399 PInPersist fin(*iFile); 402 400 vector<string> vec = fin.GetNameTags(); … … 710 708 //JEC 29/10/11 add ON-OFF/OFF 711 709 TMatrix<r_4> diffOnOffOvOff_noCalib(diffOnOff_noCalib,false); //do not share data 712 TMatrix<r_4> aSpecOffFi tltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);710 TMatrix<r_4> aSpecOffFiltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); 713 711 sa_size_t halfWidth = 35; //number of freq. bin for the 1/2 width of the filtering window 714 medianFiltering(aSpecOff,halfWidth,aSpecOffFi tltered);715 716 diffOnOffOvOff_noCalib.Div(aSpecOffFi tltered); //division in place712 medianFiltering(aSpecOff,halfWidth,aSpecOffFiltered); 713 714 diffOnOffOvOff_noCalib.Div(aSpecOffFiltered); //division in place 717 715 meanDiffONOFFovOFF_noCalib += diffOnOffOvOff_noCalib; 718 716 … … 723 721 meanDiffONOFF_perCycleCalib += diffOnOff_perCycleCalib; 724 722 725 //726 723 totalNumberCycles++; 724 727 725 //Fill NTuple 728 726 xnt[0] = totalNumberCycles; … … 798 796 //JEC 18/11/11 follow up the 1400-1420MHz OFF only 799 797 TMatrix<r_4> aSpecOffovOff(aSpecOff,false); 800 aSpecOffovOff.Div(aSpecOffFi tltered);798 aSpecOffovOff.Div(aSpecOffFiltered); 801 799 802 800 TVector<r_4> meanInRange_1410a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS); … … 808 806 809 807 TMatrix<r_4> aSpecOnovOff(aSpecOn,false); 810 aSpecOnovOff.Div(aSpecOffFi tltered);808 aSpecOnovOff.Div(aSpecOffFiltered); 811 809 812 810 TVector<r_4> meanInRange_1410a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS); … … 886 884 887 885 tag = "onoffevol"; 888 fos << PPFNameTag(tag) << onoffevolution; 886 fos << PPFNameTag(tag) << onoffevolution; 887 889 888 }//end of save 890 889 … … 912 911 TMatrix<r_4> meanONovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); // ON/Filtered_OFF 913 912 TMatrix<r_4> meanOFFovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); // OFF/Filtered_OFF 914 915 913 //Tuple only for RAW things to follow 916 914 static const int NINFO=11; … … 953 951 for (iFile = listOfFiles.begin(); iFile != iFileEnd; ++iFile) { 954 952 if (para.debuglev_>90){ 955 } 953 cout << "load file <" << *iFile << ">" << endl; 954 } 955 956 vector<string> tokens; 957 split(*iFile,"_",tokens); 958 string dateOfRun = tokens[1]; 959 if (para.debuglev_>90){ 960 cout << "date <" << dateOfRun << ">" << endl; 961 } 962 vector<string> tokens2; 963 split(tokens[2],".",tokens2); 964 string srcLower = tokens2[0]; 965 956 966 PInPersist fin(*iFile); 957 967 vector<string> vec = fin.GetNameTags(); … … 983 993 regexp(iSpec->c_str(),matchstr.c_str(),&b,&e); 984 994 if (para.debuglev_>90){ 985 cout << " sp actra <" << *iSpec << ">" << endl;995 cout << " spectra <" << *iSpec << ">" << endl; 986 996 cout << " cycle " << iSpec->substr(b) << endl; 987 997 } … … 1011 1021 for (list<int>::iterator ic=commonCycles.begin(); ic!=commonCycles.end(); ++ic){ 1012 1022 1023 // AST 28.11.11 remove non-calibrated cycles for Abell1205 and Abell2440 1024 if ( *ic == 1 && srcLower == "abell1205" ) { 1025 if ( dateOfRun == "20110502" || dateOfRun == "20110607" || dateOfRun == "20110818" ) { 1026 cout << "Skipping non-calibrated cycle " << *ic << endl; 1027 continue; 1028 } 1029 } else if ( *ic == 1 && srcLower == "abell2440" ) { 1030 if ( dateOfRun == "20110516" ) { 1031 cout << "Skipping non-calibrated cycle " << *ic << endl; 1032 continue; 1033 } 1034 } else if ( *ic == 3 && srcLower == "abell1205" ) { 1035 if ( dateOfRun == "20110810" ) { 1036 cout << "Skipping non-calibrated cycle " << *ic << endl; 1037 continue; 1038 } 1039 } 1040 1013 1041 string ppftag; 1014 1042 //load ON phase … … 1026 1054 //Perform the difference ON-OFF 1027 1055 TMatrix<r_4> diffOnOff_noCalib = aSpecOn - aSpecOff; 1056 1028 1057 meanDiffONOFF_noCalib += diffOnOff_noCalib; 1029 1058 1030 1059 //JEC 29/10/11 add ON-OFF/OFF 1031 1060 TMatrix<r_4> diffOnOffOvOff_noCalib(diffOnOff_noCalib,false); //do not share data 1032 TMatrix<r_4> aSpecOffFi tltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);1061 TMatrix<r_4> aSpecOffFiltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); 1033 1062 sa_size_t halfWidth = 35; //number of freq. bin for the 1/2 width of the filtering window 1034 medianFiltering(aSpecOff,halfWidth,aSpecOffFi tltered);1035 1036 diffOnOffOvOff_noCalib.Div(aSpecOffFi tltered); //division in place1063 medianFiltering(aSpecOff,halfWidth,aSpecOffFiltered); 1064 1065 diffOnOffOvOff_noCalib.Div(aSpecOffFiltered); //division in place 1037 1066 meanDiffONOFFovOFF_noCalib += diffOnOffOvOff_noCalib; 1038 1039 1067 1040 1068 //JEC 15/11/11 add ON/OFF and OFF/OFF 1041 1069 TMatrix<r_4> onOvOff(aSpecOn,false); 1042 onOvOff.Div(aSpecOffFi tltered);1070 onOvOff.Div(aSpecOffFiltered); 1043 1071 meanONovOFF_noCalib += onOvOff; 1044 1072 1045 1073 TMatrix<r_4> offOvOff(aSpecOff,false); 1046 offOvOff.Div(aSpecOffFi tltered);1074 offOvOff.Div(aSpecOffFiltered); 1047 1075 meanOFFovOFF_noCalib += offOvOff; 1048 1076 1049 1077 totalNumberCycles++; 1050 1051 1078 1052 1079 //Fill NTuple 1053 1080 xnt[0] = totalNumberCycles; 1054 1081 1055 1056 1082 //Follow up arround the 1420.4MHz Freq. 1057 1083 TVector<r_4> meanInRange_1420Freq_noCalib(NUMBER_OF_CHANNELS); … … 1079 1105 //JEC 18/11/11 follow up the 1400-1420MHz OFF only 1080 1106 TMatrix<r_4> aSpecOffovOff(aSpecOff,false); 1081 aSpecOffovOff.Div(aSpecOffFi tltered);1107 aSpecOffovOff.Div(aSpecOffFiltered); 1082 1108 1083 1109 TVector<r_4> meanInRange_1410a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS); … … 1088 1114 1089 1115 TMatrix<r_4> aSpecOnovOff(aSpecOn,false); 1090 aSpecOnovOff.Div(aSpecOffFi tltered);1116 aSpecOnovOff.Div(aSpecOffFiltered); 1091 1117 1092 1118 TVector<r_4> meanInRange_1410a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS); … … 1136 1162 1137 1163 tag = "onoffevol"; 1138 fos << PPFNameTag(tag) << onoffevolution; 1164 fos << PPFNameTag(tag) << onoffevolution; 1139 1165 1140 1166 }//end save -
BAORadio/AmasNancay/trunk/plotRawDiffOnOff.pic
r610 r611 3 3 set ncycles $2 4 4 5 c++compile meanSigSum6 c++link meanSigSum.so meanSigSum5 #c++compile meanSigSum 6 #c++link meanSigSum.so meanSigSum 7 7 8 8 set defatt "font=helvetica,bold,20 fixedfonsize" … … 20 20 nbins1420 = floor((${f1420high}-${f1420low})/250*8192) 21 21 22 23 openppf rawOnOffDiff_${source}-${ncycles}Cycles.ppf22 set path "/sps/baoradio/AmasNancay/AST" 23 openppf ${path}/rawOnOffDiff_${source}-${ncycles}Cycles.ppf 24 24 25 25 … … 226 226 n/pl onoffevol.(onRaw0f14101415-offRaw0f14101415)%cycle ! ! "blue cpts notit nsta" 227 227 n/pl onoffevol.(onRaw1f14101415-offRaw1f14101415)%cycle ! ! "same red cpts notit nsta" 228 settitle "Raw (O FF-ON)/OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt228 settitle "Raw (ON-OFF)/OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 229 229 setaxelabels "cycle" "DI [1410,1415]MHz Raw (a.u)" $axedefatt 230 230 … … 301 301 plot2d s0 (n/8192)*250+1250 val n>0 "cpts blue notit nsta" 302 302 plot2d s1 (n/8192)*250+1250 val n>0 "same cpts red notit nsta" 303 settitle " Raw (ON-OFF)/OFF ${source} ${ncycles} cycles, Ch 0 (blue) Ch 1 (red)"303 settitle "plot Raw (ON-OFF)/OFF ${source} ${ncycles} cycles, Ch 0 (blue) Ch 1 (red)" 304 304 setaxelabels "Freq. (MHz)" "I (a.u)" 305 305 -
BAORadio/AmasNancay/trunk/proc_specmfib.sh
r600 r611 567 567 568 568 #save results to Irods (use option -f ONLY to force override) 569 emptydir=`$FIND ${ curDir}/${dirName} -maxdepth 1 -type f | $WC -l`569 emptydir=`$FIND ${outDir} -maxdepth 1 -type f | $WC -l` 570 570 if [ "<${emptydir}>" != "<0>" ]; then 571 571 iput -v -K -r -f ${curDir}/${dirName} /baoradio/data/AmasNancay/${localpath} 572 572 else 573 $ECHO "DEBUG (${scriptName}): directory ${ curDir}/${dirName} is empty"573 $ECHO "DEBUG (${scriptName}): directory ${outDir} is empty" 574 574 fi 575 575 #
Note: See TracChangeset
for help on using the changeset viewer.