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 &, edm::ConsumesCollector &&)
void	run (edm::Event &, 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_
const float	deltaPhi1High_
const float	deltaPhi1Low_
const float	deltaPhi2B_
const float	deltaPhi2F_
const float	deltaZ1WithVertex_
const double	dPhi1Coef1_
const double	dPhi1Coef2_
const bool	dynamicPhiRoad_
const float	highPtThresh_
const std::vector< const TrajectorySeedCollection * > *	initialSeedCollectionVector_ = nullptr
const float	lowPtThresh_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_
edm::ESWatcher< IdealMagneticFieldRecord >	magneticFieldWatcher_
PixelHitMatcher	matcher_
const float	nSigmasDeltaZ1_
const float	phiMax2B_
const float	phiMax2F_
const float	phiMin2B_
const float	phiMin2F_
const float	sizeWindowENeg_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_
edm::ESWatcher< TrackerDigiGeometryRecord >	trackerGeometryWatcher_
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 46 of file ElectronSeedGenerator.h.

Member Typedef Documentation

ConstRecHitPointer

typedef TransientTrackingRecHit::ConstRecHitPointer ElectronSeedGenerator::ConstRecHitPointer

Definition at line 54 of file ElectronSeedGenerator.h.

PRecHitContainer

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

Definition at line 53 of file ElectronSeedGenerator.h.

RecHitContainer

typedef TransientTrackingRecHit::RecHitContainer ElectronSeedGenerator::RecHitContainer

Definition at line 56 of file ElectronSeedGenerator.h.

RecHitPointer

typedef TransientTrackingRecHit::RecHitPointer ElectronSeedGenerator::RecHitPointer

Definition at line 55 of file ElectronSeedGenerator.h.

Constructor & Destructor Documentation

ElectronSeedGenerator()

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

Definition at line 128 of file ElectronSeedGenerator.cc.

    : dynamicPhiRoad_(pset.getParameter<bool>("dynamicPhiRoad")),
      verticesTag_(ts.token_vtx),
      beamSpotTag_(ts.token_bs),
      magFieldToken_{cc.esConsumes()},
      trackerGeometryToken_{cc.esConsumes()},
      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.),
      // 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")),
      matcher_(pset.getParameter<double>("ePhiMin1"),
               pset.getParameter<double>("ePhiMax1"),
               phiMin2B_,
               phiMax2B_,
               phiMin2F_,
               phiMax2F_,
               pset.getParameter<double>("z2MaxB"),
               pset.getParameter<double>("r2MaxF"),
               pset.getParameter<double>("rMaxI"),
               useRecoVertex_) {}

Member Function Documentation

run()

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

Definition at line 172 of file ElectronSeedGenerator.cc.

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

                                                                 {
  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();
}

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 197 of file ElectronSeedGenerator.cc.

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

Referenced by run().

                                                                                  {
  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;
    }

    matcher_.set1stLayer(-deltaPhi1 * sizeWindowENeg_, deltaPhi1 * (1. - sizeWindowENeg_));
    matcher_.set2ndLayer(-deltaPhi2B_ / 2., deltaPhi2B_ / 2., -deltaPhi2F_ / 2., deltaPhi2F_ / 2.);
  }

  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);

    matcher_.set1stLayerZRange(myZmin1, myZmax1);

    // try electron
    auto elePixelSeeds = matcher_(*initialSeedCollectionVector_, clusterPos, vertexPos, clusterEnergy, -1.);
    seedsFromTrajectorySeeds(elePixelSeeds, caloCluster, out, false);
    // try positron
    auto posPixelSeeds = matcher_(*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_;
      }

      matcher_.set1stLayerZRange(myZmin1, myZmax1);

      // try electron
      auto elePixelSeeds = matcher_(*initialSeedCollectionVector_, clusterPos, vertexPos, clusterEnergy, -1.);
      seedsFromTrajectorySeeds(elePixelSeeds, caloCluster, out, false);
      // try positron
      auto posPixelSeeds = matcher_(*initialSeedCollectionVector_, clusterPos, vertexPos, clusterEnergy, 1.);
      seedsFromTrajectorySeeds(posPixelSeeds, caloCluster, out, true);
    }
  }
}

setupES()

void ElectronSeedGenerator::setupES

(

const edm::EventSetup &

setup

)

Definition at line 166 of file ElectronSeedGenerator.cc.

References edm::ESWatcher< T >::check(), magFieldToken_, magneticFieldWatcher_, matcher_, PixelHitMatcher::setES(), singleTopDQM_cfi::setup, trackerGeometryToken_, and trackerGeometryWatcher_.

                                                              {
  if (magneticFieldWatcher_.check(setup) || trackerGeometryWatcher_.check(setup)) {
    matcher_.setES(setup.getData(magFieldToken_), setup.getData(trackerGeometryToken_));
  }
}

Member Data Documentation

beamSpotTag_

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

private

Definition at line 75 of file ElectronSeedGenerator.h.

Referenced by run().

deltaPhi1High_

const float ElectronSeedGenerator::deltaPhi1High_

private

Definition at line 88 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

deltaPhi1Low_

const float ElectronSeedGenerator::deltaPhi1Low_

private

Definition at line 87 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 84 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

dPhi1Coef1_

const double ElectronSeedGenerator::dPhi1Coef1_

private

Definition at line 92 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

dPhi1Coef2_

const double ElectronSeedGenerator::dPhi1Coef2_

private

Definition at line 91 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

dynamicPhiRoad_

const bool ElectronSeedGenerator::dynamicPhiRoad_

private

Definition at line 72 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

highPtThresh_

const float ElectronSeedGenerator::highPtThresh_

private

Definition at line 82 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

initialSeedCollectionVector_

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

private

Definition at line 94 of file ElectronSeedGenerator.h.

Referenced by run(), and seedsFromThisCluster().

lowPtThresh_

const float ElectronSeedGenerator::lowPtThresh_

private

Definition at line 81 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

magFieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> ElectronSeedGenerator::magFieldToken_

private

Definition at line 76 of file ElectronSeedGenerator.h.

Referenced by setupES().

magneticFieldWatcher_

edm::ESWatcher<IdealMagneticFieldRecord> ElectronSeedGenerator::magneticFieldWatcher_

private

Definition at line 78 of file ElectronSeedGenerator.h.

Referenced by setupES().

matcher_

PixelHitMatcher ElectronSeedGenerator::matcher_

private

Definition at line 106 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster(), and setupES().

nSigmasDeltaZ1_

const float ElectronSeedGenerator::nSigmasDeltaZ1_

private

Definition at line 83 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.

sizeWindowENeg_

const float ElectronSeedGenerator::sizeWindowENeg_

private

Definition at line 85 of file ElectronSeedGenerator.h.

Referenced by seedsFromThisCluster().

trackerGeometryToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> ElectronSeedGenerator::trackerGeometryToken_

private

Definition at line 77 of file ElectronSeedGenerator.h.

Referenced by setupES().

trackerGeometryWatcher_

edm::ESWatcher<TrackerDigiGeometryRecord> ElectronSeedGenerator::trackerGeometryWatcher_

private

Definition at line 79 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 74 of file ElectronSeedGenerator.h.

Referenced by run().