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: