GlobalHaloData

Macros

Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97
#ifndef DATAFORMATS_METRECO_GLOBALHALODATA_H
#define DATAFORMATS_METRECO_GLOBALHALODATA_H
/*
  [class]:  GlobalHaloData
  [authors]: R. Remington, The University of Florida
  [description]: Container class to store global beam halo data synthesized from EcalHaloData, HcalHaloData, and CSCHaloData. Also stores some variables relevant to MET for possible corrections.
  [date]: October 15, 2009
*/
#include "DataFormats/METReco/interface/EcalHaloData.h"
#include "DataFormats/METReco/interface/HcalHaloData.h"
#include "DataFormats/METReco/interface/CSCHaloData.h"
#include "DataFormats/METReco/interface/PhiWedge.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/MET.h"

namespace reco {
  class GlobalHaloData {
  public:
    // Constructor
    GlobalHaloData();
    // Destructor
    ~GlobalHaloData() {}

    //A good cut-variable to isolate halo events with no overlapping physics from collisions
    float METOverSumEt() const { return METOverSumEt_; }

    //Correction to CaloMET x-component
    float DeltaMEx() const { return dMEx_; }

    //Correction to CaloMET y-component
    float DeltaMEy() const { return dMEy_; }

    //Correction to SumEt
    float DeltaSumEt() const { return dSumEt_; }

    //Get CaloMET Object corrected for BeamHalo
    reco::CaloMET GetCorrectedCaloMET(const reco::CaloMET& RawMET) const;

    std::vector<PhiWedge>& GetMatchedHcalPhiWedges() { return HcalPhiWedges; }
    const std::vector<PhiWedge>& GetMatchedHcalPhiWedges() const { return HcalPhiWedges; }

    std::vector<PhiWedge>& GetMatchedEcalPhiWedges() { return EcalPhiWedges; }
    const std::vector<PhiWedge>& GetMatchedEcalPhiWedges() const { return EcalPhiWedges; }

    edm::RefVector<EcalRecHitCollection>& GetEBRechits() { return ecalebrhRefs; }
    const edm::RefVector<EcalRecHitCollection>& GetEBRechits() const { return ecalebrhRefs; }

    edm::RefVector<EcalRecHitCollection>& GetEERechits() { return ecaleerhRefs; }
    const edm::RefVector<EcalRecHitCollection>& GetEERechits() const { return ecaleerhRefs; }

    edm::RefVector<HBHERecHitCollection>& GetHBHERechits() { return hbherhRefs; }
    const edm::RefVector<HBHERecHitCollection>& GetHBHERechits() const { return hbherhRefs; }

    bool GetSegmentIsHBCaloMatched() const { return segmentisHBcalomatched; }
    bool GetSegmentIsHECaloMatched() const { return segmentisHEcalomatched; }
    bool GetSegmentIsEBCaloMatched() const { return segmentisEBcalomatched; }
    bool GetSegmentIsEECaloMatched() const { return segmentisEEcalomatched; }

    bool GetHaloPatternFoundEB() const { return halopatternfoundEB; }
    bool GetHaloPatternFoundEE() const { return halopatternfoundEE; }
    bool GetHaloPatternFoundHB() const { return halopatternfoundHB; }
    bool GetHaloPatternFoundHE() const { return halopatternfoundHE; }

    //Setters
    void SetMETOverSumEt(float x) { METOverSumEt_ = x; }
    void SetMETCorrections(float x, float y) {
      dMEx_ = x;
      dMEy_ = y;
    }

    void SetSegmentIsHBCaloMatched(bool b) { segmentisHBcalomatched = b; }
    void SetSegmentIsHECaloMatched(bool b) { segmentisHEcalomatched = b; }
    void SetSegmentIsEBCaloMatched(bool b) { segmentisEBcalomatched = b; }
    void SetSegmentIsEECaloMatched(bool b) { segmentisEEcalomatched = b; }
    void SetHaloPatternFoundEB(bool b) { halopatternfoundEB = b; }
    void SetHaloPatternFoundEE(bool b) { halopatternfoundEE = b; }
    void SetHaloPatternFoundHB(bool b) { halopatternfoundHB = b; }
    void SetHaloPatternFoundHE(bool b) { halopatternfoundHE = b; }

  private:
    float METOverSumEt_;
    float dMEx_;
    float dMEy_;
    float dSumEt_;

    std::vector<PhiWedge> HcalPhiWedges;
    std::vector<PhiWedge> EcalPhiWedges;

    bool segmentisEBcalomatched, segmentisEEcalomatched, segmentisHBcalomatched, segmentisHEcalomatched;
    bool halopatternfoundEB, halopatternfoundEE, halopatternfoundHB, halopatternfoundHE;

    edm::RefVector<EcalRecHitCollection> ecalebrhRefs;
    edm::RefVector<EcalRecHitCollection> ecaleerhRefs;
    edm::RefVector<HBHERecHitCollection> hbherhRefs;
  };
}  // namespace reco
#endif