EcalCleaningAlgo.h

Go to the documentation of this file.

#ifndef __EcalCleaningAlgo_h_
#define __EcalCleaningAlgo_h_

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include <vector>

class DetId;

class EcalCleaningAlgo {
public:
  EcalCleaningAlgo(const edm::ParameterSet& p);

  EcalRecHit::Flags checkTopology(const DetId& id, const EcalRecHitCollection& rhs);

  void setFlags(EcalRecHitCollection& rhs);

private:
  float e4e1(const DetId& id, const EcalRecHitCollection& rhs);

  float e6e2(const DetId& id, const EcalRecHitCollection& rhs);

  float recHitE(const DetId id, const EcalRecHitCollection& recHits, bool useTimingInfo);

  bool isNearCrack(const DetId& detid);

  const std::vector<DetId> neighbours(const DetId& id);

  float ignoreOutOfTimeThresh_;

  // Parameters for tolopogical cut
  // mark anomalous if e> cThreshold &&  e4e1> a*log10(e1e1)+b
  float cThreshold_barrel_;
  float cThreshold_endcap_;
  float e4e1_a_barrel_;
  float e4e1_b_barrel_;
  float e4e1_a_endcap_;
  float e4e1_b_endcap_;
  // when calculating e4/e1, ignore hits below this threshold
  float e4e1Treshold_barrel_;
  float e4e1Treshold_endcap_;
  float tightenCrack_e1_single_;
  float tightenCrack_e4e1_single_;
  float cThreshold_double_;
  float tightenCrack_e1_double_;
  float tightenCrack_e6e2_double_;
  float e6e2thresh_;
};

#endif  // __EcalCleaningAlgo_h_

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: