# G2p Analysis Minutes

Minutes of the weekly analysis meetings

## Contents

- 1 7/1/2015
- 2 6/24/2015
- 3 6/17/2015
- 4 6/10/2015
- 5 5/27/2015
- 6 5/20/2015
- 7 5/13/2015
- 8 5/6/2015
- 9 4/29/2015
- 10 4/22/2015
- 11 4/15/2015
- 12 4/8/2015
- 13 4/1/2015
- 14 3/25/2015
- 15 3/18/2015
- 16 3/11/2015
- 17 3/4/2015
- 18 2/25/2015
- 19 2/18/2015
- 20 2/11/2015
- 21 2/4/2015
- 22 1/28/2015
- 23 1/21/2015
- 24 1/14/2015
- 25 1/7/2015
- 26 July-Dec 2014
- 27 Jan-June 2014
- 28 June-Dec 2013
- 29 Jan-May 2013
- 30 April-Dec 2012
- 31 Jan-March 2012
- 32 July-Dec 2011
- 33 Jan-June 2011

## 7/1/2015

Present: Chao, Vince, Jie, JP, Min, Melissa

By Phone: Pengjia

**Feature Presentations:**

- Min
- Gave an update on acceptance studies. She showed a summary table of each energy

setting and septum configuration. There hasn't been any problems in the 2.2 GeV, 5T

longitudinal and 1.7 GeV, 2.5T transverse settings, but some settings (2.2 GeV, 2.5T

transverse, 1.2 GeV, 2.5T transverse and 2.2 GeV, 5T transverse) have a discrepancy

in the focal plane between data and simulation. She compared the data divided by the

Mott XS to the simulation results without XS, and found that the simulation results

of θ and φ are narrower then the data. JP commented that, while including the Mott

XS might change the shape of the distribution, it shouldn't change the boundary. She

work on figuring out the cause for this discrepancy before calculating the acceptance.

More details can be seen in her slides here.

- Gave an update on acceptance studies. She showed a summary table of each energy

- Pengjia
- Comparing asymmetries determined from the MAID model to those calculated from data.

Last time he showed a factor of 6 difference between the two results. He determined the

scattering angle and Q^{2}by fitting the data and compared them to the quantities calculated

from MAID. He also showed the results of calculating the quantities g_{1}, g_{2}, F_{1}, F_{2}, A_{L}, A_{T},

dXS_{L}, dXS_{T}, XS_{tot}andXS_{mott}using MAID. He hasn't included radiative corrections in his

calculations, but this probably won't account for the factor of 6. More details can be seen

in his slides here.

- Comparing asymmetries determined from the MAID model to those calculated from data.

## 6/24/2015

Present: JP, Min, Chao, Jie, Jixie

By Phone: Vince, Toby, Ellie, Pengjia, Karl

**Feature Presentations:**

- Toby
- Presented some asymmetry plots which Ryan will show in the Chiral Dynamics conference.

The plots can be seen here. All statistics has been included in these plots along with

the dilution and the polarizaiton results. The calculated asymmetries is compared with fully

radiated MAID model. Toby will continue working on the physics asymmetries of all the

kinematics settings.

- Presented some asymmetry plots which Ryan will show in the Chiral Dynamics conference.

- Jie
- Gave an update on his study of the acceptance. He studied the multiple scattering effect

of the incoming beam. The standard deviation of this effect is ~0.07mm for position and

~0.7mrad for angle. He also studied the influence of the beam position to the acceptance.

There was a large dependence on y_tg (horizontal direction) but not on x_tg (vertical direction).

2-D plots shows that the acceptance is not uniform along the vertical direction. JP suggests

that the horizontal target field may be a reason of this effect. More details can be seen in

his slides here.

- Gave an update on his study of the acceptance. He studied the multiple scattering effect

## 6/17/2015

Present: JP, Min, Jie, Chao, Melissa

By Phone: Alex, Ryan, Toby, Ellie, Pengjia, Vince

**Feature Presentations:**

- Ryan
- Working on inelastic polarized radiative corrections, based on the formalism by Akushevich,

Ilyichev and Shumeiko. Since there is no polarized proton data yet, he is using models to

understand the systematic error on the radiative corrections. He used the MAID 2007 model to

generate the polarized cross sections, and went through the same iterative process he used for

the unpolarized radiative corrections. To do the systematic study he took the MAID 2007 cross

sections, radiated them using the POLRAD code, unfolded/radiatively corrected the cross sections,

then compared them to what he started with; the difference is attributed to systematic error. In

general, the ratio of unfolded/radiated cross section is close to 1, but there is a jump when Δσ

crosses zero. JP suggested an alternate method of doing the corrections is the correct on the spin

up XS and spin down XS separately, since neither will be zero. Ryan will look into this method as

a cross check. More details can be seen in his slides here. - Working on his talk for Chiral Dynamics, he will send out his slides to the collaboration later

this week.

- Working on inelastic polarized radiative corrections, based on the formalism by Akushevich,

- Toby
- Working on fixing the scattering angle dependence on the cross section. His original method

was to plot the calculated scattering angle, take the central value and plot it versus HRS

momentum. Fitting this with an exponential function (originally provided by Jixie), he could

then use the fit parameters in the Bosted model to calculate the cross section. The problem

with this method is it forces the scattering angle across kinematics to be continuous, but our

data is not continuous. He generated the model at all scattering angles within a momentum

setting and combined all resulting models with a weighted average. This causes the model to be

discontinuous over an energy setting, but this method of generating the model may be necessary

since the continuous model doesn't represent our data well. JP cautioned that the acceptance

could have an effect on the scattering angle, so we must be cautious using this data until

optics/acceptance studies are finalized. More details can be seen in his slides here.

- Working on fixing the scattering angle dependence on the cross section. His original method

- Min
- Gave an update on her acceptance study. She divided the data by the Mott cross section before

comparing it to the simulation results for an empty dilution run (just helium). For the scattering

angle she is using the survey result, which has an uncertainty of 0.7mrad. JP commented that we must

be cautious in how we determine the scattering angle, particularly the effect of the septum magnet

on the angle. For next time, Min will work on tuning the apertures to have better agreement between

data/simulation. More details can be seen in her slides here.

- Gave an update on her acceptance study. She divided the data by the Mott cross section before

- Chao
- Gave a status update on optics, specifically the offset problem he has discussed previously. He

treats the matrix elements separately and uses the center hole to determine the offset and 2 first

order matrix elements. He then uses the beam position scan data to determine matrix elements that

are not related to x_{fp}and delta scan data to determine the matrix elements related to x_{fp}.

However, using the updated matrix (for the 1.7 GeV data) still has a 3mm offset. JP questioned

whether the effect of the septum is being taken into account properly. Once we understand how to

treat the septum well, we'll know whether to attribute problems to the septum field, target field,

etc. Chao will go back and look into the effect of the septum before moving forward on the offset

problem. More details can be seen in his slides here.

- Gave a status update on optics, specifically the offset problem he has discussed previously. He

**General Discussion:**

- The abstract deadline for DNP is July 1st.

## 6/10/2015

Present: JP, Min, Chao, Jie, Kalyan, Melissa

By Phone: Vince, Ryan, Toby, Ellie, Pengjia, Karl

**Feature Presentations:**

- Jie
- Gave an update on his study of the change in acceptance/yields due to shifts in the

beam position using simulation. Last time he showed the number of events blocked by

different components such as the collimators, and magnet entrance/exit planes. He looked

at the change in acceptance/yields moving along BPM Y and X and found that there was

not a large dependence on Y, but the dependence on X was significant. In addition,

removing the Q1 exit plane aperture while moving along BPM X showed a large effect on

acceptance/yield ratio on the positive x side. He also looked at the effect of extending

X to +30mm. He found that moving the target ~10mm gives the largest acceptance at the

Q1 exit endplane. There was some confusion about coordinate definitions, Jie will clear

this up for next time. More details can be seen in his slides here.

- Gave an update on his study of the change in acceptance/yields due to shifts in the

## 5/27/2015

Present: JP, Min, Jie, Chao, Jixie, Melissa

By Phone: Pengjia, Ryan, Toby, Ellie, Alexandre

**Feature Presentations:**

- Pengjia
- Comparing calculated asymmetries from data to asymmetries calculated from the MAID model.

He looked at the 2.2 GeV, 5T settings (both longitudinal and transverse) and used a dilution

factor of 0.15. He needed to scale the MAID model as it was 6X larger than the asymmetries

calculated from data. He also included a plot of Q² vs W for each setting. JP commented

that the distribution was a little strange, since Q² and W are directly correlated the

distribution shouldn't decrease or stay flat. More details can be seen in his slides here.

- Comparing calculated asymmetries from data to asymmetries calculated from the MAID model.

- Min
- Gave an update on acceptance studies, specifically the 3rd septum setting. In this septum

configuration, the entire second coil is gone. She looked at groups of runs where the beam

position was stable to determine the uncertainty of the focal plane measurement. She compared

the y and φ at the center sieve hole at the focal plane determined from data to the SNAKE

result. The center sieve hole was defined by survey results. JP suggested using the no-target

field setting as a starting point to determine the offset, then look at the transverse effects.

More details can be seen in her slides here.

- Gave an update on acceptance studies, specifically the 3rd septum setting. In this septum

**General Discussion:**

- There will be no meeting next week due to the Users Group Meeting.

## 5/20/2015

Present: JP, Min, Jie, Jixie, Melissa

By Phone: Pengjia, Ryan, Toby, Kalyan

**Feature Presentations:**

- Pengjia
- Calculated asymmetries using the MAID model. He used the total cross section output from

MAID 2007, then summed the cross sections from two channels to extract the asymmetry. He

calculated the longitudinal/transverse asymmetries and g1,g2,F1 and F2 for several different

values of Q². Ryan pointed out that the calculated asymmetries cannot be compared directly

to our data, since our data is not at constant Q², and the MAID model is not radiated. More

details can be seen in his slides here.

- Calculated asymmetries using the MAID model. He used the total cross section output from

- Jie
- Investigating the field overlap of the septum and Q1 fields in the SNAKE model. Currently,

no overlap is assumed between the two fields. He looked at different points in the space

between the septum and Q1 entrance to check the field from each magnet. Although the Q1

field goes to zero close to the septum, the septum field extends to the entrance of Q1. The

current septum field map is too small to cover the full area in the x,z direction, so a larger

field map is needed. JP also suggested looking at how this will change the particle

trajectories; Min will work on estimating this for next time. More details can be seen in

Jie's slides here.

- Investigating the field overlap of the septum and Q1 fields in the SNAKE model. Currently,

## 5/13/2015

Present: JP, Chao, Min, Jie, Melissa

By Phone: Pengjia, Ryan, Toby, Karl

**General Discussion:**

- Jie posted a draft of a technote on radiation effects in the g2psim package here.

- Toby posted a draft of a technote on dilution analysis here.

**Feature Presentations:**

- Melissa
- Looking into packing fraction runs for the 1.1 GeV setting, taken with the short ammonia cell.

She doubled checked the normalization constants used to calculate the yields and didn't find

anything strange. She also looked at the s1 and Cherenkov channels for each of the four runs,

and found some small differences, but this is probably due to the fact that the runs were taken

at significantly different beam currents. Two of the runs (5197,5198) were taken with a beam

current of 20nA (compared to 40-50nA that the other runs were taken at), which means the BPMs

are not reliable for these runs, so the beam position may be significantly different for these

runs. It's possible there is some effect due to the high rate with a beam current of 40-50 nA

that is not being accounted for. More details can be seen in her slides here.

- Looking into packing fraction runs for the 1.1 GeV setting, taken with the short ammonia cell.

## 5/6/2015

Present: Chao, Min, Jie, Melissa

By Phone: Pengjia, Ryan, Ellie

**Feature Presentations:**

- Chao
- Gave a status update on optics analysis. He found that using the 2.2 GeV 5T longitudinal matrix

to replay the 1.7 GeV 2.5T transverse data resulted in an overall offset. Because of this, it is

useful to do a full check of the relation between the reconstructed kinematics and the beam position.

If there is a correlation between the beam position and reconstructed kinematics, a linear fit will

be used to determine the correction. The fitting result for φ for the two different settings is added

directly to the first order matrix element P0000 (constant term). These matrices will be used to

replay the two data sets, and then the fitting procedure can be applied to improve higher order terms.

Next he will work on checking other energy settings. More details can be seen in his slides here.

- Gave a status update on optics analysis. He found that using the 2.2 GeV 5T longitudinal matrix

- Min
- Showed an updated on focal plane matching. For the 2.2 GeV, 0T setting with the "good" septum

(484816), the data and simulation match well at the focal plane. However, data from the same setting

but with the target field at 2.5T do not match as well. She tried adjusting the beam position to find

a better match, but the correction needed was larger than the uncertainty on the beam position, so

JP suggested looking at the discrepancy starting from the central hole. By adjusting the beam position

and dipole field, she was able to match data from the central sieve hole to simulation results, within

uncertainties. Next she will work on fitting the other sieve holes to match the data. More details

can be seen in her slides here.

- Showed an updated on focal plane matching. For the 2.2 GeV, 0T setting with the "good" septum

## 4/29/2015

Present: JP, Chao, Min, Jie, Melissa

By Phone: Alex, Pengjia, Toby, Ryan, Ellie, Karl

**Feature Presentations:**

- Pengjia
- Looked at the effect of different raster cuts on the asymmetry. He looked at the 2.2 GeV 5T

settings, both longitudinal and transverse configurations. First, he placed a cut on the center

of the raster (50% of total radius), which didn't have a large effect on the asymmetry. He also

looked at the effect of splitting the beam spot into 4 quadrants; for some asymmetries there was

a considerable difference between the quadrants. He also found that the yield is larger on the

left side of the beam spot. More details can be seen in his slides here.

- Looked at the effect of different raster cuts on the asymmetry. He looked at the 2.2 GeV 5T

- Melissa
- Showed the results of packing fraction analysis for the 1.1 GeV setting. There is some

discrepancy in the yields for some materials, and the ratio of the rate/current is not consistent

in some cases. She showed the effect of applying acceptance and PID cuts to each run; the

makes the largest impact is the cut on single track events, which cuts out ~50% of the total

events in some cases. She will continue trying to figure out the cause of the variance in the

yields. More details can be seen in her slides here.

- Showed the results of packing fraction analysis for the 1.1 GeV setting. There is some

## 4/22/2015

Present: JP, Chao, Min, Jie, Melissa

By Phone: Toby, Ryan, Ellie, Karl, Vince

**Feature Presentations:**

- Chao
- Gave an update on optics analysis. Using data from the single foil target does not provide enough

constraints for optics analysis, but there are two other methods to additionally constrain the data.

One way is to look at the relation between the kinematics and beam position for beam position scan

runs, and the other way is to include data from the aluminum window into the fitting procedure to

directly correct the reconstruction matrix. The second method is what was tested for this week. To

improve the fitting routine, he first fits the k=0 matrix elements without the Aluminum window data,

then fixes the other matrix elements and fits the k=l matrix elements after adding in the aluminum

data. He found that he couldn't get a good fit using the aluminum window data, even if only the y-fp

matrix elements are allowed to vary. JP commented that, since y and φ are coupled, they cannot be

fit separately. It's possible that a higher order polynomial is necessary for the fit. More details

can be seen in his slides here.

- Gave an update on optics analysis. Using data from the single foil target does not provide enough

- Toby
- Put together a summary of the data quality checks on the yields (detailed here), including the

number of useable runs and percentage of useable statistics in each kinematic setting. Overall, the 5T

settings look the best, with the exception of the 3.3 GeV setting, which has several settings with yield

drifts that appear to be unrelated to beam position changes. Karl suggested checking the stability of the

dipole for settings with bad runs. More details can be seen in his slides here.

- Put together a summary of the data quality checks on the yields (detailed here), including the

**General Discussion:**

- Pengjia posted a draft of a technote on beam charge measurements here, and would appreciate feedback.

## 4/15/2015

Present: Jixie, Min, Chao, JP, Jie, Melissa

By Phone: Vince, Pengjia, Toby, Ryan, Ellie, Karl

**Feature Presentations:**

- Jie
- Gave a simulation update, specifically looking at acceptance and yields. Last time, he used a point

beam a moved the position from 0-10 in x and y to see the effect on the acceptance and the yields. This

time, he looked at the acceptance at each endplane using a point beam at (0,0) and (-10,0). One big

difference is the number of events blocked by the Q1 exit for these two beam positions. He went step by

step through the simulation to compare the number of events that reach each enplane to the number of

events that pass through the plane. He found that the shift of the beam position by 1cm caused a drop of

~20% in the acceptance. JP and Vince commented that this drop seems too large. It's possible the septum

is also playing a role in the decreased acceptance, not just Q1. More details can be seen in his slides here

- Gave a simulation update, specifically looking at acceptance and yields. Last time, he used a point

- Melissa
- Gave an update on packing fraction analysis. She extracted pf values for numerous elastic runs for

several settings. For settings where there is a drift in the yields, the pf also varies. The most

variation is seen in the 2.2 GeV, 2.5T setting. The raster size was changed part way through this

setting, so it may not be possible to compare runs with different rater sizes without first correcting

for this difference. She will continue working on trying to understand the discrepancies in the yields,

and finish extracting the values for the 1.1 GeV setting, for the regular and short ammonia cells. More

details can be seen in her slides here.

- Gave an update on packing fraction analysis. She extracted pf values for numerous elastic runs for

- Ryan
- Gave an update on radiative corrections. He took a closer look at the systematic error on the

inelastic radiative corrections, which are different now that he has removed the interpolation from the

unfolding procedure. He replaced the interpolation with a direct call to the Bosted Model. Comparing

the two methods, he found that they are consistent with each other. The interpolation will be the

largest contribution to the systematic error. He also wanted to test if there was a faster way to

radiate a model. There are two ways to do the integration within the Stein equation, Romberg or Simpson

integration. The Romberg method is more accurate, but takes longer. The results from using both methods

to radiation the Bosted model are nearly identical, but the Simpson method is orders of magnitude faster,

particularly for finer ν binning. Next he will work on setting up polarized radiative corrections.

More details can be seen in his slides here.

- Gave an update on radiative corrections. He took a closer look at the systematic error on the

- Toby
- Gave an updated on dilution analysis. He showed a central fit to the scattering angle compared to the

model prediction done by Jixie. This is used to scale the angle between different run types, such as

scaling carbon dilution to nitrogen production. Since the optics isn't complete yet, this scattering angle

is just a place holder until we better understand the scattering angle. At large ν, the dilution factor

(and packing fraction) starts to drop off, which he thought could be due to a problem with the scattering

angle correction. JP commented that this might not be the case, since this is a ratio, but due to the fact

that at low Q squared the DIS cross section is small. This will probably result in a large statistical

uncertainty on the dilution factor. Next Toby will extract the dilution factor for the 1.1 GeV and 1.7 GeV

settings, and finish writing a technote on this analysis. More details can be seen in his slides here.

- Gave an updated on dilution analysis. He showed a central fit to the scattering angle compared to the

## 4/8/2015

Present: Chao, Min Kalyan, Jie, Melissa

By Phone: Alex, Vince, Toby, Ryan, Karl, Ellie, Pengjia

**Feature Presentations:**

- Chao
- Gave an update on optics analysis. Previously he showed that there is an issue with the horizontal beam

position changing the dp reconstruction. Since there aren't enough constraints to do the reconstruction, it

is necessary to do a full check of the relations between the beam position and the reconstructed kinematics.

He looked at delta scan runs and beam position scans to check the correlation between the beam position

and reconstructed variables. To correct for these effects, he used a linear fit to determine a correction factor

for θ, φ and dp. To test this method he used a production run, and placed cuts on the raster pattern in

different areas around the center of the beam spot. While the "old" reconstructed dp variable shifted based

on the location of the cut, the "corrected" dp variable was stable. Next he will work on applying this method

to other settings. More details can be seen in his slides here.

- Gave an update on optics analysis. Previously he showed that there is an issue with the horizontal beam

## 4/1/2015

Present: Min, Chao, Kalyan, JP, Jie, Jixie, Melissa

By Phone: Ryan, Toby, Pengjia, Vince, Alex

**Feature Presentations:**

- Melissa
- Gave an update on packing fraction analysis. She started extracting the pf for the 2.5T settings, but the biggest

lingering issue is the discrepancy seen in the yields for elastic pf/production runs. Several suggestions were made

to check for a possible problem, such as the septum field drifting over time, looking at the raster pattern, and

checking the scalar rates and normalization values (specifically the charge). More details can be seen in her slides here.

- Gave an update on packing fraction analysis. She started extracting the pf for the 2.5T settings, but the biggest

- Ryan
- Gave an update on automated model tuning for determining the carbon/nitrogen cross section ratio. Previously, the

tune was optimized in the order dip-region, inelastic, quasielastic. For this time he swapped the order to inelastic,

quasielastic, dip-region. The final overall tune represents an average of each settings parameters. He found that

the best average reduced chi-squared was found after the second iteration, and the fit is good to within 10%. Next

he will look at saGDH nitrogen data, but is waiting on updated cross sections from Vince. More details can be seen

in his slides here.

- Gave an update on automated model tuning for determining the carbon/nitrogen cross section ratio. Previously, the

- Min
- Gave an update on acceptance studies. Last time she compared the data from a delta scan in the 2.2 GeV, 2.5T

setting (with good septum) with simulation results and found some differences. This time she adjusted the beam

position in the simulation to match the data, which was quite different from the original beam position. Kalyan

suggested looking at just the center hole of the sieve slit to determine the offset. It seems the offset is seen

in the data, but the simulation shifts the events in the opposite direction. Min/Chao have a plan to check each

part of the simulation, and Pengjia will check the beam position calibration for this setting. More details can

be seen in her slides here.

- Gave an update on acceptance studies. Last time she compared the data from a delta scan in the 2.2 GeV, 2.5T

- Jie
- Gave a simulation update, specifically checking the target center compared to the position given by the BPM at

the target. His goal was to test if the target was centered, whether the BPM could have an offset? Or if the

target was centered and the BPM position was correct, could the raster size be scaled incorrectly? Looking at

several runs he saw a shift of 1cm in the target center. However, as the runs were at different energy settings,

the beam position wouldn't necessarily be the same. Kalyan suggested looking at the carbon hole runs nearest to

these runs as a check. For next time he will focus on simulating the shifts we see in the yields. More details

can be seen in his slides here.

- Gave a simulation update, specifically checking the target center compared to the position given by the BPM at

## 3/25/2015

Present: Jie, Min, Vince, JP, Chao

By Phone: Ryan, Karl, Toby, Alexandre, Pengjia

**Feature Presentations:**

- Jie
- Gave an update on simulation package. The elastic radiation tail include soft photon contribution is compared with

the result from the predication calculated using Mo/Tsai. The tail matches at ΔE~6MeV. He also looked at the effect

of the beam position on the acceptance. The simulation covers the full angle acceptance to test the aperture cut

in the SNAKE model. The simulation shows that φ acceptance has large correlation with the horizontal beam position.

JP suggests the aperture setting in the model may be not sufficient to represent the real situation. Chao suggests

that to use the simulation to determine which aperture stops most particles. More details can be seen in his slides here.

- Gave an update on simulation package. The elastic radiation tail include soft photon contribution is compared with

- Chao
- Gave an update on the optics calibration. As Jie points out before, the reconstructed dp does not agree with the

prediction very well when horizontal beam position changes. The beam position scan runs are used to fit a correction

to the dp reconstruction matrix. The correction could correct the dp reconstruction result for optics runs. However, the

reason of this correction is still not clear. After discussion, the suggestion is to add the correction back to the dp matrix

and redo the fit of higher order terms. More details can be seen in his slides here.

- Gave an update on the optics calibration. As Jie points out before, the reconstructed dp does not agree with the

- Toby
- Gave an update on dilution study. He showed dilution results for all completed settings. The tech note of the dilution

study is in progress. He also mentioned an alternative way to calculate the packing fraction. The method only works in

DIS region but a 0th order polynomial fit at large ν can be used to find the packing fraction everywhere. All packing

fractions seem about 5-10% larger than typical values. More details can be seen in his slides here.

- Gave an update on dilution study. He showed dilution results for all completed settings. The tech note of the dilution

- Pengjia
- Presented some preliminary results of asymmetry and yield for RHRS. He compared the right arm asymmetry and yield

with left arm for two kinematic settings. The asymmetry results agree quite well but the yields has some deviation for

the 2.2GeV, 5T longitudinal field setting, perhaps caused by acceptance effect. More details can be seen in his slides here.

- Presented some preliminary results of asymmetry and yield for RHRS. He compared the right arm asymmetry and yield

## 3/18/2015

Present: Jixie, Chao, Jie, Min, Kalyan, JP, Melissa

By Phone: Ryan, Karl, Toby, Alexandre, Vince, Pengjia

**Feature Presentations:**

- Melissa
- Gave an update on packing fraction analysis. For each setting, she tried to identify ammonia runs that could

be used for the analysis, ideally runs that were taken in close proximity to the dilution (dummy/empty) runs, or

at least where the beam position is similar. In addition, the fitting routine (to understand the level of

contamination from the 2nd peak) was adjusted for each energy setting. The 2.2 GeV, 5T settings (longitudinal

and transverse) have ammonia runs taken close in time to the runs to helium. The 1.7 GeV setting has some

variation in the ammonia yields, but they are not too large. The 1.2 GeV energy setting has significant variation

in the elastic yields; she will look more into these runs to see if there is an obvious reason for the variation.

She also showed the fit results for the different energy settings; JP suggested checking the hydrogen fit, as it

seems a bit wide. More details can be seen in her slides here.

- Gave an update on packing fraction analysis. For each setting, she tried to identify ammonia runs that could

- Ryan
- Working on an automated model tune, which will be used to determine the ratio of carbon to nitrogen for dilution

analysis. Last time, he showed an automated tune of the Bosted model, which assumed that the quasi-elastic and

inelastic channels are independent one another and can be optimized independently. The caveat to this is that the

tune is only valid in the Q^{2}range and for the target type of the input data. For this time, he updated the

automated tune procedure to include an adjustment for the "dip" region. The Bosted model parametrization of the

dip region is a Gaussian that is calculated as part of the inelastic channel. Including the dip region improves the fit

to carbon data for each individual setting, however the overall fit it worse. JP suggested that fitting the dip region

first may be the problem; Ryan will try to swapping the order of the fitting routine to see if it improves the fit.

More details can be seen in his slides here.

- Working on an automated model tune, which will be used to determine the ratio of carbon to nitrogen for dilution

- Min
- Gave an update on her acceptance study. To improve the simulation, she included data from the windows into the

calibration and updated the results for the delta scan runs. This new function was applied to both horizontal and

vertical beam scan runs. This improved the agreement between the data and simulation, especially for the horizontal

beam scan run. She also looked at beam scan runs with the target field on, this time adding an additional constant

on y and φ. For the fit of the focal plane, JP cautioned not to use too high-order of a fit, as it can become

non-physical. For next time, she will adjust the fitting function. More details can be seen in her slides here.

- Gave an update on her acceptance study. To improve the simulation, she included data from the windows into the

- Pengjia
- Looked at the effect of raster cuts on the yield. He looked at several kinematic settings with consistent yields

and a applied a raster cut with a diameter of 6mm. He found that applying the raster cut had a negative effect;

meaning the spread in the yields became larger after the cut was applied. He also showed the x and y positions at

BPM A and B for each set of runs. In some cases, the x-position at BPMA is stable, but it varies at BPMB, which

could be the source of the variation in the calculated scattering angle. It may be more useful to do a relative

study, (similar to what Jie has shown previously) for individual runs than trying to compare separate runs. More

details can be seen in his slides here.

- Looked at the effect of raster cuts on the yield. He looked at several kinematic settings with consistent yields

## 3/11/2015

Present: Chao, Jie, Min, JP, Melissa

By Phone: Vince, Toby, Ryan, Karl, Alexandre

**Feature Presentations:**

- Jie
- Working on updating the energy loss model in the simulation package. Previously it was found that the energy loss

model in g2psim does not match the prediction calculated using the Mo/Tsai formulation. He found that the soft photon

contribution was missing. This term could have a large contribution, up to 20%. He also looked at the effect of

changing the scattering angle on the dp distribution. Over an 80mrad change in the scattering angle, he saw only a

~0.16% change in the dp distribution. He also looked at 5 different areas of the beam spot, to see how the dp changed

in that region. The simulation results show no change in the dp distribution in the different regions, but the data

shows that the dp distribution shifts as you move around the beam spot. More details can be seen in his slides here.

- Working on updating the energy loss model in the simulation package. Previously it was found that the energy loss

**General Discussion:**

- Ryan sent around an abstract for the upcoming Chiral Dynamics workshop, comments are appreciated.
- JP sent around the latest version of Pengjia's BPM paper, feedback is appreciated.

## 3/4/2015

Present: Chao, Jie, Min, JP, Melissa

By Phone: Vince, Toby, Ryan, Karl, Pengjia

**Feature Presentations:**

- Toby
- Gave an update on dilution analysis. He finished tuning the Bosted model in order to obtain cross section ratios used

for radiation length matching for dilution runs (for example, scaling a helium run to match the radiation length of a carbon

run). He showed the calculated background over the full kinematic range for the 2.2 GeV, 5T, Transverse setting. The

dilution factor can then be extracted using the ratio of the background yield to the total yield. Eventually he will be able

to update the dilution factor using material specific packing fraction values (when they are ready), and an improved tune

of the XS model. Next he will work on other settings and compiling a technote. More details can be seen in his slides here.

- Gave an update on dilution analysis. He finished tuning the Bosted model in order to obtain cross section ratios used

- Ryan
- Working on a method to automate model tuning for the QFS and Bosted models. For dilution analysis it is necessary

to scale carbon data as a replacement for nitrogen yields. Since there is no nitrogen and carbon cross sections available

from the same experiment, it is necessary to use to model tune of saGDH data for nitrogen, and a model tune of data

from UVa archives for carbon. Previously, this tuning was done "by hand", with a reduced χ^{2}to determine the

quality of the fit. A problem with this method is that it could result in inconsistent results for different people tuning the

same data set; an automated method could solve this problem. JP commented that a possible problem with an automated

method is that the physics is simplified. He showed results for tuning the Bosted model to carbon data, next he will work

on saGDH data, and then will repeat the procedure for the QFS model. More details, including a description of the code

structure, can be seen in his slides here.

- Working on a method to automate model tuning for the QFS and Bosted models. For dilution analysis it is necessary

- Melissa
- Gave an update on packing fraction analysis. Previously, she described an update to her method to include cross section

input. Elastic form factor models were used to determine the cross section ratio σN/σHe. To obtain radiated cross sections

two methods were used. First, using the g2psim energy loss model, and second, using the Mo/Tsai formalism to calculate

the correction to the elastic peak. There is a discrepancy between these methods, specifically that the results suggest that

the correction is smaller than it should be. This will have to be resolved to get an accurate value for the cross section ratio.

She is also working on fitting the data for dilution and production runs for other energy settings to extract packing fraction

values for other materials. More details can be seen in her slides here.

- Gave an update on packing fraction analysis. Previously, she described an update to her method to include cross section

- Min
- Gave an update on her acceptance study. Last time she found that the focal plane variables for production runs didn't

match well with the simulation. Production runs have a rastered beam, as opposed to the point beam used in optics

calibration runs, so this effect must be corrected for. For this week, she looked at events at the target window in the focal

plane, and used these events to improve the fitting at the virtual plane. She also looked at events from the carbon foil at

the focal plane, and saw a reasonable match. Next she will work on further improving the fitting and incorporating the

results into g2psim. More details can be seen in her slides here.

- Gave an update on her acceptance study. Last time she found that the focal plane variables for production runs didn't

## 2/25/2015

Present: Chao, Jie, Min, JP, Melissa

By Phone: Toby, Ryan, Ellie, Pengjia

**Feature Presentations:**

- Jie
- Working on understanding the yield drifts due to changes in the beam position. He used beam position information

as input to the simulation, and applied loose acceptance cuts to determine the yield. He looked at several momentum

settings that had a drift in yields >3% in the 1.7 and 1.1 GeV energy settings. In general he found that the effect was

opposite of what is seen in the data. JP suggested doing the study independent of data, so as not to bias the results,

and to show the dependence of the yield on each parameter (x,y,θ,φ) individually. He also suggested checking to

make sure the acceptance effects are included properly, as it is difficult to get them to agree at the boundaries. Jie will

keep working on this, and will start working on the septum SNAKE model for the 403216 septum configuration. More

details can be seen in his slides here.

- Working on understanding the yield drifts due to changes in the beam position. He used beam position information

- Pengjia
- Trying to determine if applying a raster cut to the data will fix the spread in yields. The difficulty is in determining

the total charge after this cut is applied. To determine the charge, he used the calibrated charge info from both the

happex and fastbus, and applied the same raster cut to each. He was then able to use the happex raster cut info to

get the charge from the happex bcm. JP pointed out that, since the happex and fastbus aren't synchronized, this

method may only work if the beam is stable, so it is important to quantify an uncertainty on the charge. He suggested

trying this method again with a much smaller raster cut, so that whatever effect this method will be magnified. Pengjia

will first look at a setting where the yields are stable to see the effect of this cut, then will move on to settings that have

issues. The raster cut library and testcode is available on the work disk, see his slides here for more details.

- Trying to determine if applying a raster cut to the data will fix the spread in yields. The difficulty is in determining

- Chao
- Posted the draft of a technote on target field mapping here. Any comments would be appreciated!

## 2/18/2015

Present: Chao, Jie, Kalyan, Melissa

By Phone: Vince, Ryan, Toby, Min

**Feature Presentations:**

- Min
- Gave an update on her acceptance study. She showed a comparison of simulation results with data from a production

run. As was suggested previously, she expanded the acceptance range when generating events. She found that the

focal plane variables still don't match very well. One possibility is that the production runs have a raster beam, but the

optics runs she had looked at previously used a point beam, so this effect needs to be corrected for somehow. More

details can be seen in her slides here.

- Gave an update on her acceptance study. She showed a comparison of simulation results with data from a production

**General Discussion:**

- Chao mentioned that a problem was found in the most recent BPM package, and it may have caused a problem

in the replay. Toby will discuss this with Pengjia offline.

## 2/11/2015

Present: Chao, Jie, Min, Melissa

By Phone: Pengjia, Toby, Ryan, Ellie, Vince

**Feature Presentations:**

- Toby
- Gave an update on dilution analysis. The dilution factor is define as one minus the ratio of background yield over

total yield. The background yield includes contributions from nitrogen, helium and the foil endcap on the target cell.

To extract the dilution factor he must tune the Bosted model (for helium and carbon) to match cross section data,

radiate the simulation results and take ratios to obtain necessary scaling constants, and then apply these scaling

constants to calculate the background yield. For tuning the Bosted model, he is currently using cross section data

from the UVa quasi-elastic database. The kinematics do not match ours exactly, but this will work until g2p cross

sections are available. He used a χ-squared minimization routine to determine the scaling factor for each kinematic.

He also showed the resulting cross sections ratios, which are used to determine the scaling factors α, β and γ.

Next he will apply the cross section ratios to yields from dilution runs to calculate the background yields and dilution

factor for the 2.2 GeV, 5T transverse setting, after which he will move on to other settings. More details can be seen

in his slides here.

- Gave an update on dilution analysis. The dilution factor is define as one minus the ratio of background yield over

## 2/4/2015

Present: Chao, JP, Jie, Min, Melissa

By Phone: Toby, Ryan, Pengjia, Karl

**Feature Presentations:**

- Chao
- Gave a status update on optics analysis. For the 2.2 GeV, 5T settings (longitudinal and transverse) there was

limited optics data taken. Instead, the 1.7 GeV optics matrix matrix was used for this setting, with a 0th order

correction for the longitudinal setting, and a 1st order correction for the transverse setting. He showed the

longitudinal setting as an example, applying the 1.7 GeV matrix and including the 0th order correction to center

the sieve holes. JP commented that, if the offset is large, there may be higher order effects to consider. There

are no optics runs for the 3.3 GeV setting, so he plans to use g2psim to determine the optics matrix. He is

working on an optics technote, and a technote summarizing the last target field measurement. More details can

be seen in his slides here.

- Gave a status update on optics analysis. For the 2.2 GeV, 5T settings (longitudinal and transverse) there was

- Min
- Gave an update on her acceptance study. For this time, she improved the fitting at the virtual plane (for both

x/y and , and θ/φ) and projected the results onto the focal plane. The corrections on the forward transport

function were incorporated into g2psim. For the reverse transport function, the fit was done using all dp scan runs.

The simulation does not match with the data at the edges; this is due to the boundary conditions in the simulation.

Min will try expanding the boundary to see if it better matches the data. More details can be seen in her slides here.

- Gave an update on her acceptance study. For this time, she improved the fitting at the virtual plane (for both

- Melissa
- Updated the data quality check for the 2.2 GeV, 2.5T setting to include RHRS runs. In general the problem

settings on the RHRS are consistent with the LHRS. She also gave an update on the packing fraction analysis;

there was a problem in the previous definition of the packing fraction. The updated method requires some input

from cross section models. The resulting packing fraction is slightly larger than the result from the previous

method. Next she will work on extracting packing fraction values for other settings. More details can be seen

in her slides here.

- Updated the data quality check for the 2.2 GeV, 2.5T setting to include RHRS runs. In general the problem

- Jie
- Calculated the energy loss using two methods; using his energy loss model (Monte Carlo, step by step), and

using a radiated cross section model (from Stein). He used a fixed initial energy and fixed scattering angle for

a carbon target, and didn't apply any acceptance cuts. He compared the dp spectra for the two methods and

found that his energy loss model agrees well with the radiation model. Next he will work on a full simulation

for a beam dependence study (once the acceptance is ready). More details can be seen here.

- Calculated the energy loss using two methods; using his energy loss model (Monte Carlo, step by step), and

**General Discussion:**

- Starting next week we will change the meeting start time to 9:30am.

## 1/28/2015

Present: Chao, JP, Jie, Melissa

By Phone: Vince, Ryan, Toby, Min, Pengjia, Ellie, Alexandre

**Feature Presentations:**

- Ryan
- Gave an update on data quality checks for the 1.1 GeV setting. He found that the spread in yields ranges

from 0.01% - 6.1%, with an average of ~2%. For this setting we also took runs on a short target cell (material

14). For settings with a large spread, there was a correlation between beam drift and yield drift. Overall, this

setting seems ok. Next he will move back to working on polarized radiative corrections. More details can be

in his slides here, and a detailed summary

of the data quality check can be seen here.

- Gave an update on data quality checks for the 1.1 GeV setting. He found that the spread in yields ranges

- Pengjia
- Gave an update on data quality checks for the 1.7 GeV setting. For settings with a large spread in yields,

he looked at the correlation with the calculated scattering angle. He found that in some cases the change in

yield was proportional to the change in scattering angle, but for some settings it was inverse proportional.

He also showed an example of a setting where the yield was drifting, but there was no change in the scattering

angle. He tried applying raster cuts to see the effect on the yield, but first needs to find a way to calculate the

charge while including the raster cut. There was some discussion how this could be done; it may be possible to

use the Happex DAQ, as it is triggered by helicity. More details can be seen in his slides here.

- Gave an update on data quality checks for the 1.7 GeV setting. For settings with a large spread in yields,

**General Discussion:**

- Toby will start replaying RHRS runs today, hopefully they will be done by the end of the week.

## 1/21/2015

Present: Jie, Min, Chao, JP, Kalyan, Melissa

By Phone: Pengjia, Vince, Ryan, Toby, Karl

**Feature Presentations:**

- Melissa
- Gave a summary of data quality checks for the 2.2 GeV, 2.5T, transverse setting. There are several momentum

settings with a significant spread (>2%) in the yields. In most cases, the spread in yields seems to be correlated

with a change in beam position, which often corresponds to beam down time during the run period. Examples of

problem settings can be seen in her slides here, and a summary of all settings is available here.

- Gave a summary of data quality checks for the 2.2 GeV, 2.5T, transverse setting. There are several momentum

- Chao
- Gave an update on the status of optics for the RHRS. For the longitudinal and transverse settings at 5T, only a

few optics runs were taken, so the plan is to use the 1.7 GeV optics matrix for these settings. He also included a

summary of the RMS value of each kinematic variable (δ, θ and φ) for each setting, and commented that

the RMS value for y_tg is around 3-4mm, but this is a rough estimation as the fitting is not as good for this variable.

Next he will test 2.5T transverse optics matrix on the 5T settings to see if it works. His slides are available here. - The simulation package has been updated to cover the RHRS and include some upstream geometries. The most

recent version of g2psim is available through github here: https://github.com/asymmetry/g2psim

- Gave an update on the status of optics for the RHRS. For the longitudinal and transverse settings at 5T, only a

- Min
- Gave an update on her acceptance study. For this time, she combined all dp scan runs in order to do the

corrections at the virtual plane, which was then projected onto the focal plane. The agreement at the virtual

plane looks good, but there is still some disagreement at the focal plane. Next she will fit the reverse transport

functions and compare target plane variables. More details are available in her slides here.

- Gave an update on her acceptance study. For this time, she combined all dp scan runs in order to do the

## 1/14/2015

Present: Chao, Jie, JP, Melissa

By Phone: Pengjia, Ryan, Toby, Ellie, Karl, Min, Vince, Alexandre

**Feature Presentations:**

- Jie
- Gave an update on the dp simulation, specifically trying to understand why the simulation result for elastic dp was

narrower than the data. Comparing the data (with just PID cuts) to the simulation result (using full acceptance, and

a raster diameter of 2cm), the simulation is narrower than the data. However, if he applies acceptance cuts to the data

and simulation, and uses a beam size of 4mm (which matches with the data), the results show a better match between

data and simulation. He also tried dividing the beam spot into 5 different zones, and looked at the resulting dp

distribution for each. He did see a shift in the central value of dp between zones, however this shift was not seen in the

simulation. He also compared the shift in θ and φ for the 5 zones. JP commented that the θ distribution in

the data is symmetric and centered around zero, while the simulation result is not. It is possible that this is due to the

fact that Jie is using the old SNAKE model; he will get the updated version from Min and see if it makes a difference.

More details can be see in his slides here.

- Gave an update on the dp simulation, specifically trying to understand why the simulation result for elastic dp was

- Pengjia
- Posted his false asymmetry results to the ELOG.
- Working on data quality check for the 1.7 GeV, 2.5T setting. It looks like the spread in yields in some settings is

strongly dependent on the scattering angle. JP suggested looking at beam position or other raw variables, since the

scattering angle is a calculated value and depends on optics. It seems most settings with problems correspond to a

break in time when the runs were taken, which may suggest a problem with the spectrometer. Pengjia will look into

this next. More details can be seen in his slides here.

- Toby
- Gave an update on his yield drift study for the 3.3 GeV setting. For this time, he looked at the target field (central

NMR value) over the 2.34 GeV momentum setting, a period of about 6 hours. While there are some jumps in the signal,

there isn't a consistent drift like is seen in the yields. However, there is a period where the signal jumps significantly; he

will look into a possible cause for this. He also looked at the scaler BCM counts and trigger rates for these runs. There

was no noticeable drift in the scaler rates, but the beam current was not stable during this time.. He also computed the

"scaler yield" by dividing the total triggers by the BCM counts and saw a drift similar to what is seen in the normalized

standard yield. Previous to these runs being taken, there was a flood in the hall resulting in a power outage. It's possible

that the drifting is a result of systems not being stable while these runs were taken. More details can be seen in his

slides here.

- Gave an update on his yield drift study for the 3.3 GeV setting. For this time, he looked at the target field (central

- Ryan
- Working on a data quality check for the 1.1 GeV, 2.5T setting. In general the setting looks pretty good; the maximum

spread is 10% and the minimum is 0.01%, with an average of ~2% spread. Fewer runs were taken in this setting, and the

runs were taken consecutively, resulting in smaller spreads overall. In addition, there was data taken at this setting at

3 different momentum settings, all looking at elastic, and all on the same material. Next he will try to determine the

problem for the settings with a large spread (~10%). More details can be seen in his slides here.

- Working on a data quality check for the 1.1 GeV, 2.5T setting. In general the setting looks pretty good; the maximum

## 1/7/2015

Present: Chao, Jie, JP, Min, Kalyan, Jixie, Melissa

By Phone: Yunxiu, Vince, Pengjia, Toby, Ryan, Ellie, Karl, Todd, Nilanga, Haiyan

**Feature Presentations:**

Each student gave a short update on their analysis projects since the last collaboration meeting on November 14th:

- Melissa
- Method for packing fraction analysis is complete; the packing fraction values for each material still need to be

extracted. This will probably take ~1 week, but the yield discrepancies need to be understood before it can be

completed. She is also working on data quality checks for yields, specifically the 2.2 GeV, 2.5T setting. After these

projects are complete she will work on the PbPt check using elastic asymmetries. More details can be seen in her

slides here.

- Method for packing fraction analysis is complete; the packing fraction values for each material still need to be

- Toby
- Focusing on data quality checks of yields for production runs. He has completed the 2.2 GeV 5T settings, both

longitudinal and transverse, and the 3.3 GeV setting. The 2.2 GeV 5T settings have few problems, but the 3.3 GeV

has many momentum settings where the yields drift over time. He sees a similar drift in the left and right arm data.

He is still trying to figure out the cause for the drift. In addition, he is working on tuning the Bosted model to helium

and carbon dilution runs to simulate the background for dilution analysis, with a goal of finishing the dilution analysis

by mid-march. More details can be seen in his slides here. - Information on data quality checks can be found on the wiki here.

- Focusing on data quality checks of yields for production runs. He has completed the 2.2 GeV 5T settings, both

- Jie
- Looked at the effect of changing the scattering angle on the cross section for carbon, nitrogen, helium and hydrogen.

This may help explain the discrepancy we see in the yields, but is problem not the only cause. He is also working on

finishing the dp simulation, which should be done soon, and will begin working with Min to learn more about the acceptance

study. More details can be seen in his slides here.

- Looked at the effect of changing the scattering angle on the cross section for carbon, nitrogen, helium and hydrogen.

- Chao
- Gave an update on the status of optics analysis. He is finishing up the 2nd iteration of the optics calibration; the LHRS

is complete and the RHRS has one setting to go (1.158 GeV). For the 2.2 GeV, 5T, transverse setting, there is no full

dp scan, so the longitudinal data taken at that energy setting will be used. Once the 2nd iteration of the calibration is

complete, he and Min will compile a technote on the optics study. He has also modified the geometry part of the g2psim

package to make it more configurable. More details can be seen in his slides here.

- Gave an update on the status of optics analysis. He is finishing up the 2nd iteration of the optics calibration; the LHRS

- Min
- Working on the acceptance study. In the focal plane, there seems to be good agreement between simulation and data.

In the target plane, in the θ vs φ plot, there seems to be a shift in φ for all dp settings; it seems to be a linear

correlation with dp. Next she will work on a correction for φ-target vs dp and compare target plane variables. To finish the

good septum, straight through setting it will probably take 1-2 months, then an additional 1-2 months for the other 2 septum

settings. More details can be see in her slides here.

- Working on the acceptance study. In the focal plane, there seems to be good agreement between simulation and data.

- Ryan
- He has completed a draft for a PRC paper on unpolarized He3 cross section data. He still needs to complete the He3 radiative

correction analysis for saGDH, but is waiting on updated nitrogen cross sections for saGDH. He is also working on the carbon/

nitrogen ratio for saGDH, but is waiting on updated analysis from Vince before continuing. He is also working on data quality

checks for production runs for the 1.1 GeV setting, and is starting to look into polarized radiative corrections. More details can

be seen in his slides here.

- He has completed a draft for a PRC paper on unpolarized He3 cross section data. He still needs to complete the He3 radiative

- Pengjia
- Calculated the false asymmetry by comparing the asymmetry before and after applying detector cuts. He found the false

asymmetry to be small; less than 200 ppm for all energy settings. He's also working on a data quality check for the 1.7 GeV

data set. More details can be seen in his slides here.

- Calculated the false asymmetry by comparing the asymmetry before and after applying detector cuts. He found the false

**General Discussion:**

- Yields stability/overall data quality check is the main issue at this time.
- Progress of acceptance study looks good, but could become an analysis bottleneck.