MuonTag.cc

Go to the documentation of this file.

#include "MuonTag.h"
#include "LafKernel/PhyEvent/RecRpcReadout.h"
#include "LafKernel/PhyEvent.h"
#include "LafKernel/PhyEvent/RecReadout.h"
#include "LafKernel/ToolFactory.h"

DECLARE_TOOL(MuonTag);

MuonTag::MuonTag(const std::string& name)
    : ITagBase(name)
{
  setOption("ApplyVeto",m_veto=false);
  setOption("AdMuonThreshold",m_adMuonThr);
  setOption("ShowerMuonThreshold",m_showerMuonThr);
  setOption("IwsNPmtThreshold",m_iwsNpmtThr);
  setOption("OwsNPmtThreshold",m_owsNpmtThr);

  OptionParser::setOption(name, "tagRPCMuon", m_tagRPCMuon=true);
  OptionParser::setOption(name, "tagADMuon", m_tagADMuon=true);
}

bool MuonTag::tag(PhyEvent *event)
{
  if(!event->m_good) {
    return false;
  }

  if(event->m_det==5) { // IWS
    if(event->m_nPmt>m_iwsNpmtThr) {
      event->m_muonTag=PhyEvent::kWpMuon;
      return m_veto;
    }
  } else if (event->m_det==6) { // OWS
    if(event->m_nPmt>m_owsNpmtThr) {
      event->m_muonTag=PhyEvent::kWpMuon;
      return m_veto;
    }
  } else if (event->m_det==7) {
    if(m_tagRPCMuon && ((event->m_trigType&0x00040000) != 0) ) {  // RPC 3/4
      RecRpcReadout* recRpc = event->get<RecRpcReadout>("Event/Rec/RpcSimple");
      if(recRpc->eventType == 3 || recRpc->maxHitsInLayer() < 3)  
            event->m_muonTag=PhyEvent::kRpcMuon;
      return m_veto;
    }
  } else {
    if(event->energy()>m_adMuonThr && m_tagADMuon) {
      event->m_muonTag=PhyEvent::kAdMuon;
      if(event->energy()>m_showerMuonThr) {
        event->m_muonTag=PhyEvent::kShowerMuon;
      }
      return m_veto;
    }
  }
  return false;
}