TFile *_file0 = TFile::Open("CREX-Good-miniruns_basic.root") // Get us_avg global average agg->Draw("rcdb_sign*cor_asym_us_avg_mean:cor_asym_us_avg_mean_error:run_number+0.1*minirun_n","rcdb_arm_flag==0","goff") TGraphErrors* tge = new TGraphErrors(8534,agg->GetV3(),agg->GetV1(),0,agg->GetV2()) tge = new TGraphErrors(8406,agg->GetV3(),agg->GetV1(),0,agg->GetV2()) tge->Fit("pol0") /**************************************** Minimizer is Linear Chi2 = 8445.61 NDf = 8405 p0 = 2.09018e-06 +/- 8.60713e-08*/ gStyle->SetOptFit(1) tge->Draw("ap") TCanvas * c2 = new TCanvas() // Fractional error 86.0/2090.0 // Get us_avg + usr + usl global average agg->Draw("rcdb_sign*((rcdb_arm_flag==0)*cor_asym_us_avg_mean+(rcdb_arm_flag==1)*cor_asym_usr_mean+(rcdb_arm_flag==2)*cor_asym_usl_mean):(rcdb_arm_flag==0)*cor_asym_us_avg_mean_error+(rcdb_arm_flag==1)*cor_asym_usr_mean_error+(rcdb_arm_flag==2)*cor_asym_usl_mean_error:run_number+0.1*minirun_n","","goff") tge = new TGraphErrors(8534,agg->GetV3(),agg->GetV1(),0,agg->GetV2()) tge->Fit("pol0") /***************************************** Minimizer is Linear Chi2 = 8579.31 NDf = 8533 p0 = 2.09199e-06 +/- 8.57467e-08*/ tge->Draw("ap") // Compare error bars from us_avg only to us_avg + USR + USL data -> 7.5 ppb quad error difference, 2 ppb mean difference (i.e. it is well within the range it should be) sqrt(pow(8.60713e-08,2)-pow(8.57467e-08,2)) TCanvas * c3 = new TCanvas() // Get the eigen-dit vs. regression difference (multipletwise), self weighted global average agg->Draw("rcdb_sign*((rcdb_arm_flag==0)*asym_us_avg_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==1)*asym_usr_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==2)*asym_usl_eigen_dit_vs_cor_diff_mean):(rcdb_arm_flag==0)*asym_us_avg_eigen_dit_vs_cor_diff_mean_error+(rcdb_arm_flag==1)*asym_usr_eigen_dit_vs_cor_diff_mean_error+(rcdb_arm_flag==2)*asym_usl_eigen_dit_vs_cor_diff_mean_error:run_number+0.1*minirun_n","","goff") tge = new TGraphErrors(8534,agg->GetV3(),agg->GetV1(),0,agg->GetV2()) tge->Fit("pol0") /**************************************** Minimizer is Linear Chi2 = 8667.04 NDf = 8533 p0 = -5.07888e-09 +/- 6.52643e-09*/ tge->Draw("ap") TCanvas * c4 = new TCanvas() // Get the main-det weighting factor agg->Draw("pow(((rcdb_arm_flag==0)*asym_us_avg_eigen_dit_vs_cor_diff_mean_error+(rcdb_arm_flag==1)*asym_usr_eigen_dit_vs_cor_diff_mean_error+(rcdb_arm_flag==2)*asym_usl_eigen_dit_vs_cor_diff_mean_error)/((rcdb_arm_flag==0)*cor_asym_us_avg_mean_error+(rcdb_arm_flag==1)*cor_asym_usr_mean_error+(rcdb_arm_flag==2)*cor_asym_usl_mean_error),2):(rcdb_arm_flag==0)*cor_asym_us_avg_mean_error+(rcdb_arm_flag==1)*cor_asym_usr_mean_error+(rcdb_arm_flag==2)*cor_asym_usl_mean_error:run_number+0.1*minirun_n","","goff") tge = new TGraphErrors(8534,agg->GetV3(),agg->GetV1(),0,agg->GetV2()) tge->Fit("pol0") /***************************************** * Minimizer is Linear * Chi2 = 9.36869e+10 * NDf = 8533 * p0 = 0.0128333 +/- 8.57467e-08*/ // Get the eigen-dit vs. regression difference (multipletwise), main-det weighted global average agg->Draw("rcdb_sign*((rcdb_arm_flag==0)*asym_us_avg_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==1)*asym_usr_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==2)*asym_usl_eigen_dit_vs_cor_diff_mean):sqrt(0.0128333)*((rcdb_arm_flag==0)*cor_asym_us_avg_mean_error+(rcdb_arm_flag==1)*cor_asym_usr_mean_error+(rcdb_arm_flag==2)*cor_asym_usl_mean_error):run_number+0.1*minirun_n","","goff") tge = new TGraphErrors(8534,agg->GetV3(),agg->GetV1(),0,agg->GetV2()) tge->Fit("pol0") /***************************************** * Minimizer is Linear * Chi2 = 8573.29 * NDf = 8533 * p0 = 7.01723e-09 +/- 9.71374e-09*/ tge->Draw("ap") // Look explicitly at the 4 sigma outliers... etc. agg->Scan("run_number:minirun_n:rcdb_slug:rcdb_sign*((rcdb_arm_flag==0)*asym_us_avg_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==1)*asym_usr_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==2)*asym_usl_eigen_dit_vs_cor_diff_mean):sqrt(0.0128333)*((rcdb_arm_flag==0)*cor_asym_us_avg_mean_error+(rcdb_arm_flag==1)*cor_asym_usr_mean_error+(rcdb_arm_flag==2)*cor_asym_usl_mean_error):run_number+0.1*minirun_n","abs(((rcdb_sign*((rcdb_arm_flag==0)*asym_us_avg_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==1)*asym_usr_eigen_dit_vs_cor_diff_mean+(rcdb_arm_flag==2)*asym_usl_eigen_dit_vs_cor_diff_mean))-7.017e-9)/(sqrt(0.0128333)*((rcdb_arm_flag==0)*cor_asym_us_avg_mean_error+(rcdb_arm_flag==1)*cor_asym_usr_mean_error+(rcdb_arm_flag==2)*cor_asym_usl_mean_error)))>4","") /**************************************************************************************** * * Row * run_numbe * minirun_n * rcdb_slug * rcdb_sign * sqrt(0.01 * run_numbe * Cameron's Comment * ************************************************************************************ * * 116 * 5459 * 3 * 102 * -6.18e-06 * 9.369e-07 * 5459.3 * Slow 1mm 12X energy drift * * 131 * 5461 * 0 * 102 * -4.24e-06 * 9.404e-07 * 5461 * Small amount of 12X energy noise * * 384 * 5772 * 0 * 109 * 3.374e-06 * 6.757e-07 * 5772 * Questionable: Short run, residual L R correlation, medium fast 12X energy drift, bad 12YM wire? * * 391 * 5775 * 0 * 109 * 4.482e-06 * 9.473e-07 * 5775 * Slow 12X energy drift, rcdb claims "no BMOD", residual L R correlation * * 393 * 5775 * 2 * 109 * -4.04e-06 * 9.335e-07 * 5775.2 * "" * * 400 * 5776 * 4 * 109 * -5.06e-06 * 7.856e-07 * 5776.4 * Slow 12X energy drift, residual L R correlation * * 593 * 5814 * 5 * 113 * -3.78e-06 * 8.606e-07 * 5814.5 * Looks normal... * * 1057 * 6177 * 0 * 128 * 5.199e-06 * 1.045e-06 * 6177 * Shifts during FFB on/off... looks good * * 1516 * 6516 * 1 * 138 * -4.68e-06 * 9.412e-07 * 6516.1 * Shifts during FFB on/off... looks good * * 1540 * 6520 * 1 * 138 * -3.93e-06 * 8.362e-07 * 6520.1 * Target BPM drifts * * 1542 * 6526 * 0 * 138 * 4.005e-06 * 9.394e-07 * 6526 * BPM 12 shift, otherwise good * * 1857 * 6608 * 0 * 142 * 1.276e-05 * 1.200e-06 * 6608 * Questionable: Half minirun, really noisy beam positions * * 1897 * 6632 * 5 * 142 * 3.186e-06 * 7.816e-07 * 6632.5 * Shifts during FFB on/off... looks good * * 2028 * 6644 * 7 * 144 * -3.26e-06 * 8.161e-07 * 6644.7 * Looks normal... * * 2162 * 6711 * 0 * 146 * 3.849e-06 * 9.422e-07 * 6711 * Shifts during FFB on/off... some residual L R correlation, looks good * * 3597 * 7069 * 8 * 162 * 3.927e-06 * 8.500e-07 * 7069.8 * Looks normal... * * 4008 * 7153 * 4 * 167 * -5.38e-06 * 1.008e-06 * 7153.4 * Jump in 12X during FFB off * * 4579 * 7262 * 1 * 174 * -3.98e-06 * 9.111e-07 * 7262.1 * Slightly messy target BPMs... looks good * * 4693 * 7294 * 1 * 175 * 4.468e-06 * 9.239e-07 * 7294.1 * Slightly messy target BPms and a 12X jump during FFB off * * 4714 * 7296 * 5 * 176 * 3.803e-06 * 9.250e-07 * 7296.5 * Looks normal * * 4717 * 7296 * 8 * 176 * 4.068e-06 * 7.997e-07 * 7296.8 * Looks normal * * 4725 * 7297 * 7 * 176 * -4.94e-06 * 9.281e-07 * 7297.7 * A bit of 12X noise... looks good * * 5012 * 7350 * 1 * 179 * -5.11e-06 * 9.418e-07 * 7350.1 * Messy target BPMs * * 5018 * 7351 * 1 * 180 * -4.17e-06 * 9.231e-07 * 7351.1 * Looks normal * * 5019 * 7351 * 2 * 180 * 4.437e-06 * 9.223e-07 * 7351.2 * "" * * 5026 * 7352 * 2 * 180 * 5.148e-06 * 9.307e-07 * 7352.2 * FFB off 12X jumps * * 5108 * 7368 * 7 * 180 * 4.039e-06 * 9.271e-07 * 7368.7 * Some noisy target BPms * * 5124 * 7370 * 5 * 181 * -3.96e-06 * 9.312e-07 * 7370.5 * Looks good * * 5127 * 7371 * 0 * 181 * 5.321e-06 * 9.246e-07 * 7371 * Looks good * * 5129 * 7371 * 2 * 181 * -4.69e-06 * 9.679e-07 * 7371.2 * "" * * 5131 * 7371 * 4 * 181 * 6.124e-06 * 9.346e-07 * 7371.4 * "" * * 5138 * 7372 * 2 * 181 * -3.98e-06 * 9.284e-07 * 7372.2 * A bit of BPM noise * * 5160 * 7374 * 7 * 181 * 4.275e-06 * 9.322e-07 * 7374.7 * Noisy target BPMs * * 5161 * 7374 * 8 * 181 * 3.494e-06 * 6.656e-07 * 7374.8 * "" * * 5172 * 7376 * 1 * 181 * -4.02e-06 * 9.230e-07 * 7376.1 * Noisy target BPMs * * 5191 * 7378 * 3 * 181 * -4.09e-06 * 9.174e-07 * 7378.3 * Noisy target BPMs * * 5283 * 7402 * 5 * 183 * 3.781e-06 * 9.142e-07 * 7402.5 * Noisy target BPMs (all prior runs seem to have nosiy 4eX) * * 5547 * 7700 * 0 * 188 * 6.416e-06 * 1.306e-06 * 7700 * Pretty short run (half minirun) and bmod running data only * * 5611 * 7719 * 2 * 189 * -4.10e-06 * 9.178e-07 * 7719.2 * Large BPM Y noise * * 5613 * 7719 * 4 * 189 * -4.36e-06 * 9.273e-07 * 7719.4 * "" * * 5623 * 7735 * 1 * 189 * 3.854e-06 * 9.142e-07 * 7735.1 * Looks normal * * 5666 * 7759 * 0 * 190 * 2.679e-06 * 6.668e-07 * 7759 * A bit noisy in 12X, amidst 12X shifts (cut out already) * * 6095 * 7874 * 2 * 196 * 4.489e-06 * 9.213e-07 * 7874.2 * Residual L R correlations, noisy target and 12X BPMs * * 6150 * 7897 * 4 * 197 * -3.49e-06 * 8.317e-07 * 7897.4 * Some jumps in BPMs * * 6199 * 7905 * 3 * 197 * -5.18e-06 * 9.133e-07 * 7905.3 * Noisy BPM 12X, large L R residual correlation * * 6231 * 7909 * 3 * 198 * -4.24e-06 * 9.218e-07 * 7909.3 * Noisy BPM 12X, and shifts around with hairs too * * 6235 * 7910 * 0 * 198 * 7.474e-06 * 9.502e-07 * 7910 * Noisy BPM 12X, large L R residual correlation * * 6242 * 7910 * 7 * 198 * 4.238e-06 * 7.619e-07 * 7910.7 * "" * * 6934 * 8082 * 0 * 206 * 1.188e-05 * 1.316e-06 * 8082 * Super short (half minirun) run * * 7170 * 8143 * 1 * 208 * 4.092e-06 * 9.276e-07 * 8143.1 * Just a bit noisy in the PBMs... looks good * * 7173 * 8143 * 4 * 208 * 5.249e-06 * 9.311e-07 * 8143.4 * "" * * 7175 * 8143 * 6 * 208 * -4.60e-06 * 9.296e-07 * 8143.6 * "" * * 7176 * 8143 * 7 * 208 * 3.616e-06 * 7.111e-07 * 8143.7 * ""*/