#include <DPGAnalysis/SiStripTools/src/MultiplicityTimeCorrelations.cc>

Inheritance diagram for MultiplicityTimeCorrelations:

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

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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

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

void	endJob () override

Private Attributes
edm::EDGetTokenT< APVCyclePhaseCollection >	_apvphasecollToken

std::map< unsigned int, int >	_binmax

std::vector< int >	_dbxbins

std::map< int, std::map< unsigned int, TH1F * > >	_dbxhistos

DigiBXCorrHistogramMaker< EventWithHistory >	_digibxcorrhmevent

EventWithHistoryFilter	_evfilter

edm::InputTag	_hecollection

edm::EDGetTokenT< EventWithHistory >	_hecollectionToken

int	_highorbit

int	_loworbit

int	_mindbx

int	_mintrpltdbx

edm::EDGetTokenT< std::map< unsigned int, int > >	_multiplicityMapToken

std::map< unsigned int, std::string >	_subdets

SiStripTKNumbers	_trnumb

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

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 60 of file MultiplicityTimeCorrelations.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 114 of file MultiplicityTimeCorrelations.cc.

References _binmax, _dbxbins, _dbxhistos, _digibxcorrhmevent, _evfilter, _hecollection, _mindbx, _mintrpltdbx, _subdets, _trnumb, edm::ParameterSet::addUntrackedParameter(), newFWLiteAna::bin, DigiBXCorrHistogramMaker< T >::book(), edm::EDConsumerBase::consumesCollector(), edm::ParameterSet::getUntrackedParameter(), createfilelist::int, SummaryClient_cfi::labels, LogDebug, TFileDirectory::make(), TFileService::mkdir(), SiStripTKNumbers::nstrips(), EventWithHistoryFilter::set(), and AlCaHLTBitMon_QueryRunRegistry::string.

     : _digibxcorrhmevent(iConfig, consumesCollector()),
       _evfilter(),
       _hecollection(iConfig.getParameter<edm::InputTag>("historyProduct")),
       _hecollectionToken(consumes<EventWithHistory>(_hecollection)),
       _apvphasecollToken(consumes<APVCyclePhaseCollection>(iConfig.getParameter<edm::InputTag>("apvPhaseCollection"))),
       _multiplicityMapToken(
           mayConsume<std::map<unsigned int, int> >(iConfig.getParameter<edm::InputTag>("multiplicityMap"))),
       _subdets(),
       _binmax(),
       _loworbit(iConfig.getUntrackedParameter<int>("lowedgeOrbit")),
       _highorbit(iConfig.getUntrackedParameter<int>("highedgeOrbit")),
       _mindbx(iConfig.getUntrackedParameter<int>("minDBX")),
       _mintrpltdbx(iConfig.getUntrackedParameter<int>("minTripletDBX")),
       _trnumb(),
       _dbxbins(iConfig.getUntrackedParameter<std::vector<int> >("dbxBins")) {
   //now do what ever initialization is needed
 
   // configure the filter
 
   edm::ParameterSet filterConfig;
   filterConfig.addUntrackedParameter<edm::InputTag>("historyProduct", _hecollection);
   if (_mindbx > 0) {
     std::vector<int> dbxrange;
     dbxrange.push_back(_mindbx + 1);
     dbxrange.push_back(-1);
     filterConfig.addUntrackedParameter<std::vector<int> >("dbxRange", dbxrange);
   }
   if (_mintrpltdbx > 0) {
     std::vector<int> dbxtrpltrange;
     dbxtrpltrange.push_back(_mintrpltdbx + 1);
     dbxtrpltrange.push_back(-1);
     filterConfig.addUntrackedParameter<std::vector<int> >("dbxTripletRange", dbxtrpltrange);
   }
 
   _evfilter.set(filterConfig, consumesCollector());
 
   //
 
   edm::Service<TFileService> tfserv;
 
   // create map of labels
 
   std::vector<edm::ParameterSet> wantedsubds(
       iConfig.getUntrackedParameter<std::vector<edm::ParameterSet> >("wantedSubDets"));
 
   for (std::vector<edm::ParameterSet>::iterator ps = wantedsubds.begin(); ps != wantedsubds.end(); ++ps) {
     _subdets[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<std::string>("detLabel");
     _binmax[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<int>("binMax");
   }
   std::map<int, std::string> labels;
 
   for (std::map<unsigned int, std::string>::const_iterator subd = _subdets.begin(); subd != _subdets.end(); ++subd) {
     labels[int(subd->first)] = subd->second;
   }
 
   //
 
   _digibxcorrhmevent.book("EventProcs", labels);
 
   TFileDirectory subdbxbin = tfserv->mkdir("DBXDebugging");
 
   for (std::vector<int>::const_iterator bin = _dbxbins.begin(); bin != _dbxbins.end(); bin++) {
     char hname[200];
     char htitle[200];
 
     edm::LogInfo("DBXHistosBinMaxValue") << "Setting bin max values";
 
     for (std::map<unsigned int, std::string>::const_iterator subd = _subdets.begin(); subd != _subdets.end(); ++subd) {
       if (_binmax.find(subd->first) == _binmax.end()) {
         edm::LogVerbatim("DBXHistosNotConfiguredBinMax")
             << "Bin max for " << subd->second << " not configured: " << _trnumb.nstrips(int(subd->first)) << " used";
         _binmax[subd->first] = _trnumb.nstrips(int(subd->first));
       }
 
       edm::LogVerbatim("DBXHistosBinMaxValue") << "Bin max for " << subd->second << " is " << _binmax[subd->first];
 
       sprintf(hname, "sumn%sdigi_%d", subd->second.c_str(), *bin);
       sprintf(htitle, "%s digi multiplicity at DBX = %d", subd->second.c_str(), *bin);
       LogDebug("DBXDebug") << "creating histogram " << hname << " " << htitle;
       _dbxhistos[*bin][subd->first] =
           subdbxbin.make<TH1F>(hname, htitle, 1000, 0., _binmax[subd->first] / (20 * 1000) * 1000);
       _dbxhistos[*bin][subd->first]->GetXaxis()->SetTitle("Number of Digis");
     }
     /*
     sprintf(hname,"sumntkdigi_%d",*bin);
     sprintf(htitle,"TK digi multiplicity at DBX = %d",*bin);
     LogDebug("DBXDebug") << "creating histogram " << hname << " " << htitle;
     _dbxtkhistos[*bin]= subdbxbin.make<TH1F>(hname,htitle,1000,0.,_trnumb.nstrips(0)/(20*1000)*1000);
     _dbxtkhistos[*bin]->GetXaxis()->SetTitle("Number of Digis");
 
     sprintf(hname,"sumntibdigi_%d",*bin);
     sprintf(htitle,"TIB digi multiplicity at DBX = %d",*bin);
     LogDebug("DBXDebug") << "creating histogram " << hname << " " << htitle;
     _dbxtibhistos[*bin]= subdbxbin.make<TH1F>(hname,htitle,1000,0.,_trnumb.nstrips(SiStripDetId::TIB)/(20*1000)*1000);
     _dbxtibhistos[*bin]->GetXaxis()->SetTitle("Number of Digis");
 
     sprintf(hname,"sumntiddigi_%d",*bin);
     sprintf(htitle,"TID digi multiplicity at DBX = %d",*bin);
     LogDebug("DBXDebug") << "creating histogram " << hname << " " << htitle;
     _dbxtidhistos[*bin]= subdbxbin.make<TH1F>(hname,htitle,1000,0.,_trnumb.nstrips(SiStripDetId::TID)/(20*1000)*1000);
     _dbxtidhistos[*bin]->GetXaxis()->SetTitle("Number of Digis");
 
     sprintf(hname,"sumntobdigi_%d",*bin);
     sprintf(htitle,"TOB digi multiplicity at DBX = %d",*bin);
     LogDebug("DBXDebug") << "creating histogram " << hname << " " << htitle;
     _dbxtobhistos[*bin]= subdbxbin.make<TH1F>(hname,htitle,1000,0.,_trnumb.nstrips(SiStripDetId::TOB)/(20*1000)*1000);
     _dbxtobhistos[*bin]->GetXaxis()->SetTitle("Number of Digis");
 
     sprintf(hname,"sumntecdigi_%d",*bin);
     sprintf(htitle,"TEC digi multiplicity at DBX = %d",*bin);
     LogDebug("DBXDebug") << "creating histogram " << hname << " " << htitle;
     _dbxtechistos[*bin]= subdbxbin.make<TH1F>(hname,htitle,1000,0.,_trnumb.nstrips(SiStripDetId::TEC)/(20*1000)*1000);
     _dbxtechistos[*bin]->GetXaxis()->SetTitle("Number of Digis");
     */
   }
 }

MultiplicityTimeCorrelations::~MultiplicityTimeCorrelations ( )

override

Definition at line 232 of file MultiplicityTimeCorrelations.cc.

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

Member Function Documentation

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

overrideprivate

Definition at line 242 of file MultiplicityTimeCorrelations.cc.

References _apvphasecollToken, _dbxhistos, _digibxcorrhmevent, _evfilter, _hecollectionToken, _highorbit, _loworbit, _multiplicityMapToken, _subdets, EventWithHistory::deltaBX(), EventWithHistory::depth(), DigiBXCorrHistogramMaker< T >::fill(), edm::Event::getByToken(), hcalSimParameters_cfi::he, VarParsing::mult, edm::EventBase::orbitNumber(), and EventWithHistoryFilter::selected().

                                                                                               {
   using namespace edm;
 
   // get Phase
 
   Handle<APVCyclePhaseCollection> apvphase;
   iEvent.getByToken(_apvphasecollToken, apvphase);
 
   // get HE
 
   Handle<EventWithHistory> he;
   iEvent.getByToken(_hecollectionToken, he);
 
   // check if the event is selected
 
   if ((_loworbit < 0 || iEvent.orbitNumber() >= _loworbit) && (_highorbit < 0 || iEvent.orbitNumber() <= _highorbit)) {
     if (_evfilter.selected(iEvent, iSetup)) {
       //Compute digi multiplicity
       /*
       int ntkdigi=0;
       int ntibdigi=0;
       int ntiddigi=0;
       int ntobdigi=0;
       int ntecdigi=0;
       */
       Handle<std::map<unsigned int, int> > mults;
       iEvent.getByToken(_multiplicityMapToken, mults);
 
       // create map of digi multiplicity
 
       std::map<int, int> digimap;
       for (std::map<unsigned int, int>::const_iterator mult = mults->begin(); mult != mults->end(); ++mult) {
         if (_subdets.find(mult->first) != _subdets.end())
           digimap[int(mult->first)] = mult->second;
       }
 
       _digibxcorrhmevent.fill(*he, digimap, apvphase);
 
       // fill debug histos
 
       if (he->depth() != 0) {
         long long dbx = he->deltaBX();
 
         if (_dbxhistos.find(dbx) != _dbxhistos.end()) {
           for (std::map<unsigned int, int>::const_iterator ndigi = mults->begin(); ndigi != mults->end(); ++ndigi) {
             _dbxhistos[dbx][ndigi->first]->Fill(ndigi->second);
           }
         }
         if (_dbxhistos.find(-1) != _dbxhistos.end()) {
           for (std::map<unsigned int, int>::const_iterator ndigi = mults->begin(); ndigi != mults->end(); ++ndigi) {
             _dbxhistos[-1][ndigi->first]->Fill(ndigi->second);
           }
         }
         /*
     if(_dbxtkhistos.find(dbx)!=_dbxtkhistos.end()) {
       _dbxtkhistos[dbx]->Fill(ntkdigi);
     }
     if(_dbxtkhistos.find(-1)!=_dbxtkhistos.end()) {
       _dbxtkhistos[-1]->Fill(ntkdigi);
     }
 
     if(_dbxtibhistos.find(dbx)!=_dbxtibhistos.end()) {
       _dbxtibhistos[dbx]->Fill(ntibdigi);
     }
     if(_dbxtibhistos.find(-1)!=_dbxtibhistos.end()) {
       _dbxtibhistos[-1]->Fill(ntibdigi);
     }
 
     if(_dbxtidhistos.find(dbx)!=_dbxtidhistos.end()) {
       _dbxtidhistos[dbx]->Fill(ntiddigi);
     }
     if(_dbxtidhistos.find(-1)!=_dbxtidhistos.end()) {
       _dbxtidhistos[-1]->Fill(ntiddigi);
     }
     if(_dbxtobhistos.find(dbx)!=_dbxtobhistos.end()) {
       _dbxtobhistos[dbx]->Fill(ntobdigi);
     }
     if(_dbxtobhistos.find(-1)!=_dbxtobhistos.end()) {
       _dbxtobhistos[-1]->Fill(ntobdigi);
     }
     if(_dbxtechistos.find(dbx)!=_dbxtechistos.end()) {
       _dbxtechistos[dbx]->Fill(ntecdigi);
     }
     if(_dbxtechistos.find(-1)!=_dbxtechistos.end()) {
       _dbxtechistos[-1]->Fill(ntecdigi);
     }
     */
       }
     }
   }
 }

void MultiplicityTimeCorrelations::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 335 of file MultiplicityTimeCorrelations.cc.

References _highorbit, _loworbit, and LogDebug.

                                             {
   LogDebug("IntegerDebug") << " int max and min " << INT_MIN << " " << INT_MAX;
   LogDebug("IntegerDebug") << " uint max and min " << UINT_MAX;
   LogDebug("IntegerDebug") << " long max and min " << LONG_MIN << " " << LONG_MAX;
   LogDebug("IntegerDebug") << " ulong max and min " << ULONG_MAX;
   LogDebug("IntegerDebug") << " long long max and min " << LLONG_MIN << " " << LLONG_MAX;
   LogDebug("IntegerDebug") << " u long long max and min " << ULLONG_MAX;
 
   edm::LogInfo("MultiplicityTimeCorrelations")
       << " Correlation studies performed only in the orbit # range " << _loworbit << " " << _highorbit;
 }

void MultiplicityTimeCorrelations::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 347 of file MultiplicityTimeCorrelations.cc.

References _digibxcorrhmevent, DigiBXCorrHistogramMaker< T >::beginRun(), and edm::RunBase::run().

                                                                                            {
   _digibxcorrhmevent.beginRun(iRun.run());
 }

void MultiplicityTimeCorrelations::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 351 of file MultiplicityTimeCorrelations.cc.

References DEFINE_FWK_MODULE.

351 {}

Member Data Documentation

edm::EDGetTokenT<APVCyclePhaseCollection> MultiplicityTimeCorrelations::_apvphasecollToken

private

Definition at line 87 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze().

std::map<unsigned int, int> MultiplicityTimeCorrelations::_binmax

private

Definition at line 90 of file MultiplicityTimeCorrelations.cc.

Referenced by MultiplicityTimeCorrelations().

std::vector<int> MultiplicityTimeCorrelations::_dbxbins

private

Definition at line 100 of file MultiplicityTimeCorrelations.cc.

Referenced by MultiplicityTimeCorrelations().

std::map<int, std::map<unsigned int, TH1F*> > MultiplicityTimeCorrelations::_dbxhistos

private

Definition at line 76 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze(), and MultiplicityTimeCorrelations().

DigiBXCorrHistogramMaker<EventWithHistory> MultiplicityTimeCorrelations::_digibxcorrhmevent

private

Definition at line 73 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze(), beginRun(), and MultiplicityTimeCorrelations().

EventWithHistoryFilter MultiplicityTimeCorrelations::_evfilter

private

Definition at line 74 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze(), and MultiplicityTimeCorrelations().

edm::InputTag MultiplicityTimeCorrelations::_hecollection

private

Definition at line 85 of file MultiplicityTimeCorrelations.cc.

Referenced by MultiplicityTimeCorrelations().

edm::EDGetTokenT<EventWithHistory> MultiplicityTimeCorrelations::_hecollectionToken

private

Definition at line 86 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze().

int MultiplicityTimeCorrelations::_highorbit

private

Definition at line 93 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze(), and beginJob().

int MultiplicityTimeCorrelations::_loworbit

private

Definition at line 92 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze(), and beginJob().

int MultiplicityTimeCorrelations::_mindbx

private

Definition at line 95 of file MultiplicityTimeCorrelations.cc.

Referenced by MultiplicityTimeCorrelations().

int MultiplicityTimeCorrelations::_mintrpltdbx

private

Definition at line 96 of file MultiplicityTimeCorrelations.cc.

Referenced by MultiplicityTimeCorrelations().

edm::EDGetTokenT<std::map<unsigned int, int> > MultiplicityTimeCorrelations::_multiplicityMapToken

private

Definition at line 88 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze().

std::map<unsigned int, std::string> MultiplicityTimeCorrelations::_subdets

private

Definition at line 89 of file MultiplicityTimeCorrelations.cc.

Referenced by analyze(), and MultiplicityTimeCorrelations().

SiStripTKNumbers MultiplicityTimeCorrelations::_trnumb

private

Definition at line 98 of file MultiplicityTimeCorrelations.cc.

Referenced by MultiplicityTimeCorrelations().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation