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 final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const 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
bool	wantsStreamRuns () const
	~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(), EgHLTOffHistBins_cfi::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().