#include <PUAlgoBase.h>

Inheritance diagram for l1tpf_impl::PUAlgoBase:

Public Types
enum	VertexAlgo { VertexAlgo::Old, VertexAlgo::TP, VertexAlgo::External }
	global operations More...

Public Member Functions
virtual void	doPUGlobals (const std::vector< Region > &rs, float npu, std::vector< float > &globals) const =0

virtual void	doVertexing (std::vector< Region > &rs, VertexAlgo algo, float &vz) const

	PUAlgoBase (const edm::ParameterSet &)

virtual const std::vector < std::string > &	puGlobalNames () const

virtual void	runChargedPV (Region &r, float z0) const

virtual void	runNeutralsPU (Region &r, float npu, const std::vector< float > &globals) const =0

virtual	~PUAlgoBase ()

Protected Attributes
int	debug_

float	etaCharged_

bool	vtxAdaptiveCut_

float	vtxRes_

Detailed Description

Definition at line 9 of file PUAlgoBase.h.

Member Enumeration Documentation

enum l1tpf_impl::PUAlgoBase::VertexAlgo

strong

global operations

Enumerator
Old
TP
External

Definition at line 15 of file PUAlgoBase.h.

15 { Old, TP, External };

DigiToRawDM_cff.TP

string TP

Definition: DigiToRawDM_cff.py:51

Constructor & Destructor Documentation

PUAlgoBase::PUAlgoBase ( const edm::ParameterSet & iConfig )

Definition at line 7 of file PUAlgoBase.cc.

     : debug_(iConfig.getUntrackedParameter<int>("debug", 0)),
       etaCharged_(iConfig.getParameter<double>("etaCharged")),
       vtxRes_(iConfig.getParameter<double>("vtxRes")),
       vtxAdaptiveCut_(iConfig.getParameter<bool>("vtxAdaptiveCut")) {}

PUAlgoBase::~PUAlgoBase ( )

virtual

Definition at line 13 of file PUAlgoBase.cc.

13 {}

Member Function Documentation

virtual void l1tpf_impl::PUAlgoBase::doPUGlobals	(	const std::vector< Region > &	rs,
		float	npu,
		std::vector< float > &	globals
	)		const

pure virtual

Implemented in l1tpf_impl::LinearizedPuppiAlgo, and l1tpf_impl::PuppiAlgo.

void PUAlgoBase::doVertexing	(	std::vector< Region > &	rs,
		VertexAlgo	algo,
		float &	vz
	)		const

virtual

Definition at line 28 of file PUAlgoBase.cc.

References funct::abs(), b, central, External, l1tpf_impl::PropagatedTrack::fromPV, l1tpf_impl::InputTrack::hwVtxEta, l1tpf_impl::InputTrack::hwZ0, SiStripPI::max, min(), Old, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, TP, vtxAdaptiveCut_, and vtxRes_.

                                                                                       {
   int lNBins = int(40. / vtxRes_);
   if (algo == VertexAlgo::TP)
     lNBins *= 3;
   std::unique_ptr<TH1F> h_dz(new TH1F("h_dz", "h_dz", lNBins, -20, 20));
   if (algo != VertexAlgo::External) {
     for (const Region &r : rs) {
       for (const PropagatedTrack &p : r.track) {
         if (rs.size() > 1) {
           if (!r.fiducialLocal(p.floatVtxEta(), p.floatVtxPhi()))
             continue;  // skip duplicates
         }
         h_dz->Fill(p.floatDZ(), std::min(p.floatPt(), 50.f));
       }
     }
   }
   switch (algo) {
     case VertexAlgo::External:
       break;
     case VertexAlgo::Old: {
       int imaxbin = h_dz->GetMaximumBin();
       pvdz = h_dz->GetXaxis()->GetBinCenter(imaxbin);
     }; break;
     case VertexAlgo::TP: {
       float max = 0;
       int bmax = -1;
       for (int b = 1; b <= lNBins; ++b) {
         float sum3 = h_dz->GetBinContent(b) + h_dz->GetBinContent(b + 1) + h_dz->GetBinContent(b - 1);
         if (bmax == -1 || sum3 > max) {
           max = sum3;
           bmax = b;
         }
       }
       pvdz = h_dz->GetXaxis()->GetBinCenter(bmax);
     }; break;
   }
   int16_t iZ0 = round(pvdz * InputTrack::Z0_SCALE);
   int16_t iDZ = round(1.5 * vtxRes_ * InputTrack::Z0_SCALE);
   int16_t iDZ2 = vtxAdaptiveCut_ ? round(4.0 * vtxRes_ * InputTrack::Z0_SCALE) : iDZ;
   for (Region &r : rs) {
     for (PropagatedTrack &p : r.track) {
       bool central = std::abs(p.hwVtxEta) < InputTrack::VTX_ETA_1p3;
       if (r.relativeCoordinates)
         central =
             (std::abs(r.globalAbsEta(p.floatVtxEta())) < 1.3);  // FIXME could make a better integer implementation
       p.fromPV = (std::abs(p.hwZ0 - iZ0) < (central ? iDZ : iDZ2));
     }
   }
 }

const std::vector< std::string > & PUAlgoBase::puGlobalNames ( ) const

virtual

Reimplemented in l1tpf_impl::LinearizedPuppiAlgo, and l1tpf_impl::PuppiAlgo.

Definition at line 78 of file PUAlgoBase.cc.

                                                             {
   static const std::vector<std::string> empty_;
   return empty_;
 }

void PUAlgoBase::runChargedPV	(	Region &	r,
		float	z0
	)		const

virtual

Definition at line 15 of file PUAlgoBase.cc.

References funct::abs(), Reference_intrackfit_cff::barrel, l1tpf_impl::PFParticle::chargedPV, l1tpf_impl::Region::globalAbsEta(), l1tpf_impl::PFParticle::hwId, l1tpf_impl::InputTrack::hwVtxEta, l1tpf_impl::InputTrack::hwZ0, AlCaHLTBitMon_ParallelJobs::p, l1tpf_impl::Region::pf, l1tpf_impl::Region::relativeCoordinates, l1tpf_impl::PFParticle::track, vtxAdaptiveCut_, and vtxRes_.

                                                        {
   int16_t iZ0 = round(z0 * InputTrack::Z0_SCALE);
   int16_t iDZ = round(1.5 * vtxRes_ * InputTrack::Z0_SCALE);
   int16_t iDZ2 = vtxAdaptiveCut_ ? round(4.0 * vtxRes_ * InputTrack::Z0_SCALE) : iDZ;
   for (PFParticle &p : r.pf) {
     bool barrel = std::abs(p.track.hwVtxEta) < InputTrack::VTX_ETA_1p3;
     if (r.relativeCoordinates)
       barrel =
           (std::abs(r.globalAbsEta(p.track.floatVtxEta())) < 1.3);  // FIXME could make a better integer implementation
     p.chargedPV = (p.hwId <= 1 && std::abs(p.track.hwZ0 - iZ0) < (barrel ? iDZ : iDZ2));
   }
 }