HCalDQM.h

#ifndef _DQM_HCAL_DQM_H_
#define _DQM_HCAL_DQM_H_
//----------//
//   STL    //
//----------//
#include <algorithm>
 
//----------//
//   ROOT   //
//----------//
#include "TVector3.h"
 
/*~~~~~~~~~~~~~~~*/
/*   Framework   */
/*~~~~~~~~~~~~~~~*/
#include <map>
 
#include "DetDescr/HcalGeometry.h"
#include "DetDescr/HcalID.h"
#include "Framework/Configure/Parameters.h"
#include "Framework/Event.h"
#include "Framework/EventFile.h"
#include "Framework/EventProcessor.h"
#include "Hcal/Event/HcalHit.h"
#include "Hcal/Event/HcalVetoResult.h"
#include "SimCore/Event/SimCalorimeterHit.h"
#include "Tools/AnalysisUtils.h"
namespace dqm {
 
class HCalDQM : public framework::Analyzer {
 public:
  HCalDQM(const std::string &name, framework::Process &process);
 
  ~HCalDQM() {}
 
  void configure(framework::config::Parameters &parameters) override;
 
  void analyze(const framework::Event &event) override;
 
  bool skipHit(const ldmx::HcalID &id) {
    const auto section{id.section()};
    return (section != section_ && section_ != -1);
  }
  void analyzeRecHits(const std::vector<ldmx::HcalHit> &hits);
  void analyzeSimHits(const std::vector<ldmx::SimCalorimeterHit> &hits);
 
  bool hitPassesVeto(const ldmx::HcalHit &hit, int section) {
    if (hit.getPE() < pe_veto_threshold || hit.getTime() > max_hit_time_) {
      return true;
    }
    if (section == ldmx::HcalID::HcalSection::BACK && hit.getMinPE() < 1) {
      return true;
    }
    return false;
  }
 
 private:
  std::string sim_coll_name_;
 
  std::string sim_pass_name_;
  std::string rec_coll_name_;
 
  std::string rec_pass_name_;
 
  std::string veto_name_;
 
  std::string veto_pass_;
 
  // Veto threshold for photo-electrons
  float pe_veto_threshold;
  int section_;
  double max_hit_time_;
};
 
}  // namespace dqm
 
#endif  // _DQM_HCAL_DQM_H_