CTPPSPixelClusterProducer.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "RecoPPS/Local/interface/CTPPSPixelClusterProducer.h"

CTPPSPixelClusterProducer::CTPPSPixelClusterProducer(const edm::ParameterSet &conf)
    : tokenCTPPSPixelDigi_(consumes<edm::DetSetVector<CTPPSPixelDigi> >(conf.getParameter<edm::InputTag>("tag"))),
      tokenCTPPSPixelAnalysisMask_(esConsumes()),
      tokenGainCalib_(esConsumes()),
      verbosity_(conf.getUntrackedParameter<int>("RPixVerbosity")),
      clusterizer_(conf) {
  produces<edm::DetSetVector<CTPPSPixelCluster> >();
}

CTPPSPixelClusterProducer::~CTPPSPixelClusterProducer() {}

void CTPPSPixelClusterProducer::fillDescriptions(edm::ConfigurationDescriptions &descriptions) {
  edm::ParameterSetDescription desc;
  desc.addUntracked<int>("RPixVerbosity", 0);
  desc.add<edm::InputTag>("tag", edm::InputTag("ctppsPixelDigis"));
  desc.add<int>("SeedADCThreshold", 2);
  desc.add<int>("ADCThreshold", 2);
  desc.add<double>("ElectronADCGain", 135.0);
  desc.add<int>("VCaltoElectronGain", 50);
  desc.add<int>("VCaltoElectronOffset", -411);
  desc.add<bool>("doSingleCalibration", false);
  descriptions.add("ctppsPixelClusters", desc);
}

void CTPPSPixelClusterProducer::produce(edm::Event &iEvent, const edm::EventSetup &iSetup) {
  edm::Handle<edm::DetSetVector<CTPPSPixelDigi> > rpd;
  iEvent.getByToken(tokenCTPPSPixelDigi_, rpd);

  edm::DetSetVector<CTPPSPixelCluster> output;

  if (!rpd->empty()) {
    // get analysis mask to mask channels
    const auto &mask = iSetup.getData(tokenCTPPSPixelAnalysisMask_);
    // get calibration DB
    const auto &gainCalibrations = iSetup.getData(tokenGainCalib_);
    // run clusterisation
    run(*rpd, output, mask, gainCalibrations);
  }
  // write output
  iEvent.put(std::make_unique<edm::DetSetVector<CTPPSPixelCluster> >(output));
}

void CTPPSPixelClusterProducer::run(const edm::DetSetVector<CTPPSPixelDigi> &input,
                                    edm::DetSetVector<CTPPSPixelCluster> &output,
                                    const CTPPSPixelAnalysisMask &mask,
                                    const CTPPSPixelGainCalibrations &gainCalibration) {
  for (const auto &ds_digi : input) {
    edm::DetSet<CTPPSPixelCluster> &ds_cluster = output.find_or_insert(ds_digi.id);
    clusterizer_.buildClusters(ds_digi.id, ds_digi.data, ds_cluster.data, &gainCalibration, &mask);

    if (verbosity_) {
      unsigned int cluN = 0;
      for (std::vector<CTPPSPixelCluster>::iterator iit = ds_cluster.data.begin(); iit != ds_cluster.data.end();
           iit++) {
        edm::LogInfo("CTPPSPixelClusterProducer") << "Cluster " << ++cluN << " avg row " << (*iit).avg_row()
                                                  << " avg col " << (*iit).avg_col() << " ADC.size " << (*iit).size();
      }
    }
  }
}

DEFINE_FWK_MODULE(CTPPSPixelClusterProducer);