SusyNtuple/JetSelector.h

//  -*- c++ -*-
#ifndef SUSY_JETSELECTOR_H
#define SUSY_JETSELECTOR_H

#include "SusyNtuple/SusyNtSys.h"
#include "SusyNtuple/AnalysisType.h"
#include "SusyNtuple/SusyDefs.h" // JetVector

namespace Susy {

class Jet;


class JetSelector
{
public:

    static JetSelector* build(const AnalysisType &a, bool verbose);

    JetSelector(); 
    virtual ~JetSelector() {}; 
    JetSelector& setSystematic(const NtSys::SusyNtSys&); 

    virtual bool isB(const Jet* jet);
    virtual bool isCentralLight(const Jet* jet);
    virtual bool isCentralB(const Jet* jet);
    virtual bool isForward(const Jet* jet);
    virtual bool passJvt(const Jet* jet);

    virtual bool isBaseline(const Jet* jet); 
    virtual bool isSignal(const Jet* jet); 

    bool verbose() const { return m_verbose; }
    JetSelector& setVerbose(bool v) { m_verbose = v; return *this; }
    virtual size_t count_CL_jets(const JetVector &jets) /*const*/;
    virtual size_t count_CB_jets(const JetVector &jets) /*const*/;
    virtual size_t count_F_jets(const JetVector &jets) /*const*/;


    static float mv2c20_70efficiency() { return -0.0436; } 
    static float mv2c20_77efficiency() { return -0.4434; } 
    static float mv2c20_80efficiency() { return -0.5911; } 


    static float overlapRemovalBtagEffWP() { return mv2c20_80efficiency(); }
private:
    NtSys::SusyNtSys m_systematic; 
    bool m_verbose; 
    }; // class JetSelector


//----------------------------------------------------------
//
// End generic selector, begin analysis-specific ones
//
//----------------------------------------------------------

class JetSelector_2Lep : public JetSelector
{
};
class JetSelector_3Lep : public JetSelector
{
};

class JetSelector_2LepWH : public JetSelector
{
};

class JetSelector_SS3L : public JetSelector
{
    virtual bool isBaseline(const Jet* jet);
    virtual bool isSignal(const Jet* jet);
    virtual bool isB(const Jet* jet);
};

class JetSelector_Stop2L : public JetSelector
{
    virtual bool isBaseline(const Jet* jet);
    virtual bool isSignal(const Jet* jet);
};

} // Susy

#endif