MCShower.h

Go to the documentation of this file.

#ifndef MCSHOWER_H
#define MCSHOWER_H

// LArSoft libraries
#include "lardataobj/MCBase/MCStep.h"

// Supporting tools
#include "nusimdata/SimulationBase/MCTruth.h"

// ROOT libraries
#include <TVector3.h>

// Standard C/C++ libraries
#include <string>
#include <utility> // std::move()
#include <vector>

namespace sim {

  class MCShower {

  public:
    MCShower() { ClearData(); }

    virtual ~MCShower();

    virtual void Clear() { ClearData(); }

    //--- Getters ---//

    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; }

    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 fAncestorProcess; }
    const MCStep& AncestorStart() const { return fAncestorStart; }
    const MCStep& AncestorEnd() const { return fAncestorEnd; }

    const MCStep& DetProfile() const { return fDetProfile; }

    const std::vector<unsigned int>& DaughterTrackID() const { return fDaughterTrackID; }

    double Charge(size_t plane) const;
    double dQdx(size_t plane) const;

    const std::vector<double>& Charge() const { return fPlaneCharge; }
    const std::vector<double>& dQdx() const { return fdQdx; }

    double dEdx() const { return fdEdx; }
    const TVector3& StartDir() const { return fStartDir; }

    //--- Setters ---//
    void Origin(simb::Origin_t o) { fOrigin = o; }

    void PdgCode(int id) { fPDGCode = id; }
    void TrackID(unsigned int id) { fTrackID = id; }
    void Process(const std::string& name) { fProcess = name; }
    void Start(const MCStep& s) { fStart = s; }
    void End(const MCStep& s) { fEnd = s; }
    void StartDir(const TVector3& sdir) { fStartDir = sdir; }

    void MotherPdgCode(int id) { fMotherPDGCode = id; }
    void MotherTrackID(unsigned int id) { fMotherTrackID = id; }
    void MotherProcess(const std::string& name) { fMotherProcess = name; }
    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& name) { fAncestorProcess = name; }
    void AncestorStart(const MCStep& s) { fAncestorStart = s; }
    void AncestorEnd(const MCStep& s) { fAncestorEnd = s; }

    void DetProfile(const MCStep& s) { fDetProfile = s; }

    void DaughterTrackID(const std::vector<unsigned int>& id_v) { fDaughterTrackID = id_v; }

    void Charge(const std::vector<double>& q) { fPlaneCharge = q; }
    void Charge(std::vector<double>&& q) { fPlaneCharge = std::move(q); }

    void dQdx(const std::vector<double>& dqdx) { fdQdx = dqdx; }
    void dQdx(std::vector<double>&& dqdx) { fdQdx = std::move(dqdx); }

    void dEdx(double dedx) { fdEdx = dedx; }
    void dEdxRAD(double dedx) { fdEdx_radial = dedx; }

  protected:
    simb::Origin_t fOrigin; 

    int fPDGCode;          
    unsigned int fTrackID; 
    std::string fProcess;  
    MCStep fStart;         
    MCStep fEnd;           
    TVector3 fStartDir;    

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

    int fAncestorPDGCode;          
    unsigned int fAncestorTrackID; 
    std::string fAncestorProcess;  
    MCStep fAncestorStart;         
    MCStep fAncestorEnd;           

    std::vector<unsigned int> fDaughterTrackID; 
    MCStep fDetProfile;                         
    double fdEdx;                               
    double fdEdx_radial;                        

    std::vector<double> fPlaneCharge; 
    std::vector<double> fdQdx;        

    void ClearData();

  }; // class MCShower

} // namespace sim

#endif    // MCSHOWER_H
 // end of doxygen group