EcalTBMCInfoProducer Class Reference

#include <EcalTBMCInfoProducer.h>

Inheritance diagram for EcalTBMCInfoProducer:

Public Member Functions
	EcalTBMCInfoProducer (const edm::ParameterSet &ps)
	Constructor. More...
void	produce (edm::Event &event, const edm::EventSetup &eventSetup) override
	Produce digis out of raw data. More...
	~EcalTBMCInfoProducer () override
	Destructor. More...
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
double	beamEta
double	beamPhi
double	beamTheta
double	beamXoff
double	beamYoff
int	crysNumber
ROOT::Math::Rotation3D *	fromCMStoTB
edm::EDGetTokenT< edm::HepMCProduct >	GenVtxToken
double	partXhodo
double	partYhodo
EcalTBCrystalMap *	theTestMap

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 22 of file EcalTBMCInfoProducer.h.

Constructor & Destructor Documentation

EcalTBMCInfoProducer()

EcalTBMCInfoProducer::EcalTBMCInfoProducer ( const edm::ParameterSet &  ps )

explicit

Constructor.

Definition at line 22 of file EcalTBMCInfoProducer.cc.

                                                                    {
  produces<PEcalTBInfo>();
 
  edm::FileInPath CrystalMapFile = ps.getParameter<edm::FileInPath>("CrystalMapFile");
  GenVtxToken = consumes<edm::HepMCProduct>(edm::InputTag("moduleLabelVtx", "source"));
  double fMinEta = ps.getParameter<double>("MinEta");
  double fMaxEta = ps.getParameter<double>("MaxEta");
  double fMinPhi = ps.getParameter<double>("MinPhi");
  double fMaxPhi = ps.getParameter<double>("MaxPhi");
  beamEta = (fMaxEta + fMinEta) / 2.;
  beamPhi = (fMaxPhi + fMinPhi) / 2.;
  beamTheta = 2.0 * atan(exp(-beamEta));
  beamXoff = ps.getParameter<double>("BeamMeanX");
  beamYoff = ps.getParameter<double>("BeamMeanX");
 
  string fullMapName = CrystalMapFile.fullPath();
  theTestMap = new EcalTBCrystalMap(fullMapName);
  crysNumber = 0;
 
  double deltaEta = 999.;
  double deltaPhi = 999.;
  for (int cryIndex = 1; cryIndex <= EcalTBCrystalMap::NCRYSTAL; ++cryIndex) {
    double eta = 0;
    double phi = 0.;
    theTestMap->findCrystalAngles(cryIndex, eta, phi);
    if (fabs(beamEta - eta) < deltaEta && fabs(beamPhi - phi) < deltaPhi) {
      deltaEta = fabs(beamEta - eta);
      deltaPhi = fabs(beamPhi - phi);
      crysNumber = cryIndex;
    } else if (fabs(beamEta - eta) < deltaEta && fabs(beamPhi - phi) > deltaPhi) {
      if (fabs(beamPhi - phi) < 0.017) {
        deltaEta = fabs(beamEta - eta);
        deltaPhi = fabs(beamPhi - phi);
        crysNumber = cryIndex;
      }
    } else if (fabs(beamEta - eta) > deltaEta && fabs(beamPhi - phi) < deltaPhi) {
      if (fabs(beamEta - eta) < 0.017) {
        deltaEta = fabs(beamEta - eta);
        deltaPhi = fabs(beamPhi - phi);
        crysNumber = cryIndex;
      }
    }
  }
 
  edm::LogInfo("EcalTBInfo") << "Initialize TB MC ECAL info producer with parameters: \n"
                             << "Crystal map file:  " << CrystalMapFile << "\n"
                             << "Beam average eta = " << beamEta << "\n"
                             << "Beam average phi = " << beamPhi << "\n"
                             << "Corresponding to crystal number = " << crysNumber << "\n"
                             << "Beam X offset =    " << beamXoff << "\n"
                             << "Beam Y offset =    " << beamYoff;
 
  // rotation matrix to move from the CMS reference frame to the test beam one
 
  double xx = -cos(beamTheta) * cos(beamPhi);
  double xy = -cos(beamTheta) * sin(beamPhi);
  double xz = sin(beamTheta);
 
  double yx = sin(beamPhi);
  double yy = -cos(beamPhi);
  double yz = 0.;
 
  double zx = sin(beamTheta) * cos(beamPhi);
  double zy = sin(beamTheta) * sin(beamPhi);
  double zz = cos(beamTheta);
 
  fromCMStoTB = new ROOT::Math::Rotation3D(xx, xy, xz, yx, yy, yz, zx, zy, zz);
 
  // random number
  edm::Service<edm::RandomNumberGenerator> rng;
  if (!rng.isAvailable()) {
    throw cms::Exception("Configuration") << "EcalTBMCInfoProducer requires the RandomNumberGeneratorService\n"
                                             "which is not present in the configuration file.  You must add the "
                                             "service\n"
                                             "in the configuration file or remove the modules that require it.";
  }
}

References funct::cos(), ecalTB2006H4_GenSimDigiReco_cfg::CrystalMapFile, spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, PVValHelper::eta, Exception, JetChargeProducer_cfi::exp, edm::ParameterSet::getParameter(), edm::Service< T >::isAvailable(), EcalTBCrystalMap::NCRYSTAL, funct::sin(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, geometryCSVtoXML::xz, geometryCSVtoXML::yy, geometryCSVtoXML::yz, and geometryCSVtoXML::zz.

~EcalTBMCInfoProducer()

EcalTBMCInfoProducer::~EcalTBMCInfoProducer ( )

override

Destructor.

Definition at line 100 of file EcalTBMCInfoProducer.cc.

{ delete theTestMap; }

Member Function Documentation

produce()

void EcalTBMCInfoProducer::produce ( edm::Event &  event, const edm::EventSetup &  eventSetup  )

override

Produce digis out of raw data.

Definition at line 102 of file EcalTBMCInfoProducer.cc.

                                                                                   {
  edm::Service<edm::RandomNumberGenerator> rng;
  CLHEP::HepRandomEngine *engine = &rng->getEngine(event.streamID());
 
  unique_ptr<PEcalTBInfo> product(new PEcalTBInfo());
 
  // Fill the run information
 
  product->setCrystal(crysNumber);
 
  product->setBeamDirection(beamEta, beamPhi);
  product->setBeamOffset(beamXoff, beamYoff);
 
  // Compute the event x,y vertex coordinates in the beam reference system
  // e.g. in the place orthogonal to the beam average direction
 
  partXhodo = partYhodo = 0.;
 
  edm::Handle<edm::HepMCProduct> GenEvt;
  event.getByToken(GenVtxToken, GenEvt);
 
  const HepMC::GenEvent *Evt = GenEvt->GetEvent();
  HepMC::GenEvent::vertex_const_iterator Vtx = Evt->vertices_begin();
 
  math::XYZPoint eventCMSVertex((*Vtx)->position().x(), (*Vtx)->position().y(), (*Vtx)->position().z());
 
  LogDebug("EcalTBInfo") << "Generated vertex position = " << eventCMSVertex.x() << " " << eventCMSVertex.y() << " "
                         << eventCMSVertex.z();
 
  math::XYZPoint eventTBVertex = (*fromCMStoTB) * eventCMSVertex;
 
  LogDebug("EcalTBInfo") << "Rotated vertex position   = " << eventTBVertex.x() << " " << eventTBVertex.y() << " "
                         << eventTBVertex.z();
 
  partXhodo = eventTBVertex.x();
  partYhodo = eventTBVertex.y();
 
  product->setBeamPosition(partXhodo, partYhodo);
 
  // Asynchronous phase shift
  double thisPhaseShift = CLHEP::RandFlat::shoot(engine);
 
  product->setPhaseShift(thisPhaseShift);
  LogDebug("EcalTBInfo") << "Asynchronous Phaseshift = " << thisPhaseShift;
 
  // store the object in the framework event
 
  event.put(std::move(product));
}

References edm::RandomNumberGenerator::getEngine(), edm::HepMCProduct::GetEvent(), LogDebug, and eostools::move().

Member Data Documentation

beamEta

double EcalTBMCInfoProducer::beamEta

private

Definition at line 34 of file EcalTBMCInfoProducer.h.

beamPhi

double EcalTBMCInfoProducer::beamPhi

private

Definition at line 35 of file EcalTBMCInfoProducer.h.

beamTheta

double EcalTBMCInfoProducer::beamTheta

private

Definition at line 36 of file EcalTBMCInfoProducer.h.

beamXoff

double EcalTBMCInfoProducer::beamXoff

private

Definition at line 40 of file EcalTBMCInfoProducer.h.

beamYoff

double EcalTBMCInfoProducer::beamYoff

private

Definition at line 41 of file EcalTBMCInfoProducer.h.

crysNumber

int EcalTBMCInfoProducer::crysNumber

private

Definition at line 38 of file EcalTBMCInfoProducer.h.

fromCMStoTB

ROOT::Math::Rotation3D* EcalTBMCInfoProducer::fromCMStoTB

private

Definition at line 48 of file EcalTBMCInfoProducer.h.

GenVtxToken

edm::EDGetTokenT<edm::HepMCProduct> EcalTBMCInfoProducer::GenVtxToken

private

Definition at line 50 of file EcalTBMCInfoProducer.h.

partXhodo

double EcalTBMCInfoProducer::partXhodo

private

Definition at line 43 of file EcalTBMCInfoProducer.h.

partYhodo

double EcalTBMCInfoProducer::partYhodo

private

Definition at line 44 of file EcalTBMCInfoProducer.h.

theTestMap

EcalTBCrystalMap* EcalTBMCInfoProducer::theTestMap

private

Definition at line 46 of file EcalTBMCInfoProducer.h.