PFCand_AssoMapAlgos.cc

Go to the documentation of this file.

#include "CommonTools/RecoUtils/interface/PFCand_AssoMapAlgos.h"

#include "TrackingTools/IPTools/interface/IPTools.h"

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

/*************************************************************************************/
/* dedicated constructor for the algorithms                                          */
/*************************************************************************************/

PFCand_AssoMapAlgos::PFCand_AssoMapAlgos(const edm::ParameterSet& iConfig, edm::ConsumesCollector&& iC)
    : PF_PU_AssoMapAlgos(iConfig, iC) {
  input_MaxNumAssociations_ = iConfig.getParameter<int>("MaxNumberOfAssociations");

  token_VertexCollection_ = iC.consumes<VertexCollection>(iConfig.getParameter<InputTag>("VertexCollection"));

  token_BeamSpot_ = iC.consumes<BeamSpot>(iConfig.getParameter<InputTag>("BeamSpot"));
}

/*************************************************************************************/
/* get all needed collections at the beginning                                       */
/*************************************************************************************/

void PFCand_AssoMapAlgos::GetInputCollections(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  PF_PU_AssoMapAlgos::GetInputCollections(iEvent, iSetup);

  //get the offline beam spot
  iEvent.getByToken(token_BeamSpot_, beamspotH);

  //get the input vertex collection
  iEvent.getByToken(token_VertexCollection_, vtxcollH);

  iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
}

/*************************************************************************************/
/* create the pf candidate to vertex association and the inverse map                 */
/*************************************************************************************/
std::pair<std::unique_ptr<PFCandToVertexAssMap>, std::unique_ptr<VertexToPFCandAssMap>>
PFCand_AssoMapAlgos::createMappings(edm::Handle<reco::PFCandidateCollection> pfCandH, const edm::EventSetup& iSetup) {
  unique_ptr<PFCandToVertexAssMap> pfcand2vertex(new PFCandToVertexAssMap(vtxcollH, pfCandH));
  unique_ptr<VertexToPFCandAssMap> vertex2pfcand(new VertexToPFCandAssMap(pfCandH, vtxcollH));

  int num_vertices = vtxcollH->size();
  if (num_vertices < input_MaxNumAssociations_)
    input_MaxNumAssociations_ = num_vertices;
  vector<VertexRef> vtxColl_help;
  if (input_MaxNumAssociations_ == 1)
    vtxColl_help = CreateVertexVector(vtxcollH);

  for (unsigned i = 0; i < pfCandH->size(); i++) {
    PFCandidateRef candref(pfCandH, i);

    if (input_MaxNumAssociations_ > 1)
      vtxColl_help = CreateVertexVector(vtxcollH);

    VertexPfcQuality VtxPfcQual;

    TrackRef PFCtrackref = candref->trackRef();

    if (PFCtrackref.isNull()) {
      for (int assoc_ite = 0; assoc_ite < input_MaxNumAssociations_; ++assoc_ite) {
        int quality = -1 - assoc_ite;

        // Insert the best vertex and the pair of track and the quality of this association in the map
        pfcand2vertex->insert(vtxColl_help.at(0), make_pair(candref, quality));
        vertex2pfcand->insert(candref, make_pair(vtxColl_help.at(0), quality));

        //cleanup only if multiple iterations are made
        if (input_MaxNumAssociations_ > 1)
          PF_PU_AssoMapAlgos::EraseVertex(vtxColl_help, vtxColl_help.at(0));
      }

    } else {
      TransientTrack transtrk(PFCtrackref, &(*bFieldH));
      transtrk.setBeamSpot(*beamspotH);
      transtrk.setES(iSetup);

      for (int assoc_ite = 0; assoc_ite < input_MaxNumAssociations_; ++assoc_ite) {
        VertexStepPair assocVtx = FindAssociation(PFCtrackref, vtxColl_help, bFieldH, iSetup, beamspotH, assoc_ite);
        int step = assocVtx.second;
        double distance = (IPTools::absoluteImpactParameter3D(transtrk, *(assocVtx.first))).second.value();

        int quality = DefineQuality(assoc_ite, step, distance);

        // Insert the best vertex and the pair of track and the quality of this association in the map
        pfcand2vertex->insert(assocVtx.first, make_pair(candref, quality));
        vertex2pfcand->insert(candref, make_pair(assocVtx.first, quality));

        //cleanup only if multiple iterations are made
        if (input_MaxNumAssociations_ > 2)
          PF_PU_AssoMapAlgos::EraseVertex(vtxColl_help, assocVtx.first);
      }

    }  //check PFCtrackref.isNull

  }  //i on pfCandH

  return {std::move(pfcand2vertex), std::move(vertex2pfcand)};
}

/*************************************************************************************/
/* create the pf candidate to vertex association map                                 */
/*************************************************************************************/

std::unique_ptr<PFCandToVertexAssMap> PFCand_AssoMapAlgos::CreatePFCandToVertexMap(
    edm::Handle<reco::PFCandidateCollection> pfCandH, const edm::EventSetup& iSetup) {
  return createMappings(pfCandH, iSetup).first;
}

/*************************************************************************************/
/* create the vertex to pf candidate association map                                 */
/*************************************************************************************/

std::unique_ptr<VertexToPFCandAssMap> PFCand_AssoMapAlgos::CreateVertexToPFCandMap(
    edm::Handle<reco::PFCandidateCollection> pfCandH, const edm::EventSetup& iSetup) {
  return createMappings(pfCandH, iSetup).second;
}

/*************************************************************************************/
/* create the vertex to pf candidate association map                                 */
/*************************************************************************************/

std::unique_ptr<PFCandToVertexAssMap> PFCand_AssoMapAlgos::SortPFCandAssociationMap(
    PFCandToVertexAssMap* pfcvertexassInput, edm::EDProductGetter const* getter) {
  //create a new PFCandVertexAssMap for the Output which will be sorted
  unique_ptr<PFCandToVertexAssMap> pfcvertexassOutput(new PFCandToVertexAssMap(getter));

  //Create and fill a vector of pairs of vertex and the summed (pT)**2 of the pfcandidates associated to the vertex
  VertexPtsumVector vertexptsumvector;

  //loop over all vertices in the association map
  for (PFCandToVertexAssMap::const_iterator assomap_ite = pfcvertexassInput->begin();
       assomap_ite != pfcvertexassInput->end();
       assomap_ite++) {
    const VertexRef assomap_vertexref = assomap_ite->key;
    const PFCandQualityPairVector pfccoll = assomap_ite->val;

    float ptsum = 0;

    PFCandidateRef pfcandref;

    //get the pfcandidates associated to the vertex and calculate the pT**2
    for (unsigned int pfccoll_ite = 0; pfccoll_ite < pfccoll.size(); pfccoll_ite++) {
      pfcandref = pfccoll[pfccoll_ite].first;
      int quality = pfccoll[pfccoll_ite].second;

      if ((quality <= 2) && (quality != -1))
        continue;

      double man_pT = pfcandref->pt();
      if (man_pT > 0.)
        ptsum += man_pT * man_pT;
    }

    vertexptsumvector.push_back(make_pair(assomap_vertexref, ptsum));
  }

  while (!vertexptsumvector.empty()) {
    VertexRef vertexref_highestpT;
    float highestpT = 0.;
    int highestpT_index = 0;

    for (unsigned int vtxptsumvec_ite = 0; vtxptsumvec_ite < vertexptsumvector.size(); vtxptsumvec_ite++) {
      if (vertexptsumvector[vtxptsumvec_ite].second > highestpT) {
        vertexref_highestpT = vertexptsumvector[vtxptsumvec_ite].first;
        highestpT = vertexptsumvector[vtxptsumvec_ite].second;
        highestpT_index = vtxptsumvec_ite;
      }
    }

    //loop over all vertices in the association map
    for (PFCandToVertexAssMap::const_iterator assomap_ite = pfcvertexassInput->begin();
         assomap_ite != pfcvertexassInput->end();
         assomap_ite++) {
      const VertexRef assomap_vertexref = assomap_ite->key;
      const PFCandQualityPairVector pfccoll = assomap_ite->val;

      //if the vertex from the association map the vertex with the highest pT
      //insert all associated pfcandidates in the output Association Map
      if (assomap_vertexref == vertexref_highestpT)
        for (unsigned int pfccoll_ite = 0; pfccoll_ite < pfccoll.size(); pfccoll_ite++)
          pfcvertexassOutput->insert(assomap_vertexref, pfccoll[pfccoll_ite]);
    }

    vertexptsumvector.erase(vertexptsumvector.begin() + highestpT_index);
  }

  return pfcvertexassOutput;
}