#include <LambdaCSelector.h>

Inheritance diagram for LambdaCSelector:

Public Member Functions
	LambdaCSelector ()

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 LambdaCSelector.h.

Constructor & Destructor Documentation

◆ LambdaCSelector()

LambdaCSelector::LambdaCSelector ( )

Definition at line 10 of file LambdaCSelector.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("MbcMinCut",        m_minMbc  = 2.2 );
        m_propMgr.declareProperty("MbcMaxCut",        m_maxMbc  = 10.0 );
        m_propMgr.declareProperty("DeltaEMinCut",     m_minDeltaE = -0.15 );
        m_propMgr.declareProperty("DeltaEMaxCut",     m_maxDeltaE = +0.15 );
        m_propMgr.declareProperty("DeltaMassMinCut",  m_minDeltaMass = -0.15 );
        m_propMgr.declareProperty("DeltaMassMaxCut",  m_maxDeltaMass = +0.15 );
 
        jobSvc->setMyProperties("LambdaCSelector", &m_propMgr);
}

Member Function Documentation

◆ operator()()

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

virtual

Calculating mbc and deltaE

Calculating invariant mass difference

Implements DCSelectionFunction< CDDecay >.

Definition at line 31 of file LambdaCSelector.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 chargedDMass = 2.28646;
        double deltaMass = theD.mass() - chargedDMass;
 
        if(m_useDeltaMassCuts && ((deltaMass <= m_minDeltaMass)||(deltaMass >= m_maxDeltaMass))) return false;
 
 
        return true;
}