PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo Class Reference

#include <PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h>

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup)
void	clear ()
TrajectorySeedCollection *	getTrajectorySeedCollection ()
void	init ()
	PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo (const edm::ParameterSet &, edm::ConsumesCollector &&iC)
	~PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo ()

Private Types
typedef Regions::const_iterator	IR
typedef std::vector < TrackingRegion * >	Regions

Private Member Functions
bool	inspect (const TrackingRegion &region)
void	loop ()

Private Attributes
const edm::ParameterSet	_conf
edm::ParameterSet	creatorPSet
const edm::EventSetup *	myEsetup
const edm::Event *	myEvent
PrintRecoObjects	po
reco::Vertex	primaryVertex
edm::ParameterSet	QuadCutPSet
Regions	regions
TrajectorySeedCollection *	seedCollection
edm::ParameterSet	SeedComparitorPSet
std::stringstream	ss
ClusterChecker	theClusterCheck
std::unique_ptr < CombinedHitQuadrupletGeneratorForPhotonConversion >	theHitsGenerator
std::unique_ptr < GlobalTrackingRegionProducerFromBeamSpot >	theRegionProducer
SeedForPhotonConversionFromQuadruplets *	theSeedCreator
bool	theSilentOnClusterCheck
edm::EDGetTokenT < reco::VertexCollection >	token_vertex
reco::VertexCollection	vertexCollection
edm::Handle < reco::VertexCollection >	vertexHandle

Detailed Description

Definition at line 30 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Member Typedef Documentation

typedef Regions::const_iterator PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::IR

private

Definition at line 73 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

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

private

Definition at line 72 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

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

  :_conf(conf),seedCollection(0),
   creatorPSet(conf.getParameter<edm::ParameterSet>("SeedCreatorPSet")),
   theClusterCheck(conf.getParameter<edm::ParameterSet>("ClusterCheckPSet"), iC),
   SeedComparitorPSet(conf.getParameter<edm::ParameterSet>("SeedComparitorPSet")),
   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)),
   theRegionProducer(new GlobalTrackingRegionProducerFromBeamSpot(conf.getParameter<edm::ParameterSet>("RegionFactoryPSet"), iC)) {

  token_vertex      = iC.consumes<reco::VertexCollection>(_conf.getParameter<edm::InputTag>("primaryVerticesTag"));

  init();  
}

PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::~PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo

(

)

Definition at line 38 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

                                                                                                                      {
}

Member Function Documentation

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

Definition at line 53 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

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

Referenced by PhotonConversionTrajectorySeedProducerFromQuadruplets::produce().

                                                           {

  myEsetup = &setup;
  myEvent = &event;

  if(seedCollection!=0)
    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 ;
  }

  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();
  BOOST_FOREACH(TrajectorySeed tjS,*seedCollection){
    po.print(ss, tjS);
  }
  edm::LogInfo("debugTrajSeedFromQuadruplets") << ss.str();
  //-------------------------------------------------
#endif
}

void PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::clear

(

void

)

Definition at line 42 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

References NULL, and theSeedCreator.

Referenced by PhotonConversionTrajectorySeedProducerFromQuadruplets::endRun().

       {
  if(theSeedCreator!=NULL)
    delete theSeedCreator;
}

TrajectorySeedCollection* PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::getTrajectorySeedCollection ( )

inline

Definition at line 41 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

References seedCollection.

Referenced by PhotonConversionTrajectorySeedProducerFromQuadruplets::produce().

{return seedCollection;}

void PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::init

(

void

)

Definition at line 48 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

References creatorPSet, and theSeedCreator.

Referenced by PhotonConversionTrajectorySeedProducerFromQuadruplets::beginRun(), and PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo().

      {
  theSeedCreator    = new SeedForPhotonConversionFromQuadruplets(creatorPSet);
}

bool PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::inspect ( const TrackingRegion &  region )

private

Definition at line 135 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

References cppFunctionSkipper::exception, myEsetup, myEvent, QuadCutPSet, seedCollection, SeedComparitorPSet, OrderedSeedingHits::size(), ss, theHitsGenerator, theSeedCreator, SeedForPhotonConversionFromQuadruplets::trajectorySeed(), and cms::Exception::what().

Referenced by loop().

                                       {

  const OrderedSeedingHits & hitss = theHitsGenerator->run(region, *myEvent, *myEsetup);
  
  unsigned int nHitss =  hitss.size();

#ifdef debugTSPFSLA 
  ss << "\n nHitss " << nHitss << "\n";
#endif

  if (seedCollection->empty()) seedCollection->reserve(nHitss/2); // 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)

  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, *myEsetup, ss, quadVector, SeedComparitorPSet, 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;
}

void PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::loop ( )

private

Definition at line 100 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.cc.

References inspect(), reco::Vertex::position(), primaryVertex, PtMinSelector_cfg::ptMin, ss, funct::true, vertexHandle, reco::Vertex::xError(), and reco::Vertex::zError().

Referenced by analyze().

      {

  
  ss.str("");
  
  float ptMin=0.1;
  
  BOOST_FOREACH(const reco::Vertex 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
}

Member Data Documentation

const edm::ParameterSet PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::_conf

private

Definition at line 59 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze(), and PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo().

edm::ParameterSet PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::creatorPSet

private

Definition at line 62 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by init().

const edm::EventSetup* PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::myEsetup

private

Definition at line 81 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze(), and inspect().

const edm::Event* PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::myEvent

private

Definition at line 82 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze(), and inspect().

PrintRecoObjects PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::po

private

Definition at line 84 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze().

reco::Vertex PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::primaryVertex

private

Definition at line 78 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by loop().

edm::ParameterSet PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::QuadCutPSet

private

Definition at line 64 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by inspect().

Regions PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::regions

private

Definition at line 74 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze().

TrajectorySeedCollection* PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::seedCollection

private

Definition at line 61 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze(), getTrajectorySeedCollection(), and inspect().

edm::ParameterSet PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::SeedComparitorPSet

private

Definition at line 64 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by inspect().

std::stringstream PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::ss

private

Definition at line 86 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze(), inspect(), and loop().

ClusterChecker PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::theClusterCheck

private

Definition at line 63 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze().

std::unique_ptr<CombinedHitQuadrupletGeneratorForPhotonConversion> PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::theHitsGenerator

private

Definition at line 67 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by inspect().

std::unique_ptr<GlobalTrackingRegionProducerFromBeamSpot> PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::theRegionProducer

private

Definition at line 69 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze().

SeedForPhotonConversionFromQuadruplets* PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::theSeedCreator

private

Definition at line 68 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by clear(), init(), and inspect().

bool PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::theSilentOnClusterCheck

private

Definition at line 65 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze().

edm::EDGetTokenT<reco::VertexCollection> PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::token_vertex

private

Definition at line 79 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze(), and PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo().

reco::VertexCollection PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::vertexCollection

private

Definition at line 77 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

edm::Handle<reco::VertexCollection> PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::vertexHandle

private

Definition at line 76 of file PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo.h.

Referenced by analyze(), and loop().