SeedFilter Class Reference

#include <SeedFilter.h>

Classes
struct	Tokens

Public Member Functions
	SeedFilter (const edm::ParameterSet &conf, const Tokens &tokens, edm::ConsumesCollector &iC)
void	seeds (edm::Event &, const edm::EventSetup &, const reco::SuperClusterRef &, TrajectorySeedCollection *)
	~SeedFilter ()

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotTag_
SeedGeneratorFromRegionHits *	combinatorialSeedGenerator
double	deltaEta_
double	deltaPhi_
double	deta_
double	dphi_
double	dr_
double	halflength_
RectangularEtaPhiTrackingRegion::UseMeasurementTracker	hitsfactoryMode_
edm::EDGetTokenT< MeasurementTrackerEvent >	measurementTrackerToken_
double	originradius_
double	pt_
double	ptmin_
edm::ESHandle< MagneticField >	theMagField
bool	useZvertex_
edm::EDGetTokenT< std::vector< reco::Vertex > >	vertexSrc_
double	vertexz_

Detailed Description

originally Matteo Sani's SubSeedGenerator

Definition at line 32 of file SeedFilter.h.

Constructor & Destructor Documentation

SeedFilter::SeedFilter	(	const edm::ParameterSet &	conf,
		const Tokens &	tokens,
		edm::ConsumesCollector &	iC
	)

Definition at line 35 of file SeedFilter.cc.

References edm::ConsumesCollector::consumes(), beamerCreator::create(), Exception, reco::get(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), RectangularEtaPhiTrackingRegion::intToUseMeasurementTracker(), RectangularEtaPhiTrackingRegion::kNever, ckfInOutTracksFromConversions_cfi::MeasurementTrackerEvent, AlCaHLTBitMon_QueryRunRegistry::string, tmp, SeedFilter::Tokens::token_bs, and SeedFilter::Tokens::token_vtx.

 {
  edm::LogInfo("EtaPhiRegionSeedFactory") << "Enter the EtaPhiRegionSeedFactory";
  edm::ParameterSet regionPSet = conf.getParameter<edm::ParameterSet>("RegionPSet");

  ptmin_        = regionPSet.getParameter<double>("ptMin");
  originradius_ = regionPSet.getParameter<double>("originRadius");
  halflength_   = regionPSet.getParameter<double>("originHalfLength");
  deltaEta_     = regionPSet.getParameter<double>("deltaEtaRegion");
  deltaPhi_     = regionPSet.getParameter<double>("deltaPhiRegion");
  useZvertex_   = regionPSet.getParameter<bool>("useZInVertex");
  vertexSrc_    = tokens.token_vtx;

  // setup orderedhits setup (in order to tell seed generator to use pairs/triplets, which layers)
  edm::ParameterSet hitsfactoryPSet = conf.getParameter<edm::ParameterSet>("OrderedHitsFactoryPSet");
  std::string hitsfactoryName = hitsfactoryPSet.getParameter<std::string>("ComponentName");

  // HIT FACTORY MODE (JRV)
  // -1 = look for existing hit collections for both pixels and strips
  // 0 = look for pixel hit collections and build strip hits on demand (regional unpacking)
  // 1 = build both pixel and strip hits on demand (regional unpacking)
  //
  // FIXME: when next time altering the configuration of this class,
  // please consider changing the type of useOnDemandTracker to a
  // string corresponding to the UseMeasurementTracker enumeration
  int tmp = hitsfactoryPSet.getUntrackedParameter<int>("useOnDemandTracker");
  if(tmp < -1 || tmp > 1)
    throw cms::Exception("Configuration") << "SeedFilter: useOnDemandTracker must be -1, 0, or 1; got " << tmp;
  hitsfactoryMode_ = RectangularEtaPhiTrackingRegion::intToUseMeasurementTracker(tmp);

  // get orderd hits generator from factory
  OrderedHitsGenerator*  hitsGenerator = OrderedHitsGeneratorFactory::get()->create(hitsfactoryName, hitsfactoryPSet, iC);

  // start seed generator
  edm::ParameterSet seedCreatorPSet = conf.getParameter<edm::ParameterSet>("SeedCreatorPSet");
  std::string seedCreatorType = seedCreatorPSet.getParameter<std::string>("ComponentName");

  combinatorialSeedGenerator = new SeedGeneratorFromRegionHits(hitsGenerator,nullptr,
                                   SeedCreatorFactory::get()->create(seedCreatorType, seedCreatorPSet)
                                           );
  beamSpotTag_ = tokens.token_bs; ;
  if(hitsfactoryMode_ != RectangularEtaPhiTrackingRegion::UseMeasurementTracker::kNever) {
    measurementTrackerToken_ = iC.consumes<MeasurementTrackerEvent>(conf.getParameter<edm::InputTag>("measurementTrackerEvent"));
  }
 }

SeedFilter::~SeedFilter

(

)

Definition at line 83 of file SeedFilter.cc.

                        {
  delete combinatorialSeedGenerator;
}

Member Function Documentation

void SeedFilter::seeds	(	edm::Event &	e,
		const edm::EventSetup &	setup,
		const reco::SuperClusterRef &	scRef,
		TrajectorySeedCollection *	output
	)

Definition at line 87 of file SeedFilter.cc.

References funct::abs(), TrackingRecHit::all, FTSFromVertexToPointFactory::get(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, LogDebug, ckfInOutTracksFromConversions_cfi::MeasurementTrackerEvent, FreeTrajectoryState::momentum(), reco::BeamSpot::position(), edm::Handle< T >::product(), alignCSCRings::r, TrackingRecHit::recHits(), TrackingRecHit::sharesInput(), fftjetvertexadder_cfi::sigmaZ, reco::BeamSpot::sigmaZ(), reco::BeamSpot::sigmaZ0Error(), mathSSE::sqrt(), and electrons_cff::vertices.

                                                                                                                                  {

  setup.get<IdealMagneticFieldRecord>().get(theMagField);
  auto seedColl = std::make_unique<TrajectorySeedCollection>();

  GlobalPoint vtxPos;
  double deltaZVertex;

  //edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  //e.getByType(recoBeamSpotHandle);

  // gets its position
  //const reco::BeamSpot::Point& BSPosition = recoBeamSpotHandle->position();
  //vtxPos = GlobalPoint(BSPosition.x(), BSPosition.y(), BSPosition.z());
  //double sigmaZ = recoBeamSpotHandle->sigmaZ();
  //double sigmaZ0Error = recoBeamSpotHandle->sigmaZ0Error();
  //double sq=sqrt(sigmaZ*sigmaZ+sigmaZ0Error*sigmaZ0Error);
  //deltaZVertex = 3*sq; //halflength_;

  // get the primary vertex (if any)
  reco::VertexCollection vertices;
  edm::Handle<reco::VertexCollection> h_vertices;
  if (e.getByToken(vertexSrc_, h_vertices)) {
    vertices = *(h_vertices.product());
  } else {
      LogDebug("SeedFilter") << "SeedFilter::seeds"
                  << "No vertex collection found: using beam-spot";
  }

  if (!vertices.empty() && useZvertex_) {
    vtxPos = GlobalPoint(vertices.front().x(), vertices.front().y(), vertices.front().z());
    deltaZVertex = halflength_;

  } else {
    edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
    //e.getByType(recoBeamSpotHandle);
    e.getByToken(beamSpotTag_,recoBeamSpotHandle);
    // gets its position
    const reco::BeamSpot::Point& BSPosition = recoBeamSpotHandle->position();
    double sigmaZ = recoBeamSpotHandle->sigmaZ();
    double sigmaZ0Error = recoBeamSpotHandle->sigmaZ0Error();
    double sq=sqrt(sigmaZ*sigmaZ+sigmaZ0Error*sigmaZ0Error);
    vtxPos = GlobalPoint(BSPosition.x(), BSPosition.y(), BSPosition.z());
    deltaZVertex = 3*sq;
  }

  const MeasurementTrackerEvent *measurementTracker = nullptr;
  if(!measurementTrackerToken_.isUninitialized()) {
    edm::Handle<MeasurementTrackerEvent> hmte;
    e.getByToken(measurementTrackerToken_, hmte);
    measurementTracker = hmte.product();
  }

  //seeds selection
  float energy = scRef->energy();

  const GlobalPoint clusterPos(scRef->position().x(), scRef->position().y(), scRef->position().z());

  //===============================================
  // EtaPhiRegion and seeds for electron hypothesis
  //===============================================

  TrackCharge aCharge = -1 ;
  FreeTrajectoryState fts = FTSFromVertexToPointFactory::get(*theMagField, clusterPos, vtxPos, energy, aCharge);

  RectangularEtaPhiTrackingRegion etaphiRegionMinus(fts.momentum(),
                                                    vtxPos,
                                                    ptmin_,
                                                    originradius_,
                                                    deltaZVertex,
                            //                                                    deltaEta_,
                            //                                                    deltaPhi_,-1);
                            RectangularEtaPhiTrackingRegion::Margin(std::abs(deltaEta_),std::abs(deltaEta_)),
                            RectangularEtaPhiTrackingRegion::Margin(std::abs(deltaPhi_),std::abs(deltaPhi_))
                            ,hitsfactoryMode_,
                            true, /*default in header*/
                            measurementTracker
                            );
  combinatorialSeedGenerator->run(*seedColl, etaphiRegionMinus, e, setup);

  for (unsigned int i = 0; i<seedColl->size(); ++i)
    output->push_back((*seedColl)[i]);

  seedColl->clear();

  //===============================================
  // EtaPhiRegion and seeds for positron hypothesis
  //===============================================

  TrackCharge aChargep = 1 ;
  fts = FTSFromVertexToPointFactory::get(*theMagField, clusterPos, vtxPos, energy, aChargep);

  RectangularEtaPhiTrackingRegion etaphiRegionPlus(fts.momentum(),
                                                   vtxPos,
                                                   ptmin_,
                                                   originradius_,
                                                   deltaZVertex,
                           //                                                   deltaEta_,
                           //                                                   deltaPhi_,-1);
                            RectangularEtaPhiTrackingRegion::Margin(std::abs(deltaEta_),std::abs(deltaEta_)),
                            RectangularEtaPhiTrackingRegion::Margin(std::abs(deltaPhi_),std::abs(deltaPhi_))
                            ,hitsfactoryMode_
                           ,true, /*default in header*/
                           measurementTracker
                           );

  combinatorialSeedGenerator->run(*seedColl, etaphiRegionPlus, e, setup);

  for (unsigned int i = 0; i<seedColl->size(); ++i){
    TrajectorySeed::range r = (*seedColl)[i].recHits();
    // first Hit
    TrajectorySeed::const_iterator it = r.first;
    const TrackingRecHit* a1 = &(*it);
    // now second Hit
    it++;
    const TrackingRecHit* a2 = &(*it);
    // check if the seed is already in the collection
    bool isInCollection = false;
    for(unsigned int j=0; j<output->size(); ++j) {
      TrajectorySeed::range r = (*seedColl)[i].recHits();
      // first Hit
      TrajectorySeed::const_iterator it = r.first;
      const TrackingRecHit* b1 = &(*it);
      // now second Hit
      it++;
      const TrackingRecHit* b2 = &(*it);
      if ((b1->sharesInput(a1, TrackingRecHit::all)) && (b2->sharesInput(a2, TrackingRecHit::all))) {
    isInCollection = true;
    break;
      }
    }
    if (!isInCollection)
      output->push_back((*seedColl)[i]);
  }

  seedColl->clear();
}

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> SeedFilter::beamSpotTag_

private

Definition at line 62 of file SeedFilter.h.

SeedGeneratorFromRegionHits* SeedFilter::combinatorialSeedGenerator

private

Definition at line 48 of file SeedFilter.h.

double SeedFilter::deltaEta_

private

Definition at line 53 of file SeedFilter.h.

double SeedFilter::deltaPhi_

private

Definition at line 53 of file SeedFilter.h.

double SeedFilter::deta_

private

Definition at line 51 of file SeedFilter.h.

double SeedFilter::dphi_

private

Definition at line 51 of file SeedFilter.h.

double SeedFilter::dr_

private

Definition at line 51 of file SeedFilter.h.

double SeedFilter::halflength_

private

Definition at line 53 of file SeedFilter.h.

RectangularEtaPhiTrackingRegion::UseMeasurementTracker SeedFilter::hitsfactoryMode_

private

Definition at line 60 of file SeedFilter.h.

edm::EDGetTokenT<MeasurementTrackerEvent> SeedFilter::measurementTrackerToken_

private

Definition at line 64 of file SeedFilter.h.

double SeedFilter::originradius_

private

Definition at line 53 of file SeedFilter.h.

double SeedFilter::pt_

private

Definition at line 51 of file SeedFilter.h.

double SeedFilter::ptmin_

private

Definition at line 53 of file SeedFilter.h.

edm::ESHandle<MagneticField> SeedFilter::theMagField

private

Definition at line 58 of file SeedFilter.h.

bool SeedFilter::useZvertex_

private

Definition at line 54 of file SeedFilter.h.

edm::EDGetTokenT<std::vector<reco::Vertex> > SeedFilter::vertexSrc_

private

Definition at line 56 of file SeedFilter.h.

double SeedFilter::vertexz_

private

Definition at line 53 of file SeedFilter.h.