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

Inheritance diagram for L1CaloTowerTreeProducer:

Public Member Functions
	L1CaloTowerTreeProducer (const edm::ParameterSet &)

	~L1CaloTowerTreeProducer () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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
double	ecalLSB_

edm::EDGetTokenT < EcalTrigPrimDigiCollection >	ecalToken_

edm::Service< TFileService >	fs_

edm::EDGetTokenT < HcalTrigPrimDigiCollection >	hcalToken_

edm::EDGetTokenT < l1t::CaloClusterBxCollection >	l1ClusterToken_

edm::EDGetTokenT < l1t::CaloTowerBxCollection >	l1TowerToken_

unsigned	maxCaloTP_

unsigned	maxL1Cluster_

unsigned	maxL1Tower_

bool	storeCaloClusters_

TTree *	tree_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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 () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

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 ( const edm::ParameterSet & iConfig )

explicit

Definition at line 98 of file L1CaloTowerTreeProducer.cc.

References caloTPData_, ecalLSB_, ecalToken_, fs_, edm::ParameterSet::getUntrackedParameter(), hcalToken_, edm::InputTag::instance(), l1CaloClusterData_, l1CaloTowerData_, l1ClusterToken_, l1TowerToken_, edm::InputTag::label(), TFileService::make(), maxCaloTP_, maxL1Cluster_, maxL1Tower_, or, storeCaloClusters_, AlCaHLTBitMon_QueryRunRegistry::string, and tree_.

                                                                                {
   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"));
 
   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);
 
   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);
 
   // 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 ( )

override

Definition at line 129 of file L1CaloTowerTreeProducer.cc.

                                                   {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 139 of file L1CaloTowerTreeProducer.cc.

                                                                                          {
   // do calo towers
   caloTPData_->Reset();
 
   edm::ESHandle<CaloTPGTranscoder> decoder;
   iSetup.get<CaloTPGRecord>().get(decoder);
 
   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();
 }