Raster/BPM calibration

• To determine raster calibration [[ ]]

in onlana directory/BPM run get_rast_const.C. An edit the file db_rb.dat and replace the three lines Insert non-formatted text [Raster] -23.0 18.3 18.3 0 0 1633 1643 -7.345 -2.214 0.0 -0.00475791 -0.00278243 1.63825e-06 1.83396e-06 0. 0. <- -0.00186778 -0.00286444 2.23071e-06 2.14824e-06 0. 0. <- -0.000647539 -0.0029319 2.50367e-06 2.30482e-06 0. 0. <- ! bpma lines

by the one given from the script.

Bull's eye scan run [[ http://www.jlab.org/~adaq/halog/html/0703_archive/070303205118.html Halog entry about bull's eyes scan]]

Harps scans/coda runs

Position X

Harp X A

sigma X A

Position Y A

Harp Y A

sigma Y A

1078	-4	2	1733	30	2182	72.31
1079	-4	2	1733	30	2182	72.31
1080	-2	-1	1733	30	2182	72.31
1081	-6	-1	1733	30	2182	72.31
1082	-2	5	1733	30	2182	72.31
1083	-6	5	1733	30	2182	72.31

Raster size runs

Raster runs

1075	Raster Size 3X3
1076	Raster Size 7X9
1077	Raster Size 8X11
1078	Raster Size 9X11

---

• from SKxceb

Double_t right_optic_dp_correction = 0;
 Double_t right_optic_phi_correction = 0; 
 Double_t left_optic_dp_correction = 0;
 Double_t left_optic_phi_correction = 0; 
 Double_t raster_dp_correction = 0.;
 Double_t raster_th_correction = 0;
 Double_t raster_ph_correction = 0;
 //Double_t DeltaDPCorr = 5.18;
 Double_t DeltaDPCorr = 5.18;
 //Double_t DeltaThCorr = 0.61;
 //Double_t DeltaThCorr = 5.124;
 Double_t DeltaThCorr = 5.124;
 Double_t DeltaPhCorr = 0.; 
 Double_t right_dp, right_th, right_ph, left_dp, left_th, left_ph;
 //Double_t right_angle = -0.1047198;
 //Double_t left_angle = 0.1047198;
 Double_t right_angle = -0.513884;
 Double_t left_angle = 0.76683;

if(Flag_Raster_Correction == 1)

 {
    // Raster Corrections
    
  raster_dp_correction = - ReactPt_L_y/DeltaDPCorr;
  raster_th_correction = - ReactPt_L_y*DeltaThCorr - 1700.*ReactPt_L_y*ReactPt_L_y;
  raster_ph_correction = - ReactPt_L_y*DeltaPhCorr; 
 }

 right_dp = R_gold_dp + right_optic_dp_correction + raster_dp_correction;
 right_th = R_gold_th + raster_th_correction;
 right_ph = R_gold_ph + right_optic_phi_correction + raster_ph_correction;
 left_dp = L_gold_dp + left_optic_dp_correction + raster_dp_correction;
 left_th = L_gold_th + raster_th_correction;
 left_ph = L_gold_ph + left_optic_phi_correction + raster_ph_correction;

• Raster/current calibration

T->Draw("(rb.Raster.rawcur.y-2150)/900*3+2:6+2*(rb.Raster.rawcur.x-2250)/1150")

<SRC="http://www.jlab.org/~adaq/halog/html/0703_archive/070304080440.1.gif">

• Raster correction bases on current

[[1]]

 T->Draw("(R.gold.dp+1.)*1340*(1.+2.*1343/11177.928*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5)*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5))/(1.+2.*1343/11177.928*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Rp1");

T->Draw("(R.gold.dp+1.+ ((rb.Raster.rawcur.y-2150)/900*3+2)/5180+ 0.0041*R.gold.ph*1.+1.*0.1*R.gold.ph**2+0.*4.912*R.gold.dp**2)*1340*(1.+2.*1343/11177.928*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5))/(1.+2.*1343/11177.928*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Rp1c(2000,1200,1400.)");