EcalTBMCInfoProducer Class Reference

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=default
	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
const double	beamEta
const double	beamPhi
const double	beamTheta
const double	beamXoff
const double	beamYoff
int	crysNumber
const double	fMaxEta
const double	fMaxPhi
const double	fMinEta
const double	fMinPhi
std::unique_ptr< ROOT::Math::Rotation3D >	fromCMStoTB
const edm::EDGetTokenT< edm::HepMCProduct >	GenVtxToken
double	partXhodo
double	partYhodo
std::unique_ptr< 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 33 of file EcalTBMCInfoProducer.cc.

Constructor & Destructor Documentation

EcalTBMCInfoProducer()

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

explicit

Constructor.

Definition at line 67 of file EcalTBMCInfoProducer.cc.

References beamEta, beamPhi, beamTheta, beamXoff, beamYoff, funct::cos(), crysNumber, ecalTB2006H4_GenSimDigiReco_cfg::CrystalMapFile, spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, PVValHelper::eta, Exception, fromCMStoTB, edm::ParameterSet::getParameter(), EcalTBCrystalMap::NCRYSTAL, phi, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, theTestMap, geometryCSVtoXML::xx, geometryCSVtoXML::xy, geometryCSVtoXML::xz, geometryCSVtoXML::yy, geometryCSVtoXML::yz, and geometryCSVtoXML::zz.

    : fMinEta(ps.getParameter<double>("MinEta")),
      fMaxEta(ps.getParameter<double>("MaxEta")),
      fMinPhi(ps.getParameter<double>("MinPhi")),
      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")),
      GenVtxToken(consumes<edm::HepMCProduct>(edm::InputTag("moduleLabelVtx", "source"))) {
  produces<PEcalTBInfo>();

  edm::FileInPath CrystalMapFile = ps.getParameter<edm::FileInPath>("CrystalMapFile");

  std::string fullMapName = CrystalMapFile.fullPath();
  theTestMap = std::make_unique<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::LogVerbatim("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 = std::make_unique<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.";
  }
}

~EcalTBMCInfoProducer()

EcalTBMCInfoProducer::~EcalTBMCInfoProducer ( )

overridedefault

Destructor.

Member Function Documentation

produce()

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

override

Produce digis out of raw data.

Definition at line 145 of file EcalTBMCInfoProducer.cc.

References beamEta, beamPhi, beamXoff, beamYoff, crysNumber, GenVtxToken, edm::RandomNumberGenerator::getEngine(), edm::HepMCProduct::GetEvent(), LogDebug, eostools::move(), partXhodo, and partYhodo.

                                                                                   {
  edm::Service<edm::RandomNumberGenerator> rng;
  CLHEP::HepRandomEngine *engine = &rng->getEngine(event.streamID());

  std::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.;

  const edm::Handle<edm::HepMCProduct> &GenEvt = event.getHandle(GenVtxToken);

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

Member Data Documentation

beamEta

const double EcalTBMCInfoProducer::beamEta

private

Definition at line 49 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer(), and produce().

beamPhi

const double EcalTBMCInfoProducer::beamPhi

private

Definition at line 50 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer(), and produce().

beamTheta

const double EcalTBMCInfoProducer::beamTheta

private

Definition at line 51 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer().

beamXoff

const double EcalTBMCInfoProducer::beamXoff

private

Definition at line 52 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer(), and produce().

beamYoff

const double EcalTBMCInfoProducer::beamYoff

private

Definition at line 53 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer(), and produce().

crysNumber

int EcalTBMCInfoProducer::crysNumber

private

Definition at line 57 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer(), and produce().

fMaxEta

const double EcalTBMCInfoProducer::fMaxEta

private

Definition at line 46 of file EcalTBMCInfoProducer.cc.

fMaxPhi

const double EcalTBMCInfoProducer::fMaxPhi

private

Definition at line 48 of file EcalTBMCInfoProducer.cc.

fMinEta

const double EcalTBMCInfoProducer::fMinEta

private

Definition at line 45 of file EcalTBMCInfoProducer.cc.

fMinPhi

const double EcalTBMCInfoProducer::fMinPhi

private

Definition at line 47 of file EcalTBMCInfoProducer.cc.

fromCMStoTB

std::unique_ptr<ROOT::Math::Rotation3D> EcalTBMCInfoProducer::fromCMStoTB

private

Definition at line 64 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer().

GenVtxToken

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

private

Definition at line 55 of file EcalTBMCInfoProducer.cc.

Referenced by produce().

partXhodo

double EcalTBMCInfoProducer::partXhodo

private

Definition at line 59 of file EcalTBMCInfoProducer.cc.

Referenced by produce().

partYhodo

double EcalTBMCInfoProducer::partYhodo

private

Definition at line 60 of file EcalTBMCInfoProducer.cc.

Referenced by produce().

theTestMap

std::unique_ptr<EcalTBCrystalMap> EcalTBMCInfoProducer::theTestMap

private

Definition at line 62 of file EcalTBMCInfoProducer.cc.

Referenced by EcalTBMCInfoProducer().