LocalEptoPiPiEta3PiSelector.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 "EventModel/EventModel.h"
#include "EventModel/EventHeader.h"
#include "EvtRecEvent/EvtRecEvent.h"
#include "EvtRecEvent/EvtRecTrack.h"
#include "ParticleID/ParticleID.h"
#include "DTagAlg/LocalEptoPiPiEta3PiSelector.h"
 
LocalEptoPiPiEta3PiSelector::LocalEptoPiPiEta3PiSelector()
{
  IJobOptionsSvc* jobSvc;
  Gaudi::svcLocator()->service("JobOptionsSvc", jobSvc);
 
  PropertyMgr m_propMgr;
 
  //Declare the properties
  m_propMgr.declareProperty("minMassEptoPiPiEta3Pi",                            m_minMass = 0.938);
  m_propMgr.declareProperty("maxMassEptoPiPiEta3Pi",                                m_maxMass = 0.978);
   
   
  jobSvc->setMyProperties("LocalEptoPiPiEta3PiSelector", &m_propMgr);
}
 
bool LocalEptoPiPiEta3PiSelector::operator() (CDDecay& aEp) {
 
  //aEp.setUserTag(1);
 
  int userTag_pion1=aEp.decay().child(0).userTag();
  int userTag_pion2=aEp.decay().child(1).userTag();
  int userTag_Eta=aEp.decay().child(2).userTag();
 
  if(userTag_pion1==1 && userTag_pion2==1 && userTag_Eta==1 ) aEp.setUserTag(1);
  else aEp.setUserTag(2);
 
  double mass = aEp.mass();
  if(mass >= m_minMass && mass<= m_maxMass)
    return true;
  else
    return false;
}
 
 
LocalEptoPiPiEta3PiSelector eptoPiPiEta3PiSelector;