#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

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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, ModuleToResolverIndicies 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
void	produce (edm::Event &, const edm::EventSetup &) override

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
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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
	)

overrideprivate

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().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation