;ò žM9Jc@s]dkZdklZdkTd„Zd„Zddd„Zd„Zd „Zd „Z dS( N(ssqrt(s*c Cs¯|iƒ}|iƒ|iƒ|iƒf\}}xjt|ƒD]\}xSt|ƒD]E}|i ||ƒt |ƒ|i ||ƒ}|i|||ƒqTWqAW|iƒSdS(N(sh1sCloneshdiffsResets GetNbinsXs GetNbinsYsnxsnysrangesibxsibys GetBinContentsfloatsash2scs SetBinContent( sh1sh2sasnxscsibyshdiffsnysibx((sC/home/jbueno/Reports/ts/results/nice_residuals_plots/accessories.pys add_histoss    ,cCsÈ|dt|ƒjo|djn otSn||dj otSn||}t|ƒ|dj otSnt|d|dƒ}|d|jo|djn otSntSdS(Niiiiii( scvalssabssctsFalsespspzssqrtsptsTrue(spsctscvalssptspz((sC/home/jbueno/Reports/ts/results/nice_residuals_plots/accessories.pysin_fids- 'f22.998999999999999f27.001000000000001c Csõ|iƒ}|iƒ|iƒ|iƒf\}} x°t|ƒD]¢}|i ƒi |ƒ}||jo |jnoqAnx_t| ƒD]Q}|iƒi |ƒ} t|| |ƒo#|i|||i||ƒƒqŽqŽWqAW|iƒSdS(N(sh1sClonesh2sResets GetNbinsXs GetNbinsYsnxsnysrangesibxsGetXaxiss GetBinCenterspspcut_minspcut_maxsibysGetYaxissctsin_fidsfids SetBinContents GetBinContent( sh1sfidspcut_minspcut_maxsibxspsh2sibysnxsnysct((sC/home/jbueno/Reports/ts/results/nice_residuals_plots/accessories.pysfid_onlys    +cCsÑtdddddƒ}|iƒidƒxxt|iƒƒD]d}x[t|iƒƒD]G}|i ||ƒ}d|jo djno|i |ƒqWqWWq>WdG|i ƒGd G|iƒGH|iƒSdS( Nssimpleidif0.01iùf0.0f1.0smean=s; RMS=(sTH1DshsGetXaxiss SetNdivisionssrangesh1s GetNbinsXsibxs GetNbinsYsibys GetBinContentscsFillsGetMeansGetRMSsClone(sh1scshsibxsiby((sC/home/jbueno/Reports/ts/results/nice_residuals_plots/accessories.pys simple_histo/scCs%|iƒiƒ}|iƒxÞt|iƒƒD]Ê}g}x[t|i ƒƒD]G}|i ||ƒ}d|jo djno|i |ƒqNqNWt|ƒdjoM|i|t|ƒt|ƒƒ|i|ti|ƒtt|ƒƒƒq/q/W|idƒ|idƒ|iƒSdS(Nf0.0f1.0iif0.0050000000000000001(sh1s ProjectionXsClonesh2sResetsranges GetNbinsXsibxsxs GetNbinsYsibys GetBinContentscsappendslens SetBinContentssums SetBinErrorsnumpysstdssqrtsSetMarkerStyles SetMaximum(sh1scsh2sibxsibysx((sC/home/jbueno/Reports/ts/results/nice_residuals_plots/accessories.pys project_vs_p:s   1  cCs%|iƒiƒ}|iƒxÞt|iƒƒD]Ê}g}x[t|i ƒƒD]G}|i ||ƒ}d|jo djno|i |ƒqNqNWt|ƒdjoM|i|t|ƒt|ƒƒ|i|ti|ƒtt|ƒƒƒq/q/W|idƒ|idƒ|iƒSdS(Nf0.0f1.0iif0.0050000000000000001(sh1s ProjectionYsClonesh2sResetsranges GetNbinsYsibysxs GetNbinsXsibxs GetBinContentscsappendslens SetBinContentssums SetBinErrorsnumpysstdssqrtsSetMarkerStyles SetMaximum(sh1scsh2sibxsibysx((sC/home/jbueno/Reports/ts/results/nice_residuals_plots/accessories.pys project_vs_ctKs   1  ( snumpysmathssqrtsROOTs add_histossin_fidsfid_onlys simple_histos project_vs_ps project_vs_ct(s project_vs_ctsnumpyssqrts simple_histosin_fids project_vs_ps add_histossfid_only((sC/home/jbueno/Reports/ts/results/nice_residuals_plots/accessories.pys?s