L1AnalysisGCT.cc

Go to the documentation of this file.

#include "L1Trigger/L1TNtuples/interface/L1AnalysisGCT.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

L1Analysis::L1AnalysisGCT::L1AnalysisGCT() : verbose_(false) {}

L1Analysis::L1AnalysisGCT::L1AnalysisGCT(bool verbose) { verbose_ = verbose; }

L1Analysis::L1AnalysisGCT::~L1AnalysisGCT() {}
void L1Analysis::L1AnalysisGCT::SetJet(const edm::Handle<L1GctJetCandCollection> l1CenJets,
                                       const edm::Handle<L1GctJetCandCollection> l1ForJets,
                                       const edm::Handle<L1GctJetCandCollection> l1TauJets,
                                       const edm::Handle<L1GctJetCandCollection> l1IsoTauJets) {
  // 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;
      }
    }
  }
}

void L1Analysis::L1AnalysisGCT::SetES(const edm::Handle<L1GctEtMissCollection> l1EtMiss,
                                      const edm::Handle<L1GctHtMissCollection> l1HtMiss,
                                      const edm::Handle<L1GctEtHadCollection> l1EtHad,
                                      const edm::Handle<L1GctEtTotalCollection> l1EtTotal) {
  // 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) {
  //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::SetEm(const edm::Handle<L1GctEmCandCollection> l1IsoEm,
                                      const edm::Handle<L1GctEmCandCollection> l1NonIsoEm) {
  // 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());
  }
}