Changeset 606 for BAORadio/AmasNancay/trunk/mergeAnaFiles.cc
- Timestamp:
- Nov 21, 2011, 10:47:14 AM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/mergeAnaFiles.cc
r605 r606 305 305 ,"onoffCycle01420side","onoffCycle11420side" 306 306 ,"onoffRaw0f14001420","onoffRaw1f14001420" 307 ,"offRaw0f140014 15","offRaw1f14001415"308 ,"onRaw0f140014 15","onRaw1f14001415"307 ,"offRaw0f14001420","offRaw1f14001420" 308 ,"onRaw0f14001420","onRaw1f14001420" 309 309 }; 310 310 NTuple onoffevolution(NINFO,onffTupleName); … … 326 326 //1400-1420Mhz 327 327 sa_size_t ch1400 = freqToChan(1400); 328 sa_size_t ch1405 = freqToChan(1400); 328 329 sa_size_t ch1415 = freqToChan(1415); 329 330 sa_size_t ch1420 = freqToChan(1420); … … 754 755 755 756 //JEC 18/11/11 follow up the 1400-1420MHz OFF only 756 TVector<r_4> meanInRange_1400a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS); 757 meanInRange(aSpecOff,ch1400,ch1415,meanInRange_1400a1415Freq_OFF_noCalib); 758 759 xnt[21] = meanInRange_1400a1415Freq_OFF_noCalib(0); 760 xnt[22] = meanInRange_1400a1415Freq_OFF_noCalib(1); 761 762 TVector<r_4> meanInRange_1400a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS); 763 meanInRange(aSpecOn,ch1400,ch1415,meanInRange_1400a1415Freq_ON_noCalib); 764 765 xnt[23] = meanInRange_1400a1415Freq_ON_noCalib(0); 766 xnt[24] = meanInRange_1400a1415Freq_ON_noCalib(1); 757 TMatrix<r_4> aSpecOffovOff(aSpecOff,false); 758 aSpecOffovOff.Div(aSpecOffFitltered); 759 760 TVector<r_4> meanInRange_1400a1420Freq_OFF_noCalib(NUMBER_OF_CHANNELS); 761 meanInRange(aSpecOffovOff,ch1400,ch1420,meanInRange_1400a1420Freq_OFF_noCalib); 762 763 xnt[21] = meanInRange_1400a1420Freq_OFF_noCalib(0); 764 xnt[22] = meanInRange_1400a1420Freq_OFF_noCalib(1); 765 766 // TMatrix<r_4> aSpecOnFitltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ); 767 // medianFiltering(aSpecOn,halfWidth,aSpecOnFitltered); 768 TMatrix<r_4> aSpecOnovOff(aSpecOn,false); 769 aSpecOnovOff.Div(aSpecOffFitltered); 770 771 TVector<r_4> meanInRange_1400a1420Freq_ON_noCalib(NUMBER_OF_CHANNELS); 772 meanInRange(aSpecOnovOff,ch1400,ch1420,meanInRange_1400a1420Freq_ON_noCalib); 773 774 xnt[23] = meanInRange_1400a1420Freq_ON_noCalib(0); 775 xnt[24] = meanInRange_1400a1420Freq_ON_noCalib(1); 767 776 768 777 //store infos to Ntuple
Note: See TracChangeset
for help on using the changeset viewer.