MCTrack.h

Go to the documentation of this file.

#ifndef MCTRACK_H
#define MCTRACK_H

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

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