TrigBphys_v1.h

Go to the documentation of this file.

// Dear emacs, this is -*- c++ -*-
 
/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
// $Id: TrigBphys_v1.h 696633 2015-09-24 14:59:50Z jwalder $
#ifndef XAODTRIGBPHYS_VERSIONS_TRIGBPHYS_V1_H
#define XAODTRIGBPHYS_VERSIONS_TRIGBPHYS_V1_H
 
// System include(s):
extern "C" {
#   include <stdint.h>
}
#include <vector>
#include <string>
 
// xAOD include(s):
#include "AthContainers/AuxElement.h"
#include "xAODTracking/TrackParticle.h"
#include "xAODTracking/TrackParticleContainer.h"
 
// athena includes
#include "AthLinks/ElementLink.h"
#include "AthContainers/DataVector.h"
 
namespace xAOD {
 
    // forward declarations, so that we can use TrigBphysContainer_v1 in TrigBphys_v1
    class TrigBphys_v1;
    typedef DataVector< xAOD::TrigBphys_v1 > TrigBphysContainer_v1;

    class TrigBphys_v1 : public SG::AuxElement {
 
    public:
        enum pType {
           PHIKK     = 0,
           DSPHIPI   = 1,
           BMUMU     = 2,
           BMUMUX    = 3,
           JPSIEE    = 4,
           JPSIMUMU  = 5,
           MULTIMU   = 6,
           BKMUMU    = 7,
           BDKSTMUMU = 8,
           BSPHIMUMU = 9,
           LBLMUMU   = 10,
           BCDSMUMU  = 11,
           // added JW: 21/09/15 for 2016 edm running
           LAMBDAPIP = 12, // Lambda^0 -> pi- p
           BCDPMUMU  = 13, // B_c -> D+ mu mu
           BCDSTMUMU = 14, // B_c -> D*+ mu mu
           OMEGALAMK = 15, // Omega- -> Lambda^0 K-
           XILAMPI   = 16, // Xi- -> Lambda^0 pi-
           DSTDZPI   = 17, // D*+ -> D0 pi+
           DZMUKX    = 18, // D0 -> mu+ K- X
           DZKPI     = 19, // D0 -> K- pi+
           DPKPIPI   = 20, // D+ -> K- pi+ pi+
           BCPIMUMU  = 21, // B_c+ -> pi+ mu mu
           LBPQMUMU  = 22, // Lambda_b0 -> P_c+(-> J/psi p) K-, pentaquark states
           B2D0MUX   = 23, // B+ -> mu+ nu_mu anti-D0(-> K+ pi-)
           BD2DMMUX  = 24, // B0 -> mu+ nu_mu D-(->K+ pi- pi-)
           BD2DSTMUX = 25, // B0 -> mu+ nu_mu D*-(->anti-D0(->K+ pi-) pi-)
           BS2DSMUX  = 26, // B_s0 -> mu+ nu_mu D_s-(->phi(->K+ K-) pi-)
           LB2LCMUX  = 27, // anti-Lambda_b0 -> mu+ nu_mu anti-Lambda_c-(-> anti-p K+ pi-)
           BHH       = 28, // B0 or B_s0 -> h+ h-, h = pi, K
           UNKNOWNPTYPE = 9999
        };

        enum levelType{
           UNKOWNLEVEL = 0,
           L2          = 2,
           EF          = 3,
           HLT         = 4,
           LASTLEVEL
        };

        TrigBphys_v1();
 
        //**clone method, essentially a replacement for a copy constructor*/
        void clone(const TrigBphys_v1 &rhs);
        
        //** initialise with defaults. Attempt to call this after creating a particle */
        void initialise();
 
        //** initialise with given parameters */
        void initialise(uint32_t roi, float eta, float phi, float pt,
                        pType particleType,
                        float mass, levelType level);
        void initialise(uint32_t roi, float eta, float phi, pType particleType,
                        float mass, levelType level);
        void initialise(uint32_t roi, float eta, float phi, float pt,
                        pType particleType,
                        float mass, const ElementLink<xAOD::TrigBphysContainer_v1>& sDecay,
                        levelType level);
        void initialise(uint32_t roi, float eta, float phi, pType particleType,
                        float mass, const ElementLink<xAOD::TrigBphysContainer_v1>& sDecay,
                        levelType level);
 
        bool operator==(const TrigBphys_v1& rhs) const;
        bool operator!=(const TrigBphys_v1& rhs) const {return !((*this)==rhs);}

        uint32_t roiId()    const;
        pType particleType() const;
        levelType level() const;
        float eta()      const;
        float phi()      const;
        float pt()      const;
        float mass()      const;
        float fitmass()  const;
        float fitchi2()  const;
        int fitndof()  const;
        float fitx()  const;
        float fity()  const;
        float fitz()  const;

        float lxy     () const;
        float lxyError() const;
        float tau     () const;
        float tauError() const;
        float fitmassError() const;
        
        const TrigBphys_v1* secondaryDecay() const;
        const ElementLink< TrigBphysContainer_v1 >& secondaryDecayLink() const;

        const std::vector< ElementLink< xAOD::TrackParticleContainer > >&
        trackParticleLinks() const;
        void setTrackParticleLinks(
           const std::vector< ElementLink< TrackParticleContainer > >& links );
        size_t nTrackParticles() const;
        const xAOD::TrackParticle* trackParticle( size_t i ) const;

        void setEta (float);
        void setPhi (float);
        void setMass(float);
        void setPt(float);

        void setRoiId(uint32_t id);
        void setParticleType(pType type);
        void setLevel(levelType type);
        void setFitmass(float FitMass);
        void setFitchi2(float FitChi2);
        void setFitndof(int FitNdof);
        void setFitx(float FitX);
        void setFity(float FitY);
        void setFitz(float FitZ);
        void addTrackParticleLink(const ElementLink<xAOD::TrackParticleContainer>& track);
        void setSecondaryDecayLink(const ElementLink<xAOD::TrigBphysContainer_v1> & link);

        void setLxy     (float v);
        void setLxyError(float v);
        void setTau     (float v);
        void setTauError(float v);
        void setFitmassError(float v);
 
        
        const std::vector< uint32_t > & vecRoiIds() const;
        void setVecRoiIds(const std::vector< uint32_t >& roiIds );
        size_t nVecRoiIds() const;
        uint32_t vecRoiId( size_t i ) const;
        void addVecRoiId(uint32_t roiId);
 
        
        
        const TrigBphys_v1* lowerChain() const;
        const ElementLink< TrigBphysContainer_v1 >& lowerChainLink() const;
        void setLowerChainLink(const ElementLink<xAOD::TrigBphysContainer_v1> & link);
        
        const std::vector< ElementLink< xAOD::IParticleContainer > >&
        particleLinks() const;
        void setParticleLinks(const std::vector< ElementLink< IParticleContainer > >& links );
        size_t nParticles() const;
        const xAOD::IParticle* particle( size_t i ) const;
        void addParticleLink(const ElementLink<xAOD::IParticleContainer>& particle);
        
        TrigBphys_v1& operator=(const TrigBphys_v1 & rhs) = default;
    private:
        //** Copy constructor */
        TrigBphys_v1(const TrigBphys_v1 & rhs) = delete;
 
    }; // TrigBphys
 
} // namespace AOD
 
#endif // XAODTRIGBPHYS_VERSIONS_TRIGBPHYS_V1_H