SimAnalyzerMinbias Class Reference

Inheritance diagram for SimAnalyzerMinbias:

Classes
struct	myInfo

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	beginRun (const edm::Run &r, const edm::EventSetup &iSetup) override
void	endJob () override
void	endRun (const edm::Run &r, const edm::EventSetup &iSetup) override
	SimAnalyzerMinbias (const edm::ParameterSet &)
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
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 &&)=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::vector< ModuleDescription const * > &modules, 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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
int	cells
int	depth
edm::Service< TFileService >	fs_
int	ieta
int	iphi
float	mom0_MB
float	mom1_MB
float	mom2_MB
float	mom3_MB
float	mom4_MB
int	mydet
std::map< HcalDetId, myInfo >	myMap_
int	mysubd
TTree *	myTree_
double	timeCut_
edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_
edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hcal_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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 40 of file SimAnalyzerMinbias.cc.

Constructor & Destructor Documentation

SimAnalyzerMinbias()

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

explicit

Definition at line 70 of file SimAnalyzerMinbias.cc.

                                                                     {
  usesResource(TFileService::kSharedResource);
  timeCut_ = iConfig.getUntrackedParameter<double>("TimeCut", 500);
 
  // get token names of modules, producing object collections
  tok_evt_ = consumes<edm::HepMCProduct>(edm::InputTag("generator"));
  tok_hcal_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "HcalHits"));
 
  edm::LogVerbatim("AnalyzerMB") << "Use Time cut of " << timeCut_ << " ns";
}

References edm::ParameterSet::getUntrackedParameter(), TFileService::kSharedResource, timeCut_, tok_evt_, and tok_hcal_.

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 137 of file SimAnalyzerMinbias.cc.

                                                                             {
  edm::LogVerbatim("AnalyzerMB") << " Start SimAnalyzerMinbias::analyze " << iEvent.id().run() << ":"
                                 << iEvent.id().event();
 
  edm::Handle<edm::HepMCProduct> evtMC;
  iEvent.getByToken(tok_evt_, evtMC);
  if (!evtMC.isValid()) {
    edm::LogWarning("AnalyzerMB") << "no HepMCProduct found";
  } else {
    const HepMC::GenEvent* myGenEvent = evtMC->GetEvent();
    edm::LogVerbatim("AnalyzerMB") << "Event with " << myGenEvent->particles_size() << " particles + "
                                   << myGenEvent->vertices_size() << " vertices";
  }
 
  edm::Handle<edm::PCaloHitContainer> hcalHits;
  iEvent.getByToken(tok_hcal_, hcalHits);
  if (!hcalHits.isValid()) {
    edm::LogWarning("AnalyzerMB") << "Error! can't get HcalHits product!";
    return;
  }
 
  const edm::PCaloHitContainer* HitHcal = hcalHits.product();
  std::map<HcalDetId, double> hitMap;
  for (std::vector<PCaloHit>::const_iterator hcalItr = HitHcal->begin(); hcalItr != HitHcal->end(); ++hcalItr) {
    double time = hcalItr->time();
    if (time < timeCut_) {
      double energyhit = hcalItr->energy();
      HcalDetId hid = HcalDetId(hcalItr->id());
      std::map<HcalDetId, double>::iterator itr1 = hitMap.find(hid);
      if (itr1 == hitMap.end()) {
        hitMap[hid] = 0;
        itr1 = hitMap.find(hid);
      }
      itr1->second += energyhit;
    }
  }
  edm::LogVerbatim("AnalyzerMB") << "extract information of " << hitMap.size() << " towers from " << HitHcal->size()
                                 << " hits";
 
  for (std::map<HcalDetId, double>::const_iterator hcalItr = hitMap.begin(); hcalItr != hitMap.end(); ++hcalItr) {
    HcalDetId hid = hcalItr->first;
    double energyhit = hcalItr->second;
    std::map<HcalDetId, myInfo>::iterator itr1 = myMap_.find(hid);
    if (itr1 == myMap_.end()) {
      myInfo info;
      myMap_[hid] = info;
      itr1 = myMap_.find(hid);
    }
    itr1->second.theMB0++;
    itr1->second.theMB1 += energyhit;
    itr1->second.theMB2 += (energyhit * energyhit);
    itr1->second.theMB3 += (energyhit * energyhit * energyhit);
    itr1->second.theMB4 += (energyhit * energyhit * energyhit * energyhit);
    edm::LogVerbatim("AnalyzerMB") << "ID " << hid << " with energy " << energyhit;
  }
}

References edm::HepMCProduct::GetEvent(), iEvent, info(), edm::HandleBase::isValid(), myMap_, edm::Handle< T >::product(), ntuplemaker::time, timeCut_, tok_evt_, and tok_hcal_.

beginJob()

void SimAnalyzerMinbias::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 87 of file SimAnalyzerMinbias.cc.

                                  {
  myTree_ = fs_->make<TTree>("SimJet", "SimJet 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("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");
 
  myMap_.clear();
}

References cells, depth, fs_, ieta, iphi, TFileService::make(), mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, mydet, myMap_, mysubd, and myTree_.

beginRun()

void SimAnalyzerMinbias::beginRun ( const edm::Run &  r, const edm::EventSetup &  iSetup  )

inlineoverride

Definition at line 47 of file SimAnalyzerMinbias.cc.

{};

endJob()

void SimAnalyzerMinbias::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 106 of file SimAnalyzerMinbias.cc.

                                {
  cells = 0;
  for (std::map<HcalDetId, myInfo>::const_iterator itr = myMap_.begin(); itr != myMap_.end(); ++itr) {
    mysubd = itr->first.subdet();
    depth = itr->first.depth();
    iphi = itr->first.iphi();
    ieta = itr->first.ieta();
    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;
      cells++;
 
      edm::LogVerbatim("AnalyzerMB") << " Result=  " << mysubd << " " << ieta << " " << iphi << " mom0  " << mom0_MB
                                     << " mom1 " << mom1_MB << " mom2 " << mom2_MB << " mom3 " << mom3_MB << " mom4 "
                                     << mom4_MB;
      myTree_->Fill();
    }
  }
  edm::LogVerbatim("AnalyzerMB") << "cells " << cells;
}

References cells, depth, ieta, info(), iphi, mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, myMap_, mysubd, and myTree_.

Referenced by o2olib.O2ORunMgr::executeJob().

endRun()

void SimAnalyzerMinbias::endRun ( const edm::Run &  r, const edm::EventSetup &  iSetup  )

inlineoverride

Definition at line 48 of file SimAnalyzerMinbias.cc.

{};

fillDescriptions()

void SimAnalyzerMinbias::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 81 of file SimAnalyzerMinbias.cc.

                                                                                    {
  edm::ParameterSetDescription desc;
  desc.addUntracked<double>("TimeCut", 500);
  descriptions.add("simAnalyzerMinbias", desc);
}

References edm::ConfigurationDescriptions::add(), and submitPVResolutionJobs::desc.

Member Data Documentation

cells

int SimAnalyzerMinbias::cells

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

depth

int SimAnalyzerMinbias::depth

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

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

fs_

edm::Service<TFileService> SimAnalyzerMinbias::fs_

private

Definition at line 48 of file SimAnalyzerMinbias.cc.

Referenced by beginJob().

ieta

int SimAnalyzerMinbias::ieta

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

iphi

int SimAnalyzerMinbias::iphi

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom0_MB

float SimAnalyzerMinbias::mom0_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom1_MB

float SimAnalyzerMinbias::mom1_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom2_MB

float SimAnalyzerMinbias::mom2_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom3_MB

float SimAnalyzerMinbias::mom3_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom4_MB

float SimAnalyzerMinbias::mom4_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mydet

int SimAnalyzerMinbias::mydet

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob().

myMap_

std::map<HcalDetId, myInfo> SimAnalyzerMinbias::myMap_

private

Definition at line 63 of file SimAnalyzerMinbias.cc.

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

mysubd

int SimAnalyzerMinbias::mysubd

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

myTree_

TTree* SimAnalyzerMinbias::myTree_

private

Definition at line 54 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

timeCut_

double SimAnalyzerMinbias::timeCut_

private

Definition at line 53 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().

tok_evt_

edm::EDGetTokenT<edm::HepMCProduct> SimAnalyzerMinbias::tok_evt_

private

Definition at line 64 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().

tok_hcal_

edm::EDGetTokenT<edm::PCaloHitContainer> SimAnalyzerMinbias::tok_hcal_

private

Definition at line 65 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().