FastPrimaryVertexWithWeightsProducer Class Reference

#include <RecoBTag/FastPrimaryVertexWithWeightsProducer/src/FastPrimaryVertexWithWeightsProducer.cc>

Inheritance diagram for FastPrimaryVertexWithWeightsProducer:

Public Member Functions
	FastPrimaryVertexWithWeightsProducer (const edm::ParameterSet &)
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken
edm::EDGetTokenT< SiPixelClusterCollectionNew >	clustersToken
edm::EDGetTokenT< edm::View< reco::Jet > >	jetsToken
const bool	m_barrel
const edm::InputTag	m_clusters
const double	m_EC_weight
const bool	m_endCap
const double	m_maxDeltaPhi
const double	m_maxDeltaPhi_EC
const double	m_maxJetEta
const double	m_maxJetEta_EC
const double	m_maxSizeX
const double	m_maxSizeY_q
const double	m_maxZ
const double	m_minJetEta_EC
const double	m_minJetPt
const double	m_minSizeY_q
const int	m_njets
const double	m_peakSizeY_q
const double	m_PixelCellHeightOverWidth
std::string	m_pixelCPE
const bool	m_ptWeighting
const double	m_ptWeighting_offset
const double	m_ptWeighting_slope
const double	m_weight_charge_down
const double	m_weight_charge_peak
const double	m_weight_charge_up
const double	m_weight_dPhi
const double	m_weight_dPhi_EC
const double	m_weight_rho_up
const double	m_weight_SizeX1
const double	m_weightCut_step2
const double	m_weightCut_step3
const double	m_zClusterSearchArea_step2
const double	m_zClusterSearchArea_step3
const double	m_zClusterWidth_step1
const double	m_zClusterWidth_step2
const double	m_zClusterWidth_step3

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: The FastPrimaryVertex is an algorithm used to find the primary vertex of an event . It takes the pixel clusters compabible with a jet and project it to the beamSpot with the eta-angle of the jet. The peak on the z-projected clusters distribution is our FastPrimaryVertex.

Implementation: [Notes on implementation]

Definition at line 56 of file FastPrimaryVertexWithWeightsProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 125 of file FastPrimaryVertexWithWeightsProducer.cc.

References beamSpotToken, clustersToken, edm::ParameterSet::getParameter(), and jetsToken.

                                                                                                        :
  m_maxZ(iConfig.getParameter<double>("maxZ")),
  m_pixelCPE(iConfig.getParameter<std::string>("pixelCPE")),

  m_njets(iConfig.getParameter<int>("njets")),
  m_maxJetEta(iConfig.getParameter<double>("maxJetEta")),
  m_minJetPt(iConfig.getParameter<double>("minJetPt")),

  m_barrel(iConfig.getParameter<bool>("barrel")),
  m_maxSizeX(iConfig.getParameter<double>("maxSizeX")),
  m_maxDeltaPhi(iConfig.getParameter<double>("maxDeltaPhi")),
  m_weight_charge_down(iConfig.getParameter<double>("weight_charge_down")),
  m_weight_charge_up(iConfig.getParameter<double>("weight_charge_up")),
  m_PixelCellHeightOverWidth(iConfig.getParameter<double>("PixelCellHeightOverWidth")),
  m_minSizeY_q(iConfig.getParameter<double>("minSizeY_q")),
  m_maxSizeY_q(iConfig.getParameter<double>("maxSizeY_q")),
  
  m_weight_dPhi(iConfig.getParameter<double>("weight_dPhi")),
  m_weight_SizeX1(iConfig.getParameter<double>("weight_SizeX1")),
  m_weight_rho_up(iConfig.getParameter<double>("weight_rho_up")),
  m_weight_charge_peak(iConfig.getParameter<double>("weight_charge_peak")),
  m_peakSizeY_q(iConfig.getParameter<double>("peakSizeY_q")),

  m_endCap(iConfig.getParameter<bool>("endCap")),
  m_minJetEta_EC(iConfig.getParameter<double>("minJetEta_EC")),
  m_maxJetEta_EC(iConfig.getParameter<double>("maxJetEta_EC")),
  m_maxDeltaPhi_EC(iConfig.getParameter<double>("maxDeltaPhi_EC")),
  m_EC_weight(iConfig.getParameter<double>("EC_weight")),
  m_weight_dPhi_EC(iConfig.getParameter<double>("weight_dPhi_EC")),

  m_zClusterWidth_step1(iConfig.getParameter<double>("zClusterWidth_step1")),

  m_zClusterWidth_step2(iConfig.getParameter<double>("zClusterWidth_step2")),
  m_zClusterSearchArea_step2(iConfig.getParameter<double>("zClusterSearchArea_step2")),
  m_weightCut_step2(iConfig.getParameter<double>("weightCut_step2")),

  m_zClusterWidth_step3(iConfig.getParameter<double>("zClusterWidth_step3")),
  m_zClusterSearchArea_step3(iConfig.getParameter<double>("zClusterSearchArea_step3")),
  m_weightCut_step3(iConfig.getParameter<double>("weightCut_step3")),

  m_ptWeighting(iConfig.getParameter<bool>("ptWeighting")),
  m_ptWeighting_slope(iConfig.getParameter<double>("ptWeighting_slope")),
  m_ptWeighting_offset(iConfig.getParameter<double>("ptWeighting_offset"))
{

  clustersToken = consumes<SiPixelClusterCollectionNew>(iConfig.getParameter<edm::InputTag>("clusters"));
  beamSpotToken = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpot"));
  jetsToken = consumes<edm::View<reco::Jet> >(iConfig.getParameter<edm::InputTag>("jets"));

  produces<reco::VertexCollection>();
  produces<float>();

}

Member Function Documentation

void FastPrimaryVertexWithWeightsProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 180 of file FastPrimaryVertexWithWeightsProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag> ("clusters",edm::InputTag("hltSiPixelClusters"));
  desc.add<edm::InputTag> ("beamSpot",edm::InputTag("hltOnlineBeamSpot"));
  desc.add<edm::InputTag> ("jets",edm::InputTag("hltCaloJetL1FastJetCorrected"));
  desc.add<std::string>("pixelCPE","hltESPPixelCPEGeneric");
  desc.add<double>("maxZ",19.0);
  desc.add<int>("njets",999);
  desc.add<double>("maxJetEta",2.6);
  desc.add<double>("minJetPt",40.);
  desc.add<bool>("barrel",true);
  desc.add<double>("maxSizeX",2.1);
  desc.add<double>("maxDeltaPhi",0.21);
  desc.add<double>("PixelCellHeightOverWidth",1.8);
  desc.add<double>("weight_charge_down",11.*1000.);
  desc.add<double>("weight_charge_up",190.*1000.);
  desc.add<double>("maxSizeY_q",2.0);
  desc.add<double>("minSizeY_q",-0.6);
  desc.add<double>("weight_dPhi",0.13888888);
  desc.add<double>("weight_SizeX1",0.88);
  desc.add<double>("weight_rho_up",22.);
  desc.add<double>("weight_charge_peak",22.*1000.);
  desc.add<double>("peakSizeY_q",1.0);
  desc.add<bool>("endCap",true);
  desc.add<double>("minJetEta_EC",1.3);
  desc.add<double>("maxJetEta_EC",2.6);
  desc.add<double>("maxDeltaPhi_EC",0.14);
  desc.add<double>("EC_weight",0.008);
  desc.add<double>("weight_dPhi_EC",0.064516129);
  desc.add<double>("zClusterWidth_step1",2.0);
  desc.add<double>("zClusterWidth_step2",0.65);
  desc.add<double>("zClusterSearchArea_step2",3.0);
  desc.add<double>("weightCut_step2",0.05);
  desc.add<double>("zClusterWidth_step3",0.3);
  desc.add<double>("zClusterSearchArea_step3",0.55);
  desc.add<double>("weightCut_step3",0.1);
  desc.add<bool>("ptWeighting",false); // <---- newMethod 
  desc.add<double>("ptWeighting_slope",1/20.);
  desc.add<double>("ptWeighting_offset",-1);
  descriptions.add("fastPrimaryVertexWithWeightsProducer",desc);
}

void FastPrimaryVertexWithWeightsProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

if it is a cluster to project

Definition at line 224 of file FastPrimaryVertexWithWeightsProducer.cc.

References funct::abs(), ecalDrivenElectronSeedsParameters_cff::beamSpot, beamSpotToken, edm::View< T >::begin(), edmNew::DetSetVector< T >::begin(), ALCARECOTkAlJpsiMuMu_cff::charge, SiPixelCluster::charge(), clustersToken, DEFINE_FWK_MODULE, hiPixelPairStep_cff::deltaPhi, MillePedeFileConverter_cfg::e, edm::View< T >::end(), edmNew::DetSetVector< T >::end(), PVValHelper::eta, JetChargeProducer_cfi::exp, FindPeakFastPV(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, TrackerGeometry::idToDet(), TrackerGeometry::idToDetUnit(), fwrapper::jets, jetsToken, PixelClusterParameterEstimator::localParametersV(), m_barrel, m_EC_weight, m_endCap, m_maxDeltaPhi, m_maxDeltaPhi_EC, m_maxJetEta, m_maxJetEta_EC, m_maxSizeX, m_maxSizeY_q, m_maxZ, m_minJetEta_EC, m_minJetPt, m_minSizeY_q, m_njets, m_peakSizeY_q, m_PixelCellHeightOverWidth, m_pixelCPE, m_ptWeighting, m_ptWeighting_offset, m_ptWeighting_slope, m_weight_charge_down, m_weight_charge_peak, m_weight_charge_up, m_weight_dPhi, m_weight_dPhi_EC, m_weight_rho_up, m_weight_SizeX1, m_weightCut_step2, m_weightCut_step3, m_zClusterSearchArea_step2, m_zClusterSearchArea_step3, m_zClusterWidth_step1, m_zClusterWidth_step2, m_zClusterWidth_step3, eostools::move(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::phi(), HiRegitMuonDetachedTripletStep_cff::pixelClusters, GeomDet::position(), createTree::pp, edm::Handle< T >::product(), edm::ESHandle< T >::product(), EnergyCorrector::pt, edm::Event::put(), rho, objectSelection_cff::selectedJets, edmNew::DetSet< T >::size(), SiPixelCluster::sizeX(), SiPixelCluster::sizeY(), RecoTauValidation_cfi::sizeY, mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), mixOne_simraw_on_sim_cfi::tracker, findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), reco::BeamSpot::x(), reco::BeamSpot::x0(), PV3DBase< T, PVType, FrameType >::y(), reco::BeamSpot::y(), reco::BeamSpot::y0(), z, and PV3DBase< T, PVType, FrameType >::z().

{

   using namespace edm;
   using namespace reco;
   using namespace std;

   const float barrel_lenght=30; //half lenght of the pixel barrel 30 cm

   //get pixel cluster
   Handle<SiPixelClusterCollectionNew> cH;
   iEvent.getByToken(clustersToken,cH);
   const SiPixelClusterCollectionNew & pixelClusters = *cH.product();

   //get jets
   Handle<edm::View<reco::Jet> > jH;
   iEvent.getByToken(jetsToken,jH);
   const edm::View<reco::Jet> & jets = *jH.product();

   vector<const reco::Jet*> selectedJets;
   int countjet=0;   
   for(edm::View<reco::Jet>::const_iterator it = jets.begin() ; it != jets.end() && countjet<m_njets ; it++)
   {
   if( //select jets used in barrel or endcap pixel cluster projections
    ((it->pt() >= m_minJetPt) && std::abs(it->eta()) <= m_maxJetEta) || //barrel 
    ((it->pt() >= m_minJetPt) && std::abs(it->eta()) <= m_maxJetEta_EC && std::abs(it->eta()) >= m_minJetEta_EC) //endcap 
    )
    {
      selectedJets.push_back(&(*it));
      countjet++;
    } 
   }
  
   //get PixelClusterParameterEstimator
   edm::ESHandle<PixelClusterParameterEstimator> pe; 
   const PixelClusterParameterEstimator * pp ;
   iSetup.get<TkPixelCPERecord>().get(m_pixelCPE , pe );  
   pp = pe.product();

   //get beamSpot
   edm::Handle<BeamSpot> beamSpot;
   iEvent.getByToken(beamSpotToken,beamSpot);
 
   //get TrackerGeometry
   edm::ESHandle<TrackerGeometry> tracker;
   iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
   const TrackerGeometry * trackerGeometry = tracker.product();

// PART I: get z-projections with z-weights
   std::vector<float> zProjections;
   std::vector<float> zWeights;
   int jet_count=0;
   for(vector<const reco::Jet*>::iterator jit = selectedJets.begin() ; jit != selectedJets.end() ; jit++)
   {//loop on selected jets
     float px=(*jit)->px();
     float py=(*jit)->py();
     float pz=(*jit)->pz();
     float pt=(*jit)->pt();
     float eta=(*jit)->eta();
     float jetZOverRho = (*jit)->momentum().Z()/(*jit)->momentum().Rho();
     float pt_weight= pt*m_ptWeighting_slope +m_ptWeighting_offset;
     for(SiPixelClusterCollectionNew::const_iterator it = pixelClusters.begin() ; it != pixelClusters.end() ; it++) //Loop on pixel modules with clusters
     {//loop on pixel modules
        DetId id = it->detId();
        const edmNew::DetSet<SiPixelCluster> & detset  = (*it);
        Point3DBase<float, GlobalTag> modulepos=trackerGeometry->idToDet(id)->position();
        float zmodule = modulepos.z() - ((modulepos.x()-beamSpot->x0())*px+(modulepos.y()-beamSpot->y0())*py)/pt * pz/pt; 
        if ((std::abs(deltaPhi((*jit)->momentum().Phi(),modulepos.phi()))< m_maxDeltaPhi*2)&&(std::abs(zmodule)<(m_maxZ+barrel_lenght))){//if it is a compatible module
        for(size_t j = 0 ; j < detset.size() ; j ++) 
        {//loop on pixel clusters on this module
      const SiPixelCluster & aCluster =  detset[j];
          if(//it is a cluster to project
          (// barrel
          m_barrel &&
          std::abs(modulepos.z())<barrel_lenght &&
          pt>=m_minJetPt &&
          jet_count<m_njets &&
          std::abs(eta)<=m_maxJetEta &&
          aCluster.sizeX() <= m_maxSizeX &&
          aCluster.sizeY()  >= m_PixelCellHeightOverWidth*std::abs(jetZOverRho)+m_minSizeY_q &&
          aCluster.sizeY() <= m_PixelCellHeightOverWidth*std::abs(jetZOverRho)+m_maxSizeY_q
          )
          ||
          (// EC
          m_endCap &&
          std::abs(modulepos.z())>barrel_lenght &&
          pt > m_minJetPt &&
          jet_count<m_njets &&
          std::abs(eta) <= m_maxJetEta_EC &&
          std::abs(eta) >= m_minJetEta_EC &&
          aCluster.sizeX() <= m_maxSizeX 
          )   
          ){
            Point3DBase<float, GlobalTag> v = trackerGeometry->idToDet(id)->surface().toGlobal(pp->localParametersV( aCluster,( *trackerGeometry->idToDetUnit(id)))[0].first) ;
            GlobalPoint v_bs(v.x()-beamSpot->x0(),v.y()-beamSpot->y0(),v.z());
            if(//it pass DeltaPhi(Jet,Cluster) requirements
          (m_barrel && std::abs(modulepos.z())<barrel_lenght && std::abs(deltaPhi((*jit)->momentum().Phi(),v_bs.phi())) <= m_maxDeltaPhi ) || //barrel
          (m_endCap && std::abs(modulepos.z())>barrel_lenght && std::abs(deltaPhi((*jit)->momentum().Phi(),v_bs.phi())) <= m_maxDeltaPhi_EC )  //EC
            )
            {
              float z = v.z() - ((v.x()-beamSpot->x0())*px+(v.y()-beamSpot->y0())*py)/pt * pz/pt;   //calculate z-projection
              if(std::abs(z) < m_maxZ)
              {
            zProjections.push_back(z); //add z-projection in zProjections
        float weight=0;
        //calculate zWeight
            if(std::abs(modulepos.z())<barrel_lenght)
            { //barrel
                //calculate weight_sizeY
            float sizeY=aCluster.sizeY();
            float sizeY_up =  m_PixelCellHeightOverWidth*std::abs(jetZOverRho)+m_maxSizeY_q;
            float sizeY_peak = m_PixelCellHeightOverWidth*std::abs(jetZOverRho)+m_peakSizeY_q; 
            float sizeY_down = m_PixelCellHeightOverWidth*std::abs(jetZOverRho)+m_minSizeY_q;
            float weight_sizeY_up = (sizeY_up-sizeY)/(sizeY_up-sizeY_peak);
            float weight_sizeY_down = (sizeY-sizeY_down)/(sizeY_peak-sizeY_down);
            weight_sizeY_down = weight_sizeY_down *(weight_sizeY_down>0)*(weight_sizeY_down<1);
            weight_sizeY_up = weight_sizeY_up *(weight_sizeY_up>0)*(weight_sizeY_up<1);
            float weight_sizeY = weight_sizeY_up + weight_sizeY_down;                   

                //calculate weight_rho
                float rho = sqrt(v_bs.x()*v_bs.x() +  v_bs.y()*v_bs.y());
            float weight_rho = ((m_weight_rho_up - rho)/m_weight_rho_up);
                
                //calculate weight_dPhi
                float weight_dPhi = exp(-  std::abs(deltaPhi((*jit)->momentum().Phi(),v_bs.phi()))/m_weight_dPhi );
                
                //calculate weight_sizeX1
            float weight_sizeX1= (aCluster.sizeX()==2) + (aCluster.sizeX()==1)*m_weight_SizeX1;

                //calculate weight_charge
            float charge=aCluster.charge();
            float weightCluster_up = (m_weight_charge_up-charge)/(m_weight_charge_up-m_weight_charge_peak);
            float weightCluster_down = (charge-m_weight_charge_down)/(m_weight_charge_peak-m_weight_charge_down);
            weightCluster_down = weightCluster_down *(weightCluster_down>0)*(weightCluster_down<1);
            weightCluster_up = weightCluster_up *(weightCluster_up>0)*(weightCluster_up<1);
            float weight_charge = weightCluster_up +  weightCluster_down;
        
                //calculate the final weight
            weight = weight_dPhi * weight_sizeY  * weight_rho * weight_sizeX1 * weight_charge ;
            }
            else if(std::abs(modulepos.z())>barrel_lenght) // EC
            {// EC
                //calculate weight_dPhi
                float weight_dPhi = exp(-  std::abs(deltaPhi((*jit)->momentum().Phi(),v_bs.phi())) /m_weight_dPhi_EC );
                //calculate the final weight
            weight=  m_EC_weight*(weight_dPhi) ;                
            }
           if(m_ptWeighting)    weight=weight*pt_weight; 
            zWeights.push_back(weight); //add the weight to zWeights
          } 
        }//if it pass DeltaPhi(Jet,Cluster) requirements
          }
    }//loop on pixel clusters on this module
      }//if it is a compatible module
    }//loop on pixel modules
   jet_count++;
   }//loop on selected jets
   
  //order zProjections and zWeights by z
  std::multimap<float,float>  zWithW;
  size_t i=0;
  for(i=0;i<zProjections.size();i++) zWithW.insert(std::pair<float,float>(zProjections[i],zWeights[i])); 
  i=0;
  for(std::multimap<float,float>::iterator it=zWithW.begin(); it!=zWithW.end(); it++,i++) { zProjections[i]=it->first; zWeights[i]=it->second; } //order zProjections and zWeights by z
   

  //calculate zWeightsSquared
  std::vector<float> zWeightsSquared;
  for(std::vector<float>::iterator it=zWeights.begin();it!=zWeights.end();it++) {zWeightsSquared.push_back((*it)*(*it));}

  //do multi-step peak searching
  float res_step1 = FindPeakFastPV(  zProjections, zWeights, 0.0, m_zClusterWidth_step1, 999.0, -1.0);    
  float res_step2 = FindPeakFastPV(  zProjections, zWeights, res_step1, m_zClusterWidth_step2, m_zClusterSearchArea_step2, m_weightCut_step2);  
  float res_step3 = FindPeakFastPV(  zProjections, zWeightsSquared, res_step2, m_zClusterWidth_step3, m_zClusterSearchArea_step3, m_weightCut_step3*m_weightCut_step3);  

  float centerWMax=res_step3;
  //End of PART II
  
  //Make the output
  float res=0; 
  if(zProjections.size() > 2) 
  {
     res=centerWMax;
     Vertex::Error e; 
     e(0, 0) = 0.0015 * 0.0015;
     e(1, 1) = 0.0015 * 0.0015;
     e(2, 2) = 1.5 * 1.5;
     Vertex::Point p(beamSpot->x(res), beamSpot->y(res), res);
     Vertex thePV(p, e, 1, 1, 0);
     auto pOut = std::make_unique<reco::VertexCollection>();
     pOut->push_back(thePV);
     iEvent.put(std::move(pOut));
   } else
   {
     Vertex::Error e;
     e(0, 0) = 0.0015 * 0.0015;
     e(1, 1) = 0.0015 * 0.0015;
     e(2, 2) = 1.5 * 1.5;
     Vertex::Point p(beamSpot->x(res), beamSpot->y(res), res);
     Vertex thePV(p, e, 0, 0, 0);
     auto pOut = std::make_unique<reco::VertexCollection>();
     pOut->push_back(thePV);
     iEvent.put(std::move(pOut));
   }

//Finally, calculate the zClusterQuality as Sum(weights near the fastPV)/sqrt(Sum(total weights)) [a kind of zCluster significance]

 

  const float half_width_peak=1;
  float nWeightedTot=0;
  float nWeightedTotPeak=0;
  for(std::vector<float>::iterator it = zProjections.begin();it!=zProjections.end(); it++)
  {
    nWeightedTot+=zWeights[it-zProjections.begin()]; 
    if((res-half_width_peak)<=(*it) && (*it)<=(res+half_width_peak))
    {
        nWeightedTotPeak+=zWeights[it-zProjections.begin()];    
    }
  }

auto zClusterQuality = std::make_unique<float>();
*zClusterQuality=-1;
if(nWeightedTot!=0)
{
   *zClusterQuality=nWeightedTotPeak / sqrt(nWeightedTot/(2*half_width_peak)); // where 30 is the beam spot lenght 
   iEvent.put(std::move(zClusterQuality));
}
else
   iEvent.put(std::move(zClusterQuality));

}

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> FastPrimaryVertexWithWeightsProducer::beamSpotToken

private

Definition at line 70 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by FastPrimaryVertexWithWeightsProducer(), and produce().

edm::EDGetTokenT<SiPixelClusterCollectionNew> FastPrimaryVertexWithWeightsProducer::clustersToken

private

Definition at line 69 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by FastPrimaryVertexWithWeightsProducer(), and produce().

edm::EDGetTokenT<edm::View<reco::Jet> > FastPrimaryVertexWithWeightsProducer::jetsToken

private

Definition at line 71 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by FastPrimaryVertexWithWeightsProducer(), and produce().

const bool FastPrimaryVertexWithWeightsProducer::m_barrel

private

Definition at line 77 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const edm::InputTag FastPrimaryVertexWithWeightsProducer::m_clusters

private

Definition at line 67 of file FastPrimaryVertexWithWeightsProducer.cc.

const double FastPrimaryVertexWithWeightsProducer::m_EC_weight

private

Definition at line 102 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const bool FastPrimaryVertexWithWeightsProducer::m_endCap

private

Definition at line 96 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_maxDeltaPhi

private

Definition at line 79 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_maxDeltaPhi_EC

private

Definition at line 99 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_maxJetEta

private

Definition at line 75 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_maxJetEta_EC

private

Definition at line 98 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_maxSizeX

private

Definition at line 78 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_maxSizeY_q

private

Definition at line 84 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_maxZ

private

Definition at line 66 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_minJetEta_EC

private

Definition at line 97 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_minJetPt

private

Definition at line 76 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_minSizeY_q

private

Definition at line 83 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const int FastPrimaryVertexWithWeightsProducer::m_njets

private

Definition at line 74 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_peakSizeY_q

private

Definition at line 93 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_PixelCellHeightOverWidth

private

Definition at line 82 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

std::string FastPrimaryVertexWithWeightsProducer::m_pixelCPE

private

Definition at line 68 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const bool FastPrimaryVertexWithWeightsProducer::m_ptWeighting

private

Definition at line 120 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_ptWeighting_offset

private

Definition at line 122 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_ptWeighting_slope

private

Definition at line 121 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weight_charge_down

private

Definition at line 80 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weight_charge_peak

private

Definition at line 92 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weight_charge_up

private

Definition at line 81 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weight_dPhi

private

Definition at line 89 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weight_dPhi_EC

private

Definition at line 103 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weight_rho_up

private

Definition at line 91 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weight_SizeX1

private

Definition at line 90 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weightCut_step2

private

Definition at line 112 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_weightCut_step3

private

Definition at line 117 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_zClusterSearchArea_step2

private

Definition at line 111 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_zClusterSearchArea_step3

private

Definition at line 116 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_zClusterWidth_step1

private

Definition at line 107 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_zClusterWidth_step2

private

Definition at line 110 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().

const double FastPrimaryVertexWithWeightsProducer::m_zClusterWidth_step3

private

Definition at line 115 of file FastPrimaryVertexWithWeightsProducer.cc.

Referenced by produce().