da/d7b/FastPrimaryVertexWithWeightsProducer_8cc_source.html

 // -*- C++ -*-

 //

 // Package:    FastPrimaryVertexWithWeightsProducer

 // Class:      FastPrimaryVertexWithWeightsProducer

 //

 //

 // Original Author:  Silvio DONATO

 //         Created:  Wed Dec 18 10:05:40 CET 2013

 //

 //


 // system include files

 #include <memory>

 #include <vector>

 #include <math.h>

 // user include files

 #include "FWCore/Framework/interface/Frameworkfwd.h"

 #include "FWCore/Framework/interface/stream/EDProducer.h"

 #include "FWCore/Framework/interface/ESHandle.h"

 #include "FWCore/Framework/interface/Event.h"

 #include "FWCore/Framework/interface/MakerMacros.h"

 #include "FWCore/ParameterSet/interface/ParameterSet.h"

 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"

 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"


 #include "DataFormats/JetReco/interface/CaloJet.h"

 #include "DataFormats/JetReco/interface/CaloJetCollection.h"

 #include "DataFormats/Math/interface/deltaPhi.h"

 #include "DataFormats/VertexReco/interface/VertexFwd.h"

 #include "DataFormats/VertexReco/interface/Vertex.h"

 #include "DataFormats/TrackReco/interface/Track.h"

 #include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"

 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"

 #include "DataFormats/TrackingRecHit/interface/TrackingRecHitFwd.h"

 #include "DataFormats/BeamSpot/interface/BeamSpot.h"


 #include "RecoLocalTracker/ClusterParameterEstimator/interface/PixelClusterParameterEstimator.h"

 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"

 #include "RecoLocalTracker/Records/interface/TkPixelCPERecord.h"


 #include "RecoPixelVertexing/PixelVertexFinding/interface/FindPeakFastPV.h"


 #include "FWCore/Utilities/interface/InputTag.h"


 using namespace std;


 class FastPrimaryVertexWithWeightsProducer : public edm::stream::EDProducer<> {

    public:

       explicit FastPrimaryVertexWithWeightsProducer(const edm::ParameterSet&);

       static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);


    private:

       virtual void produce(edm::Event&, const edm::EventSetup&);


   const double m_maxZ;      // Use only pixel clusters with |z| < maxZ

   const edm::InputTag m_clusters;   // PixelClusters InputTag

   std::string m_pixelCPE;   // PixelCPE (PixelClusterParameterEstimator)

   edm::EDGetTokenT<SiPixelClusterCollectionNew> clustersToken;

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

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


 // PARAMETERS USED IN THE BARREL PIXEL CLUSTERS PROJECTION

   const int m_njets;            // Use only the first njets

   const double m_maxJetEta;     // Use only jets with |eta| < maxJetEta

   const double m_minJetPt;      // Use only jets with Pt > minJetPt

   const bool m_barrel;      // Use clusters from pixel endcap

   const double m_maxSizeX;      // Use only pixel clusters with sizeX <= maxSizeX

   const double m_maxDeltaPhi;       // Use only pixel clusters with DeltaPhi(Jet,Cluster) < maxDeltaPhi

   const double m_weight_charge_down;    // Use only pixel clusters with ClusterCharge > weight_charge_down

   const double m_weight_charge_up;  // Use only pixel clusters with ClusterCharge < weight_charge_up

   const double m_PixelCellHeightOverWidth;//It is the ratio between pixel cell height and width along z coordinate about 285µm/150µm=1.9

   const double m_minSizeY_q;        // Use only pixel clusters with sizeY > PixelCellHeightOverWidth * |jetZOverRho| + minSizeY_q

   const double m_maxSizeY_q;        // Use only pixel clusters with sizeY < PixelCellHeightOverWidth * |jetZOverRho| + maxSizeY_q


 // PARAMETERS USED TO WEIGHT THE BARREL PIXEL CLUSTERS

   // The cluster weight is defined as weight = weight_dPhi * weight_sizeY  * weight_rho * weight_sizeX1 * weight_charge


   const double m_weight_dPhi;       // used in weight_dPhi = exp(-|DeltaPhi(JetCluster)|/m_weight_dPhi)

   const double  m_weight_SizeX1;    // used in weight_SizeX1 = (ClusterSizeX==2)*1+(ClusterSizeX==1)*m_weight_SizeX1;

   const double m_weight_rho_up;     // used in weight_rho = (m_weight_rho_up - ClusterRho)/m_weight_rho_up

   const double m_weight_charge_peak;    // Give the maximum weight_charge for a cluster with Charge = m_weight_charge_peak

   const double m_peakSizeY_q;       // Give the maximum weight_sizeY for a cluster with sizeY = PixelCellHeightOverWidth * |jetZOverRho| + peakSizeY_q


 // PARAMETERS USED IN THE ENDCAP PIXEL CLUSTERS PROJECTION

   const bool m_endCap;      // Use clusters from pixel endcap

   const double m_minJetEta_EC;  // Use only jets with |eta| > minJetEta_EC

   const double m_maxJetEta_EC;  // Use only jets with |eta| < maxJetEta_EC

   const double m_maxDeltaPhi_EC;    // Use only pixel clusters with DeltaPhi(Jet,Cluster) < maxDeltaPhi_EC


 // PARAMETERS USED TO WEIGHT THE ENDCAP PIXEL CLUSTERS

   const double m_EC_weight;     // In EndCap the weight is defined as weight = m_EC_weight*(weight_dPhi)

   const double m_weight_dPhi_EC;    // Used in weight_dPhi = exp(-|DeltaPhi|/m_weight_dPhi_EC )


 // PARAMETERS USED TO FIND THE FASTPV AS PEAK IN THE Z-PROJECTIONS DISTRIBUTION

   // First Iteration: look for a cluster with a width = m_zClusterWidth_step1

   const double m_zClusterWidth_step1;          // cluster width in step1


   // Second Iteration: use only z-projections with weight > weightCut_step2 and look for a cluster with a width = m_zClusterWidth_step2, within of weightCut_step2 of the previous result

   const double m_zClusterWidth_step2;   // cluster width in step2

   const double m_zClusterSearchArea_step2;  // cluster width in step2

   const double m_weightCut_step2;       // minimum z-projections weight required in step2


   // Third Iteration: use only z-projections with weight > weightCut_step3 and look for a cluster with a width = m_zClusterWidth_step3, within of weightCut_step3 of the previous result

   const double m_zClusterWidth_step3;   // cluster width in step3

   const double m_zClusterSearchArea_step3;  // cluster width in step3

   const double m_weightCut_step3;       // minimum z-projections weight required in step3


   // use the jetPt weighting

   const bool m_ptWeighting;             // use weight=weight*pt_weigth ?;

   const double m_ptWeighting_slope;     // pt_weigth= pt*ptWeighting_slope +m_ptWeighting_offset;

   const double m_ptWeighting_offset;        //

 };


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

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


 }


 void

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

 {

   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)

 {


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

      std::auto_ptr<reco::VertexCollection> pOut(new reco::VertexCollection());

      pOut->push_back(thePV);

      iEvent.put(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);

      std::auto_ptr<reco::VertexCollection> pOut(new reco::VertexCollection());

      pOut->push_back(thePV);

      iEvent.put(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()];

     }

   }


 std::auto_ptr<float > zClusterQuality(new float());

 *zClusterQuality=-1;

 if(nWeightedTot!=0)

 {

    *zClusterQuality=nWeightedTotPeak / sqrt(nWeightedTot/(2*half_width_peak)); // where 30 is the beam spot lenght

    iEvent.put(zClusterQuality);

 }

 else

    iEvent.put(zClusterQuality);


 }


 //define this as a plug-in

 DEFINE_FWK_MODULE(FastPrimaryVertexWithWeightsProducer);

math.h

edm::ParameterSet::getParameter
T getParameter(std::string const &) const

FastPrimaryVertexWithWeightsProducer::m_endCap
const bool m_endCap
Definition: FastPrimaryVertexWithWeightsProducer.cc:96

i
int i
Definition: DBlmapReader.cc:9

TrackerDigiGeometryRecord.h

SiPixelCluster::charge
float charge() const
Definition: SiPixelCluster.h:145

patRefSel_refMuJets_cfi.selectedJets
tuple selectedJets
Definition: patRefSel_refMuJets_cfi.py:70

edmNew::DetSetVector< SiPixelCluster >::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:118

SiPixelRawToDigiRegional_cfi.beamSpot
tuple beamSpot
Definition: SiPixelRawToDigiRegional_cfi.py:11

create_public_lumi_plots.exp
tuple exp
Definition: create_public_lumi_plots.py:1087

FastPrimaryVertexWithWeightsProducer::m_weightCut_step2
const double m_weightCut_step2
Definition: FastPrimaryVertexWithWeightsProducer.cc:112

TrackingRecHitFwd.h

FastPrimaryVertexWithWeightsProducer::m_weightCut_step3
const double m_weightCut_step3
Definition: FastPrimaryVertexWithWeightsProducer.cc:117

FastPrimaryVertexWithWeightsProducer::m_minJetPt
const double m_minJetPt
Definition: FastPrimaryVertexWithWeightsProducer.cc:76

createTree.pp
tuple pp
Definition: createTree.py:15

FastPrimaryVertexWithWeightsProducer::m_weight_charge_down
const double m_weight_charge_down
Definition: FastPrimaryVertexWithWeightsProducer.cc:80

FastPrimaryVertexWithWeightsProducer::m_weight_dPhi_EC
const double m_weight_dPhi_EC
Definition: FastPrimaryVertexWithWeightsProducer.cc:103

FastPrimaryVertexWithWeightsProducer::m_maxJetEta_EC
const double m_maxJetEta_EC
Definition: FastPrimaryVertexWithWeightsProducer.cc:98

FastPrimaryVertexWithWeightsProducer::m_weight_charge_up
const double m_weight_charge_up
Definition: FastPrimaryVertexWithWeightsProducer.cc:81

FastPrimaryVertexWithWeightsProducer::m_zClusterWidth_step3
const double m_zClusterWidth_step3
Definition: FastPrimaryVertexWithWeightsProducer.cc:115

edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449

DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17

SiPixelRecHit.h

FastPrimaryVertexWithWeightsProducer::m_maxDeltaPhi_EC
const double m_maxDeltaPhi_EC
Definition: FastPrimaryVertexWithWeightsProducer.cc:99

FastPrimaryVertexWithWeightsProducer::beamSpotToken
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken
Definition: FastPrimaryVertexWithWeightsProducer.cc:70

Event.h

CaloJetCollection.h

rho
Definition: DDAxes.h:10

MakerMacros.h

edm::Handle
Definition: AssociativeIterator.h:47

PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69

FastPrimaryVertexWithWeightsProducer::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: FastPrimaryVertexWithWeightsProducer.cc:180

PV3DBase::y
T y() const
Definition: PV3DBase.h:63

reco::Vertex::Error
math::Error< dimension >::type Error
covariance error matrix (3x3)
Definition: Vertex.h:43

reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9

HLT_25ns14e33_v1_cff.weight_dPhi
tuple weight_dPhi
Definition: HLT_25ns14e33_v1_cff.py:27385

FastPrimaryVertexWithWeightsProducer::m_zClusterWidth_step2
const double m_zClusterWidth_step2
Definition: FastPrimaryVertexWithWeightsProducer.cc:110

FastPrimaryVertexWithWeightsProducer::m_weight_charge_peak
const double m_weight_charge_peak
Definition: FastPrimaryVertexWithWeightsProducer.cc:92

FastPrimaryVertexWithWeightsProducer::m_PixelCellHeightOverWidth
const double m_PixelCellHeightOverWidth
Definition: FastPrimaryVertexWithWeightsProducer.cc:82

VertexFwd.h

FastPrimaryVertexWithWeightsProducer::FastPrimaryVertexWithWeightsProducer
FastPrimaryVertexWithWeightsProducer(const edm::ParameterSet &)
Definition: FastPrimaryVertexWithWeightsProducer.cc:125

findQualityFiles.v
v
Definition: findQualityFiles.py:177

eta
T eta() const
Definition: Basic3DVectorLD.h:177

dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107

HLT_25ns14e33_v1_cff.InputTag
tuple InputTag
Definition: HLT_25ns14e33_v1_cff.py:44427

edm::EDGetTokenT
Definition: EDGetToken.h:32

FastPrimaryVertexWithWeightsProducer::m_weight_dPhi
const double m_weight_dPhi
Definition: FastPrimaryVertexWithWeightsProducer.cc:89

FastPrimaryVertexWithWeightsProducer::m_maxSizeX
const double m_maxSizeX
Definition: FastPrimaryVertexWithWeightsProducer.cc:78

Frameworkfwd.h

TkPixelCPERecord.h

detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15

FastPrimaryVertexWithWeightsProducer::produce
virtual void produce(edm::Event &, const edm::EventSetup &)
Definition: FastPrimaryVertexWithWeightsProducer.cc:224

AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255

edm::ParameterSetDescription
Definition: ParameterSetDescription.h:50

FastPrimaryVertexWithWeightsProducer
Definition: FastPrimaryVertexWithWeightsProducer.cc:56

ParameterSet.h

EDProducer.h

PixelClusterParameterEstimator.h

FastPrimaryVertexWithWeightsProducer::m_ptWeighting_slope
const double m_ptWeighting_slope
Definition: FastPrimaryVertexWithWeightsProducer.cc:121

edm::View< reco::Jet >

iEvent
int iEvent
Definition: GenABIO.cc:230

reco::Vertex
Definition: Vertex.h:34

ParameterSetDescription.h

FastPrimaryVertexWithWeightsProducer::m_weight_SizeX1
const double m_weight_SizeX1
Definition: FastPrimaryVertexWithWeightsProducer.cc:90

edm::Event::put
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113

mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48

edm::ESHandle
Definition: DTSurvey.h:22

fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:42

PV3DBase::z
T z() const
Definition: PV3DBase.h:64

FastPrimaryVertexWithWeightsProducer::m_peakSizeY_q
const double m_peakSizeY_q
Definition: FastPrimaryVertexWithWeightsProducer.cc:93

funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22

ESHandle.h

j
int j
Definition: DBlmapReader.cc:9

edmNew::DetSetVector< SiPixelCluster >

reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:39

patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69

FastPrimaryVertexWithWeightsProducer::m_maxJetEta
const double m_maxJetEta
Definition: FastPrimaryVertexWithWeightsProducer.cc:75

edm::EventSetup
Definition: EventSetup.h:44

FastPrimaryVertexWithWeightsProducer::m_clusters
const edm::InputTag m_clusters
Definition: FastPrimaryVertexWithWeightsProducer.cc:67

edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:89

FastPrimaryVertexWithWeightsProducer::m_weight_rho_up
const double m_weight_rho_up
Definition: FastPrimaryVertexWithWeightsProducer.cc:91

edmNew::DetSet
Definition: DetSetNew.h:21

FastPrimaryVertexWithWeightsProducer::m_njets
const int m_njets
Definition: FastPrimaryVertexWithWeightsProducer.cc:74

Vertex.h

EnergyCorrector.pt
int pt
Definition: EnergyCorrector.py:45

FastPrimaryVertexWithWeightsProducer::jetsToken
edm::EDGetTokenT< edm::View< reco::Jet > > jetsToken
Definition: FastPrimaryVertexWithWeightsProducer.cc:71

DetId
Definition: DetId.h:18

FastPrimaryVertexWithWeightsProducer::m_EC_weight
const double m_EC_weight
Definition: FastPrimaryVertexWithWeightsProducer.cc:102

FastPrimaryVertexWithWeightsProducer::m_zClusterSearchArea_step3
const double m_zClusterSearchArea_step3
Definition: FastPrimaryVertexWithWeightsProducer.cc:116

SiPixelRawToDigiRegional_cfi.deltaPhi
tuple deltaPhi
Definition: SiPixelRawToDigiRegional_cfi.py:9

edm::stream::EDProducer
Definition: EDProducer.h:33

TkPixelCPERecord
Definition: TkPixelCPERecord.h:16

edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55

FastPrimaryVertexWithWeightsProducer::m_barrel
const bool m_barrel
Definition: FastPrimaryVertexWithWeightsProducer.cc:77

edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:65

PixelClusterParameterEstimator
Definition: PixelClusterParameterEstimator.h:15

alignCSCRings.e
list e
Definition: alignCSCRings.py:90

SiPixelCluster::sizeY
int sizeY() const
Definition: SiPixelCluster.h:142

FastPrimaryVertexWithWeightsProducer::m_ptWeighting
const bool m_ptWeighting
Definition: FastPrimaryVertexWithWeightsProducer.cc:120

AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152

Point3DBase< float, GlobalTag >

SiPixelCluster
Pixel cluster – collection of neighboring pixels above threshold.
Definition: SiPixelCluster.h:33

edm::InputTag
Definition: InputTag.h:17

edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:81

FastPrimaryVertexWithWeightsProducer::m_maxDeltaPhi
const double m_maxDeltaPhi
Definition: FastPrimaryVertexWithWeightsProducer.cc:79

deltaPhi.h

InputTag.h

TrackerGeometry.h

BeamSpot.h

FastPrimaryVertexWithWeightsProducer::m_zClusterSearchArea_step2
const double m_zClusterSearchArea_step2
Definition: FastPrimaryVertexWithWeightsProducer.cc:111

FastPrimaryVertexWithWeightsProducer::m_ptWeighting_offset
const double m_ptWeighting_offset
Definition: FastPrimaryVertexWithWeightsProducer.cc:122

FastPrimaryVertexWithWeightsProducer::m_maxZ
const double m_maxZ
Definition: FastPrimaryVertexWithWeightsProducer.cc:66

FastPrimaryVertexWithWeightsProducer::m_pixelCPE
std::string m_pixelCPE
Definition: FastPrimaryVertexWithWeightsProducer.cc:68

edm::ParameterSet
Definition: ParameterSet.h:35

FastPrimaryVertexWithWeightsProducer::m_maxSizeY_q
const double m_maxSizeY_q
Definition: FastPrimaryVertexWithWeightsProducer.cc:84

FindPeakFastPV.h

ConfigurationDescriptions.h

FindPeakFastPV
float FindPeakFastPV(const std::vector< float > &zProjections, const std::vector< float > &zWeights, const float oldVertex, const float m_zClusterWidth, const float m_zClusterSearchArea, const float m_weightCut)
Definition: FindPeakFastPV.h:12

histoStyle.weight
int weight
Definition: histoStyle.py:50

edmNew::DetSet::size
size_type size() const
Definition: DetSetNew.h:86

CaloJet.h

edm::Event
Definition: Event.h:59

HLT_25ns14e33_v1_cff.pixelClusters
tuple pixelClusters
Definition: HLT_25ns14e33_v1_cff.py:4488

Track.h

TrackerDigiGeometryRecord
Definition: TrackerDigiGeometryRecord.h:32

SiPixelCluster.h

SiPixelCluster::sizeX
int sizeX() const
Definition: SiPixelCluster.h:139

FastPrimaryVertexWithWeightsProducer::clustersToken
edm::EDGetTokenT< SiPixelClusterCollectionNew > clustersToken
Definition: FastPrimaryVertexWithWeightsProducer.cc:69

PV3DBase::x
T x() const
Definition: PV3DBase.h:62

TrackerGeometry
Definition: TrackerGeometry.h:29

FastPrimaryVertexWithWeightsProducer::m_zClusterWidth_step1
const double m_zClusterWidth_step1
Definition: FastPrimaryVertexWithWeightsProducer.cc:107

FastPrimaryVertexWithWeightsProducer::m_minSizeY_q
const double m_minSizeY_q
Definition: FastPrimaryVertexWithWeightsProducer.cc:83

edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:27

FastPrimaryVertexWithWeightsProducer::m_minJetEta_EC
const double m_minJetEta_EC
Definition: FastPrimaryVertexWithWeightsProducer.cc:97