Reconstruction, July 21, 2009

From GlueXWiki
Jump to: navigation, search

Time: 11:00 EST/10:00 CST


Documents to Review

Tentative Agenda

  1. Physics Simulation Meetings
  2. Generating Background Events (Pythia: 124 μb @ 9GeV beam energy)
  3. Reconstruction: Physics Reactions

Discussion

  1. Blake looked at his bggen folders to find that a BGGEN hddm file is ~1kB/event (10^7 events is ~10GB) and an HDGEANT hddm file (from filtered BGGEN events) is ~27kB/event. Filtering the data for neutral energy greater than 6GeV culled 80% of events. The final data set is after runnning though HDGEANT ~55GB for 2x10^6 events therefore an unfiltered set of 10^7 events would have a data storage size of 275GB.
  2. Reconstruction


Reconstruction

  • It is required that there are 4 or 5(~5% of accepted events) photons showers (none are due to charged particles - comparing swum tracks to shower positions), and 1 proton. This gives 6(30) possible combinations.
  • Pi0 and Eta candidates are chosen as the eta-pi0 pair with the largest probability (from the TwoGammaFit )
  • Before χ^2 cuts ~50% of the generated signal events can be reconstructed. A cut on χ^2 < 3 leaves approximately 39% of generated events,
  • 10^7 Pythia hadronic background events were generated and then filtered by requiring that the total neutral energy be greater than 6GeV. 50,000 signal a2(1318) events were generated. This gives a ratio of cross-sections being 200:1 where the Pythia cross section at 9 GeV is 124 μb giving a cross section for the a2 to be ~0.6 μb. The 4 photon, 1 proton requirement leaves 12,690 background events and 24,870 signal events (5% more including 5 photon events).
  • Protons (acceptance and smearing) are handle via HDParSim, a parametric Monte Carlo used as a plugin with the analysis code



π0 mass: the lines in black are Pythia + signal events, red are cutting on χ^2 < 3, green is Pythia background only cut on χ^2 < 3.

PythiaA2 pi0 MASS.png PythiaA2 pi0 Prob.png


η mass: the lines in black are Pythia + signal events, red are cutting on χ^2 < 3,, green is Pythia background only cut on χ^2 < 3.

PythiaA2 eta MASS.png PythiaA2 eta Prob.png

resonance mass

  • a2(1318) resonance with angular distribution from IU generator.
  • Prob(a2) = Prob(π0)*Prob(η)

PythiaA2 res MASS.png PythiaA2 res Prob.png


Total transverse momentum (proton + 4γ)

PythiaA2 Pt.png

Gottfreid-Jackson angle vs resonance mass

π0 and η: χ^2 < 3

PythiaA2 res cosGJMASS.png

Minutes