# Basic ideas of follows convention of Hall A Tech Notes JLAB-TN-02-012: # http://hallaweb.jlab.org/publications/Technotes/files/2002/02-012.pdf ######################################################### # Common Section for THaOptics ######################################################### # source of beam information. # this item could be overwriten by # 1. parameter in the constructor function # 2. any call of function THaOptics::SetBeamName BeamDetectorName = B # OLD INPUTS: # shift of target coordinate system origin, in the hall coordinate system # rotation of the whole spectrometer will be read out from db_run.dat # target_x = 0.0 # target_y = 0.0 # target_z = 0.0 # the shift of Transport Coordinate System Origin in Target Coordinate System trans_x = -0.37 trans_y = 0.0 trans_z = 2.4 # the rotation for a vector expressed in Transport Coordinate System # to be translate into Target Coordinate System # # The rotation is described by a series of angles ( unit is degree ) as below. # # The rotation take the order of x -> y-> z trans_rot_x = -0.0 trans_rot_y = -25.0 trans_rot_z = 0.0 #0st order corection ThShift = -0.04945 PhShift = 0.0 PShift = 0.0 PathLenShift= 0.0 ######################################################### # Additional Section for THaOpticsAGen ######################################################### # description of this part could be found at # http://hallaweb.jlab.org/experiment/E02-013/offline/BBdatabase.txt effbendplanetheta = 10.0 effbendz = 1.85 effbendx = 0.0 effbendy = 0.0 ciDef = 0.283 ciDefXbend = 0.126 cThetaTarg = 0.0 cY = 0.10 cPhi = 0.0 ca = 0.0 cV = 1.0 cVX = 0.012 cVXp = 0.0 cVY = 0.05 cVYp = 0.0 ######################################################### # Additional Section for THaOpticsHRS ######################################################### # format is described in # http://hallaweb.jlab.org/publications/Technotes/files/2002/02-012.pdf [matrix] t 1 0 0 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 y 1 0 0 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 p 1 0 0 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 D 0 0 0 -0.08749 -0.67627 1.91361 -3.11562 -6.33754 11.5908 D 1 0 0 2.80163 -8.35861 17.7113 17.3145 -56.5428 0.00000 D 2 0 0 11.6524 -46.0247 7.91907 119.282 0.00000 0.00000 D 3 0 0 44.3530 -77.9449 -99.1704 0.00000 0.00000 0.00000 D 4 0 0 79.0134 -4.96614 0.00000 0.00000 0.00000 0.00000 D 5 0 0 42.3238 0.00000 0.00000 0.00000 0.00000 0.00000 D 0 0 0 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 T 0 0 0 0.01714 0.55084 -0.06026 -0.07013 0.04245 T 0 0 1 -0.00797 -0.23834 0.09164 -0.14835 0.00000 T 0 0 2 -0.07563 0.00000 0.00000 0.00000 0.00000 T 0 1 0 0.00721 0.08022 -0.13420 0.20116 0.00000 T 0 1 1 0.06494 0.00000 0.27701 0.00000 0.00000 T 0 2 0 0.05660 0.00000 -0.38301 0.00000 0.00000 T 0 2 1 -0.35652 0.00000 0.00000 0.00000 0.00000 T 0 3 0 0.21441 0.00000 0.00000 0.00000 0.00000 T 1 0 0 -0.48367 -0.04962 0.17132 0.15176 0.00000 T 1 0 1 0.19989 0.00000 0.70217 0.00000 0.00000 T 1 1 0 -0.06061 0.00000 -0.38277 0.00000 0.00000 T 1 3 0 0.28484 0.00000 0.00000 0.00000 0.00000 T 2 0 0 0.04598 -0.14871 -0.41783 0.00000 0.00000 T 2 0 1 0.00000 -1.38744 0.00000 0.00000 0.00000 T 2 1 0 0.00000 0.26537 0.00000 0.00000 0.00000 T 3 0 0 0.11351 0.00000 0.00000 0.00000 0.00000 T 3 0 1 0.49991 0.00000 0.00000 0.00000 0.00000 T 4 0 0 0.28425 0.00000 0.00000 0.00000 0.00000 P 0 0 0 0.00458 0.00904 -0.03227 0.00085 0.00000 P 0 0 1 0.95166 -0.12289 1.09276 -1.39597 0.00000 P 0 0 2 0.40387 -1.40952 -1.16763 0.00000 0.00000 P 0 0 3 0.00000 -6.86876 0.00000 0.00000 0.00000 P 0 1 0 0.10867 -0.16805 -1.53930 2.06102 0.00000 P 0 1 1 -0.56894 3.47503 3.62306 0.00000 0.00000 P 0 1 2 0.00000 -15.6968 0.00000 0.00000 0.00000 P 0 2 0 0.27395 -1.93917 -1.50013 0.00000 0.00000 P 0 2 1 0.00000 53.7498 0.00000 0.00000 0.00000 P 0 3 0 1.13541 -20.6053 0.00000 0.00000 0.00000 P 1 0 0 -0.04127 -0.00523 0.33872 0.00000 0.00000 P 1 0 1 0.90995 -0.41927 -10.9548 0.00000 0.00000 P 1 0 2 0.00000 12.0085 0.00000 0.00000 0.00000 P 1 1 0 -0.89483 0.00000 13.8273 0.00000 0.00000 P 1 1 1 -0.22486 -16.9376 0.00000 0.00000 0.00000 P 1 2 0 0.43072 4.53636 0.00000 0.00000 0.00000 P 1 2 1 -28.2799 0.00000 0.00000 0.00000 0.00000 P 2 0 0 0.01202 -0.25563 0.00000 0.00000 0.00000 P 2 0 1 0.17365 13.0352 0.00000 0.00000 0.00000 P 2 0 2 -6.23897 0.00000 0.00000 0.00000 0.00000 P 2 1 0 -0.54670 -13.7950 0.00000 0.00000 0.00000 P 2 1 1 15.9660 0.00000 0.00000 0.00000 0.00000 P 2 2 0 -8.13482 0.00000 0.00000 0.00000 0.00000 P 3 0 1 -6.19519 0.00000 0.00000 0.00000 0.00000 P 3 1 0 5.77147 0.00000 0.00000 0.00000 0.00000 Y 0 0 0 0.02318 0.00000 0.03641 -0.13060 Y 0 0 1 -2.76580 -0.24199 -1.03855 5.63012 Y 0 0 2 -2.33569 4.94169 0.00000 0.00000 Y 0 0 3 19.9915 -67.6181 0.00000 0.00000 Y 0 1 0 0.91701 0.35879 1.91131 -8.24815 Y 0 1 1 1.98481 -9.50365 0.00000 0.00000 Y 0 1 2 -30.8046 177.272 0.00000 0.00000 Y 0 2 0 -0.63517 4.41392 -0.22726 0.00000 Y 0 2 1 20.2199 -168.771 0.00000 0.00000 Y 0 3 0 -6.60272 49.9945 0.00000 0.00000 Y 1 0 0 0.03889 0.00000 0.00000 0.00000 Y 1 0 1 -0.61325 0.00000 0.24360 0.00000 Y 1 1 0 0.87685 0.00000 0.00000 0.00000 Y 1 1 1 2.25219 0.00000 0.00000 0.00000 Y 1 1 2 -28.7562 0.00000 0.00000 0.00000 Y 1 2 0 -1.35994 0.00000 0.00000 0.00000 Y 1 2 1 47.3592 0.00000 0.00000 0.00000 Y 2 0 1 -0.87017 0.00000 0.00000 0.00000 Y 2 1 0 1.61122 0.00000 0.00000 0.00000 Y 2 2 0 4.62632 0.00000 0.00000 0.00000 ######################################################### # Additional Section for THaOpticsAnalytical ######################################################### # P=q*B*R=K*I*R, in which, # P is calculated momentum # q is charge of the incoming particle # B is the magnetic field # R is bending radius # K is an constant related to the magnetic # I is the current in coil # # there are two method of scaling bending radius to momentum: # 1. given current I and constant K # 2. given radius of central momentum # this class will try both of above method. # Priority of method 1 is higher # K<=0 will a push this class to use method 2 # setup the constant K here, as discussed above # units: GeV/A/m/c K=7.96996e-4 # The radius of curve for a particle with central momentum to move in the average magnetic field. # CurveRadiusOfCentralMomentum>0 in meter to be valid. # in case CurveRadiusOfCentralMomentum<=0, code will be automatically calculated one according to definition. CurveRadiusOfCentralMomentum= -1 # see http://www.jlab.org/~jinhuang/BigBiteDoc/THaOpticsAnalytical.html for reference # no need to change if old parameters working fine # model defined fitted parameters MagneticFieldFrontShift = 0.0 MagneticFieldTopCrossing = 1.773 MagneticFieldRearSlopeAngle = 20 #CurveRadiusOfCentralMomentum = 1.57 ######################################################### # Additional Section for THaOpticsNoFeild #########################################################