FakeTBHodoscopeRawInfoProducer.cc

Go to the documentation of this file.

/*
 * \file FakeTBHodoscopeRawInfoProducer.cc
 *
 *
 */
 
#include "SimG4CMS/EcalTestBeam/interface/FakeTBHodoscopeRawInfoProducer.h"
 
using namespace cms;
using namespace std;
 
FakeTBHodoscopeRawInfoProducer::FakeTBHodoscopeRawInfoProducer(const edm::ParameterSet &ps) {
  ecalTBInfo_ = consumes<PEcalTBInfo>(edm::InputTag("EcalTBInfoLabel", "SimEcalTBG4Object"));
  produces<EcalTBHodoscopeRawInfo>();
 
  theTBHodoGeom_ = new EcalTBHodoscopeGeometry();
}
 
FakeTBHodoscopeRawInfoProducer::~FakeTBHodoscopeRawInfoProducer() { delete theTBHodoGeom_; }
 
void FakeTBHodoscopeRawInfoProducer::produce(edm::Event &event, const edm::EventSetup &eventSetup) {
  unique_ptr<EcalTBHodoscopeRawInfo> product(new EcalTBHodoscopeRawInfo());
 
  // get the vertex information from the event
 
  edm::Handle<PEcalTBInfo> theEcalTBInfo;
  event.getByToken(ecalTBInfo_, theEcalTBInfo);
 
  double partXhodo = theEcalTBInfo->evXbeam();
  double partYhodo = theEcalTBInfo->evYbeam();
 
  LogDebug("EcalTBHodo") << "TB frame vertex (X,Y) for hodoscope simulation: \n"
                         << "x = " << partXhodo << " y = " << partYhodo;
 
  // for each hodoscope plane determine the fibre number corresponding
  // to the event vertex coordinates in the TB reference frame
  // plane 0/2 = x plane 1/3 = y
 
  int nPlanes = (int)theTBHodoGeom_->getNPlanes();
  product->setPlanes(nPlanes);
 
  for (int iPlane = 0; iPlane < nPlanes; ++iPlane) {
    float theCoord = (float)partXhodo;
    if (iPlane == 1 || iPlane == 3)
      theCoord = (float)partYhodo;
 
    vector<int> firedChannels = theTBHodoGeom_->getFiredFibresInPlane(theCoord, iPlane);
    unsigned int nChannels = firedChannels.size();
 
    EcalTBHodoscopePlaneRawHits planeHit(nChannels);
    for (unsigned int i = 0; i < nChannels; ++i) {
      planeHit.addHit(firedChannels[i]);
    }
 
    product->setPlane((unsigned int)iPlane, planeHit);
  }
 
  LogDebug("EcalTBHodo") << (*product);
 
  event.put(std::move(product));
}