LambdaCSelector.cxx

Go to the documentation of this file.

#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/IJobOptionsSvc.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/PropertyMgr.h"
 
#include "CLHEP/Vector/LorentzVector.h"
#include "DTagAlg/LambdaCSelector.h"
 
 
LambdaCSelector::LambdaCSelector( ){
 
        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);
}
 
bool LambdaCSelector::operator() (CDDecay& theD) {
 
        /// 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;
}
 
 
LambdaCSelector lambdaCSelector;