#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

	MinBias (const edm::ParameterSet &)

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

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

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

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
bool	allowMissingInputs_

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::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 &)

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

Definition at line 35 of file MinBias.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 10 of file MinBias.cc.

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

                                                : 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

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

override

Definition at line 107 of file MinBias.cc.

References allowMissingInputs_, edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, hbheLabel_, hbheToken_, photonIsolationHIProducer_cfi::hf, hfLabel_, hfToken_, photonIsolationHIProducer_cfi::ho, hoLabel_, hoToken_, triggerObjects_cff::id, edm::HandleBase::isValid(), 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);
       }
     }
   }
   
 }

void cms::MinBias::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 22 of file MinBias.cc.

References depth, eta, ieta, iphi, TFileService::make(), 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");
 
 }

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

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 39 of file MinBias.cc.

References geo_, edm::EventSetup::get(), CaloGeometry::getValidDetIds(), DetId::Hcal, triggerObjects_cff::id, edm::ESHandle< T >::product(), theFillDetMap0_, theFillDetMap1_, theFillDetMap2_, theFillDetMap3_, and theFillDetMap4_.

                                                                       {
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   geo_ = pG.product();
   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.;
     }
   }
 }

void cms::MinBias::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 56 of file MinBias.cc.

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

                      {
   const HcalGeometry* hgeo = (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;
   */
 }

Member Data Documentation

bool cms::MinBias::allowMissingInputs_

private

Definition at line 53 of file MinBias.h.

Referenced by analyze(), and MinBias().

int cms::MinBias::depth

private

Definition at line 58 of file MinBias.h.

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

float cms::MinBias::eta

private

Definition at line 59 of file MinBias.h.

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

const CaloGeometry* cms::MinBias::geo_

private

Definition at line 61 of file MinBias.h.

Referenced by beginRun(), and endJob().

std::string cms::MinBias::hbheLabel_

private

Definition at line 48 of file MinBias.h.

Referenced by analyze(), and MinBias().

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

private

Definition at line 49 of file MinBias.h.

Referenced by analyze(), and MinBias().

std::string cms::MinBias::hfLabel_

private

Definition at line 48 of file MinBias.h.

Referenced by analyze(), and MinBias().

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

private

Definition at line 51 of file MinBias.h.

Referenced by analyze(), and MinBias().

std::string cms::MinBias::hoLabel_

private

Definition at line 48 of file MinBias.h.

Referenced by analyze(), and MinBias().

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

private

Definition at line 50 of file MinBias.h.

Referenced by analyze(), and MinBias().

int cms::MinBias::ieta

private

Definition at line 58 of file MinBias.h.

Referenced by beginJob(), and endJob().

int cms::MinBias::iphi

private

Definition at line 58 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom1

private

Definition at line 60 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom2

private

Definition at line 60 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom3

private

Definition at line 60 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom4

private

Definition at line 60 of file MinBias.h.

Referenced by beginJob(), and endJob().

int cms::MinBias::mydet

private

Definition at line 58 of file MinBias.h.

Referenced by beginJob(), and endJob().

int cms::MinBias::mysubd

private

Definition at line 58 of file MinBias.h.

Referenced by beginJob(), and endJob().

TTree* cms::MinBias::myTree_

private

Definition at line 56 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::occup

private

Definition at line 60 of file MinBias.h.

float cms::MinBias::phi

private

Definition at line 59 of file MinBias.h.

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

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

private

Definition at line 63 of file MinBias.h.

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

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

private

Definition at line 64 of file MinBias.h.

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

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

private

Definition at line 65 of file MinBias.h.

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

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

private

Definition at line 66 of file MinBias.h.

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

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

private

Definition at line 67 of file MinBias.h.

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation