Difference between revisions of "CDC B field"
(→56% Ar, 44% CO2, 1200A) |
|||
(54 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
+ | __TOC__ | ||
+ | |||
+ | === B field in CDC volume === | ||
+ | |||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
|+Magnetic field | |+Magnetic field | ||
Line 7: | Line 11: | ||
At large drift radius, change in B field of 0.1T causes change in drift time of ~ 5ns. This would cause approx 30um error in drift radius. Resolution in 0T is 80um-100um at large drift radius. | At large drift radius, change in B field of 0.1T causes change in drift time of ~ 5ns. This would cause approx 30um error in drift radius. Resolution in 0T is 80um-100um at large drift radius. | ||
+ | |||
Plan: | Plan: | ||
Line 32: | Line 37: | ||
|[[Image:cdc_bfield_100um_tdiff.png|thumb|x250px|Drift radius change 100um]] | |[[Image:cdc_bfield_100um_tdiff.png|thumb|x250px|Drift radius change 100um]] | ||
|} | |} | ||
+ | |||
+ | |||
+ | <hr/> | ||
+ | === Finding drift time difference wrt 1.3T === | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Drift times from Garfield | ||
+ | |[[Image:cdc_bfield_orig.png|thumb|x250px|Garfield drift times]] | ||
+ | |[[Image:cdc_bfield_origz2.png|thumb|x250px|Zoomed in at small radius]] | ||
+ | |[[Image:cdc_bfield_origz.png|thumb|x250px|Zoomed in at large radius]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Drift times from Garfield, edited | ||
+ | |[[Image:cdc_bfield_edited.png|thumb|x250px|Zoomed in at small radius]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Individual fits to drift time difference [0]*radius**[1] | ||
+ | |[[Image:cdc_bfield_indiv_fits.png|thumb|x250px|Fit f(r) to each curve]] | ||
+ | |[[Image:cdc_bfield_fitpar0.png|thumb|x250px|Fit of param 0 values]] | ||
+ | |[[Image:cdc_bfield_fitpar1.png|thumb|x250px|Fit of param 1 values]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+One fit to drift time difference, excluded 1.5T data from fit | ||
+ | |[[Image:cdc_bfield_one_fit.png|thumb|x250px|One fit f(B,r) for all curves (except 1.5T)]] | ||
+ | |} | ||
+ | |||
+ | Max error between fitted function and time difference is 5.5ns at radius 7.6mm in 1.5T dataset | ||
+ | |||
+ | |||
+ | <hr/> | ||
+ | |||
+ | === Finding drift time difference wrt 1.3T (more data) === | ||
+ | |||
+ | Repeated the above with more datasets | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Drift times from Garfield | ||
+ | |[[Image:cdc_bfield2_origz.png|thumb|x250px|Zoomed in at large radius]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Drift times from Garfield, edited | ||
+ | |[[Image:cdc_bfield2_editedz.png|thumb|x250px|Zoomed in at large radius]] | ||
+ | |[[Image:cdc_bfield2_edited.png|thumb|x250px|]] | ||
+ | |[[Image:cdc_bfield2_edited_zz.png|thumb|x250px|Zoomed in at large radius, key reversed]] | ||
+ | |} | ||
+ | |||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Individual fits to drift time difference [0]*driftradius**[1] | ||
+ | |[[Image:cdc_bfield2_indiv_fits.png|thumb|x250px|Fit f(driftradius) to each curve]] | ||
+ | |[[Image:cdc_bfield2_fitpar0.png|thumb|x250px|Fit of param 0 values]] | ||
+ | |[[Image:cdc_bfield2_fitpar1.png|thumb|x250px|Fit of param 1 values]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+One fit to drift time difference, excluded 1.5T data from fit | ||
+ | |[[Image:cdc_bfield2_one_fit.png|thumb|x250px|One fit f(B,driftradius) for all curves]] | ||
+ | |} | ||
+ | |||
+ | Max error using one fit f(B,driftradius) for all curves is 2.4ns at 7.6mm for 1.7T | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | <pre> | ||
+ | Param 0 fit status 0 function 10.72 + -19.17x + 10.68x^2 + -1.74x^3 | ||
+ | Param 1 fit status 0 function 9.88 + -8.22x + 2.06x^2 | ||
+ | |||
+ | 1.35T Global fit av abs err 0.4ns max abs err -1.2ns at radius 5.8mm | ||
+ | 1.40T Global fit av abs err 0.2ns max abs err -0.5ns at radius 4.5mm | ||
+ | 1.45T Global fit av abs err 0.3ns max abs err 2.0ns at radius 7.6mm | ||
+ | 1.50T Global fit av abs err 0.3ns max abs err 0.6ns at radius 2.4mm | ||
+ | 1.55T Global fit av abs err 0.4ns max abs err -1.1ns at radius 6.4mm | ||
+ | 1.60T Global fit av abs err 0.3ns max abs err 1.1ns at radius 7.6mm | ||
+ | 1.65T Global fit av abs err 0.4ns max abs err 1.2ns at radius 7.0mm | ||
+ | 1.70T Global fit av abs err 0.4ns max abs err 2.4ns at radius 7.6mm | ||
+ | 1.75T Global fit av abs err 0.5ns max abs err -1.1ns at radius 6.1mm | ||
+ | 1.80T Global fit av abs err 0.3ns max abs err 1.4ns at radius 7.6mm | ||
+ | 1.85T Global fit av abs err 0.6ns max abs err -1.4ns at radius 6.2mm | ||
+ | 1.90T Global fit av abs err 0.3ns max abs err 1.2ns at radius 7.6mm | ||
+ | </pre> | ||
+ | |||
+ | |||
+ | === Finding B === | ||
+ | |||
+ | |||
+ | |||
+ | <div style="background-color:gray"> | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
Line 99: | Line 196: | ||
|} | |} | ||
− | < | + | <pre> |
+ | Mostly backward angle z (16 to 86cm) pol3 fit parameter fit parameters as function of r | ||
+ | param [0] [1] r [2] r^2 [3] r^3 | ||
+ | 0 1.045e+00 2.290e-04 1.972e-05 5.338e-07 | ||
+ | 1 1.518e-02 1.009e-04 -2.982e-06 7.920e-08 | ||
+ | 2 -1.157e-04 -3.444e-06 1.134e-07 -2.934e-09 | ||
+ | 3 3.657e-07 2.528e-08 -8.925e-10 2.196e-11 | ||
+ | Mostly forward angle z (84 to 168cm) pol3 fit parameter fit parameters as function of r | ||
+ | param [0] [1] r [2] r^2 [3] r^3 | ||
+ | 0 1.451e+00 -4.060e-02 1.916e-03 -3.024e-05 | ||
+ | 1 2.443e-03 1.017e-03 -4.767e-05 7.464e-07 | ||
+ | 2 1.570e-05 -8.240e-06 3.862e-07 -5.966e-09 | ||
+ | 3 -8.188e-08 2.168e-08 -1.019e-09 1.552e-11 | ||
+ | |||
+ | Double_t lo_z0[4]={1.045E+00, 2.290E-04, 1.972E-05, 5.338E-07}; | ||
+ | Double_t lo_z1[4]={1.518E-02, 1.009E-04, -2.982E-06, 7.920E-08}; | ||
+ | Double_t lo_z2[4]={-1.157E-04, -3.444E-06, 1.134E-07, -2.934E-09}; | ||
+ | Double_t lo_z3[4]={3.657E-07, 2.528E-08, -8.925E-10, 2.196E-11}; | ||
+ | |||
+ | Double_t hi_z0[4]={1.451E+00, -4.060E-02, 1.916E-03, -3.024E-05}; | ||
+ | Double_t hi_z1[4]={2.443E-03, 1.017E-03, -4.767E-05, 7.464E-07}; | ||
+ | Double_t hi_z2[4]={1.570E-05, -8.240E-06, 3.862E-07, -5.966E-09}; | ||
+ | Double_t hi_z3[4]={-8.188E-08, 2.168E-08, -1.019E-09, 1.552E-11}; | ||
+ | </pre> | ||
+ | |||
+ | so for z<85cm, B ~ (1.045 + 2.290e-04r + 1.972e-05r^2 + 5.338e-07r^3) + (1.518e-02 + 1.009e-04r - 2.982e-06r^2 + 7.920e-08r^3)z + (-1.157e-04 - 3.444e-06r + 1.134e-07r^2 - 2.934e-09r^3)z^2 | ||
+ | + ( 3.657e-07 + 2.528e-08r - 8.925e-10r^2 + 2.196e-11r^3)z^3 with B in T and z and r in cm | ||
+ | |||
+ | Max difference between fit and original histogram from bfield.root is 0.05T. | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fit function - original Btot histogram |
− | |[[Image: | + | |[[Image:cdc_bfield_bdiffs.png|thumb|x190px| Fit - histo (T)]] |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
+ | </div> | ||
+ | |||
+ | The above was for Btot. Repeated for Bz, params below. Max err is 0.0449T at z=22cm, r=55cm | ||
+ | Fits are in bfield/bfits | ||
+ | |||
+ | |||
+ | <pre> | ||
+ | |||
+ | Mostly backward angle z (16 to 86cm) pol3 fit parameters | ||
+ | r (cm) [0] [1] z [2] z^2 [3] z^3 | ||
+ | 10.0 1.046e+00 1.608e-02 -1.432e-04 5.577e-07 | ||
+ | 20.0 1.051e+00 1.699e-02 -1.666e-04 7.033e-07 | ||
+ | 30.0 1.059e+00 1.865e-02 -2.103e-04 9.811e-07 | ||
+ | 35.0 1.065e+00 1.987e-02 -2.435e-04 1.196e-06 | ||
+ | 40.0 1.072e+00 2.139e-02 -2.856e-04 1.472e-06 | ||
+ | 45.0 1.081e+00 2.330e-02 -3.399e-04 1.833e-06 | ||
+ | 50.0 1.091e+00 2.573e-02 -4.099e-04 2.304e-06 | ||
+ | 53.0 1.103e+00 2.880e-02 -5.001e-04 2.916e-06 | ||
+ | 57.0 1.109e+00 3.057e-02 -5.527e-04 3.274e-06 | ||
+ | |||
+ | Mostly forward angle z (84 to 168cm) pol3 fit parameters | ||
+ | r (cm) [0] [1] z [2] z^2 [3] z^3 | ||
+ | 10.0 1.200e+00 8.750e-03 -3.532e-05 5.202e-08 | ||
+ | 20.0 1.183e+00 9.148e-03 -3.779e-05 5.612e-08 | ||
+ | 30.0 1.124e+00 1.055e-02 -4.755e-05 7.711e-08 | ||
+ | 35.0 1.057e+00 1.213e-02 -5.904e-05 1.037e-07 | ||
+ | 40.0 9.488e-01 1.472e-02 -7.826e-05 1.496e-07 | ||
+ | 45.0 7.734e-01 1.892e-02 -1.100e-04 2.272e-07 | ||
+ | 50.0 4.956e-01 2.561e-02 -1.613e-04 3.550e-07 | ||
+ | 53.0 6.308e-02 3.610e-02 -2.429e-04 5.608e-07 | ||
+ | 57.0 -2.243e-01 4.310e-02 -2.977e-04 7.005e-07 | ||
+ | |||
+ | Mostly backward angle z (16 to 86cm) pol3 fit parameter fit parameters as function of r | ||
+ | param [0] [1] r [2] r^2 [3] r^3 | ||
+ | 0 1.045e+00 1.338e-04 4.523e-06 2.368e-07 | ||
+ | 1 1.513e-02 1.160e-04 -3.189e-06 1.053e-07 | ||
+ | 2 -1.146e-04 -3.751e-06 1.237e-07 -3.429e-09 | ||
+ | 3 3.595e-07 2.689e-08 -9.586e-10 2.461e-11 | ||
+ | |||
+ | Mostly forward angle z (84 to 168cm) pol3 fit parameter fit parameters as function of r | ||
+ | param [0] [1] r [2] r^2 [3] r^3 | ||
+ | 0 1.459e+00 -4.190e-02 1.966e-03 -3.085e-05 | ||
+ | 1 2.268e-03 1.043e-03 -4.871e-05 7.577e-07 | ||
+ | 2 1.686e-05 -8.411e-06 3.930e-07 -6.034e-09 | ||
+ | 3 -8.439e-08 2.205e-08 -1.033e-09 1.565e-11 | ||
+ | |||
+ | Double_t lo_z0[4]={1.0445E+00, 1.3383E-04, 4.5227E-06, 2.3683E-07}; | ||
+ | Double_t lo_z1[4]={1.5132E-02, 1.1596E-04, -3.1894E-06, 1.0531E-07}; | ||
+ | Double_t lo_z2[4]={-1.1458E-04, -3.7512E-06, 1.2367E-07, -3.4290E-09}; | ||
+ | Double_t lo_z3[4]={3.5954E-07, 2.6887E-08, -9.5863E-10, 2.4608E-11}; | ||
+ | |||
+ | Double_t hi_z0[4]={1.4594E+00, -4.1900E-02, 1.9661E-03, -3.0849E-05}; | ||
+ | Double_t hi_z1[4]={2.2682E-03, 1.0434E-03, -4.8708E-05, 7.5771E-07}; | ||
+ | Double_t hi_z2[4]={1.6860E-05, -8.4111E-06, 3.9302E-07, -6.0335E-09}; | ||
+ | Double_t hi_z3[4]={-8.4393E-08, 2.2049E-08, -1.0333E-09, 1.5649E-11}; | ||
+ | |||
+ | </pre> | ||
+ | |||
+ | === B fits (part 2) === | ||
+ | |||
+ | Bz fit parameters | ||
+ | |||
+ | <pre> | ||
+ | /* Max err 0.067 at z 22.0 r 55.0 */ | ||
+ | Double_t lo_z0[4]={1.0440E+00, 2.8431E-04, 2.2198E-05, 4.4318E-07}; | ||
+ | Double_t lo_z1[4]={1.5235E-02, 8.4686E-05, -2.0227E-06, 6.4062E-08}; | ||
+ | Double_t lo_z2[4]={-1.1729E-04, -2.9396E-06, 8.2057E-08, -2.4271E-09}; | ||
+ | Double_t lo_z3[4]={3.7660E-07, 2.1683E-08, -6.6856E-10, 1.8324E-11}; | ||
+ | |||
+ | Double_t hi_z0[4]={1.4139E+00, -3.3108E-02, 1.5528E-03, -2.5238E-05}; | ||
+ | Double_t hi_z1[4]={3.3934E-03, 8.2637E-04, -3.8508E-05, 6.2072E-07}; | ||
+ | Double_t hi_z2[4]={7.7830E-06, -6.6706E-06, 3.1153E-07, -4.9497E-09}; | ||
+ | Double_t hi_z3[4]={-6.0532E-08, 1.7501E-08, -8.2145E-10, 1.2855E-11}; | ||
+ | </pre> | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fit function - original Bz histogram |
− | |[[Image: | + | |[[Image:cdc_bfield_bdiffs2.png|thumb|x190px| Fit - histo (T)]] |
|} | |} | ||
+ | |||
+ | === B fits (part 3) === | ||
+ | |||
+ | Used pol2 fit for r<41cm and pairs of pol3 fits (as before) for r>=41cm | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fit function - original Bz histogram |
− | |[[Image: | + | |[[Image:cdc_bfield_bdiffs3.png|thumb|x190px| Fit - histo (T)]] |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
+ | Used pol2 fit for r<=30cm and pairs of pol3 fits (as before) for r>30cm | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fit function - original Bz histogram |
− | |[[Image: | + | |[[Image:cdc_bfield_bdiffs3_30.png|thumb|x190px| Fit - histo (T)]] |
|} | |} | ||
− | + | Used pol2 fit for r<=30cm and pairs of pol2 fits for r>30cm | |
+ | {| border="0" cellpadding="2" | ||
+ | |+ Fit function - original Bz histogram | ||
+ | |[[Image:cdc_bfield_bdiffs3_30_pol2.png|thumb|x190px| Fit - histo (T)]] | ||
+ | |} | ||
− | |||
− | + | Used pol2 fit for r<=25cm and pairs of pol2 fits for r>25cm | |
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fits of B field (z) for fixed radius |
− | |[[Image: | + | |[[Image:cdc_bfield_fit4_r10cm.png|thumb|x90px|10cm]] |
+ | |[[Image:cdc_bfield_fit4_r15cm.png|thumb|x90px|15cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r20cm.png|thumb|x90px|20cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r25cm.png|thumb|x90px|25cm]] | ||
+ | |- | ||
+ | |[[Image:cdc_bfield_fit4_r30cm.png|thumb|x90px|30cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r35cm.png|thumb|x90px|35cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r40cm.png|thumb|x90px|40cm]] | ||
+ | |- | ||
+ | |[[Image:cdc_bfield_fit4_r42_5cm.png|thumb|x90px|42.5cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r45cm.png|thumb|x90px|45cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r47_5cm.png|thumb|x90px|47.5cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r50cm.png|thumb|x90px|50cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r52_5cm.png|thumb|x90px|52.5cm]] | ||
+ | |[[Image:cdc_bfield_fit4_r55cm.png|thumb|x90px|55cm]] | ||
|} | |} | ||
+ | |||
+ | |||
+ | Fits to fit parameters above | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fit params for r <= 25cm, as func (r) |
− | |[[Image: | + | |[[Image:cdc_bfield_fit4_lowr_par0.png|thumb|x90px| param 0]] |
− | |[[Image: | + | |[[Image:cdc_bfield_fit4_lowr_par1.png|thumb|x90px| param 1 (r)]] |
− | |[[Image: | + | |[[Image:cdc_bfield_fit4_lowr_par2.png|thumb|x90px| param 2 (r^2)]] |
|} | |} | ||
+ | {| border="0" cellpadding="2" | ||
+ | |+ Fit params for mostly backward angles, small z, as func (r) | ||
+ | |[[Image:cdc_bfield_fit4_lowz_par0.png|thumb|x90px| param 0]] | ||
+ | |[[Image:cdc_bfield_fit4_lowz_par1.png|thumb|x90px| param 1 (r)]] | ||
+ | |[[Image:cdc_bfield_fit4_lowz_par2.png|thumb|x90px| param 2 (r^2)]] | ||
+ | |[[Image:cdc_bfield_fit4_lowz_par3.png|thumb|x90px| param 3 (r^3)]] | ||
+ | |} | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fit params for mostly forward angles, large z, as func (r) |
− | |[[Image: | + | |[[Image:cdc_bfield_fit4_hiz_par0.png|thumb|x90px| param 0]] |
− | |[[Image: | + | |[[Image:cdc_bfield_fit4_hiz_par1.png|thumb|x90px| param 1 (r)]] |
− | |[[Image: | + | |[[Image:cdc_bfield_fit4_hiz_par2.png|thumb|x90px| param 2 (r^2)]] |
+ | |[[Image:cdc_bfield_fit4_hiz_par3.png|thumb|x90px| param 3 (r^3)]] | ||
|} | |} | ||
{| border="0" cellpadding="2" | {| border="0" cellpadding="2" | ||
− | |+ | + | |+ Fit function - original Bz histogram |
− | |[[Image: | + | |[[Image:cdc_bfield_bdiffs4_25.png|thumb|x190px| Fit - histo (T)]] |
|} | |} | ||
− | + | This gives results comparable to the original fits with slightly simpler polynomials | |
+ | <pre> | ||
+ | Double_t lo_r0[3]={1.1746E+00, -1.1176E-04, 9.1798E-05}; | ||
+ | Double_t lo_r1[3]={8.6423E-03, 4.5702E-06, -1.4617E-06}; | ||
+ | Double_t lo_r2[3]={-2.5830E-05, -2.6167E-08, 5.9721E-09}; | ||
+ | Double_t lo_z0[3]={1.0753E+00, -2.0457E-03, 7.8510E-05}; | ||
+ | Double_t lo_z1[3]={2.0352E-02, -2.7862E-04, 6.4074E-06}; | ||
+ | Double_t lo_z2[3]={-3.0851E-04, 1.0707E-05, -2.3605E-07}; | ||
+ | Double_t lo_z3[3]={1.8161E-06, -8.1165E-08, 1.7311E-09}; | ||
+ | Double_t hi_z0[3]={-6.1452E-01, 1.1061E-01, -1.7751E-03}; | ||
+ | Double_t hi_z1[3]={5.3402E-02, -2.7137E-03, 4.3399E-05}; | ||
+ | Double_t hi_z2[3]={-3.9186E-04, 2.1597E-05, -3.4203E-07}; | ||
+ | Double_t hi_z3[3]={9.7935E-07, -5.6001E-08, 8.7685E-10}; | ||
+ | if (r<=25) { | ||
− | + | bfit0 = lo_r0[0] + lo_r0[1]*r + lo_r0[2]*r*r; | |
− | + | bfit1 = lo_r1[0] + lo_r1[1]*r + lo_r1[2]*r*r; | |
− | + | bfit2 = lo_r2[0] + lo_r2[1]*r + lo_r2[2]*r*r; | |
+ | bfit3 = 0; | ||
− | 1.35T Global fit av abs err 0. | + | } else { |
− | 1.40T Global fit av abs err 0. | + | |
+ | if (z<85) { | ||
+ | |||
+ | bfit0 = lo_z0[0] + lo_z0[1]*r + lo_z0[2]*r*r; | ||
+ | bfit1 = lo_z1[0] + lo_z1[1]*r + lo_z1[2]*r*r; | ||
+ | bfit2 = lo_z2[0] + lo_z2[1]*r + lo_z2[2]*r*r; | ||
+ | bfit3 = lo_z3[0] + lo_z3[1]*r + lo_z3[2]*r*r; | ||
+ | |||
+ | } else { | ||
+ | |||
+ | bfit0 = hi_z0[0] + hi_z0[1]*r + hi_z0[2]*r*r; | ||
+ | bfit1 = hi_z1[0] + hi_z1[1]*r + hi_z1[2]*r*r; | ||
+ | bfit2 = hi_z2[0] + hi_z2[1]*r + hi_z2[2]*r*r; | ||
+ | bfit3 = hi_z3[0] + hi_z3[1]*r + hi_z3[2]*r*r; | ||
+ | |||
+ | } | ||
+ | } | ||
+ | |||
+ | bfit = bfit0 + bfit1*z + bfit2*z*z + bfit3*z*z*z; | ||
+ | |||
+ | </pre> | ||
+ | |||
+ | |||
+ | === Radial field === | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Magnetic field | ||
+ | |[[Image:cdc_Bx_1350A.png|thumb|x400px|CDC space Bx (T) vs r and z for 1350A]] | ||
+ | |[[Image:cdc_Bz_1350A.png|thumb|x400px|CDC space Bz (T) vs r and z for 1350A]] | ||
+ | |} | ||
+ | |||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Drift time differences for Bx = -0.1T, -0.2T, -0.3T and various Bz | ||
+ | |[[Image:cdc_bx_timediff_0_1T.png|thumb|x190px|Drift time for Bx = -0.1T - drift time for Bx=0T]] | ||
+ | |[[Image:cdc_bx_timediff_0_2T.png|thumb|x190px|Drift time for Bx = -0.2T - drift time for Bx=0T]] | ||
+ | |[[Image:cdc_bx_timediff_0_3T.png|thumb|x190px|Drift time for Bx = -0.3T - drift time for Bx=0T]] | ||
+ | |} | ||
+ | |||
+ | === Drift time difference due to radial field === | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Drift time difference due to radial field | ||
+ | |[[Image:cdc_Bx_dtmap4.png|thumb|x400px|Max drift time difference (ns) due to Bx for 1350A]] | ||
+ | |} | ||
+ | |||
+ | 3ns diff at center occurs with Bz=1.65T, Bx=0.05T. Max drift times for Bz=1.65T are 665ns (Bx=0), 662ns (Bx=0.05T), 663ns (Bx=0.1T) | ||
+ | |||
+ | === Drift time correction for variation in Bz === | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Drift times from Garfield | ||
+ | |[[Image:cdc_bfield_dt_vs_t.png|thumb|x250px|Difference in drift time (for same driftradius) between higher field and 1.3T]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Linear fits to time differences | ||
+ | |[[Image:cdc_bfield_dt_fit_1_35.png|thumb|x190px|1.35T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_40.png|thumb|x190px|1.40T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_45.png|thumb|x190px|1.45T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_50.png|thumb|x190px|1.50T]] | ||
+ | |- | ||
+ | |[[Image:cdc_bfield_dt_fit_1_55.png|thumb|x190px|1.55T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_60.png|thumb|x190px|1.60T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_65.png|thumb|x190px|1.65T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_70.png|thumb|x190px|1.70T]] | ||
+ | |- | ||
+ | |[[Image:cdc_bfield_dt_fit_1_75.png|thumb|x190px|1.75T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_80.png|thumb|x190px|1.80T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_85.png|thumb|x190px|1.85T]] | ||
+ | |[[Image:cdc_bfield_dt_fit_1_90.png|thumb|x190px|1.90T]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Fit of param 0 values | ||
+ | |[[Image:cdc_bfield_dt_fitparam0.png|thumb|x250px|Fit of param 0 values]] | ||
+ | |} | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Global fit | ||
+ | |[[Image:cdc_bfield_dt_globalfit.png|thumb|x250px|One fit f(B,r) for all curves]] | ||
+ | |} | ||
+ | |||
+ | |||
+ | <pre> | ||
+ | Global fit function (-0.10 + 0.08 B)* drift-time | ||
+ | 1.35T Global fit av abs err 0.3ns max abs err 1.8ns at radius 7.6mm | ||
+ | 1.40T Global fit av abs err 0.3ns max abs err -0.7ns at radius 4.5mm | ||
1.45T Global fit av abs err 0.3ns max abs err 2.0ns at radius 7.6mm | 1.45T Global fit av abs err 0.3ns max abs err 2.0ns at radius 7.6mm | ||
1.50T Global fit av abs err 0.3ns max abs err 0.6ns at radius 2.4mm | 1.50T Global fit av abs err 0.3ns max abs err 0.6ns at radius 2.4mm | ||
− | 1.55T Global fit av abs err 0.4ns max abs err - | + | 1.55T Global fit av abs err 0.4ns max abs err -0.9ns at radius 6.3mm |
− | 1.60T Global fit av abs err 0. | + | 1.60T Global fit av abs err 0.4ns max abs err 1.1ns at radius 7.6mm |
− | 1.65T Global fit av abs err 0.4ns max abs err 1. | + | 1.65T Global fit av abs err 0.4ns max abs err 1.1ns at radius 7.0mm |
− | 1.70T Global fit av abs err 0. | + | 1.70T Global fit av abs err 0.5ns max abs err 1.6ns at radius 7.6mm |
− | 1.75T Global fit av abs err 0. | + | 1.75T Global fit av abs err 0.8ns max abs err -1.6ns at radius 7.6mm |
− | 1.80T Global fit av abs err 0. | + | 1.80T Global fit av abs err 0.6ns max abs err 1.6ns at radius 3.3mm |
− | 1.85T Global fit av abs err | + | 1.85T Global fit av abs err 1.1ns max abs err 2.1ns at radius 2.9mm |
− | 1.90T Global fit av abs err 0. | + | 1.90T Global fit av abs err 0.9ns max abs err 2.1ns at radius 3.5mm |
</pre> | </pre> | ||
+ | |||
+ | === The plan (1350A, 50% Ar, 50% CO2) === | ||
+ | |||
+ | # Measure drift time dt, radius r, axial z | ||
+ | # Look up B(r,z) | ||
+ | # Calculate drift time correction dtc = (-0.10 + 0.08 * B)* dt | ||
+ | # Find equivalent time for 1.3T field dte = dt - dtc | ||
+ | # Look up drift radius for time dte in Garfield table for 1.3T | ||
+ | # Ta-daa! | ||
+ | |||
+ | |||
+ | === 1200A, 50% Ar, 50% CO2 === | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Magnetic field maps of CDC region for solenoid current 1200A | ||
+ | |[[Image:cdc_Bz_1200.png|thumb|x400px|Axial field at 1200A]] | ||
+ | |[[Image:cdc_Br_1200.png|thumb|x400px|Radial field at 1200A]] | ||
+ | |} | ||
+ | |||
+ | New function is dtc = (-0.077 + 0.068 * B)* dt and base tables are for 1.1T.<br/> | ||
+ | |||
+ | # Measure drift time dt, radius r, axial z | ||
+ | # Look up B(r,z) | ||
+ | # Calculate drift time correction dtc = (-0.077 + 0.068 * B)* dt | ||
+ | # Find equivalent time for 1.1T field dte = dt - dtc | ||
+ | # Look up drift radius for time dte in Garfield table for 1.1T | ||
+ | |||
+ | <pre> | ||
+ | Global fit param 0 fit status 0 function -0.077 + 0.068x | ||
+ | 1.15T Global fit av abs err 0.2ns max abs err 0.7ns at radius 5.3mm | ||
+ | 1.20T Global fit av abs err 0.7ns max abs err 2.7ns at radius 7.6mm | ||
+ | 1.25T Global fit av abs err 0.5ns max abs err -1.7ns at radius 7.0mm | ||
+ | 1.30T Global fit av abs err 0.4ns max abs err 1.1ns at radius 4.1mm | ||
+ | 1.35T Global fit av abs err 0.5ns max abs err -1.0ns at radius 7.6mm | ||
+ | 1.40T Global fit av abs err 0.7ns max abs err -2.7ns at radius 7.4mm | ||
+ | 1.45T Global fit av abs err 0.7ns max abs err -1.7ns at radius 7.6mm | ||
+ | 1.50T Global fit av abs err 0.7ns max abs err -1.6ns at radius 7.0mm | ||
+ | 1.55T Global fit av abs err 1.0ns max abs err 1.9ns at radius 3.6mm | ||
+ | 1.60T Global fit av abs err 1.2ns max abs err -4.8ns at radius 7.6mm | ||
+ | 1.65T Global fit av abs err 1.1ns max abs err 2.1ns at radius 3.6mm | ||
+ | 1.70T Global fit av abs err 1.2ns max abs err -2.4ns at radius 7.6mm | ||
+ | 1.75T Global fit av abs err 1.5ns max abs err 2.9ns at radius 4.4mm | ||
+ | </pre> | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Global fit for 1200A, 50% Ar, 50% CO2 | ||
+ | |[[Image:cdc_bfield_dt_globalfit_1200_50_50.png|thumb|x250px|One fit f(B,r) for all curves]] | ||
+ | |} | ||
+ | |||
+ | |||
+ | === 1200A, 56% Ar, 44% CO2 === | ||
+ | |||
+ | New function is dtc = (-0.103 + 0.092 * B)* dt and base tables are for 1.1T.<br/> | ||
+ | |||
+ | # Measure drift time dt, radius r, axial z | ||
+ | # Look up B(r,z) | ||
+ | # Calculate drift time correction dtc = (-0.103 + 0.092 * B)* dt | ||
+ | # Find equivalent time for 1.1T field dte = dt - dtc | ||
+ | # Look up drift radius for time dte in Garfield table for 1.1T | ||
+ | |||
+ | <pre> | ||
+ | Global fit function errors: | ||
+ | Global fit param 0 fit status 0 function -0.103 + 0.092x | ||
+ | 1.15T Global fit av abs err 0.7ns max abs err 2.7ns at radius 7.6mm | ||
+ | 1.20T Global fit av abs err 0.3ns max abs err -1.0ns at radius 7.3mm | ||
+ | 1.25T Global fit av abs err 0.2ns max abs err 0.5ns at radius 3.8mm | ||
+ | 1.30T Global fit av abs err 0.4ns max abs err -1.4ns at radius 7.2mm | ||
+ | 1.35T Global fit av abs err 0.6ns max abs err -1.9ns at radius 7.6mm | ||
+ | 1.40T Global fit av abs err 0.4ns max abs err 1.1ns at radius 3.2mm | ||
+ | 1.45T Global fit av abs err 0.7ns max abs err -2.0ns at radius 6.9mm | ||
+ | 1.50T Global fit av abs err 0.8ns max abs err -2.2ns at radius 7.6mm | ||
+ | 1.55T Global fit av abs err 0.7ns max abs err 1.4ns at radius 3.4mm | ||
+ | 1.60T Global fit av abs err 0.8ns max abs err -1.7ns at radius 6.9mm | ||
+ | 1.65T Global fit av abs err 1.2ns max abs err 2.5ns at radius 4.3mm | ||
+ | 1.70T Global fit av abs err 1.1ns max abs err 2.4ns at radius 3.8mm | ||
+ | 1.75T Global fit av abs err 1.3ns max abs err -3.2ns at radius 7.6mm | ||
+ | </pre> | ||
+ | |||
+ | {| border="0" cellpadding="2" | ||
+ | |+Global fit for 1200A, 56% Ar, 44% CO2 | ||
+ | |[[Image:cdc_bfield_dt_globalfit_1200_56_44.png|thumb|x250px|One fit f(B,r) for all curves]] | ||
+ | |} |
Latest revision as of 09:49, 7 April 2015
Contents
- 1 B field in CDC volume
- 2 Finding drift time difference wrt 1.3T
- 3 Finding drift time difference wrt 1.3T (more data)
- 4 Finding B
- 5 B fits (part 2)
- 6 B fits (part 3)
- 7 Radial field
- 8 Drift time difference due to radial field
- 9 Drift time correction for variation in Bz
- 10 The plan (1350A, 50% Ar, 50% CO2)
- 11 1200A, 50% Ar, 50% CO2
- 12 1200A, 56% Ar, 44% CO2
B field in CDC volume
CDC volume occupies z=17 to 167 and r=10 to 60
At large drift radius, change in B field of 0.1T causes change in drift time of ~ 5ns. This would cause approx 30um error in drift radius. Resolution in 0T is 80um-100um at large drift radius.
Plan:
- measure drift time, r, z
- find B field at r, z
- use time-diff function(B) to calculate what drift time would be in 1.3T
- use garfield tables for 1.3T to find drift distance
Finding drift time difference wrt 1.3T
Max error between fitted function and time difference is 5.5ns at radius 7.6mm in 1.5T dataset
Finding drift time difference wrt 1.3T (more data)
Repeated the above with more datasets
Max error using one fit f(B,driftradius) for all curves is 2.4ns at 7.6mm for 1.7T
Param 0 fit status 0 function 10.72 + -19.17x + 10.68x^2 + -1.74x^3 Param 1 fit status 0 function 9.88 + -8.22x + 2.06x^2 1.35T Global fit av abs err 0.4ns max abs err -1.2ns at radius 5.8mm 1.40T Global fit av abs err 0.2ns max abs err -0.5ns at radius 4.5mm 1.45T Global fit av abs err 0.3ns max abs err 2.0ns at radius 7.6mm 1.50T Global fit av abs err 0.3ns max abs err 0.6ns at radius 2.4mm 1.55T Global fit av abs err 0.4ns max abs err -1.1ns at radius 6.4mm 1.60T Global fit av abs err 0.3ns max abs err 1.1ns at radius 7.6mm 1.65T Global fit av abs err 0.4ns max abs err 1.2ns at radius 7.0mm 1.70T Global fit av abs err 0.4ns max abs err 2.4ns at radius 7.6mm 1.75T Global fit av abs err 0.5ns max abs err -1.1ns at radius 6.1mm 1.80T Global fit av abs err 0.3ns max abs err 1.4ns at radius 7.6mm 1.85T Global fit av abs err 0.6ns max abs err -1.4ns at radius 6.2mm 1.90T Global fit av abs err 0.3ns max abs err 1.2ns at radius 7.6mm
Finding B
Mostly backward angle z (16 to 86cm) pol3 fit parameters radius (cm) [0] [1] z [2] z^2 [3] z^3 10 1.049e+00 1.597e-02 -1.419e-04 5.523e-07 20 1.062e+00 1.660e-02 -1.616e-04 6.818e-07 30 1.084e+00 1.771e-02 -1.979e-04 9.257e-07 35 1.100e+00 1.851e-02 -2.252e-04 1.112e-06 40 1.119e+00 1.951e-02 -2.596e-04 1.352e-06 45 1.142e+00 2.077e-02 -3.040e-04 1.664e-06 50 1.169e+00 2.238e-02 -3.613e-04 2.072e-06 53 1.200e+00 2.448e-02 -4.361e-04 2.607e-06 57 1.217e+00 2.572e-02 -4.800e-04 2.922e-06 Mostly forward angle z (84 to 168cm) pol3 fit parameters radius (cm) [0] [1] z [2] z^2 [3] z^3 10 1.200e+00 8.748e-03 -3.533e-05 5.208e-08 20 1.184e+00 9.145e-03 -3.786e-05 5.646e-08 30 1.126e+00 1.056e-02 -4.785e-05 7.824e-08 35 1.059e+00 1.216e-02 -5.957e-05 1.055e-07 40 9.507e-01 1.477e-02 -7.907e-05 1.524e-07 45 7.758e-01 1.898e-02 -1.111e-04 2.311e-07 50 4.998e-01 2.567e-02 -1.627e-04 3.600e-07 53 7.226e-02 3.611e-02 -2.441e-04 5.662e-07 57 -2.104e-01 4.304e-02 -2.987e-04 7.055e-07
Fits to fit parameters above as function of r
Mostly backward angle z (16 to 86cm) pol3 fit parameter fit parameters as function of r param [0] [1] r [2] r^2 [3] r^3 0 1.045e+00 2.290e-04 1.972e-05 5.338e-07 1 1.518e-02 1.009e-04 -2.982e-06 7.920e-08 2 -1.157e-04 -3.444e-06 1.134e-07 -2.934e-09 3 3.657e-07 2.528e-08 -8.925e-10 2.196e-11 Mostly forward angle z (84 to 168cm) pol3 fit parameter fit parameters as function of r param [0] [1] r [2] r^2 [3] r^3 0 1.451e+00 -4.060e-02 1.916e-03 -3.024e-05 1 2.443e-03 1.017e-03 -4.767e-05 7.464e-07 2 1.570e-05 -8.240e-06 3.862e-07 -5.966e-09 3 -8.188e-08 2.168e-08 -1.019e-09 1.552e-11 Double_t lo_z0[4]={1.045E+00, 2.290E-04, 1.972E-05, 5.338E-07}; Double_t lo_z1[4]={1.518E-02, 1.009E-04, -2.982E-06, 7.920E-08}; Double_t lo_z2[4]={-1.157E-04, -3.444E-06, 1.134E-07, -2.934E-09}; Double_t lo_z3[4]={3.657E-07, 2.528E-08, -8.925E-10, 2.196E-11}; Double_t hi_z0[4]={1.451E+00, -4.060E-02, 1.916E-03, -3.024E-05}; Double_t hi_z1[4]={2.443E-03, 1.017E-03, -4.767E-05, 7.464E-07}; Double_t hi_z2[4]={1.570E-05, -8.240E-06, 3.862E-07, -5.966E-09}; Double_t hi_z3[4]={-8.188E-08, 2.168E-08, -1.019E-09, 1.552E-11};
so for z<85cm, B ~ (1.045 + 2.290e-04r + 1.972e-05r^2 + 5.338e-07r^3) + (1.518e-02 + 1.009e-04r - 2.982e-06r^2 + 7.920e-08r^3)z + (-1.157e-04 - 3.444e-06r + 1.134e-07r^2 - 2.934e-09r^3)z^2 + ( 3.657e-07 + 2.528e-08r - 8.925e-10r^2 + 2.196e-11r^3)z^3 with B in T and z and r in cm
Max difference between fit and original histogram from bfield.root is 0.05T.
The above was for Btot. Repeated for Bz, params below. Max err is 0.0449T at z=22cm, r=55cm Fits are in bfield/bfits
Mostly backward angle z (16 to 86cm) pol3 fit parameters r (cm) [0] [1] z [2] z^2 [3] z^3 10.0 1.046e+00 1.608e-02 -1.432e-04 5.577e-07 20.0 1.051e+00 1.699e-02 -1.666e-04 7.033e-07 30.0 1.059e+00 1.865e-02 -2.103e-04 9.811e-07 35.0 1.065e+00 1.987e-02 -2.435e-04 1.196e-06 40.0 1.072e+00 2.139e-02 -2.856e-04 1.472e-06 45.0 1.081e+00 2.330e-02 -3.399e-04 1.833e-06 50.0 1.091e+00 2.573e-02 -4.099e-04 2.304e-06 53.0 1.103e+00 2.880e-02 -5.001e-04 2.916e-06 57.0 1.109e+00 3.057e-02 -5.527e-04 3.274e-06 Mostly forward angle z (84 to 168cm) pol3 fit parameters r (cm) [0] [1] z [2] z^2 [3] z^3 10.0 1.200e+00 8.750e-03 -3.532e-05 5.202e-08 20.0 1.183e+00 9.148e-03 -3.779e-05 5.612e-08 30.0 1.124e+00 1.055e-02 -4.755e-05 7.711e-08 35.0 1.057e+00 1.213e-02 -5.904e-05 1.037e-07 40.0 9.488e-01 1.472e-02 -7.826e-05 1.496e-07 45.0 7.734e-01 1.892e-02 -1.100e-04 2.272e-07 50.0 4.956e-01 2.561e-02 -1.613e-04 3.550e-07 53.0 6.308e-02 3.610e-02 -2.429e-04 5.608e-07 57.0 -2.243e-01 4.310e-02 -2.977e-04 7.005e-07 Mostly backward angle z (16 to 86cm) pol3 fit parameter fit parameters as function of r param [0] [1] r [2] r^2 [3] r^3 0 1.045e+00 1.338e-04 4.523e-06 2.368e-07 1 1.513e-02 1.160e-04 -3.189e-06 1.053e-07 2 -1.146e-04 -3.751e-06 1.237e-07 -3.429e-09 3 3.595e-07 2.689e-08 -9.586e-10 2.461e-11 Mostly forward angle z (84 to 168cm) pol3 fit parameter fit parameters as function of r param [0] [1] r [2] r^2 [3] r^3 0 1.459e+00 -4.190e-02 1.966e-03 -3.085e-05 1 2.268e-03 1.043e-03 -4.871e-05 7.577e-07 2 1.686e-05 -8.411e-06 3.930e-07 -6.034e-09 3 -8.439e-08 2.205e-08 -1.033e-09 1.565e-11 Double_t lo_z0[4]={1.0445E+00, 1.3383E-04, 4.5227E-06, 2.3683E-07}; Double_t lo_z1[4]={1.5132E-02, 1.1596E-04, -3.1894E-06, 1.0531E-07}; Double_t lo_z2[4]={-1.1458E-04, -3.7512E-06, 1.2367E-07, -3.4290E-09}; Double_t lo_z3[4]={3.5954E-07, 2.6887E-08, -9.5863E-10, 2.4608E-11}; Double_t hi_z0[4]={1.4594E+00, -4.1900E-02, 1.9661E-03, -3.0849E-05}; Double_t hi_z1[4]={2.2682E-03, 1.0434E-03, -4.8708E-05, 7.5771E-07}; Double_t hi_z2[4]={1.6860E-05, -8.4111E-06, 3.9302E-07, -6.0335E-09}; Double_t hi_z3[4]={-8.4393E-08, 2.2049E-08, -1.0333E-09, 1.5649E-11};
B fits (part 2)
Bz fit parameters
/* Max err 0.067 at z 22.0 r 55.0 */ Double_t lo_z0[4]={1.0440E+00, 2.8431E-04, 2.2198E-05, 4.4318E-07}; Double_t lo_z1[4]={1.5235E-02, 8.4686E-05, -2.0227E-06, 6.4062E-08}; Double_t lo_z2[4]={-1.1729E-04, -2.9396E-06, 8.2057E-08, -2.4271E-09}; Double_t lo_z3[4]={3.7660E-07, 2.1683E-08, -6.6856E-10, 1.8324E-11}; Double_t hi_z0[4]={1.4139E+00, -3.3108E-02, 1.5528E-03, -2.5238E-05}; Double_t hi_z1[4]={3.3934E-03, 8.2637E-04, -3.8508E-05, 6.2072E-07}; Double_t hi_z2[4]={7.7830E-06, -6.6706E-06, 3.1153E-07, -4.9497E-09}; Double_t hi_z3[4]={-6.0532E-08, 1.7501E-08, -8.2145E-10, 1.2855E-11};
B fits (part 3)
Used pol2 fit for r<41cm and pairs of pol3 fits (as before) for r>=41cm
Used pol2 fit for r<=30cm and pairs of pol3 fits (as before) for r>30cm
Used pol2 fit for r<=30cm and pairs of pol2 fits for r>30cm
Used pol2 fit for r<=25cm and pairs of pol2 fits for r>25cm
Fits to fit parameters above
This gives results comparable to the original fits with slightly simpler polynomials
Double_t lo_r0[3]={1.1746E+00, -1.1176E-04, 9.1798E-05}; Double_t lo_r1[3]={8.6423E-03, 4.5702E-06, -1.4617E-06}; Double_t lo_r2[3]={-2.5830E-05, -2.6167E-08, 5.9721E-09}; Double_t lo_z0[3]={1.0753E+00, -2.0457E-03, 7.8510E-05}; Double_t lo_z1[3]={2.0352E-02, -2.7862E-04, 6.4074E-06}; Double_t lo_z2[3]={-3.0851E-04, 1.0707E-05, -2.3605E-07}; Double_t lo_z3[3]={1.8161E-06, -8.1165E-08, 1.7311E-09}; Double_t hi_z0[3]={-6.1452E-01, 1.1061E-01, -1.7751E-03}; Double_t hi_z1[3]={5.3402E-02, -2.7137E-03, 4.3399E-05}; Double_t hi_z2[3]={-3.9186E-04, 2.1597E-05, -3.4203E-07}; Double_t hi_z3[3]={9.7935E-07, -5.6001E-08, 8.7685E-10}; if (r<=25) { bfit0 = lo_r0[0] + lo_r0[1]*r + lo_r0[2]*r*r; bfit1 = lo_r1[0] + lo_r1[1]*r + lo_r1[2]*r*r; bfit2 = lo_r2[0] + lo_r2[1]*r + lo_r2[2]*r*r; bfit3 = 0; } else { if (z<85) { bfit0 = lo_z0[0] + lo_z0[1]*r + lo_z0[2]*r*r; bfit1 = lo_z1[0] + lo_z1[1]*r + lo_z1[2]*r*r; bfit2 = lo_z2[0] + lo_z2[1]*r + lo_z2[2]*r*r; bfit3 = lo_z3[0] + lo_z3[1]*r + lo_z3[2]*r*r; } else { bfit0 = hi_z0[0] + hi_z0[1]*r + hi_z0[2]*r*r; bfit1 = hi_z1[0] + hi_z1[1]*r + hi_z1[2]*r*r; bfit2 = hi_z2[0] + hi_z2[1]*r + hi_z2[2]*r*r; bfit3 = hi_z3[0] + hi_z3[1]*r + hi_z3[2]*r*r; } } bfit = bfit0 + bfit1*z + bfit2*z*z + bfit3*z*z*z;
Radial field
Drift time difference due to radial field
3ns diff at center occurs with Bz=1.65T, Bx=0.05T. Max drift times for Bz=1.65T are 665ns (Bx=0), 662ns (Bx=0.05T), 663ns (Bx=0.1T)
Drift time correction for variation in Bz
Global fit function (-0.10 + 0.08 B)* drift-time 1.35T Global fit av abs err 0.3ns max abs err 1.8ns at radius 7.6mm 1.40T Global fit av abs err 0.3ns max abs err -0.7ns at radius 4.5mm 1.45T Global fit av abs err 0.3ns max abs err 2.0ns at radius 7.6mm 1.50T Global fit av abs err 0.3ns max abs err 0.6ns at radius 2.4mm 1.55T Global fit av abs err 0.4ns max abs err -0.9ns at radius 6.3mm 1.60T Global fit av abs err 0.4ns max abs err 1.1ns at radius 7.6mm 1.65T Global fit av abs err 0.4ns max abs err 1.1ns at radius 7.0mm 1.70T Global fit av abs err 0.5ns max abs err 1.6ns at radius 7.6mm 1.75T Global fit av abs err 0.8ns max abs err -1.6ns at radius 7.6mm 1.80T Global fit av abs err 0.6ns max abs err 1.6ns at radius 3.3mm 1.85T Global fit av abs err 1.1ns max abs err 2.1ns at radius 2.9mm 1.90T Global fit av abs err 0.9ns max abs err 2.1ns at radius 3.5mm
The plan (1350A, 50% Ar, 50% CO2)
- Measure drift time dt, radius r, axial z
- Look up B(r,z)
- Calculate drift time correction dtc = (-0.10 + 0.08 * B)* dt
- Find equivalent time for 1.3T field dte = dt - dtc
- Look up drift radius for time dte in Garfield table for 1.3T
- Ta-daa!
1200A, 50% Ar, 50% CO2
New function is dtc = (-0.077 + 0.068 * B)* dt and base tables are for 1.1T.
- Measure drift time dt, radius r, axial z
- Look up B(r,z)
- Calculate drift time correction dtc = (-0.077 + 0.068 * B)* dt
- Find equivalent time for 1.1T field dte = dt - dtc
- Look up drift radius for time dte in Garfield table for 1.1T
Global fit param 0 fit status 0 function -0.077 + 0.068x 1.15T Global fit av abs err 0.2ns max abs err 0.7ns at radius 5.3mm 1.20T Global fit av abs err 0.7ns max abs err 2.7ns at radius 7.6mm 1.25T Global fit av abs err 0.5ns max abs err -1.7ns at radius 7.0mm 1.30T Global fit av abs err 0.4ns max abs err 1.1ns at radius 4.1mm 1.35T Global fit av abs err 0.5ns max abs err -1.0ns at radius 7.6mm 1.40T Global fit av abs err 0.7ns max abs err -2.7ns at radius 7.4mm 1.45T Global fit av abs err 0.7ns max abs err -1.7ns at radius 7.6mm 1.50T Global fit av abs err 0.7ns max abs err -1.6ns at radius 7.0mm 1.55T Global fit av abs err 1.0ns max abs err 1.9ns at radius 3.6mm 1.60T Global fit av abs err 1.2ns max abs err -4.8ns at radius 7.6mm 1.65T Global fit av abs err 1.1ns max abs err 2.1ns at radius 3.6mm 1.70T Global fit av abs err 1.2ns max abs err -2.4ns at radius 7.6mm 1.75T Global fit av abs err 1.5ns max abs err 2.9ns at radius 4.4mm
1200A, 56% Ar, 44% CO2
New function is dtc = (-0.103 + 0.092 * B)* dt and base tables are for 1.1T.
- Measure drift time dt, radius r, axial z
- Look up B(r,z)
- Calculate drift time correction dtc = (-0.103 + 0.092 * B)* dt
- Find equivalent time for 1.1T field dte = dt - dtc
- Look up drift radius for time dte in Garfield table for 1.1T
Global fit function errors: Global fit param 0 fit status 0 function -0.103 + 0.092x 1.15T Global fit av abs err 0.7ns max abs err 2.7ns at radius 7.6mm 1.20T Global fit av abs err 0.3ns max abs err -1.0ns at radius 7.3mm 1.25T Global fit av abs err 0.2ns max abs err 0.5ns at radius 3.8mm 1.30T Global fit av abs err 0.4ns max abs err -1.4ns at radius 7.2mm 1.35T Global fit av abs err 0.6ns max abs err -1.9ns at radius 7.6mm 1.40T Global fit av abs err 0.4ns max abs err 1.1ns at radius 3.2mm 1.45T Global fit av abs err 0.7ns max abs err -2.0ns at radius 6.9mm 1.50T Global fit av abs err 0.8ns max abs err -2.2ns at radius 7.6mm 1.55T Global fit av abs err 0.7ns max abs err 1.4ns at radius 3.4mm 1.60T Global fit av abs err 0.8ns max abs err -1.7ns at radius 6.9mm 1.65T Global fit av abs err 1.2ns max abs err 2.5ns at radius 4.3mm 1.70T Global fit av abs err 1.1ns max abs err 2.4ns at radius 3.8mm 1.75T Global fit av abs err 1.3ns max abs err -3.2ns at radius 7.6mm