l1tpf_impl::PUAlgoBase Class Reference

#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

VertexAlgo

enum l1tpf_impl::PUAlgoBase::VertexAlgo

strong

global operations

Enumerator
Old
TP
External

Definition at line 15 of file PUAlgoBase.h.

{ Old, TP, External };

Constructor & Destructor Documentation

PUAlgoBase()

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::~PUAlgoBase ( )

virtual

Definition at line 13 of file PUAlgoBase.cc.

{}

Member Function Documentation

doPUGlobals()

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.

doVertexing()

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, ReggeGribovPartonMC_AdvancedParameters_cfi::bmax, central, External, createfilelist::int, SiStripPI::max, SiStripPI::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));
    }
  }
}

puGlobalNames()

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

runChargedPV()

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

virtual

Definition at line 15 of file PUAlgoBase.cc.

References funct::abs(), Reference_intrackfit_cff::barrel, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, vtxAdaptiveCut_, vtxRes_, and HLTMuonOfflineAnalyzer_cfi::z0.

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

runNeutralsPU()

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

pure virtual

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

Member Data Documentation

debug_

int l1tpf_impl::PUAlgoBase::debug_

protected

Definition at line 27 of file PUAlgoBase.h.

Referenced by l1tpf_impl::PuppiAlgo::computePuppiMedRMS(), l1tpf_impl::LinearizedPuppiAlgo::computePuppiWeights(), l1tpf_impl::PuppiAlgo::computePuppiWeights(), and l1tpf_impl::PuppiAlgo::PuppiAlgo().

etaCharged_

float l1tpf_impl::PUAlgoBase::etaCharged_

protected

Definition at line 28 of file PUAlgoBase.h.

Referenced by l1tpf_impl::PuppiAlgo::computePuppiMedRMS(), l1tpf_impl::LinearizedPuppiAlgo::computePuppiWeights(), and l1tpf_impl::PuppiAlgo::computePuppiWeights().

vtxAdaptiveCut_

bool l1tpf_impl::PUAlgoBase::vtxAdaptiveCut_

protected

Definition at line 29 of file PUAlgoBase.h.

Referenced by doVertexing(), and runChargedPV().

vtxRes_

float l1tpf_impl::PUAlgoBase::vtxRes_

protected

Definition at line 28 of file PUAlgoBase.h.

Referenced by doVertexing(), and runChargedPV().