HcalVetoProcessor.h

Go to the documentation of this file.

 
#ifndef __HCAL_HCAL_VETO_PROCESSOR_H__
#define __HCAL_HCAL_VETO_PROCESSOR_H__
 
//----------------//
//   C++ StdLib   //
//----------------//
#include <string>
 
//----------//
//   LDMX   //
//----------//
#include "Event/HcalHit.h"
#include "Framework/Configure/Parameters.h"
#include "Framework/EventProcessor.h"
#include "Tracking/Event/Track.h"
 
namespace hcal {
 
class HcalVetoProcessor : public framework::Producer {
 public:
  HcalVetoProcessor(const std::string& name, framework::Process& process);
 
  virtual ~HcalVetoProcessor() = default;
 
  void configure(framework::config::Parameters& parameters) override;
 
  void produce(framework::Event& event) override;
 
  std::vector<float> trackProp(const ldmx::Tracks& tracks,
                               ldmx::TrackStateType ts_type,
                               const std::string& ts_title);
 
 private:
  double total_pe_threshold_{5};
 
  float max_time_{50};  // ns
 
  float back_min_pe_{1};
 
  /*
   * A hit representing the case where we never reach the maxPE condition. This
   * is rare but can happen which previously would record uninitialized memory.
   * Stored inside of the producer so that it will have a valid lifetime when we
   * persist it.
   *
   * See https://github.com/LDMX-Software/Hcal/issues/58 for details
   *
   * Ideally, this would just be stored as a part of the HcalVetoResult, but
   * changing that would be breaking change so for now we work around it like
   * this.
   *
   * It contains nonsense values but since they are predictable, they are harder
   * to mistake for real hits_. See constructor for the actual values.
   */
  ldmx::HcalHit default_max_hit_;
 
  std::string output_coll_name_;
  std::string input_hit_coll_name_;
  std::string input_hit_pass_name_;
  std::string track_pass_name_;
  bool exclude_recoil_ele_;
  std::string track_collection_;
  float dr_from_recoil_max_;
  bool inverse_skim_{false};
};  // HcalVetoProcessor
}  // namespace hcal
 
#endif  // HCAL_HCALVETOPROCESSOR_H_