## Kaon LT Meetings » LTsep_mtg_22aug26.txt

1 |
LT Separation Meeting - 26/08/22 |
---|---|

2 | |

3 |
Stephen Kay |

4 |
Garth Huber |

5 |
Vijay Kumar |

6 |
Ali Usman |

7 |
Muhammad Junaid |

8 | |

9 |
Garth has a set of slides - https://redmine.jlab.org/attachments/download/1593/LT_sep_iterations.pdf |

10 |
See also Bill's slides - https://redmine.jlab.org/attachments/download/1594/Kaon_LT_tutorial_18nov28.pdf |

11 |
Start from Bill's code - https://github.com/billlee77/omega_analysis |

12 | |

13 |
- Prior to LT separation, need *final* normalised yields |

14 |
- Most of the work involved is actually in getting to this point! |

15 |
- If things are converging, iteration steps should be quick |

16 |
- This includes diamond cut already |

17 |
- One could analyse the variations across only the larger, high epsilon diamond |

18 |
- All efficincies included |

19 |
- LT |

20 |
- FADC DT |

21 |
- Cryotarget |

22 |
- Yield corrections |

23 |
- ALL tested for reliablity over a wide range and applied |

24 |
- All kinematic offsets determined and finalised |

25 |
- Need one thing kept constant between iterations |

26 |
- The normalised yields and distributions |

27 |
- Cross section varies across experimental acceptance |

28 |
- Need to choose a functional form that will reasonably account for this |

29 |
- Don't know this in advance |

30 |
- Make a choice, start iteration process with it |

31 |
- Can check previous analyses for guidance on forms to try |

32 |
- Can be different for each diamond plot |

33 |
- So long as it it well understood for each individual diamond |

34 |
- On diamond plot, t bins slice across diagonally |

35 |
- W and Q2 different for each t bin |

36 |
- Results quote different Q2/W value for each t bin |

37 |
- t min varies across the diamond |

38 |
- t min varies even within a t bin |

39 |
- In SIMC, replace physics_pion.f with physics_iterate.f |

40 |
- Need to make some good guesses of intial paramter values for your first iteration |

41 |
- Need some Q2 dependence, t dependence etc |

42 |
- Will likely have to go back and change these functions a lot |

43 |
- Each Q2/W should be done separately |

44 |
- Can't expect the procedure to work globally |

45 |
- *Keep notes organised!* |

46 |
- *Keep all output!* |

47 |
- Store each iteration in its own directory |

48 |
- Step 1 - SIMC distributions |

49 |
- Run SIMC for large #events |

50 |
- Generate spectrometer and physics variables using functional form and fitpar |

51 |
- Do this setting by setting for a given Q2/W diamond |

52 |
- E.g. Left/Centre/Right all different |

53 |
- Compare to data |

54 |
- Step 2 - Combine Left/Centre/Right settings at high and low epsilon for each W/Q2/t/phi/epsilon bin for data and MC |

55 |
- Get statistical errors for each bin |

56 |
- Get yield, error in yield for data and simulation |

57 |
- Calculate ratio and error in ratio for each bin |

58 |
- Process enough events to minimise SIMC statistical error! |

59 |
- Shouldn't be dominated by simulation statistical error! |

60 |
- Step 3 - Calculate average kinematics |

61 |
- Mean data values of W/Q2/Theta/epsilon for each t bin at high an low epsilon |

62 |
- Need this for data and MC |

63 |
- Values will differ between high and low epsilon slghtly |

64 |
- Will also change slightly as the model is iterated |

65 |
- Step 4a - Inspect and understand data closely |

66 |
- Examine variation on phi between data and MC within each t bin |

67 |
- Deviations between data/MC are usually indicated as wiggles in the ratio |

68 |
- Want R ~ 1 across a broad kinematic range |

69 |
- Some plots earlier on might look ok, only by examining the ratio very closely in individual bins can some differences be found |

70 |
- How things vary should guide what changes in the next iteration |

71 |
- E.g. Is it the interference terms that are off? |

72 |
- Step 4b - Even closer examination |

73 |
- Each t bin (for one epsilon), binned in theta |

74 |
- Interference terms depend upon theta |

75 |
- Interference terms vanish in parallel kinematics |

76 |
- Explicitly chose a functional form that makes this true |

77 |
- In lowest bin, should have minimal interference terms |

78 |
- Wiggles get larger as theta increases |

79 |
- Phi distributions for each t bin, subdivided into 8 theta bins |

80 |
- A lot of plots |

81 |
- More than you probably want to put in thesis |

82 |
- BUT a tech note/report is encouraged |

83 |
- Even here, you probably want to show only a few plots as an example |

84 |
- Step 5 - Calculate unseparated cross section |

85 |
- Using fitpar for the iteration, evaluate the model at *average* kinematics of the data for each t-bin |

86 |
- Calculate using formula and fit parameters, calculate the cross section at the average values of W/Q2/Theta and epsilon for each bin |

87 |
- Read over Blok PRC 78 (2008) 045202 carefully |

88 |
- Each t/phi bin at both high and low epsilon |

89 |
- The weight in SIMC is NOT the unseparated two fold cross section |

90 |
- It is the five fold cross section (d5sigma) |

91 |
- Fitting procedure iterated until the experimental cross section changed by less than a prescribed amount (~1%) |

92 |
- Step 6 - Fit Rosenbluth Eqn |

93 |
- Each t bin is fit separately |

94 |
- Fit result gives L, T, LT, TT cross sections for each t bin |

95 |
- Step 7 - Fit L, T, LT, TT values to get new fitpar |

96 |
- Compare t bins with each other |

97 |
- Fit with physics_iterate.f functions to give next iteration model parameters |

98 |
- Now, we're back at the top. Feed these back into step 1 |

99 |
- Repeat steps 1-7 until separated cros sections are stable |

100 |
- *Don't re-run SIMC* |

101 |
- Recalculate weight for each event |

102 |
- Any theses older than Tanja's are likely using different procedures |

103 |
- Tanja's procedures were improved on earlier ones |

104 |
- Some slight modifications after this |

105 |
- Only Marco and Bill used Tanja's procedure |

106 |
- Evaluating if fit equations are ok |

107 |
- Procedure usually works ok, but for some kinematics, pi-/pi+ sigma wouldn't converge |

108 |
- Compare fitpar from different Q2/W and see if they were slowly varying |

109 |
- If not, could use their variation as a pointer to an alternative functional form variety to try |

110 |
- Things will behave in a consistent manner if the correct solution is found |

111 |
- General pointers |

112 |
- Read Blok paper carefully |

113 |
- Review Bill's slides |

114 |
- A single Q2/W iteration should only take 1-2 hours |

115 |
- Fitpar should converge in a few iterations to give 0.5<R<2 |

116 |
- The main work is getting R to be acceptably flat |