#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

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

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::vector< ModuleDescription const * > &modules, 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	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

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 64 of file L1CaloTowerTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

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

Definition at line 139 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

Definition at line 154 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 ( 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();
 
 }