BeamHaloPairGenerator Class Reference

#include <BeamHaloPairGenerator.h>

Inheritance diagram for BeamHaloPairGenerator:

Public Member Functions
	BeamHaloPairGenerator (const edm::ParameterSet &conf, edm::ConsumesCollector &iC)
const OrderedSeedingHits &	run (const TrackingRegion &region, const edm::Event &ev, const edm::EventSetup &es) override
	~BeamHaloPairGenerator () override
Public Member Functions inherited from OrderedHitsGenerator
virtual void	clear ()
	OrderedHitsGenerator ()
virtual	~OrderedHitsGenerator ()

Private Attributes
OrderedHitPairs	hitPairs
double	theMaxTheta
edm::EDGetTokenT< SeedingLayerSetsHits >	theSeedingLayerToken

Additional Inherited Members
Public Attributes inherited from OrderedHitsGenerator
unsigned int	theMaxElement

Detailed Description

Definition at line 16 of file BeamHaloPairGenerator.h.

Constructor & Destructor Documentation

BeamHaloPairGenerator()

BeamHaloPairGenerator::BeamHaloPairGenerator	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	iC
	)

Definition at line 8 of file BeamHaloPairGenerator.cc.

References edm::ParameterSet::getParameter(), funct::sin(), and theMaxTheta.

    : theSeedingLayerToken(iC.consumes<SeedingLayerSetsHits>(conf.getParameter<edm::InputTag>("LayerSrc"))) {
  edm::LogInfo("CtfSpecialSeedGenerator|BeamHaloPairGenerator") << "Constructing BeamHaloPairGenerator";
  theMaxTheta = conf.getParameter<double>("maxTheta");
  theMaxTheta = fabs(sin(theMaxTheta));
}

~BeamHaloPairGenerator()

BeamHaloPairGenerator::~BeamHaloPairGenerator ( )

inlineoverride

Definition at line 19 of file BeamHaloPairGenerator.h.

{};

Member Function Documentation

run()

const OrderedSeedingHits & BeamHaloPairGenerator::run ( const TrackingRegion &  region, const edm::Event &  ev, const edm::EventSetup &  es  )

overridevirtual

Implements OrderedHitsGenerator.

Definition at line 15 of file BeamHaloPairGenerator.cc.

References ztail::d, MillePedeFileConverter_cfg::e, hitPairs, hgcalTopologyTester_cfi::layers, LogDebug, eostools::ls(), HLT_2022v12_cff::region, funct::sin(), theMaxTheta, and theSeedingLayerToken.

                                                                              {
  hitPairs.clear();
  edm::Handle<SeedingLayerSetsHits> hlayers;
  e.getByToken(theSeedingLayerToken, hlayers);
  const SeedingLayerSetsHits& layers = *hlayers;
  if (layers.numberOfLayersInSet() != 2)
    throw cms::Exception("CtfSpecialSeedGenerator")
        << "You are using " << layers.numberOfLayersInSet() << " layers in set instead of 2 ";
  for (SeedingLayerSetsHits::SeedingLayerSet ls : layers) {
    auto innerHits = region.hits(ls[0]);
    auto outerHits = region.hits(ls[1]);

    for (auto iOuterHit = outerHits.begin(); iOuterHit != outerHits.end(); iOuterHit++) {
      for (auto iInnerHit = innerHits.begin(); iInnerHit != innerHits.end(); iInnerHit++) {
        //do something in there... if necessary
        SeedingHitSet::ConstRecHitPointer crhpi = &(**iInnerHit);
        SeedingHitSet::ConstRecHitPointer crhpo = &(**iOuterHit);
        GlobalVector d = crhpo->globalPosition() - crhpi->globalPosition();
        double ABSsinDtheta = fabs(sin(d.theta()));
        LogDebug("BeamHaloPairGenerator")
            << "position1: " << crhpo->globalPosition() << " position2: " << crhpi->globalPosition()
            << " |sin(Dtheta)|: " << ABSsinDtheta << ((ABSsinDtheta > theMaxTheta) ? " skip" : " keep");

        if (ABSsinDtheta > theMaxTheta) {
          continue;
        }

        hitPairs.push_back(OrderedHitPair(crhpi, crhpo));
      }
    }
  }
  return hitPairs;
}

Member Data Documentation

hitPairs

OrderedHitPairs BeamHaloPairGenerator::hitPairs

private

Definition at line 24 of file BeamHaloPairGenerator.h.

Referenced by run().

theMaxTheta

double BeamHaloPairGenerator::theMaxTheta

private

Definition at line 25 of file BeamHaloPairGenerator.h.

Referenced by BeamHaloPairGenerator(), and run().

theSeedingLayerToken

edm::EDGetTokenT<SeedingLayerSetsHits> BeamHaloPairGenerator::theSeedingLayerToken

private

Definition at line 23 of file BeamHaloPairGenerator.h.

Referenced by run().