This is just to give a very brief update on how the simulations are going.

Some reminder on features:
  Ability to simulate multiple particles in events
  Creates simulated Qt output based eithier on energy deposition or optical photon tracking
  Can create QT data for entire event, or just the effect of individual particles
  Stores simulated information on particles that we choose for future comparison

New Features:
  Added event number based on macro file names
  propogated event number to final output file for comparision

With this implimented we can now start to cross reference event by event various quantities in reconstruction vs simulation. As an example, I've included a comparision between the total photonic energy in the event (including pi0-yy decays) vs E12 of the reconstruction.

Some notes on this:
1st, I'd like to start off by saying that this plot is more of an example of the kinds of things we can do, rather than any usable results at this stage. I haven't implimented any checks to see if pi0 decays both hit the detector, checked on if additional photons are within cluster or not, made cut on photon number, etc...

With that said here's a listing of various important quantities for the simulation
  ~1000 pythia events that had a pi0 with pt>=2 aimed at FMS(eta cut)
  keeping all final state particles (minus more exotic ones such as neutrinos) within same eta cut
  simulated QT based on photon tracking
  cluster size for reconstruction 200 miliradians
  ESim is the sum of all saved trajectories with 0 mass (right now any primary photons plus photons from primary nuetral meson decays)

Future work:
  1). Look into fast way to link entries of final output tree to simulation tree (attached pdf made by looping through and checking event numbers)
  2). Learn about/impliment triggering.