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
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)
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	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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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

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, ieta, HcalDetId::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 = 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 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().