#include <NeutralDSelector.h>

Inheritance diagram for NeutralDSelector:

Public Member Functions
	NeutralDSelector ()

void	setebeam (double ebeam)

void	setbeta (Hep3Vector beta)

bool	operator() (CDDecay &theD)

Public Member Functions inherited from DCSelectionFunction< CDDecay >
	DCSelectionFunction ()

virtual	~DCSelectionFunction ()

virtual bool	operator() (CDDecay &)=0

bool	operator() (CDDecay &iArg) const

Detailed Description

Definition at line 7 of file NeutralDSelector.h.

Constructor & Destructor Documentation

◆ NeutralDSelector()

NeutralDSelector::NeutralDSelector ( )

Definition at line 10 of file NeutralDSelector.cxx.

                                   {
  
  IJobOptionsSvc* jobSvc;
  Gaudi::svcLocator()->service("JobOptionsSvc", jobSvc);
  
  PropertyMgr m_propMgr;
  
  m_propMgr.declareProperty("UseMbcCuts",        m_useMbcCuts       = true );
  m_propMgr.declareProperty("UseDeltaECuts",     m_useDeltaECuts    = true );
  m_propMgr.declareProperty("UseDeltaMassCuts",  m_useDeltaMassCuts = false );
 
  m_propMgr.declareProperty("DMbcMinCut",        m_minMbc  = 1.83 );
  m_propMgr.declareProperty("DMbcMaxCut",        m_maxMbc  = 10.0 );
  m_propMgr.declareProperty("DDeltaEMinCut",     m_minDeltaE = -0.1 );
  m_propMgr.declareProperty("DDeltaEMaxCut",     m_maxDeltaE = +0.1 );
  m_propMgr.declareProperty("DDeltaMassMinCut",  m_minDeltaMass = -0.085 );
  m_propMgr.declareProperty("DDeltaMassMaxCut",  m_maxDeltaMass = +0.085 );
 
  jobSvc->setMyProperties("NeutralDSelector", &m_propMgr);
}

Member Function Documentation

◆ operator()()

bool NeutralDSelector::operator() ( CDDecay & theD )

virtual

Calculating mbc and deltaE

Calculating invariant mass difference

Implements DCSelectionFunction< CDDecay >.

Definition at line 31 of file NeutralDSelector.cxx.

                                                {
  
  /// Calculating mbc and deltaE
  HepLorentzVector p4( theD.momentum(), theD.energy() );
  p4.boost(-m_beta);
  double mbc2 = m_ebeam*m_ebeam - p4.v().mag2();
  double mbc = mbc2 > 0 ? sqrt( mbc2 ) : -10;
  double deltaE = p4.t() - m_ebeam;
 
  if(m_useMbcCuts&&((mbc <= m_minMbc)||(mbc >= m_maxMbc))) return false;
  if(m_useDeltaECuts&&((deltaE <= m_minDeltaE)||(deltaE >= m_maxDeltaE))) return false;
 
  /// Calculating invariant mass difference
  double neutralDMass = 1.8648;
  double deltaMass = theD.mass() - neutralDMass;
 
  if(m_useDeltaMassCuts &&
     ((deltaMass <= m_minDeltaMass)||(deltaMass >= m_maxDeltaMass))) return false;
 
 
  return true;
}