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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 23 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 35 of file EcalTBMCInfoProducer.h.

beamPhi

double EcalTBMCInfoProducer::beamPhi

private

Definition at line 36 of file EcalTBMCInfoProducer.h.

beamTheta

double EcalTBMCInfoProducer::beamTheta

private

Definition at line 37 of file EcalTBMCInfoProducer.h.

beamXoff

double EcalTBMCInfoProducer::beamXoff

private

Definition at line 41 of file EcalTBMCInfoProducer.h.

beamYoff

double EcalTBMCInfoProducer::beamYoff

private

Definition at line 42 of file EcalTBMCInfoProducer.h.

crysNumber

int EcalTBMCInfoProducer::crysNumber

private

Definition at line 39 of file EcalTBMCInfoProducer.h.

fromCMStoTB

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

private

Definition at line 49 of file EcalTBMCInfoProducer.h.

GenVtxToken

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

private

Definition at line 51 of file EcalTBMCInfoProducer.h.

partXhodo

double EcalTBMCInfoProducer::partXhodo

private

Definition at line 44 of file EcalTBMCInfoProducer.h.

partYhodo

double EcalTBMCInfoProducer::partYhodo

private

Definition at line 45 of file EcalTBMCInfoProducer.h.

theTestMap

EcalTBCrystalMap* EcalTBMCInfoProducer::theTestMap

private

Definition at line 47 of file EcalTBMCInfoProducer.h.