CastorDigiAnalyzer.cc

Go to the documentation of this file.

#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "SimCalorimetry/CastorSim/plugins/CastorDigiAnalyzer.h"
#include <iostream>

CastorDigiAnalyzer::CastorDigiAnalyzer(edm::ParameterSet const &conf)
    : hitReadoutName_("CastorHits"),
      simParameterMap_(),
      castorHitAnalyzer_("CASTORDigi", 1., &simParameterMap_, &castorFilter_),
      castorDigiStatistics_("CASTORDigi", 3, 10., 6., 0.1, 0.5, castorHitAnalyzer_),
      castorDigiCollectionTag_(conf.getParameter<edm::InputTag>("castorDigiCollectionTag")) {}

namespace CastorDigiAnalyzerImpl {
  template <class Collection>
  void analyze(edm::Event const &e, CastorDigiStatistics &statistics, edm::InputTag &tag) {
    edm::Handle<Collection> digis;
    e.getByLabel(tag, digis);
    if (!digis.isValid()) {
      edm::LogError("CastorDigiAnalyzer") << "Could not find Castor Digi Container ";
    } else {
      for (unsigned i = 0; i < digis->size(); ++i) {
        statistics.analyze((*digis)[i]);
      }
    }
  }
}  // namespace CastorDigiAnalyzerImpl

void CastorDigiAnalyzer::analyze(edm::Event const &e, edm::EventSetup const &c) {
  //  edm::Handle<edm::PCaloHitContainer> hits;
  edm::Handle<CrossingFrame<PCaloHit>> castorcf;

  e.getByLabel("mix", "g4SimHitsCastorFI", castorcf);

  // access to SimHits
  std::unique_ptr<MixCollection<PCaloHit>> hits(new MixCollection<PCaloHit>(castorcf.product()));
  //  if (hits.isValid()) {
  castorHitAnalyzer_.fillHits(*hits);
  CastorDigiAnalyzerImpl::analyze<CastorDigiCollection>(e, castorDigiStatistics_, castorDigiCollectionTag_);
}