Difference between revisions of "Raster/BPM calibration"

From Hall A Wiki
Jump to: navigation, search
 
(20 intermediate revisions by 2 users not shown)
Line 1: Line 1:
−
Bull's eye scan run
+
*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]]
 +
{| border="1"
 +
|+ 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 size runs  
−
{!
+
 
−
! run  !! comment !!         
+
{| border="1"
 +
|+ Raster runs
 
|-
 
|-
−
| 1075 || Raster Size 3X3 ||
+
| 1075 || Raster Size 3X3  
 
|-
 
|-
−
| 1076 || Raster Size 7X9||
+
| 1076 || Raster Size 7X9
 
|-
 
|-
−
| 1077 || Raster Size 8X11||
+
| 1077 || Raster Size 8X11
 
|-
 
|-
−
| 1078 || Raster Size 9X11 ||
+
| 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
 +
[[http://www.jlab.org/~adaq/halog/html/0703_archive/070304080440.html]]
 +
 +
  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.)");

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]]

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.)");

Retrieved from "https://hallaweb.jlab.org/wiki/index.php?title=Raster/BPM_calibration&oldid=1637"