Changeset 611

Timestamp:

Nov 28, 2011, 10:04:36 AM (13 years ago)

Author:

torrento

Message:

Add pic analysis files and emptydir function for proc_specmfib

Location:

BAORadio/AmasNancay/trunk

Files:

• dabchrom_plot.pic (added)
• emptydir.sh (added)
• etude_dab.pic (added)
• etude_gain_medfiltmtx.pic (added)
• etude_mergeAna.pic (modified) (12 diffs)
• mergeAnaFiles.cc (modified) (15 diffs)
• plotRawDiffOnOff.pic (modified) (4 diffs)
• proc_specmfib.sh (modified) (1 diff)

BAORadio/AmasNancay/trunk/etude_mergeAna.pic

@@ +1 @@
+
+set path "/sps/baoradio/AmasNancay/AST"
 
 set source  $1
@@ -9 +11 @@
 
 set fcalib 1410
-set f1420low 1420.2
-set f1420high 1420.6
-
-
+
+# Abell 85, 500 cycles
+if ( ${source} == "abell85" ) then
+  set f1420low 1420.2
+  set f1420high 1420.6
+  set packVal "1 10 25 50 100"
+endif
+
+# Abell 1205, 600 cycles
+if ( ${source} == "abell1205" ) then
+  set f1420low 1420.1
+  set f1420high 1420.6
+  set packVal "1 10 25 50 100 120"
+endif
+
+# Abell 2440, 210 cycles
+if ( ${source} == "abell2440" ) then
+  set f1420low 1420.4
+  set f1420high 1420.7
+  set packVal "1 10 25 42"
+endif
+
+echo "f1420low=${f1420low}, f1420high=${f1420high}"
+echo "Packing values="${packVal}
 
 nbins1420 = floor((${f1420high}-${f1420low})/250*8192)
-
 echo $nbins1420
-
 
 clearscript diffOnOff
@@ -26 +46 @@
 clearscript sigmaRaw1400a1420
 clearscript evolAt1410a1415OffOn
+
 ############################################################################
 defscript evolAtCalibFreq
@@ -90 +111 @@
 
 
-graphicatt "xylimits=0,500,-0.01,0.01"
+graphicatt "xylimits=0,${ncycles},-0.02,0.02"
 
 newwin 1 2
@@ -159 +180 @@
 line2vec vecCh1All $linCh1All
 
-#packing values only valid for 500 cycles...
-set packVal "1 10 25 50 100"
 set sigmaCh0 ""
 set errsigCh0 ""
 set sigmaCh1 ""
 set errsigCh1 ""
-
 
 foreach ipack ( $packVal )
@@ -236 +254 @@
 line2vec vecCh1All $linCh1All
 
-#packing values only valid for 500 cycles...
-set packVal "1 10 25 50 100"
 set sigmaCh0 ""
 set errsigCh0 ""
@@ -329 +345 @@
 
 
-graphicatt "xylimits=0,500,-0.01,0.01"
+graphicatt "xylimits=0,${ncycles},-0.01,0.01"
 
 newwin 1 2
@@ -394 +410 @@
 
 
-graphicatt "xylimits=0,500,-0.1,0.1"
+graphicatt "xylimits=0,${ncycles},-0.1,0.1"
 
 newwin 1 2
@@ -487 +503 @@
 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
 
-
-
 graphicatt "xylimits=1250,1500,-0.01,0.01"
 del snocal0 snocal1
@@ -499 +513 @@
 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
 
-del sruncal0 sruncal1
-objaoper meanPerRunCalib row 0 sruncal0
-objaoper meanPerRunCalib row 1 sruncal1
-newwin 1 1
-plot2d sruncal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
-plot2d sruncal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
-settitle "ON-OFF Calib Run  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
-setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
-
-newwin 1 1
-del scyclecal0 scyclecal1
-objaoper meanPerCycleCalib row 0 scyclecal0
-objaoper meanPerCycleCalib row 1 scyclecal1
-plot2d scyclecal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
-plot2d scyclecal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
-settitle "ON-OFF Calib Cycle  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
-setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
+#del sruncal0 sruncal1
+#objaoper meanPerRunCalib row 0 sruncal0
+#objaoper meanPerRunCalib row 1 sruncal1
+#newwin 1 1
+#plot2d sruncal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
+#plot2d sruncal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
+#settitle "ON-OFF Calib Run  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
+#setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
+
+#newwin 1 1
+#del scyclecal0 scyclecal1
+#objaoper meanPerCycleCalib row 0 scyclecal0
+#objaoper meanPerCycleCalib row 1 scyclecal1
+#plot2d scyclecal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
+#plot2d scyclecal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
+#settitle "ON-OFF Calib Cycle  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
+#setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
 
 
@@ -526 +540 @@
 #use minuscule for source label
 
-
-
-
-openppf onoffsurvey_${source}-${ncycles}Cycles.ppf
+openppf ${path}/onoffsurvey_${source}-${ncycles}Cycles.ppf
 print onoffevol
+
+#openppf ${path}/rawOnOffDiff_${source}-${ncycles}Cycles.ppf
+#ls
 
 set defatt "font=helvetica,bold,20 fixedfonsize"
@@ -546 +560 @@
 #evolAt1420SideFreq
 #IntAt1420
+#sigmaRaw1420Side
+#sigmaRaw1400a1420
+#evolAt1400a1420OffOn

BAORadio/AmasNancay/trunk/mergeAnaFiles.cc

@@ -330 +330 @@
 
   //mean ON-OFF over the list of cycles
-  TMatrix<r_4> meanDiffONOFFovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);          //set to 0 
+  TMatrix<r_4> meanDiffONOFFovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);     //set to 0 
   TMatrix<r_4> meanDiffONOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);          //set to 0 
   TMatrix<r_4> meanDiffONOFF_perRunCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);      //set to 0
@@ -397 +397 @@
     string srcLower = tokens2[0];
 
-
-    
     PInPersist fin(*iFile);
     vector<string> vec = fin.GetNameTags();
@@ -710 +708 @@
       //JEC 29/10/11 add ON-OFF/OFF
       TMatrix<r_4> diffOnOffOvOff_noCalib(diffOnOff_noCalib,false); //do not share data
-      TMatrix<r_4> aSpecOffFitltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);
+      TMatrix<r_4> aSpecOffFiltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);
       sa_size_t halfWidth = 35; //number of freq. bin for the 1/2 width of the filtering window
-      medianFiltering(aSpecOff,halfWidth,aSpecOffFitltered);
-      
-      diffOnOffOvOff_noCalib.Div(aSpecOffFitltered); //division in place
+      medianFiltering(aSpecOff,halfWidth,aSpecOffFiltered);
+      
+      diffOnOffOvOff_noCalib.Div(aSpecOffFiltered); //division in place
       meanDiffONOFFovOFF_noCalib += diffOnOffOvOff_noCalib;
 
@@ -723 +721 @@
       meanDiffONOFF_perCycleCalib += diffOnOff_perCycleCalib;
 
-      //
       totalNumberCycles++;
+
       //Fill NTuple
       xnt[0] = totalNumberCycles;
@@ -798 +796 @@
       //JEC 18/11/11 follow up the 1400-1420MHz OFF only
       TMatrix<r_4> aSpecOffovOff(aSpecOff,false);
-      aSpecOffovOff.Div(aSpecOffFitltered);
+      aSpecOffovOff.Div(aSpecOffFiltered);
 
       TVector<r_4> meanInRange_1410a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS);
@@ -808 +806 @@
 
       TMatrix<r_4> aSpecOnovOff(aSpecOn,false);
-      aSpecOnovOff.Div(aSpecOffFitltered);
+      aSpecOnovOff.Div(aSpecOffFiltered);
 
       TVector<r_4> meanInRange_1410a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS);
@@ -886 +884 @@
     
     tag = "onoffevol";
-    fos << PPFNameTag(tag) << onoffevolution;    
+    fos << PPFNameTag(tag) << onoffevolution;
+ 
   }//end of save
 
@@ -912 +911 @@
   TMatrix<r_4> meanONovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); //  ON/Filtered_OFF
   TMatrix<r_4> meanOFFovOFF_noCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); // OFF/Filtered_OFF
-
   //Tuple only for RAW things to follow 
   static const int NINFO=11;
@@ -953 +951 @@
   for (iFile = listOfFiles.begin(); iFile != iFileEnd; ++iFile) {
     if (para.debuglev_>90){
-    }
+      cout << "load file <" << *iFile << ">" << endl;
+    }
+
+    vector<string> tokens;
+    split(*iFile,"_",tokens);
+    string dateOfRun = tokens[1];
+    if (para.debuglev_>90){
+      cout << "date <" << dateOfRun << ">" << endl;
+    }
+    vector<string> tokens2;
+    split(tokens[2],".",tokens2);
+    string srcLower = tokens2[0];
+
     PInPersist fin(*iFile);
     vector<string> vec = fin.GetNameTags();
@@ -983 +993 @@
         regexp(iSpec->c_str(),matchstr.c_str(),&b,&e);
         if (para.debuglev_>90){
-          cout << " spactra <" << *iSpec << ">" << endl;
+          cout << " spectra <" << *iSpec << ">" << endl;
           cout << " cycle " << iSpec->substr(b) << endl;
         }
@@ -1011 +1021 @@
     for (list<int>::iterator ic=commonCycles.begin(); ic!=commonCycles.end(); ++ic){
       
+      // AST 28.11.11 remove non-calibrated cycles for Abell1205 and Abell2440
+      if ( *ic == 1 && srcLower == "abell1205" ) {
+          if ( dateOfRun == "20110502" || dateOfRun == "20110607" || dateOfRun == "20110818" ) {
+            cout << "Skipping non-calibrated cycle " << *ic << endl;
+            continue;
+          }
+      } else if ( *ic == 1 && srcLower == "abell2440" ) {
+          if ( dateOfRun == "20110516" ) {
+            cout << "Skipping non-calibrated cycle " << *ic << endl;
+            continue;
+          }
+      } else if ( *ic == 3 && srcLower == "abell1205" ) {
+          if ( dateOfRun == "20110810" ) {
+            cout << "Skipping non-calibrated cycle " << *ic << endl;
+            continue;
+          }
+      }
+
       string ppftag;
       //load ON phase
@@ -1026 +1054 @@
       //Perform the difference ON-OFF
       TMatrix<r_4> diffOnOff_noCalib = aSpecOn - aSpecOff;
+
       meanDiffONOFF_noCalib += diffOnOff_noCalib;
       
       //JEC 29/10/11 add ON-OFF/OFF
       TMatrix<r_4> diffOnOffOvOff_noCalib(diffOnOff_noCalib,false); //do not share data
-      TMatrix<r_4> aSpecOffFitltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);
+      TMatrix<r_4> aSpecOffFiltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);
       sa_size_t halfWidth = 35; //number of freq. bin for the 1/2 width of the filtering window
-      medianFiltering(aSpecOff,halfWidth,aSpecOffFitltered);
-      
-      diffOnOffOvOff_noCalib.Div(aSpecOffFitltered); //division in place
+      medianFiltering(aSpecOff,halfWidth,aSpecOffFiltered);
+      
+      diffOnOffOvOff_noCalib.Div(aSpecOffFiltered); //division in place
       meanDiffONOFFovOFF_noCalib += diffOnOffOvOff_noCalib;
-
 
       //JEC 15/11/11 add ON/OFF and OFF/OFF
       TMatrix<r_4> onOvOff(aSpecOn,false);
-      onOvOff.Div(aSpecOffFitltered);
+      onOvOff.Div(aSpecOffFiltered);
       meanONovOFF_noCalib += onOvOff;
       
       TMatrix<r_4> offOvOff(aSpecOff,false);
-      offOvOff.Div(aSpecOffFitltered);
+      offOvOff.Div(aSpecOffFiltered);
       meanOFFovOFF_noCalib += offOvOff;
 
       totalNumberCycles++;
-
 
       //Fill NTuple
       xnt[0] = totalNumberCycles;
  
-
       //Follow up arround the 1420.4MHz Freq.
       TVector<r_4> meanInRange_1420Freq_noCalib(NUMBER_OF_CHANNELS);
@@ -1079 +1105 @@
       //JEC 18/11/11 follow up the 1400-1420MHz OFF only
       TMatrix<r_4> aSpecOffovOff(aSpecOff,false);
-      aSpecOffovOff.Div(aSpecOffFitltered);
+      aSpecOffovOff.Div(aSpecOffFiltered);
 
       TVector<r_4> meanInRange_1410a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS);
@@ -1088 +1114 @@
 
       TMatrix<r_4> aSpecOnovOff(aSpecOn,false);
-      aSpecOnovOff.Div(aSpecOffFitltered);
+      aSpecOnovOff.Div(aSpecOffFiltered);
 
       TVector<r_4> meanInRange_1410a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS);
@@ -1136 +1162 @@
 
     tag = "onoffevol";
-    fos << PPFNameTag(tag) << onoffevolution;    
+    fos << PPFNameTag(tag) << onoffevolution;
 
   }//end save

BAORadio/AmasNancay/trunk/plotRawDiffOnOff.pic

@@ -3 +3 @@
 set ncycles $2
 
-c++compile meanSigSum
-c++link    meanSigSum.so meanSigSum
+#c++compile meanSigSum
+#c++link    meanSigSum.so meanSigSum
 
 set defatt "font=helvetica,bold,20 fixedfonsize"
@@ -20 +20 @@
 nbins1420 = floor((${f1420high}-${f1420low})/250*8192)
 
-
-openppf rawOnOffDiff_${source}-${ncycles}Cycles.ppf
+set path "/sps/baoradio/AmasNancay/AST"
+openppf ${path}/rawOnOffDiff_${source}-${ncycles}Cycles.ppf
 
 
@@ -226 +226 @@
 n/pl onoffevol.(onRaw0f14101415-offRaw0f14101415)%cycle ! ! "blue cpts notit nsta"
 n/pl onoffevol.(onRaw1f14101415-offRaw1f14101415)%cycle ! ! "same red cpts notit nsta"
-settitle "Raw (OFF-ON)/OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
+settitle "Raw (ON-OFF)/OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
 setaxelabels "cycle" "DI [1410,1415]MHz Raw (a.u)" $axedefatt
 
@@ -301 +301 @@
 plot2d s0 (n/8192)*250+1250 val n>0 "cpts blue notit nsta"
 plot2d s1 (n/8192)*250+1250 val n>0 "same cpts red notit nsta"
-settitle "Raw (ON-OFF)/OFF ${source} ${ncycles} cycles, Ch 0 (blue) Ch 1 (red)"
+settitle "plot Raw (ON-OFF)/OFF ${source} ${ncycles} cycles, Ch 0 (blue) Ch 1 (red)"
 setaxelabels "Freq. (MHz)" "I (a.u)"
 

BAORadio/AmasNancay/trunk/proc_specmfib.sh

@@ -567 +567 @@
 
 #save results to Irods (use option -f ONLY to force override)
-    emptydir=`$FIND ${curDir}/${dirName} -maxdepth 1 -type f | $WC -l`
+    emptydir=`$FIND ${outDir} -maxdepth 1 -type f | $WC -l`
     if [ "<${emptydir}>" != "<0>" ]; then
       iput -v -K -r -f ${curDir}/${dirName} /baoradio/data/AmasNancay/${localpath}
     else
-      $ECHO "DEBUG (${scriptName}): directory ${curDir}/${dirName} is empty"
+      $ECHO "DEBUG (${scriptName}): directory ${outDir} is empty"
     fi
 #