MuellerNavelet.h

Go to the documentation of this file.

#pragma once
 
#include "Core/Objects/interface/Jet.h"
 
#include <algorithm>
#include <list>
#include <optional>
#include <vector>
 
namespace DAS::MN {
 
struct Obs2Jets {
 
    Obs2Jets ( const DAS::RecJet & Fwd, const DAS::RecJet & Bwd );
    Obs2Jets ( const DAS::GenJet & Fwd, const DAS::GenJet & Bwd );
 
    float PtAve () const;
    float PtFWD () const;
    float PtBWD () const;
    float PtRatMN () const;
    float DEta () const;
    float DPhi () const;
    float CosDPhi (int n 
            ) const;
    size_t size () const;
 
    const float weightFWD, 
                weightBWD; 
    const float weight;    
 
private:
    std::vector<DAS::FourVector> LeadingJets; 
};
 
struct ObsMiniJets {
 
    ObsMiniJets ( const std::vector<DAS::RecJet> & recJets );
    ObsMiniJets ( const std::vector<DAS::GenJet> & genJets );
 
    float PtAveMini () const;
    float DEtaAveMini () const;
    float REtaAveMini () const;
    float RPtExpDEta () const;
    size_t size () const;
 
    const float weight;
 
private:
    std::vector<DAS::FourVector> MiniJets; 
};// end of class ObsMiniJets
 
template<typename Jet>
std::optional<std::pair<Jet,Jet>> GetMNJet (std::vector<Jet> jets, std::function<bool(Jet&)> ptcut=[](Jet& jet) {return jet.p4.Pt() < 35 ;});
 
template<typename Jet>
std::optional<std::pair<Jet, Jet>> GetMNJet (std::vector <Jet> jets, std::function<bool(Jet&)> ptcut)
{
    jets.erase(remove_if(jets.begin(), jets.end(), ptcut), jets.end());
    if (jets.size() < 2) return std::nullopt;
    auto result = minmax_element( jets.begin(), jets.end() , [](Jet& j1, Jet& j2)
            {return j1.p4.Eta() > j2.p4.Eta();} );  
 
    if (result.first->p4.Eta() <= 0 || result.second->p4.Eta() >= 0) return std::nullopt;
    return std::make_optional<std::pair<Jet,Jet>>(*result.first, *result.second);
}
 
template<typename Jet>
std::vector<Jet> GetMiniJets (std::vector<Jet> jets, const std::pair<Jet,Jet>& MNJets, std::function<bool(Jet&)> ptcut=[](Jet& jet) {return jet.p4.Pt()<20;});
 
template<typename Jet>
std::vector<Jet> GetMiniJets (std::vector<Jet> jets, const std::pair<Jet,Jet>& MNJets, std::function<bool(Jet&)> ptcut)
{
    jets.erase(std::remove_if(jets.begin(), jets.end(), ptcut), jets.end());
    jets.erase(std::remove_if(jets.begin(), jets.end(),
               [MNJets](Jet& jet){return jet.p4.Eta() >= MNJets.first.p4.Eta();}), jets.end());
 
    jets.erase(std::remove_if(jets.begin(), jets.end(),
               [MNJets](Jet& jet){return jet.p4.Eta() <= MNJets.second.p4.Eta();}), jets.end());
 
    std::sort( jets.begin(), jets.end() , [](Jet& j1, Jet& j2)
            {return j1.p4.Eta() > j2.p4.Eta();} );  
 
    return jets;
}
 
} // end of namespace DAS::MN