Funs.h

Go to the documentation of this file.

#ifndef __FUNS_H
#define __FUNS_H

#include "jbltools/sfh/FillIterator.h"
#include "jbltools/sfh/BaseCut.h"
#include "jbltools/sfh/FloatFun.h"

#include <TClonesArray.h>
#include "Marana/TTools.h"

class NSelJetsFun: public FloatFun {
  public:
    NSelJetsFun (TClonesArray* jet_, 
                 const std::string& name_      
                ) 
      : FloatFun (name_), jet (jet_)
      {}
    virtual float operator() () const {
      return (jet < 0) ? 0 : float(jet->GetEntries());
    }
  protected:
    // Data members
    TClonesArray* jet;
    virtual ~NSelJetsFun() {}
};

class PtSelJet1Fun: public FloatFun {
  public:
    PtSelJet1Fun (TClonesArray* jet_, 
                  const std::string& name_      
                 ) 
      : FloatFun (name_), jet (jet_)
      {}
    virtual float operator() () const {
      float result = 0;
      if (jet) {
        for (int i = 0; i < jet->GetEntries(); i++) {
          TMarBody *Jet = (TMarBody*) jet->At(i);
          if (Jet->GetPt() > result) result = Jet->GetPt();
        }
      }
      return result;
    }
  protected:
    // Data members
    TClonesArray* jet;
    virtual ~PtSelJet1Fun() {}
};

class ThSelJet1Fun: public FloatFun {
  public:
    ThSelJet1Fun (TClonesArray* jet_, 
                  const std::string& name_      
                 ) 
      : FloatFun (name_), jet (jet_)
      {}
    virtual float operator() () const {
      float result = 0;
      if (jet) {
        float pt = 0;
        for (int i = 0; i < jet->GetEntries(); i++) {
          TMarBody *Jet = (TMarBody*) jet->At(i);
          if (Jet->GetPt() > pt) {
            pt = Jet->GetPt();
            result = Jet->GetTheta();
          }
        }
      }
      return result;
    }
  protected:
    // Data members
    TClonesArray* jet;
    virtual ~ThSelJet1Fun() {}
};

class TimimgCut: public BaseCut {
  public:
    TimimgCut (TTools* tools_, 
               const std::string& name_      
                 ) 
      : BaseCut (name_), tools (tools_)
      {}
    virtual bool operator() () const {
      return (tools) ? tools->Is2004CalibratedTimingOK() : false;
    }
  protected:
    // Data members
    TTools* tools;
    virtual ~TimimgCut() {}
};


#endif /* #ifndef __FUNS_H */ 

---