ChargedDSelector.cxx

Go to the documentation of this file.

#include "GaudiKernel/Bootstrap.h"
 
#include "ChargedDSelector.h"
 
DECLARE_COMPONENT( ChargedDSelector )
 
ChargedDSelector::ChargedDSelector( const std::string& type, const std::string& name,
                                    const IInterface* parent )
    : base_class( type, name, parent ) {
  declareProperty( "UseMbcCuts", m_useMbcCuts = true );
  declareProperty( "UseDeltaECuts", m_useDeltaECuts = true );
  declareProperty( "UseDeltaMassCuts", m_useDeltaMassCuts = false );
 
  declareProperty( "DMbcMinCut", m_minMbc = 1.83 );
  declareProperty( "DMbcMaxCut", m_maxMbc = 10.0 );
  declareProperty( "DDeltaEMinCut", m_minDeltaE = -0.1 );
  declareProperty( "DDeltaEMaxCut", m_maxDeltaE = +0.1 );
  declareProperty( "DDeltaMassMinCut", m_minDeltaMass = -0.085 );
  declareProperty( "DDeltaMassMaxCut", m_maxDeltaMass = +0.085 );
}
 
bool ChargedDSelector::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 = 1.8696;
  double deltaMass    = theD.mass() - chargedDMass;
 
  if ( m_useDeltaMassCuts &&
       ( ( deltaMass <= m_minDeltaMass ) || ( deltaMass >= m_maxDeltaMass ) ) )
    return false;
 
  return true;
}