HLTTauRefCombiner.cc

Go to the documentation of this file.

#include "HLTriggerOffline/Tau/interface/HLTTauRefCombiner.h"
#include "Math/GenVector/VectorUtil.h"



using namespace edm;
using namespace std;

HLTTauRefCombiner::HLTTauRefCombiner(const edm::ParameterSet& iConfig)
{
  inputColl_   = iConfig.getParameter< std::vector<InputTag> >("InputCollections");
  matchDeltaR_ = iConfig.getParameter<double>("MatchDeltaR");
  outName_     = iConfig.getParameter<string>("OutputCollection");

  produces<LorentzVectorCollection> ( outName_);
}

HLTTauRefCombiner::~HLTTauRefCombiner(){ }

void HLTTauRefCombiner::produce(edm::Event& iEvent, const edm::EventSetup& iES)
{
    auto_ptr<LorentzVectorCollection> out_product(new LorentzVectorCollection);

    //Create The Handles..
    std::vector< Handle<LorentzVectorCollection> > handles;

    bool allCollectionsExist = true;
    //Map the Handles to the collections if all collections exist
    for(size_t i = 0;i<inputColl_.size();++i)
      {
    edm::Handle<LorentzVectorCollection> tmp;
    if(iEvent.getByLabel(inputColl_[i],tmp))
      {
        handles.push_back(tmp);
      }
    else
      {
        allCollectionsExist = false;
      }

      }

    //The reference output object collection will be the first one..
    if(allCollectionsExist)
      {
    //loop on the first collection
    for(size_t i = 0; i < (handles[0])->size();++i)
      {
        bool MatchedObj = true;

        //get reference Vector
        const LorentzVector lvRef = (*(handles[0]))[i];
        
        //Loop on all other collections and match
            for(size_t j = 1; j < handles.size();++j)
          {
            if(!match(lvRef,*(handles[j])))
              MatchedObj = false;

          }
        
        //If the Object is Matched Everywhere store
        if(MatchedObj)
          {
            out_product->push_back(lvRef);
          }


      }

    //Put product to file
    iEvent.put(out_product,outName_);

      }
    
    


}



bool 
HLTTauRefCombiner::match(const LorentzVector& lv,const LorentzVectorCollection& lvcol)
{
 bool matched=false;

 if(lvcol.size()>0)
  for(LorentzVectorCollection::const_iterator it = lvcol.begin();it!=lvcol.end();++it)
   {
          double delta = ROOT::Math::VectorUtil::DeltaR(lv,*it);
      if(delta<matchDeltaR_)
        {
          matched=true;
         
        }
   }



 return matched;
}