#include <PositionAtECalEntranceComputer.h>

Public Member Functions
void	beginEvent (const edm::EventSetup &)

reco::Candidate::Point	operator() (const reco::Candidate *particle, bool &success) const

	PositionAtECalEntranceComputer (edm::ConsumesCollector &&, bool isPhase2=false)

	PositionAtECalEntranceComputer (edm::ConsumesCollector &, bool isPhase2=false)

	~PositionAtECalEntranceComputer ()

Private Attributes
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	bField_esToken_

double	bField_z_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeo_esToken_

float	hgcalFace_z_

bool	isPhase2_

hgcal::RecHitTools	recHitTools_

Static Private Attributes
static constexpr float	ecalBarrelEndcapEtaBorder_ = 1.479

static constexpr float	hgcalHfEtaBorder_ = 3.0

Detailed Description

Extrapolate particle (charged or neutral) to ECAL entrance, in order to compute the distance of the tau to ECAL cracks and/or dead ECAL channels

Authors: Fabio Colombo, Christian Veelken

Definition at line 25 of file PositionAtECalEntranceComputer.h.

Constructor & Destructor Documentation

◆ PositionAtECalEntranceComputer() [1/2]

PositionAtECalEntranceComputer::PositionAtECalEntranceComputer	(	edm::ConsumesCollector &&	cc,
		bool	isPhase2 = `false`
	)

Definition at line 11 of file PositionAtECalEntranceComputer.cc.

     : bField_esToken_(cc.esConsumes<MagneticField, IdealMagneticFieldRecord>()),
       caloGeo_esToken_(cc.esConsumes<CaloGeometry, CaloGeometryRecord>()),
       bField_z_(-1.),
       isPhase2_(isPhase2) {}

◆ PositionAtECalEntranceComputer() [2/2]

PositionAtECalEntranceComputer::PositionAtECalEntranceComputer	(	edm::ConsumesCollector &	cc,
		bool	isPhase2 = `false`
	)

Definition at line 17 of file PositionAtECalEntranceComputer.cc.

     : bField_esToken_(cc.esConsumes<MagneticField, IdealMagneticFieldRecord>()),
       caloGeo_esToken_(cc.esConsumes<CaloGeometry, CaloGeometryRecord>()),
       bField_z_(-1.),
       isPhase2_(isPhase2) {}

◆ ~PositionAtECalEntranceComputer()

PositionAtECalEntranceComputer::~PositionAtECalEntranceComputer ( )

Definition at line 23 of file PositionAtECalEntranceComputer.cc.

23 {}

Member Function Documentation

◆ beginEvent()

void PositionAtECalEntranceComputer::beginEvent ( const edm::EventSetup & es )

Definition at line 25 of file PositionAtECalEntranceComputer.cc.

                                                                        {
   bField_z_ = es.getData(bField_esToken_).inTesla(GlobalPoint(0., 0., 0.)).z();
   if (isPhase2_) {
     recHitTools_.setGeometry(es.getData(caloGeo_esToken_));
     hgcalFace_z_ = recHitTools_.getPositionLayer(1).z();  // HGCal 1st layer
   }
 }

References bField_esToken_, bField_z_, caloGeo_esToken_, edm::EventSetup::getData(), hgcal::RecHitTools::getPositionLayer(), hgcalFace_z_, isPhase2_, recHitTools_, hgcal::RecHitTools::setGeometry(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by AntiElectronDeadECAL::beginEvent(), TauDiscriminationAgainstElectronMVA6< TauType, TauDiscriminator, ElectronType >::beginEvent(), and pat::PATTauProducer::produce().

◆ operator()()

reco::Candidate::Point PositionAtECalEntranceComputer::operator()	(	const reco::Candidate *	particle,
		bool &	success
	)		const

Definition at line 33 of file PositionAtECalEntranceComputer.cc.

                                                                                        {
   assert(bField_z_ != -1.);
   reco::Candidate::Point position;
   BaseParticlePropagator propagator = BaseParticlePropagator(
       RawParticle(particle->p4(),
                   math::XYZTLorentzVector(particle->vertex().x(), particle->vertex().y(), particle->vertex().z(), 0.),
                   particle->charge()),
       0.,
       0.,
       bField_z_);
   if (!isPhase2_ || std::abs(particle->eta()) < ecalBarrelEndcapEtaBorder_) {  // ECal
     propagator.propagateToEcalEntrance(false);
   } else {  // HGCal
     if (std::abs(particle->vertex().z()) >= hgcalFace_z_) {
       success = false;
       return position;
     }
     propagator.setPropagationConditions(152.6, hgcalFace_z_, false);
     propagator.propagate();
   }
   if (propagator.getSuccess() != 0) {
     position = propagator.particle().vertex().Vect();
     success = (std::abs(position.eta()) <= hgcalHfEtaBorder_);
   } else {
     success = false;
   }
   return position;
 }