---

TrigObjects.h

Go to the documentation of this file.

//================================================
//
// A set of function objects used to seed the trigger word bits
//
// ===============================================
//
#ifndef ATLFAST_TRIGOBJECTS_H
#define ATLFAST_TRIGOBJECTS_H

// Gaudi includes
#include "GaudiKernel/SmartDataPtr.h"

#include <cmath>
#include <cstdlib>


//--------------------------------------------------------
//
// Function objects for sort
//

namespace TrigObj {
  const float etaMaxElec    = xxx.x;
  const float etaMaxPhoton  = xxx.x;
  const float etaMaxMuon    = xxx.x;
  const float etaMaxEMu_E   = xxx.x; //eta covrge El-Mu diff from e or mu alone
  const float etaMaxEMu_M   = xxx.x; //lumi dep in ATL++, but values degenerate
  const float ptMinEM1      = xxx.x;
  const float ptMinEM2Elec  = xxx.x;
  const float ptMinEM2Phot  = xxx.x;
  const float ptMinMU1Isol  = xxx.x;
  const float ptMinMU2Isol  = xxx.x;
  const float ptMinMU1NIsol = xxx.x;
  const float ptMinMU2NIsol = xxx.x;
  const float ptMinEMUElec  = xxx.x;
  const float ptMinEMUMuon  = xxx.x;

  class IsolatedElectron {
  public:
    IsolatedElectron() m_em1(0), m_em2(0), m_emu(0){}
    bool operator() ( ReconstructedParticle* a) { 
      if( abs(a->eta()) <= etaMaxElec) return false;
      float pt=a->pt();
      if( pt> ptMinEM1) ++m_em1;
      if( pt> ptMinEM2) ++m_em2;
      if( pt> ptMinEMU) ++m_emu;
      return true;                        // return value not used!
    }
    bool em1(){return m_em1>1;}
    bool em2(){return m_em2;}
    int  emu(){return m_emu;}
  private:
    bool m_em1;
    int m_em2;
    int m_emu;
  };
  
  class IsolatedPhoton {
  public:
    Photon() m_ph1(0), m_em2(0) {}
    bool operator() ( ReconstructedParticle* a) { 
      if( abs(a->eta()) <= etaMaxPhoton) return false;
      float pt=a->pt();
      if( pt> ptMinPH1) ++m_ph1;
      if( pt> ptMinEM2) ++m_em2;
      return true;                        // return value not used!
    }
    bool em1(){return m_ph1>1;}
    int  em2(){return m_em2;}
  private:
    int m_ph1;
    int m_em2;
  };
//
  class IsolatedMuon {
  public:
    IsolatedMuon(int lumi) m_lumi(lumi), m_mu1(0), m_mu2(0) {}
    bool operator() ( ReconstructedParticle* a) { 
      if( abs(a->eta()) <= etaMaxMuon) return;
      float pt=a->pt();
      if(               pt> ptMinMU1Isol) ++m_mu1;
      if( highLumi() && pt> ptMinMU2Isol) ++m_mu2;
      return true;                        // return value not used!
    }
    bool mu1(){return m_mu1>1;}
    bool mu2(){return m_mu2>1;}
  private:
    bool highLumi(){return m_lumi==2;}
    int m_lumi;
    int m_mu1;
    int m_mu2;
  };
  class NonIsolatedMuonTrig {
  public:
    NonIsolatedMuon(int lumi) m_lumi(lumi), m_mu1(0), m_mu2(0) {}
    bool operator() ( ReconstructedParticle* a) { 
      etaMuon=abs(a->eta())
        if( etaMuon <= etaMaxMuon){
          float pt=a->pt();
          if(               pt> ptMinMU1NIsol) ++m_mu1;
          if( highlumi() && pt> ptMinMU2NIsol) ++m_mu2;
          return true;                        // return value not used!
        }
      if(etaMuon<=etaMaxEMu_M) ++m_emu; 
        
        
      }
      
      bool mu1(){return m_mu1>1;}
      bool mu2(){return m_mu2>2;}
    private:
      bool highLumi(){return m_lumi==2;}
      int m_lumi;
      int m_mu1;
      int m_mu2;
    };
}

---