TtSemiLRSignalSelObservables.h

Go to the documentation of this file.

#ifndef TtSemiLRSignalSelObservables_h
#define TtSemiLRSignalSelObservables_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 "TVector3.h"
#include "TVectorD.h"

#include "TMatrix.h"
#include "TMatrixDSymEigen.h"
#include "TMatrixDSym.h"
#include "TMatrixTSym.h"

#include "AnalysisDataFormats/TopObjects/interface/TtSemiEvtSolution.h"
#include "TopQuarkAnalysis/TopTools/interface/MEzCalculator.h"
#include "DataFormats/PatCandidates/interface/Jet.h"

const double PI=3.14159265;

class TtSemiLRSignalSelObservables{
  
 public:
  
  TtSemiLRSignalSelObservables();
  ~TtSemiLRSignalSelObservables();  
  
  void  operator()(TtSemiEvtSolution&, const std::vector<pat::Jet>&);
  
 private:
  
  std::vector<std::pair<unsigned int,double> > evtselectVarVal;
  
  // compare two jets in ET
  struct CompareET {
    bool operator()( pat::Jet j1, pat::Jet j2 ) const
    {
      return j1.et() > j2.et();
    }
  };
  
  CompareET EtComparator;
  
  // compare two jets in bdisc
  struct CompareBdisc {
    bool operator()( pat::Jet j1, 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;
};

#endif