L1AnalysisL1Extra.cc

Go to the documentation of this file.

#include "L1Trigger/L1TNtuples/interface/L1AnalysisL1Extra.h"

L1Analysis::L1AnalysisL1Extra::L1AnalysisL1Extra()
{
}

L1Analysis::L1AnalysisL1Extra::~L1AnalysisL1Extra()
{

}

void L1Analysis::L1AnalysisL1Extra::SetIsoEm(const edm::Handle<l1extra::L1EmParticleCollection> isoEm, unsigned maxL1Extra)
{
  for(l1extra::L1EmParticleCollection::const_iterator it=isoEm->begin(); it!=isoEm->end() && l1extra_.nIsoEm<maxL1Extra; it++){
      
      l1extra_.isoEmEt .push_back(it->et());
      l1extra_.isoEmEta.push_back(it->eta());
      l1extra_.isoEmPhi.push_back(it->phi());
      l1extra_.isoEmBx .push_back(it->bx());
      l1extra_.nIsoEm++;
    }
}


void L1Analysis::L1AnalysisL1Extra::SetNonIsoEm(const edm::Handle<l1extra::L1EmParticleCollection> nonIsoEm, unsigned maxL1Extra)
{
  for(l1extra::L1EmParticleCollection::const_iterator it=nonIsoEm->begin(); it!=nonIsoEm->end() && l1extra_.nNonIsoEm<maxL1Extra; it++){
      
      l1extra_.nonIsoEmEt .push_back(it->et());
      l1extra_.nonIsoEmEta.push_back(it->eta());
      l1extra_.nonIsoEmPhi.push_back(it->phi());
      l1extra_.nonIsoEmBx .push_back(it->bx());
      l1extra_.nNonIsoEm++;
    }
}

void L1Analysis::L1AnalysisL1Extra::SetCenJet(const edm::Handle<l1extra::L1JetParticleCollection> cenJet, unsigned maxL1Extra)
{
//      std::cout << "Filling L1 Extra cenJets" << maxL1Extra << " " << cenJet->size() << std::endl;      
 
  for(l1extra::L1JetParticleCollection::const_iterator it=cenJet->begin(); it!=cenJet->end() && l1extra_.nCenJets<maxL1Extra; it++){
      //printf("L1CenJet (et,eta,phi,bx,) (%f,%f,%f,%d) \n",it->et(),it->eta(),it->phi(),it->bx() );
//      std::cout << "L1 CenJets et,eta,phi,bx = " << it->et() << ", " << it->eta() <<", " <<it->phi() <<", " << it->bx() << std::endl;
      l1extra_.cenJetEt .push_back(it->et());
      l1extra_.cenJetEta.push_back(it->eta());
      l1extra_.cenJetPhi.push_back(it->phi());
      l1extra_.cenJetBx .push_back(it->bx());
      l1extra_.nCenJets++;
    }
}

void L1Analysis::L1AnalysisL1Extra::SetFwdJet(const edm::Handle<l1extra::L1JetParticleCollection> fwdJet, unsigned maxL1Extra)
{
      //std::cout << "Filling L1 Extra fwdJets" << std::endl;      
   for(l1extra::L1JetParticleCollection::const_iterator it=fwdJet->begin(); it!=fwdJet->end() && l1extra_.nFwdJets<maxL1Extra; it++){
      
      //printf("L1fwdJet (et,eta,phi,bx,) (%f,%f,%f,%d)\n",it->et(),it->eta(),it->phi(),it->bx() );
      l1extra_.fwdJetEt .push_back(it->et());
      l1extra_.fwdJetEta.push_back(it->eta());
      l1extra_.fwdJetPhi.push_back(it->phi());
      l1extra_.fwdJetBx .push_back(it->bx());
      l1extra_.nFwdJets++;
    }
}

void L1Analysis::L1AnalysisL1Extra::SetTauJet(const edm::Handle<l1extra::L1JetParticleCollection> tauJet, unsigned maxL1Extra)
{
      //std::cout << "Filling L1 Extra tauJets" << std::endl;      
   for(l1extra::L1JetParticleCollection::const_iterator it=tauJet->begin(); it!=tauJet->end() && l1extra_.nTauJets<maxL1Extra; it++){
      
     // printf("L1tauJet (et,eta,phi,bx,) (%f,%f,%f,%d)\n",it->et(),it->eta(),it->phi(),it->bx() );
      l1extra_.tauJetEt .push_back(it->et());
      l1extra_.tauJetEta.push_back(it->eta());
      l1extra_.tauJetPhi.push_back(it->phi());
      l1extra_.tauJetBx .push_back(it->bx());
      l1extra_.nTauJets++;
    }
}

void L1Analysis::L1AnalysisL1Extra::SetIsoTauJet(const edm::Handle<l1extra::L1JetParticleCollection> isoTauJet, unsigned maxL1Extra)
{
  // std::cout << "Filling L1 Extra isoTauJets" << std::endl;      
   for(l1extra::L1JetParticleCollection::const_iterator it=isoTauJet->begin(); it!=isoTauJet->end() && l1extra_.nIsoTauJets<maxL1Extra; it++){
      
     // printf("L1isoTauJet (et,eta,phi,bx,) (%f,%f,%f,%d)\n",it->et(),it->eta(),it->phi(),it->bx() );
      l1extra_.isoTauJetEt .push_back(it->et());
      l1extra_.isoTauJetEta.push_back(it->eta());
      l1extra_.isoTauJetPhi.push_back(it->phi());
      l1extra_.isoTauJetBx .push_back(it->bx());
      l1extra_.nIsoTauJets++;
    }
}


void L1Analysis::L1AnalysisL1Extra::SetMuon(const edm::Handle<l1extra::L1MuonParticleCollection> muon, unsigned maxL1Extra)
{
  for(l1extra::L1MuonParticleCollection::const_iterator it=muon->begin(); it!=muon->end() && l1extra_.nMuons<maxL1Extra; it++){
      
      l1extra_.muonEt .push_back( it->et());
      l1extra_.muonEta.push_back(it->eta());
      l1extra_.muonPhi.push_back(it->phi());
      l1extra_.muonChg.push_back(it->charge());
      l1extra_.muonIso.push_back(it->isIsolated());
      l1extra_.muonMip.push_back(it->isMip());
      l1extra_.muonFwd.push_back(it->isForward());
      l1extra_.muonRPC.push_back(it->isRPC());
      l1extra_.muonBx .push_back(it->bx());
      l1extra_.muonQuality .push_back(it->gmtMuonCand().quality());
        
//      std::cout << "gmtmuon cand: pt " << it->gmtMuonCand().ptValue() 
//                  << "; ptExtra " << it->et() 
//                  << "; qual " << it->gmtMuonCand().quality() 
//                  << std::endl;
      l1extra_.nMuons++;
    }
}

void L1Analysis::L1AnalysisL1Extra::SetMet(const edm::Handle<l1extra::L1EtMissParticleCollection> mets)
{
  for(l1extra::L1EtMissParticleCollection::const_iterator it=mets->begin(); it!=mets->end(); it++) {
    l1extra_.et.    push_back( it->etTotal() ); 
    l1extra_.met.   push_back( it->et() );
    l1extra_.metPhi.push_back( it->phi() );
    l1extra_.metBx. push_back( it->bx() );
    l1extra_.nMet++;
  }
}

void L1Analysis::L1AnalysisL1Extra::SetMht(const edm::Handle<l1extra::L1EtMissParticleCollection> mhts)
{
  for(l1extra::L1EtMissParticleCollection::const_iterator it=mhts->begin(); it!=mhts->end(); it++) {
    l1extra_.ht.    push_back( it->etTotal() );
    l1extra_.mht.   push_back( it->et() );
    l1extra_.mhtPhi.push_back( it->phi() );
    l1extra_.mhtBx. push_back( it->bx() );
    l1extra_.nMht++;
  }
}

void L1Analysis::L1AnalysisL1Extra::SetHFring(const edm::Handle<l1extra::L1HFRingsCollection> hfRings)
{  
   l1extra_.hfEtSum.resize(4);
   l1extra_.hfBitCnt.resize(4);
   l1extra_.hfBx.resize(4);
  
    for(unsigned int i=0; i<4; ++i)  
    {
      if (hfRings->empty()) continue;

      l1extra_.hfEtSum[i] = (hfRings->begin()->hfEtSum((l1extra::L1HFRings::HFRingLabels) i));
      l1extra_.hfBitCnt[i] = (hfRings->begin()->hfBitCount((l1extra::L1HFRings::HFRingLabels) i));
      l1extra_.hfBx[i] = hfRings->begin()->bx();
      }
}