Changeset 611 for BAORadio/AmasNancay/trunk/mergeAnaFiles.cc

Timestamp:

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

Author:

torrento

Message:

Add pic analysis files and emptydir function for proc_specmfib

File:

• BAORadio/AmasNancay/trunk/mergeAnaFiles.cc (modified) (15 diffs)

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