L1CaloTowerTreeProducer Class Reference

#include <L1TriggerDPG/L1Ntuples/src/L1CaloTowerTreeProducer.cc>

Inheritance diagram for L1CaloTowerTreeProducer:

Public Member Functions
	L1CaloTowerTreeProducer (const edm::ParameterSet &)
	~L1CaloTowerTreeProducer () override=default
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Public Attributes
L1Analysis::L1AnalysisCaloTPDataFormat *	caloTPData_
L1Analysis::L1AnalysisL1CaloClusterDataFormat *	l1CaloClusterData_
L1Analysis::L1AnalysisL1CaloTowerDataFormat *	l1CaloTowerData_

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob (void) override
void	endJob () override

Private Attributes
const edm::ESGetToken< CaloTPGTranscoder, CaloTPGRecord >	decoderToken_
const double	ecalLSB_
const edm::EDGetTokenT< EcalTrigPrimDigiCollection >	ecalToken_
edm::Service< TFileService >	fs_
const edm::EDGetTokenT< HcalTrigPrimDigiCollection >	hcalToken_
edm::EDGetTokenT< l1t::CaloClusterBxCollection >	l1ClusterToken_
const edm::EDGetTokenT< l1t::CaloTowerBxCollection >	l1TowerToken_
const unsigned	maxCaloTP_
const unsigned	maxL1Cluster_
const unsigned	maxL1Tower_
bool	storeCaloClusters_
TTree *	tree_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: Produce L1 Extra tree

Implementation:

Definition at line 62 of file L1CaloTowerTreeProducer.cc.

Constructor & Destructor Documentation

L1CaloTowerTreeProducer()

L1CaloTowerTreeProducer::L1CaloTowerTreeProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 99 of file L1CaloTowerTreeProducer.cc.

References caloTPData_, StandaloneTrackMonitor_cfi::clusterTag, fs_, edm::ParameterSet::getUntrackedParameter(), TFileService::kSharedResource, l1CaloClusterData_, l1CaloTowerData_, l1ClusterToken_, TFileService::make(), or, storeCaloClusters_, AlCaHLTBitMon_QueryRunRegistry::string, and tree_.

    : ecalLSB_(iConfig.getUntrackedParameter<double>("ecalLSB", 0.5)),
      maxCaloTP_(iConfig.getUntrackedParameter<unsigned int>("maxCaloTP", 5760)),
      maxL1Tower_(iConfig.getUntrackedParameter<unsigned int>("maxL1Tower", 5760)),
      maxL1Cluster_(iConfig.getUntrackedParameter<unsigned int>("maxL1Cluster", 5760)),
      ecalToken_(consumes<EcalTrigPrimDigiCollection>(iConfig.getUntrackedParameter<edm::InputTag>("ecalToken"))),
      hcalToken_(consumes<HcalTrigPrimDigiCollection>(iConfig.getUntrackedParameter<edm::InputTag>("hcalToken"))),
      l1TowerToken_(consumes<l1t::CaloTowerBxCollection>(iConfig.getUntrackedParameter<edm::InputTag>("l1TowerToken"))),
      decoderToken_(esConsumes<CaloTPGTranscoder, CaloTPGRecord>()) {
  edm::InputTag clusterTag = iConfig.getUntrackedParameter<edm::InputTag>("l1ClusterToken");
  storeCaloClusters_ = true;
  if (clusterTag.label() == std::string("") or clusterTag.label() == std::string("none"))
    storeCaloClusters_ = false;

  if (clusterTag.instance() != std::string(""))
    l1ClusterToken_ = consumes<l1t::CaloClusterBxCollection>(clusterTag);

  usesResource(TFileService::kSharedResource);

  // set up output
  tree_ = fs_->make<TTree>("L1CaloTowerTree", "L1CaloTowerTree");
  tree_->Branch("CaloTP", "L1Analysis::L1AnalysisCaloTPDataFormat", &caloTPData_, 32000, 3);
  tree_->Branch("L1CaloTower", "L1Analysis::L1AnalysisL1CaloTowerDataFormat", &l1CaloTowerData_, 32000, 3);

  if (storeCaloClusters_)
    tree_->Branch("L1CaloCluster", "L1Analysis::L1AnalysisL1CaloClusterDataFormat", &l1CaloClusterData_, 32000, 3);

  caloTPData_ = new L1Analysis::L1AnalysisCaloTPDataFormat();
  l1CaloTowerData_ = new L1Analysis::L1AnalysisL1CaloTowerDataFormat();
  l1CaloClusterData_ = new L1Analysis::L1AnalysisL1CaloClusterDataFormat();
}

~L1CaloTowerTreeProducer()

L1CaloTowerTreeProducer::~L1CaloTowerTreeProducer ( )

overridedefault

Member Function Documentation

analyze()

void L1CaloTowerTreeProducer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 136 of file L1CaloTowerTreeProducer.cc.

References funct::abs(), BXVector< T >::begin(), caloTPData_, decoderToken_, ecalLSB_, ecalToken_, L1Analysis::L1AnalysisCaloTPDataFormat::ecalTPCaliphi, L1Analysis::L1AnalysisCaloTPDataFormat::ecalTPcompEt, L1Analysis::L1AnalysisCaloTPDataFormat::ecalTPet, L1Analysis::L1AnalysisCaloTPDataFormat::ecalTPfineGrain, L1Analysis::L1AnalysisCaloTPDataFormat::ecalTPieta, L1Analysis::L1AnalysisCaloTPDataFormat::ecalTPiphi, BXVector< T >::end(), l1tnanotables_cff::et, L1Analysis::L1AnalysisL1CaloClusterDataFormat::et, L1Analysis::L1AnalysisL1CaloTowerDataFormat::et, L1Analysis::L1AnalysisL1CaloClusterDataFormat::eta, L1Analysis::L1AnalysisL1CaloTowerDataFormat::eta, BXVector< T >::getFirstBX(), edm::EventSetup::getHandle(), BXVector< T >::getLastBX(), CaloTPGTranscoder::hcaletValue(), hcalToken_, L1Analysis::L1AnalysisCaloTPDataFormat::hcalTPCaliphi, L1Analysis::L1AnalysisCaloTPDataFormat::hcalTPcompEt, L1Analysis::L1AnalysisCaloTPDataFormat::hcalTPet, L1Analysis::L1AnalysisCaloTPDataFormat::hcalTPfineGrain, L1Analysis::L1AnalysisCaloTPDataFormat::hcalTPieta, L1Analysis::L1AnalysisCaloTPDataFormat::hcalTPiphi, L1Analysis::L1AnalysisL1CaloTowerDataFormat::iem, L1Analysis::L1AnalysisL1CaloClusterDataFormat::iet, L1Analysis::L1AnalysisL1CaloTowerDataFormat::iet, hcalRecHitTable_cff::ieta, L1Analysis::L1AnalysisL1CaloClusterDataFormat::ieta, L1Analysis::L1AnalysisL1CaloTowerDataFormat::ieta, iEvent, L1Analysis::L1AnalysisL1CaloTowerDataFormat::ihad, hcalRecHitTable_cff::iphi, L1Analysis::L1AnalysisL1CaloClusterDataFormat::iphi, L1Analysis::L1AnalysisL1CaloTowerDataFormat::iphi, L1Analysis::L1AnalysisL1CaloClusterDataFormat::iqual, L1Analysis::L1AnalysisL1CaloTowerDataFormat::iqual, L1Analysis::L1AnalysisL1CaloTowerDataFormat::iratio, edm::HandleBase::isValid(), l1CaloClusterData_, l1CaloTowerData_, l1ClusterToken_, l1TowerToken_, L1Analysis::L1AnalysisL1CaloClusterDataFormat::nCluster, L1Analysis::L1AnalysisCaloTPDataFormat::nECALTP, L1Analysis::L1AnalysisCaloTPDataFormat::nHCALTP, L1Analysis::L1AnalysisL1CaloTowerDataFormat::nTower, L1Analysis::L1AnalysisL1CaloClusterDataFormat::phi, L1Analysis::L1AnalysisL1CaloTowerDataFormat::phi, edm::Handle< T >::product(), L1Analysis::L1AnalysisL1CaloClusterDataFormat::Reset(), L1Analysis::L1AnalysisL1CaloTowerDataFormat::Reset(), L1Analysis::L1AnalysisCaloTPDataFormat::Reset(), storeCaloClusters_, and tree_.

                                                                                         {
  // do calo towers
  caloTPData_->Reset();

  edm::ESHandle<CaloTPGTranscoder> decoder;
  decoder = iSetup.getHandle(decoderToken_);

  edm::Handle<EcalTrigPrimDigiCollection> ecalTPs;
  edm::Handle<HcalTrigPrimDigiCollection> hcalTPs;

  iEvent.getByToken(ecalToken_, ecalTPs);
  iEvent.getByToken(hcalToken_, hcalTPs);

  if (ecalTPs.isValid()) {
    for (const auto& itr : *(ecalTPs.product())) {
      short ieta = (short)itr.id().ieta();
      //      unsigned short absIeta = (unsigned short) abs(ieta);
      //      short sign = ieta/absIeta;

      unsigned short cal_iphi = (unsigned short)itr.id().iphi();
      unsigned short iphi = (72 + 18 - cal_iphi) % 72;
      unsigned short compEt = itr.compressedEt();
      double et = ecalLSB_ * compEt;
      unsigned short fineGrain = (unsigned short)itr.fineGrain();

      if (compEt > 0) {
        caloTPData_->ecalTPieta.push_back(ieta);
        caloTPData_->ecalTPCaliphi.push_back(cal_iphi);
        caloTPData_->ecalTPiphi.push_back(iphi);
        caloTPData_->ecalTPet.push_back(et);
        caloTPData_->ecalTPcompEt.push_back(compEt);
        caloTPData_->ecalTPfineGrain.push_back(fineGrain);
        caloTPData_->nECALTP++;
      }
    }

  } else {
    edm::LogWarning("L1TNtuple") << "ECAL TPs not found, branch will not be filled";
  }

  if (hcalTPs.isValid()) {
    for (auto itr : (*hcalTPs.product())) {
      int ver = itr.id().version();
      short ieta = (short)itr.id().ieta();
      unsigned short absIeta = (unsigned short)abs(ieta);
      //      short sign = ieta/absIeta;

      unsigned short cal_iphi = (unsigned short)itr.id().iphi();
      unsigned short iphi = (72 + 18 - cal_iphi) % 72;

      unsigned short compEt = itr.SOI_compressedEt();
      double et = decoder->hcaletValue(itr.id(), itr.t0());

      unsigned short fineGrain = (unsigned short)itr.SOI_fineGrain();

      if (compEt > 0 && (absIeta < 29 || ver == 1)) {
        caloTPData_->hcalTPieta.push_back(ieta);
        caloTPData_->hcalTPCaliphi.push_back(cal_iphi);
        caloTPData_->hcalTPiphi.push_back(iphi);
        caloTPData_->hcalTPet.push_back(et);
        caloTPData_->hcalTPcompEt.push_back(compEt);
        caloTPData_->hcalTPfineGrain.push_back(fineGrain);
        caloTPData_->nHCALTP++;
      }
    }

  } else {
    edm::LogWarning("L1TNtuple") << "HCAL TPs not found, branch will not be filled";
  }

  // do L1 towers
  l1CaloTowerData_->Reset();

  edm::Handle<l1t::CaloTowerBxCollection> l1Towers;
  iEvent.getByToken(l1TowerToken_, l1Towers);

  if (l1Towers.isValid()) {
    for (int ibx = l1Towers->getFirstBX(); ibx <= l1Towers->getLastBX(); ++ibx) {
      for (auto itr = l1Towers->begin(ibx); itr != l1Towers->end(ibx); ++itr) {
        if (itr->hwPt() <= 0)
          continue;

        //  l1CaloTowerData_->bx.push_back( ibx );
        l1CaloTowerData_->et.push_back(itr->pt());
        l1CaloTowerData_->eta.push_back(itr->eta());
        l1CaloTowerData_->phi.push_back(itr->phi());
        l1CaloTowerData_->iet.push_back(itr->hwPt());
        l1CaloTowerData_->ieta.push_back(itr->hwEta());
        l1CaloTowerData_->iphi.push_back(itr->hwPhi());
        l1CaloTowerData_->iem.push_back(itr->hwEtEm());
        l1CaloTowerData_->ihad.push_back(itr->hwEtHad());
        l1CaloTowerData_->iratio.push_back(itr->hwEtRatio());
        l1CaloTowerData_->iqual.push_back(itr->hwQual());

        l1CaloTowerData_->nTower++;
      }
    }

  } else {
    edm::LogWarning("L1TNtuple") << "L1 Calo Towerss not found, branch will not be filled";
  }

  // do L1 clusters
  if (storeCaloClusters_) {
    l1CaloClusterData_->Reset();

    edm::Handle<l1t::CaloClusterBxCollection> l1Clusters;
    if (!l1ClusterToken_.isUninitialized())
      iEvent.getByToken(l1ClusterToken_, l1Clusters);

    if (l1Clusters.isValid()) {
      for (int ibx = l1Clusters->getFirstBX(); ibx <= l1Clusters->getLastBX(); ++ibx) {
        for (auto itr = l1Clusters->begin(ibx); itr != l1Clusters->end(ibx); ++itr) {
          if (itr->hwPt() <= 0)
            continue;

          //    l1CaloClusterData_->bx.push_back( ibx );
          l1CaloClusterData_->et.push_back(itr->pt());
          l1CaloClusterData_->eta.push_back(itr->eta());
          l1CaloClusterData_->phi.push_back(itr->phi());
          l1CaloClusterData_->iet.push_back(itr->hwPt());
          l1CaloClusterData_->ieta.push_back(itr->hwEta());
          l1CaloClusterData_->iphi.push_back(itr->hwPhi());
          l1CaloClusterData_->iqual.push_back(itr->hwQual());

          l1CaloClusterData_->nCluster++;
        }
      }

    } else {
      edm::LogWarning("L1TNtuple") << "L1 Calo Clusters not found, branch will not be filled";
    }
  }

  tree_->Fill();
}

beginJob()

void L1CaloTowerTreeProducer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 274 of file L1CaloTowerTreeProducer.cc.

{}

endJob()

void L1CaloTowerTreeProducer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 277 of file L1CaloTowerTreeProducer.cc.

{}

Member Data Documentation

caloTPData_

L1Analysis::L1AnalysisCaloTPDataFormat* L1CaloTowerTreeProducer::caloTPData_

Definition at line 73 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze(), and L1CaloTowerTreeProducer().

decoderToken_

const edm::ESGetToken<CaloTPGTranscoder, CaloTPGRecord> L1CaloTowerTreeProducer::decoderToken_

private

Definition at line 87 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze().

ecalLSB_

const double L1CaloTowerTreeProducer::ecalLSB_

private

Definition at line 78 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze().

ecalToken_

const edm::EDGetTokenT<EcalTrigPrimDigiCollection> L1CaloTowerTreeProducer::ecalToken_

private

Definition at line 84 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze().

fs_

edm::Service<TFileService> L1CaloTowerTreeProducer::fs_

private

Definition at line 91 of file L1CaloTowerTreeProducer.cc.

Referenced by L1CaloTowerTreeProducer().

hcalToken_

const edm::EDGetTokenT<HcalTrigPrimDigiCollection> L1CaloTowerTreeProducer::hcalToken_

private

Definition at line 85 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze().

l1CaloClusterData_

L1Analysis::L1AnalysisL1CaloClusterDataFormat* L1CaloTowerTreeProducer::l1CaloClusterData_

Definition at line 75 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze(), and L1CaloTowerTreeProducer().

l1CaloTowerData_

L1Analysis::L1AnalysisL1CaloTowerDataFormat* L1CaloTowerTreeProducer::l1CaloTowerData_

Definition at line 74 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze(), and L1CaloTowerTreeProducer().

l1ClusterToken_

edm::EDGetTokenT<l1t::CaloClusterBxCollection> L1CaloTowerTreeProducer::l1ClusterToken_

private

Definition at line 88 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze(), and L1CaloTowerTreeProducer().

l1TowerToken_

const edm::EDGetTokenT<l1t::CaloTowerBxCollection> L1CaloTowerTreeProducer::l1TowerToken_

private

Definition at line 86 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze().

maxCaloTP_

const unsigned L1CaloTowerTreeProducer::maxCaloTP_

private

Definition at line 79 of file L1CaloTowerTreeProducer.cc.

maxL1Cluster_

const unsigned L1CaloTowerTreeProducer::maxL1Cluster_

private

Definition at line 81 of file L1CaloTowerTreeProducer.cc.

maxL1Tower_

const unsigned L1CaloTowerTreeProducer::maxL1Tower_

private

Definition at line 80 of file L1CaloTowerTreeProducer.cc.

storeCaloClusters_

bool L1CaloTowerTreeProducer::storeCaloClusters_

private

Definition at line 96 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze(), and L1CaloTowerTreeProducer().

tree_

TTree* L1CaloTowerTreeProducer::tree_

private

Definition at line 94 of file L1CaloTowerTreeProducer.cc.

Referenced by analyze(), and L1CaloTowerTreeProducer().