MpTrigger.cc

Go to the documentation of this file.

#include "MuonProphet.h"

/*
 * In this version, I am going to skip all muon geant4 simulation.
 * Have a quick check with geant4 for neutron and heavy isotopes only.
 *
 * Zhe Wang, 3 Dec. 2009
 */

// Triggered by RPC?
MpTriggerStat MuonProphet::triggeredByRpc(HepMC::GenParticle * /*pMuon*/)
{
  /* old one
  if(m_crossRpc) {
    return MpTriggerStat::kTriggered;
  } else {
    return MpTriggerStat::kNeedSim;
  }
  */
  
  if(m_crossRpc) {
    return MpTriggerStat::kTriggered;
  } else {
    return MpTriggerStat::kFarAway;
  }
  //
}

// Triggered by water pool?
MpTriggerStat MuonProphet::triggeredByWaterPool(HepMC::GenParticle * /*pMuon*/)
{
  /* old one
  if(m_crossWater) {
    if(m_trkLengthInWater > m_trkLengthInWaterThres) {
      return MpTriggerStat::kTriggered;
    }
  }

  return MpTriggerStat::kNeedSim;
  */

  if(m_crossWater) {
    if(m_trkLengthInWater > m_trkLengthInWaterThres) {    // 100% triggered
      return MpTriggerStat::kTriggered;
    } else {
      if(m_rnd.Rndm()<0.5) {   // 50% - 50%  triggered or inefficent
        return MpTriggerStat::kTriggered;
      } else {
        return MpTriggerStat::kIneffi;
      }
    } 
  } else {
    if(m_rnd.Rndm()<0.05) {   // 5% triggered  95% faraway
      return MpTriggerStat::kTriggered;
    } else {
      return MpTriggerStat::kFarAway;
    }
  }
}