27 : inputJetTag_(iConfig.getParameter<
edm::
InputTag>(
"inputJetTag")) {
28 produces<std::vector<double>>();
46 std::unique_ptr<std::vector<double>>
result(
new std::vector<double>);
49 std::vector<TLorentzVector> myJets;
50 reco::CaloJetCollection::const_iterator jetIt;
51 for (jetIt = calojet_handle->begin(); jetIt != calojet_handle->end(); jetIt++) {
53 j.SetPtEtaPhiE(jetIt->pt(), jetIt->eta(), jetIt->phi(), jetIt->energy());
67 std::vector<double> ETs;
71 for (
unsigned int i = 0;
i <
jets.size();
i++) {
72 if (
jets[
i].Et() > 50. && fabs(
jets[
i].Eta()) < 2.5)
73 ETs.push_back(
jets[
i].Et());
90 std::vector<double>
diff(1 << (ETs.size() - 1), 0.);
91 for (
unsigned i = 0;
i <
diff.size();
i++)
92 for (
unsigned j = 0;
j < ETs.size();
j++)
93 diff[
i] += ETs[
j] * (1 - 2 * (
int(
i >>
j) & 1));
94 std::vector<double>::const_iterator it;
96 for (it =
diff.begin(); it !=
diff.end(); it++)
108 for (
unsigned int i = 0;
i <
jets.size();
i++) {
109 if (
jets[
i].Et() > 50. && fabs(
jets[
i].Eta()) < 2.5)
118 for (
unsigned int i = 0;
i <
jets.size();
i++) {
119 if (
jets[
i].Et() > 50. && fabs(
jets[
i].Eta()) < 2.5)
T getParameter(std::string const &) const
static double CalcMHT(const std::vector< TLorentzVector > &)
void produce(edm::StreamID, edm::Event &, const edm::EventSetup &) const override
static double alphaT(const double, const double, const double)
edm::InputTag inputJetTag_
AlphaTVarProducer(const edm::ParameterSet &)
#define DEFINE_FWK_MODULE(type)
edm::EDGetTokenT< reco::CaloJetCollection > inputJetTagToken_
static double deltaHt(const std::vector< double > &)
static double CalcHT(const std::vector< TLorentzVector > &)
double CalcAlphaT(const std::vector< TLorentzVector > &) const