Inheritance diagram for AnalyzerMinbias:

Classes
struct	myInfo

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

	AnalyzerMinbias (const edm::ParameterSet &)

void	beginJob () override

void	endJob () override

	~AnalyzerMinbias () 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 Member Functions
void	analyzeHcal (const HcalRespCorrs *myRecalib, const HBHERecHitCollection &HithbheNS, const HBHERecHitCollection &HithbheMB, const HFRecHitCollection &HithfNS, const HFRecHitCollection &HithfMB, int algoBit, bool fill)

Private Attributes
int	cells

int	depth

float	eta

std::string	fOutputFileName

TH1D *	h_Noise [4]

TH1D *	h_Signal [4]

std::string	hcalfile_

TFile *	hOutputFile

int	ieta

bool	ignoreL1_

int	iphi

float	mom0_Diff

float	mom0_MB

float	mom0_Noise

float	mom1_Diff

float	mom1_MB

float	mom1_Noise

float	mom2_Diff

float	mom2_MB

float	mom2_Noise

float	mom3_Diff

float	mom3_MB

float	mom3_Noise

float	mom4_Diff

float	mom4_MB

float	mom4_Noise

int	mydet

std::map< std::pair< int, HcalDetId >, myInfo >	myMap_

std::ofstream *	myout_hcal

int	mysubd

TTree *	myTree

float	occup

float	phi

double	rnnum

double	rnnumber

bool	runNZS_

bool	theRecalib_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbheNormal_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbherecoMB_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbherecoNoise_

edm::EDGetTokenT< HFRecHitCollection >	tok_hfrecoMB_

edm::EDGetTokenT< HFRecHitCollection >	tok_hfrecoNoise_

edm::EDGetTokenT< L1GlobalTriggerObjectMapRecord >	tok_hltL1GtMap_

edm::EDGetTokenT< HORecHitCollection >	tok_horecoMB_

edm::EDGetTokenT< HORecHitCollection >	tok_horecoNoise_

edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_respCorr_

int	trigbit

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 59 of file AnalyzerMinbias.cc.

Constructor & Destructor Documentation

◆ AnalyzerMinbias()

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

explicit

Definition at line 114 of file AnalyzerMinbias.cc.

                                                                {
   // get name of output file with histogramms
   fOutputFileName = iConfig.getUntrackedParameter<std::string>("HistOutFile");
  
   // get token names of modules, producing object collections
   tok_hbherecoMB_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"));
   tok_horecoMB_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputMB"));
   tok_hfrecoMB_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"));
  
   tok_hbherecoNoise_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputNoise"));
   tok_horecoNoise_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputNoise"));
   tok_hfrecoNoise_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputNoise"));
  
   theRecalib_ = iConfig.getParameter<bool>("Recalib");
   ignoreL1_ = iConfig.getUntrackedParameter<bool>("IgnoreL1", true);
   runNZS_ = iConfig.getUntrackedParameter<bool>("RunNZS", true);
  
   tok_hbheNormal_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
   tok_hltL1GtMap_ = consumes<L1GlobalTriggerObjectMapRecord>(edm::InputTag("hltL1GtObjectMap"));
  
   tok_respCorr_ = esConsumes<HcalRespCorrs, HcalRespCorrsRcd>();
 }

References fOutputFileName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ignoreL1_, runNZS_, AlCaHLTBitMon_QueryRunRegistry::string, theRecalib_, tok_hbheNormal_, tok_hbherecoMB_, tok_hbherecoNoise_, tok_hfrecoMB_, tok_hfrecoNoise_, tok_hltL1GtMap_, tok_horecoMB_, tok_horecoNoise_, and tok_respCorr_.

◆ ~AnalyzerMinbias()

AnalyzerMinbias::~AnalyzerMinbias ( )

override

Definition at line 137 of file AnalyzerMinbias.cc.

137 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 242 of file AnalyzerMinbias.cc.

                                                                                  {
   rnnum = (float)iEvent.run();
   const HcalRespCorrs* myRecalib = nullptr;
   if (theRecalib_) {
     myRecalib = &iSetup.getData(tok_respCorr_);
   }  // theRecalib
  
   edm::Handle<HBHERecHitCollection> hbheNormal;
   iEvent.getByToken(tok_hbheNormal_, hbheNormal);
   if (!hbheNormal.isValid()) {
     edm::LogVerbatim("AnalyzerMB") << " hbheNormal failed";
   } else {
     edm::LogVerbatim("AnalyzerMB") << " The size of the normal collection " << hbheNormal->size();
   }
  
   edm::Handle<HBHERecHitCollection> hbheNS;
   iEvent.getByToken(tok_hbherecoNoise_, hbheNS);
   if (!hbheNS.isValid()) {
     edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbheNoise product!";
     return;
   }
   const HBHERecHitCollection HithbheNS = *(hbheNS.product());
   edm::LogVerbatim("AnalyzerMB") << "HBHE NS size of collection " << HithbheNS.size();
   if (runNZS_ && HithbheNS.size() != 5184) {
     edm::LogWarning("AnalyzerMB") << "HBHE NS problem " << rnnum << " size " << HithbheNS.size();
     return;
   }
  
   edm::Handle<HBHERecHitCollection> hbheMB;
   iEvent.getByToken(tok_hbherecoMB_, hbheMB);
   if (!hbheMB.isValid()) {
     edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe product!";
     return;
   }
   const HBHERecHitCollection HithbheMB = *(hbheMB.product());
   edm::LogVerbatim("AnalyzerMB") << "HBHE MB size of collection " << HithbheMB.size();
   if (runNZS_ && HithbheMB.size() != 5184) {
     edm::LogWarning("AnalyzerMB") << "HBHE problem " << rnnum << " size " << HithbheMB.size();
     return;
   }
  
   edm::Handle<HFRecHitCollection> hfNS;
   iEvent.getByToken(tok_hfrecoNoise_, hfNS);
   if (!hfNS.isValid()) {
     edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hfNoise product!";
     return;
   }
   const HFRecHitCollection HithfNS = *(hfNS.product());
   edm::LogVerbatim("AnalyzerMB") << "HF NS size of collection " << HithfNS.size();
   if (runNZS_ && HithfNS.size() != 1728) {
     edm::LogWarning("AnalyzerMB") << "HF NS problem " << rnnum << " size " << HithfNS.size();
     return;
   }
  
   edm::Handle<HFRecHitCollection> hfMB;
   iEvent.getByToken(tok_hfrecoMB_, hfMB);
   if (!hfMB.isValid()) {
     edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf product!";
     return;
   }
   const HFRecHitCollection HithfMB = *(hfMB.product());
   edm::LogVerbatim("AnalyzerMB") << "HF MB size of collection " << HithfMB.size();
   if (runNZS_ && HithfMB.size() != 1728) {
     edm::LogWarning("AnalyzerMB") << "HF problem " << rnnum << " size " << HithfMB.size();
     return;
   }
  
   if (ignoreL1_) {
     analyzeHcal(myRecalib, HithbheNS, HithbheMB, HithfNS, HithfMB, 1, true);
   } else {
     edm::Handle<L1GlobalTriggerObjectMapRecord> gtObjectMapRecord;
     iEvent.getByToken(tok_hltL1GtMap_, gtObjectMapRecord);
     if (gtObjectMapRecord.isValid()) {
       const std::vector<L1GlobalTriggerObjectMap>& objMapVec = gtObjectMapRecord->gtObjectMap();
       int ii(0);
       bool ok(false), fill(true);
       for (std::vector<L1GlobalTriggerObjectMap>::const_iterator itMap = objMapVec.begin(); itMap != objMapVec.end();
            ++itMap, ++ii) {
         bool resultGt = (*itMap).algoGtlResult();
         if (resultGt == 1) {
           ok = true;
           int algoBit = (*itMap).algoBitNumber();
           analyzeHcal(myRecalib, HithbheNS, HithbheMB, HithfNS, HithfMB, algoBit, fill);
           fill = false;
           std::string algoNameStr = (*itMap).algoName();
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("AnalyzerMB") << "Trigger[" << ii << "] " << algoNameStr << " bit " << algoBit << " entered";
 #endif
         }
       }
       if (!ok)
         edm::LogVerbatim("AnalyzerMB") << "No passed L1 Triggers";
     }
   }
 }

References analyzeHcal(), ntuplemaker::fill, dqmMemoryStats::float, edm::EventSetup::getData(), L1GlobalTriggerObjectMapRecord::gtObjectMap(), iEvent, ignoreL1_, cuy::ii, edm::HandleBase::isValid(), convertSQLiteXML::ok, edm::Handle< T >::product(), rnnum, runNZS_, edm::SortedCollection< T, SORT >::size(), AlCaHLTBitMon_QueryRunRegistry::string, theRecalib_, tok_hbheNormal_, tok_hbherecoMB_, tok_hbherecoNoise_, tok_hfrecoMB_, tok_hfrecoNoise_, tok_hltL1GtMap_, and tok_respCorr_.

◆ analyzeHcal()

void AnalyzerMinbias::analyzeHcal	(	const HcalRespCorrs *	myRecalib,
		const HBHERecHitCollection &	HithbheNS,
		const HBHERecHitCollection &	HithbheMB,
		const HFRecHitCollection &	HithfNS,
		const HFRecHitCollection &	HithfMB,
		int	algoBit,
		bool	fill
	)

private

Definition at line 338 of file AnalyzerMinbias.cc.

                                              {
   // Noise part for HB HE
   std::map<std::pair<int, HcalDetId>, myInfo> tmpMap;
   tmpMap.clear();
  
   for (HBHERecHitCollection::const_iterator hbheItr = HithbheNS.begin(); hbheItr != HithbheNS.end(); hbheItr++) {
     // Recalibration of energy
     float icalconst = 1.;
     DetId mydetid = hbheItr->id().rawId();
     if (theRecalib_)
       icalconst = myRecalib->getValues(mydetid)->getValue();
  
     HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
     double energyhit = aHit.energy();
  
     DetId id = (*hbheItr).detid();
     HcalDetId hid = HcalDetId(id);
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
     if (itr1 == myMap_.end()) {
       myInfo info;
       myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
       itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
     }
     itr1->second.theNS0++;
     itr1->second.theNS1 += energyhit;
     itr1->second.theNS2 += (energyhit * energyhit);
     itr1->second.theNS3 += (energyhit * energyhit * energyhit);
     itr1->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
     itr1->second.runcheck = rnnum;
     if (fill)
       h_Noise[hid.subdet() - 1]->Fill(energyhit);
  
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
     if (itr2 == tmpMap.end()) {
       myInfo info;
       tmpMap[std::pair<int, HcalDetId>(algoBit, hid)] = info;
       itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
     }
     itr2->second.theNS0++;
     itr2->second.theNS1 += energyhit;
     itr2->second.theNS2 += (energyhit * energyhit);
     itr2->second.theNS3 += (energyhit * energyhit * energyhit);
     itr2->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
     itr2->second.runcheck = rnnum;
  
   }  // HBHE_NS
  
   // Signal part for HB HE
  
   for (HBHERecHitCollection::const_iterator hbheItr = HithbheMB.begin(); hbheItr != HithbheMB.end(); hbheItr++) {
     // Recalibration of energy
     float icalconst = 1.;
     DetId mydetid = hbheItr->id().rawId();
     if (theRecalib_)
       icalconst = myRecalib->getValues(mydetid)->getValue();
  
     HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
     double energyhit = aHit.energy();
  
     DetId id = (*hbheItr).detid();
     HcalDetId hid = HcalDetId(id);
  
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
  
     if (itr1 == myMap_.end()) {
       myInfo info;
       myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
       itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
     }
     itr1->second.theMB0++;
     itr1->second.theDif0 = 0;
     itr1->second.theMB1 += energyhit;
     itr1->second.theMB2 += (energyhit * energyhit);
     itr1->second.theMB3 += (energyhit * energyhit * energyhit);
     itr1->second.theMB4 += (energyhit * energyhit * energyhit * energyhit);
     itr1->second.runcheck = rnnum;
     float mydiff = 0.0;
     if (itr2 != tmpMap.end()) {
       mydiff = energyhit - (itr2->second.theNS1);
       itr1->second.theDif0++;
       itr1->second.theDif1 += mydiff;
       itr1->second.theDif2 += (mydiff * mydiff);
       if (fill)
         h_Signal[hid.subdet() - 1]->Fill(mydiff);
     }
   }  // HBHE_MB
  
   // HF
  
   for (HFRecHitCollection::const_iterator hbheItr = HithfNS.begin(); hbheItr != HithfNS.end(); hbheItr++) {
     // Recalibration of energy
     float icalconst = 1.;
     DetId mydetid = hbheItr->id().rawId();
     if (theRecalib_)
       icalconst = myRecalib->getValues(mydetid)->getValue();
  
     HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
     double energyhit = aHit.energy();
     // Remove PMT hits
     if (fabs(energyhit) > 40.)
       continue;
     DetId id = (*hbheItr).detid();
     HcalDetId hid = HcalDetId(id);
  
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
  
     if (itr1 == myMap_.end()) {
       myInfo info;
       myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
       itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
     }
     itr1->second.theNS0++;
     itr1->second.theNS1 += energyhit;
     itr1->second.theNS2 += (energyhit * energyhit);
     itr1->second.theNS3 += (energyhit * energyhit * energyhit);
     itr1->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
     itr1->second.runcheck = rnnum;
     if (fill)
       h_Noise[hid.subdet() - 1]->Fill(energyhit);
  
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
     if (itr2 == tmpMap.end()) {
       myInfo info;
       tmpMap[std::pair<int, HcalDetId>(algoBit, hid)] = info;
       itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
     }
     itr2->second.theNS0++;
     itr2->second.theNS1 += energyhit;
     itr2->second.theNS2 += (energyhit * energyhit);
     itr2->second.theNS3 += (energyhit * energyhit * energyhit);
     itr2->second.theNS4 += (energyhit * energyhit * energyhit * energyhit);
     itr2->second.runcheck = rnnum;
  
   }  // HF_NS
  
   // Signal part for HF
  
   for (HFRecHitCollection::const_iterator hbheItr = HithfMB.begin(); hbheItr != HithfMB.end(); hbheItr++) {
     // Recalibration of energy
     float icalconst = 1.;
     DetId mydetid = hbheItr->id().rawId();
     if (theRecalib_)
       icalconst = myRecalib->getValues(mydetid)->getValue();
     HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
  
     double energyhit = aHit.energy();
     // Remove PMT hits
     if (fabs(energyhit) > 40.)
       continue;
  
     DetId id = (*hbheItr).detid();
     HcalDetId hid = HcalDetId(id);
  
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
     std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr2 = tmpMap.find(std::pair<int, HcalDetId>(algoBit, hid));
  
     if (itr1 == myMap_.end()) {
       myInfo info;
       myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
       itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
     }
     itr1->second.theMB0++;
     itr1->second.theDif0 = 0;
     itr1->second.theMB1 += energyhit;
     itr1->second.theMB2 += (energyhit * energyhit);
     itr1->second.theMB3 += (energyhit * energyhit * energyhit);
     itr1->second.theMB4 += (energyhit * energyhit * energyhit * energyhit);
     itr1->second.runcheck = rnnum;
     float mydiff = 0.0;
     if (itr2 != tmpMap.end()) {
       mydiff = energyhit - (itr2->second.theNS1);
       itr1->second.theDif0++;
       itr1->second.theDif1 += mydiff;
       itr1->second.theDif2 += (mydiff * mydiff);
       if (fill)
         h_Signal[hid.subdet() - 1]->Fill(mydiff);
     }
   }
 }

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), ntuplemaker::fill, HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), h_Noise, h_Signal, info(), myMap_, DetId::rawId(), rnnum, HcalDetId::subdet(), and theRecalib_.

Referenced by analyze().

◆ beginJob()

void AnalyzerMinbias::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 139 of file AnalyzerMinbias.cc.

                                {
   std::string det[4] = {"HB", "HE", "HO", "HF"};
   char name[80], title[80];
   for (int subd = 0; subd < 4; ++subd) {
     sprintf(name, "Noise_%s", det[subd].c_str());
     sprintf(title, "Energy Distribution for Noise in %s", det[subd].c_str());
     h_Noise[subd] = new TH1D(name, title, 100, -10., 10.);
     sprintf(name, "Signal_%s", det[subd].c_str());
     sprintf(title, "Energy Distribution for Signal in %s", det[subd].c_str());
     h_Signal[subd] = new TH1D(name, title, 100, -10., 10.);
   }
  
   hOutputFile = new TFile(fOutputFileName.c_str(), "RECREATE");
   myTree = new TTree("RecJet", "RecJet Tree");
   myTree->Branch("mydet", &mydet, "mydet/I");
   myTree->Branch("mysubd", &mysubd, "mysubd/I");
   myTree->Branch("cells", &cells, "cells");
   myTree->Branch("depth", &depth, "depth/I");
   myTree->Branch("ieta", &ieta, "ieta/I");
   myTree->Branch("iphi", &iphi, "iphi/I");
   myTree->Branch("eta", &eta, "eta/F");
   myTree->Branch("phi", &phi, "phi/F");
   myTree->Branch("mom0_MB", &mom0_MB, "mom0_MB/F");
   myTree->Branch("mom1_MB", &mom1_MB, "mom1_MB/F");
   myTree->Branch("mom2_MB", &mom2_MB, "mom2_MB/F");
   myTree->Branch("mom3_MB", &mom3_MB, "mom3_MB/F");
   myTree->Branch("mom4_MB", &mom4_MB, "mom4_MB/F");
   myTree->Branch("mom0_Noise", &mom0_Noise, "mom0_Noise/F");
   myTree->Branch("mom1_Noise", &mom1_Noise, "mom1_Noise/F");
   myTree->Branch("mom2_Noise", &mom2_Noise, "mom2_Noise/F");
   myTree->Branch("mom3_Noise", &mom3_Noise, "mom3_Noise/F");
   myTree->Branch("mom4_Noise", &mom4_Noise, "mom4_Noise/F");
   myTree->Branch("mom0_Diff", &mom0_Diff, "mom0_Diff/F");
   myTree->Branch("mom1_Diff", &mom1_Diff, "mom1_Diff/F");
   myTree->Branch("mom2_Diff", &mom2_Diff, "mom2_Diff/F");
   myTree->Branch("occup", &occup, "occup/F");
   myTree->Branch("trigbit", &trigbit, "trigbit/I");
   myTree->Branch("rnnumber", &rnnumber, "rnnumber/D");
  
   myMap_.clear();
 }

References cells, depth, eta, fOutputFileName, h_Noise, h_Signal, hOutputFile, ieta, iphi, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom3_MB, mom3_Noise, mom4_MB, mom4_Noise, mydet, myMap_, mysubd, myTree, Skims_PA_cff::name, occup, phi, rnnumber, AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, and trigbit.

◆ endJob()

void AnalyzerMinbias::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 183 of file AnalyzerMinbias.cc.

                              {
   int ii = 0;
   for (std::map<std::pair<int, HcalDetId>, myInfo>::const_iterator itr = myMap_.begin(); itr != myMap_.end(); ++itr) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("AnalyzerMB") << "Fired trigger bit number " << itr->first.first;
 #endif
     myInfo info = itr->second;
     if (info.theMB0 > 0) {
       mom0_MB = info.theMB0;
       mom1_MB = info.theMB1;
       mom2_MB = info.theMB2;
       mom3_MB = info.theMB3;
       mom4_MB = info.theMB4;
       mom0_Noise = info.theNS0;
       mom1_Noise = info.theNS1;
       mom2_Noise = info.theNS2;
       mom3_Noise = info.theNS3;
       mom4_Noise = info.theNS4;
       mom0_Diff = info.theDif0;
       mom1_Diff = info.theDif1;
       mom2_Diff = info.theDif2;
       rnnumber = info.runcheck;
       trigbit = itr->first.first;
       mysubd = itr->first.second.subdet();
       depth = itr->first.second.depth();
       ieta = itr->first.second.ieta();
       iphi = itr->first.second.iphi();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("AnalyzerMB") << " Result=  " << trigbit << " " << mysubd << " " << ieta << " " << iphi
                                      << " mom0  " << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB << " mom3 "
                                      << mom3_MB << " mom4 " << mom4_MB << " mom0_Noise " << mom0_Noise << " mom1_Noise "
                                      << mom1_Noise << " mom2_Noise " << mom2_Noise << " mom3_Noise " << mom3_Noise
                                      << " mom4_Noise " << mom4_Noise << " mom0_Diff " << mom0_Diff << " mom1_Diff "
                                      << mom1_Diff << " mom2_Diff " << mom2_Diff;
 #endif
       myTree->Fill();
       ii++;
     }
   }
   cells = ii;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("AnalyzerMB") << "cells " << cells;
 #endif
   hOutputFile->Write();
   hOutputFile->cd();
   myTree->Write();
   for (int i = 0; i < 4; i++) {
     h_Noise[i]->Write();
     h_Signal[i]->Write();
   }
   hOutputFile->Close();
 }

References cells, depth, h_Noise, h_Signal, hOutputFile, mps_fire::i, ieta, cuy::ii, info(), iphi, genParticles_cff::map, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom3_MB, mom3_Noise, mom4_MB, mom4_Noise, myMap_, mysubd, myTree, rnnumber, and trigbit.

Referenced by o2olib.O2ORunMgr::executeJob().

Member Data Documentation

◆ cells

int AnalyzerMinbias::cells

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ depth

int AnalyzerMinbias::depth

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by dqm_interfaces.DirID::__eq__(), dqm_interfaces.DirID::__repr__(), beginJob(), endJob(), and dqm_interfaces.DirWalkerDB::run().

◆ eta

float AnalyzerMinbias::eta

private

Definition at line 100 of file AnalyzerMinbias.cc.

Referenced by Particle.Particle::__str__(), beginJob(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

◆ fOutputFileName

std::string AnalyzerMinbias::fOutputFileName

private

Definition at line 89 of file AnalyzerMinbias.cc.

Referenced by AnalyzerMinbias(), and beginJob().

◆ h_Noise

TH1D* AnalyzerMinbias::h_Noise[4]

private

Definition at line 93 of file AnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and endJob().

◆ h_Signal

TH1D * AnalyzerMinbias::h_Signal[4]

private

Definition at line 93 of file AnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and endJob().

◆ hcalfile_

std::string AnalyzerMinbias::hcalfile_

private

Definition at line 89 of file AnalyzerMinbias.cc.

◆ hOutputFile

TFile* AnalyzerMinbias::hOutputFile

private

Definition at line 91 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ ieta

int AnalyzerMinbias::ieta

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ ignoreL1_

bool AnalyzerMinbias::ignoreL1_

private

Definition at line 94 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ iphi

int AnalyzerMinbias::iphi

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom0_Diff

float AnalyzerMinbias::mom0_Diff

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom0_MB

float AnalyzerMinbias::mom0_MB

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom0_Noise

float AnalyzerMinbias::mom0_Noise

private

Definition at line 102 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom1_Diff

float AnalyzerMinbias::mom1_Diff

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom1_MB

float AnalyzerMinbias::mom1_MB

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom1_Noise

float AnalyzerMinbias::mom1_Noise

private

Definition at line 102 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom2_Diff

float AnalyzerMinbias::mom2_Diff

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom2_MB

float AnalyzerMinbias::mom2_MB

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom2_Noise

float AnalyzerMinbias::mom2_Noise

private

Definition at line 102 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom3_Diff

float AnalyzerMinbias::mom3_Diff

private

Definition at line 103 of file AnalyzerMinbias.cc.

◆ mom3_MB

float AnalyzerMinbias::mom3_MB

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom3_Noise

float AnalyzerMinbias::mom3_Noise

private

Definition at line 102 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom4_Diff

float AnalyzerMinbias::mom4_Diff

private

Definition at line 103 of file AnalyzerMinbias.cc.

◆ mom4_MB

float AnalyzerMinbias::mom4_MB

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom4_Noise

float AnalyzerMinbias::mom4_Noise

private

Definition at line 102 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mydet

int AnalyzerMinbias::mydet

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by beginJob().

◆ myMap_

std::map<std::pair<int, HcalDetId>, myInfo> AnalyzerMinbias::myMap_

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and endJob().

◆ myout_hcal

std::ofstream* AnalyzerMinbias::myout_hcal

private

Definition at line 90 of file AnalyzerMinbias.cc.

◆ mysubd

int AnalyzerMinbias::mysubd

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ myTree

TTree* AnalyzerMinbias::myTree

private

Definition at line 92 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ occup

float AnalyzerMinbias::occup

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by beginJob().

◆ phi

float AnalyzerMinbias::phi

private

Definition at line 100 of file AnalyzerMinbias.cc.

Referenced by Particle.Particle::__str__(), beginJob(), and ntupleDataFormat.Track::phiPull().

◆ rnnum

double AnalyzerMinbias::rnnum

private

Definition at line 95 of file AnalyzerMinbias.cc.

Referenced by analyze(), and analyzeHcal().

◆ rnnumber

double AnalyzerMinbias::rnnumber

private

Definition at line 98 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ runNZS_

bool AnalyzerMinbias::runNZS_

private

Definition at line 94 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ theRecalib_

bool AnalyzerMinbias::theRecalib_

private

Definition at line 94 of file AnalyzerMinbias.cc.

Referenced by analyze(), analyzeHcal(), and AnalyzerMinbias().

◆ tok_hbheNormal_

edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbheNormal_

private

Definition at line 109 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ tok_hbherecoMB_

edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbherecoMB_

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ tok_hbherecoNoise_

edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbherecoNoise_

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ tok_hfrecoMB_

edm::EDGetTokenT<HFRecHitCollection> AnalyzerMinbias::tok_hfrecoMB_

private

Definition at line 107 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ tok_hfrecoNoise_

edm::EDGetTokenT<HFRecHitCollection> AnalyzerMinbias::tok_hfrecoNoise_

private

Definition at line 107 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ tok_hltL1GtMap_

edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> AnalyzerMinbias::tok_hltL1GtMap_

private

Definition at line 110 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ tok_horecoMB_

edm::EDGetTokenT<HORecHitCollection> AnalyzerMinbias::tok_horecoMB_

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by AnalyzerMinbias().

◆ tok_horecoNoise_

edm::EDGetTokenT<HORecHitCollection> AnalyzerMinbias::tok_horecoNoise_

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by AnalyzerMinbias().

◆ tok_respCorr_

edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> AnalyzerMinbias::tok_respCorr_

private

Definition at line 111 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

◆ trigbit

int AnalyzerMinbias::trigbit

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AnalyzerMinbias()

◆ ~AnalyzerMinbias()

Member Function Documentation

◆ analyze()

◆ analyzeHcal()

◆ beginJob()

◆ endJob()

Member Data Documentation

◆ cells

◆ depth

◆ eta

◆ fOutputFileName

◆ h_Noise

◆ h_Signal

◆ hcalfile_

◆ hOutputFile

◆ ieta

◆ ignoreL1_

◆ iphi

◆ mom0_Diff

◆ mom0_MB

◆ mom0_Noise

◆ mom1_Diff

◆ mom1_MB

◆ mom1_Noise

◆ mom2_Diff

◆ mom2_MB

◆ mom2_Noise

◆ mom3_Diff

◆ mom3_MB

◆ mom3_Noise

◆ mom4_Diff

◆ mom4_MB

◆ mom4_Noise

◆ mydet

◆ myMap_

◆ myout_hcal

◆ mysubd

◆ myTree

◆ occup

◆ phi

◆ rnnum

◆ rnnumber

◆ runNZS_

◆ theRecalib_

◆ tok_hbheNormal_

◆ tok_hbherecoMB_

◆ tok_hbherecoNoise_

◆ tok_hfrecoMB_

◆ tok_hfrecoNoise_

◆ tok_hltL1GtMap_

◆ tok_horecoMB_

◆ tok_horecoNoise_

◆ tok_respCorr_

◆ trigbit