Changeset 607 for BAORadio/AmasNancay/trunk
- Timestamp:
- Nov 21, 2011, 3:35:31 PM (13 years ago)
- Location:
- BAORadio/AmasNancay/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/etude_mergeAna.pic
r606 r607 25 25 clearscript sigmaRaw1420Side 26 26 clearscript sigmaRaw1400a1420 27 clearscript evolAt14 00a1420OffOn27 clearscript evolAt1410a1415OffOn 28 28 ############################################################################ 29 29 defscript evolAtCalibFreq … … 298 298 endscript 299 299 ########################################### 300 defscript evolAt14 00a1420OffOn301 302 graphicatt "xylimits=0, 500,0.99,1.02"303 304 newwin 1 1 305 n/pl onoffevol.offRaw0f14 001420%cycle ! ! "blue cpts notit nsta"306 n/pl onoffevol.offRaw1f14 001420%cycle ! ! "same red cpts notit nsta"307 n/pl onoffevol.onRaw0f14 001420%cycle ! ! "same cyan cpts notit nsta"308 n/pl onoffevol.onRaw1f14 001420%cycle ! ! "same orange cpts notit nsta"300 defscript evolAt1410a1415OffOn 301 302 graphicatt "xylimits=0,${ncycles},0.99,1.02" 303 304 newwin 1 1 305 n/pl onoffevol.offRaw0f14101415%cycle ! ! "blue cpts notit nsta" 306 n/pl onoffevol.offRaw1f14101415%cycle ! ! "same red cpts notit nsta" 307 n/pl onoffevol.onRaw0f14101415%cycle ! ! "same cyan cpts notit nsta" 308 n/pl onoffevol.onRaw1f14101415%cycle ! ! "same orange cpts notit nsta" 309 309 settitle "Raw (OFF,ON)/OFF ${source} Ch 0 (blue/cyan) Ch 1 (red/orange) ${ncycles}cycles " ' ' $defatt 310 setaxelabels "cycle" "I[14 00,1420]MHz Raw (a.u)" $axedefatt310 setaxelabels "cycle" "I[1410,1415]MHz Raw (a.u)" $axedefatt 311 311 312 graphicatt "xylimits=0,${ncycles},-0.005,0.005" 313 newwin 1 1 314 n/pl onoffevol.(onRaw0f14101415-offRaw0f14101415)%cycle ! ! "blue cpts notit nsta" 315 n/pl onoffevol.(onRaw1f14101415-offRaw1f14101415)%cycle ! ! "same red cpts notit nsta" 316 settitle "Raw (OFF-ON)/OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 317 setaxelabels "cycle" "DI [1410,1415]MHz Raw (a.u)" $axedefatt 318 319 312 320 endscript 313 321 ########################################### -
BAORadio/AmasNancay/trunk/mergeAnaFiles.cc
r606 r607 304 304 ,"onoffRun01420side","onoffRun11420side" 305 305 ,"onoffCycle01420side","onoffCycle11420side" 306 ,"onoffRaw0f14 001420","onoffRaw1f14001420"307 ,"offRaw0f14 001420","offRaw1f14001420"308 ,"onRaw0f14 001420","onRaw1f14001420"306 ,"onoffRaw0f14101415","onoffRaw1f14101415" 307 ,"offRaw0f14101415","offRaw1f14101415" 308 ,"onRaw0f14101415","onRaw1f14101415" 309 309 }; 310 310 NTuple onoffevolution(NINFO,onffTupleName); … … 327 327 sa_size_t ch1400 = freqToChan(1400); 328 328 sa_size_t ch1405 = freqToChan(1400); 329 sa_size_t ch1410 = freqToChan(1410); 329 330 sa_size_t ch1415 = freqToChan(1415); 330 331 sa_size_t ch1420 = freqToChan(1420); … … 758 759 aSpecOffovOff.Div(aSpecOffFitltered); 759 760 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); 761 TVector<r_4> meanInRange_1410a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS); 762 meanInRange(aSpecOffovOff,ch1410,ch1415,meanInRange_1410a1415Freq_OFF_noCalib); 763 764 xnt[21] = meanInRange_1410a1415Freq_OFF_noCalib(0); 765 xnt[22] = meanInRange_1410a1415Freq_OFF_noCalib(1); 766 768 767 TMatrix<r_4> aSpecOnovOff(aSpecOn,false); 769 768 aSpecOnovOff.Div(aSpecOffFitltered); 770 769 771 TVector<r_4> meanInRange_14 00a1420Freq_ON_noCalib(NUMBER_OF_CHANNELS);772 meanInRange(aSpecOnovOff,ch14 00,ch1420,meanInRange_1400a1420Freq_ON_noCalib);773 774 xnt[23] = meanInRange_14 00a1420Freq_ON_noCalib(0);775 xnt[24] = meanInRange_14 00a1420Freq_ON_noCalib(1);770 TVector<r_4> meanInRange_1410a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS); 771 meanInRange(aSpecOnovOff,ch1410,ch1415,meanInRange_1410a1415Freq_ON_noCalib); 772 773 xnt[23] = meanInRange_1410a1415Freq_ON_noCalib(0); 774 xnt[24] = meanInRange_1410a1415Freq_ON_noCalib(1); 776 775 777 776 //store infos to Ntuple … … 852 851 //Used like: 853 852 // 854 void meanRawDiffOnOffCycles() throw(string) { 853 void 854 meanRawDiffOnOffCycles() throw(string) { 855 855 list<string> listOfFiles; 856 856 string directoryName;
Note: See TracChangeset
for help on using the changeset viewer.