Changeset 608 for BAORadio

Timestamp:

Nov 22, 2011, 11:33:35 AM (13 years ago)

Author:

campagne

Message:

setup pure RAW data analysis, and uniformaize the use of (ON-OFF)/OFF (jec)

Location:

BAORadio/AmasNancay/trunk

Files:

• mergeAnaFiles.cc (modified) (4 diffs)
• plotRawDiffOnOff.pic (added)

BAORadio/AmasNancay/trunk/mergeAnaFiles.cc

@@ -326 +326 @@
   //1400-1420Mhz
   sa_size_t ch1400 = freqToChan(1400);
-  sa_size_t ch1405 = freqToChan(1400);
+  //  sa_size_t ch1405 = freqToChan(1400);
   sa_size_t ch1410 = freqToChan(1410);
   sa_size_t ch1415 = freqToChan(1415);
@@ -869 +869 @@
   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;
+  char* onffTupleName[NINFO]={"cycle"
+                              ,"onoffRaw01420","onoffRaw11420"
+                              ,"onoffRaw01420side","onoffRaw11420side"
+                              ,"onoffRaw0f14101415","onoffRaw1f14101415"
+                              ,"offRaw0f14101415","offRaw1f14101415"
+                              ,"onRaw0f14101415","onRaw1f14101415"
+  };
+  NTuple onoffevolution(NINFO,onffTupleName);
+  r_4 xnt[NINFO];
+
+  //Lower and Higher freq.  just arround 1420.4MHz Freq. bin to perform mean follow up
+  sa_size_t ch1420Low  = freqToChan(1420.4-0.2);
+  sa_size_t ch1420High = freqToChan(1420.4+0.2);
+
+  //Lower and Higher freq. on the sides of 1420.4Mhz Freq. bin to perform mean follow up
+  sa_size_t ch1420aLow  = freqToChan(1418);
+  sa_size_t ch1420aHigh = freqToChan(1419);
+  sa_size_t ch1420bLow  = freqToChan(1422);
+  sa_size_t ch1420bHigh = freqToChan(1423);
+  
+  //1400-1420Mhz
+  sa_size_t ch1400 = freqToChan(1400);
+  //  sa_size_t ch1405 = freqToChan(1400);
+  sa_size_t ch1410 = freqToChan(1410);
+  sa_size_t ch1415 = freqToChan(1415);
+  sa_size_t ch1420 = freqToChan(1420);
+
+  if (para.debuglev_>0){
+    cout << "freq. band for follow up [" <<  ch1420Low << ", "<< ch1420High << "]" << endl;
+    cout << "freq. band for follow up [" <<  ch1420aLow << ", "<< ch1420aHigh << "]" << endl;
+    cout << "freq. band for follow up [" <<  ch1420bLow << ", "<< ch1420bHigh << "]" << endl;
+  }
+
   int totalNumberCycles=0; //total number of cycles
 
@@ -970 +1005 @@
       totalNumberCycles++;
 
+
+      //Fill NTuple
+      xnt[0] = totalNumberCycles;
+ 
+
+      //Follow up arround the 1420.4MHz Freq.
+      TVector<r_4> meanInRange_1420Freq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOffovOff_noCalib,ch1420Low,ch1420High,meanInRange_1420Freq_noCalib);
+      xnt[1] = meanInRange_1420Freq_noCalib(0);
+      xnt[2] = meanInRange_1420Freq_noCalib(1);
+
+
+      //Follow up below the 1420.4MHz Freq.
+      TVector<r_4> meanInRange_1420aFreq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOffovOff_noCalib,ch1420aLow,ch1420aHigh,meanInRange_1420aFreq_noCalib);
+      TVector<r_4> meanInRange_1420bFreq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOffovOff_noCalib,ch1420bLow,ch1420bHigh,meanInRange_1420bFreq_noCalib);
+
+      xnt[3] = (meanInRange_1420aFreq_noCalib(0) + meanInRange_1420bFreq_noCalib(0))/2.;
+      xnt[4] = (meanInRange_1420aFreq_noCalib(1) + meanInRange_1420bFreq_noCalib(1))/2.;
+
+
+      //JEC 25/10/11 follow 1400-1420 MHz bande protege et n'inclue pas le 1420.4Mhz de la Galaxie
+      TVector<r_4> meanInRange_1400a1420Freq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOffovOff_noCalib,ch1400,ch1420,meanInRange_1400a1420Freq_noCalib);
+      xnt[5] = meanInRange_1400a1420Freq_noCalib(0);
+      xnt[6] = meanInRange_1400a1420Freq_noCalib(1);
+
+      //JEC 18/11/11 follow up the 1400-1420MHz OFF only
+      TMatrix<r_4> aSpecOffovOff(aSpecOff,false);
+      aSpecOffovOff.Div(aSpecOffFitltered);
+
+      TVector<r_4> meanInRange_1410a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(aSpecOffovOff,ch1410,ch1415,meanInRange_1410a1415Freq_OFF_noCalib);
+
+      xnt[7] = meanInRange_1410a1415Freq_OFF_noCalib(0);
+      xnt[8] = meanInRange_1410a1415Freq_OFF_noCalib(1);
+
+      TMatrix<r_4> aSpecOnovOff(aSpecOn,false);
+      aSpecOnovOff.Div(aSpecOffFitltered);
+
+      TVector<r_4> meanInRange_1410a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(aSpecOnovOff,ch1410,ch1415,meanInRange_1410a1415Freq_ON_noCalib);
+
+      xnt[9] = meanInRange_1410a1415Freq_ON_noCalib(0);
+      xnt[10] = meanInRange_1410a1415Freq_ON_noCalib(1);
+      
+      //store infos to Ntuple
+      onoffevolution.Fill(xnt);
+
       //Quit if enough
       if (totalNumberCycles >= para.maxNumberCycles_) break;    
@@ -1005 +1090 @@
     tag = "meanOffovOffNoCalib";
     fos << PPFNameTag(tag) << meanOFFovOFF_noCalib;
+
+    tag = "onoffevol";
+    fos << PPFNameTag(tag) << onoffevolution;    
 
   }//end save