Vertexer.h

#ifndef TRACKING_RECO_VERTEXER_H_
#define TRACKING_RECO_VERTEXER_H_
 
//--- Framework ---//
#include "Framework/Configure/Parameters.h"
#include "Framework/EventProcessor.h"
 
// --- Tracking --- //
#include "Tracking/Event/Track.h"
#include "Tracking/Sim/BFieldXYZUtils.h"
#include "Tracking/Sim/TrackingUtils.h"
 
// --- ACTS --- //
 
// Utils and definitions
 
#include "Acts/Definitions/Common.hpp"
#include "Acts/Definitions/Units.hpp"
 
// Vertexing
 
#include "Acts/Vertexing/FullBilloirVertexFitter.hpp"
#include "Acts/Vertexing/HelicalTrackLinearizer.hpp"
#include "Acts/Vertexing/Vertex.hpp"
 
// Magfield
 
#include "Acts/MagneticField/ConstantBField.hpp"
#include "Acts/MagneticField/MagneticFieldContext.hpp"
#include "Acts/MagneticField/MagneticFieldProvider.hpp"
 
// Propagator
 
#include "Acts/Propagator/EigenStepper.hpp"
#include "Acts/Propagator/Propagator.hpp"
//#include "Acts/Propagator/Navigator.hpp"
//#include "Acts/Propagator/StandardAborters.hpp"
 
// Geometry
#include "Acts/Surfaces/PerigeeSurface.hpp"
 
// Propagator with void navigator
using VoidPropagator = Acts::Propagator<Acts::EigenStepper<>>;
 
namespace tracking {
namespace reco {
 
class Vertexer : public framework::Producer {
 public:
  Vertexer(const std::string& name, framework::Process& process);
 
  virtual ~Vertexer() = default;
 
  void onProcessStart() override;
  void onProcessEnd() override;
 
  void configure(framework::config::Parameters& parameters) override;
 
  void produce(framework::Event& event) override;
 
  void TaggerRecoilMonitoring(const std::vector<ldmx::Track>& tagger_tracks,
                              const std::vector<ldmx::Track>& recoil_tracks);
 
 private:
  Acts::GeometryContext gctx_;
  Acts::MagneticFieldContext bctx_;
 
  int nevents_{0};
  int nvertices_{0};
  int nreconstructable_{0};
  std::shared_ptr<InterpolatedMagneticField3> sp_interpolated_bField_;
  std::shared_ptr<Acts::ConstantBField> bField_;
 
  std::string field_map_;
 
  std::string trk_c_name_1{"TaggerTracks"};
  std::string trk_c_name_2{"RecoilTracks"};
  std::shared_ptr<VoidPropagator> propagator_;
 
  // Monitoring histograms
  TH1F* h_delta_d0;
  TH1F* h_delta_z0;
  TH1F* h_delta_p;
 
  TH1F* h_delta_phi;
  TH1F* h_delta_theta;
 
  TH2F* h_delta_d0_vs_recoil_p;
  TH2F* h_delta_z0_vs_recoil_p;
 
  TH2F* h_td0_vs_rd0;
  TH2F* h_tz0_vs_rz0;
};
 
}  // namespace reco
}  // namespace tracking
 
#endif