5#ifndef LArTBH6BeamInfo_H
6#define LArTBH6BeamInfo_H
10#include "Identifier/Identifier.h"
33 virtual StatusCode
execute()
override;
34 virtual StatusCode
finalize()
override;
37 typedef std::vector<double>
dVect;
44 double &a1,
double &a2,
double &
chi2,
dVect &residual);
47 Gaudi::Property<bool>
m_Primary{
this,
"PrimaryTrackOnly",
true};
48 Gaudi::Property<int>
m_pcode{
this,
"PrimaryParticle", 999};
55 std::vector< SG::ReadHandle< AthenaHitsVector<LArG4H6FrontHit> > >
m_hitcoll;
Handle class for reading from StoreGate.
Handle class for recording to StoreGate.
AthAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor with parameters:
Gaudi::Property< int > m_pcode
virtual StatusCode execute() override
virtual StatusCode finalize() override
Gaudi::Property< bool > m_Primary
LArTBH6BeamInfo(const std::string &name, ISvcLocator *pSvcLocator)
virtual ~LArTBH6BeamInfo()=default
bool fitVect(const dVect &vec_x, const dVect &vec_xz, const dVect &vec_ex, double &a1, double &a2, double &chi2, dVect &residual)
Fit data to the function u = a1 + a2*w.
std::vector< SG::ReadHandle< AthenaHitsVector< LArG4H6FrontHit > > > m_hitcoll
virtual StatusCode initialize() override
SG::WriteHandle< TBTrack > m_track
std::vector< double > dVect
SG::ReadHandle< TBEventInfo > m_theEventInfo
Gaudi::Property< std::vector< std::string > > m_HitsCollNames
double chi2(TH1 *h0, TH1 *h1)