NtpMCStdHep.h

Go to the documentation of this file.

// NtpMCStdHep
//
// NtpMCStdHep is an ntuple class to hold stdhep data for a single
// particle which can be either a primary or a secondary generated
// during particle transport through the detector. 
//
// Based on Roy Lee's CandEventSR/StdHepSR class
//

#ifndef NTPMCSTDHEP_H
#define NTPMCSTDHEP_H

#include "TMath.h"
#include "TObject.h"
#include "MCNtuple/NtpMCStdHepHit.h" // data member

class NtpMCStdHep;
std::ostream &operator << (std::ostream& os, const NtpMCStdHep& sh);

class NtpMCStdHep : public TObject {

 public:
  // Constructors & Destructors
  NtpMCStdHep();
  virtual ~NtpMCStdHep() {}

  // State testing methods
  virtual std::ostream& Print(std::ostream& os,
                              const Option_t* option="",
                              std::string printpreface="",
                              const TClonesArray* stdheparray=0) const;
  virtual void Print(const Option_t* option="") const;

  // particle momentum
  Float_t P() const {return TMath::Sqrt(p4[0]*p4[0]+p4[1]*p4[1]+p4[2]*p4[2]);}
  Float_t Charge() const; // particle charge
  std::string ParticleName() const; // particle name
       
 public:

  // Ntuple is treated like a C-struct with public data members and
  // rule-breaking field data members not prefaced by "f" and all
  // lowercase, by popular demand.
 
  UInt_t index;               // index of this stdhep entry in stdhep array
  Short_t mc;                 // index of associated mc primary in mc array
  Int_t parent[2];            // indices of first/last parent
  Int_t child[2];             // indices of first/last child 
  Int_t IstHEP;               // status code(see list in UserManual MC chapter)
  Int_t IdHEP;                // particle ID (PDG standard)
  Float_t mass;               // particle mass (GeV)
  Float_t p4[4];              // momentum 4-vector (GeV)
  Float_t vtx[4];             // vertex 4-vector ((x,y,z)(m),t(sec))
  Int_t ndethit;              // number of det digihits produced by this track
  NtpMCStdHepHit dethit[2];   // record info about first/last detector digihit
 
  ClassDef(NtpMCStdHep,8)
};
  
#endif // NTPMCSTDHEP_H

---