CTPPSLocalTrackLiteProducer Class Reference

Distills the essential track data from all RPs. More...

Inheritance diagram for CTPPSLocalTrackLiteProducer:

Public Member Functions
	CTPPSLocalTrackLiteProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &)

Private Attributes
edm::EDGetTokenT< edm::DetSetVector< CTPPSDiamondLocalTrack > >	diamondTrackToken_
bool	includeDiamonds_
bool	includePixels_
bool	includeStrips_
edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelLocalTrack > >	pixelTrackToken_
double	pixelTrackTxMax_
double	pixelTrackTxMin_
double	pixelTrackTyMax_
double	pixelTrackTyMin_
edm::EDGetTokenT< edm::DetSetVector< TotemRPLocalTrack > >	siStripTrackToken_
double	timingTrackTMax_
double	timingTrackTMin_

Static Private Attributes
static float	HPTDC_TIME_SLICE_WIDTH = 25.
	HPTDC time slice width, in ns. More...

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

Distills the essential track data from all RPs.

Definition at line 31 of file CTPPSLocalTrackLiteProducer.cc.

Constructor & Destructor Documentation

CTPPSLocalTrackLiteProducer::CTPPSLocalTrackLiteProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 57 of file CTPPSLocalTrackLiteProducer.cc.

References diamondTrackToken_, edm::ParameterSet::getParameter(), pixelTrackToken_, siStripTrackToken_, and ctppsLocalTrackLiteProducer_cff::tagPixelTrack.

    : includeStrips_(iConfig.getParameter<bool>("includeStrips")),
      includeDiamonds_(iConfig.getParameter<bool>("includeDiamonds")),
      includePixels_(iConfig.getParameter<bool>("includePixels")),
      pixelTrackTxMin_(iConfig.getParameter<double>("pixelTrackTxMin")),
      pixelTrackTxMax_(iConfig.getParameter<double>("pixelTrackTxMax")),
      pixelTrackTyMin_(iConfig.getParameter<double>("pixelTrackTyMin")),
      pixelTrackTyMax_(iConfig.getParameter<double>("pixelTrackTyMax")),
      timingTrackTMin_(iConfig.getParameter<double>("timingTrackTMin")),
      timingTrackTMax_(iConfig.getParameter<double>("timingTrackTMax")) {
  auto tagSiStripTrack = iConfig.getParameter<edm::InputTag>("tagSiStripTrack");
  if (!tagSiStripTrack.label().empty())
    siStripTrackToken_ = consumes<edm::DetSetVector<TotemRPLocalTrack>>(tagSiStripTrack);

  auto tagDiamondTrack = iConfig.getParameter<edm::InputTag>("tagDiamondTrack");
  if (!tagDiamondTrack.label().empty())
    diamondTrackToken_ = consumes<edm::DetSetVector<CTPPSDiamondLocalTrack>>(tagDiamondTrack);

  auto tagPixelTrack = iConfig.getParameter<edm::InputTag>("tagPixelTrack");
  if (!tagPixelTrack.label().empty())
    pixelTrackToken_ = consumes<edm::DetSetVector<CTPPSPixelLocalTrack>>(tagPixelTrack);

  produces<CTPPSLocalTrackLiteCollection>();
}

Member Function Documentation

void CTPPSLocalTrackLiteProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descr )

static

Definition at line 237 of file CTPPSLocalTrackLiteProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, HLT_2018_cff::InputTag, and edm::ParameterDescriptionNode::setComment().

                                                                                      {
  edm::ParameterSetDescription desc;

  // By default: all includeXYZ flags set to false.
  // The includeXYZ are switched on when the "ctpps_2016" era is declared in
  // python config, see:
  // RecoCTPPS/TotemRPLocal/python/ctppsLocalTrackLiteProducer_cff.py

  desc.add<bool>("includeStrips", false)->setComment("whether tracks from Si strips should be included");
  desc.add<edm::InputTag>("tagSiStripTrack", edm::InputTag("totemRPLocalTrackFitter"))
      ->setComment("input TOTEM strips' local tracks collection to retrieve");

  desc.add<bool>("includeDiamonds", false)->setComment("whether tracks from diamonds strips should be included");
  desc.add<edm::InputTag>("tagDiamondTrack", edm::InputTag("ctppsDiamondLocalTracks"))
      ->setComment("input diamond detectors' local tracks collection to retrieve");

  desc.add<bool>("includePixels", false)->setComment("whether tracks from pixels should be included");
  desc.add<edm::InputTag>("tagPixelTrack", edm::InputTag("ctppsPixelLocalTracks"))
      ->setComment("input pixel detectors' local tracks collection to retrieve");
  desc.add<double>("timingTrackTMin", -12.5)->setComment("minimal track time selection for timing detectors, in ns");
  desc.add<double>("timingTrackTMax", +12.5)->setComment("maximal track time selection for timing detectors, in ns");

  desc.add<double>("pixelTrackTxMin", -10.0);
  desc.add<double>("pixelTrackTxMax", 10.0);
  desc.add<double>("pixelTrackTyMin", -10.0);
  desc.add<double>("pixelTrackTyMax", 10.0);

  descr.add("ctppsLocalTrackLiteDefaultProducer", desc);
}

void CTPPSLocalTrackLiteProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &    )

override

Definition at line 84 of file CTPPSLocalTrackLiteProducer.cc.

References diamondTrackToken_, edm::Event::getByToken(), HPTDC_TIME_SLICE_WIDTH, includeDiamonds_, includePixels_, includeStrips_, invalid, edm::EDGetTokenT< T >::isUninitialized(), eostools::move(), pixelTrackToken_, pixelTrackTxMax_, pixelTrackTxMin_, pixelTrackTyMax_, pixelTrackTyMin_, edm::Event::put(), siStripTrackToken_, and timingTrackTMax_.

                                                                                 {
  // prepare output
  auto pOut = std::make_unique<CTPPSLocalTrackLiteCollection>();

  //----- TOTEM strips

  // get input from Si strips
  if (includeStrips_) {
    edm::Handle<edm::DetSetVector<TotemRPLocalTrack>> inputSiStripTracks;
    iEvent.getByToken(siStripTrackToken_, inputSiStripTracks);

    // process tracks from Si strips
    for (const auto &rpv : *inputSiStripTracks) {
      const uint32_t rpId = rpv.detId();
      for (const auto &trk : rpv) {
        if (!trk.isValid())
          continue;

        float roundedX0 = MiniFloatConverter::reduceMantissaToNbitsRounding<14>(trk.x0());
        float roundedX0Sigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.x0Sigma());
        float roundedY0 = MiniFloatConverter::reduceMantissaToNbitsRounding<13>(trk.y0());
        float roundedY0Sigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.y0Sigma());
        float roundedTx = MiniFloatConverter::reduceMantissaToNbitsRounding<11>(trk.tx());
        float roundedTxSigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.txSigma());
        float roundedTy = MiniFloatConverter::reduceMantissaToNbitsRounding<11>(trk.ty());
        float roundedTySigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.tySigma());
        float roundedChiSquaredOverNDF = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.chiSquaredOverNDF());

        pOut->emplace_back(rpId,  // detector info
                                  // spatial info
                           roundedX0,
                           roundedX0Sigma,
                           roundedY0,
                           roundedY0Sigma,
                           // angular info
                           roundedTx,
                           roundedTxSigma,
                           roundedTy,
                           roundedTySigma,
                           // reconstruction info
                           roundedChiSquaredOverNDF,
                           CTPPSpixelLocalTrackReconstructionInfo::invalid,
                           trk.numberOfPointsUsedForFit(),
                           // timing info
                           0.,
                           0.);
      }
    }
  }

  //----- diamond detectors

  if (includeDiamonds_) {
    // get input from diamond detectors
    edm::Handle<edm::DetSetVector<CTPPSDiamondLocalTrack>> inputDiamondTracks;
    iEvent.getByToken(diamondTrackToken_, inputDiamondTracks);

    // process tracks from diamond detectors
    for (const auto &rpv : *inputDiamondTracks) {
      const unsigned int rpId = rpv.detId();
      for (const auto &trk : rpv) {
        if (!trk.isValid())
          continue;

        const float abs_time = trk.time() + trk.ootIndex() * HPTDC_TIME_SLICE_WIDTH;
        if (abs_time < timingTrackTMin_ || abs_time > timingTrackTMax_)
          continue;

        float roundedX0 = MiniFloatConverter::reduceMantissaToNbitsRounding<16>(trk.x0());
        float roundedX0Sigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.x0Sigma());
        float roundedY0 = MiniFloatConverter::reduceMantissaToNbitsRounding<13>(trk.y0());
        float roundedY0Sigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.y0Sigma());
        float roundedT = MiniFloatConverter::reduceMantissaToNbitsRounding<16>(abs_time);
        float roundedTSigma = MiniFloatConverter::reduceMantissaToNbitsRounding<13>(trk.timeSigma());

        pOut->emplace_back(rpId,  // detector info
                                  // spatial info
                           roundedX0,
                           roundedX0Sigma,
                           roundedY0,
                           roundedY0Sigma,
                           // angular info
                           0.,
                           0.,
                           0.,
                           0.,
                           // reconstruction info
                           0.,
                           CTPPSpixelLocalTrackReconstructionInfo::invalid,
                           trk.numberOfPlanes(),
                           // timing info
                           roundedT,
                           roundedTSigma);
      }
    }
  }

  //----- pixel detectors

  if (includePixels_) {
    edm::Handle<edm::DetSetVector<CTPPSPixelLocalTrack>> inputPixelTracks;
    if (!pixelTrackToken_.isUninitialized()) {
      iEvent.getByToken(pixelTrackToken_, inputPixelTracks);

      // process tracks from pixels
      for (const auto &rpv : *inputPixelTracks) {
        const uint32_t rpId = rpv.detId();
        for (const auto &trk : rpv) {
          if (!trk.isValid())
            continue;
          if (trk.tx() > pixelTrackTxMin_ && trk.tx() < pixelTrackTxMax_ && trk.ty() > pixelTrackTyMin_ &&
              trk.ty() < pixelTrackTyMax_) {
            float roundedX0 = MiniFloatConverter::reduceMantissaToNbitsRounding<16>(trk.x0());
            float roundedX0Sigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.x0Sigma());
            float roundedY0 = MiniFloatConverter::reduceMantissaToNbitsRounding<13>(trk.y0());
            float roundedY0Sigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.y0Sigma());
            float roundedTx = MiniFloatConverter::reduceMantissaToNbitsRounding<11>(trk.tx());
            float roundedTxSigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.txSigma());
            float roundedTy = MiniFloatConverter::reduceMantissaToNbitsRounding<11>(trk.ty());
            float roundedTySigma = MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.tySigma());
            float roundedChiSquaredOverNDF =
                MiniFloatConverter::reduceMantissaToNbitsRounding<8>(trk.chiSquaredOverNDF());

            pOut->emplace_back(rpId,  // detector info
                                      // spatial info
                               roundedX0,
                               roundedX0Sigma,
                               roundedY0,
                               roundedY0Sigma,
                               // angular info
                               roundedTx,
                               roundedTxSigma,
                               roundedTy,
                               roundedTySigma,
                               // reconstruction info
                               roundedChiSquaredOverNDF,
                               trk.recoInfo(),
                               trk.numberOfPointsUsedForFit(),
                               // timing info
                               0.,
                               0.);
          }
        }
      }
    }
  }

  // save output to event
  iEvent.put(std::move(pOut));
}

Member Data Documentation

edm::EDGetTokenT<edm::DetSetVector<CTPPSDiamondLocalTrack> > CTPPSLocalTrackLiteProducer::diamondTrackToken_

private

Definition at line 46 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by CTPPSLocalTrackLiteProducer(), and produce().

float CTPPSLocalTrackLiteProducer::HPTDC_TIME_SLICE_WIDTH = 25.

staticprivate

HPTDC time slice width, in ns.

Definition at line 40 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

bool CTPPSLocalTrackLiteProducer::includeDiamonds_

private

Definition at line 45 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

bool CTPPSLocalTrackLiteProducer::includePixels_

private

Definition at line 48 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

bool CTPPSLocalTrackLiteProducer::includeStrips_

private

Definition at line 42 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

edm::EDGetTokenT<edm::DetSetVector<CTPPSPixelLocalTrack> > CTPPSLocalTrackLiteProducer::pixelTrackToken_

private

Definition at line 49 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by CTPPSLocalTrackLiteProducer(), and produce().

double CTPPSLocalTrackLiteProducer::pixelTrackTxMax_

private

Definition at line 51 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

double CTPPSLocalTrackLiteProducer::pixelTrackTxMin_

private

Definition at line 51 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

double CTPPSLocalTrackLiteProducer::pixelTrackTyMax_

private

Definition at line 51 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

double CTPPSLocalTrackLiteProducer::pixelTrackTyMin_

private

Definition at line 51 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

edm::EDGetTokenT<edm::DetSetVector<TotemRPLocalTrack> > CTPPSLocalTrackLiteProducer::siStripTrackToken_

private

Definition at line 43 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by CTPPSLocalTrackLiteProducer(), and produce().

double CTPPSLocalTrackLiteProducer::timingTrackTMax_

private

Definition at line 52 of file CTPPSLocalTrackLiteProducer.cc.

Referenced by produce().

double CTPPSLocalTrackLiteProducer::timingTrackTMin_

private

Definition at line 52 of file CTPPSLocalTrackLiteProducer.cc.