NtpMCTruth.h

Go to the documentation of this file.

// NtpMCTruth
//
// NtpMCTruth is an ntuple class to hold mc truth data.
//
// Based on Roy Lee's CandEventSR/MCTruthSR class
//

#ifndef NTPMCTRUTH_H
#define NTPMCTRUTH_H

#include "TObject.h"       // base class
#include "NtpMCFluxInfo.h" // data member
#include "NtpMCFluxWgt.h"  // data member

class NtpMCTruth;
std::ostream &operator << (std::ostream& os, const NtpMCTruth& mc);

class NtpMCTruth : public TObject {

 public:
  // Constructors & Destructors
  NtpMCTruth();
  virtual ~NtpMCTruth();

  // State testing methods
  virtual std::ostream& Print(std::ostream& os) const;
  virtual void Print(const Option_t* option = "") const;

 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.

  UShort_t index; // index of this mctruth entry in mc array
  Int_t stdhep[2]; // index range in stdhep array of associated particles
  Int_t inu;      // PDG id of neutrino, nu_e=12, nu_mu=14, nu_tau=16, anti=-
  Int_t inunoosc; // PDG id of un-oscillated neutrino
  Int_t itg;      // PDG id of target nucleon, p=2212, n=2112
  Int_t iboson;   // PDG id of exchange boson, Z0=23, W+=24
  // resonance code: 1001=QE, 1002=resonance prod., 1003=deep inelastic
  // 1004=coherent pion production, 1005=inverse muon decay
  Int_t iresonance; // resonance code
  Int_t iaction;    // 0=NC 1=CC (2=QE 3=QE+pi)
  Int_t istruckq;   // PDG id of struck quark
  // for Daikon (w/ neugen3 v3_5_0) iflags can be interpreted as:
  //  0: non-DIS,  1: old KNO hadronization, 2: modified KNO, 3: charm
  // 11: JETSET/string frag, 12: JETSET/cluster frag, 13: JETSET/other frag
  Int_t iflags;     // generator flags
  Int_t ndigu;      // total number of digits in u-view
  Int_t ndigv;      // total number of digits in v-view
  Float_t tphu;     // summed pulse height u-view
  Float_t tphv;     // summed pulse height v-view
  Float_t a;        // nucleus A
  Float_t z;        // nucleus Z
  Float_t sigma;    // total cross section for this interaction
  Float_t sigmadiff; // differential cross section
  Float_t x;        // Bj x
  Float_t y;        // y
  Float_t q2;       // Q**2
  Float_t w2;       // W**2
  Float_t emfrac;   // fraction of "shower" energy which is EM
  Float_t vtxx;     // vertex x-pos in meters (was cm previous to 8/31/04)     
  Float_t vtxy;     // vertex y-pos in meters (was cm previous to 8/31/04)
  Float_t vtxz;     // vertex z-pos in meters (was cm previous to 8/31/04)
  Float_t p4neu[4]; // neutrino momentum+energy (GeV) 4-vector
  Float_t p4neunoosc[4]; // unoscillated neutrino 4-vector
  // p4tgt Note: Fermi motion may mean that target is not at rest in lab frame 
  Float_t p4tgt[4];     // target nucleon 4-vector
  Float_t p4shw[4];     // shower (hadronic+em) 
  Float_t p4mu1[4];     // primary muon
  Float_t p4mu2[4];     // secondary muon (e.g. from charm production & decay)
  Float_t p4el1[4];     // first electron
  Float_t p4el2[4];     // highest E second electron
  Float_t p4tau[4];     // highest E tau (may have decayed)

  NtpMCFluxInfo flux;     // flux info associated with primary
  NtpMCFluxWgt fluxwgt;   // flux weight associated with primary

  ClassDef(NtpMCTruth,7)
};

#endif // NTPMCTRUTH_H

---