ElectronSeedGenerator Class Reference

#include <ElectronSeedGenerator.h>

Classes
struct	Tokens

Public Types
typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer
typedef edm::OwnVector< TrackingRecHit >	PRecHitContainer
typedef TransientTrackingRecHit::RecHitContainer	RecHitContainer
typedef TransientTrackingRecHit::RecHitPointer	RecHitPointer

Public Member Functions
	ElectronSeedGenerator (const edm::ParameterSet &, const Tokens &)
void	run (edm::Event &, const edm::EventSetup &setup, const reco::SuperClusterRefVector &, const std::vector< const TrajectorySeedCollection * > &seedsV, reco::ElectronSeedCollection &)
void	setupES (const edm::EventSetup &setup)

Private Member Functions
void	seedsFromThisCluster (edm::Ref< reco::SuperClusterCollection > seedCluster, reco::BeamSpot const &beamSpot, std::vector< reco::Vertex > const *vertices, reco::ElectronSeedCollection &out)

Private Attributes
const edm::EDGetTokenT< reco::BeamSpot >	beamSpotTag_
unsigned long long	cacheIDCkfComp_ = 0
unsigned long long	cacheIDMagField_ = 0
unsigned long long	cacheIDTrkGeom_ = 0
const float	deltaPhi1High_
const float	deltaPhi1Low_
const float	deltaPhi2B_
const float	deltaPhi2F_
const float	deltaZ1WithVertex_
const double	dPhi1Coef1_
const double	dPhi1Coef2_
const bool	dynamicPhiRoad_
PixelHitMatcher	electronMatcher_
const float	highPtThresh_
const std::vector< const TrajectorySeedCollection * > *	initialSeedCollectionVector_ = nullptr
const float	lowPtThresh_
edm::ESHandle< MagneticField >	magField_
const float	nSigmasDeltaZ1_
const float	phiMax2B_
const float	phiMax2F_
const float	phiMin2B_
const float	phiMin2F_
PixelHitMatcher	positronMatcher_
std::unique_ptr< PropagatorWithMaterial >	propagator_
const float	sizeWindowENeg_
edm::ESHandle< TrackerGeometry >	trackerGeometry_
const bool	useRecoVertex_
const edm::EDGetTokenT< std::vector< reco::Vertex > >	verticesTag_

Detailed Description

Class to generate the trajectory seed from two hits in the pixel detector which have been found compatible with an ECAL cluster.

Author

U.Berthon, C.Charlot, LLR Palaiseau

Version

1st Version May 30, 2006

Description: Top algorithm producing ElectronSeeds, ported from ORCA

Implementation: future redesign...

Definition at line 41 of file ElectronSeedGenerator.h.

Member Typedef Documentation

ConstRecHitPointer

typedef TransientTrackingRecHit::ConstRecHitPointer ElectronSeedGenerator::ConstRecHitPointer

Definition at line 49 of file ElectronSeedGenerator.h.

PRecHitContainer

typedef edm::OwnVector<TrackingRecHit> ElectronSeedGenerator::PRecHitContainer

Definition at line 48 of file ElectronSeedGenerator.h.

RecHitContainer

typedef TransientTrackingRecHit::RecHitContainer ElectronSeedGenerator::RecHitContainer

Definition at line 51 of file ElectronSeedGenerator.h.

RecHitPointer

typedef TransientTrackingRecHit::RecHitPointer ElectronSeedGenerator::RecHitPointer

Definition at line 50 of file ElectronSeedGenerator.h.

Constructor & Destructor Documentation

ElectronSeedGenerator()

ElectronSeedGenerator::ElectronSeedGenerator	(	const edm::ParameterSet &	pset,
		const Tokens &	ts
	)

Definition at line 131 of file ElectronSeedGenerator.cc.

    : dynamicPhiRoad_(pset.getParameter<bool>("dynamicPhiRoad")),
      verticesTag_(ts.token_vtx),
      beamSpotTag_(ts.token_bs),
      lowPtThresh_(pset.getParameter<double>("LowPtThreshold")),
      highPtThresh_(pset.getParameter<double>("HighPtThreshold")),
      nSigmasDeltaZ1_(pset.getParameter<double>("nSigmasDeltaZ1")),
      deltaZ1WithVertex_(pset.getParameter<double>("deltaZ1WithVertex")),
      sizeWindowENeg_(pset.getParameter<double>("SizeWindowENeg")),
      deltaPhi1Low_(pset.getParameter<double>("DeltaPhi1Low")),
      deltaPhi1High_(pset.getParameter<double>("DeltaPhi1High")),
      // so that deltaPhi1 = dPhi1Coef1_ + dPhi1Coef2_/clusterEnergyT
      dPhi1Coef2_(dynamicPhiRoad_ ? (deltaPhi1Low_ - deltaPhi1High_) / (1. / lowPtThresh_ - 1. / highPtThresh_) : 0.),
      dPhi1Coef1_(dynamicPhiRoad_ ? deltaPhi1Low_ - dPhi1Coef2_ / lowPtThresh_ : 0.),
      propagator_(nullptr),
      // use of reco vertex
      useRecoVertex_(pset.getParameter<bool>("useRecoVertex")),
      // new B/F configurables
      deltaPhi2B_(pset.getParameter<double>("DeltaPhi2B")),
      deltaPhi2F_(pset.getParameter<double>("DeltaPhi2F")),
      phiMin2B_(-pset.getParameter<double>("PhiMax2B")),
      phiMin2F_(-pset.getParameter<double>("PhiMax2F")),
      phiMax2B_(pset.getParameter<double>("PhiMax2B")),
      phiMax2F_(pset.getParameter<double>("PhiMax2F")),
      electronMatcher_(pset.getParameter<double>("ePhiMin1"),
                       pset.getParameter<double>("ePhiMax1"),
                       phiMin2B_,
                       phiMax2B_,
                       phiMin2F_,
                       phiMax2F_,
                       -pset.getParameter<double>("z2MaxB"),
                       pset.getParameter<double>("z2MaxB"),
                       -pset.getParameter<double>("r2MaxF"),
                       pset.getParameter<double>("r2MaxF"),
                       -pset.getParameter<double>("rMaxI"),
                       pset.getParameter<double>("rMaxI"),
                       useRecoVertex_),
      positronMatcher_(pset.getParameter<double>("pPhiMin1"),
                       pset.getParameter<double>("pPhiMax1"),
                       phiMin2B_,
                       phiMax2B_,
                       phiMin2F_,
                       phiMax2F_,
                       -pset.getParameter<double>("z2MaxB"),
                       pset.getParameter<double>("z2MaxB"),
                       -pset.getParameter<double>("r2MaxF"),
                       pset.getParameter<double>("r2MaxF"),
                       -pset.getParameter<double>("rMaxI"),
                       pset.getParameter<double>("rMaxI"),
                       useRecoVertex_) {}

Member Function Documentation

run()

void ElectronSeedGenerator::run	(	edm::Event &	e,
		const edm::EventSetup &	setup,
		const reco::SuperClusterRefVector &	sclRefs,
		const std::vector< const TrajectorySeedCollection * > &	seedsV,
		reco::ElectronSeedCollection &	out
	)

Definition at line 205 of file ElectronSeedGenerator.cc.

                                                                 {
  initialSeedCollectionVector_ = &seedsV;
 
  // get the beamspot from the Event:
  auto const &beamSpot = e.get(beamSpotTag_);
 
  // if required get the vertices
  std::vector<reco::Vertex> const *vertices = nullptr;
  if (useRecoVertex_)
    vertices = &e.get(verticesTag_);
 
  for (unsigned int i = 0; i < sclRefs.size(); ++i) {
    // Find the seeds
    LogDebug("ElectronSeedGenerator") << "new cluster, calling seedsFromThisCluster";
    seedsFromThisCluster(sclRefs[i], beamSpot, vertices, out);
  }
 
  LogDebug("ElectronSeedGenerator") << ": For event " << e.id();
  LogDebug("ElectronSeedGenerator") << "Nr of superclusters after filter: " << sclRefs.size()
                                    << ", no. of ElectronSeeds found  = " << out.size();
}

References pwdgSkimBPark_cfi::beamSpot, beamSpotTag_, MillePedeFileConverter_cfg::e, mps_fire::i, initialSeedCollectionVector_, LogDebug, MillePedeFileConverter_cfg::out, seedsFromThisCluster(), edm::RefVector< C, T, F >::size(), useRecoVertex_, pwdgSkimBPark_cfi::vertices, and verticesTag_.

seedsFromThisCluster()

void ElectronSeedGenerator::seedsFromThisCluster ( edm::Ref< reco::SuperClusterCollection >  seedCluster, reco::BeamSpot const &  beamSpot, std::vector< reco::Vertex > const *  vertices, reco::ElectronSeedCollection &  out  )

private

Definition at line 231 of file ElectronSeedGenerator.cc.

                                                                                  {
  float clusterEnergy = seedCluster->energy();
  GlobalPoint clusterPos(seedCluster->position().x(), seedCluster->position().y(), seedCluster->position().z());
  reco::ElectronSeed::CaloClusterRef caloCluster(seedCluster);
 
  if (dynamicPhiRoad_) {
    float clusterEnergyT = clusterEnergy / cosh(EleRelPoint(clusterPos, beamSpot.position()).eta());
 
    float deltaPhi1;
    if (clusterEnergyT < lowPtThresh_) {
      deltaPhi1 = deltaPhi1Low_;
    } else if (clusterEnergyT > highPtThresh_) {
      deltaPhi1 = deltaPhi1High_;
    } else {
      deltaPhi1 = dPhi1Coef1_ + dPhi1Coef2_ / clusterEnergyT;
    }
 
    float ephimin1 = -deltaPhi1 * sizeWindowENeg_;
    float ephimax1 = deltaPhi1 * (1. - sizeWindowENeg_);
    float pphimin1 = -deltaPhi1 * (1. - sizeWindowENeg_);
    float pphimax1 = deltaPhi1 * sizeWindowENeg_;
 
    float phimin2B = -deltaPhi2B_ / 2.;
    float phimax2B = deltaPhi2B_ / 2.;
    float phimin2F = -deltaPhi2F_ / 2.;
    float phimax2F = deltaPhi2F_ / 2.;
 
    electronMatcher_.set1stLayer(ephimin1, ephimax1);
    positronMatcher_.set1stLayer(pphimin1, pphimax1);
    electronMatcher_.set2ndLayer(phimin2B, phimax2B, phimin2F, phimax2F);
    positronMatcher_.set2ndLayer(phimin2B, phimax2B, phimin2F, phimax2F);
  }
 
  if (!useRecoVertex_)  // here use the beam spot position
  {
    double sigmaZ = beamSpot.sigmaZ();
    double sigmaZ0Error = beamSpot.sigmaZ0Error();
    double sq = sqrt(sigmaZ * sigmaZ + sigmaZ0Error * sigmaZ0Error);
    double myZmin1 = beamSpot.position().z() - nSigmasDeltaZ1_ * sq;
    double myZmax1 = beamSpot.position().z() + nSigmasDeltaZ1_ * sq;
 
    GlobalPoint vertexPos;
    ele_convert(beamSpot.position(), vertexPos);
 
    electronMatcher_.set1stLayerZRange(myZmin1, myZmax1);
    positronMatcher_.set1stLayerZRange(myZmin1, myZmax1);
 
    // try electron
    auto elePixelSeeds = electronMatcher_(*initialSeedCollectionVector_, clusterPos, vertexPos, clusterEnergy, -1.);
    seedsFromTrajectorySeeds(elePixelSeeds, caloCluster, out, false);
    // try positron
    auto posPixelSeeds = positronMatcher_(*initialSeedCollectionVector_, clusterPos, vertexPos, clusterEnergy, 1.);
    seedsFromTrajectorySeeds(posPixelSeeds, caloCluster, out, true);
 
  } else if (vertices)  // here we use the reco vertices
  {
    for (auto const &vertex : *vertices) {
      GlobalPoint vertexPos(vertex.position().x(), vertex.position().y(), vertex.position().z());
      double myZmin1, myZmax1;
      if (vertexPos.z() == beamSpot.position().z()) {  // in case vetex not found
        double sigmaZ = beamSpot.sigmaZ();
        double sigmaZ0Error = beamSpot.sigmaZ0Error();
        double sq = sqrt(sigmaZ * sigmaZ + sigmaZ0Error * sigmaZ0Error);
        myZmin1 = beamSpot.position().z() - nSigmasDeltaZ1_ * sq;
        myZmax1 = beamSpot.position().z() + nSigmasDeltaZ1_ * sq;
      } else {  // a vertex has been recoed
        myZmin1 = vertex.position().z() - deltaZ1WithVertex_;
        myZmax1 = vertex.position().z() + deltaZ1WithVertex_;
      }
 
      electronMatcher_.set1stLayerZRange(myZmin1, myZmax1);
      positronMatcher_.set1stLayerZRange(myZmin1, myZmax1);
 
      // try electron
      auto elePixelSeeds = electronMatcher_(*initialSeedCollectionVector_, clusterPos, vertexPos, clusterEnergy, -1.);
      seedsFromTrajectorySeeds(elePixelSeeds, caloCluster, out, false);
      // try positron
      auto posPixelSeeds = positronMatcher_(*initialSeedCollectionVector_, clusterPos, vertexPos, clusterEnergy, 1.);
      seedsFromTrajectorySeeds(posPixelSeeds, caloCluster, out, true);
    }
  }
}

References pwdgSkimBPark_cfi::beamSpot, deltaPhi1High_, deltaPhi1Low_, deltaPhi2B_, deltaPhi2F_, deltaZ1WithVertex_, dPhi1Coef1_, dPhi1Coef2_, dynamicPhiRoad_, ele_convert(), electronMatcher_, PVValHelper::eta, highPtThresh_, initialSeedCollectionVector_, lowPtThresh_, nSigmasDeltaZ1_, MillePedeFileConverter_cfg::out, positronMatcher_, PixelHitMatcher::set1stLayer(), PixelHitMatcher::set1stLayerZRange(), PixelHitMatcher::set2ndLayer(), fftjetvertexadder_cfi::sigmaZ, sizeWindowENeg_, mathSSE::sqrt(), useRecoVertex_, bphysicsOniaDQM_cfi::vertex, and pwdgSkimBPark_cfi::vertices.

Referenced by run().

setupES()

void ElectronSeedGenerator::setupES

(

const edm::EventSetup &

setup

)

Definition at line 182 of file ElectronSeedGenerator.cc.

                                                              {
  // get records if necessary (called once per event)
  bool tochange = false;
 
  if (cacheIDMagField_ != setup.get<IdealMagneticFieldRecord>().cacheIdentifier()) {
    setup.get<IdealMagneticFieldRecord>().get(magField_);
    cacheIDMagField_ = setup.get<IdealMagneticFieldRecord>().cacheIdentifier();
    propagator_ = std::make_unique<PropagatorWithMaterial>(alongMomentum, .000511, &(*magField_));
    tochange = true;
  }
 
  if (cacheIDTrkGeom_ != setup.get<TrackerDigiGeometryRecord>().cacheIdentifier()) {
    cacheIDTrkGeom_ = setup.get<TrackerDigiGeometryRecord>().cacheIdentifier();
    setup.get<TrackerDigiGeometryRecord>().get(trackerGeometry_);
    tochange = true;  //FIXME
  }
 
  if (tochange) {
    electronMatcher_.setES(magField_.product(), trackerGeometry_.product());
    positronMatcher_.setES(magField_.product(), trackerGeometry_.product());
  }
}

References alongMomentum, edm::eventsetup::EventSetupRecord::cacheIdentifier(), cacheIDMagField_, cacheIDTrkGeom_, electronMatcher_, get, magField_, positronMatcher_, edm::ESHandle< T >::product(), propagator_, PixelHitMatcher::setES(), singleTopDQM_cfi::setup, and trackerGeometry_.

Member Data Documentation

beamSpotTag_

const edm::EDGetTokenT<reco::BeamSpot> ElectronSeedGenerator::beamSpotTag_

private

Definition at line 70 of file ElectronSeedGenerator.h.

Referenced by run().

cacheIDCkfComp_

unsigned long long ElectronSeedGenerator::cacheIDCkfComp_ = 0

private

Definition at line 93 of file ElectronSeedGenerator.h.

cacheIDMagField_

unsigned long long ElectronSeedGenerator::cacheIDMagField_ = 0

private

Definition at line 92 of file ElectronSeedGenerator.h.

Referenced by setupES().

cacheIDTrkGeom_

unsigned long long ElectronSeedGenerator::cacheIDTrkGeom_ = 0

private

Definition at line 94 of file ElectronSeedGenerator.h.

Referenced by setupES().

deltaPhi1High_

const float ElectronSeedGenerator::deltaPhi1High_

private

Definition at line 79 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

deltaPhi1Low_

const float ElectronSeedGenerator::deltaPhi1Low_

private

Definition at line 78 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

deltaPhi2B_

const float ElectronSeedGenerator::deltaPhi2B_

private

Definition at line 98 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

deltaPhi2F_

const float ElectronSeedGenerator::deltaPhi2F_

private

Definition at line 99 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

deltaZ1WithVertex_

const float ElectronSeedGenerator::deltaZ1WithVertex_

private

Definition at line 75 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

dPhi1Coef1_

const double ElectronSeedGenerator::dPhi1Coef1_

private

Definition at line 83 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

dPhi1Coef2_

const double ElectronSeedGenerator::dPhi1Coef2_

private

Definition at line 82 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

dynamicPhiRoad_

const bool ElectronSeedGenerator::dynamicPhiRoad_

private

Definition at line 68 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

electronMatcher_

PixelHitMatcher ElectronSeedGenerator::electronMatcher_

private

Definition at line 106 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster(), and setupES().

highPtThresh_

const float ElectronSeedGenerator::highPtThresh_

private

Definition at line 73 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

initialSeedCollectionVector_

const std::vector<const TrajectorySeedCollection*>* ElectronSeedGenerator::initialSeedCollectionVector_ = nullptr

private

Definition at line 85 of file ElectronSeedGenerator.h.

Referenced by run(), and seedsFromThisCluster().

lowPtThresh_

const float ElectronSeedGenerator::lowPtThresh_

private

Definition at line 72 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

magField_

edm::ESHandle<MagneticField> ElectronSeedGenerator::magField_

private

Definition at line 87 of file ElectronSeedGenerator.h.

Referenced by setupES().

nSigmasDeltaZ1_

const float ElectronSeedGenerator::nSigmasDeltaZ1_

private

Definition at line 74 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

phiMax2B_

const float ElectronSeedGenerator::phiMax2B_

private

Definition at line 103 of file ElectronSeedGenerator.h.

phiMax2F_

const float ElectronSeedGenerator::phiMax2F_

private

Definition at line 104 of file ElectronSeedGenerator.h.

phiMin2B_

const float ElectronSeedGenerator::phiMin2B_

private

Definition at line 101 of file ElectronSeedGenerator.h.

phiMin2F_

const float ElectronSeedGenerator::phiMin2F_

private

Definition at line 102 of file ElectronSeedGenerator.h.

positronMatcher_

PixelHitMatcher ElectronSeedGenerator::positronMatcher_

private

Definition at line 107 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster(), and setupES().

propagator_

std::unique_ptr<PropagatorWithMaterial> ElectronSeedGenerator::propagator_

private

Definition at line 89 of file ElectronSeedGenerator.h.

Referenced by setupES().

sizeWindowENeg_

const float ElectronSeedGenerator::sizeWindowENeg_

private

Definition at line 76 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

trackerGeometry_

edm::ESHandle<TrackerGeometry> ElectronSeedGenerator::trackerGeometry_

private

Definition at line 88 of file ElectronSeedGenerator.h.

Referenced by setupES().

useRecoVertex_

const bool ElectronSeedGenerator::useRecoVertex_

private

Definition at line 96 of file ElectronSeedGenerator.h.

Referenced by run(), and seedsFromThisCluster().

verticesTag_

const edm::EDGetTokenT<std::vector<reco::Vertex> > ElectronSeedGenerator::verticesTag_

private

Definition at line 69 of file ElectronSeedGenerator.h.

Referenced by run().