JetCalibrationMethods.cc

Go to the documentation of this file.

// JetCalibrationMethods.cc
// Authors: Inga Bucinskaite, UIC
//
// This file should contain the different algorithms used for Jet Calibration
//

#include "L1Trigger/L1TCalorimeter/interface/JetCalibrationMethods.h"
#include <vector>

namespace l1t {

  bool compareCorrJets (l1t::Jet m, l1t::Jet n){
    return (m.hwPt() < n.hwPt() );
  }

  void JetCalibration1(std::vector<l1t::Jet> * uncalibjets,
               std::vector<double> jetSF,
               std::vector<l1t::Jet> * jets,
               bool applyJetCalibration,
               double jetLSB) 
  {

    for (std::vector<l1t::Jet>::const_iterator uncalibjet = uncalibjets->begin(); uncalibjet != uncalibjets->end(); ++uncalibjet){

      if (!applyJetCalibration) {
    l1t::Jet corrjets = *uncalibjet;
    jets->push_back(corrjets);
    continue;
      }
    
      int jetPt = (uncalibjet->hwPt())*jetLSB;  // correction factors are parameterized as functions of physical pt
      int jetPhi = uncalibjet->hwPhi();
      int jetEta = uncalibjet->hwEta();
      int jetQual = uncalibjet->hwQual();
      double jpt = 0.0;
      
      double alpha = jetSF[2*jetEta + 0]; //Scale factor (See jetSF_cfi.py)
      double gamma = ((jetSF[2*jetEta + 1])); //Offset
      
      jpt = jetPt*alpha+gamma;
      unsigned int corjetET =(int) (jpt/jetLSB);
      
      ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > *jetLorentz =
    new ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> >();
      l1t::Jet corrjets(*jetLorentz, corjetET, jetEta, jetPhi, jetQual);
      
      jets->push_back(corrjets);
      
      
    }
    
    std::sort(jets->begin(), jets->end(), compareCorrJets);
    std::reverse(jets->begin(), jets->end());
  }
}