OverburdenCalculator.h

Go to the documentation of this file.

#ifndef OVERBURDENCALCULATOR__H
#define OVERBURDENCALCULATOR__H

#include <TVector3.h>
#include "art/Framework/Services/Registry/ServiceMacros.h"
#include "art/Framework/Services/Registry/ActivityRegistry.h"
#include "fhiclcpp/ParameterSet.h"
#include "RecoBase/Track.h"
#include "TGeoManager.h"

namespace red{
class OverburdenCalculator{
    public:
        OverburdenCalculator(const fhicl::ParameterSet& p, art::ActivityRegistry& reg);
    
        double trace (double x, double y, double z, double px, double py, double pz) const;
        double trace (const TVector3 start, const TVector3 dir) const {
            return trace(start.X(),start.Y(),start.Z(),dir.X(),dir.Y(),dir.Z());
        }
        double traceBack (const rb::Track& t)const {
            return trace(t.Stop(),-t.Dir());
        }
        void onBeginRun(const art::Run& r);
    protected:
        TGeoManager* fGeoManager = nullptr;
        std::string fExcludeVolume = "";
};

}
DECLARE_ART_SERVICE(red::OverburdenCalculator, LEGACY)
#endif