EcalTimeDigiProducer.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Event.h"
#include "SimCalorimetry/EcalSimProducers/interface/EcalTimeDigiProducer.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalTimeMapDigitizer.h"

#include "CalibFormats/CaloObjects/interface/CaloSamples.h"

#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "Geometry/EcalAlgo/interface/EcalEndcapGeometry.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/Common/interface/Handle.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "SimGeneral/MixingModule/interface/PileUpEventPrincipal.h"


//#define ecal_time_debug 1

EcalTimeDigiProducer::EcalTimeDigiProducer( const edm::ParameterSet& params, edm::ProducerBase& mixMod,edm::ConsumesCollector& sumes) :
   DigiAccumulatorMixMod(),
   m_EBdigiCollection ( params.getParameter<std::string>("EBtimeDigiCollection") ) ,
   m_EEdigiCollection ( params.getParameter<std::string>("EEtimeDigiCollection") ) ,
   m_hitsProducerTagEB  ( params.getParameter<edm::InputTag>("hitsProducerEB"    ) ) ,
   m_hitsProducerTagEE  ( params.getParameter<edm::InputTag>("hitsProducerEE"    ) ) ,
   m_hitsProducerTokenEB  ( sumes.consumes<std::vector<PCaloHit> >( m_hitsProducerTagEB) ) ,
   m_hitsProducerTokenEE  ( sumes.consumes<std::vector<PCaloHit> >( m_hitsProducerTagEE) ) ,
   m_timeLayerEB     (  params.getParameter<int> ("timeLayerBarrel") ),
   m_timeLayerEE     (  params.getParameter<int> ("timeLayerEndcap") ),
   m_Geometry          ( nullptr ) 
{
   mixMod.produces<EcalTimeDigiCollection>(m_EBdigiCollection);
   mixMod.produces<EcalTimeDigiCollection>(m_EEdigiCollection);

   m_BarrelDigitizer = new EcalTimeMapDigitizer(EcalBarrel);
   m_EndcapDigitizer = new EcalTimeMapDigitizer(EcalEndcap);

#ifdef ecal_time_debug
   std::cout << "[EcalTimeDigiProducer]::Create EB " << m_EBdigiCollection << " and EE " << m_EEdigiCollection << " collections and digitizers" << std::endl;
#endif

   m_BarrelDigitizer->setTimeLayerId(m_timeLayerEB);
   m_EndcapDigitizer->setTimeLayerId(m_timeLayerEE);
}

EcalTimeDigiProducer::~EcalTimeDigiProducer() 
{
}

void
EcalTimeDigiProducer::initializeEvent(edm::Event const& event, edm::EventSetup const& eventSetup) {
   checkGeometry( eventSetup );
   //    checkCalibrations( event, eventSetup );
   // here the methods to clean the maps
   m_BarrelDigitizer->initializeMap();
   m_EndcapDigitizer->initializeMap();
}

void
EcalTimeDigiProducer::accumulateCaloHits(HitsHandle const& ebHandle, HitsHandle const& eeHandle, int bunchCrossing) {
  // accumulate the simHits and do the averages in a given layer per bunch crossing
  if(ebHandle.isValid()) {
    m_BarrelDigitizer->add(*ebHandle.product(), bunchCrossing);
  }

  if(eeHandle.isValid()) {
    m_EndcapDigitizer->add(*eeHandle.product(), bunchCrossing);
  }

}

void
EcalTimeDigiProducer::accumulate(edm::Event const& e, edm::EventSetup const& eventSetup) {
  // Step A: Get Inputs
  edm::Handle<std::vector<PCaloHit> > ebHandle;
  e.getByToken(m_hitsProducerTokenEB, ebHandle);

  edm::Handle<std::vector<PCaloHit> > eeHandle;
  e.getByToken(m_hitsProducerTokenEE, eeHandle);

#ifdef ecal_time_debug
  std::cout << "[EcalTimeDigiProducer]::Accumulate Hits HS  event" << std::endl;
#endif

  accumulateCaloHits(ebHandle, eeHandle, 0);
}

void
EcalTimeDigiProducer::accumulate(PileUpEventPrincipal const& e, edm::EventSetup const& eventSetup, edm::StreamID const&) {
  edm::Handle<std::vector<PCaloHit> > ebHandle;
  e.getByLabel(m_hitsProducerTagEB, ebHandle);

  edm::Handle<std::vector<PCaloHit> > eeHandle;
  e.getByLabel(m_hitsProducerTagEE, eeHandle);

#ifdef ecal_time_debug
  std::cout << "[EcalTimeDigiProducer]::Accumulate Hits for BC " << e.bunchCrossing() << std::endl;
#endif
  accumulateCaloHits(ebHandle, eeHandle, e.bunchCrossing());
}

void 
EcalTimeDigiProducer::finalizeEvent(edm::Event& event, edm::EventSetup const& eventSetup) {
   std::unique_ptr<EcalTimeDigiCollection> barrelResult   ( new EcalTimeDigiCollection() ) ;
   std::unique_ptr<EcalTimeDigiCollection> endcapResult   ( new EcalTimeDigiCollection() ) ;

#ifdef ecal_time_debug
   std::cout << "[EcalTimeDigiProducer]::finalizeEvent" << std::endl;
#endif

   // here basically just put everything in the final collections
   m_BarrelDigitizer->run( *barrelResult ) ;

#ifdef ecal_time_debug
   std::cout << "[EcalTimeDigiProducer]::EB Digi size " <<  barrelResult->size() << std::endl;
#endif

   edm::LogInfo("TimeDigiInfo") << "EB time Digis: " << barrelResult->size() ;

   m_EndcapDigitizer->run( *endcapResult ) ;

#ifdef ecal_time_debug
   std::cout << "[EcalTimeDigiProducer]::EE Digi size " <<  endcapResult->size() << std::endl;
#endif

    edm::LogInfo("TimeDigiInfo") << "EE Digis: " << endcapResult->size() ;

#ifdef ecal_time_debug
    std::cout << "[EcalTimeDigiProducer]::putting collections into the event " << std::endl;
#endif

    event.put( std::move(barrelResult),    m_EBdigiCollection ) ;
    event.put( std::move(endcapResult),    m_EEdigiCollection ) ;
}

void 
EcalTimeDigiProducer::checkGeometry( const edm::EventSetup & eventSetup ) 
{
   // TODO find a way to avoid doing this every event
   edm::ESHandle<CaloGeometry>               hGeometry   ;
   eventSetup.get<CaloGeometryRecord>().get( hGeometry ) ;

   const CaloGeometry* pGeometry = &*hGeometry;

   if( pGeometry != m_Geometry )
   {
      m_Geometry = pGeometry;
      updateGeometry();
   }
}

void
EcalTimeDigiProducer::updateGeometry() 
{
   m_BarrelDigitizer->setGeometry(
      m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalBarrel    ) ) ;
   m_EndcapDigitizer->setGeometry(
      m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalEndcap    ) ) ;
}