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 <vector>
#  include <utility> // std::move()

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