EcalMipTrackingFeatures.cxx

 
#include "DQM/EcalMipTrackingFeatures.h"
 
#include <cmath>
 
#include "Ecal/Event/EcalMipResult.h"
#include "Ecal/Event/EcalVetoResult.h"
 
namespace dqm {
 
void EcalMipTrackingFeatures::configure(framework::config::Parameters& ps) {
  ecal_mip_name_ = ps.get<std::string>("ecal_mip_name");
  ecal_mip_pass_ = ps.get<std::string>("ecal_mip_pass");
  // Maybe create two objects for the second recoil information
  ecal_veto_name_ = ps.get<std::string>("ecal_veto_name");
  ecal_veto_pass_ = ps.get<std::string>("ecal_veto_pass");
 
  return;
}
 
void EcalMipTrackingFeatures::analyze(const framework::Event& event) {
  const auto& mip_result{
      event.getObject<ldmx::EcalMipResult>(ecal_mip_name_, ecal_mip_pass_)};
  const auto& veto{
      event.getObject<ldmx::EcalVetoResult>(ecal_veto_name_, ecal_veto_pass_)};
 
  histograms_.fill("n_straight_tracks", mip_result.getNStraightTracks());
  histograms_.fill("n_linreg_segments", mip_result.getNLinRegTracks());
  histograms_.fill("first_near_photon_layer", mip_result.getFirstNearPhLayer());
  histograms_.fill("n_tracking_hits", veto.getNHitsInPhotonTerritory());
  histograms_.fill("ep_ang", veto.getEPAng());
  histograms_.fill("ep_sep", veto.getEPSep());
  auto recoil_mom = veto.getRecoilMomentum();
  histograms_.fill("recoil_pz", recoil_mom[2]);
  histograms_.fill("recoil_pt", std::sqrt(recoil_mom[0] * recoil_mom[0] +
                                          recoil_mom[1] * recoil_mom[1]));
  histograms_.fill("recoil_x", veto.getRecoilX());
  histograms_.fill("recoil_y", veto.getRecoilY());
 
  return;
}
 
}  // namespace dqm
 
DECLARE_ANALYZER(dqm::EcalMipTrackingFeatures);