Changeset 577


Ignore:
Timestamp:
Oct 19, 2011, 2:46:27 PM (13 years ago)
Author:
campagne
Message:

improve analysis (jec)

Location:
BAORadio/AmasNancay/trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • BAORadio/AmasNancay/trunk/etude_mergeAna.pic

    r563 r577  
     1defscript evolAtCalibFreq
    12
    2 set source  $1
    3 set ncycles $2
    43
    5 openppf onoffsurvey_${source}-${ncycles}Cycles.ppf
    6 print onoffevol
     4#Follow up of the intensity arround the Calib frequency in different calibration conditions
     5# NO calibration
     6# Calibration coeff. mean per Run
     7# Calibration coeff. per cycles
    78
    8 set defatt "font=helvetica,bold,20 fixedfonsize"
    9 set axedefatt  "font=helvetica,bold,30 fixedfonsize grid"
    10 graphicatt "$defatt"
    11 setaxesatt "$axedefatt"
     9set fcalib $1
     10
    1211
    1312newwin 1 2
    … …  
    1514settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
    1615n/pl onoffevol.onoffRaw1%cycle ! ! "red cpts notit nsta"
    17 setaxelabels "cycle" "I[1346MHz/6.25MHz] (a.u)" $axedefatt
     16setaxelabels "cycle" "I[${fcalib}MHz/6.25MHz] Raw (a.u)" $axedefatt
    1817
    1918newwin 1 2
    … …  
    2120settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
    2221n/proj hRaw1 onoffevol.onoffRaw1 ! ! ! "red cpts notit"
    23 
     22settitle "I[${fcalib}MHz/6.25MHz] Raw (a.u)"
    2423
    2524
    … …  
    2928settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
    3029n/pl onoffevol.onoffRun1%cycle ! ! "red cpts notit nsta"
    31 setaxelabels "cycle" "I[1346MHz/6.25MHz] (a.u)" $axedefatt
     30setaxelabels "cycle" "I[${fcalib}MHz/6.25MHz] Run (a.u)" $axedefatt
     31
    3232
    3333newwin 1 2
    … …  
    3535settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
    3636n/proj hRun1 onoffevol.onoffRun1 ! ! ! "red cpts notit"
     37settitle "I[${fcalib}MHz/6.25MHz] Run (a.u)"
    3738
    3839
    … …  
    4142settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
    4243n/pl onoffevol.onoffCycle1%cycle ! ! "red cpts notit nsta"
    43 setaxelabels "cycle" "I[1346MHz/6.25MHz] (a.u)" $axedefatt
     44setaxelabels "cycle" "I[${fcalib}MHz/6.25MHz] Cycle (a.u)" $axedefatt
    4445
    4546newwin 1 2
    … …  
    4748settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
    4849n/proj hCycle1 onoffevol.onoffCycle1 ! ! ! "red cpts notit"
     50settitle "I[${fcalib}MHz/6.25MHz] Cycle (a.u)"
     51
     52
     53endscript
     54###########################################
     55defscript evolAt1420Freq
     56
     57
     58#Follow up of the intensity arround the (1420.4 +/- 2) MHz frequency in different calibration conditions
     59# NO calibration
     60# Calibration coeff. mean per Run
     61# Calibration coeff. per cycles
     62
     63
     64graphicatt "xylimits=0,500,-0.01,0.01"
     65
     66newwin 1 2
     67n/pl onoffevol.onoffRaw01420%cycle ! ! "blue cpts notit nsta"
     68settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
     69n/pl onoffevol.onoffRaw11420%cycle ! ! "red cpts notit nsta"
     70setaxelabels "cycle" "I[1420.4MHz +/- 2MHz] Raw (a.u)" $axedefatt
     71
     72
     73newwin 1 2
     74n/pl onoffevol.onoffRun01420%cycle ! ! "blue cpts notit nsta"
     75settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
     76n/pl onoffevol.onoffRun11420%cycle ! ! "red cpts notit nsta"
     77setaxelabels "cycle" "I[1420.4MHz +/- 2MHz] Run (a.u)" $axedefatt
     78
     79
     80newwin 1 2
     81n/pl onoffevol.onoffCycle01420%cycle ! ! "blue cpts notit nsta"
     82settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
     83n/pl onoffevol.onoffCycle11420%cycle ! ! "red cpts notit nsta"
     84setaxelabels "cycle" "I[1420.4MHz +/- 2MHz] Cycle (a.u)" $axedefatt
     85
     86
     87newh1d hSigmaRaw0 -0.02 0.02 50
     88newh1d hSigmaRaw1 -0.02 0.02 50
     89newh1d hSigmaRun0 -0.02 0.02 50
     90newh1d hSigmaRun1 -0.02 0.02 50
     91newh1d hSigmaCycl0 -0.02 0.02 50
     92newh1d hSigmaCycl1 -0.02 0.02 50
     93
     94graphicatt ""
     95
     96newwin 1 1
     97n/proj hSigmaRaw0 onoffevol.onoffRaw01420 ! ! ! "blue cpts notit"
     98n/proj hSigmaRaw1 onoffevol.onoffRaw11420 ! ! ! "same red cpts notit"
     99settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
     100setaxelabels "(a.u)" "I[1420.4MHz +/- 2MHz] Raw" $axedefatt
     101
     102newwin 1 1
     103n/proj hSigmaRun0 onoffevol.onoffRun01420 ! ! ! "blue cpts notit"
     104n/proj hSigmaRun1 onoffevol.onoffRun11420 ! ! ! "same red cpts notit"
     105settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
     106setaxelabels "(a.u)" "I[1420.4MHz +/- 2MHz] Run" $axedefatt
     107
     108
     109newwin 1 1
     110n/proj hSigmaCycl0 onoffevol.onoffCycle01420 ! ! ! "blue cpts notit"
     111n/proj hSigmaCycl1 onoffevol.onoffCycle11420 ! ! ! "same red cpts notit"
     112settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
     113setaxelabels "(a.u)" "I[1420.4MHz +/- 2MHz] Cycle" $axedefatt
     114
     115
     116endscript
     117##################################################
     118defscript diffOnOff
     119
     120#
     121#Display the ON-OFF over the whole freq. range in diffrent Calib. conditions
     122# NO calibration
     123# Calibration coeff. mean per Run
     124# Calibration coeff. per cycles
     125#
    49126
    50127
    … …  
    59136settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
    60137setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
    61 
    62 
    63138
    64139del sruncal0 sruncal1
    … …  
    79154settitle "ON-OFF Calib Cycle  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
    80155setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
     156
     157
     158endscript
     159
     160
     161
     162
     163#Main
     164#use minuscule for source label
     165
     166set source  $1
     167set ncycles $2
     168set fcalib $3
     169
     170openppf onoffsurvey_${source}-${ncycles}Cycles.ppf
     171print onoffevol
     172
     173set defatt "font=helvetica,bold,20 fixedfonsize"
     174set axedefatt  "font=helvetica,bold,30 fixedfonsize grid"
     175
     176set defatt ""
     177set axedefatt  "grid"
     178
     179graphicatt ""
     180setaxesatt ""
     181
     182
     183
     184#diffOnOff
     185
     186
     187evolAt1420Freq

  • BAORadio/AmasNancay/trunk/mergeAnaFiles.cc

    r563 r577  
    3838//Usage
    3939//
    40 // ./Objs/mergeAnaFiles -act meanCalibBAODiffOnOff -inPath /sps/baoradio/AmasNancay/JEC/ -src Abell85 -ppfFile dataRaw -debug 1 -mxcycles 500
    41 //
     40//./Objs/mergeAnaFiles -act meanCalibBAODiffOnOff -inPath /sps/baoradio/AmasNancay/JEC/ -src Abell85 -ppfFile dataRaw -debug 1 -mxcycles 500
     41//./Objs/mergeAnaFiles -src Abell85 -act meanCalibBAODiffOnOff -mxcycles 500 -calibfreq 1410 -inPath /sps/baoradio/AmasNancay/JEC/ -ppfFile dataRaw -debug 1 >& mergeAna-500.log
    4242//
    4343//
    … …  
    267267  TMatrix<r_4> meanDiffONOFF_perRunCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);      //set to 0
    268268  TMatrix<r_4> meanDiffONOFF_perCycleCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);    //set to 0
    269   char* onffTupleName[7]={"cycle",
    270                           "onoffRaw0","onoffRaw1",
    271                           "onoffRun0","onoffRun1",
    272                           "onoffCycle0","onoffCycle1"};
    273   NTuple onoffevolution(7,onffTupleName);
    274   r_4 xnt[7];
    275 
    276   //Lower and Higher freq. bin to perform mean follow up
    277   sa_size_t chLow  = freqToChan(para.rcalibFreq_ - (para.rcalibBandFreq_*0.5));
    278   sa_size_t chHigh = freqToChan(para.rcalibFreq_ + (para.rcalibBandFreq_*0.5));
     269  static const int NINFO=19;
     270  char* onffTupleName[NINFO]={"cycle"
     271                          ,"onoffRaw0","onoffRaw1"
     272                          ,"onoffRun0","onoffRun1"
     273                          ,"onoffCycle0","onoffCycle1"
     274                          ,"onoffRaw01420","onoffRaw11420"
     275                          ,"onoffRun01420","onoffRun11420"
     276                          ,"onoffCycle01420","onoffCycle11420"
     277                          ,"onoffRaw01420side","onoffRaw11420side"
     278                          ,"onoffRun01420side","onoffRun11420side"
     279                          ,"onoffCycle01420side","onoffCycle11420side"
     280  };
     281  NTuple onoffevolution(NINFO,onffTupleName);
     282  r_4 xnt[NINFO];
     283
     284  //Lower and Higher freq. arround the Calibration Freq. bin to perform mean follow up
     285  sa_size_t chCalibLow  = freqToChan(para.rcalibFreq_ - (para.rcalibBandFreq_*0.5));
     286  sa_size_t chCalibHigh = freqToChan(para.rcalibFreq_ + (para.rcalibBandFreq_*0.5));
     287  //Lower and Higher freq.  just arround 1420.4MHz Freq. bin to perform mean follow up
     288  sa_size_t ch1420Low  = freqToChan(1420);
     289  sa_size_t ch1420High = freqToChan(1421);
     290
     291  //Lower and Higher freq. on the sides of 1420.4Mhz Freq. bin to perform mean follow up
     292  sa_size_t ch1420aLow  = freqToChan(1418);
     293  sa_size_t ch1420aHigh = freqToChan(1419);
     294  sa_size_t ch1420bLow  = freqToChan(1422);
     295  sa_size_t ch1420bHigh = freqToChan(1423);
     296
     297
    279298  if (para.debuglev_>90){
    280     cout << "freq. band for follow up [" <<  chLow << ", "<< chHigh << "]" << endl;
     299    cout << "freq. band for follow up [" <<  chCalibLow << ", "<< chCalibHigh << "]" << endl;
     300    cout << "freq. band for follow up [" <<  ch1420Low << ", "<< ch1420High << "]" << endl;
    281301  }
    282302 
    … …  
    619639      //Fill NTuple
    620640      xnt[0] = totalNumberCycles;
     641 
     642      //Follow up arround the Calibration Freq.
     643      TVector<r_4> meanInRange_CalibFreq_noCalib(NUMBER_OF_CHANNELS);
     644      meanInRange(diffOnOff_noCalib,chCalibLow,chCalibHigh,meanInRange_CalibFreq_noCalib);
     645      xnt[1] = meanInRange_CalibFreq_noCalib(0);
     646      xnt[2] = meanInRange_CalibFreq_noCalib(1);
     647
     648      TVector<r_4> meanInRange_CalibFreq_perRunCalib(NUMBER_OF_CHANNELS);
     649      meanInRange(diffOnOff_perRunCalib,chCalibLow,chCalibHigh,meanInRange_CalibFreq_perRunCalib);
     650      xnt[3] = meanInRange_CalibFreq_perRunCalib(0);
     651      xnt[4] = meanInRange_CalibFreq_perRunCalib(1);
     652
     653      TVector<r_4> meanInRange_CalibFreq_perCycleCalib(NUMBER_OF_CHANNELS);
     654      meanInRange(diffOnOff_perCycleCalib,chCalibLow,chCalibHigh,meanInRange_CalibFreq_perCycleCalib);
     655      xnt[5] = meanInRange_CalibFreq_perCycleCalib(0);
     656      xnt[6] = meanInRange_CalibFreq_perCycleCalib(1);
     657
     658
     659      //Follow up arround the 1420.4MHz Freq.
     660      TVector<r_4> meanInRange_1420Freq_noCalib(NUMBER_OF_CHANNELS);
     661      meanInRange(diffOnOff_noCalib,ch1420Low,ch1420High,meanInRange_1420Freq_noCalib);
     662      xnt[7] = meanInRange_1420Freq_noCalib(0);
     663      xnt[8] = meanInRange_1420Freq_noCalib(1);
     664
     665      TVector<r_4> meanInRange_1420Freq_perRunCalib(NUMBER_OF_CHANNELS);
     666      meanInRange(diffOnOff_perRunCalib,ch1420Low,ch1420High,meanInRange_1420Freq_perRunCalib);
     667      xnt[9] = meanInRange_1420Freq_perRunCalib(0);
     668      xnt[10] = meanInRange_1420Freq_perRunCalib(1);
     669
     670      TVector<r_4> meanInRange_1420Freq_perCycleCalib(NUMBER_OF_CHANNELS);
     671      meanInRange(diffOnOff_perCycleCalib,ch1420Low,ch1420High,meanInRange_1420Freq_perCycleCalib);
     672      xnt[11] = meanInRange_1420Freq_perCycleCalib(0);
     673      xnt[12] = meanInRange_1420Freq_perCycleCalib(1);
     674
     675
     676      //Follow up below the 1420.4MHz Freq.
     677      TVector<r_4> meanInRange_1420aFreq_noCalib(NUMBER_OF_CHANNELS);
     678      meanInRange(diffOnOff_noCalib,ch1420aLow,ch1420aHigh,meanInRange_1420aFreq_noCalib);
     679      TVector<r_4> meanInRange_1420bFreq_noCalib(NUMBER_OF_CHANNELS);
     680      meanInRange(diffOnOff_noCalib,ch1420bLow,ch1420bHigh,meanInRange_1420bFreq_noCalib);
     681
     682      xnt[13] = (meanInRange_1420aFreq_noCalib(0) + meanInRange_1420bFreq_noCalib(0))/2.;
     683      xnt[14] = (meanInRange_1420aFreq_noCalib(1) + meanInRange_1420bFreq_noCalib(1))/2.;
     684
     685
     686      TVector<r_4> meanInRange_1420aFreq_perRun(NUMBER_OF_CHANNELS);
     687      meanInRange(diffOnOff_perRunCalib,ch1420aLow,ch1420aHigh,meanInRange_1420aFreq_perRun);
     688      TVector<r_4> meanInRange_1420bFreq_perRun(NUMBER_OF_CHANNELS);
     689      meanInRange(diffOnOff_perRunCalib,ch1420bLow,ch1420bHigh,meanInRange_1420bFreq_perRun);
     690
     691      xnt[15] = (meanInRange_1420aFreq_perRun(0) + meanInRange_1420bFreq_perRun(0))/2.;
     692      xnt[16] = (meanInRange_1420aFreq_perRun(1) + meanInRange_1420bFreq_perRun(1))/2.;
     693
     694
     695      TVector<r_4> meanInRange_1420aFreq_perCycle(NUMBER_OF_CHANNELS);
     696      meanInRange(diffOnOff_perCycleCalib,ch1420aLow,ch1420aHigh,meanInRange_1420aFreq_perCycle);
     697      TVector<r_4> meanInRange_1420bFreq_perCycle(NUMBER_OF_CHANNELS);
     698      meanInRange(diffOnOff_perCycleCalib,ch1420bLow,ch1420bHigh,meanInRange_1420bFreq_perCycle);
     699
     700      xnt[17] = (meanInRange_1420aFreq_perCycle(0) + meanInRange_1420bFreq_perCycle(0))/2.;
     701      xnt[18] = (meanInRange_1420aFreq_perCycle(1) + meanInRange_1420bFreq_perCycle(1))/2.;
     702
    621703     
    622       TVector<r_4> meanInRange_noCalib(NUMBER_OF_CHANNELS);
    623       meanInRange(diffOnOff_noCalib,chLow,chHigh,meanInRange_noCalib);
    624       xnt[1] = meanInRange_noCalib(0);
    625       xnt[2] = meanInRange_noCalib(1);
    626 
    627       TVector<r_4> meanInRange_perRunCalib(NUMBER_OF_CHANNELS);
    628       meanInRange(diffOnOff_perRunCalib,chLow,chHigh,meanInRange_perRunCalib);
    629       xnt[3] = meanInRange_perRunCalib(0);
    630       xnt[4] = meanInRange_perRunCalib(1);
    631 
    632       TVector<r_4> meanInRange_perCycleCalib(NUMBER_OF_CHANNELS);
    633       meanInRange(diffOnOff_perCycleCalib,chLow,chHigh,meanInRange_perCycleCalib);
    634       xnt[5] = meanInRange_perCycleCalib(0);
    635       xnt[6] = meanInRange_perCycleCalib(1);
    636      
     704      //store infos to Ntuple
    637705      onoffevolution.Fill(xnt);
    638706
    … …  
    9521020           << " -calibfreq <number> (cf. freq. used by calibration operation)\n"
    9531021           << " -calibbandfreq <number> (cf. band of freq. used by calibration operation)\n"
    954            << " -src [Abell85]\n -inPath [.]|<top_root_dir of the ppf file>\n"
     1022           << " -src [Abell85]\n"
     1023           << " -inPath [.]|<top_root_dir of the ppf file>\n"
    9551024           << " (ex. /sps/baoradio/AmasNancay/JEC/\n "
    9561025           << " -outPath [.]|<dir of the output> \n"
Note: See TracChangeset for help on using the changeset viewer.