#include <PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h>

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup)

TrajectorySeedCollection *	getTrajectorySeedCollection ()

	PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

	~PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo ()

Private Types
typedef Regions::const_iterator	IR

typedef std::vector< std::unique_ptr< TrackingRegion > >	Regions

Private Member Functions
bool	inspect (const TrackingRegion &region)

void	loop ()

Private Attributes
const edm::ParameterSet	_conf

const edm::EventSetup *	myEsetup

const edm::Event *	myEvent

PrintRecoObjects	po

reco::Vertex	primaryVertex

edm::ParameterSet	QuadCutPSet

Regions	regions

TrajectorySeedCollection *	seedCollection

std::stringstream	ss

ClusterChecker	theClusterCheck

std::unique_ptr< CombinedHitQuadrupletGeneratorForPhotonConversion >	theHitsGenerator

std::unique_ptr< GlobalTrackingRegionProducerFromBeamSpot >	theRegionProducer

std::unique_ptr< SeedForPhotonConversionFromQuadruplets >	theSeedCreator

bool	theSilentOnClusterCheck

edm::EDGetTokenT< reco::VertexCollection >	token_vertex

reco::VertexCollection	vertexCollection

edm::Handle< reco::VertexCollection >	vertexHandle

Detailed Description

Definition at line 28 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Member Typedef Documentation

◆ IR

typedef Regions::const_iterator PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::IR

private

Definition at line 63 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

◆ Regions

typedef std::vector<std::unique_ptr<TrackingRegion> > PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::Regions

private

Definition at line 62 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Constructor & Destructor Documentation

◆ PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo()

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

Definition at line 21 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

     : _conf(conf),
       seedCollection(nullptr),
       theClusterCheck(conf.getParameter<edm::ParameterSet>("ClusterCheckPSet"), iC),
       QuadCutPSet(conf.getParameter<edm::ParameterSet>("QuadCutPSet")),
       theSilentOnClusterCheck(conf.getParameter<edm::ParameterSet>("ClusterCheckPSet")
                                   .getUntrackedParameter<bool>("silentClusterCheck", false)),
       theHitsGenerator(new CombinedHitQuadrupletGeneratorForPhotonConversion(
           conf.getParameter<edm::ParameterSet>("OrderedHitsFactoryPSet"), iC)),
       theSeedCreator(
           new SeedForPhotonConversionFromQuadruplets(conf.getParameter<edm::ParameterSet>("SeedCreatorPSet"),
                                                      iC,
                                                      conf.getParameter<edm::ParameterSet>("SeedComparitorPSet"))),
       theRegionProducer(
           new GlobalTrackingRegionProducerFromBeamSpot(conf.getParameter<edm::ParameterSet>("RegionFactoryPSet"), iC)) {
   token_vertex = iC.consumes<reco::VertexCollection>(_conf.getParameter<edm::InputTag>("primaryVerticesTag"));
 }

References _conf, edm::ParameterSet::getParameter(), and token_vertex.

◆ ~PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo()

PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::~PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo ( )

Definition at line 40 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

40 {

41 }

Member Function Documentation

◆ analyze()

void PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

Definition at line 43 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

                                                                                                     {
   myEsetup = &setup;
   myEvent = &event;
  
   if (seedCollection != nullptr)
     delete seedCollection;
  
   seedCollection = new TrajectorySeedCollection();
  
   size_t clustsOrZero = theClusterCheck.tooManyClusters(event);
   if (clustsOrZero) {
     if (!theSilentOnClusterCheck)
       edm::LogError("TooManyClusters") << "Found too many clusters (" << clustsOrZero << "), bailing out.\n";
     return;
   }
  
   // Why is the regions variable (and theRegionProducer) needed at all?
   regions = theRegionProducer->regions(event, setup);
  
   event.getByToken(token_vertex, vertexHandle);
   if (!vertexHandle.isValid()) {
     edm::LogError("PhotonConversionFinderFromTracks")
         << "Error! Can't get the product primary Vertex Collection "
         << _conf.getParameter<edm::InputTag>("primaryVerticesTag") << "\n";
     return;
   }
  
   //Do the analysis
   loop();
  
 #ifdef mydebug_knuenz
   std::cout << "Running PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo" << std::endl;
 #endif
  
 #ifdef debugTSPFSLA
   std::stringstream ss;
   ss.str("");
   ss << "\n++++++++++++++++++\n";
   ss << "seed collection size " << seedCollection->size();
   for (auto const& tjS : *seedCollection) {
     po.print(ss, tjS);
   }
   edm::LogInfo("debugTrajSeedFromQuadruplets") << ss.str();
   //-------------------------------------------------
 #endif
 }

References _conf, gather_cfg::cout, edmPickEvents::event, edm::ParameterSet::getParameter(), edm::HandleBase::isValid(), loop(), myEsetup, myEvent, po, PrintRecoObjects::print(), regions, seedCollection, singleTopDQM_cfi::setup, ss, theClusterCheck, theRegionProducer, theSilentOnClusterCheck, token_vertex, ClusterChecker::tooManyClusters(), and vertexHandle.

◆ getTrajectorySeedCollection()

TrajectorySeedCollection* PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::getTrajectorySeedCollection ( )

inline

Definition at line 34 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

34 { return seedCollection; }

References seedCollection.

◆ inspect()

bool PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::inspect ( const TrackingRegion & region )

private

Definition at line 116 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

                                                                                                     {
   const OrderedSeedingHits& hitss = theHitsGenerator->run(region, *myEvent, *myEsetup);
  
   unsigned int nHitss = hitss.size();
  
 #ifdef debugTSPFSLA
   ss << "\n nHitss " << nHitss << "\n";
 #endif
  
   // don't do multiple reserves in the case of multiple regions: it would make things even worse
   // as it will cause N re-allocations instead of the normal log(N)/log(2)
   if (seedCollection->empty())
     seedCollection->reserve(nHitss / 2);
  
   unsigned int iHits = 0, jHits = 1;
  
   //
   // Trivial arbitration
   //
   // Vector to store old quad values
   std::vector<Quad> quadVector;
  
   for (; iHits < nHitss && jHits < nHitss; iHits += 2, jHits += 2) {
 #ifdef debugTSPFSLA
     //    ss << "\n iHits " << iHits << " " << jHits << "\n";
 #endif
     //phits is the pair of hits for the positron
     const SeedingHitSet& phits = hitss[iHits];
     //mhits is the pair of hits for the electron
     const SeedingHitSet& mhits = hitss[jHits];
  
     try {
       //FIXME (modify the interface of the seed generator if needed)
       //passing the region, that is centered around the primary vertex
       theSeedCreator->trajectorySeed(
           *seedCollection, phits, mhits, region, *myEvent, *myEsetup, ss, quadVector, QuadCutPSet);
     } catch (cms::Exception& er) {
       edm::LogError("SeedingConversion") << " Problem in the Quad Seed creator " << er.what() << std::endl;
     } catch (std::exception& er) {
       edm::LogError("SeedingConversion") << " Problem in the Quad Seed creator " << er.what() << std::endl;
     }
   }
   quadVector.clear();
   return true;
 }

References cppFunctionSkipper::exception, myEsetup, myEvent, phits, QuadCutPSet, HLT_FULL_cff::region, seedCollection, OrderedSeedingHits::size(), ss, theHitsGenerator, theSeedCreator, and cms::Exception::what().

Referenced by loop().

◆ loop()

void PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::loop ( )

private

Definition at line 91 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

                                                                      {
   ss.str("");
  
   float ptMin = 0.1;
  
   for (auto const& primaryVertex : *vertexHandle) {
     //FIXME currnetly using the first primaryVertex, should loop on the promaryVertexes
     GlobalTrackingRegion region(
         ptMin,
         GlobalPoint(primaryVertex.position().x(), primaryVertex.position().y(), primaryVertex.position().z()),
         primaryVertex.xError() * 10,
         primaryVertex.zError() * 10,
         true);
  
 #ifdef debugTSPFSLA
     ss << "[PrintRegion] " << region.print() << std::endl;
 #endif
  
     inspect(region);
   }
 #ifdef debugTSPFSLA
   edm::LogInfo("debugTrajSeedFromQuadruplets") << ss.str();
 #endif
 }