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 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	Mjj
	)

Definition at line 80 of file L1TJetsMatching.h.

References mps_fire::i, convertSQLitetoXML_cfg::output, hiDetachedQuadStep_cff::pt1, and hiDetachedQuadStep_cff::pt2.

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() < pt1) && (myJet1.pt() > pt2) && (myJet2.pt() > pt2))
            {
                
                const T &  myJetTest = (pfMatchedJets)[0];
                if (myJetTest.pt()>pt1){
                    output.second.push_back(myJet1);
                    output.second.push_back(myJet2);
                    output.second.push_back(myJetTest);
                    
                }
            }
            
        }
        
        return output;
        
    }