ParticleFlowDQM.h

#ifndef DQM_PARTICLEFLOWDQM_H
#define DQM_PARTICLEFLOWDQM_H
 
#include <string>
#include <vector>
 
#include "Framework/Configure/Parameters.h"
#include "Framework/EventProcessor.h"
#include "Recon/Event/PFCandidate.h"
 
namespace dqm {
 
class ParticleFlowDQM : public framework::Analyzer {
 public:
  ParticleFlowDQM(const std::string& name, framework::Process& process)
      : Analyzer(name, process) {}
  ~ParticleFlowDQM() override = default;
 
  void configure(framework::config::Parameters& ps) override;
  void analyze(const framework::Event& event) override;
 
 private:
  void analyzeTruth(const framework::Event& event,
                    std::vector<ldmx::PFCandidate>& cands);
  bool reconstructable(int pdg, double px, double py, double pz,
                       double ke) const;
  int interactionRegion(const std::vector<double>& end) const;
 
  std::string pf_cand_coll_name_{"PFCandidates"};
  std::string pf_cand_pass_name_{""};
  std::string ecal_sp_coll_name_{"EcalScoringPlaneHits"};
  std::string ecal_sp_pass_name_{""};
  std::string sim_particles_coll_name_{"SimParticles"};
  std::string sim_particles_pass_name_{""};
  // minimum visible energy of a truth particle [MeV]
  double truth_min_energy_{50.};
  // neutral candidate to truth match window at the Ecal [mm]
  double neutral_match_dist_{30.};
  // max energy ratio either way for a neutral match
  double neutral_match_max_ratio_{2.};
  // Hcal only candidate to truth match cone from the target [deg]
  double neutral_hcal_match_angle_{30.};
  // reconstructable truth: max angle to the beam axis [deg], min |p| for
  // charged pi and for p, K [MeV], min KE for neutral hadrons and pi0 [MeV]
  double reco_max_theta_{80.};
  double reco_min_p_pion_{100.};
  double reco_min_p_heavy_{800.};
  double reco_min_ke_neutral_{1000.};
  double reco_min_ke_pi0_{2000.};
  // regions for neutral hadron end points [mm]: Ecal box up to ecal_back_z
  // within ecal_half_xy, side Hcal around it, back Hcal behind it
  double ecal_back_z_{760.};
  double ecal_half_xy_{300.};
  double hcal_max_z_{4000.};
  double hcal_max_xy_{1600.};
};
 
}  // namespace dqm
 
#endif