GnumiTree.h

Go to the documentation of this file.

#ifndef GNUMI_H
#define GNUMI_H

class TTree;
class TH2D;
class GnumiTree {
    TTree* tree;

public:
    int run, evtno;
    float ndxdz, ndydz, npz, nenergy; // random decay neutrino
    float ndxdznea,ndydznea,nenergyn,nwtnear; // near detector neutrino
    float ndxdzfar,ndydzfar,nenergyf,nwtfar; //  far  detector neutrino
    int norig;                  // 1=primary hadron, 2=secondary interaction, 3=muon
    int ndecay;                 // decay type
    int ntype;                  // neutrino type
    float vx,vy,vz;             // Vertex of hadron decay (cm)
    float pdpx,pdpy,pdpz;       // Momentum of neutrino parent at decay point
    float ppdxdz,ppdydz,pppz,ppenergy; // direction and energy of parent at parent production
    int ppmedium, ptype;        // GEANT medium of parent production, parent type (GEANT pid)
    float ppvx, ppvy, ppvz;     // Parent production vertex
    float muparpx;              // If parent is muon
    float muparpy;              //    then vector of hadron parent of muon
    float muparpz;              // else if muon produced with neutrino  
    float mupare;               //    then muon vector
                                // else zero
    float necm;                 // E_nu in parent cm
    float nimpwt;               // Weight of event (from rejecting low P hadrons)
    float xpoint,ypoint,zpoint; // special run values selectable in gbeam.input
    float tvx,tvy,tvz;          // Target exit point or in target production point
    float tpx,tpy,tpz;          // Ibid but momentum;
    int tptype;                 // particle type exiting target
    
    void SetTree(TTree& tree);
    void paint_flux(TH2D* hist, double npot, double zloc=1039.8e2);
};


#endif  // GNUMI_H

---