View source for Minutes July2015 to Dec2015

==12/23/2015==


Present: Min, Jixie, JP <br>
By Phone: Jie, Chao, Ryan <br>

'''Feature Presentations:'''

*Jie
Gave an update on his yield study using simulation. He studied the effect of the raster size cut on the yield. The yield changes with the area of the raster cut size as the prediction but still shows some correlation. He was suggested to study if there is some high order effect which could be corrected. After applied the raster cut, half of the runs with the yield problem could be resolved to a reasonable 3% range. More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/jie/2015_12_20_raster_plots/yields_update_2015_12_20_backup.pdf here].

*Min
Gave an update on her very preliminary calculation results of the asymmetries and structure functions of the 1.7 GeV setting. She applied dilution effects and the radiative correction to the calculation. She also provided a systematic uncertainty estimation. More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/asymmetry/12232015_g2p.pdf here].

'''General Discussion:'''

*The meeting time will be changed to '''10am Tuesday''' next year.


==12/16/2015==


Present: Min, Jixie, Melissa <br>
By Phone: Jie, Chao, Ryan, Toby, Ellie <br>

'''Feature Presentations:'''

*Toby
**Showed updated dilution results for the 2.2 GeV, 5T longitudinal setting, using the <br>most updated optics matrix.  He used a weighted average of the packing fraction <br>values for material 17; there is some drift in the yields for this setting, and therefore <br> some drift in the packing fraction.  He showed the effect of applying 4 different<br> acceptance cuts, and found that the structure of the dilution factor was stable.  The <br>dilution factor is very dependent on the packing fraction.  The method for dilution <br>analysis is complete, but until the yield drifts are completely understood (and final <br>packing fraction values can be extracted), the dilution analysis cannot be finalized.  <br>He also started looking at the transverse setting.  Before applying any acceptance cuts,<br> he found there was already good agreement between the data and the Bosted prediction.<br>  It may not be necessary to apply acceptance cuts at this setting, but it's not clear why <br>there is no acceptance effect. Ryan suggested that the acceptance effect might be <br>minimized since the scattering angle is slightly larger for this configuration.More details <br>can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/diltalk_121615.pdf here].

'''General Discussion:'''

*Chao has finished updating the optics databases for all of the LHRS settings.  He will <br>update the DB on ifarm so that Toby can replay the data.


==12/9/2015==


Present: JP, Jixie, Melissa <br>
By Phone: Vince, Jie, Chao, Toby, Ellie, Ryan, Min <br>

'''Feature Presentations:'''

*Jie
**Compiled a table summarizing the results from his yields study using simulation. In the <br>simulation, the beam position and raster size are taken into account.  He also simulated <br>the change in the acceptance from run to run.  He found he was able to reproduce the shift <br>in the yields with simulation for about 50% of the runs.  There was some question about the <br>sign of the beam position, it's possible there is a difference in the definition of the coordinate <br>systems used by Pengjia and included in the simulation.  Jie will double check with Pengjia <br>to make sure the sign is correct.  He is also summarizing his results into a technote.


==12/2/2015==


Present: JP, Min, Melissa <br>
By Phone: Vince, Jie, Chao, Toby, Ellie, Ryan, Karl, Pengjia <br>

'''Feature Presentations:'''

*Toby
**Gave an update on his acceptance and yield study using the updated optics database. He <br>looked at the helium run for one setting in the longitudinal configuration with the smallest p0, <br>since it is way past the delta the dp distribution should be fairly flat.  Without applying any <br>acceptance cuts, a sawtooth pattern is seen in the data, while the Bosted model prediction is <br>much flatter. By applying a cut on dp vs phi, he was able to find a region where the dp distribution <br>matches well with the Bosted prediction.  Ellie suggested using a minimization procedure, instead <br>of selecting the cut by eye, to improve the results further.  Karl cautioned that we don't want to bias<br> our result; we don't necessarily want the distribution to be flat, we just want it to be continuous.  JP <br>commented that it is hard to do this without a full acceptance study; Karl suggested varying phi by <br>1-2% to see the effect. More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/updatedDB_acc2.pdf here].

*Min
**Showed very preliminary results for asymmetries for the 1.7 GeV setting.  She applied loose acceptance <br>cuts on the focal plane, but will update the results later with cuts on the target plane.  She is currently <br>using a constant value for the dilution factor, and hasn't applied radiative corrections.  She will work on <br>applying these corrections next.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/asymmetry/12022015_asym.pdf here].

*Pengjia
**Showed the yield spread for different momentum settings, when different sized raster cuts are applied.  JP <br>suggested showing just one momentum setting and showing the effect of the raster cuts for each run.  He <br>also suggested starting with a setting that has stable yields, to confirm that the yields are still stable after <br>the raster cuts are applied.  More details can be seen in his slide [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/12012015/raster_cut_reliability.pdf here].

'''General Discussion:'''

*Jie will give the g2p update at the Hall A Collaboration meeting in January.


==11/25/2015==


Present: JP, Kalyan, Min, Melissa <br>
By Phone: Vince, Jie, Toby, Ryan, Karl, Chao <br>

'''Feature Presentations:'''

*Melissa
**Showed updated "very preliminary" results for g1/g2.  She tried to estimate the uncertainty on the <br>MAID model by comparing it with EG1B data for g1.  She looked at the difference between the <br>model prediction and the EG1B results for their lowest values of Q2, and assigned a constant <br>uncertainty of 40% on the MAID model.  JP suggested breaking the uncertainty up into different <br>regions, since the data in the delta region matches well with the model prediction but there is a <br>larger discrepancy for higher W.  Karl questioned whether the enhancement in the longitudinal <br>asymmetry at large W is real, this will have to be considered carefully for the final analysis. More<br> details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/g1_g2/g1g2_11_25.pdf here].

*Chao
**Gave an update on the issue with the optics database.  Toby showed last week that applying a <br>cut to phi resulted in structure in the dp distribution.  Chao scanned the optics matrix to find the <br>element causing this issue, and found that the x^2 and x^3 elements of the phi tensor were causing <br>the over-fitting of the matrix.  He showed the dp distribution without these two elements and the <br>unusual structure seen when a cut is placed on phi was gone. JP/Vince commented that it is troubling <br>that the x^2 dependence is causing a problem.  They will discuss is more during the optics meeting. <br>More details can be seen in Chao's slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20151125/Chao_WeeklyMeeting_11252015.pdf here].


==11/18/2015==


Present: Min, JP, Melissa <br>
By Phone: Vince, Jie, Toby, Ryan, Karl, Ellie, Chao <br>

'''Feature Presentations:'''

*Toby
**Gave an update on using the newest optics reconstruction matrix.  He showed an example of one<br>run, located way past the delta resonance, so the dp distribution should be fairly flat.  He tried <br>applying the phi cuts he has shown previously to remove acceptance and scattering angle effects,<br> but making a cut on phi revealed a strange structure in the dp spectrum.  He showed the phi vs. dp <br>reconstruction both at the septum plane (i.e. without target effects) and after being reconstructed <br>back to the target.  The strange structure exists in both distributions, so it is probably not a target <br>field effect.  He is currently working with Chao to resolve this; Chao is doing a comparison of the <br>old and new databases to determine which element is causing this effect.  More details can be seen <br>in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/updatedDB_acc.pdf here].
**He has updated the target polarization uncertainties in the mysql.  He will send out an email to the <br>collaboration with more details about this.

*Pengjia
**Showed the results of a study to check the effect of applying raster cuts on the yield.  He showed <br>several examples of groups of runs where the yield was drifting.  In some cases, the raster cuts <br>got rid of the drift in the yields, or at least made the drift significantly smaller. For one setting, the <br>raster cut did not resolve the change in yields, so there is likely something else causing the drift.  <br>JP suggested doing some sanity check to confirm that the method works, such as starting with <br>a setting that has consistent yields and applying the raster cuts to see the effect.  More details <br>can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/11182015/raster_cut_library.pdf here].


==11/11/2015==

Present: JP, Melissa <br>
By Phone: Vince, Toby, Ryan, Chao, Min, Karl <br>

'''Feature Presentations:'''

*Melissa
**Showed "very preliminary" results for g1/g2.  She is using a dilution factor determined from <br>the Bosted model, and using the average value of the packing fraction for each material. For the <br>radiative corrections, she is using a model prediction for the asymmetry where the polarized piece <br>comes from the MAID model and the unpolarized cross section is the Bosted model.  The difference <br>between the model prediction ''without'' radiative effects and the model prediction ''with'' radiative effects <br>is used as an additive correction to the asymmetry.  There was some concern over the large effect <br>of the R.C. at high W.  Ryan commented that the final R.C will be large, so they will need to be done<br>very carefully.  A cross check with EG1 data will be helpful.  In order to match the kinematics of the<br>longitudinal and transverse settings, she scaled the asymmetry by the ratio of the unpolarized cross<br>section at the two different kinematic settings.  JP commented that this is not valid, and suggested <br>using model input for the asymmetries, instead of trying to match the kinematics between the data<br>taken in the longitudinal and transverse settings.  Next she will work on updating systematics for the <br>asymmetries and structure functions.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/g1_g2/g1g2_11_10.pdf here].

*Chao
**Gave an update on optics analysis.  Comparisons of the data with simulation revealed that there <br>were offsets in the reconstructed kinematics.  This correction could be fixed using a finely-tuned<br>fitting procedure.  The longitudinal setting has already been fixed, today he showed the results <br>for the 1.7 GeV, transverse setting.  He showed the results of doing the calibration; most of the <br>holes near the center have a difference of less than 0.5 MeV between data and simulation.  He is <br>currently working on finishing the other transverse settings.  More details can be found in his <br>slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20151111/Chao_WeeklyMeeting_11112015.pdf here].


==11/4/2015==

Present: Min, Jixie, Melissa <br>
By Phone: Chao, Vince, Toby, Ryan, Alex, Karl <br>

'''Feature Presentations:'''

*Ryan
**Gave an overview of g2p elastic tail and angle reconstruction.  The elastic tail is highly sensitive<br>to the scattering angle, so his goal is to develop a method to include these variations in the elastic <br>tail subtraction.  He started by looking at the simplest case, helium dilution runs.  Using his method, <br>he looked at the elastic tail across a series of different momentum settings, andfound that the tail <br>was not continuous, which he assumed was due to variations in the scattering angle.  He did a pseudo-<br>simulation using the Bosted model and varying the scattering angle to confirm the "jump" seen in the <br>elastic tail.  However, when he looked at the normalized yields for this same data, he did not see the <br>same discontinuity. This suggests there might be a problem with the scattering angle reconstruction. <br>Chao mentioned that the reconstruction code has been updated (at least for the longitudinal setting) <br>the last replay; he will provide this code to Toby so he can do another replay.  Details of the method <br>and cuts used can be seen in his slides [https://userweb.jlab.org/~rbziel/Weekly_Meetings/g2p_11415.pdf here].

*Jie
**Gave an update on his simulation study to understand the beam dependence on the yields.  He used the <br>beam position and raster size as input to g2psim, and included radiative and acceptance effects. In <br>general, he found there is a much larger dependence on the beam position using a rastered beam <br>compared to a point beam.  He summarized his results for various different momentum settings and found <br>that, in some cases, his results matched what we saw in the data, but not in all cases.  Vince suggested <br>testing how much the beam position would need to be changed (for a given momentum setting) to match <br>the simulation with the data.  The details of his study can be found in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/jie/yields_update_2015_11_04.pdf here], he is also compiling <br>a technote on this analysis.


==10/28/2015==

Present: JP, Min, Kalyan, Melissa <br>
By Phone: Alexandre, Vince, Jie, Ryan, Toby, Ellie, Karl <br>

'''Feature Presentations:'''

*Melissa
**Gave an overview of contributions to the systematic uncertainty for the asymmetry.  There are <br>several contributions which are very small (<0.1%) including the charge asymmetry, livetime <br>asymmetry and pion contamination. The target polarization uncertainty was surprisingly small.<br>The values currently listed in the mysql DB for this quantity appear to be incorrect; Toby will<br>update these values.  For the contribution from the dilution factor and radiative corrections, she<br> assigned 5% as a placeholder until a better estimate could be made.  Karl/Ryan commented <br>that 5% uncertainty on the radiative corrections may be too small, 10-15% would be a better<br>estimate at this point.  The tables shown for the longitudinal and transverse settings represent<br>an overestimate of the uncertainty; in the long-term this could be looked at more carefully to<br>decrease the overall systematic uncertainty. More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/Asym/Asym_10_28.pdf here].


'''General Discussion:'''
*There has been a request from the Hall A CC to have an update on g2p at the winter collaboration <br>meeting.  Ellie and Jie are both interested, they will check the dates of the meetings to see if they <br>are available.


==10/14/2015==

Present: JP, Min, Jixie, Kalyan, Melissa <br>
By Phone: Alexandre, Vince, Jie, Todd, Haiyan, Pengjia, Karl, Toby, Ryan, Xiaochao, Nilanga, Chao <br>

'''Feature Presentations:'''

*Karl
**Went over the goals of the collaboration meeting:
***Discuss physics topics for students graduating soon
***Identify major remaining physics topics
***Establish a policy for the release of plots
**Suggested a policy for collaboration release of physics results, including specific requirements <br>for moving from "very preliminary" to "preliminary" to "final" results.  Details can be seen in his <br>slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/slifer/2015/slifer_g2pcollab_101415.pdf here], and they will be circulated to the g2p mailing list with some small modifications.

*Pengjia
**Will be defending next Wednesday (Oct. 21st).  He showed a "very preliminary" physics <br> asymmetry for the 2.2 GeV 5T setting (both longitudinal and transverse) with a comparison to <br> a radiated model prediction.  He used this result to extract a preliminary value for g1 and g2,<br>using the Bosted model as input for the cross section.  This result does not include radiative<br>corrections.  More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/graduate_2015.pdf here].

*Min
**Showed an overview of status of acceptance study; the overall goal is to match the simulation<br>results to the data, so that the simulation can be used to calculate the acceptance.  She<br>described the various attempts to correct the simulation results to match the data.  Adjusting the<br> apertures based on optics data works well, but the other two methods (adjusting the apertures <br>based on dilution data and adjusting the acceptance cuts after removing the Mott XS from the <br>data) only work well around the target run.  In addition, the simulated yield has a different shape <br>than the data, so the simulation must be adjusted to match this.  She will go through this process <br>again once Chao is finished with the optics correction.  She is planning to graduate next January, <br>and would like to include a cross section difference for one kinematic setting in her thesis.  More <br>details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/collab/10142015_collab_Min.pdf here].

*Melissa
**Gave an overview of packing fraction analysis.  Of the 5 total settings there are 3 that are reasonably <br>stable with some small drifts in yields, and 2 settings that have larger problems.  During the 2.2 GeV, <br>2.5T setting the raster size was changed, which could possibly be corrected by applying a cut on<br>the center of the raster.  The 1.1 GeV setting has large variations in the elastic yields; some of these <br>runs were taken at low currents (8-20nA), so the beam position information is not reliable.  She also <br>showed the results of a pion asymmetry study.  Although the pion asymmetry is very large, the <br>residual pion contamination after PID cuts are applied is very small, so the correction is not large. <br>She has scheduled her defense for December 4th, and would like to include a preliminary result <br>for g1,g2 in her thesis.  To include a fully corrected result she will need preliminary radiative <br>corrections and dilution results, which she will work with Ryan and Toby to get.  More details can <br>be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/talks/MCummings_CollabMeeting_20151013.pdf here].

*Toby
**Showed a preliminary dilution value that he had shown previously, and described several issues <br>with the result.  These issues include that the model assumes a constant scattering angle, and does <br>not include some acceptance effects that are seen in the data.  He has been working on correcting <br>these issues by incorporating a scattering angle from data into the model and applying acceptance <br>cuts to the data to remove acceptance effects.  He showed the result of applying the acceptance cut <br>across all kinematics for the 5T longitudinal setting; including the acceptance cut gets rid of some of <br>the strange structure seen in the yields, such as the sawtooth pattern seen at high nu.  The downside <br>of this method is that the acceptance cut causes a large loss in statistics.  He is currently working on <br>trying to re-incorporate some of the lost statistics from the acceptance cut.  More details can be seen <br>in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/collaborationupdate.pdf here].

*Jie
**Gave an update on his study using simulation to understand the drifts in yields.  He looked at several <br>momentum settings in the 2.2 GeV, 2.5T setting that had variation in the yields.  He varied the beam <br>position along the horizontal and vertical directions using both a point beam and rastered beam.  He <br>generated a table summarizing the results using a Monte Carlo method, cross checking with averaged <br>acceptance.  His results were similar to his previous test using only the Bosted model; some results <br>agree with what was seen in the data, some have the opposite sign.  He will write up his method and <br>results in a technote.  He is planning to graduate in summer of 2016.  More details can be seen in his <br>slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/jie/2015_10_09_simulation/2015_10_14_collaboration_meeting.pdf here].

*Chao
**Gave an update on optics analysis.  Comparisons of simulation with data revealed that the <br>reconstructed kinematics have some offset. A correction was added to the reconstructed kinematics <br>(vs beam position) and he is currently working on merging this correction into the database with a <br>finely tuned fitting procedure.  The problem is that the horizontal beam position shifts the dp <br>reconstruction.  He showed the results of applying the new optics matrix which corrects for this shift. <br> He will apply this method to all kinematic settings and study the uncertainty of the optics.  Once this <br>is done he will write a NIM paper for g2p optics. He plans to graduate in February/March and would <br>like to include the cross section difference for one or two settings with tuned acceptance and model <br>inputs to obtain a physics result.  More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20151014/Chao_WeeklyMeeting_10142015.pdf here].

*Ryan
**Since the last collaboration meeting he has improved the model tuning procedure by using an <br>automated procedure instead of tuning the model by hand.  He has tested this using carbon data <br>from the UVa database and will run it for nitrogen data once he has updated saGDH N2 cross <br>sections from Vince. He also updated the systematic error study for inelastic unpolarized radiative <br>corrections by separating the integration/iteration error from the interpolation error.  He has the <br>polarized radiative correction code running, and has completed a systematic uncertainty study on <br>the inelastic radiative corrections.  He is also working on evaluating physics quantities using models, <br>specifically the hyperfine splitting calculation.  He plans to graduate in summer/fall of 2016 and would <br>like to include full radiative corrections on g2p polarized cross sections in his thesis.  More details can<br> be seen in his slides [https://userweb.jlab.org/~rbziel/Weekly_Meetings/g2p_101315.pdf here].


'''General Discussion:'''
*Important to do cross checks of results between students (ex. asymmetries)
*Must be careful to include details of any physics results (whether radiative corrections are included,<br>dilution factor, model input instead of data, etc.)


==9/30/2015==

Present: JP, Kalyan, Jixie, Min <br>
By Phone: Pengjia, Toby, Ryan, Vince, Jie, Chao, Karl <br>

'''Feature Presentations:'''

*Toby
**Showed an update on his acceptance/dilution study.  He continued with the study of the cuts on the <br>reconstructed &theta; vs &phi; distribution at the target.  He scanned several different &theta; vs &phi; cuts to see how <br>the acceptance effects the scattering angle and resulting yields. It turns out the small angle region is <br>influenced by the acceptance a lot but the center region gives a good agreement to the model. He <br>also loosed the center cut to &theta; ~ ± 40 mrad and &phi; ~ ± 10 mrad, the data still agreed with the model. <br>He then recalculated the dilution factor for 2.254 GeV longitudinal setting. More details on his <br>acceptance cuts and method can be found in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/diltalk_093015.pdf here].


==9/23/2015==

Present: JP, Kalyan, Min, Melissa <br>
By Phone: Pengjia, Toby, Ryan, Vince, Jie, Chao, Karl <br>

'''Feature Presentations:'''

*Pengjia
**Showed the structure of his code to create a radiated model of the MAID and Bosted <br>models, based on Ryan's radiative correction work.  He showed the breakdown of the <br>different components to the radiated model (internal, external and elastic tail).  He used<br> the radiated model to calculate an asymmetry and cross section differences, which he <br> compared with data.  Finally, he combined asymmetries calculated from data with the <br>Bosted model cross section to calculate a value for g1 and g2, which he compared <br>with values of g1/g2 determined using the MAID model.  More details can be seen in <br>his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/09222015/asym_g12.pdf here].

*Min
**Gave an update on acceptance studies.  Last time, she found that by adjusting the<br>apertures in the simulation she could match the data/simulation results for dp (for one <br>momentum setting). She tested this method for other momentum settings and found <br>that it did not work as well.  She also showed a comparison of the yield spectrum for<br>helium runs in the 1.7 GeV setting; the shape of the two distributions are quite <br>different, which might suggest that the model deviates significantly form the data.  In<br>the simulated results she saw several "bumps" in the yield, which correspond to the<br>location of resonances.  She will update these results when the optics corrections are <br>complete.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/acceptance/09232015.pdf here].


'''General Discussion:'''
*We would like to schedule another collaboration meeting in the near future.  The <br>meeting is tentatively scheduled for Wednesday October 14th during our usual meeting <br>time (10am-12pm).
*The purpose of the meeting is to:
**Discuss student thesis topics (particularly any updates from the last <br>collaboration meeting)
**Establish a procedure for presentation/release of collaboration plots/results
*Each student should prepare a few slides (~5 minutes) to show analysis progress


==9/16/2015==

Present: Min, JP, Jixie, Melissa <br>
By Phone: Vince, Jie, Chao, Pengjia, Alex, Ryan, Toby, Ellie <br>

'''Feature Presentations:'''

*Jie
**Gave an update on his simulation study using MUDIFI.  Last time he showed that he uses a<br>least square fit to go from endplane to endplane.  For the g2p simulation, it's used to go from<br>(x,y,&theta;,&phi;,dp) in the target plane to (x,y,&theta;,&phi;,dp) in the focal plane. SNAKE produces (x,y,&theta;,&phi;,dp)<br>data in both planes and a polynomial function can be used to determine the relationship<br>between the two planes.  First the maximum power for each variable is determined, symmetries<br>are included (if applicable), then the priority of each function must be determined. A Gram-<br>Schmidt transformation is used to do the minimization, and the final step is to decide whether<br>to accept the function f<sub>i</sub>.  JP/Vince asked how the priority for each function is chosen, since the<br>choice could introduce some ordering dependence.  Since we don't have enough data to<br>constrain, JP suggested introducing some physics.  He also suggested getting rid of dp in the<br>focal plane (since it is really just dp in the target plane) so that there is one less dimension to<br>fit.  More details can be seen in Jie's slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/jie/2015_09_15_simulation/2015_09_15_Simulation_update.pdf here].

*Ryan
**Gave an update on evaluating physics quantities using models.  His goal was to do a<br>comparison of the hyperfine splitting calculation in C. Carlson's 2006/2008 paper using MAID<br>and Hall B models.  Previously he showed a comparison using the Hall B model, where he saw<br>good agreement between his calculation and Carlson's result.  For this time he also did a<br>comparison using the MAID model.  There is a large difference between the two models for the<br>value of the integrand of &Delta;<sub>2</sub>, over the full range of Q<sup>2</sup>. JP suggested comparing the value of g2<br>between the two models, instead of &Delta;<sub>2</sub>.  More details can be seen in his slides [https://userweb.jlab.org/~rbziel/Weekly_Meetings/g2p_091615.pdf here].


'''General Discussion:'''
*Toby
**Last time he showed a comparison between data and cross sections calculated using the<br>Bosted model (seen [http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/acceptance_scatangle.pdf here]).  The model was not continuous for one momentum setting.  It <br>turns out the BPM was not recording for that setting, so he re-calculated the XS using<br>values for &theta; and &phi; from surrounding runs and found that the result was now consistent with<br>the neighboring momentum settings.  He also played around with the acceptance cuts he<br>showed, and found a region using the "extreme cut" where the sawtooth pattern go away<br>and the yields become flat.


==9/9/2015==

Present: Min, JP, Jixie, Melissa <br>
By Phone: Vince, Jie, Chao, Ryan, Toby, Karl <br>

'''Feature Presentations:'''

*Jie
**Gave an overview of his simulation study using SNAKE and MUDIFI.  Using SNAKE, the<br>idea is to trace the particle trajectory in a series of free boxes, then an endplane can be<br>defined within the free boxes to determine the trajectory from endplane to endplane. In<br>g2psim this will be used to trace the trajectory between the sieve plane and the focal<br>plane.  There are two important issues to consider, the magnetic field in the free boxes,<br>and the transportation method between the two endplanes.  A Runge-Kutta method is <br>used for the transportation, and an adaptive method (which changes with field size and<br>gradient) is used to determine the step size.  For the magnetic field, a 3D Tosca field<br>map is used for the septum, and an analytic model is used for the spectrometer magnets. <br>  He also described the method for choosing the transportation function to determine the <br>25 matrix terms for each kinematic variable (up to 4th order).  JP asked how to determine <br>which order is more important, and cautioned to be careful not to over fit the data.  He and<br>Jie will discuss this more offline.  More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/jie/2015_09_09_simulation_update/2015_09_09_Simulation%20update.pdf here].

*Chao
**Gave an optics status update.  The current problem is that changes in the horizontal beam<br>position changes the dp reconstruction.  This was fixed by applying a linear correction in<br>BPM-x and BPM-y to &theta;, &phi; and dp.  Applying this correction made the data/simulation <br>agree, but it is still necessary to understand the correction.  Using dp scan data with beam<br>position scan data still leaves an offset in &theta; and &phi;.  He also tried two other methods to do<br>the calibration.  Method 1 uses the central hole to determine a constant offset, uses the<br>horizontal beam position scan data to determine all the matrix elements related to y, and<br>uses dp scan runs to fit the additional matrix elements.  Method 2 uses beam position<br>scan data to determine the matrix elements that are not related to x<sub>fp</sub>.  Neither of these<br>methods works very well, even when the fitting order of the matrix elements are adjusted.<br>He will continue to check more fitting procedures.  More details can be found in his slides <br>[http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20150909/Chao_WeeklyMeeting_09092015.pdf here].

*Pengjia
**Gave an update on comparing asymmetries from data/models for the 2.2 GeV, 5T <br>longitudinal setting.  He first radiated the Bosted model XS; he showed the break down of<br>the different radiative corrections including the elastic tail (internal and external) and the<br>radiative effects from inelastic scattering.  He showed a comparison with Ryan's results<br>which had good agreement, though there is some deviation when he used different proton<br>form factors (provided by Moshe).  In order to compare the asymmetry result with data, he<br>needs to do the radiative corrections for the polarized MAID model, which he will continue<br>to work on.  More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/09092015/asym_20150909.pdf here].


==9/2/2015==

Present: JP, Min, Kalyan, Melissa <br>
By Phone: Vince, Alex, Jie, Chao, Ryan, Toby, Ellie, Karl <br>

'''Feature Presentations:'''

*Melissa
**Gave an update on determining the correction to asymmetries for pion contamination.  <br>Previously, she used the 2.2 GeV, 5T longitudinal setting as an example, for this time <br>she used the transverse configuration, as it has the largest pion contamination.  She <br>modified her method for differentiating between electrons and pions; previously she <br>identified pions as events that do not trigger the Cherenkov, this time she used a cut <br>below the single photo electron peak to select pions.  JP suggested using stricter cuts <br>for selecting pions when determining the pion asymmetry, to make sure it is a clean <br>sample.  Although this setting has a large pion contamination, the number of pion<br> events that remain after lead glass cuts are applied is low, so the overall correction <br>is very small.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/Pion/PionAsym_08_31.pdf here].

*Min
**Gave an update on acceptance tuning.  She looked at results from the 1.7 GeV setting <br>and found that the height of the peaks is different between data and simulation for <br>different holes in the sieve slit.  Looking at the &theta; vs &phi; distribution, she noticed on the<br>positive &theta; side there were fewer holes in simulation than in data, and on the negative <br> &theta; side there were more holes in simulation compared to data.  She tried adjusting the<br>aperture cuts by changing the physical location of the Q1, Q2 and septum planes to <br>uncover events that were previously blocked.  After the adjustment, the data/simulation <br>results agreed much better.  She showed a similar procedure for a dilution run and was<br>able to see better agreement between data/simulation after the aperture adjustments <br>were made.  JP/Vince suggested being cautious when adjusting the aperture cuts, as<br> the apertures should be well defined already.  They suggested adjusting things that are<br> less well known first, such as the septum field. More details can be seen in her slides<br>[http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/acceptance/09022015.pdf here].

*Ryan
**Working on generating physics quantities using models to eventually compare with data.<br>So far, he has incorporated the MAID2007 and Hall B models.  He calculated the <br>hyperfine splitting quantities &Delta;1 and &Delta;2 using the Hall B EG1 model and compared them<br> to the results in C. Carlson's 2006 paper.  His results agreed very well, he only saw a<br> small deviation for &Delta;2 in the lowest Q<sup>2</sup> range.  To check that nothing strange was going <br> on in the low Q<sup>2</sup> range, he calculated &Delta;2 for various small ranges around zero, to ensure<br>that the integration method is well behaved.  The Carlson paper also includes calculations <br>using the Simula model, he would like to do a check of this method as well, but is unsure <br>where to find the Simula model.  Since he was able to reproduce the hyperfine splitting <br>result, he is confident in his calculation method; it should be straightforward to calculate <br>other spin structure function integrals.  More details on his method can be found in his <br>slides [https://userweb.jlab.org/~rbziel/Weekly_Meetings/g2p_09215.pdf here].

*Toby
**Showed an update on his acceptance/scattering angle study.  He applied cuts of varying<br>strictness on the reconstructed &theta; vs &phi; distribution at the target.  He then calculated the <br>real scattering angle before/after scattering for each cut.  Using these scattering angles <br> as inputs for the Bosted model he was able to calculate the XS and compare it to data.<br>He showed the comparison between data/Bosted model for each acceptance cut.  He <br>hoped that the more extreme acceptance cuts on the data would flatten the yield <br>distribution, but it actually just enhanced the sawtooth pattern.  JP questioned why the <br>Bosted prediction for one momentum bin (the green one) was so much lower than the<br> neighboring ones.  Toby said the only input that goes into the model is the scattering angle,<br> so it must be considerably different for this momentum setting. He also commented that<br>the model prediction is scaled to match the data, so the scaling enhances the jump in<br>the model.  More details on his acceptance cuts and method can be found in his slides<br>[http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/acceptance_scatangle.pdf here].


==8/26/2015==

Present: Jixie, Alex, Min, Kalyan, JP, Chao, Melissa <br>
By Phone: Vince, Jie, Ryan, Toby, Karl <br>

'''Feature Presentations:'''

*Chao
**Gave an update on optics studies.  Previously, he showed a comparison between optics data <br>and simulation where he cut on individual sieve holes.  For this week, he included a cut on y <br>and dp that he didn't include previously.  In the plot of &theta; vs y, JP pointed out that there were <br>additional tails on some sieve holes, Chao says these tails go away when the &phi; cut is <br>included. He showed a 2D plot of &theta; vs &phi;, along with the 1D projection of these variables.  JP<br>commented that the data is more smeared out than the simulation, so it is difficult to compare <br>data with simulation results.  It's possible that some of the smearing comes from events that <br>punch through the sieve slit and are not being cut out (the simulation currently assumes these <br>events are stopped by the sieve).  Looking at the center hole, the data/simulation agree <br>reasonably well, but the agreement gets worse for holes on the edge of the sieve slit.  Chao<br>will continue to look into this.  More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20150826/Chao_WeeklyMeeting_08262015.pdf here].
**He also mentioned that ~1 month ago he reported a problem with the [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20150708/Chao_WeeklyMeeting_07082015.pdf BPM simulation].<br>He found it was actually a problem with the simulation itself, which has now been updated <br> to correct the problem.


'''General Discussion:'''
*Pengjia has submitted his BPM NIM paper.


==8/19/2015==

Present: Chao, Jie, Min, Jixie, Kalyan, Melissa <br>
By Phone: Vince, Ryan, Toby, Karl, Pengjia <br>

'''Feature Presentations:'''

*Melissa
**Showed a method for determining the correction to asymmetries for pion contamination. The <br>measured asymmetry can be written in terms of the electron and pion asymmetries, multiplied <br>by the fraction of events of that particle.  To differentiate between electrons and pions, a cut is<br>made on events that trigger the Cherenkov detector (electrons) and events that do not trigger<br>the Cherenkov (pions).  Kalyan commented that this method of selecting pions may be too blind,<br> and suggested making a 2D cut on E/p vs Cherenkov instead.  To determine the fraction of<br>electrons/pions, "good event" cuts are applied (loose acceptance cuts, single track events and<br>pion rejector cuts) and the number of electrons/pions that survive the cuts are counted.  Although<br>the pion asymmetry is much larger than the electron asymmetry, the number of residual pions<br>is very small, so the overall correction is small.  For this example the 2.2 GeV, 5T longitudinal<br>setting was used, but the correction may be larger for the 5T transverse settings, which she <br>will do next.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/Pion/PionAsym_08_19.pdf here].

*Min 
**Gave an update on acceptance studies.  She used the 1.7 GeV, 2.5T setting (with 3rd septum <br>configuration) to test the acceptance cuts using the yield ratio of data to simulation.  She <br>showed a plot of target &theta; vs &phi; for both data and simulation, and applied a graphical cut to get<br>rid of events on the edges.  Looking at the resulting 1D plots shows reasonable agreement<br>between &theta; and &phi;.  Vince pointed out that there seems to be some additional structure in the<br>dp distribution, and suggested expanding the &delta; cut out to &pm; 4% to see if there is more <br>hidden structure.  Min said she has seen this in other settings as well and will look into it.  She<br>also tried adjusting the cuts on &theta;, &phi; and &delta; and looked at how the ratio of yields from data/<br>simulation changed; the ratio starts to drop off as the cuts get wider.  Next she will work on<br>calculating the yields using W bins.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/acceptance/08192015.pdf here].


==8/12/2015==

Present: Jie, Chao, Kalyan, Min, Jixie, Melissa <br>
By Phone: Pengjia, Ryan, Toby, Karl, Ellie, Alex <br>

'''Feature Presentations:'''

*Chao
**Gave an update on optics studies.  Previously, Min has shown comparisons of data/simulation<br> where she looked at different columns of sieve holes.  Chao continued this study by cutting on <br>individual sieve holes.  Looking at the center sieve hole, the &phi; distribution looks reasonable, but <br>the &theta; distribution has some issues.  Specifically, the data is wider on the left side (negative &theta;) <br>than the simulation.  He showed results for several sieve holes, including one that was in the <br>row as the center hole, but in the farthest column on the left side of the sieve slit.  In this case <br>the &theta; distribution had similar width for both data and simulation, but the height of the simulation<br> distribution was larger.  This may suggest that something is wrong with the acceptance in the <br>simulation.  The problem with &theta; seems to be systematic; for most of the sieve holes there are <br>more events on the left side of the data than in the simulation.  Kalyan asked whether the <br>thickness of the sieve slit was taken into account, Chao says it was included in the simulation <br>and optimizer.  Jixie suggested comparing the dp distributions for each sieve hole. More details <br>can be seen in Chao's slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20150812/Chao_WeeklyMeeting_08122015.pdf here].


==7/29/2015==

Present: Chao, Min, JP, Melissa <br>
By Phone: Alex, Pengjia <br>

'''Feature Presentations:'''

*Min
**Gave an update on acceptance studies.  Previously, Chao showed a correction for the beam<br>position in g2psim.  Min checked this correction by looking at optics runs.  On the data she <br>applied a cut on dp (to select elastic events) and on the focal plane (to get rid of junk <br> events).  On the simulation, only an aperture cut was applied.  She compared the data and<br>simulation results for &theta; and &phi; both before and after drifting from the sieve slit to the target.<br>The effect of the target field on &theta; is not very large, though there is a shift down in angle.<br>However, the shape of the &theta; distribution is different between data and simulation.  The<br>effect of the target field on &phi; is significant; it's not clear why this is the case.  To try to<br>determine the reason for the discrepancies, she applied a cut along each column of sieve<br>slit holes in &phi; and compared the data and simulation results for &theta;,&phi; and dp.  Column 2<br>seems to have the best agreement (for dp), columns 1 and 2 have a shift while columns<br>3-6 have a difference in width.  Chao suggested also cutting along each row of sieve slit<br>holes to better study the effect of the drift on &phi;.  Min will work on identifying the differences<br>between data and simulation, specifically whether the boundaries are different between<br>the two.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/07272015/07292015.pdf here].


==7/22/2015==

Present: JP, Min, Jixie, Melissa <br>
By Phone: Alex, Ryan, Toby, Ellie, Karl, Pengjia <br>

'''Feature Presentations:'''

*Melissa
**Showed various checks for packing fraction runs that have large discrepancies, including <br>beam position, beam energy, spectrometer central momentum setting, septum current, and<br>the location of anneals.  For both the 2.2 and 1.1 GeV settings, there are significant shifts in <br>the x and y beam position (2mm or more).  She showed the average and standard deviation <br> of epics readings for the LHRS septum current for each run.  For several runs, there are small <br>number of epics readings that drift from the average current.  It's not clear if this is the <br>current is actually changing or if it is a readout problem.  She will look into the yield vs time<br>to try to determine this.  JP suggested making a raster cut on the center of the beam to get<br> rid of any potential scraping.  Pengjia says there is a way to cut on raster size, but there will<br> be an additional BPM uncertainty.  Melissa will discuss this more offline with Pengjia.  More<br>details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/Elastic/Elastic_07_22.pdf here].

*Pengjia
**Gave an update on calculating asymmetries/cross sections using the MAID model.  This time, <br>he used the Bosted model to get the unpolarized cross sections for the proton.  Using unradiated <br>cross section models, he still saw a factor of 6 difference when compared with data.  Ryan also <br>calculated asymmetries using unradiated models, and the results were consistent with Pengjia's, <br>so it seems the radiative effects contribute significantly to the factor of 6 difference.  More <br>details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/07222015/MAID_update_20150722.pdf here].


==7/15/2015==

Present: JP, Min, Chao, Jixie, Melissa <br>
By Phone: Ryan, Toby, Ellie, Karl, Vince, Pengjia <br>

'''Feature Presentations:'''

*Chao
**Doing a check of simulation/data results for the 2.2 GeV, longitudinal setting in order <br>to check that the optics calibration matrix will work for all momentum settings (not <br>just elastic).  He showed a comparison for several different momentum settings.<br>  Vince suggested simulating more events, since the statistics start to drop off for <br>lower momentum settings.  In general, the data is broader then the result from<br>simulation.  In the &theta; distribution, the "tails" seen in the data are not symmetric. JP<br> suggested carefully selecting cuts to remove any background from the data to <br>ensure a good comparison with simulation.  More details can be seen in his <br>slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20150715/Chao_WeeklyMeeting_07152015.pdf here].

*Pengjia
**More discussion on Pengjia's [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/07082015/MAID_update_20150708.pdf slides] from last week.  Pengjia has seen a factor of <br>six difference between his asymmetry calculated from data compared to the <br>asymmetry determined using the MAID model.  Ryan did not see this factor of six <br>in his study, but in addition to using radiated XS models, he is using a combination <br>of the Bosted and MAID models (Pengjia is using just the MAID model).  Karl <br>commented that the XS from MAID is not very good for our kinematics.  Pengjia <br>also showed an estimated NH3 XS.  JP commented that it was surprising that the <br>delta peak wasn't visible.  Pengjia will look into this for next time.


==7/8/2015==

Present: JP, Min, Chao, Melissa <br>
By Phone: Toby, Vince <br>

'''Feature Presentations:'''

*Pengjia
**Working on comparing asymmetries determined from MAID to those calculated from <br>data; last time he showed a factor of 6 difference between the two results.  This <br>time, he showed two different methods for calculating the differential XS from the <br>virtual photon XS.  The first method calculates the differential XS directly from <br>the virtual photon XS, and the second method calculates the differential XS using <br>F<sub>1</sub>, F<sub>2</sub>, g<sub>1</sub> and g<sub>2</sub>.  He also showed two different methods for calculating the <br>asymmetry from the virtual photon XS; one method uses A<sub>1</sub> and A<sub>2</sub>, while the other <br>is calculated from the virtual photon XS.  There was a discrepancy between the two <br>methods, but he found that there was actually a mistake in one of the equations. <br>While the results from the two methods agree with each other, they are still a factor <br>of 6 larger than the results from data.  More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/07082015/MAID_update_20150708.pdf here].

*Toby
**Showed an update of his scattering angle study for model reconstruction.  He is <br>using BPM and optics data to determine the scattering angle.  He showed an example <br>of the calculated scattering angle; JP questioned whether the range was really from <br>0-16 deg.  Toby pointed out that there really aren't many events at the boundaries, <br>the distribution is peaked between 4.5-9 deg.  He then uses this scattering angle as<br>input for the Bosted model.  At low &nu; the simulation results match reasonably well<br>with the data, but the discrepancy starts to become larger as &nu; increases. He took <br>a closer look at the reconstructed scattering angle for each central momentum, and <br>found that a second peak is visible for smaller values of p0. If this second peak was <br>real, the simulation should recreate the data, so this could suggest a problem with <br>the reconstruction.  However, Chao pointed out that, for this energy setting, the optics <br>calibration is not complete (currently using the longitudinal optics matrix), so the <br>reconstructed &theta; and &phi; may not be correct yet.  For next time, Toby will test this <br>analysis using the 1.7 GeV setting, which has calibrated optics.  More details can <br> be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/toby/g2p%20meetings/scatangle_talk.pdf here].

*Melissa
**Gave a summary of packing fraction analysis.  For some settings, the variation in the <br>yields (and P<sub>f</sub>) is due to a fluctuation in the beam position. This will hopefully be<br>resolved by Jie's beam position/acceptance study.  For the 2.2 GeV, 2.5T transverse setting<br> there is also an issue of a shift in &nu; between runs.  For the 1.1 GeV, 2.5T transverse <br>setting, the yields vary depending on the beam current used for each run.  JP suggested <br>a few things to check as the possible cause for these discrepancies including drift in the <br>BCM calibration, anneals of the target material, drift in field (HRS or septum), etc.  She <br>has posted a technote draft [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/Technotes/PF_technote_draft.pdf here], feedback would be appreciated.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/melissac/Elastic/Elastic_07_08.pdf here].

*Min
**Gave an update on her acceptance study.  Last time, she showed that the simulation results <br>are more narrow than the data.  For this time, she tried cutting on just the center hole, <br>and found that the dp distribution from data better matched the simulation.  JP suggested <br>cutting on each of the different holes individually to determine which hole causes the dp <br>distribution to be wider.  More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/07082015/07062015.pdf here].

*Chao
**Gave an update on optics analysis, currently checking whether simulation results match <br>with data, starting with optics data.  First he checked whether the geometry was correct, <br>some deviation was found in the recent optics meeting.  The position of the BPMs is hard-<br>coded into the simulation, and the drifting algorithm is used to drift the electrons <br>backward to the BPMs to simulate the readout of BPM A and B.  He compared the &theta; <br>distribution before drifting in the target field using a fixed and non-fixed BPM location. <br>The results using a fixed BPM match well with the data.  He is currently working on <br>checking this for all momentum settings, and will follow up with Pengjia about the BPM <br>problem. More details can be seen in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/chao/20150708/Chao_WeeklyMeeting_07082015.pdf here].


==7/1/2015==

Present: Chao, Vince, Jie, JP, Min, Melissa <br>
By Phone: Pengjia <br>

'''Feature Presentations:'''

*Min
**Gave an update on acceptance studies.  She showed a summary table of each energy <br>setting and septum configuration.  There hasn't been any problems in the 2.2 GeV, 5T <br>longitudinal and 1.7 GeV, 2.5T transverse settings, but some settings (2.2 GeV, 2.5T <br> transverse, 1.2 GeV, 2.5T transverse and 2.2 GeV, 5T transverse) have a discrepancy <br>in the focal plane between data and simulation.  She compared the data divided by the <br>Mott XS to the simulation results without XS, and found that the simulation results <br>of &theta; and &phi; are narrower then the data.  JP commented that, while including the Mott<br>XS might change the shape of the distribution, it shouldn't change the boundary.  She <br>work on figuring out the cause for this discrepancy before calculating the acceptance. <br>More details can be seen in her slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/mhuang/07012015/07012015.pdf here].

*Pengjia
**Comparing asymmetries determined from the MAID model to those calculated from data.<br>Last time he showed a factor of 6 difference between the two results.  He determined the <br>scattering angle and Q<sup>2</sup> by fitting the data and compared them to the quantities calculated <br>from MAID.  He also showed the results of calculating the quantities g<sub>1</sub>, g<sub>2</sub>, F<sub>1</sub>, F<sub>2</sub>, A<sub>L</sub>, A<sub>T</sub>,<br> dXS<sub>L</sub>, dXS<sub>T</sub>, XS<sub>tot</sub> andXS<sub>mott</sub> using MAID.  He hasn't included radiative corrections in his <br>calculations, but this probably won't account for the factor of 6.  More details can be seen <br>in his slides [http://hallaweb.jlab.org/experiment/g2p/collaborators/pzhu/07012015/MAID_20150701.pdf here].