BeamHaloPairGenerator Class Reference

#include <BeamHaloPairGenerator.h>

Inheritance diagram for BeamHaloPairGenerator:

Public Member Functions
	BeamHaloPairGenerator (const edm::ParameterSet &conf, edm::ConsumesCollector &iC)
virtual const OrderedSeedingHits &	run (const TrackingRegion &region, const edm::Event &ev, const edm::EventSetup &es)
virtual	~BeamHaloPairGenerator ()
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	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	iC
	)

Definition at line 10 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));
} 

virtual BeamHaloPairGenerator::~BeamHaloPairGenerator ( )

inlinevirtual

Definition at line 19 of file BeamHaloPairGenerator.h.

{};

Member Function Documentation

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

virtual

Implements OrderedHitsGenerator.

Definition at line 19 of file BeamHaloPairGenerator.cc.

References edm::Event::getByToken(), hitPairs, TrackingRegion::hits(), LayerTriplets::layers(), LogDebug, python.rootplot.utilities::ls(), SeedingLayerSetsHits::numberOfLayersInSet(), funct::sin(), theMaxTheta, theSeedingLayerToken, and PV3DBase< T, PVType, FrameType >::theta().

                                                                    {
    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(e, es, ls[0]);
        auto outerHits  = region.hits(e, es, 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

OrderedHitPairs BeamHaloPairGenerator::hitPairs

private

Definition at line 25 of file BeamHaloPairGenerator.h.

Referenced by run().

double BeamHaloPairGenerator::theMaxTheta

private

Definition at line 26 of file BeamHaloPairGenerator.h.

Referenced by BeamHaloPairGenerator(), and run().

edm::EDGetTokenT<SeedingLayerSetsHits> BeamHaloPairGenerator::theSeedingLayerToken

private

Definition at line 24 of file BeamHaloPairGenerator.h.

Referenced by run().