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 18 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 6 of file L1AnalysisGCT.cc.

References verbose, and verbose_.

{ verbose_ = verbose; }

L1Analysis::L1AnalysisGCT::~L1AnalysisGCT

(

)

Definition at line 8 of file L1AnalysisGCT.cc.

{}

Member Function Documentation

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

inline

Definition at line 41 of file L1AnalysisGCT.h.

References gct_.

{ return &gct_; }

void L1Analysis::L1AnalysisGCT::Reset ( )

inline

Definition at line 39 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 160 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 84 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 134 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 9 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 45 of file L1AnalysisGCT.h.

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

bool L1Analysis::L1AnalysisGCT::verbose_

private

Definition at line 44 of file L1AnalysisGCT.h.

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