#include <EcalTPGAnalyzer.h>

Inheritance diagram for EcalTPGAnalyzer:

Classes
struct	EcalTPGVariables

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

	EcalTPGAnalyzer (const edm::ParameterSet &)

void	endRun (edm::Run const &, edm::EventSetup const &) override

	~EcalTPGAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	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)

Private Attributes
const bool	allowTP_

const edm::InputTag	digiCollectionEB_

const edm::InputTag	digiCollectionEE_

const edm::EDGetTokenT< EBDigiCollection >	ebDigiToken_

const edm::ESGetToken< CaloSubdetectorGeometry, EcalBarrelGeometryRecord >	ebGeometryToken_

const edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_

const edm::ESGetToken< CaloSubdetectorGeometry, EcalEndcapGeometryRecord >	eeGeometryToken_

const EcalTrigTowerConstituentsMap *	eTTmap_

const edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	eTTMapToken_

TFile *	file_

const std::string	gtRecordCollectionTag_

const edm::ESGetToken< L1GtTriggerMask, L1GtTriggerMaskAlgoTrigRcd >	l1GtMaskToken_

const edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	l1GtReadoutRecordToken_

const bool	print_

const CaloSubdetectorGeometry *	theBarrelGeometry_

const CaloSubdetectorGeometry *	theEndcapGeometry_

const edm::InputTag	tpCollection_

const edm::InputTag	tpEmulatorCollection_

const edm::EDGetTokenT< EcalTrigPrimDigiCollection >	tpEmulToken_

const edm::EDGetTokenT< EcalTrigPrimDigiCollection >	tpToken_

TTree *	tree_

EcalTPGVariables	treeVariables_

const bool	useEE_

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

Definition at line 42 of file EcalTPGAnalyzer.h.

Constructor & Destructor Documentation

◆ EcalTPGAnalyzer()

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

explicit

Definition at line 27 of file EcalTPGAnalyzer.cc.

     : tpCollection_(iConfig.getParameter<edm::InputTag>("TPCollection")),
       tpEmulatorCollection_(iConfig.getParameter<edm::InputTag>("TPEmulatorCollection")),
       digiCollectionEB_(iConfig.getParameter<edm::InputTag>("DigiCollectionEB")),
       digiCollectionEE_(iConfig.getParameter<edm::InputTag>("DigiCollectionEE")),
       gtRecordCollectionTag_(iConfig.getParameter<std::string>("GTRecordCollection")),
       l1GtReadoutRecordToken_(consumes<L1GlobalTriggerReadoutRecord>(edm::InputTag(gtRecordCollectionTag_))),
       tpToken_(consumes<EcalTrigPrimDigiCollection>(tpCollection_)),
       tpEmulToken_(consumes<EcalTrigPrimDigiCollection>(tpEmulatorCollection_)),
       ebDigiToken_(consumes<EBDigiCollection>(digiCollectionEB_)),
       eeDigiToken_(consumes<EEDigiCollection>(digiCollectionEE_)),
       eTTMapToken_(esConsumes<edm::Transition::BeginRun>()),
       ebGeometryToken_(esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", "EcalBarrel"))),
       eeGeometryToken_(esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", "EcalEndcap"))),
       l1GtMaskToken_(esConsumes()),
       allowTP_(iConfig.getParameter<bool>("ReadTriggerPrimitives")),
       useEE_(iConfig.getParameter<bool>("UseEndCap")),
       print_(iConfig.getParameter<bool>("Print")) {
   // file
   file_ = new TFile("ECALTPGtree.root", "RECREATE");
   file_->cd();
 
   // tree
   tree_ = new TTree("EcalTPGAnalysis", "EcalTPGAnalysis");
 
   tree_->Branch("runNb", &treeVariables_.runNb, "runNb/i");                                                //
   tree_->Branch("evtNb", &treeVariables_.evtNb, "evtNb/i");                                                //
   tree_->Branch("bxNb", &treeVariables_.bxNb, "bxNb/i");                                                   //
   tree_->Branch("orbitNb", &treeVariables_.orbitNb, "orbitNb/i");                                          //
   tree_->Branch("nbOfActiveTriggers", &treeVariables_.nbOfActiveTriggers, "nbOfActiveTriggers/i");         //
   tree_->Branch("activeTriggers", treeVariables_.activeTriggers, "activeTriggers[nbOfActiveTriggers]/I");  //
 
   tree_->Branch("nbOfTowers", &treeVariables_.nbOfTowers, "nbOfTowers/i");             //
   tree_->Branch("ieta", treeVariables_.ieta, "ieta[nbOfTowers]/I");                    //
   tree_->Branch("iphi", treeVariables_.iphi, "iphi[nbOfTowers]/I");                    //
   tree_->Branch("nbOfXtals", treeVariables_.nbOfXtals, "nbOfXtals[nbOfTowers]/I");     //
   tree_->Branch("rawTPData", treeVariables_.rawTPData, "rawTPData[nbOfTowers]/I");     //
   tree_->Branch("rawTPEmul1", treeVariables_.rawTPEmul1, "rawTPEmul1[nbOfTowers]/I");  //
   tree_->Branch("rawTPEmul2", treeVariables_.rawTPEmul2, "rawTPEmul2[nbOfTowers]/I");  //
   tree_->Branch("rawTPEmul3", treeVariables_.rawTPEmul3, "rawTPEmul3[nbOfTowers]/I");  //
   tree_->Branch("rawTPEmul4", treeVariables_.rawTPEmul4, "rawTPEmul4[nbOfTowers]/I");  //
   tree_->Branch("rawTPEmul5", treeVariables_.rawTPEmul5, "rawTPEmul5[nbOfTowers]/I");  //
   tree_->Branch("eRec", treeVariables_.eRec, "eRec[nbOfTowers]/F");                    //
 }

◆ ~EcalTPGAnalyzer()

EcalTPGAnalyzer::~EcalTPGAnalyzer ( )

override

Definition at line 72 of file EcalTPGAnalyzer.cc.

References file_, and tree_.

                                   {
   file_->cd();
   tree_->Write();
   file_->Close();
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 87 of file EcalTPGAnalyzer.cc.

                                                                                  {
   using namespace edm;
   using namespace std;
 
   if (print_)
     std::cout << "===========" << iEvent.id() << std::endl;
 
   map<EcalTrigTowerDetId, towerEner> mapTower;
   map<EcalTrigTowerDetId, towerEner>::iterator itTT;
 
   // get Evts info
 
   treeVariables_.runNb = iEvent.id().run();
   treeVariables_.evtNb = iEvent.id().event();
   treeVariables_.bxNb = iEvent.bunchCrossing();
   treeVariables_.orbitNb = iEvent.orbitNumber();
 
   // get L1 info
 
   edm::Handle<L1GlobalTriggerReadoutRecord> gtRecord;
   iEvent.getByToken(l1GtReadoutRecordToken_, gtRecord);
   DecisionWord dWord = gtRecord->decisionWord();  // this will get the decision word *before* masking disabled bits
 
   const auto& l1GtTmAlgo = iSetup.getData(l1GtMaskToken_);
   std::vector<unsigned int> triggerMaskAlgoTrig = l1GtTmAlgo.gtTriggerMask();
 
   // apply masks on algo
   int iDaq = 0;
   int iBit = -1;
   treeVariables_.nbOfActiveTriggers = 0;
   for (std::vector<bool>::iterator itBit = dWord.begin(); itBit != dWord.end(); ++itBit) {
     iBit++;
     int maskBit = triggerMaskAlgoTrig[iBit] & (1 << iDaq);
     if (maskBit)
       *itBit = false;
     if (*itBit) {
       treeVariables_.activeTriggers[treeVariables_.nbOfActiveTriggers] = iBit;
       treeVariables_.nbOfActiveTriggers++;
     }
   }
 
   // Get TP data
 
   edm::Handle<EcalTrigPrimDigiCollection> tp;
   iEvent.getByToken(tpToken_, tp);
   if (print_)
     std::cout << "TP collection size=" << tp.product()->size() << std::endl;
 
   for (unsigned int i = 0; i < tp.product()->size(); i++) {
     EcalTriggerPrimitiveDigi d = (*(tp.product()))[i];
     const EcalTrigTowerDetId TPtowid = d.id();
     towerEner tE;
     tE.iphi_ = TPtowid.iphi();
     tE.ieta_ = TPtowid.ieta();
     tE.tpgADC_ = d[0].raw();
     mapTower[TPtowid] = tE;
   }
 
   // Get Emulators TP
 
   edm::Handle<EcalTrigPrimDigiCollection> tpEmul;
   iEvent.getByToken(tpEmulToken_, tpEmul);
   if (print_)
     std::cout << "TPEmulator collection size=" << tpEmul.product()->size() << std::endl;
 
   for (unsigned int i = 0; i < tpEmul.product()->size(); i++) {
     EcalTriggerPrimitiveDigi d = (*(tpEmul.product()))[i];
     const EcalTrigTowerDetId TPtowid = d.id();
     itTT = mapTower.find(TPtowid);
     if (itTT != mapTower.end())
       for (int j = 0; j < 5; j++)
         (itTT->second).tpgEmul_[j] = d[j].raw();
   }
 
   // Get nb of crystals read out
 
   // Get EB xtal digi inputs
   edm::Handle<EBDigiCollection> digiEB;
   iEvent.getByToken(ebDigiToken_, digiEB);
 
   for (unsigned int i = 0; i < digiEB.product()->size(); i++) {
     const EBDataFrame& df = (*(digiEB.product()))[i];
     const EBDetId& id = df.id();
     const EcalTrigTowerDetId towid = id.tower();
     itTT = mapTower.find(towid);
     if (itTT != mapTower.end())
       (itTT->second).nbXtal_++;
   }
 
   if (useEE_) {
     // Get EE xtal digi inputs
     edm::Handle<EEDigiCollection> digiEE;
     iEvent.getByToken(eeDigiToken_, digiEE);
     for (unsigned int i = 0; i < digiEE.product()->size(); i++) {
       const EEDataFrame& df = (*(digiEE.product()))[i];
       const EEDetId& id = df.id();
       const EcalTrigTowerDetId towid = eTTmap_->towerOf(id);
       itTT = mapTower.find(towid);
       if (itTT != mapTower.end())
         (itTT->second).nbXtal_++;
     }
   }
 
   // Get rechits
 
   //... to be completed
 
   // fill tree
 
   treeVariables_.nbOfTowers = mapTower.size();
   int towerNb = 0;
   for (itTT = mapTower.begin(); itTT != mapTower.end(); ++itTT) {
     treeVariables_.ieta[towerNb] = (itTT->second).ieta_;
     treeVariables_.iphi[towerNb] = (itTT->second).iphi_;
     treeVariables_.nbOfXtals[towerNb] = (itTT->second).nbXtal_;
     treeVariables_.rawTPData[towerNb] = (itTT->second).tpgADC_;
     treeVariables_.rawTPEmul1[towerNb] = (itTT->second).tpgEmul_[0];
     treeVariables_.rawTPEmul2[towerNb] = (itTT->second).tpgEmul_[1];
     treeVariables_.rawTPEmul3[towerNb] = (itTT->second).tpgEmul_[2];
     treeVariables_.rawTPEmul4[towerNb] = (itTT->second).tpgEmul_[3];
     treeVariables_.rawTPEmul5[towerNb] = (itTT->second).tpgEmul_[4];
     treeVariables_.eRec[towerNb] = (itTT->second).eRec_;
     towerNb++;
   }
 
   tree_->Fill();
 }

◆ beginRun()

void EcalTPGAnalyzer::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	evtSetup
	)

override

Definition at line 78 of file EcalTPGAnalyzer.cc.

References ebGeometryToken_, eeGeometryToken_, eTTmap_, eTTMapToken_, edm::EventSetup::getData(), theBarrelGeometry_, and theEndcapGeometry_.

                                                                            {
   // geometry
   eTTmap_ = &evtSetup.getData(eTTMapToken_);
   theBarrelGeometry_ = &evtSetup.getData(ebGeometryToken_);
   theEndcapGeometry_ = &evtSetup.getData(eeGeometryToken_);
 }

◆ endRun()

void EcalTPGAnalyzer::endRun	(	edm::Run const &	,
		edm::EventSetup const &	evtSetup
	)

override

Definition at line 85 of file EcalTPGAnalyzer.cc.

85 {}

Member Data Documentation

◆ allowTP_

const bool EcalTPGAnalyzer::allowTP_

private

Definition at line 96 of file EcalTPGAnalyzer.h.

◆ digiCollectionEB_

const edm::InputTag EcalTPGAnalyzer::digiCollectionEB_

private

Definition at line 81 of file EcalTPGAnalyzer.h.

◆ digiCollectionEE_

const edm::InputTag EcalTPGAnalyzer::digiCollectionEE_

private

Definition at line 82 of file EcalTPGAnalyzer.h.

◆ ebDigiToken_

const edm::EDGetTokenT<EBDigiCollection> EcalTPGAnalyzer::ebDigiToken_

private

Definition at line 88 of file EcalTPGAnalyzer.h.

Referenced by analyze().

◆ ebGeometryToken_

const edm::ESGetToken<CaloSubdetectorGeometry, EcalBarrelGeometryRecord> EcalTPGAnalyzer::ebGeometryToken_

private

Definition at line 92 of file EcalTPGAnalyzer.h.

Referenced by beginRun().

◆ eeDigiToken_

const edm::EDGetTokenT<EEDigiCollection> EcalTPGAnalyzer::eeDigiToken_

private

Definition at line 89 of file EcalTPGAnalyzer.h.

Referenced by analyze().

◆ eeGeometryToken_

const edm::ESGetToken<CaloSubdetectorGeometry, EcalEndcapGeometryRecord> EcalTPGAnalyzer::eeGeometryToken_

private

Definition at line 93 of file EcalTPGAnalyzer.h.

Referenced by beginRun().

◆ eTTmap_

const EcalTrigTowerConstituentsMap* EcalTPGAnalyzer::eTTmap_

private

Definition at line 102 of file EcalTPGAnalyzer.h.

Referenced by analyze(), and beginRun().

◆ eTTMapToken_

const edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> EcalTPGAnalyzer::eTTMapToken_

private

Definition at line 91 of file EcalTPGAnalyzer.h.

Referenced by beginRun().

◆ file_

TFile* EcalTPGAnalyzer::file_

private

Definition at line 75 of file EcalTPGAnalyzer.h.

Referenced by EcalTPGAnalyzer(), and ~EcalTPGAnalyzer().

◆ gtRecordCollectionTag_

const std::string EcalTPGAnalyzer::gtRecordCollectionTag_

private

Definition at line 83 of file EcalTPGAnalyzer.h.

◆ l1GtMaskToken_

const edm::ESGetToken<L1GtTriggerMask, L1GtTriggerMaskAlgoTrigRcd> EcalTPGAnalyzer::l1GtMaskToken_

private

Definition at line 94 of file EcalTPGAnalyzer.h.

Referenced by analyze().

◆ l1GtReadoutRecordToken_

const edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> EcalTPGAnalyzer::l1GtReadoutRecordToken_

private

Definition at line 85 of file EcalTPGAnalyzer.h.

Referenced by analyze().

◆ print_

const bool EcalTPGAnalyzer::print_

private

Definition at line 98 of file EcalTPGAnalyzer.h.

Referenced by analyze().

◆ theBarrelGeometry_

const CaloSubdetectorGeometry* EcalTPGAnalyzer::theBarrelGeometry_

private

Definition at line 100 of file EcalTPGAnalyzer.h.

Referenced by beginRun().

◆ theEndcapGeometry_

const CaloSubdetectorGeometry* EcalTPGAnalyzer::theEndcapGeometry_

private

Definition at line 101 of file EcalTPGAnalyzer.h.

Referenced by beginRun().

◆ tpCollection_

const edm::InputTag EcalTPGAnalyzer::tpCollection_

private

Definition at line 79 of file EcalTPGAnalyzer.h.

◆ tpEmulatorCollection_

const edm::InputTag EcalTPGAnalyzer::tpEmulatorCollection_

private

Definition at line 80 of file EcalTPGAnalyzer.h.

◆ tpEmulToken_

const edm::EDGetTokenT<EcalTrigPrimDigiCollection> EcalTPGAnalyzer::tpEmulToken_

private

Definition at line 87 of file EcalTPGAnalyzer.h.

Referenced by analyze().

◆ tpToken_

const edm::EDGetTokenT<EcalTrigPrimDigiCollection> EcalTPGAnalyzer::tpToken_

private

Definition at line 86 of file EcalTPGAnalyzer.h.

Referenced by analyze().

◆ tree_

TTree* EcalTPGAnalyzer::tree_

private

Definition at line 76 of file EcalTPGAnalyzer.h.

Referenced by analyze(), EcalTPGAnalyzer(), and ~EcalTPGAnalyzer().

◆ treeVariables_

EcalTPGVariables EcalTPGAnalyzer::treeVariables_

private

Definition at line 77 of file EcalTPGAnalyzer.h.

Referenced by analyze(), and EcalTPGAnalyzer().

◆ useEE_

const bool EcalTPGAnalyzer::useEE_

private

Definition at line 97 of file EcalTPGAnalyzer.h.

Referenced by analyze().

Classes

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalTPGAnalyzer()

◆ ~EcalTPGAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endRun()

Member Data Documentation

◆ allowTP_

◆ digiCollectionEB_

◆ digiCollectionEE_

◆ ebDigiToken_

◆ ebGeometryToken_

◆ eeDigiToken_

◆ eeGeometryToken_

◆ eTTmap_

◆ eTTMapToken_

◆ file_

◆ gtRecordCollectionTag_

◆ l1GtMaskToken_

◆ l1GtReadoutRecordToken_

◆ print_

◆ theBarrelGeometry_

◆ theEndcapGeometry_

◆ tpCollection_

◆ tpEmulatorCollection_

◆ tpEmulToken_

◆ tpToken_

◆ tree_

◆ treeVariables_

◆ useEE_