DerivationFramework::Select_onia2mumu Class Reference

#include <Select_onia2mumu.h>

Inheritance diagram for DerivationFramework::Select_onia2mumu:

Collaboration diagram for DerivationFramework::Select_onia2mumu:

Public Member Functions
	Select_onia2mumu (const std::string &t, const std::string &n, const IInterface *p)
StatusCode	initialize () override
	inirialization and finalization
virtual StatusCode	addBranches (const EventContext &ctx) const override
	: augmentation and selection Retrieved vertices are augmented with usual information.

Private Member Functions
void	ProcessVertex (xAOD::BPhysHypoHelper &, xAOD::BPhysHelper::pv_type) const

Private Attributes
ToolHandle< Trk::V0Tools >	m_v0Tools
	tools
std::string	m_hypoName
	job options
SG::ReadHandleKey< xAOD::VertexContainer >	m_inputVtxContainerName
	name of the input container name
std::vector< double >	m_trkMasses
	track mass hypotheses
double	m_massHypo
	vertex mass hypothesis
double	m_massMax
	invariant mass range
double	m_massMin
	invariant mass range
double	m_chi2Max
	max chi2 cut
double	m_lxyMin
	min lxy cut
int	m_DoVertexType
	Allows user to skip certain vertexes - bitwise test 7==all(111)
bool	m_do3d

Detailed Description

Definition at line 32 of file Select_onia2mumu.h.

Constructor & Destructor Documentation

Select_onia2mumu()

DerivationFramework::Select_onia2mumu::Select_onia2mumu	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Definition at line 25 of file Select_onia2mumu.cxx.

                           : 
    base_class(t,n,p),
    m_v0Tools("Trk::V0Tools")
  {
 
    // Declare tools    
    declareProperty("V0Tools", m_v0Tools);
 
    // Declare user-defined properties
    
    declareProperty("HypothesisName"       , m_hypoName              = "A");
    declareProperty("InputVtxContainerName", m_inputVtxContainerName = "JpsiCandidates");
    declareProperty("TrkMasses"            , m_trkMasses             = std::vector<double>(2, ParticleConstants::muonMassInMeV) );    
    declareProperty("VtxMassHypo"          , m_massHypo              = ParticleConstants::JpsiMassInMeV );                  
    declareProperty("MassMax"              , m_massMax               = 6000);                   
    declareProperty("MassMin"              , m_massMin               = 2000);                   
    declareProperty("Chi2Max"              , m_chi2Max               = 200);
    declareProperty("DoVertexType"         , m_DoVertexType          = 7);
    declareProperty("LxyMin"               , m_lxyMin                = std::numeric_limits<double>::lowest());
    declareProperty("Do3d"                 , m_do3d = false);
    
  }

Member Function Documentation

addBranches()

StatusCode DerivationFramework::Select_onia2mumu::addBranches ( const EventContext & ctx ) const

overridevirtual

: augmentation and selection Retrieved vertices are augmented with usual information.

Selection is performed and each candidate is decorated with the Char_t flag named "passed_"+name() to indicate whether if the candidate passed the selection. This flag is then used by the event selection tool and by the vertex thinning tool.

Definition at line 131 of file Select_onia2mumu.cxx.

  {
    SG::ReadHandle<xAOD::VertexContainer> oniaContainer(m_inputVtxContainerName, ctx);
    SG::auxid_set_t decor_auxids;
 
    bool doPt   = (m_DoVertexType & 1) != 0;
    bool doA0   = (m_DoVertexType & 2) != 0;
    bool doZ0   = (m_DoVertexType & 4) != 0;
    bool doZ0BA = (m_DoVertexType & 8) != 0;
    // loop over onia candidates and perform selection and augmentation
    xAOD::VertexContainer::const_iterator oniaItr = oniaContainer->begin();
    for(; oniaItr!=oniaContainer->end(); ++oniaItr) {
      // create BPhysHypoHelper
      xAOD::BPhysHypoHelper onia(m_hypoName, *oniaItr, &decor_auxids);
      if((*oniaItr)->nTrackParticles() != m_trkMasses.size())
          ATH_MSG_WARNING("Vertex has " << (*oniaItr)->nTrackParticles() << " while provided masses " << m_trkMasses.size());
      //----------------------------------------------------
      // decorate the vertex
      //----------------------------------------------------
      // a) invariant mass and error
      if( !onia.setMass(m_trkMasses) ) ATH_MSG_WARNING("Decoration onia.setMass failed");
      
      double massErr = m_v0Tools->invariantMassError(onia.vtx(), m_trkMasses);
      if( !onia.setMassErr(massErr) ) ATH_MSG_WARNING("Decoration onia.setMassErr failed");
      
      // b) proper decay time and error: 
      // retrieve the refitted PV (or the original one, if the PV refitting was turned off)
      if(doPt)   ProcessVertex(onia, xAOD::BPhysHelper::PV_MAX_SUM_PT2);
      if(doA0)   ProcessVertex(onia, xAOD::BPhysHelper::PV_MIN_A0);
      if(doZ0)   ProcessVertex(onia, xAOD::BPhysHelper::PV_MIN_Z0);
      if(doZ0BA) ProcessVertex(onia, xAOD::BPhysHelper::PV_MIN_Z0_BA);
      
      //----------------------------------------------------
      // perform the selection (i.e. flag the vertex)
      //----------------------------------------------------
      // flag the vertex indicating that it is selected by this selector
      onia.setPass(true);
      
      // now we check othe cuts. if one of them didn't pass, set the flag to 0
      // and continue to the next candidate:
      
      // 1) invariant mass cut
      if( onia.mass() < m_massMin || onia.mass() > m_massMax) {
        onia.setPass(false); // flag as failed
        continue;
      }
 
      // 2) chi2 cut
      if( onia.vtx()->chiSquared() > m_chi2Max) {
        onia.setPass(false);; // flag as failed
        continue;
      }
      // 3) lxy cut
      if( onia.lxy(xAOD::BPhysHelper::PV_MAX_SUM_PT2) < m_lxyMin) {
        onia.setPass(false);; // flag as failed
        continue;
      }
 
    } // end of loop over onia candidates
 
    // Lock the decorations we just produced.
    xAOD::VertexContainer* onia_nc ATLAS_THREAD_SAFE =
      const_cast<xAOD::VertexContainer*> (oniaContainer.cptr());
    for (SG::auxid_t auxid : decor_auxids) {
      onia_nc->lockDecoration (auxid);
    }
    
    // all OK
    return StatusCode::SUCCESS;
  }  

initialize()

StatusCode DerivationFramework::Select_onia2mumu::initialize ( )

override

inirialization and finalization

Definition at line 52 of file Select_onia2mumu.cxx.

  {
  
    ATH_MSG_DEBUG("in initialize()");
    
    // retrieve V0 tools
    CHECK( m_v0Tools.retrieve() );
    ATH_CHECK(m_inputVtxContainerName.initialize());
 
    return StatusCode::SUCCESS;
    
  }

ProcessVertex()

void DerivationFramework::Select_onia2mumu::ProcessVertex ( xAOD::BPhysHypoHelper & onia, xAOD::BPhysHelper::pv_type pv_t ) const

private

Definition at line 67 of file Select_onia2mumu.cxx.

                                                                                                  {
      constexpr float errConst = -9999999;
      const xAOD::Vertex* pv = onia.pv(pv_t); 
      if(pv) {
        // decorate the vertex. 
        // Proper decay time assuming constant mass hypothesis m_massHypo
        BPHYS_CHECK( onia.setTau( m_v0Tools->tau(onia.vtx(), pv,  m_massHypo), 
                                  pv_t,
                                  xAOD::BPhysHypoHelper::TAU_CONST_MASS) );
        // Proper decay time assuming error constant mass hypothesis m_massHypo
        BPHYS_CHECK( onia.setTauErr( m_v0Tools->tauError(onia.vtx(), pv,  m_massHypo), 
                                     pv_t,
                                     xAOD::BPhysHypoHelper::TAU_CONST_MASS) );
        
        BPHYS_CHECK( onia.setTau( m_v0Tools->tau(onia.vtx(), pv,  m_trkMasses), 
                                  pv_t,
                                  xAOD::BPhysHypoHelper::TAU_INV_MASS) );
 
        BPHYS_CHECK( onia.setTauErr( m_v0Tools->tauError(onia.vtx(), pv,  m_trkMasses), 
                                     pv_t,
                                     xAOD::BPhysHypoHelper::TAU_INV_MASS) );       
        
        //enum pv_type {PV_MAX_SUM_PT2, PV_MIN_A0, PV_MIN_Z0, PV_MIN_Z0_BA};
      }else{
      
 
        BPHYS_CHECK( onia.setTau(errConst, pv_t,
                                  xAOD::BPhysHypoHelper::TAU_CONST_MASS) );
        // Proper decay time assuming error constant mass hypothesis m_massHypo
        BPHYS_CHECK( onia.setTauErr( errConst, 
                                     pv_t,
                                     xAOD::BPhysHypoHelper::TAU_CONST_MASS) );
        
        BPHYS_CHECK( onia.setTau( errConst, 
                                  pv_t,
                                  xAOD::BPhysHypoHelper::TAU_INV_MASS) );
 
        BPHYS_CHECK( onia.setTauErr( errConst, 
                                     pv_t,
                                     xAOD::BPhysHypoHelper::TAU_INV_MASS) );        
    }
 
    if(m_do3d){
        BPHYS_CHECK( onia.setTau3d( pv ? m_v0Tools->tau3D(onia.vtx(), pv,  m_massHypo) : errConst, 
                                  pv_t,
                                  xAOD::BPhysHypoHelper::TAU_CONST_MASS) );
        // Proper decay time assuming error constant mass hypothesis m_massHypo
        BPHYS_CHECK( onia.setTau3dErr( pv ? m_v0Tools->tau3DError(onia.vtx(), pv,  m_massHypo) : errConst,
                                     pv_t,
                                     xAOD::BPhysHypoHelper::TAU_CONST_MASS) );
        
        BPHYS_CHECK( onia.setTau3d( pv ? m_v0Tools->tau3D(onia.vtx(), pv,  m_trkMasses) : errConst,
                                  pv_t,
                                  xAOD::BPhysHypoHelper::TAU_INV_MASS) );
 
        BPHYS_CHECK( onia.setTau3dErr( pv ? m_v0Tools->tau3DError(onia.vtx(), pv,  m_trkMasses) : errConst,
                                     pv_t,
                                     xAOD::BPhysHypoHelper::TAU_INV_MASS) );    
 
    }
 
  }

Member Data Documentation

m_chi2Max

double DerivationFramework::Select_onia2mumu::m_chi2Max

private

max chi2 cut

Definition at line 63 of file Select_onia2mumu.h.

m_do3d

bool DerivationFramework::Select_onia2mumu::m_do3d

private

Definition at line 66 of file Select_onia2mumu.h.

m_DoVertexType

int DerivationFramework::Select_onia2mumu::m_DoVertexType

private

Allows user to skip certain vertexes - bitwise test 7==all(111)

Definition at line 65 of file Select_onia2mumu.h.

m_hypoName

std::string DerivationFramework::Select_onia2mumu::m_hypoName

private

job options

name of the mass hypothesis. E.g. Jpis, Upsi, etc. Will be used as a prefix for decorations

Definition at line 57 of file Select_onia2mumu.h.

m_inputVtxContainerName

SG::ReadHandleKey<xAOD::VertexContainer> DerivationFramework::Select_onia2mumu::m_inputVtxContainerName

private

name of the input container name

Definition at line 58 of file Select_onia2mumu.h.

m_lxyMin

double DerivationFramework::Select_onia2mumu::m_lxyMin

private

min lxy cut

Definition at line 64 of file Select_onia2mumu.h.

m_massHypo

double DerivationFramework::Select_onia2mumu::m_massHypo

private

vertex mass hypothesis

Definition at line 60 of file Select_onia2mumu.h.

m_massMax

double DerivationFramework::Select_onia2mumu::m_massMax

private

invariant mass range

Definition at line 61 of file Select_onia2mumu.h.

m_massMin

double DerivationFramework::Select_onia2mumu::m_massMin

private

invariant mass range

Definition at line 62 of file Select_onia2mumu.h.

m_trkMasses

std::vector<double> DerivationFramework::Select_onia2mumu::m_trkMasses

private

track mass hypotheses

Definition at line 59 of file Select_onia2mumu.h.

m_v0Tools

ToolHandle<Trk::V0Tools> DerivationFramework::Select_onia2mumu::m_v0Tools

private

tools

Definition at line 53 of file Select_onia2mumu.h.

---

The documentation for this class was generated from the following files:

• Select_onia2mumu.h
• Select_onia2mumu.cxx