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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 58 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 264 of file CTPPSLocalTrackLiteProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, 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 89 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.getX0());
         float roundedX0Sigma =
             MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                 trk.getX0Sigma());
         float roundedY0 =
             MiniFloatConverter::reduceMantissaToNbitsRounding<13>(trk.getY0());
         float roundedY0Sigma =
             MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                 trk.getY0Sigma());
         float roundedTx =
             MiniFloatConverter::reduceMantissaToNbitsRounding<11>(trk.getTx());
         float roundedTxSigma =
             MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                 trk.getTxSigma());
         float roundedTy =
             MiniFloatConverter::reduceMantissaToNbitsRounding<11>(trk.getTy());
         float roundedTySigma =
             MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                 trk.getTySigma());
         float roundedChiSquaredOverNDF =
             MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                 trk.getChiSquaredOverNDF());
 
         pOut->emplace_back(rpId, // detector info
                                  // spatial info
                            roundedX0, roundedX0Sigma, roundedY0, roundedY0Sigma,
                            // angular info
                            roundedTx, roundedTxSigma, roundedTy, roundedTySigma,
                            // reconstruction info
                            roundedChiSquaredOverNDF,
                            CTPPSpixelLocalTrackReconstructionInfo::invalid,
                            trk.getNumberOfPointsUsedForFit(),
                            // 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.getT() + trk.getOOTIndex() * HPTDC_TIME_SLICE_WIDTH;
         if (abs_time < timingTrackTMin_ || abs_time > timingTrackTMax_)
           continue;
 
         float roundedX0 =
             MiniFloatConverter::reduceMantissaToNbitsRounding<16>(trk.getX0());
         float roundedX0Sigma =
             MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                 trk.getX0Sigma());
         float roundedY0 =
             MiniFloatConverter::reduceMantissaToNbitsRounding<13>(trk.getY0());
         float roundedY0Sigma =
             MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                 trk.getY0Sigma());
         float roundedT =
             MiniFloatConverter::reduceMantissaToNbitsRounding<16>(abs_time);
         float roundedTSigma =
             MiniFloatConverter::reduceMantissaToNbitsRounding<13>(
                 trk.getTSigma());
 
         pOut->emplace_back(rpId, // detector info
                                  // spatial info
                            roundedX0, roundedX0Sigma, roundedY0, roundedY0Sigma,
                            // angular info
                            0., 0., 0., 0.,
                            // reconstruction info
                            0., CTPPSpixelLocalTrackReconstructionInfo::invalid,
                            trk.getNumOfPlanes(),
                            // 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.getTx() > pixelTrackTxMin_ &&
               trk.getTx() < pixelTrackTxMax_ &&
               trk.getTy() > pixelTrackTyMin_ &&
               trk.getTy() < pixelTrackTyMax_) {
             float roundedX0 =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<16>(
                     trk.getX0());
             float roundedX0Sigma =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                     trk.getX0Sigma());
             float roundedY0 =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<13>(
                     trk.getY0());
             float roundedY0Sigma =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                     trk.getY0Sigma());
             float roundedTx =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<11>(
                     trk.getTx());
             float roundedTxSigma =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                     trk.getTxSigma());
             float roundedTy =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<11>(
                     trk.getTy());
             float roundedTySigma =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                     trk.getTySigma());
             float roundedChiSquaredOverNDF =
                 MiniFloatConverter::reduceMantissaToNbitsRounding<8>(
                     trk.getChiSquaredOverNDF());
 
             pOut->emplace_back(
                 rpId, // detector info
                       // spatial info
                 roundedX0, roundedX0Sigma, roundedY0, roundedY0Sigma,
                 // angular info
                 roundedTx, roundedTxSigma, roundedTy, roundedTySigma,
                 // reconstruction info
                 roundedChiSquaredOverNDF, trk.getRecoInfo(),
                 trk.getNumberOfPointsUsedForFit(),
                 // timing info
                 0., 0.);
           }
         }
       }
     }
   }
 
   // save output to event
   iEvent.put(std::move(pOut));
 }