cms::MinBias Class Reference

#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 ()
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
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 const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
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	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 &)
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
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 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<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 32 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_, HLT_FULL_cff::InputTag, 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  )

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 108 of file MinBias.cc.

References allowMissingInputs_, edm::Event::getByToken(), hbheLabel_, hbheToken_, hfLabel_, hfToken_, hcalSimParameters_cfi::ho, hoLabel_, hoToken_, gpuClustering::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 21 of file MinBias.cc.

References depth, eta, compareTotals::fs, 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 37 of file MinBias.cc.

References geo_, edm::EventSetup::get(), CaloGeometry::getValidDetIds(), DetId::Hcal, gpuClustering::id, edm::ESHandle< class >::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 54 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, gpuClustering::id, ieta, HcalDetId::ieta(), 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 = 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;
  */
  }

Member Data Documentation

bool cms::MinBias::allowMissingInputs_

private

Definition at line 50 of file MinBias.h.

Referenced by analyze(), and MinBias().

int cms::MinBias::depth

private

Definition at line 55 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 56 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 58 of file MinBias.h.

Referenced by beginRun(), and endJob().

std::string cms::MinBias::hbheLabel_

private

Definition at line 45 of file MinBias.h.

Referenced by analyze(), and MinBias().

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

private

Definition at line 46 of file MinBias.h.

Referenced by analyze(), and MinBias().

std::string cms::MinBias::hfLabel_

private

Definition at line 45 of file MinBias.h.

Referenced by analyze(), and MinBias().

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

private

Definition at line 48 of file MinBias.h.

Referenced by analyze(), and MinBias().

std::string cms::MinBias::hoLabel_

private

Definition at line 45 of file MinBias.h.

Referenced by analyze(), and MinBias().

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

private

Definition at line 47 of file MinBias.h.

Referenced by analyze(), and MinBias().

int cms::MinBias::ieta

private

Definition at line 55 of file MinBias.h.

Referenced by beginJob(), and endJob().

int cms::MinBias::iphi

private

Definition at line 55 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom1

private

Definition at line 57 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom2

private

Definition at line 57 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom3

private

Definition at line 57 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::mom4

private

Definition at line 57 of file MinBias.h.

Referenced by beginJob(), and endJob().

int cms::MinBias::mydet

private

Definition at line 55 of file MinBias.h.

Referenced by beginJob(), and endJob().

int cms::MinBias::mysubd

private

Definition at line 55 of file MinBias.h.

Referenced by beginJob(), and endJob().

TTree* cms::MinBias::myTree_

private

Definition at line 53 of file MinBias.h.

Referenced by beginJob(), and endJob().

float cms::MinBias::occup

private

Definition at line 57 of file MinBias.h.

float cms::MinBias::phi

private

Definition at line 56 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 60 of file MinBias.h.

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

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

private

Definition at line 61 of file MinBias.h.

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

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

private

Definition at line 62 of file MinBias.h.

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

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

private

Definition at line 63 of file MinBias.h.

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

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

private

Definition at line 64 of file MinBias.h.

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