CMSSW_4_4_3_patch1/src/SimCalorimetry/HcalSimAlgos/src/HcalShapes.cc

Go to the documentation of this file.

#include "SimCalorimetry/HcalSimAlgos/interface/HcalShapes.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "SimCalorimetry/CaloSimAlgos/interface/CaloCachedShapeIntegrator.h"
#include "CondFormats/HcalObjects/interface/HcalMCParam.h"
#include "CondFormats/HcalObjects/interface/HcalMCParams.h"
#include "CondFormats/DataRecord/interface/HcalMCParamsRcd.h"
#include "DataFormats/HcalDetId/interface/HcalGenericDetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

HcalShapes::HcalShapes()
: theMCParams(0),
  theShapes(5),
  theHcalShape(),
  theHFShape(),
  theZDCShape(),
  theSiPMShape()
{
  theShapes[0] = new CaloCachedShapeIntegrator(&theHcalShape);
  theShapes[1] = new CaloCachedShapeIntegrator(&theHcalShape);
  theShapes[2] = new CaloCachedShapeIntegrator(&theSiPMShape);
  theShapes[3] = new CaloCachedShapeIntegrator(&theHFShape);
  theShapes[4] = new CaloCachedShapeIntegrator(&theZDCShape);
}


HcalShapes::~HcalShapes()
{
  for(std::vector<const CaloVShape *>::const_iterator shapeItr = theShapes.begin();
      shapeItr != theShapes.end();  ++shapeItr)
  {
    delete *shapeItr;
  }
  theShapes.clear();
  delete theMCParams;
}


void HcalShapes::beginRun(edm::EventSetup const & es)
{
  edm::ESHandle<HcalMCParams> p;
  es.get<HcalMCParamsRcd>().get(p);
  theMCParams = new HcalMCParams(*p.product()); 
}


void HcalShapes::endRun()
{
  delete theMCParams;
  theMCParams = 0;
}


const CaloVShape * HcalShapes::shape(const DetId & detId) const
{
  if(!theMCParams) {
    edm::LogWarning("HcalShapes") << "Cannot find HCAL MC Params ";
    // try to figure the appropriate shape
    HcalGenericDetId::HcalGenericSubdetector subdet = HcalGenericDetId(detId).genericSubdet();
    if(subdet == HcalGenericDetId::HcalGenBarrel || subdet == HcalGenericDetId::HcalGenEndcap) return theShapes[0];
    else if(subdet == HcalGenericDetId::HcalGenOuter) return theShapes[2];
    else if(subdet == HcalGenericDetId::HcalGenForward) return theShapes[3];
    else if(subdet == HcalGenericDetId::HcalGenZDC) return theShapes[4];
    else return 0;
  }
  int shapeType = theMCParams->getValues(detId)->signalShape();
  return theShapes.at(shapeType);
}