TrigEgammaMonitorTagAndProbeAlgorithm.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
#ifndef TrigEgammaMonitorTagAndProbeAlgorithm_H
#define TrigEgammaMonitorTagAndProbeAlgorithm_H
 
#include "TrigEgammaMonitorAnalysisAlgorithm.h"
#include "StoreGate/ReadDecorHandleKey.h"
 
 
 
 
class TrigEgammaMonitorTagAndProbeAlgorithm: public TrigEgammaMonitorAnalysisAlgorithm 
{
 
  public:
 
    TrigEgammaMonitorTagAndProbeAlgorithm( const std::string& name, ISvcLocator* pSvcLocator );
 
    virtual ~TrigEgammaMonitorTagAndProbeAlgorithm() override;
 
    virtual StatusCode initialize() override;
    
    virtual StatusCode fillHistograms( const EventContext& ctx) const override;
 
 
  protected:
 
    bool executeTandP( const EventContext& ctx, std::vector<std::shared_ptr<const xAOD::Electron>> & ) const;
 
    void matchObjects(const std::string& trigItem,  std::vector<std::shared_ptr<const xAOD::Electron>>&, 
                      std::vector<std::pair<const xAOD::Egamma*, const TrigCompositeUtils::Decision*>> &) const; 
    
   
    
 
    std::vector<std::string> m_trigList;
    
 
  private:
 
    bool minimalTriggerRequirement () const;
    bool isTagElectron( const ToolHandle<GenericMonitoringTool>& monGroup, const xAOD::Electron *el) const;
    bool isGoodProbeElectron( const ToolHandle<GenericMonitoringTool>& monGroup, const xAOD::Electron *el, const xAOD::JetContainer *) const;
    float getPseudoLifetime(const xAOD::Electron*,const xAOD::Electron*) const;
    double simple_lxy(int ,double, double, double , double , double , double , double, double ) const;
    void dressPid(const xAOD::Electron *eg) const;
    
 
    Gaudi::Property<std::vector<std::string>> m_trigInputList{this, "TriggerList", {}};
    Gaudi::Property<std::vector<std::string>> m_tagTrigList{ this, "TagTriggerList", {}};
    Gaudi::Property<float> m_ZeeMassMin{ this, "ZeeLowerMass", 80};
    Gaudi::Property<float> m_ZeeMassMax{ this, "ZeeUpperMass", 100};
    Gaudi::Property<std::string> m_offTagTightness{ this, "OfflineTagSelector", "Tight"};
    Gaudi::Property<std::string> m_offProbeTightness{ this, "OfflineProbeSelector", "Loose"};
    Gaudi::Property<bool> m_oppositeCharge{ this, "OppositeCharge", true};
    Gaudi::Property<float> m_tagMinEt{ this, "OfflineTagMinEt", 25};
    Gaudi::Property<float> m_probeMinEt{this, "OfflineProbeMinEt", 4};
    Gaudi::Property<std::string> m_offProbeIsolation{ this, "OfflineProbeIsolation", "Loose"};
    Gaudi::Property<bool> m_rmCrack{this, "RemoveCrack", true};
    Gaudi::Property<bool> m_applyMinimalTrigger{this, "ApplyMinimalTrigger", true};
    Gaudi::Property<bool> m_applyJetNearProbeSelection{this, "ApplyJetNearProbeSelection", true};
    Gaudi::Property<bool> m_doJpsiee{this,"DoJpsiee", false};
    Gaudi::Property<std::string> m_anatype{ this, "Analysis","Zee"};
    SG::ReadHandleKey<xAOD::ElectronContainer> m_offElectronKey{ this, "ElectronKey", "Electrons", ""};
    SG::ReadHandleKey<xAOD::JetContainer> m_jetKey{ this, "JetKey" , "AntiKt4LCTopoJets", ""};
    SG::ReadDecorHandleKey<xAOD::EventInfo> m_eventInfoDecorKey{this, "LArStatusFlag", "EventInfo.larFlags", "Key for EventInfo object"};
    SG::ReadDecorHandleKey<xAOD::ElectronContainer> m_electronIsolationKeyPtCone20 { this, "PtCone20Decoration", m_offElectronKey, "ptcone20", "Decoration key for the ptcone20 isolation decoration" };
    
};
 
#endif