OverburdenCalculator.h
Go to the documentation of this file.
1 #ifndef OVERBURDENCALCULATOR__H
2 #define OVERBURDENCALCULATOR__H
3 
4 #include <TVector3.h>
8 #include "RecoBase/Track.h"
9 #include "TGeoManager.h"
10 
11 namespace red{
13  public:
15 
16  double trace (double x, double y, double z, double px, double py, double pz) const;
17  double trace (const TVector3 start, const TVector3 dir) const {
18  return trace(start.X(),start.Y(),start.Z(),dir.X(),dir.Y(),dir.Z());
19  }
20  double traceBack (const rb::Track& t)const {
21  return trace(t.Stop(),-t.Dir());
22  }
23  void onBeginRun(const art::Run& r);
24  protected:
25  TGeoManager* fGeoManager = nullptr;
27 };
28 
29 }
31 #endif
double traceBack(const rb::Track &t) const
const char * p
Definition: xmltok.h:285
#define DECLARE_ART_SERVICE(svc, scope)
Definition: ServiceMacros.h:91
A rb::Prong with full reconstructed trajectory.
Definition: Track.h:20
Definition: Run.h:31
virtual TVector3 Dir() const
Unit vector describing prong direction.
Definition: Prong.h:77
double trace(const TVector3 start, const TVector3 dir) const
z
Definition: test.py:28
OverburdenCalculator(const fhicl::ParameterSet &p, art::ActivityRegistry &reg)
TDirectory * dir
Definition: macro.C:5
TRandom3 r(0)
TVector3 Stop() const
Position of the final trajectory point.
Definition: Track.cxx:186
double trace(double x, double y, double z, double px, double py, double pz) const
enum BeamMode string