BH event in root file my working dir at jlab ifarm /work/halla/solid/zwzhao/BH/grape-dilepton_work from count to asymmetry stat error, as mentioned in Shengyin's thesis eq4.7 delta_A=sqrt(1-A*A)/sqrt(N) histogram of count are stored in root file = try1 ================================= from grape output "JLab_11_BH_muon_3fold_decaypairelectron_deg5-50" done with 3 acceptance, FA, FA+LA,FA+LA+Barrel. so just focus on FA phi angle "Phi_LH" according to Belitsky, differ from trento kind convention by pi for the 3 bins used in the EPJA paper, as defined // if(abs(xi-0.135)<0.015 && abs(-t-0.25)<0.05){ // if(abs(xip+0.06)<0.03 && abs(Q2-1.25)<0.25) {hcount1->Fill(Phi_LH*DEG,thisweight);} // if(abs(xip-0.075)<0.015 && abs(Q2-1.25)<0.25) {hcount2->Fill(Phi_LH*DEG,thisweight);} // if(abs(xip-0.075)<0.015 && abs(Q2-1.875)<0.375) {hcount3->Fill(Phi_LH*DEG,thisweight);} // } hcount1 is for black curve, hcount2 is for red curve, hcount3 is for blue curve don't use hcount2 as its counts is too low, just focus on hcount1 and hcount3 = try2 ================================== from grape output "JLab_11_BH_muon_3fold_decaypairelectron_deg5-50" use acceptance FA only phi angle "Phi_LH" is from lepton plane to recoil proton, kinda like trento convention for SIDIS at PhysRevD.70.117504 which has phi angle from lepton plane outgoing hadron 2 bins as defined // if(abs(xi-0.135)<0.015 && abs(-t-0.25)<0.05){ // if(abs(xip+0.06)<0.03 && abs(Q2-1.25)<0.25) {hcount[0]->Fill(Phi_LH*DEG,weight[windex]);} // if(abs(xip-0.075)<0.05 && abs(Q2-2.5)<1) {hcount[1]->Fill(Phi_LH*DEG,weight[windex]);} // } hcount0 in black is the same as try1 hcount1 hcount1 in red is enlarged comparing to try1 hcount3 == try3 ============================================ from grape output "JLab_11_BH_muon_4fold_nocut" no acceptance and no luminosity hcount0 shows BH total crosssection in phi (trento convention) this histogram needs to divide the bin size 12deg to get results in unit of pb/deg == try4 ============================================ from grape output "JLab_11_BH_muon_4fold_nocut" same like try3, but with cuts as follows 0.05 < Q2 < 10.0 2.0 < Qp2 < 10.0 -1.2 < t < 0.0 0.001 < xi < 0.7 "nocut" generator input has no any cut "nocut_noISR" like above, and initial state radiation is off for incoming e- "cutQ2" generator input has cut 0.05 < Q2 < 10.0 "cutQ2_noISR" like above, and initial state radiation is off for incoming e- == try5 ============================================== 2024/11/13 from grape output "JLab_11_BH_muon_3fold_decaypairelectron_deg5-50" output.root has condition of 11GeV beam,1.2e37/cm2/s,50 days,0.7 overall efficiency histogram "lepIM1_3" with acceptance condition of e- angle7-30deg mom>1, mu+mu- angle 7-50deg mom>1 histogram "lepIM1_4" with condition of last one and W>2,1e-4<-t<1,0.052,Q2+Qp2>1,W>2 == try6 ================================================== from grape output JLab_11_BH_muon_3fold_decaypairelectron_deg5-50 JLab_11_CO_muon_3fold_decaypairelectron_deg5-50 JLab_11_QED_muon_3fold_decaypairelectron_deg5-50 JLab_11_EW_muon_3fold_decaypairelectron_deg5-50 compare_grape_mode.png, compare different grape running mode GRASLL, refer to https://research.kek.jp/people/tabe/grape/CPC/node5.html so to include all channels,"QED" or "EW" should be used for rate estimation instead of "BH" from grape output "HERA_EW_muon_e27.6p920_massll5-100" compare_grape_HERA.png, compare grape EW mode with HERA data (https://doi.org/10.1016/j.physletb.2004.01.003) it has e 27.6GeV p 920GeV, two muon angle theta 20-160deg, mum.Pt > 2GeV,mup.Pt>1.75GeV instead of mu1 and mu2 mentioned in paper