Changeset 628 for BAORadio/AmasNancay/trunk
- Timestamp:
- Jan 13, 2012, 11:41:42 AM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/etude_mergeAna.pic
r613 r628 59 59 set fcalib $1 60 60 61 62 61 newwin 1 2 63 62 n/pl onoffevol.onoffRaw0%cycle ! ! "blue cpts notit nsta" … … 106 105 107 106 108 #Follow up of the intensity ar round the [1420,1421]MHz frequency in different calibration conditions107 #Follow up of the intensity around the [1420,1421]MHz frequency in different calibration conditions 109 108 # NO calibration 110 109 # Calibration coeff. mean per Run … … 134 133 setaxelabels "cycle" "I[${f1420low},${f1420high}]MHz Cycle (a.u)" $axedefatt 135 134 136 137 newh1d hSigmaRaw0 -0.01 0.01 25 138 newh1d hSigmaRaw1 -0.01 0.01 25 139 newh1d hSigmaRun0 -0.01 0.01 25 140 newh1d hSigmaRun1 -0.01 0.01 25 141 newh1d hSigmaCycl0 -0.01 0.01 25 142 newh1d hSigmaCycl1 -0.01 0.01 25 135 set low -0.02 136 set high 0.02 137 newh1d hSigmaRaw0 ${low} ${high} 25 138 newh1d hSigmaRaw1 ${low} ${high} 25 139 newh1d hSigmaRun0 ${low} ${high} 25 140 newh1d hSigmaRun1 ${low} ${high} 25 141 newh1d hSigmaCycl0 ${low} ${high} 25 142 newh1d hSigmaCycl1 ${low} ${high} 25 143 143 144 144 graphicatt "" … … 236 236 setaxelabels "num of cycles" "Sigma (a.u)" $axedefatt 237 237 238 239 240 endscript 238 endscript 239 241 240 ################################################## 242 241 defscript sigmaRaw1400a1420 243 242 244 c++compile rebining 245 c++link rebining.so dorebin 246 247 248 #col 19 = onoffRaw0f14001420 Ch 0 249 #col 20 = onoffRaw1f14001420 Ch 1 243 ## mergeAnaFiles.cc -> CalibBAODiffOnOffCycles -> onoffsurvey_[source]-[n]Cycles.ppf 244 ##col 19 = onoffRaw0f14001420 Ch 0 245 ##col 20 = onoffRaw1f14001420 Ch 1 250 246 251 247 ntcol2var onoffevol 19 linCh0All 252 248 line2vec vecCh0All $linCh0All 253 254 249 ntcol2var onoffevol 20 linCh1All 255 250 line2vec vecCh1All $linCh1All 256 251 252 graphicatt "" 257 253 set sigmaCh0 "" 258 254 set errsigCh0 "" … … 260 256 set errsigCh1 "" 261 257 258 c++compile rebining 259 c++link rebining.so dorebin 262 260 263 261 foreach ipack ( $packVal ) … … 285 283 end 286 284 287 288 285 newnt nt0 x y ey 289 286 newnt nt1 x y ey … … 296 293 set y1 ( $sigmaCh1 ) 297 294 set ey1 ( $errsigCh1 ) 298 295 299 296 for i 0:$#x0 300 297 line2nt nt0 $x0[i] $y0[i] $ey0[i] 301 298 line2nt nt1 $x1[i] $y1[i] $ey1[i] 302 299 end 303 304 300 305 301 newwin 1 1 … … 311 307 setaxelabels "num of cycles" "Sigma (a.u)" $axedefatt 312 308 313 314 315 endscript 309 newwin 1 1 310 plot2de nt1 x y 0 ey 1 "red marker=fcircle,9 notit nsta" 311 plot2de nt0 x y 0 ey 1 "same blue marker=fcircle,9 notit nsta" 312 n = ${#x0}-1 313 func $y0[0]/sqrt(x) $x0[0] $x0[n] 100 "same" 314 settitle "Sigma [1400,1420]MHz ${source} Ch 0 (blue) Ch 1 (red)" ' ' $defatt 315 setaxelabels "num of cycles" "Sigma (a.u)" $axedefatt 316 317 newwin 1 1 318 n/pl onoffevol.onoffRaw0f14001420 ! ! "cpts blue notit" 319 n/pl onoffevol.onoffRaw1f14001420 ! ! "cpts same red notit" 320 settitle "onoffRaw [1400,1420]MHz ${source} Ch 0 (blue) Ch 1 (red)" ' ' $defatt 321 newwin 1 1 322 n/pl onoffevol.onoffRaw0f14001420%cycle ! ! "cpts blue notit" 323 n/pl onoffevol.onoffRaw1f14001420%cycle ! ! "cpts same red notit" 324 settitle "onoffRaw [1400,1420]MHz ${source} Ch 0 (blue) Ch 1 (red)" ' ' $defatt 325 326 endscript 327 316 328 ########################################### 317 329 defscript evolAt1410a1415OffOn … … 410 422 # Calibration coeff. per cycles 411 423 412 413 424 graphicatt "xylimits=0,${ncycles},-0.1,0.1" 414 425 … … 455 466 setaxelabels "(a.u)" "I 1420MHz Run" $axedefatt 456 467 457 458 459 468 newwin 1 1 460 469 n/proj h1420IntCycl0 onoffevol.${nbins1420}*(onoffCycle01420-onoffCycle01420side) ! ! ! "blue cpts notit" … … 462 471 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 463 472 setaxelabels "(a.u)" "I 1420MHz Cycle" $axedefatt 464 465 466 473 467 474 endscript … … 476 483 # 477 484 478 graphicatt "" 485 # ON/OFF, OFF/OFF 486 479 487 del sonovoff0 sonovoff1 480 488 objaoper meanOnovOffNoCalib row 0 sonovoff0 … … 493 501 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 494 502 495 503 # (ON-OFF)/OFF 496 504 graphicatt "xylimits=1250,1500,-0.01,0.01" 497 505 del sonoffovoff0 sonoffovoff1 … … 504 512 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 505 513 506 graphicatt "xylimits=1250,1500,-0.01,0.01" 507 del snocal0 snocal1 508 objaoper meanNoCalib row 0 snocal0 509 objaoper meanNoCalib row 1 snocal1 510 newwin 1 1 511 plot2d snocal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta" 512 plot2d snocal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta" 513 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 514 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 514 # (ON-OFF) 515 #graphicatt "xylimits=1250,1500,-0.01,0.01" 516 #del snocal0 snocal1 517 #objaoper meanNoCalib row 0 snocal0 518 #objaoper meanNoCalib row 1 snocal1 519 #newwin 1 1 520 #plot2d snocal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta" 521 #plot2d snocal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta" 522 #settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 523 #setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt 515 524 516 525 #del sruncal0 sruncal1 … … 541 550 #use minuscule for source label 542 551 543 openppf ${path}/onoffsurvey_${source}-${ncycles}Cycles.ppf544 print onoffevol545 546 #openppf ${path}/rawOnOffDiff_${source}-${ncycles}Cycles.ppf547 #ls552 #openppf ${path}/onoffsurvey_${source}-${ncycles}Cycles.ppf 553 #print onoffevol 554 555 openppf ${path}/rawOnOffDiff_${source}-${ncycles}Cycles.ppf 556 ls 548 557 549 558 set defatt "font=helvetica,bold,20 fixedfonsize" … … 557 566 558 567 559 #diffOnOff568 diffOnOff 560 569 #evolAt1420Freq 561 570 #evolAt1420SideFreq
Note: See TracChangeset
for help on using the changeset viewer.