phys_tree/Event.h

Go to the documentation of this file.

// $Id$
//
//    File: Event.h
// Created: Wed Sep  2 20:18:06 EDT 2009
// Creator: davidl (on Darwin Amelia.local 9.8.0 i386)
//

#ifndef _Event_
#define _Event_

#include <TObject.h>
#include <TLorentzVector.h>
#include <TClonesArray.h>

// Maximum number of particles of a given type
#define MAX_PART 10

class Event:public TObject{

   public:

      Event();
      ~Event();
      void Clear(void);

      int event;
      UInt_t Npip;
      UInt_t Npim;
      UInt_t NKp;
      UInt_t NKm;
      UInt_t Nproton;
      UInt_t Nphoton;
      UInt_t Nneutron;

      TClonesArray *pip;
      TClonesArray *pim;
      TClonesArray *Kp;
      TClonesArray *Km;
      TClonesArray *proton;
      TClonesArray *photon;
      TClonesArray *neutron;

      TClonesArray *pip_match;      ///< Closest match pi+ (truth or recon, whatever is not in pip)
      TClonesArray *pim_match;      ///< Closest match pi- (truth or recon, whatever is not in pim)
      TClonesArray *Kp_match;       ///< Closest match K+ (truth or recon, whatever is not in Kp)
      TClonesArray *Km_match;       ///< Closest match K- (truth or recon, whatever is not in Km)
      TClonesArray *proton_match;   ///< Closest match proton (truth or recon, whatever is not in proton)
      TClonesArray *photon_match;   ///< Closest match photon (truth or recon, whatever is not in photon)
      TClonesArray *neutron_match;  ///< Closest match neutron (truth or recon, whatever is not in neutron)

      TLorentzVector target;  // Initial state target momentum
      TLorentzVector beam;    // Initial state target beam photon momentum
      TVector3 vertex;        // Vertex position
      TLorentzVector W;       // Final state 4-momentum of everything *except* the proton(s)

      bool all_fiducial;         // True if the following four are true
      bool all_mesons_fiducial;  // True if all pi+ and pi- are in fiducial region
      bool all_photons_fiducial; // True if all photons are in fiducial region
      bool all_neutrons_fiducial;   // True if all neutrons are in fiducial region
      bool all_protons_fiducial; // True if all protons are in fiducial region

   private:
      ClassDef(Event,1);

};

#endif // _Event_