L1TJetsMatching.h File Reference

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/TauReco/interface/PFTauFwd.h"
#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
#include "Math/GenVector/VectorUtil.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
#include "FWCore/Utilities/interface/EDMException.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "HLTrigger/HLTcore/interface/defaultModuleLabel.h"
#include "DataFormats/Math/interface/deltaR.h"
#include <map>
#include <vector>

Go to the source code of this file.

Classes
class	L1TJetsMatching< T >

Functions
template<typename T >
std::pair< std::vector< T > , std::vector< T > >	categorise (const std::vector< T > &pfMatchedJets, double pt1, double pt2, double pt3, double Mjj)

Function Documentation

template<typename T >

std::pair<std::vector<T>, std::vector<T> > categorise	(	const std::vector< T > &	pfMatchedJets,
		double	pt1,
		double	pt2,
		double	pt3,
		double	Mjj
	)

Definition at line 71 of file L1TJetsMatching.h.

References mps_fire::i, dqmiolumiharvest::j, and convertSQLitetoXML_cfg::output.

Referenced by L1TJetsMatching< T >::produce().

                                                                                       {
  std::pair<std::vector<T>, std::vector<T>> output;
  unsigned int i1 = 0;
  unsigned int i2 = 0;
  double mjj = 0;
  if (pfMatchedJets.size() > 1) {
    for (unsigned int i = 0; i < pfMatchedJets.size() - 1; i++) {
      const T& myJet1 = (pfMatchedJets)[i];

      for (unsigned int j = i + 1; j < pfMatchedJets.size(); j++) {
        const T& myJet2 = (pfMatchedJets)[j];

        const double mjj_test = (myJet1.p4() + myJet2.p4()).M();

        if (mjj_test > mjj) {
          mjj = mjj_test;
          i1 = i;
          i2 = j;
        }
      }
    }

    const T& myJet1 = (pfMatchedJets)[i1];
    const T& myJet2 = (pfMatchedJets)[i2];

    if ((mjj > Mjj) && (myJet1.pt() >= pt1) && (myJet2.pt() > pt2)) {
      output.first.push_back(myJet1);
      output.first.push_back(myJet2);
    }

    if ((mjj > Mjj) && (myJet1.pt() < pt3) && (myJet1.pt() > pt2) && (myJet2.pt() > pt2)) {
      const T& myJetTest = (pfMatchedJets)[0];
      if (myJetTest.pt() > pt3) {
        output.second.push_back(myJet1);
        output.second.push_back(myJet2);
        output.second.push_back(myJetTest);
      }
    }
  }

  return output;
}