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

Inheritance diagram for L1CaloTowerTreeProducer:

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

	~L1CaloTowerTreeProducer ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Public Attributes
L1Analysis::L1AnalysisCaloTPDataFormat *	caloTPData_

L1Analysis::L1AnalysisL1CaloClusterDataFormat *	l1CaloClusterData_

L1Analysis::L1AnalysisL1CaloTowerDataFormat *	l1CaloTowerData_

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob (void)

virtual void	endJob ()

Private Attributes
unsigned long long	ecalScaleCacheID_
	scales More...

edm::EDGetTokenT < EcalTrigPrimDigiCollection >	ecalToken_

edm::Service< TFileService >	fs_

unsigned long long	hcalScaleCacheID_

edm::EDGetTokenT < HcalTrigPrimDigiCollection >	hcalToken_

edm::EDGetTokenT < l1t::CaloClusterBxCollection >	l1ClusterToken_

edm::EDGetTokenT < l1t::CaloTowerBxCollection >	l1TowerToken_

unsigned	maxCaloTP_

unsigned	maxL1Cluster_

unsigned	maxL1Tower_

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 &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

Detailed Description

Description: Produce L1 Extra tree

Implementation:

Definition at line 66 of file L1CaloTowerTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 109 of file L1CaloTowerTreeProducer.cc.

References caloTPData_, ecalToken_, fs_, edm::ParameterSet::getUntrackedParameter(), hcalToken_, edm::InputTag::instance(), l1CaloClusterData_, l1CaloTowerData_, l1ClusterToken_, l1TowerToken_, TFileService::make(), maxCaloTP_, maxL1Cluster_, maxL1Tower_, 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");
   if (clusterTag.instance() != std::string(""))
     l1ClusterToken_ = consumes<l1t::CaloClusterBxCollection>(clusterTag);
  
   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);
   tree_->Branch("L1CaloCluster", "L1Analysis::L1AnalysisL1CaloClusterDataFormat", &l1CaloClusterData_, 32000, 3);
 
   caloTPData_ = new L1Analysis::L1AnalysisCaloTPDataFormat();
   l1CaloTowerData_ = new L1Analysis::L1AnalysisL1CaloTowerDataFormat();
   l1CaloClusterData_ = new L1Analysis::L1AnalysisL1CaloClusterDataFormat();
 
 }

L1CaloTowerTreeProducer::~L1CaloTowerTreeProducer ( )

Definition at line 137 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
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 152 of file L1CaloTowerTreeProducer.cc.

 {
 
   // do calo towers
   caloTPData_->Reset();
 
   edm::ESHandle<L1CaloEcalScale> ecalScale;
   if( iSetup.get< L1CaloEcalScaleRcd > ().cacheIdentifier() != ecalScaleCacheID_ ) {
     iSetup.get<L1CaloEcalScaleRcd>().get(ecalScale);
     ecalScaleCacheID_ = iSetup.get< L1CaloEcalScaleRcd > ().cacheIdentifier();
   }
 
   edm::ESHandle<L1CaloHcalScale> hcalScale;
   iSetup.get<L1CaloHcalScaleRcd>().get(hcalScale);
   hcalScaleCacheID_ = iSetup.get< L1CaloHcalScaleRcd > ().cacheIdentifier();
 
   edm::Handle<EcalTrigPrimDigiCollection> ecalTPs;
   edm::Handle<HcalTrigPrimDigiCollection> hcalTPs;
 
   iEvent.getByToken(ecalToken_, ecalTPs);
   iEvent.getByToken(hcalToken_, hcalTPs);
 
   if (ecalTPs.isValid()){ 
 
     for ( 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 = 0.;
       if (ecalScale.isValid()) et = ecalScale->et( compEt, absIeta, sign );
       
       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") << "HCAL TPs not found, branch will not be filled";
   }
   
   
   if (hcalTPs.isValid()) {
       
     for ( auto itr : (*hcalTPs.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;
       if (absIeta >= 29) {  // special treatment for HF
     iphi = iphi/4;
       }
       
       unsigned short compEt = itr.SOI_compressedEt();
       double et = 0.;
       if (hcalScale.isValid()) et = hcalScale->et( compEt, absIeta, sign );
       
       unsigned short fineGrain = (unsigned short) itr.SOI_fineGrain();
       
       if (compEt > 0) {
     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
   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();
 
 }