L1AnalysisRecoMet.cc

Go to the documentation of this file.

#include "L1Trigger/L1TNtuples/interface/L1AnalysisRecoMet.h"
#include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"

#include <TVector2.h>


L1Analysis::L1AnalysisRecoMet::L1AnalysisRecoMet() 
{
}

L1Analysis::L1AnalysisRecoMet::~L1AnalysisRecoMet()
{
}

void L1Analysis::L1AnalysisRecoMet::SetMet(const edm::Handle<reco::CaloMETCollection> recoMet)
{ 
  const reco::CaloMETCollection *metCol = recoMet.product();
  const reco::CaloMET theMet = metCol->front();

  recoMet_.met    = theMet.et();
  recoMet_.metPhi = theMet.phi();
  recoMet_.sumEt  = theMet.sumEt();

}

void L1Analysis::L1AnalysisRecoMet::SetHtMht(const edm::Handle<reco::CaloJetCollection> caloJets, float jetptThreshold)
{  
  float mHx = 0.;
  float mHy = 0.;

  recoMet_.Ht     = 0;
  recoMet_.mHt    = -999;
  recoMet_.mHtPhi = -999;
  
  for (reco::CaloJetCollection::const_iterator calojet = caloJets->begin(); calojet!=caloJets->end(); ++calojet)
  {
    if (calojet->pt()>jetptThreshold){
      mHx += -1.*calojet->px();
      mHy += -1.*calojet->py();
      recoMet_.Ht  += calojet->pt();
    }
  }

  TVector2 *tv2 = new TVector2(mHx,mHy);
  
  recoMet_.mHt  = tv2->Mod();
  recoMet_.mHtPhi= tv2->Phi();

}

void L1Analysis::L1AnalysisRecoMet::SetECALFlags(const edm::ESHandle<EcalChannelStatus> chStatus,
                         const edm::Handle<EcalRecHitCollection> ebRecHits,
                         const edm::Handle<EcalRecHitCollection> eeRecHits,
                                                 const EcalSeverityLevelAlgo* sevlv)
{
  int ecalFlag=0;
  
  // loop over EB rechits
  for(EcalRecHitCollection::const_iterator  rechit = ebRecHits->begin();
      rechit != ebRecHits->end();
      ++rechit){

    EBDetId eid(rechit->id());
    //rechit->recoFlag();
    //rechit->chi2();
    //rechit->outOfTimeChi2();

    //int flag = EcalSeverityLevelAlgo::severityLevel( eid, *ebRecHits, *chStatus );
    int flag = sevlv->severityLevel( eid, *ebRecHits);
    if (flag>ecalFlag) ecalFlag = flag;
  }
  
  // not clear what flags are in EE....  don't use them yet

  recoMet_.ecalFlag = ecalFlag;

}