HcalPulseContainmentManager.cc

Go to the documentation of this file.

#include "CalibCalorimetry/HcalAlgos/interface/HcalPulseContainmentManager.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "CondFormats/DataRecord/interface/HcalTimeSlewRecord.h"
#include <iostream>

HcalPulseContainmentManager::HcalPulseContainmentManager(float max_fracerror ) 
: entries_(),
  shapes_(),
  max_fracerror_(max_fracerror)
{
  hcalTimeSlew_delay_ = nullptr;
}

void HcalPulseContainmentManager::beginRun(edm::EventSetup const & es)
{
  edm::ESHandle<HcalTimeSlew> delay;
  es.get<HcalTimeSlewRecord>().get("HBHE", delay);
  hcalTimeSlew_delay_ = &*delay;

  shapes_.beginRun(es);
}

void HcalPulseContainmentManager::beginRun(const HcalDbService* conditions, const HcalTimeSlew* delay)
{
  hcalTimeSlew_delay_ = delay;

  shapes_.beginRun(conditions);
}

double HcalPulseContainmentManager::correction(const HcalDetId & detId, 
                                               int toAdd, float fixedphase_ns, double fc_ampl)
{

  /*
          int sub     = detId.subdet();
          int depth   = detId.depth();
          int inteta  = detId.ieta();
          int intphi  = detId.iphi();
          
          std::cout << "* HcalPulseContainmentManager::correction,   cell:" 
                    << " sub, ieta, iphi, depth = " 
                    << sub << "  " << inteta << "  " << intphi 
                    << "  " << depth  << "  toAaa= " << toAdd 
            <<"  phase = " << fixedphase_ns << "  ampl = " 
            << fc_ampl 
                    << std::endl;
  */

  return get(detId, toAdd, fixedphase_ns)->getCorrection(fc_ampl);
}

const HcalPulseContainmentCorrection * 
HcalPulseContainmentManager::get(const HcalDetId & detId, int toAdd, float fixedphase_ns)
{
  // const HcalPulseShape * shape = &(shapes_.shape(detId));
  const HcalPulseShape * shape = &(shapes_.shapeForReco(detId));
  for(std::vector<HcalPulseContainmentEntry>::const_iterator entryItr = entries_.begin();
      entryItr != entries_.end(); ++entryItr)
  {
    if (entryItr->shape_ == shape && entryItr->toAdd_ == toAdd && entryItr->fixedphase_ns_ == fixedphase_ns)
    {
        return &entryItr->correction_;
    }
  }

  /*
          int sub     = detId.subdet();
          int depth   = detId.depth();
          int inteta  = detId.ieta();
          int intphi  = detId.iphi();
          
          std::cout << "* HcalPulseContainmentManager::get new entry,  cell:" 
                    << " sub, ieta, iphi, depth = " 
                    << sub << "  " << inteta << "  " << intphi 
                    << "  " << depth   
                    << std::endl;
  */

  // didn't find it.  Make one.
  HcalPulseContainmentEntry entry(toAdd, fixedphase_ns, shape, HcalPulseContainmentCorrection(shape, toAdd, fixedphase_ns, max_fracerror_, hcalTimeSlew_delay_));
  entries_.push_back(entry);
  return &(entries_.back().correction_);
}