L1Analysis::L1AnalysisGCT Class Reference

#include <L1AnalysisGCT.h>

Public Member Functions
L1AnalysisGCTDataFormat *	getData ()
	L1AnalysisGCT ()
	L1AnalysisGCT (bool verbose)
void	Reset ()
void	SetEm (const edm::Handle< L1GctEmCandCollection > l1IsoEm, const edm::Handle< L1GctEmCandCollection > l1NonIsoEm)
void	SetES (const edm::Handle< L1GctEtMissCollection > l1EtMiss, const edm::Handle< L1GctHtMissCollection > l1HtMiss, const edm::Handle< L1GctEtHadCollection > l1EtHad, const edm::Handle< L1GctEtTotalCollection > l1EtTotal)
void	SetHFminbias (const edm::Handle< L1GctHFRingEtSumsCollection > l1HFSums, const edm::Handle< L1GctHFBitCountsCollection > l1HFCounts)
void	SetJet (const edm::Handle< L1GctJetCandCollection > l1CenJets, const edm::Handle< L1GctJetCandCollection > l1ForJets, const edm::Handle< L1GctJetCandCollection > l1TauJets, const edm::Handle< L1GctJetCandCollection > l1IsoTauJets)
	~L1AnalysisGCT ()

Private Attributes
L1AnalysisGCTDataFormat	gct_
bool	verbose_

Detailed Description

Definition at line 19 of file L1AnalysisGCT.h.

Constructor & Destructor Documentation

L1Analysis::L1AnalysisGCT::L1AnalysisGCT

(

)

Definition at line 4 of file L1AnalysisGCT.cc.

                                      :verbose_(false)
{
}

L1Analysis::L1AnalysisGCT::L1AnalysisGCT

(

bool

verbose

)

Definition at line 8 of file L1AnalysisGCT.cc.

References verbose, and verbose_.

{
  verbose_ = verbose;
}

L1Analysis::L1AnalysisGCT::~L1AnalysisGCT

(

)

Definition at line 13 of file L1AnalysisGCT.cc.

{

}

Member Function Documentation

L1AnalysisGCTDataFormat* L1Analysis::L1AnalysisGCT::getData ( void  )

inline

Definition at line 42 of file L1AnalysisGCT.h.

References gct_.

{return &gct_;}

void L1Analysis::L1AnalysisGCT::Reset ( )

inline

Definition at line 40 of file L1AnalysisGCT.h.

References gct_, and L1Analysis::L1AnalysisGCTDataFormat::Reset().

{gct_.Reset();}

void L1Analysis::L1AnalysisGCT::SetEm	(	const edm::Handle< L1GctEmCandCollection >	l1IsoEm,
		const edm::Handle< L1GctEmCandCollection >	l1NonIsoEm
	)

Definition at line 208 of file L1AnalysisGCT.cc.

References gct_, L1Analysis::L1AnalysisGCTDataFormat::IsoEmBx, L1Analysis::L1AnalysisGCTDataFormat::IsoEmEta, L1Analysis::L1AnalysisGCTDataFormat::IsoEmPhi, L1Analysis::L1AnalysisGCTDataFormat::IsoEmRnk, L1Analysis::L1AnalysisGCTDataFormat::IsoEmSize, Vispa.Plugins.EdmBrowser.EdmDataAccessor::ne(), L1Analysis::L1AnalysisGCTDataFormat::NonIsoEmBx, L1Analysis::L1AnalysisGCTDataFormat::NonIsoEmEta, L1Analysis::L1AnalysisGCTDataFormat::NonIsoEmPhi, L1Analysis::L1AnalysisGCTDataFormat::NonIsoEmRnk, L1Analysis::L1AnalysisGCTDataFormat::NonIsoEmSize, and verbose_.

{   
        
      // Isolated EM
      if ( verbose_ ) {
        edm::LogInfo("L1Prompt") << "L1TGCT: number of iso em cands: " 
              << l1IsoEm->size() << std::endl;
      }
      
      gct_.IsoEmSize = l1IsoEm->size();
      for (L1GctEmCandCollection::const_iterator ie=l1IsoEm->begin(); ie!=l1IsoEm->end(); ie++) {
        
        gct_.IsoEmEta.push_back(ie->regionId().ieta());
        gct_.IsoEmPhi.push_back(ie->regionId().iphi());
        gct_.IsoEmRnk.push_back(ie->rank());
        gct_.IsoEmBx.push_back(ie->bx());
      } 

      // Non-isolated EM
      if ( verbose_ ) {
        edm::LogInfo("L1Prompt") << "L1TGCT: number of non-iso em cands: " 
              << l1NonIsoEm->size() << std::endl;
      }
      gct_.NonIsoEmSize = l1NonIsoEm->size();
      
      for (L1GctEmCandCollection::const_iterator ne=l1NonIsoEm->begin(); ne!=l1NonIsoEm->end(); ne++) {
        gct_.NonIsoEmEta.push_back(ne->regionId().ieta());
        gct_.NonIsoEmPhi.push_back(ne->regionId().iphi());
        gct_.NonIsoEmRnk.push_back(ne->rank());
    gct_.NonIsoEmBx.push_back(ne->bx());
          
      } 
   
}

void L1Analysis::L1AnalysisGCT::SetES	(	const edm::Handle< L1GctEtMissCollection >	l1EtMiss,
		const edm::Handle< L1GctHtMissCollection >	l1HtMiss,
		const edm::Handle< L1GctEtHadCollection >	l1EtHad,
		const edm::Handle< L1GctEtTotalCollection >	l1EtTotal
	)

Definition at line 108 of file L1AnalysisGCT.cc.

References L1Analysis::L1AnalysisGCTDataFormat::EtHad, L1Analysis::L1AnalysisGCTDataFormat::EtHadBX, L1Analysis::L1AnalysisGCTDataFormat::EtHadSize, L1Analysis::L1AnalysisGCTDataFormat::EtMiss, L1Analysis::L1AnalysisGCTDataFormat::EtMissBX, L1Analysis::L1AnalysisGCTDataFormat::EtMissPhi, L1Analysis::L1AnalysisGCTDataFormat::EtMissSize, L1Analysis::L1AnalysisGCTDataFormat::EtTot, L1Analysis::L1AnalysisGCTDataFormat::EtTotBX, L1Analysis::L1AnalysisGCTDataFormat::EtTotSize, gct_, L1Analysis::L1AnalysisGCTDataFormat::HtMiss, L1Analysis::L1AnalysisGCTDataFormat::HtMissBX, L1Analysis::L1AnalysisGCTDataFormat::HtMissPhi, L1Analysis::L1AnalysisGCTDataFormat::HtMissSize, and verbose_.

{ 
  
  // Energy sums
  for (L1GctEtMissCollection::const_iterator etm = l1EtMiss->begin();
       etm != l1EtMiss->end();
       ++etm) {

    gct_.EtMiss.push_back( etm->et() );
    gct_.EtMissPhi.push_back( etm->phi() );
    gct_.EtMissBX.push_back( etm->bx() );
    gct_.EtMissSize++;
    
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: Et Miss " 
                   << etm->et() << ", " <<  etm->phi()
                   << ", " <<  etm->bx() << std::endl;
    }
  }

   for (L1GctHtMissCollection::const_iterator htm = l1HtMiss->begin();
       htm != l1HtMiss->end();
       ++htm) {

    gct_.HtMiss.push_back( htm->et() );
    gct_.HtMissPhi.push_back( htm->phi() );
    gct_.HtMissBX.push_back( htm->bx() );
    gct_.HtMissSize++;
    
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: Ht Miss " 
                   << htm->et() << ", " <<  htm->phi()
                   << ", " <<  htm->bx() << std::endl;
    }
  }

   for (L1GctEtHadCollection::const_iterator ht = l1EtHad->begin();
       ht != l1EtHad->end();
       ++ht) {

    gct_.EtHad.push_back( ht->et() );
    gct_.EtHadBX.push_back( ht->bx() );
    gct_.EtHadSize++;
    
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: Ht Total " 
                   << ht->et()
                   << ", " <<  ht->bx() << std::endl;
    }
  }

  for (L1GctEtTotalCollection::const_iterator ett = l1EtTotal->begin();
       ett != l1EtTotal->end();
       ++ett) {

    gct_.EtTot.push_back( ett->et() );
    gct_.EtTotBX.push_back( ett->bx() );
    gct_.EtTotSize++;
    
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: Et Total " 
                   << ett->et()
                   << ", " <<  ett->bx() << std::endl;
    }
  }

}

void L1Analysis::L1AnalysisGCT::SetHFminbias	(	const edm::Handle< L1GctHFRingEtSumsCollection >	l1HFSums,
		const edm::Handle< L1GctHFBitCountsCollection >	l1HFCounts
	)

Definition at line 179 of file L1AnalysisGCT.cc.

References gct_, L1Analysis::L1AnalysisGCTDataFormat::HFBitCountsEta, L1Analysis::L1AnalysisGCTDataFormat::HFBitCountsSize, L1Analysis::L1AnalysisGCTDataFormat::HFRingEtSumEta, L1Analysis::L1AnalysisGCTDataFormat::HFRingEtSumSize, and verbose_.

{   
    
   //Fill HF Ring Histograms
    gct_.HFRingEtSumSize=l1HFSums->size();
    int ies=0;
    for (L1GctHFRingEtSumsCollection::const_iterator hfs=l1HFSums->begin(); hfs!=l1HFSums->end(); hfs++){ 
       gct_.HFRingEtSumEta.push_back(hfs->etSum(ies));
      if ( verbose_ ) {
    edm::LogInfo("L1Prompt") << "L1NtupleProducer: HF Sums " 
                     << l1HFSums->size() << ", " << hfs->etSum(ies) << std::endl;
      }
      ies++;
    }
    
    int ibc=0;
    gct_.HFBitCountsSize=l1HFCounts->size();
    for (L1GctHFBitCountsCollection::const_iterator hfc=l1HFCounts->begin(); hfc!=l1HFCounts->end(); hfc++){ 
      gct_.HFBitCountsEta.push_back(hfc->bitCount(ibc));
      if ( verbose_ ) {
    edm::LogInfo("L1Prompt") << "L1NtupleProducer: HF Counts " 
                     << l1HFCounts->size() << ", " << hfc->bitCount(ibc) << std::endl;
      }
      ibc++;
    }
   
}

void L1Analysis::L1AnalysisGCT::SetJet	(	const edm::Handle< L1GctJetCandCollection >	l1CenJets,
		const edm::Handle< L1GctJetCandCollection >	l1ForJets,
		const edm::Handle< L1GctJetCandCollection >	l1TauJets,
		const edm::Handle< L1GctJetCandCollection >	l1IsoTauJets
	)

Definition at line 17 of file L1AnalysisGCT.cc.

References L1Analysis::L1AnalysisGCTDataFormat::CJetBx, L1Analysis::L1AnalysisGCTDataFormat::CJetEta, L1Analysis::L1AnalysisGCTDataFormat::CJetPhi, L1Analysis::L1AnalysisGCTDataFormat::CJetRnk, L1Analysis::L1AnalysisGCTDataFormat::CJetSize, L1Analysis::L1AnalysisGCTDataFormat::FJetBx, L1Analysis::L1AnalysisGCTDataFormat::FJetEta, L1Analysis::L1AnalysisGCTDataFormat::FJetPhi, L1Analysis::L1AnalysisGCTDataFormat::FJetRnk, L1Analysis::L1AnalysisGCTDataFormat::FJetSize, gct_, L1Analysis::L1AnalysisGCTDataFormat::IsoTJetBx, L1Analysis::L1AnalysisGCTDataFormat::IsoTJetEta, L1Analysis::L1AnalysisGCTDataFormat::IsoTJetPhi, L1Analysis::L1AnalysisGCTDataFormat::IsoTJetRnk, L1Analysis::L1AnalysisGCTDataFormat::IsoTJetSize, edm::HandleBase::isValid(), L1Analysis::L1AnalysisGCTDataFormat::TJetBx, L1Analysis::L1AnalysisGCTDataFormat::TJetEta, L1Analysis::L1AnalysisGCTDataFormat::TJetPhi, L1Analysis::L1AnalysisGCTDataFormat::TJetRnk, L1Analysis::L1AnalysisGCTDataFormat::TJetSize, and verbose_.

{   
  // std::cout << "L1IsoTauJets valid: " << l1IsoTauJets.isValid() <<std::endl;
  // edm::LogAbsolute("L1") << "CCLA: " << l1IsoTauJets.isValid() <<std::endl;
    // Central jets
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: number of central jets = " 
        << l1CenJets->size() << std::endl;
    }
    gct_.CJetSize= l1CenJets->size();//1
    for (L1GctJetCandCollection::const_iterator cj = l1CenJets->begin();
     cj != l1CenJets->end(); cj++) {
      gct_.CJetEta.push_back(cj->regionId().ieta());//2
      gct_.CJetPhi.push_back(cj->regionId().iphi());//3
      gct_.CJetRnk.push_back(cj->rank());//4
      gct_.CJetBx .push_back(cj->bx());//4
      if ( verbose_ ) {
    edm::LogInfo("L1Prompt") << "L1NtupleProducer: Central jet " 
          << cj->regionId().iphi() << ", " << cj->regionId().ieta()
          << ", " << cj->rank() << std::endl;
      }
    }

    // Forward jets
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: number of forward jets = " 
        << l1ForJets->size() << std::endl;
    }
    gct_.FJetSize= l1ForJets->size();//5
    for (L1GctJetCandCollection::const_iterator fj = l1ForJets->begin();
     fj != l1ForJets->end(); fj++) {
      gct_.FJetEta.push_back(fj->regionId().ieta());//6
      gct_.FJetPhi.push_back(fj->regionId().iphi());//7
      gct_.FJetRnk.push_back(fj->rank());//8
      gct_.FJetBx .push_back(fj->bx());//8
      if ( verbose_ ) {
    edm::LogInfo("L1Prompt") << "L1NtupleProducer: Forward jet " 
          << fj->regionId().iphi() << ", " << fj->regionId().ieta()
          << ", " << fj->rank() << std::endl;
      }
    }

    // Tau jets
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: number of tau jets = " 
        << l1TauJets->size() << std::endl;
    }
    gct_.TJetSize= l1TauJets->size();//9
      for (L1GctJetCandCollection::const_iterator tj = l1TauJets->begin();
     tj != l1TauJets->end(); tj++) {
      //if ( tj->rank() == 0 ) continue;
      gct_.TJetEta.push_back(tj->regionId().ieta());//10
      gct_.TJetPhi.push_back(tj->regionId().iphi());//11
      gct_.TJetRnk.push_back(tj->rank());//12
      gct_.TJetBx .push_back(tj->bx());//
      if ( verbose_ ) {
    edm::LogInfo("L1Prompt") << "L1NtupleProducer: Tau jet " 
                   << tj->regionId().iphi() << ", " << tj->regionId().ieta()
                   << ", " << tj->rank() << std::endl;
      }
    }

    // Isolated Tau Jets
    if (l1IsoTauJets.isValid()){

      if ( verbose_ ) {
    edm::LogInfo("L1Prompt") << "L1NtupleProducer: number of isoTau jets = " 
                 << l1IsoTauJets->size() << std::endl;
      }
      gct_.IsoTJetSize= l1IsoTauJets->size();//9
      for (L1GctJetCandCollection::const_iterator tj = l1IsoTauJets->begin();
       tj != l1IsoTauJets->end(); tj++) {
    //if ( tj->rank() == 0 ) continue;
    gct_.IsoTJetEta.push_back(tj->regionId().ieta());//10
    gct_.IsoTJetPhi.push_back(tj->regionId().iphi());//11
    gct_.IsoTJetRnk.push_back(tj->rank());//12
    gct_.IsoTJetBx .push_back(tj->bx());//
    if ( verbose_ ) {
      edm::LogInfo("L1Prompt") << "L1NtupleProducer: IsoTau jet " 
                   << tj->regionId().iphi() << ", " << tj->regionId().ieta()
                   << ", " << tj->rank() << std::endl;
    }
      }
      
    }
  
}

Member Data Documentation

L1AnalysisGCTDataFormat L1Analysis::L1AnalysisGCT::gct_

private

Definition at line 46 of file L1AnalysisGCT.h.

Referenced by getData(), Reset(), SetEm(), SetES(), SetHFminbias(), and SetJet().

bool L1Analysis::L1AnalysisGCT::verbose_

private

Definition at line 45 of file L1AnalysisGCT.h.

Referenced by L1AnalysisGCT(), SetEm(), SetES(), SetHFminbias(), and SetJet().