Changeset 588 for BAORadio/AmasNancay

Timestamp:

Oct 24, 2011, 5:10:52 PM (13 years ago)

Author:

campagne

Message:

improve analysis: sigma en ft du nbre de cycle (jec)

Location:

BAORadio/AmasNancay/trunk

Files:

• etude_mergeAna.pic (modified) (3 diffs)
• rebining.cc (added)

BAORadio/AmasNancay/trunk/etude_mergeAna.pic

@@ -27 +27 @@
 clearscript sigmaRaw1420Side
 ############################################################################
-
 defscript evolAtCalibFreq
 
@@ -147 +146 @@
 defscript sigmaRaw1420Side
 
-
-newh1d hSigRaw0All -0.01 0.01 25
-newh1d hSigRaw1All -0.01 0.01 25
-
-n/proj hSigRaw0All onoffevol.onoffRaw01420side
-n/proj hSigRaw1All onoffevol.onoffRaw11420side
-
-echo "Resultats sur ${ncycles} cycles (Ch0, Ch1): mean (${hSigRaw0All.mean}, ${hSigRaw1All.mean}) sigma (${hSigRaw0All.sigma}, ${hSigRaw1All.sigma})"
-
-newwin 1 1
-disp  hSigRaw0All "blue cpts notit"
-disp  hSigRaw1All "red same cpts notit"
-settitle "Mean and Sigma for 1 cycle"
-
-#Ts les 10cycles
-n10cycles = floor(${ncycles}/10)
-x10min = ${hSigRaw0All.xmin}
-x10max = ${hSigRaw0All.xmax}
-n10bin = ${hSigRaw0All.nbin}
-
+c++compile rebining
+c++link rebining.so dorebin
+
+
+
+
+#col 13 = onoffRaw01420side
+#col 14 = onoffRaw11420side
+
+ntcol2var onoffevol 13 linCh0All
+line2vec vecCh0All $linCh0All
+
+ntcol2var onoffevol 14 linCh1All
+line2vec vecCh1All $linCh1All
+
+
+set packVal "1 10 25 50 100"
+set sigmaCh0 ""
+set errsigCh0 ""
+set sigmaCh1 ""
+set errsigCh1 ""
+
+
+foreach ipack ( $packVal )
+
+#Packing per ipack
+  del invec 
+  del outvec
+  cp vecCh0All invec
+  call dorebin $ipack
+  mean0 = ${outvec.sum}/${outvec.size}
+  sigma0 = sqrt(${outvec.sumsq}/${outvec.size}-${mean0}*${mean0})
+  errsig0 = ${sigma0}/sqrt(2*${outvec.size})    
+  set sigmaCh0 "${sigmaCh0} ${sigma0}"  
+  set errsigCh0 "${errsigCh0} ${errsig0}"       
+#
+  del invec
+  del outvec
+  cp vecCh1All invec
+  call dorebin $ipack
+  mean1 = ${outvec.sum}/${outvec.size}
+  sigma1 = sqrt(${outvec.sumsq}/${outvec.size}-${mean1}*${mean1})
+  errsig1 = ${sigma0}/sqrt(2*${outvec.size})    
+  set sigmaCh1 "${sigmaCh1} ${sigma1}"  
+  set errsigCh1 "${errsigCh1} ${errsig1}"
+end
+
+
+newnt nt0 x y ey
+newnt nt1 x y ey
+
+set x0 ( $packVal )
+set y0 ( $sigmaCh0 )
+set ey0 ( $errsigCh0 )
+
+set x1 ( $packVal )
+set y1 ( $sigmaCh1 )
+set ey1 ( $errsigCh1 )
         
-for i 0:${n10cycles}
- ic = ${i}*10+1
- jc = ${ic}+9
- newh1d hSigRaw010c${i} ${x10min} ${x10max} ${n10bin}
- newh1d hSigRaw110c${i} ${x10min} ${x10max} ${n10bin}
-
- n/proj hSigRaw010c${i}  onoffevol.onoffRaw01420side (cycle>=${ic})&&(cycle<=${jc})
- cp hSigRaw010c${i} h0
- mSig010c[i] = ${h0.mean}
-
- n/proj hSigRaw110c${i}  onoffevol.onoffRaw11420side (cycle>=${ic})&&(cycle<=${jc}) 
- cp hSigRaw110c${i} h1
- mSig110c[i] = ${h1.mean}
+for i 0:$#x0
+ line2nt nt0 $x0[i] $y0[i] $ey0[i]
+ line2nt nt1 $x1[i] $y1[i] $ey1[i]
 end
 
 
-line2vec hSigRaw010c $mSig010c
-line2vec hSigRaw110c $mSig110c
-
-
-
-
-newwin 1 1
-disp  hSigRaw010c "blue cpts notit"
-disp  hSigRaw110c "red same cpts notit"
-settitle "Mean and Sigma for 10 cycles"
+newwin 1 1
+plot2de nt0 x y 0 ey 1 "blue marker=fcircle,9 notit nsta"     
+plot2de nt1 x y 0 ey 1 "same red  marker=fcircle,9 notit nsta"
+n = ${#x0}-1
+func $y0[0]/sqrt(x) $x0[0] $x0[n]  100 "same"
+settitle "Sigma [1418,1419]U[1422,1423]MHz ${source} Ch 0 (blue) Ch 1 (red)" ' ' $defatt
+setaxelabels "num of cycles" "Sigma (a.u)" $axedefatt
+
 
 
@@ -389 +414 @@
 
 graphicatt ""
-setaxesatt ""
+setaxesatt "$axedefatt"