#include <MinBias.h>

Inheritance diagram for cms::MinBias:

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

void	beginJob () override

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

void	endJob () override

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

	MinBias (const edm::ParameterSet &)

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
bool	allowMissingInputs_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryESToken_

int	depth

float	eta

const CaloGeometry *	geo_

std::string	hbheLabel_

edm::EDGetTokenT< HBHERecHitCollection >	hbheToken_

std::string	hfLabel_

edm::EDGetTokenT< HFRecHitCollection >	hfToken_

std::string	hoLabel_

edm::EDGetTokenT< HORecHitCollection >	hoToken_

int	ieta

int	iphi

float	mom1

float	mom2

float	mom3

float	mom4

int	mydet

int	mysubd

TTree *	myTree_

float	occup

float	phi

std::map< DetId, double >	theFillDetMap0_

std::map< DetId, double >	theFillDetMap1_

std::map< DetId, double >	theFillDetMap2_

std::map< DetId, double >	theFillDetMap3_

std::map< DetId, double >	theFillDetMap4_

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 24 of file MinBias.h.

Constructor & Destructor Documentation

◆ MinBias()

cms::MinBias::MinBias ( const edm::ParameterSet & iConfig )

explicit

Definition at line 9 of file MinBias.cc.

References allowMissingInputs_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hbheLabel_, hbheToken_, hfLabel_, hfToken_, hoLabel_, hoToken_, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

       : caloGeometryESToken_(esConsumes<edm::Transition::BeginRun>()), geo_(nullptr) {
     // get names of modules, producing object collections
     hbheLabel_ = iConfig.getParameter<std::string>("hbheInput");
     hoLabel_ = iConfig.getParameter<std::string>("hoInput");
     hfLabel_ = iConfig.getParameter<std::string>("hfInput");
     hbheToken_ = mayConsume<HBHERecHitCollection>(edm::InputTag(hbheLabel_));
     hoToken_ = mayConsume<HORecHitCollection>(edm::InputTag(hoLabel_));
     hfToken_ = mayConsume<HFRecHitCollection>(edm::InputTag(hfLabel_));
     allowMissingInputs_ = iConfig.getUntrackedParameter<bool>("AllowMissingInputs", false);
   }

Member Function Documentation

◆ analyze()

void cms::MinBias::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 108 of file MinBias.cc.

References allowMissingInputs_, photonIsolationHIProducer_cfi::hbhe, hbheLabel_, hbheToken_, photonIsolationHIProducer_cfi::hf, hfLabel_, hfToken_, photonIsolationHIProducer_cfi::ho, hoLabel_, hoToken_, l1ctLayer2EG_cff::id, iEvent, funct::pow(), theFillDetMap0_, theFillDetMap1_, theFillDetMap2_, theFillDetMap3_, and theFillDetMap4_.

                                                                     {
     if (!hbheLabel_.empty()) {
       edm::Handle<HBHERecHitCollection> hbhe;
       iEvent.getByToken(hbheToken_, hbhe);
       if (!hbhe.isValid()) {
         // can't find it!
         if (!allowMissingInputs_) {
           *hbhe;  // will throw the proper exception
         }
       } else {
         for (auto const& hbheItr : (HBHERecHitCollection)(*hbhe)) {
           DetId id = (hbheItr).detid();
           if (hbheItr.energy() > 0.)
             edm::LogWarning("MinBias") << " Energy = " << hbheItr.energy();
           theFillDetMap0_[id] += 1.;
           theFillDetMap1_[id] += hbheItr.energy();
           theFillDetMap2_[id] += pow(hbheItr.energy(), 2);
           theFillDetMap3_[id] += pow(hbheItr.energy(), 3);
           theFillDetMap4_[id] += pow(hbheItr.energy(), 4);
         }
       }
     }
 
     if (!hoLabel_.empty()) {
       edm::Handle<HORecHitCollection> ho;
       iEvent.getByToken(hoToken_, ho);
       if (!ho.isValid()) {
         // can't find it!
         if (!allowMissingInputs_) {
           *ho;  // will throw the proper exception
         }
       } else {
         for (auto const& hoItr : (HORecHitCollection)(*ho)) {
           DetId id = hoItr.detid();
           theFillDetMap0_[id] += 1.;
           theFillDetMap1_[id] += hoItr.energy();
           theFillDetMap2_[id] += pow(hoItr.energy(), 2);
           theFillDetMap3_[id] += pow(hoItr.energy(), 3);
           theFillDetMap4_[id] += pow(hoItr.energy(), 4);
         }
       }
     }
 
     if (!hfLabel_.empty()) {
       edm::Handle<HFRecHitCollection> hf;
       iEvent.getByToken(hfToken_, hf);
       if (!hf.isValid()) {
         // can't find it!
         if (!allowMissingInputs_) {
           *hf;  // will throw the proper exception
         }
       } else {
         for (auto const hfItr : (HFRecHitCollection)(*hf)) {
           DetId id = hfItr.detid();
           theFillDetMap0_[id] += 1.;
           theFillDetMap1_[id] += hfItr.energy();
           theFillDetMap2_[id] += pow(hfItr.energy(), 2);
           theFillDetMap3_[id] += pow(hfItr.energy(), 3);
           theFillDetMap4_[id] += pow(hfItr.energy(), 4);
         }
       }
     }
   }

◆ beginJob()

void cms::MinBias::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 21 of file MinBias.cc.

References depth, eta, compareTotals::fs, ieta, iphi, mom1, mom2, mom3, mom4, mydet, mysubd, myTree_, and phi.

                          {
     edm::Service<TFileService> fs;
     myTree_ = fs->make<TTree>("RecJet", "RecJet Tree");
     myTree_->Branch("mydet", &mydet, "mydet/I");
     myTree_->Branch("mysubd", &mysubd, "mysubd/I");
     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("mom1", &mom1, "mom1/F");
     myTree_->Branch("mom2", &mom2, "mom2/F");
     myTree_->Branch("mom3", &mom3, "mom3/F");
     myTree_->Branch("mom4", &mom4, "mom4/F");
   }

◆ beginRun()

void cms::MinBias::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	iSetup
	)

override

Definition at line 37 of file MinBias.cc.

References caloGeometryESToken_, geo_, edm::EventSetup::getData(), CaloGeometry::getValidDetIds(), DetId::Hcal, l1ctLayer2EG_cff::id, theFillDetMap0_, theFillDetMap1_, theFillDetMap2_, theFillDetMap3_, and theFillDetMap4_.

                                                                    {
     geo_ = &iSetup.getData(caloGeometryESToken_);
     std::vector<DetId> did = geo_->getValidDetIds();
 
     for (auto const& id : did) {
       if ((id).det() == DetId::Hcal) {
         theFillDetMap0_[id] = 0.;
         theFillDetMap1_[id] = 0.;
         theFillDetMap2_[id] = 0.;
         theFillDetMap3_[id] = 0.;
         theFillDetMap4_[id] = 0.;
       }
     }
   }

◆ endJob()

void cms::MinBias::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 54 of file MinBias.cc.

References depth, HcalDetId::depth(), eta, geo_, HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), HcalGeometry::getValidDetIds(), DetId::Hcal, mps_fire::i, l1ctLayer2EG_cff::id, ieta, HcalDetId::ieta(), createfilelist::int, iphi, HcalDetId::iphi(), mom1, mom2, mom3, mom4, mydet, mysubd, myTree_, phi, funct::pow(), theFillDetMap0_, theFillDetMap1_, theFillDetMap2_, theFillDetMap3_, and theFillDetMap4_.

                        {
     const HcalGeometry* hgeo = static_cast<const HcalGeometry*>(geo_->getSubdetectorGeometry(DetId::Hcal, 1));
     const std::vector<DetId>& did = hgeo->getValidDetIds();
     int i = 0;
     for (const auto& id : did) {
       //      if( id.det() == DetId::Hcal ) {
       GlobalPoint pos = hgeo->getPosition(id);
       mydet = (int)(id.det());
       mysubd = (id.subdetId());
       depth = HcalDetId(id).depth();
       ieta = HcalDetId(id).ieta();
       iphi = HcalDetId(id).iphi();
       phi = pos.phi();
       eta = pos.eta();
       if (theFillDetMap0_[id] > 0.) {
         mom1 = theFillDetMap1_[id] / theFillDetMap0_[id];
         mom2 = theFillDetMap2_[id] / theFillDetMap0_[id] - (mom1 * mom1);
         mom3 = theFillDetMap3_[id] / theFillDetMap0_[id] - 3. * mom1 * theFillDetMap2_[id] / theFillDetMap0_[id] +
                2. * pow(mom2, 3);
         mom4 = (theFillDetMap4_[id] - 4. * mom1 * theFillDetMap3_[id] + 6. * pow(mom1, 2) * theFillDetMap2_[id]) /
                    theFillDetMap0_[id] -
                3. * pow(mom1, 4);
 
       } else {
         mom1 = 0.;
         mom2 = 0.;
         mom3 = 0.;
         mom4 = 0.;
       }
       edm::LogWarning("MinBias") << " Detector " << id.rawId() << " mydet " << mydet << " " << mysubd << " " << depth
                                  << " " << ieta << " " << iphi << " " << pos.eta() << " " << pos.phi();
       edm::LogWarning("MinBias") << " Energy " << mom1 << " " << mom2 << std::endl;
       myTree_->Fill();
       i++;
       //      }
     }
     edm::LogWarning("MinBias") << " The number of CaloDet records " << did.size();
     edm::LogWarning("MinBias") << " The number of Hcal records " << i << std::endl;
 
     /*
   std::cout << "===== Start writing user histograms =====" << std::endl;
   hOutputFile->SetCompressionLevel(2);
   hOutputFile->cd();
   myTree->Write();
   hOutputFile->Close() ;
   std::cout << "===== End writing user histograms =======" << std::endl;
   */
   }

◆ endRun()

void cms::MinBias::endRun	(	edm::Run const &	,
		edm::EventSetup const &	iSetup
	)

override

Definition at line 52 of file MinBias.cc.

52 {}

Member Data Documentation

◆ allowMissingInputs_

bool cms::MinBias::allowMissingInputs_

private

Definition at line 44 of file MinBias.h.

Referenced by analyze(), and MinBias().

◆ caloGeometryESToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> cms::MinBias::caloGeometryESToken_

private

Definition at line 42 of file MinBias.h.

Referenced by beginRun().

◆ depth

int cms::MinBias::depth

private

Definition at line 49 of file MinBias.h.

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

◆ eta

float cms::MinBias::eta

private

Definition at line 50 of file MinBias.h.

Referenced by Particle.Particle::__str__(), beginJob(), datamodel.Object::DeltaR(), endJob(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

◆ geo_

const CaloGeometry* cms::MinBias::geo_

private

Definition at line 52 of file MinBias.h.

Referenced by beginRun(), and endJob().

◆ hbheLabel_

std::string cms::MinBias::hbheLabel_

private

Definition at line 38 of file MinBias.h.

Referenced by analyze(), and MinBias().

◆ hbheToken_

edm::EDGetTokenT<HBHERecHitCollection> cms::MinBias::hbheToken_

private

Definition at line 39 of file MinBias.h.

Referenced by analyze(), and MinBias().

◆ hfLabel_

std::string cms::MinBias::hfLabel_

private

Definition at line 38 of file MinBias.h.

Referenced by analyze(), and MinBias().

◆ hfToken_

edm::EDGetTokenT<HFRecHitCollection> cms::MinBias::hfToken_

private

Definition at line 41 of file MinBias.h.

Referenced by analyze(), and MinBias().

◆ hoLabel_

std::string cms::MinBias::hoLabel_

private

Definition at line 38 of file MinBias.h.

Referenced by analyze(), and MinBias().

◆ hoToken_

edm::EDGetTokenT<HORecHitCollection> cms::MinBias::hoToken_

private

Definition at line 40 of file MinBias.h.

Referenced by analyze(), and MinBias().

◆ ieta

int cms::MinBias::ieta

private

Definition at line 49 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ iphi

int cms::MinBias::iphi

private

Definition at line 49 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ mom1

float cms::MinBias::mom1

private

Definition at line 51 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ mom2

float cms::MinBias::mom2

private

Definition at line 51 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ mom3

float cms::MinBias::mom3

private

Definition at line 51 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ mom4

float cms::MinBias::mom4

private

Definition at line 51 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ mydet

int cms::MinBias::mydet

private

Definition at line 49 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ mysubd

int cms::MinBias::mysubd

private

Definition at line 49 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ myTree_

TTree* cms::MinBias::myTree_

private

Definition at line 47 of file MinBias.h.

Referenced by beginJob(), and endJob().

◆ occup

float cms::MinBias::occup

private

Definition at line 51 of file MinBias.h.

◆ phi

float cms::MinBias::phi

private

Definition at line 50 of file MinBias.h.

Referenced by Particle.Particle::__str__(), beginJob(), datamodel.Object::DeltaR(), endJob(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

◆ theFillDetMap0_

std::map<DetId, double> cms::MinBias::theFillDetMap0_

private

Definition at line 54 of file MinBias.h.

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

◆ theFillDetMap1_

std::map<DetId, double> cms::MinBias::theFillDetMap1_

private

Definition at line 55 of file MinBias.h.

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

◆ theFillDetMap2_

std::map<DetId, double> cms::MinBias::theFillDetMap2_

private

Definition at line 56 of file MinBias.h.

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

◆ theFillDetMap3_

std::map<DetId, double> cms::MinBias::theFillDetMap3_

private

Definition at line 57 of file MinBias.h.

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

◆ theFillDetMap4_

std::map<DetId, double> cms::MinBias::theFillDetMap4_

private

Definition at line 58 of file MinBias.h.

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MinBias()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ endJob()

◆ endRun()

Member Data Documentation

◆ allowMissingInputs_

◆ caloGeometryESToken_

◆ depth

◆ eta

◆ geo_

◆ hbheLabel_

◆ hbheToken_

◆ hfLabel_

◆ hfToken_

◆ hoLabel_

◆ hoToken_

◆ ieta

◆ iphi

◆ mom1

◆ mom2

◆ mom3

◆ mom4

◆ mydet

◆ mysubd

◆ myTree_

◆ occup

◆ phi

◆ theFillDetMap0_

◆ theFillDetMap1_

◆ theFillDetMap2_

◆ theFillDetMap3_

◆ theFillDetMap4_