Difference between revisions of "Raster/BPM calibration"
| (3 intermediate revisions by the same user not shown) | |||
| Line 1: | Line 1: | ||
*To determine raster calibration [[ ]] | *To determine raster calibration [[ ]] | ||
| − | in | + | in directory onlana/BPM run get_rast_const.C. |
An edit the file db_rb.dat and replace the three lines | An edit the file db_rb.dat and replace the three lines | ||
| − | + | *[Raster] | |
| − | -23.0 18.3 18.3 0 0 1633 1643 | + | *-23.0 18.3 18.3 0 0 1633 1643 |
| − | -7.345 | + | *-7.345 |
| − | -2.214 | + | *-2.214 |
| − | 0.0 | + | *0.0 |
| − | -0.00475791 -0.00278243 1.63825e-06 1.83396e-06 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.00186778 -0.00286444 2.23071e-06 2.14824e-06 0. 0. <- |
| − | -0.000647539 -0.0029319 2.50367e-06 2.30482e-06 0. 0. <- | + | *-0.000647539 -0.0029319 2.50367e-06 2.30482e-06 0. 0. <- |
| − | ! bpma lines | + | *! bpma lines |
| + | |||
| + | I changed the constants for the 2.3 GeV running, previous constants are saved in db_rb_1.3GeV.dat, all files are in $DB_DIR. | ||
| + | |||
by the one given from the script. | by the one given from the script. | ||
Latest revision as of 04:10, 5 March 2007
- To determine raster calibration [[ ]]
in directory onlana/BPM run get_rast_const.C. An edit the file db_rb.dat and replace the three lines
- [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
I changed the constants for the 2.3 GeV running, previous constants are saved in db_rb_1.3GeV.dat, all files are in $DB_DIR.
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]]
| 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
| 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.)");
