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));
}