Meeting No. 1

Converting measured CS offsets to FF:

To be able to calculate FFs from the measured CS differences in the tail region, we need a transformation function, which tells us, how differences in CS reflect differences in the FF. The corresponding function for all three energies is shown below:

01.)

To get this function I considered our simulation. For each kin. setting I ran simulation 5 times, each with different scaling of the form-factor (-2%,-1%, 0%, +1% +2%). However, when scaled, both ISR and FSR get scaled. Hence, to get the effect of the ISR only, I need to compare ISR(scaled) + FSR(unscaled) with the ISR(unscaled)+FSR(unscaled). This can be done, by looking at the simulated Vertex-Q2 histograms, fitting the peaks to identify the ISR and FSR parts and then do the comparison. The plot below shows the spectrum for 0% simulation on the left and 2% percent simulation on the right.

02.)

Analysis of Snout contribution

Following three plots show the size of the snout effect. This contribution needs to be subtracted from the measured CS ratios, in order to get proper results. The corrections are small for the 495MeV setting but very large for the 195MeV setting.

03.) 04.) 05.)

The size of the correction was determined by estimating how much snout contributes to the measured spectra, by looking at the Vertex plot after performing all the other cuts. The background was parameterized with the exponential function, using data to the left and to the right of the target. The uncertainties were estimated by comparing exponential functions with different ranges. Once determined the fit, I subtracted the background from the data. I was happy with the fit, when the "clean spectrum" matched with the flipped spectrum, assuming left/right rotational asymmetry for the clean spectrum. The following plots show the details of the analysis for the three energy settings:

195MeV (0, 5, 10, 19)

06.) 07.) 08.) 09.)

330MeV (0, 3, 5, 10)

10.) 11.) 12.) 13.)

495MeV (0, 2, 4, 7)

14.) 15.) 16.) 17.)

Analysis of the Target Frame contribution

The following three plots show the contribution of the target frame to the measured spectra for the three beam energies.

18.) 19.) 20.)

To get rid of the Tg. frame correction in the analysis we consider a theta0 cut [-3.8, -1], which represents the upper part of the acceptance. Now it is interesting to see, how the results change if different interval is considered. The easiest thing is to compare the intergral of the vertex plots for different theta0 ranges, ie. R(x) = Vertex(-3.8,x)/Vertex(-3.8,1), for different cuts on Vertex-Z. One sees, that when moving along the theta0 the ratio is not flat, as I had expected. It seems, that acceptance is still large enough that CS depends on the out-of-plane angle. To get rid of the effect we need to apply more strict cuts on Vertex-z. More strict was the cut, more consistent are the results.

21.)

Hence, I decided to use in the analysis vertex cut [0,10]mm. To estimate the residual effect of the Target frame, I compared the ratio with the one with the cut [5,10]mm (rescaled to the same number of bins). This is show on the plot below:

22.)

Then I fitted this ratio in the region between [-3.0, -2.5] in order to get correction parameter and the corresponding error. I just need to reconsider why I have chosen this range?!

Latest results:

Following two plots show the present status of the results. I have added the 195MeV data so that we can now see the full range in Q2 that we can cover. Furthermore, I have now properly considered the Snout and Frame contributions for the 495MeV setting. In the previous analyses, this was done by a simple approximation. Now I as still trying to determine the best RadCorrFac factors for the elastic settings. Once I will know them, I will run the simulation with full statistics.

23.)

24.)

Simulation with the Lupolen target:

I have started working on a simulation to see, if the solidstate (Lupolen) target can be used for the ISR experiment. The idea is measure each setting with Lupolen and Carbon targets, and then subtract the two to isolate the Hydrogen spectrum. The question is, how well can we do this. Therefore, I added a new model to the Simul++, which simulates Lupolen target (in the first order Peaking approximation). I compared my results with the data taken during the Calib2014 beam time. Unfortunately are the data available only for spectrometer A. Although the C12 peak dominates, the H spectrum is still clearly visible. Furthermore, when running at low Q2, the excited states will disappear and there will also be no QE peak at low E'.

25.)

In the next step I will improve this first attempt and consider our full Hydrogen simulation in order to get a feeling how much statistics would be needed for such experiment and what is the precision that we can expect.