Updating ZDC SMD Data Analysis for 4 Spin States

  • Hal has calculated the asymmetries as a function of phi and per slat for the 4 spin states.
  • Code is now updated to reflect the 4 spin states.
  • Tested with 2 spin states to make sure code runs and reproduces old 2 state results.  This is also a check of my uncertainty calculation for the new, more complicated asymmetries.
  • Now the fit to the phi distribution includes 2 parameters: an amplitude and phase. 
  • Also I have changed the zdcAna code, which produces the histograms of good hits, to run on the MuDsts to save the step & disk space from making zdc trees.

Important note: You must comment the lines:

for(int iew=0; iew<2; iew++){
for(int ispin=2; ispin<4; ispin++){
for(int iphi=0; iphi<nPhiBins*2; iphi++){
countsPhi[iew][ispin][iphi]=1.0;
}
}
}

 when you run with 4 spin states to get the right answer!