common.h

Go to the documentation of this file.

#pragma once
 
#include <cmath>
#include <type_traits>
#include <vector>
 
namespace DAS {
 
static const std::vector<double> binsY = {0., 0.25, .5, 0.75, 1., 1.25, 1.5, 1.75, 2., 2.25, 2.5, 2.75, 3., 3.25, 3.5, 3.75, 4., 4.25, 4.5, 4.75, 5., 5.25, 5.5, 5.75, 6., 6.25, 6.5, 6.75, 7., 7.25, 7.5, 7.75, 8., 8.25, 8.5, 8.75, 9., 9.25, 9.5};
static const std::vector<double> binsPhi = {0.,.1*M_PI, .2*M_PI, .3*M_PI, .4*M_PI, .5*M_PI, .6*M_PI, .7*M_PI, .8*M_PI, .9*M_PI, M_PI};
 
} // end of DAS namespace
 
namespace DAS::TPS {
 
template<typename T, typename = void>
struct has_p4 : std::false_type {};
template<typename T>
struct has_p4<T, std::void_t<decltype(std::declval<T>().p4.Eta())>> : std::true_type {};
template<class Jet>
std::vector<Jet> PhaseSpaceSelection(const std::vector<Jet>& jets, float minpt) {
    std::vector<Jet> tempjets;
    for (const auto& jet : jets) {
        double eta;
        if constexpr(has_p4<Jet>::value) {
            eta = jet.p4.Eta();
        } else {
            eta = jet.eta();
        }
        if (jet.CorrPt() > minpt && eta < 5.0)
            tempjets.push_back(jet);
    }
    return tempjets;
}
 
} // end of DAS::TPS namespace