SiStripClusterInfo.h

Go to the documentation of this file.

#ifndef SISTRIPCLUSTERIZER_SISTRIPCLUSTERINFO_H
#define SISTRIPCLUSTERIZER_SISTRIPCLUSTERINFO_H

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
#include <algorithm>
#include <numeric>

class SiStripNoises;
class SiStripGain;
class SiStripQuality;

class SiStripClusterInfo {
public:
  SiStripClusterInfo(const SiStripCluster& cluster,
                     const edm::EventSetup& es,
                     const int detid,
                     const std::string& qualityLabel = "");

  const SiStripCluster* cluster() const { return cluster_ptr; }

  uint32_t detId() const { return detId_; }
  uint16_t width() const { return cluster()->amplitudes().size(); }
  uint16_t firstStrip() const { return cluster()->firstStrip(); }
  float baryStrip() const { return cluster()->barycenter(); }
  uint16_t maxStrip() const { return firstStrip() + maxIndex(); }
  float variance() const;

  auto stripCharges() const -> decltype(cluster()->amplitudes()) { return cluster()->amplitudes(); }
  std::vector<float> stripGains() const;
  std::vector<float> stripNoises() const;
  std::vector<float> stripNoisesRescaledByGain() const;
  std::vector<bool> stripQualitiesBad() const;

  uint16_t charge() const { return std::accumulate(stripCharges().begin(), stripCharges().end(), uint16_t(0)); }
  uint8_t maxCharge() const { return *std::max_element(stripCharges().begin(), stripCharges().end()); }
  uint16_t maxIndex() const {
    return std::max_element(stripCharges().begin(), stripCharges().end()) - stripCharges().begin();
  }
  std::pair<uint16_t, uint16_t> chargeLR() const;

  float noise() const { return calculate_noise(stripNoises()); }
  float noiseRescaledByGain() const { return calculate_noise(stripNoisesRescaledByGain()); }

  float signalOverNoise() const { return charge() / noiseRescaledByGain(); }

  bool IsAnythingBad() const;
  bool IsApvBad() const;
  bool IsFiberBad() const;
  bool IsModuleBad() const;
  bool IsModuleUsable() const;

  std::vector<SiStripCluster> reclusterize(const edm::ParameterSet&) const;

private:
  float calculate_noise(const std::vector<float>&) const;

  const SiStripCluster* cluster_ptr;
  const edm::EventSetup& es;
  edm::ESHandle<SiStripNoises> noiseHandle;
  edm::ESHandle<SiStripGain> gainHandle;
  edm::ESHandle<SiStripQuality> qualityHandle;
  std::string qualityLabel;
  uint32_t detId_;
};

#endif