DigiCollectionProfiler.h

Go to the documentation of this file.

#ifndef DPGAnalysis_SiStripTools_DigiCollectionProfile_H
#define DPGAnalysis_SiStripTools_DigiCollectionProfile_H

#include <vector>
#include "CommonTools/UtilAlgos/interface/DetIdSelector.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/SiStripDigi/interface/SiStripDigi.h"
#include "DataFormats/SiStripDigi/interface/SiStripRawDigi.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/DetSetVectorNew.h"
#include "DataFormats/Common/interface/DetSet.h"
#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"

#include "TH1.h"
#include "TH2.h"
#include "TProfile.h"

template <class T>
class DigiCollectionProfiler {
public:
  DigiCollectionProfiler();
  DigiCollectionProfiler(const edm::ParameterSet& iConfig);
  ~DigiCollectionProfiler(){};

  void fill(edm::Handle<T> digis,
            const std::vector<TH1F*>&,
            const std::vector<TProfile*>&,
            const std::vector<TH2F*>&) const;

private:
  bool m_folded;
  bool m_want1dHisto;
  bool m_wantProfile;
  bool m_want2dHisto;

  std::vector<DetIdSelector> m_selections;
};

template <class T>
DigiCollectionProfiler<T>::DigiCollectionProfiler()
    : m_folded(false), m_want1dHisto(false), m_wantProfile(false), m_want2dHisto(false), m_selections() {}

template <class T>
DigiCollectionProfiler<T>::DigiCollectionProfiler(const edm::ParameterSet& iConfig)
    : m_folded(iConfig.getUntrackedParameter<bool>("foldedStrips", false)),
      m_want1dHisto(iConfig.getUntrackedParameter<bool>("want1dHisto", true)),
      m_wantProfile(iConfig.getUntrackedParameter<bool>("wantProfile", true)),
      m_want2dHisto(iConfig.getUntrackedParameter<bool>("want2dHisto", false))

{
  std::vector<edm::ParameterSet> selconfigs = iConfig.getParameter<std::vector<edm::ParameterSet> >("selections");

  for (std::vector<edm::ParameterSet>::const_iterator selconfig = selconfigs.begin(); selconfig != selconfigs.end();
       ++selconfig) {
    DetIdSelector selection(*selconfig);
    m_selections.push_back(selection);
  }
}

template <class T>
void DigiCollectionProfiler<T>::fill(edm::Handle<T> digis,
                                     const std::vector<TH1F*>& hist,
                                     const std::vector<TProfile*>& hprof,
                                     const std::vector<TH2F*>& hist2d) const {}

template <>
inline void DigiCollectionProfiler<edm::DetSetVector<SiStripDigi> >::fill(
    edm::Handle<edm::DetSetVector<SiStripDigi> > digis,
    const std::vector<TH1F*>& hist,
    const std::vector<TProfile*>& hprof,
    const std::vector<TH2F*>& hist2d) const {
  for (edm::DetSetVector<SiStripDigi>::const_iterator mod = digis->begin(); mod != digis->end(); mod++) {
    for (unsigned int isel = 0; isel < m_selections.size(); ++isel) {
      if (m_selections[isel].isSelected(mod->detId())) {
        TH1F* tobefilled1d = nullptr;
        TProfile* tobefilledprof = nullptr;
        TH2F* tobefilled2d = nullptr;

        if (m_want1dHisto)
          tobefilled1d = hist[isel];
        if (m_wantProfile)
          tobefilledprof = hprof[isel];
        if (m_want2dHisto)
          tobefilled2d = hist2d[isel];

        for (edm::DetSet<SiStripDigi>::const_iterator digi = mod->begin(); digi != mod->end(); digi++) {
          if (digi->adc() > 0) {
            unsigned int strip = digi->strip();
            if (m_folded)
              strip = strip % 256;
            if (tobefilled1d)
              tobefilled1d->Fill(strip);
            if (tobefilledprof)
              tobefilledprof->Fill(strip, digi->adc());
            if (tobefilled2d)
              tobefilled2d->Fill(strip, digi->adc());
          }
        }
      }
    }
  }
}

template <>
inline void DigiCollectionProfiler<edm::DetSetVector<SiStripRawDigi> >::fill(
    edm::Handle<edm::DetSetVector<SiStripRawDigi> > digis,
    const std::vector<TH1F*>& hist,
    const std::vector<TProfile*>& hprof,
    const std::vector<TH2F*>& hist2d) const {
  for (edm::DetSetVector<SiStripRawDigi>::const_iterator mod = digis->begin(); mod != digis->end(); mod++) {
    for (unsigned int isel = 0; isel < m_selections.size(); ++isel) {
      if (m_selections[isel].isSelected(mod->detId())) {
        TH1F* tobefilled1d = nullptr;
        TProfile* tobefilledprof = nullptr;
        TH2F* tobefilled2d = nullptr;

        if (m_want1dHisto)
          tobefilled1d = hist[isel];
        if (m_wantProfile)
          tobefilledprof = hprof[isel];
        if (m_want2dHisto)
          tobefilled2d = hist2d[isel];

        unsigned int istrip = 0;
        for (edm::DetSet<SiStripRawDigi>::const_iterator digi = mod->begin(); digi != mod->end(); digi++, ++istrip) {
          if (digi->adc() > 0) {
            unsigned int strip = istrip;
            if (m_folded)
              strip = strip % 256;
            if (tobefilled1d)
              tobefilled1d->Fill(strip);
            if (tobefilledprof)
              tobefilledprof->Fill(strip, digi->adc());
            if (tobefilled2d)
              tobefilled2d->Fill(strip, digi->adc());
          }
        }
      }
    }
  }
}

template <>
inline void DigiCollectionProfiler<edmNew::DetSetVector<SiStripCluster> >::fill(
    edm::Handle<edmNew::DetSetVector<SiStripCluster> > digis,
    const std::vector<TH1F*>& hist,
    const std::vector<TProfile*>& hprof,
    const std::vector<TH2F*>& hist2d) const {
  for (edmNew::DetSetVector<SiStripCluster>::const_iterator mod = digis->begin(); mod != digis->end(); mod++) {
    for (unsigned int isel = 0; isel < m_selections.size(); ++isel) {
      if (m_selections[isel].isSelected(mod->detId())) {
        TH1F* tobefilled1d = nullptr;
        TProfile* tobefilledprof = nullptr;
        TH2F* tobefilled2d = nullptr;

        if (m_want1dHisto)
          tobefilled1d = hist[isel];
        if (m_wantProfile)
          tobefilledprof = hprof[isel];
        if (m_want2dHisto)
          tobefilled2d = hist2d[isel];

        for (edmNew::DetSet<SiStripCluster>::const_iterator clus = mod->begin(); clus != mod->end(); clus++) {
          for (unsigned int digi = 0; digi < clus->amplitudes().size(); ++digi) {
            if (clus->amplitudes()[digi] > 0) {
              unsigned int strip = clus->firstStrip() + digi;
              if (m_folded)
                strip = strip % 256;
              if (tobefilled1d)
                tobefilled1d->Fill(strip);
              if (tobefilledprof)
                tobefilledprof->Fill(strip, clus->amplitudes()[digi]);
              if (tobefilled2d)
                tobefilled2d->Fill(strip, clus->amplitudes()[digi]);
            }
          }
        }
      }
    }
  }
}

#endif  // DPGAnalysis_SiStripTools_DigiCollectionProfile_H