TotemRPLocalTrackFitter Class Reference

Fits tracks trough a single RP. More...

Inheritance diagram for TotemRPLocalTrackFitter:

Public Member Functions
void	produce (edm::Event &e, const edm::EventSetup &c) override
	TotemRPLocalTrackFitter (const edm::ParameterSet &conf)
	~TotemRPLocalTrackFitter () 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
TotemRPLocalTrackFitterAlgorithm	fitter_
	The instance of the fitter module. More...
edm::ESWatcher< VeryForwardRealGeometryRecord >	geometryWatcher
	A watcher to detect geometry changes. More...
edm::EDGetTokenT< edm::DetSetVector< TotemRPUVPattern > >	patternCollectionToken
edm::InputTag	tagUVPattern
	Selection of the pattern-recognition module. More...
int	verbosity_

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

Fits tracks trough a single RP.

Definition at line 34 of file TotemRPLocalTrackFitter.cc.

Constructor & Destructor Documentation

TotemRPLocalTrackFitter::TotemRPLocalTrackFitter ( const edm::ParameterSet &  conf )

explicit

Definition at line 67 of file TotemRPLocalTrackFitter.cc.

References edm::ParameterSet::getParameter(), patternCollectionToken, and tagUVPattern.

   : verbosity_(conf.getParameter<int>("verbosity")), fitter_(conf)
{
  tagUVPattern = conf.getParameter<edm::InputTag>("tagUVPattern");
  patternCollectionToken = consumes<DetSetVector<TotemRPUVPattern>>(tagUVPattern);

  produces<DetSetVector<TotemRPLocalTrack>>();
}

TotemRPLocalTrackFitter::~TotemRPLocalTrackFitter ( )

inlineoverride

Definition at line 39 of file TotemRPLocalTrackFitter.cc.

References EnergyCorrector::c, MillePedeFileConverter_cfg::e, fillDescriptions(), and produce().

{}

Member Function Documentation

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

static

Definition at line 189 of file TotemRPLocalTrackFitter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and DEFINE_FWK_MODULE.

Referenced by ~TotemRPLocalTrackFitter().

{
  edm::ParameterSetDescription desc;

  desc.add<edm::InputTag>( "tagUVPattern", edm::InputTag( "totemRPUVPatternFinder" ) )
    ->setComment( "input U-V patterns collection to retrieve" );
  desc.add<int>( "verbosity", 0 );

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

void TotemRPLocalTrackFitter::produce ( edm::Event &  e, const edm::EventSetup &  c  )

override

Definition at line 78 of file TotemRPLocalTrackFitter.cc.

References edm::ESWatcher< T >::check(), edm::DetSetVector< T >::find_or_insert(), fitter_, TotemRPLocalTrackFitterAlgorithm::fitTrack(), geometry, geometryWatcher, edm::EventSetup::get(), edm::Event::getByToken(), TotemRPLocalTrack::getHits(), TotemRPUVPattern::getProjection(), CTPPSGeometry::getRPTranslation(), h, hfClusterShapes_cfi::hits, photons_cff::ids, input, TotemRPLocalTrack::isValid(), convertSQLitetoXML_cfg::output, listBenchmarks::pattern, patternCollectionToken, pi, TotemRPUVPattern::projU, TotemRPUVPattern::projV, edm::DetSet< T >::push_back(), edm::Event::put(), TotemRPLocalTrackFitterAlgorithm::reset(), HiIsolationCommonParameters_cff::track, and verbosity_.

Referenced by ~TotemRPLocalTrackFitter().

{
  if (verbosity_ > 5)
    LogVerbatim("TotemRPLocalTrackFitter") << ">> TotemRPLocalTrackFitter::produce";

  // get geometry
  edm::ESHandle<CTPPSGeometry> geometry;
  setup.get<VeryForwardRealGeometryRecord>().get(geometry);

  if (geometryWatcher.check(setup))
    fitter_.reset();
  
  // get input
  edm::Handle<DetSetVector<TotemRPUVPattern>> input;
  e.getByToken(patternCollectionToken, input);

  // run fit for each RP
  DetSetVector<TotemRPLocalTrack> output;

  for (const auto &rpv : *input)
  {
    CTPPSDetId rpId(rpv.detId());

    // is U-V association unique?
    unsigned int n_U=0, n_V=0;
    unsigned int idx_U=0, idx_V=0;
    for (unsigned int pi = 0; pi < rpv.size(); pi++)
    {
      const TotemRPUVPattern &pattern = rpv[pi];

      // here it would make sense to skip non-fittable patterns, but to keep the logic
      // equivalent to version 7_0_4, nothing is skipped
      /*
      if (pattern.getFittable() == false)
        continue;
      */

      switch (pattern.getProjection())
      {
        case TotemRPUVPattern::projU:
          n_U++;
          idx_U=pi;
          break;
     
        case TotemRPUVPattern::projV:
          n_V++;
          idx_V=pi;
          break;

        default:
          break;
      }   
    }

    if (n_U != 1 || n_V != 1)
    {
      if (verbosity_)
        LogVerbatim("TotemRPLocalTrackFitter")
          << ">> TotemRPLocalTrackFitter::produce > Impossible to combine U and V patterns in RP " << rpId
          << " (n_U=" << n_U << ", n_V=" << n_V << ").";

      continue;
    }
    
    // again, to follow the logic from version 7_0_4, skip the non-fittable patterns here
    if (!rpv[idx_U].getFittable() || !rpv[idx_V].getFittable())
      continue;

    // combine U and V hits
    DetSetVector<TotemRPRecHit> hits;
    for (auto &ids : rpv[idx_U].getHits())
    {
      auto &ods = hits.find_or_insert(ids.detId());
      for (auto &h : ids)
        ods.push_back(h);
    }

    for (auto &ids : rpv[idx_V].getHits())
    {
      auto &ods = hits.find_or_insert(ids.detId());
      for (auto &h : ids)
        ods.push_back(h);
    }

    // run fit
    double z0 = geometry->getRPTranslation(rpId).z();

    TotemRPLocalTrack track;
    fitter_.fitTrack(hits, z0, *geometry, track);
    
    DetSet<TotemRPLocalTrack> &ds = output.find_or_insert(rpId);
    ds.push_back(track);

    if (verbosity_ > 5)
    {
      unsigned int n_hits = 0;
      for (auto &hds : track.getHits())
        n_hits += hds.size();

      LogVerbatim("TotemRPLocalTrackFitter")
        << "    track in RP " << rpId << ": valid = " << track.isValid() << ", hits = " << n_hits;
    }
  }

  // save results
  e.put(make_unique<DetSetVector<TotemRPLocalTrack>>(output));
}

Member Data Documentation

TotemRPLocalTrackFitterAlgorithm TotemRPLocalTrackFitter::fitter_

private

The instance of the fitter module.

Definition at line 56 of file TotemRPLocalTrackFitter.cc.

Referenced by produce().

edm::ESWatcher<VeryForwardRealGeometryRecord> TotemRPLocalTrackFitter::geometryWatcher

private

A watcher to detect geometry changes.

Definition at line 53 of file TotemRPLocalTrackFitter.cc.

Referenced by produce().

edm::EDGetTokenT<edm::DetSetVector<TotemRPUVPattern> > TotemRPLocalTrackFitter::patternCollectionToken

private

Definition at line 50 of file TotemRPLocalTrackFitter.cc.

Referenced by produce(), and TotemRPLocalTrackFitter().

edm::InputTag TotemRPLocalTrackFitter::tagUVPattern

private

Selection of the pattern-recognition module.

Definition at line 48 of file TotemRPLocalTrackFitter.cc.

Referenced by TotemRPLocalTrackFitter().

int TotemRPLocalTrackFitter::verbosity_

private

Definition at line 45 of file TotemRPLocalTrackFitter.cc.

Referenced by produce().