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

Timestamp:

Oct 19, 2011, 2:46:27 PM (13 years ago)

Author:

campagne

Message:

improve analysis (jec)

File:

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

BAORadio/AmasNancay/trunk/mergeAnaFiles.cc

@@ -38 +38 @@
 //Usage
 //
-// ./Objs/mergeAnaFiles -act meanCalibBAODiffOnOff -inPath /sps/baoradio/AmasNancay/JEC/ -src Abell85 -ppfFile dataRaw -debug 1 -mxcycles 500
-//
+//./Objs/mergeAnaFiles -act meanCalibBAODiffOnOff -inPath /sps/baoradio/AmasNancay/JEC/ -src Abell85 -ppfFile dataRaw -debug 1 -mxcycles 500
+//./Objs/mergeAnaFiles -src Abell85 -act meanCalibBAODiffOnOff -mxcycles 500 -calibfreq 1410 -inPath /sps/baoradio/AmasNancay/JEC/ -ppfFile dataRaw -debug 1 >& mergeAna-500.log
 //
 //
@@ -267 +267 @@
   TMatrix<r_4> meanDiffONOFF_perRunCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);      //set to 0
   TMatrix<r_4> meanDiffONOFF_perCycleCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);    //set to 0
-  char* onffTupleName[7]={"cycle",
-                          "onoffRaw0","onoffRaw1",
-                          "onoffRun0","onoffRun1",
-                          "onoffCycle0","onoffCycle1"};
-  NTuple onoffevolution(7,onffTupleName);
-  r_4 xnt[7];
-
-  //Lower and Higher freq. bin to perform mean follow up
-  sa_size_t chLow  = freqToChan(para.rcalibFreq_ - (para.rcalibBandFreq_*0.5));
-  sa_size_t chHigh = freqToChan(para.rcalibFreq_ + (para.rcalibBandFreq_*0.5));
+  static const int NINFO=19;
+  char* onffTupleName[NINFO]={"cycle"
+                          ,"onoffRaw0","onoffRaw1"
+                          ,"onoffRun0","onoffRun1"
+                          ,"onoffCycle0","onoffCycle1"
+                          ,"onoffRaw01420","onoffRaw11420"
+                          ,"onoffRun01420","onoffRun11420"
+                          ,"onoffCycle01420","onoffCycle11420"
+                          ,"onoffRaw01420side","onoffRaw11420side"
+                          ,"onoffRun01420side","onoffRun11420side"
+                          ,"onoffCycle01420side","onoffCycle11420side"
+  };
+  NTuple onoffevolution(NINFO,onffTupleName);
+  r_4 xnt[NINFO];
+
+  //Lower and Higher freq. arround the Calibration Freq. bin to perform mean follow up
+  sa_size_t chCalibLow  = freqToChan(para.rcalibFreq_ - (para.rcalibBandFreq_*0.5));
+  sa_size_t chCalibHigh = freqToChan(para.rcalibFreq_ + (para.rcalibBandFreq_*0.5));
+  //Lower and Higher freq.  just arround 1420.4MHz Freq. bin to perform mean follow up
+  sa_size_t ch1420Low  = freqToChan(1420);
+  sa_size_t ch1420High = freqToChan(1421);
+
+  //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);
+
+
   if (para.debuglev_>90){
-    cout << "freq. band for follow up [" <<  chLow << ", "<< chHigh << "]" << endl;
+    cout << "freq. band for follow up [" <<  chCalibLow << ", "<< chCalibHigh << "]" << endl;
+    cout << "freq. band for follow up [" <<  ch1420Low << ", "<< ch1420High << "]" << endl;
   }
   
@@ -619 +639 @@
       //Fill NTuple
       xnt[0] = totalNumberCycles;
+ 
+      //Follow up arround the Calibration Freq.
+      TVector<r_4> meanInRange_CalibFreq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_noCalib,chCalibLow,chCalibHigh,meanInRange_CalibFreq_noCalib);
+      xnt[1] = meanInRange_CalibFreq_noCalib(0);
+      xnt[2] = meanInRange_CalibFreq_noCalib(1);
+
+      TVector<r_4> meanInRange_CalibFreq_perRunCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perRunCalib,chCalibLow,chCalibHigh,meanInRange_CalibFreq_perRunCalib);
+      xnt[3] = meanInRange_CalibFreq_perRunCalib(0);
+      xnt[4] = meanInRange_CalibFreq_perRunCalib(1);
+
+      TVector<r_4> meanInRange_CalibFreq_perCycleCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perCycleCalib,chCalibLow,chCalibHigh,meanInRange_CalibFreq_perCycleCalib);
+      xnt[5] = meanInRange_CalibFreq_perCycleCalib(0);
+      xnt[6] = meanInRange_CalibFreq_perCycleCalib(1);
+
+
+      //Follow up arround the 1420.4MHz Freq.
+      TVector<r_4> meanInRange_1420Freq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_noCalib,ch1420Low,ch1420High,meanInRange_1420Freq_noCalib);
+      xnt[7] = meanInRange_1420Freq_noCalib(0);
+      xnt[8] = meanInRange_1420Freq_noCalib(1);
+
+      TVector<r_4> meanInRange_1420Freq_perRunCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perRunCalib,ch1420Low,ch1420High,meanInRange_1420Freq_perRunCalib);
+      xnt[9] = meanInRange_1420Freq_perRunCalib(0);
+      xnt[10] = meanInRange_1420Freq_perRunCalib(1);
+
+      TVector<r_4> meanInRange_1420Freq_perCycleCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perCycleCalib,ch1420Low,ch1420High,meanInRange_1420Freq_perCycleCalib);
+      xnt[11] = meanInRange_1420Freq_perCycleCalib(0);
+      xnt[12] = meanInRange_1420Freq_perCycleCalib(1);
+
+
+      //Follow up below the 1420.4MHz Freq.
+      TVector<r_4> meanInRange_1420aFreq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_noCalib,ch1420aLow,ch1420aHigh,meanInRange_1420aFreq_noCalib);
+      TVector<r_4> meanInRange_1420bFreq_noCalib(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_noCalib,ch1420bLow,ch1420bHigh,meanInRange_1420bFreq_noCalib);
+
+      xnt[13] = (meanInRange_1420aFreq_noCalib(0) + meanInRange_1420bFreq_noCalib(0))/2.;
+      xnt[14] = (meanInRange_1420aFreq_noCalib(1) + meanInRange_1420bFreq_noCalib(1))/2.;
+
+
+      TVector<r_4> meanInRange_1420aFreq_perRun(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perRunCalib,ch1420aLow,ch1420aHigh,meanInRange_1420aFreq_perRun);
+      TVector<r_4> meanInRange_1420bFreq_perRun(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perRunCalib,ch1420bLow,ch1420bHigh,meanInRange_1420bFreq_perRun);
+
+      xnt[15] = (meanInRange_1420aFreq_perRun(0) + meanInRange_1420bFreq_perRun(0))/2.;
+      xnt[16] = (meanInRange_1420aFreq_perRun(1) + meanInRange_1420bFreq_perRun(1))/2.;
+
+
+      TVector<r_4> meanInRange_1420aFreq_perCycle(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perCycleCalib,ch1420aLow,ch1420aHigh,meanInRange_1420aFreq_perCycle);
+      TVector<r_4> meanInRange_1420bFreq_perCycle(NUMBER_OF_CHANNELS);
+      meanInRange(diffOnOff_perCycleCalib,ch1420bLow,ch1420bHigh,meanInRange_1420bFreq_perCycle);
+
+      xnt[17] = (meanInRange_1420aFreq_perCycle(0) + meanInRange_1420bFreq_perCycle(0))/2.;
+      xnt[18] = (meanInRange_1420aFreq_perCycle(1) + meanInRange_1420bFreq_perCycle(1))/2.;
+
       
-      TVector<r_4> meanInRange_noCalib(NUMBER_OF_CHANNELS);
-      meanInRange(diffOnOff_noCalib,chLow,chHigh,meanInRange_noCalib);
-      xnt[1] = meanInRange_noCalib(0);
-      xnt[2] = meanInRange_noCalib(1);
-
-      TVector<r_4> meanInRange_perRunCalib(NUMBER_OF_CHANNELS);
-      meanInRange(diffOnOff_perRunCalib,chLow,chHigh,meanInRange_perRunCalib);
-      xnt[3] = meanInRange_perRunCalib(0);
-      xnt[4] = meanInRange_perRunCalib(1);
-
-      TVector<r_4> meanInRange_perCycleCalib(NUMBER_OF_CHANNELS);
-      meanInRange(diffOnOff_perCycleCalib,chLow,chHigh,meanInRange_perCycleCalib);
-      xnt[5] = meanInRange_perCycleCalib(0);
-      xnt[6] = meanInRange_perCycleCalib(1);
-      
+      //store infos to Ntuple
       onoffevolution.Fill(xnt);
 
@@ -952 +1020 @@
            << " -calibfreq <number> (cf. freq. used by calibration operation)\n"
            << " -calibbandfreq <number> (cf. band of freq. used by calibration operation)\n"
-           << " -src [Abell85]\n -inPath [.]|<top_root_dir of the ppf file>\n" 
+           << " -src [Abell85]\n" 
+           << " -inPath [.]|<top_root_dir of the ppf file>\n" 
            << " (ex. /sps/baoradio/AmasNancay/JEC/\n " 
            << " -outPath [.]|<dir of the output> \n"