Inheritance diagram for SimAnalyzerMinbias:

Classes
struct	myInfo

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

virtual void	beginJob ()

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

virtual void	endJob ()

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

	SimAnalyzerMinbias (const edm::ParameterSet &)

	~SimAnalyzerMinbias ()

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

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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::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 38 of file SimAnalyzerMinbias.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 72 of file SimAnalyzerMinbias.cc.

References edm::ParameterSet::getUntrackedParameter(), HLT_FULL_cff::InputTag, timeCut_, tok_evt_, and tok_hcal_.

                                                                      {
 
   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::LogInfo("AnalyzerMB")  << "Use Time cut of " << timeCut_ << " ns";
 }

SimAnalyzerMinbias::~SimAnalyzerMinbias ( )

Definition at line 83 of file SimAnalyzerMinbias.cc.

                                         {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 148 of file SimAnalyzerMinbias.cc.

References edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), info(), edm::HandleBase::isValid(), myMap_, edm::Handle< T >::product(), reco::return(), edm::EventID::run(), timeCut_, tok_evt_, and tok_hcal_.

                                        {
   edm::LogInfo("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::LogInfo("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::LogInfo("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::LogInfo("AnalyzerMB") << "ID " << hid << " with energy " << energyhit;
   }
 }

void SimAnalyzerMinbias::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 94 of file SimAnalyzerMinbias.cc.

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

                                   {
 
   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",  &mom2_MB, "mom3_MB/F");
   myTree_->Branch("mom4_MB",  &mom4_MB, "mom4_MB/F");
 
   myMap_.clear();
 }

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 88 of file SimAnalyzerMinbias.cc.

88 {

89 }

void SimAnalyzerMinbias::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 114 of file SimAnalyzerMinbias.cc.

References cells, depth, ieta, info(), iphi, mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, myMap_, mysubd, myTree_, SimAnalyzerMinbias::myInfo::theMB0, SimAnalyzerMinbias::myInfo::theMB1, SimAnalyzerMinbias::myInfo::theMB2, SimAnalyzerMinbias::myInfo::theMB3, and SimAnalyzerMinbias::myInfo::theMB4.

                                 {
    
   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::LogInfo("AnalyzerMB") << " Result=  " << mysubd << " " << ieta << " "
                  << iphi << " mom0  " << mom0_MB << " mom1 " 
                  << mom1_MB << " mom2 " << mom2_MB << " mom3 "
                  << mom3_MB << " mom4 " << mom4_MB;
       myTree_->Fill();
     }
   }
   edm::LogInfo("AnalyzerMB") << "cells " << cells;    
 }

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 91 of file SimAnalyzerMinbias.cc.

91 {

92 }

Member Data Documentation

int SimAnalyzerMinbias::cells

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

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().

edm::Service<TFileService> SimAnalyzerMinbias::fs_

private

Definition at line 52 of file SimAnalyzerMinbias.cc.

Referenced by beginJob().

int SimAnalyzerMinbias::ieta

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int SimAnalyzerMinbias::iphi

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom0_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom1_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom2_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom3_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by endJob().

float SimAnalyzerMinbias::mom4_MB

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int SimAnalyzerMinbias::mydet

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob().

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

private

Definition at line 65 of file SimAnalyzerMinbias.cc.

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

int SimAnalyzerMinbias::mysubd

private

Definition at line 57 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

TTree* SimAnalyzerMinbias::myTree_

private

Definition at line 54 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

double SimAnalyzerMinbias::timeCut_

private

Definition at line 53 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().

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

private

Definition at line 66 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().

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

private

Definition at line 67 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().

Classes

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation