Changeset 577 for BAORadio/AmasNancay/trunk
- Timestamp:
- Oct 19, 2011, 2:46:27 PM (13 years ago)
- Location:
- BAORadio/AmasNancay/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/etude_mergeAna.pic
r563 r577 1 defscript evolAtCalibFreq 1 2 2 set source $13 set ncycles $24 3 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 7 8 8 set defatt "font=helvetica,bold,20 fixedfonsize" 9 set axedefatt "font=helvetica,bold,30 fixedfonsize grid" 10 graphicatt "$defatt" 11 setaxesatt "$axedefatt" 9 set fcalib $1 10 12 11 13 12 newwin 1 2 … … 15 14 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 16 15 n/pl onoffevol.onoffRaw1%cycle ! ! "red cpts notit nsta" 17 setaxelabels "cycle" "I[ 1346MHz/6.25MHz](a.u)" $axedefatt16 setaxelabels "cycle" "I[${fcalib}MHz/6.25MHz] Raw (a.u)" $axedefatt 18 17 19 18 newwin 1 2 … … 21 20 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 22 21 n/proj hRaw1 onoffevol.onoffRaw1 ! ! ! "red cpts notit" 23 22 settitle "I[${fcalib}MHz/6.25MHz] Raw (a.u)" 24 23 25 24 … … 29 28 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 30 29 n/pl onoffevol.onoffRun1%cycle ! ! "red cpts notit nsta" 31 setaxelabels "cycle" "I[1346MHz/6.25MHz] (a.u)" $axedefatt 30 setaxelabels "cycle" "I[${fcalib}MHz/6.25MHz] Run (a.u)" $axedefatt 31 32 32 33 33 newwin 1 2 … … 35 35 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 36 36 n/proj hRun1 onoffevol.onoffRun1 ! ! ! "red cpts notit" 37 settitle "I[${fcalib}MHz/6.25MHz] Run (a.u)" 37 38 38 39 … … 41 42 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 42 43 n/pl onoffevol.onoffCycle1%cycle ! ! "red cpts notit nsta" 43 setaxelabels "cycle" "I[ 1346MHz/6.25MHz](a.u)" $axedefatt44 setaxelabels "cycle" "I[${fcalib}MHz/6.25MHz] Cycle (a.u)" $axedefatt 44 45 45 46 newwin 1 2 … … 47 48 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 48 49 n/proj hCycle1 onoffevol.onoffCycle1 ! ! ! "red cpts notit" 50 settitle "I[${fcalib}MHz/6.25MHz] Cycle (a.u)" 51 52 53 endscript 54 ########################################### 55 defscript 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 64 graphicatt "xylimits=0,500,-0.01,0.01" 65 66 newwin 1 2 67 n/pl onoffevol.onoffRaw01420%cycle ! ! "blue cpts notit nsta" 68 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 69 n/pl onoffevol.onoffRaw11420%cycle ! ! "red cpts notit nsta" 70 setaxelabels "cycle" "I[1420.4MHz +/- 2MHz] Raw (a.u)" $axedefatt 71 72 73 newwin 1 2 74 n/pl onoffevol.onoffRun01420%cycle ! ! "blue cpts notit nsta" 75 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 76 n/pl onoffevol.onoffRun11420%cycle ! ! "red cpts notit nsta" 77 setaxelabels "cycle" "I[1420.4MHz +/- 2MHz] Run (a.u)" $axedefatt 78 79 80 newwin 1 2 81 n/pl onoffevol.onoffCycle01420%cycle ! ! "blue cpts notit nsta" 82 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 83 n/pl onoffevol.onoffCycle11420%cycle ! ! "red cpts notit nsta" 84 setaxelabels "cycle" "I[1420.4MHz +/- 2MHz] Cycle (a.u)" $axedefatt 85 86 87 newh1d hSigmaRaw0 -0.02 0.02 50 88 newh1d hSigmaRaw1 -0.02 0.02 50 89 newh1d hSigmaRun0 -0.02 0.02 50 90 newh1d hSigmaRun1 -0.02 0.02 50 91 newh1d hSigmaCycl0 -0.02 0.02 50 92 newh1d hSigmaCycl1 -0.02 0.02 50 93 94 graphicatt "" 95 96 newwin 1 1 97 n/proj hSigmaRaw0 onoffevol.onoffRaw01420 ! ! ! "blue cpts notit" 98 n/proj hSigmaRaw1 onoffevol.onoffRaw11420 ! ! ! "same red cpts notit" 99 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 100 setaxelabels "(a.u)" "I[1420.4MHz +/- 2MHz] Raw" $axedefatt 101 102 newwin 1 1 103 n/proj hSigmaRun0 onoffevol.onoffRun01420 ! ! ! "blue cpts notit" 104 n/proj hSigmaRun1 onoffevol.onoffRun11420 ! ! ! "same red cpts notit" 105 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 106 setaxelabels "(a.u)" "I[1420.4MHz +/- 2MHz] Run" $axedefatt 107 108 109 newwin 1 1 110 n/proj hSigmaCycl0 onoffevol.onoffCycle01420 ! ! ! "blue cpts notit" 111 n/proj hSigmaCycl1 onoffevol.onoffCycle11420 ! ! ! "same red cpts notit" 112 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 113 setaxelabels "(a.u)" "I[1420.4MHz +/- 2MHz] Cycle" $axedefatt 114 115 116 endscript 117 ################################################## 118 defscript 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 # 49 126 50 127 … … 59 136 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 60 137 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 61 62 63 138 64 139 del sruncal0 sruncal1 … … 79 154 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 80 155 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 156 157 158 endscript 159 160 161 162 163 #Main 164 #use minuscule for source label 165 166 set source $1 167 set ncycles $2 168 set fcalib $3 169 170 openppf onoffsurvey_${source}-${ncycles}Cycles.ppf 171 print onoffevol 172 173 set defatt "font=helvetica,bold,20 fixedfonsize" 174 set axedefatt "font=helvetica,bold,30 fixedfonsize grid" 175 176 set defatt "" 177 set axedefatt "grid" 178 179 graphicatt "" 180 setaxesatt "" 181 182 183 184 #diffOnOff 185 186 187 evolAt1420Freq -
BAORadio/AmasNancay/trunk/mergeAnaFiles.cc
r563 r577 38 38 //Usage 39 39 // 40 // 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 42 42 // 43 43 // … … 267 267 TMatrix<r_4> meanDiffONOFF_perRunCalib(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); //set to 0 268 268 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 279 298 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; 281 301 } 282 302 … … 619 639 //Fill NTuple 620 640 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 621 703 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 637 705 onoffevolution.Fill(xnt); 638 706 … … 952 1020 << " -calibfreq <number> (cf. freq. used by calibration operation)\n" 953 1021 << " -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" 955 1024 << " (ex. /sps/baoradio/AmasNancay/JEC/\n " 956 1025 << " -outPath [.]|<dir of the output> \n"
Note: See TracChangeset
for help on using the changeset viewer.