EMJetMatching.cpp

Go to the documentation of this file.

#include <iostream>
#include "caf_util/EMJetMatching.hpp"
#include "cafe/Config.hpp"
#include "kinem_util/AnglesUtil.hpp"
#include "jetcorr/CalTool.hpp"

using namespace cafe;
using namespace std;

namespace caf_util {

  EMJetMatching::EMJetMatching(const char *name) : cafe::SelectUserObjects<TMBJet>(name), _event(0), _vars("_CalDetectorPhi","_CalDetectorEta")
  {
    cafe::Config config(name);
    _jetBranch  = config.get("From",""); 
    _electronBranch = config.get("electronBranch","EMscone"); 
    _dR_Jet_EM  = config.get("DR_Jet_EM",0.4);
    _debug = config.get("Debug", 0);
  }

  bool EMJetMatching::processEvent(cafe::Event &event)
  {
    _event = &event ;
    SelectUserObjects<TMBJet>::processEvent(event) ;
    return true  ;
  }

  bool EMJetMatching::selectObject(const TMBJet &jet)
  {
    return true ;
  }

  void EMJetMatching::after(cafe::Collection<TMBJet>& input, cafe::Collection<TMBJet>& rejected)
  {
    Collection<TMBEMCluster> emcol(_event->getCollection<TMBEMCluster>(_electronBranch.c_str(),_vars));
    static const float phiConv = std::acos(-1.)/32.;

    for (Collection<TMBJet>::iterator it_jet = input.begin() ;
         it_jet != input.end() ; ++ it_jet)  {

      //      if (!it_jet->isGood()) continue;
      float dRmin = kinem::TWOPI;

      for (Collection<TMBEMCluster>::const_iterator it_em = emcol.begin() ; 
        it_em != emcol.end() ; ++ it_em) {

         float dphi = kinem::delta_phi(phiConv*it_jet->detPhi(), it_em->CalDetectorPhi());
         float deta = fabs(CalTool::EtaDToEtaDet(it_jet->detEta()) - it_em->CalDetectorEta());
         float dR = sqrt(dphi*dphi + deta*deta);
         if (dR < dRmin) dRmin = dR;

      } // for EM

      if ( dRmin < _dR_Jet_EM ) it_jet->SetFlag(it_jet->isGood(), it_jet->isL1Conf(), it_jet->isBad(),
                                        it_jet->isNoise(), true, it_jet->hasMU(), 
                                        it_jet->isSmeared(), it_jet->isSmearedMU());
      else it_jet->SetFlag(it_jet->isGood(), it_jet->isL1Conf(), it_jet->isBad(),
                                        it_jet->isNoise(), false, it_jet->hasMU(), 
                                        it_jet->isSmeared(), it_jet->isSmearedMU());

    } // for jet

  }

} // namespace

ClassImp(caf_util::EMJetMatching);
 

---