ClusterProducerFP420.h

Go to the documentation of this file.

#ifndef ClusterProducerFP420_h
#define ClusterProducerFP420_h

#include "RecoRomanPot/RecoFP420/interface/ClusterNoiseFP420.h"

#include "DataFormats/FP420Cluster/interface/ClusterFP420.h"
#include "DataFormats/FP420Digi/interface/HDigiFP420.h"

#include <vector>
#include <algorithm>
#include <cmath>

class ClusterProducerFP420 {
public:
  typedef std::vector<HDigiFP420>::const_iterator HDigiFP420Iter;

  ClusterProducerFP420(float electrode_thr, float seed_thr, float clust_thr, int max_voids)
      : theChannelThreshold(electrode_thr),
        theSeedThreshold(seed_thr),
        theClusterThreshold(clust_thr),
        max_voids_(max_voids){};

  std::vector<ClusterFP420> clusterizeDetUnit(HDigiFP420Iter begin,
                                              HDigiFP420Iter end,
                                              unsigned int detid,
                                              const ElectrodNoiseVector& vnoise) const;
  std::vector<ClusterFP420> clusterizeDetUnitPixels(HDigiFP420Iter begin,
                                                    HDigiFP420Iter end,
                                                    unsigned int detid,
                                                    const ElectrodNoiseVector& vnoise,
                                                    unsigned int xytype,
                                                    int verb) const;

  int difNarr(unsigned int xytype, HDigiFP420Iter ichannel, HDigiFP420Iter jchannel) const;
  int difWide(unsigned int xytype, HDigiFP420Iter ichannel, HDigiFP420Iter jchannel) const;

  float channelThresholdInNoiseSigma() const { return theChannelThreshold; }
  float seedThresholdInNoiseSigma() const { return theSeedThreshold; }
  float clusterThresholdInNoiseSigma() const { return theClusterThreshold; }

private:
  float theChannelThreshold;
  float theSeedThreshold;
  float theClusterThreshold;
  int max_voids_;

  bool badChannel(int channel, const std::vector<short>& badChannels) const;
};

class AboveSeed {
public:
  AboveSeed(float aseed, const ElectrodNoiseVector& vnoise) : verb(0), seed(aseed), vnoise_(vnoise){};

  bool operator()(const HDigiFP420& digi) const {
    return (!vnoise_[digi.channel()].getDisable() && digi.adc() >= seed * vnoise_[digi.channel()].getNoise());
  }

private:
  int verb;
  float seed;
  const ElectrodNoiseVector& vnoise_;
};

#endif