TtHadLRSignalSelObservables.h

Go to the documentation of this file.

#ifndef TtHadLRSignalSelObservables_h
#define TtHadLRSignalSelObservables_h

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"

#include <iostream>
#include <string>
#include <vector>

#include "TLorentzVector.h"
#include "TVector.h"
#include "TVectorD.h"
#include "TMatrix.h"
#include "TMatrixDSymEigen.h"
#include "TMatrixDSym.h"
#include "TMatrixTSym.h"

#include "AnalysisDataFormats/TopObjects/interface/TtHadEvtSolution.h"
#include "DataFormats/PatCandidates/interface/Jet.h"

const double PI = 3.14159265;

class TtHadLRSignalSelObservables {
public:
  TtHadLRSignalSelObservables();
  ~TtHadLRSignalSelObservables();

  void operator()(TtHadEvtSolution&);

private:
  // compare two jets in ET
  struct CompareET {
    bool operator()(const pat::Jet& j1, const pat::Jet& j2) const { return j1.et() > j2.et(); }
  };

  CompareET EtComparator;

  // compare two jets in bdisc
  struct CompareBdisc {
    bool operator()(const pat::Jet& j1, const pat::Jet& j2) const {
      return j1.bDiscriminator("trackCountingJetTags") > j2.bDiscriminator("trackCountingJetTags");
    }
  };

  CompareBdisc BdiscComparator;

  // compare two double
  struct CompareDouble {
    bool operator()(double j1, double j2) const { return j1 > j2; }
  };

  CompareDouble dComparator;

  std::vector<std::pair<unsigned int, double> > evtselectVarVal;
};

#endif