MCTrack.h

Go to the documentation of this file.

#ifndef MCTRACK_H
#define MCTRACK_H

#include <vector>
#include "MCStep.h"
#include "nusimdata/SimulationBase/MCTruth.h"

namespace sim{

  class MCTrack : public std::vector<sim::MCStep> {
    
  public:
    
    MCTrack() : std::vector<sim::MCStep>() {Clear();}
    
    virtual ~MCTrack(){}

    void Clear();


    simb::Origin_t      Origin   () const { return fOrigin;   }
    int                 PdgCode  () const { return fPDGCode;  }
    unsigned int        TrackID  () const { return fTrackID;  }
    const std::string&  Process  () const { return fProcess;  }
    const MCStep&       Start    () const { return fStart;    }
    const MCStep&       End      () const { return fEnd;      }
    const std::vector< std::vector<double> >& dQdx() const{ return fdQdx;}  // dQdx[# of MCSteps][# of plane]
    const std::vector<double>& dEdx() const {return fdEdx;} // dEdx[# of MCSteps]

    int                MotherPdgCode () const { return fMotherPDGCode; }
    unsigned int       MotherTrackID () const { return fMotherTrackID; }
    const std::string& MotherProcess () const { return fMotherProcess; }
    const MCStep&      MotherStart   () const { return fMotherStart;   }
    const MCStep&      MotherEnd     () const { return fMotherEnd;     }

    int                AncestorPdgCode () const { return fAncestorPDGCode; }
    unsigned int       AncestorTrackID () const { return fAncestorTrackID; }
    const std::string& AncestorProcess () const { return fMotherProcess;   }
    const MCStep&      AncestorStart   () const { return fAncestorStart;   }
    const MCStep&      AncestorEnd     () const { return fAncestorEnd;     }

    void Origin        ( simb::Origin_t o ) { fOrigin  = o;    }
    void PdgCode       ( int id           ) { fPDGCode = id;   }
    void TrackID       ( unsigned int id  ) { fTrackID = id;   }
    void Process       ( std::string name ) { fProcess = name; }
    void Start         ( const MCStep s   ) { fStart   = s;    }
    void End           ( const MCStep s   ) { fEnd     = s;    }
    void dEdx          ( const std::vector<double> s) { fdEdx = s;}
    void dQdx          ( const std::vector<std::vector<double> > s) { fdQdx = s;}


    void MotherPdgCode ( int id               ) { fMotherPDGCode = id; }
    void MotherTrackID ( unsigned int id      ) { fMotherTrackID = id; }
    void MotherProcess ( const std::string& n ) { fMotherProcess = n;  }
    void MotherStart   ( const MCStep& s      ) { fMotherStart   = s;  }
    void MotherEnd     ( const MCStep& s      ) { fMotherEnd     = s;  }

    void AncestorPdgCode ( int id               ) { fAncestorPDGCode = id; }
    void AncestorTrackID ( unsigned int id      ) { fAncestorTrackID = id; }
    void AncestorProcess ( const std::string& n ) { fAncestorProcess = n;  }
    void AncestorStart   ( const MCStep& s      ) { fAncestorStart   = s;  }
    void AncestorEnd     ( const MCStep& s      ) { fAncestorEnd     = s;  }


  protected:

    simb::Origin_t fOrigin;  
    int            fPDGCode; 
    unsigned int   fTrackID; 
    std::string    fProcess; 
    MCStep         fStart;   
    MCStep         fEnd;     
    std::vector<std::vector<double> > fdQdx; //< the G4 electron yeild per plane between each step // [N Plane][MCSteps - 1]
    std::vector<double> fdEdx;//< the G4 "ionization" energy loss between each step // [MCSteps - 1]

    int            fMotherPDGCode; 
    unsigned int   fMotherTrackID; 
    std::string    fMotherProcess; 
    MCStep         fMotherStart;   
    MCStep         fMotherEnd;     

    int            fAncestorPDGCode; 
    unsigned int   fAncestorTrackID; 
    std::string    fAncestorProcess; 
    MCStep         fAncestorStart;   
    MCStep         fAncestorEnd;     
  };
}
#endif
 // end of doxygen group