common.h

Go to the documentation of this file.

#pragma once
 
#include <vector>
#include <cmath>
 
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<class Jet>
std::vector<Jet> PhaseSpaceSelection(const std::vector<Jet>& jets, float minpt) {
    std::vector<Jet> tempjets;
    for (const auto& jet : jets) 
        if (jet.CorrPt() > minpt && jet.p4.Eta() < 5.0) 
            tempjets.push_back(jet); 
    return tempjets;
}
 
} // end of DAS::TPS namespace