EcalTPGAnalyzer Class Reference

#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 &)
	~EcalTPGAnalyzer () 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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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)

Private Attributes
bool	allowTP_
edm::InputTag	digiCollectionEB_
edm::InputTag	digiCollectionEE_
edm::ESHandle< EcalTrigTowerConstituentsMap >	eTTmap_
TFile *	file_
std::string	gtRecordCollectionTag_
bool	print_
const CaloSubdetectorGeometry *	theBarrelGeometry_
const CaloSubdetectorGeometry *	theEndcapGeometry_
edm::InputTag	tpCollection_
edm::InputTag	tpEmulatorCollection_
TTree *	tree_
EcalTPGVariables	treeVariables_
bool	useEE_

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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 35 of file EcalTPGAnalyzer.h.

Constructor & Destructor Documentation

EcalTPGAnalyzer()

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

explicit

Definition at line 50 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");
 
  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");                    //
}

References edm::ParameterSet::getParameter(), edm::hash_detail::print_(), and AlCaHLTBitMon_QueryRunRegistry::string.

~EcalTPGAnalyzer()

EcalTPGAnalyzer::~EcalTPGAnalyzer ( )

override

Definition at line 88 of file EcalTPGAnalyzer.cc.

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

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 108 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.getByLabel(edm::InputTag(gtRecordCollectionTag_), gtRecord);
  DecisionWord dWord = gtRecord->decisionWord();  // this will get the decision word *before* masking disabled bits
 
  edm::ESHandle<L1GtTriggerMask> l1GtTmAlgo;
  iSetup.get<L1GtTriggerMaskAlgoTrigRcd>().get(l1GtTmAlgo);
  std::vector<unsigned int> triggerMaskAlgoTrig = l1GtTmAlgo.product()->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.getByLabel(tpCollection_, 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.getByLabel(tpEmulatorCollection_, 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.getByLabel(digiCollectionEB_, 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.getByLabel(digiCollectionEE_, 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();
}

References gather_cfg::cout, ztail::d, L1GlobalTriggerReadoutRecord::decisionWord(), hgcalPerformanceValidation::df, edm::EventSetup::get(), edm::get(), L1GtTriggerMask::gtTriggerMask(), mps_fire::i, EcalTrigTowerDetId::ieta(), towerEner::ieta_, iEvent, EcalTrigTowerDetId::iphi(), towerEner::iphi_, dqmiolumiharvest::j, edm::hash_detail::print_(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), edm::SortedCollection< T, SORT >::size(), edm::DataFrameContainer::size(), cmsswSequenceInfo::tp, and towerEner::tpgADC_.

beginRun()

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

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 94 of file EcalTPGAnalyzer.cc.

                                                                           {
  // geometry
  ESHandle<CaloGeometry> theGeometry;
  ESHandle<CaloSubdetectorGeometry> theEndcapGeometry_handle, theBarrelGeometry_handle;
 
  evtSetup.get<CaloGeometryRecord>().get(theGeometry);
  evtSetup.get<EcalEndcapGeometryRecord>().get("EcalEndcap", theEndcapGeometry_handle);
  evtSetup.get<EcalBarrelGeometryRecord>().get("EcalBarrel", theBarrelGeometry_handle);
 
  evtSetup.get<IdealGeometryRecord>().get(eTTmap_);
  theEndcapGeometry_ = &(*theEndcapGeometry_handle);
  theBarrelGeometry_ = &(*theBarrelGeometry_handle);
}

References edm::EventSetup::get(), and edm::get().

Member Data Documentation

allowTP_

bool EcalTPGAnalyzer::allowTP_

private

Definition at line 77 of file EcalTPGAnalyzer.h.

digiCollectionEB_

edm::InputTag EcalTPGAnalyzer::digiCollectionEB_

private

Definition at line 73 of file EcalTPGAnalyzer.h.

digiCollectionEE_

edm::InputTag EcalTPGAnalyzer::digiCollectionEE_

private

Definition at line 74 of file EcalTPGAnalyzer.h.

eTTmap_

edm::ESHandle<EcalTrigTowerConstituentsMap> EcalTPGAnalyzer::eTTmap_

private

Definition at line 83 of file EcalTPGAnalyzer.h.

file_

TFile* EcalTPGAnalyzer::file_

private

Definition at line 67 of file EcalTPGAnalyzer.h.

gtRecordCollectionTag_

std::string EcalTPGAnalyzer::gtRecordCollectionTag_

private

Definition at line 75 of file EcalTPGAnalyzer.h.

print_

bool EcalTPGAnalyzer::print_

private

Definition at line 79 of file EcalTPGAnalyzer.h.

theBarrelGeometry_

const CaloSubdetectorGeometry* EcalTPGAnalyzer::theBarrelGeometry_

private

Definition at line 82 of file EcalTPGAnalyzer.h.

theEndcapGeometry_

const CaloSubdetectorGeometry* EcalTPGAnalyzer::theEndcapGeometry_

private

Definition at line 81 of file EcalTPGAnalyzer.h.

tpCollection_

edm::InputTag EcalTPGAnalyzer::tpCollection_

private

Definition at line 71 of file EcalTPGAnalyzer.h.

tpEmulatorCollection_

edm::InputTag EcalTPGAnalyzer::tpEmulatorCollection_

private

Definition at line 72 of file EcalTPGAnalyzer.h.

tree_

TTree* EcalTPGAnalyzer::tree_

private

Definition at line 68 of file EcalTPGAnalyzer.h.

treeVariables_

EcalTPGVariables EcalTPGAnalyzer::treeVariables_

private

Definition at line 69 of file EcalTPGAnalyzer.h.

useEE_

bool EcalTPGAnalyzer::useEE_

private

Definition at line 78 of file EcalTPGAnalyzer.h.