HcalNoiseAlgo.h

Go to the documentation of this file.

#ifndef _RECOMET_METALGORITHMS_HCALNOISEALGO_H_
#define _RECOMET_METALGORITHMS_HCALNOISEALGO_H_

#include "DataFormats/METReco/interface/HcalNoiseRBX.h"
#include "DataFormats/Common/interface/RefVector.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/Common/interface/RefProd.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

class CommonHcalNoiseRBXData {
  
 public:
  CommonHcalNoiseRBXData(const reco::HcalNoiseRBX& rbx, double minRecHitE, double minLowHitE, double minHighHitE,
      double TS4TS5EnergyThreshold, std::vector<std::pair<double, double> > &TS4TS5UpperCut,
      std::vector<std::pair<double, double> > &TS4TS5LowerCut, double MinRBXRechitR45E);
  ~CommonHcalNoiseRBXData() {}

  // accessors to internal variables
  inline double energy(void) const { return energy_; }
  inline double ratio(void) const { return e2ts_/e10ts_; }
  inline double e2ts(void) const { return e2ts_; }
  inline double e10ts(void) const { return e10ts_; }
  inline bool validRatio(void) const { return e10ts_!=0.0; }
  inline int numHPDHits(void) const { return numHPDHits_; }
  inline int numRBXHits(void) const { return numRBXHits_; }
  inline int numHPDNoOtherHits(void) const { return numHPDNoOtherHits_; }
  inline int numZeros(void) const { return numZeros_; }
  inline double minLowEHitTime(void) const { return minLowEHitTime_; }
  inline double maxLowEHitTime(void) const { return maxLowEHitTime_; }
  inline double lowEHitTimeSqrd(void) const { return lowEHitTimeSqrd_; }
  inline int numLowEHits(void) const { return numLowEHits_; }
  inline double minHighEHitTime(void) const { return minHighEHitTime_; }
  inline double maxHighEHitTime(void) const { return maxHighEHitTime_; }
  inline double highEHitTimeSqrd(void) const { return highEHitTimeSqrd_; }
  inline int numHighEHits(void) const { return numHighEHits_; }
  inline double RBXEMF(void) const { return RBXEMF_; }
  inline double HPDEMF(void) const { return HPDEMF_; }
  inline bool PassTS4TS5(void) const { return TS4TS5Decision_; }
  inline edm::RefVector<CaloTowerCollection> rbxTowers(void) const { return rbxtowers_; }
  inline int r45Count(void) const { return r45Count_; }
  inline double r45Fraction(void) const { return r45Fraction_; }
  inline double r45EnergyFraction(void) const { return r45EnergyFraction_; }

  bool CheckPassFilter(double Charge, double Discriminant, std::vector<std::pair<double, double> > &Cuts, int Side);

private:

  // values
  double energy_;            // RBX hadronic energy as determined by the sum of calotowers
  double e2ts_;              // pedestal subtracted charge in two peak TS for RBX
  double e10ts_;             // pedestal subtracted charge in all 10 TS for RBX
  int numHPDHits_;           // largest number of hits in an HPD in the RBX
  int numRBXHits_;           // number of hits in the RBX
  int numHPDNoOtherHits_;    // largest number of hits in an HPD when no other HPD has a hit in the RBX
  int numZeros_;             // number of ADC 0 counts in all hits in all TS in the RBX
  double minLowEHitTime_;    // minimum time found for any low energy hit in the RBX
  double maxLowEHitTime_;    // maximum time found for any low energy hit in the RBX
  double lowEHitTimeSqrd_;   // low energy hit time^2
  int numLowEHits_;          // number of low energy hits
  double minHighEHitTime_;   // minimum time found for any high energy hit in the RBX
  double maxHighEHitTime_;   // maximum time found for any high energy hit in the RBX
  double highEHitTimeSqrd_;  // high energy hit time^2
  int numHighEHits_;         // number of high energy hits
  double HPDEMF_;            // minimum electromagnetic fraction found in an HPD in the RBX
  double RBXEMF_;            // electromagnetic fraction of the RBX
  bool TS4TS5Decision_;      // if this RBX fails TS4TS5 variable or not
  edm::RefVector<CaloTowerCollection> rbxtowers_; // calotowers associated with the RBX
  int r45Count_;             // Number of rechits above some threshold flagged by R45
  double r45Fraction_;       // Fraction of rechits above some threshold flagged by R45
  double r45EnergyFraction_; // Energy fraction of rechits above some threshold

};

class HcalNoiseAlgo {

public:

  HcalNoiseAlgo(const edm::ParameterSet& iConfig);
  virtual ~HcalNoiseAlgo() {}

  // an rbx is "interesting/problematic" (i.e. is recorded to the event record)
  bool isProblematic(const CommonHcalNoiseRBXData&) const;

  // an rbx passes a noise filter
  bool passLooseNoiseFilter(const CommonHcalNoiseRBXData&) const;
  bool passTightNoiseFilter(const CommonHcalNoiseRBXData&) const;
  bool passHighLevelNoiseFilter(const CommonHcalNoiseRBXData&) const;

  // loose filter broken down into separate components
  bool passLooseRatio(const CommonHcalNoiseRBXData&) const;
  bool passLooseHits(const CommonHcalNoiseRBXData&) const;
  bool passLooseZeros(const CommonHcalNoiseRBXData&) const;
  bool passLooseTiming(const CommonHcalNoiseRBXData&) const;
  bool passLooseRBXRechitR45(const CommonHcalNoiseRBXData &) const;

  // tight filter broken down into separate components
  bool passTightRatio(const CommonHcalNoiseRBXData&) const;
  bool passTightHits(const CommonHcalNoiseRBXData&) const;
  bool passTightZeros(const CommonHcalNoiseRBXData&) const;
  bool passTightTiming(const CommonHcalNoiseRBXData&) const;
  bool passTightRBXRechitR45(const CommonHcalNoiseRBXData &) const;

  // an rbx passes an energy (or other) threshold to test a certain variable
  // for instance, the EMF cut might require that the RBX have 20 GeV of energy
  bool passRatioThreshold(const CommonHcalNoiseRBXData&) const;
  bool passZerosThreshold(const CommonHcalNoiseRBXData&) const;
  bool passEMFThreshold(const CommonHcalNoiseRBXData&) const;

  
private:

  // energy thresholds used for problematic cuts
  double pMinERatio_;         // minimum energy to apply ratio cuts
  double pMinEZeros_;         // minimum energy to apply zeros cuts
  double pMinEEMF_;           // minimum energy to apply EMF cuts

  // energy thresholds used for loose, tight and high level cuts
  double minERatio_;          // minimum energy to apply ratio cuts
  double minEZeros_;          // minimum energy to apply zeros cuts
  double minEEMF_;            // minimum energy to apply EMF cuts

  // "problematic" cuts
  // used to determine whether an RBX is stored in the EDM
  double pMinE_;              // minimum energy
  double pMinRatio_;          // minimum ratio
  double pMaxRatio_;          // maximum ratio
  int    pMinHPDHits_;        // minimum # of HPD hits
  int    pMinRBXHits_;        // minimum # of RBX hits
  int    pMinHPDNoOtherHits_; // minimum # of HPD hits with no other hits in the RBX
  int    pMinZeros_;          // minimum # of zeros
  double pMinLowEHitTime_;    // minimum low energy hit time
  double pMaxLowEHitTime_;    // maximum low energy hit time
  double pMinHighEHitTime_;   // minimum high energy hit time
  double pMaxHighEHitTime_;   // maximum high energy hit time
  double pMaxHPDEMF_;         // maximum HPD EMF
  double pMaxRBXEMF_;         // maximum RBX EMF
  int    pMinRBXRechitR45Count_;       // number of R45-flagged hits
  double pMinRBXRechitR45Fraction_;    // fraction of R45-flagged hits
  double pMinRBXRechitR45EnergyFraction_;   // energy fraction of R45-flagged hits

  // "loose" cuts
  // used to determine whether an RBX fails the loose noise cuts
  double lMinRatio_;          // minimum ratio
  double lMaxRatio_;          // maximum ratio
  int    lMinHPDHits_;        // minimum # of HPD hits
  int    lMinRBXHits_;        // minimum # of RBX hits
  int    lMinHPDNoOtherHits_; // minimum # of HPD hits with no other hits in the RBX
  int    lMinZeros_;          // minimum # of zeros
  double lMinLowEHitTime_;    // minimum low energy hit time
  double lMaxLowEHitTime_;    // maximum low energy hit time
  double lMinHighEHitTime_;   // minimum high energy hit time
  double lMaxHighEHitTime_;   // maximum high energy hit time
  std::vector<double> lMinRBXRechitR45Cuts_;

  // "tight" cuts
  // used to determine whether an RBX fails the tight noise cuts
  double tMinRatio_;          // minimum ratio
  double tMaxRatio_;          // maximum ratio
  int    tMinHPDHits_;        // minimum # of HPD hits
  int    tMinRBXHits_;        // minimum # of RBX hits
  int    tMinHPDNoOtherHits_; // minimum # of HPD hits with no other hits in the RBX
  int    tMinZeros_;          // minimum # of zeros
  double tMinLowEHitTime_;    // minimum low energy hit time
  double tMaxLowEHitTime_;    // maximum low energy hit time
  double tMinHighEHitTime_;   // minimum high energy hit time
  double tMaxHighEHitTime_;   // maximum high energy hit time
  std::vector<double> tMinRBXRechitR45Cuts_;

  // "high level" cuts
  // used to determine where an RBX fails the high level noise cuts
  double hlMaxHPDEMF_;        // maximum HPD EMF
  double hlMaxRBXEMF_;        // maximum RBX EMF
  
};

class JoinCaloTowerRefVectorsWithoutDuplicates
{
 public:
  JoinCaloTowerRefVectorsWithoutDuplicates() {}
  ~JoinCaloTowerRefVectorsWithoutDuplicates() {}

  void operator()(edm::RefVector<CaloTowerCollection>& v1, const edm::RefVector<CaloTowerCollection>& v2) const;

 private:

  // helper function to compare calotower references
  struct twrrefcomp {
    inline bool operator() ( const edm::Ref<CaloTowerCollection> & t1, const edm::Ref<CaloTowerCollection> & t2 ) const{
      return t1->id() < t2->id();
    }
  };
  typedef std::set< edm::Ref<CaloTowerCollection>, twrrefcomp> twrrefset_t;
};


#endif