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) : param_(conf), clusterizer_(conf) {
  src_ = conf.getParameter<std::string>("label");
  verbosity_ = conf.getUntrackedParameter<int>("RPixVerbosity");
 
  tokenCTPPSPixelDigi_ = consumes<edm::DetSetVector<CTPPSPixelDigi> >(edm::InputTag(src_));
 
  produces<edm::DetSetVector<CTPPSPixelCluster> >();
}
 
CTPPSPixelClusterProducer::~CTPPSPixelClusterProducer() {}
 
void CTPPSPixelClusterProducer::fillDescriptions(edm::ConfigurationDescriptions &descriptions) {
  edm::ParameterSetDescription desc;
  desc.addUntracked<int>("RPixVerbosity", 0);
  desc.add<std::string>("label", "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);
 
  // get analysis mask to mask channels
  edm::ESHandle<CTPPSPixelAnalysisMask> aMask;
 
  if (!rpd->empty())
    iSetup.get<CTPPSPixelAnalysisMaskRcd>().get(aMask);
 
  edm::DetSetVector<CTPPSPixelCluster> output;
 
  // run clusterisation
  if (!rpd->empty()) {
    // get calibration DB
    theGainCalibrationDB.getDB(iEvent, iSetup);
    run(*rpd, output, aMask.product());
  }
  // 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) {
  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, theGainCalibrationDB.getCalibs(), 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);