SimAnalyzerMinbias Class Reference

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 (const std::string &iProcessName, std::vector< const char * > &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
std::string	fOutputFileName
TFile *	hOutputFile
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 73 of file SimAnalyzerMinbias.cc.

References gather_cfg::cout, fOutputFileName, edm::ParameterSet::getUntrackedParameter(), HLT_25ns14e33_v1_cff::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, timeCut, tok_evt_, and tok_hcal_.

                                                                     {
  // get name of output file with histogramms
  fOutputFileName = iConfig.getUntrackedParameter<std::string>("HistOutFile", "simOutput.root");
  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"));

#ifdef debugLog
  std::cout << "Use Time cut of " << timeCut << " ns and store o/p in "
        << fOutputFileName << std::endl;
#endif
}

SimAnalyzerMinbias::~SimAnalyzerMinbias

(

)

Definition at line 88 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 171 of file SimAnalyzerMinbias.cc.

References gather_cfg::cout, edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), info(), edm::HandleBase::isValid(), myMap, edm::Handle< T >::product(), mathSSE::return(), edm::EventID::run(), cond::rpcobgas::time, timeCut, tok_evt_, and tok_hcal_.

                                       {
#ifdef debugLog    
  std::cout << " Start SimAnalyzerMinbias::analyze " << iEvent.id().run() 
        << ":" << iEvent.id().event() << std::endl;
#endif
  
  edm::Handle<edm::HepMCProduct> evtMC;
  iEvent.getByToken(tok_evt_, evtMC);  
  if (!evtMC.isValid()) {
    edm::LogInfo("AnalyzerMB") << "no HepMCProduct found";
#ifdef debugLog
    std::cout                  << "no HepMCProduct found" << std::endl;
#endif
  } else {
    const HepMC::GenEvent * myGenEvent = evtMC->GetEvent();
    edm::LogInfo("AnalyzerMB") << "Event with " << myGenEvent->particles_size()
                   << " particles + " << myGenEvent->vertices_size()
                   << " vertices";
#ifdef debugLog
    std::cout                  << "Event with " << myGenEvent->particles_size()
                   << " particles + " << myGenEvent->vertices_size()
                   << " vertices" << std::endl;
#endif
  }

 
  edm::Handle<edm::PCaloHitContainer> hcalHits;
  iEvent.getByToken(tok_hcal_,hcalHits);
  if (!hcalHits.isValid()) {
    edm::LogInfo("AnalyzerMB") << "Error! can't get HcalHits product!";
#ifdef debugLog
    std::cout                  << "Error! can't get HcalHits product!" 
                   << std::endl;
#endif
    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;
    }
  }
#ifdef debugLog
  std::cout << "extract information of " << hitMap.size() << " towers from "
        << HitHcal->size() << " hits" << std::endl;
#endif

  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);
#ifdef debugLog
    std::cout << "ID " << hid << " with energy " << energyhit << std::endl;
#endif
  }
}

void SimAnalyzerMinbias::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 99 of file SimAnalyzerMinbias.cc.

References cells, depth, fOutputFileName, hOutputFile, ieta, iphi, mom0_MB, mom1_MB, mom2_MB, mom4_MB, mydet, myMap, mysubd, myTree, and mathSSE::return().

                                  {
  hOutputFile   = new TFile( fOutputFileName.c_str(), "RECREATE" ) ;

  myTree = new 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();
  return ;
}

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 93 of file SimAnalyzerMinbias.cc.

                                                                               {
}

void SimAnalyzerMinbias::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 121 of file SimAnalyzerMinbias.cc.

References cells, gather_cfg::cout, depth, hOutputFile, ieta, info(), iphi, mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, myMap, mysubd, myTree, mathSSE::return(), SimAnalyzerMinbias::myInfo::theMB0, SimAnalyzerMinbias::myInfo::theMB1, SimAnalyzerMinbias::myInfo::theMB2, SimAnalyzerMinbias::myInfo::theMB3, and SimAnalyzerMinbias::myInfo::theMB4.

                                {
   
  cells = 0;
  // std::cout << " Hello me here" << std::endl;
  for (std::map<HcalDetId,myInfo>::const_iterator itr=myMap.begin(); 
       itr != myMap.end(); ++itr) {
    //  std::cout << " Hello me here" << std::endl;
    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;
#ifdef debugLog
      std::cout                  << " Result=  " << mysubd << " " << ieta << " "
                 << iphi << " mom0  " << mom0_MB << " mom1 " 
                 << mom1_MB << " mom2 " << mom2_MB << " mom3 "
                 << mom3_MB << " mom4 " << mom4_MB << std::endl;
#endif
      myTree->Fill();
    }
  }
  edm::LogInfo("AnalyzerMB") << "cells " << cells;    
#ifdef debugLog
  std::cout                  << "cells " << cells << std::endl;
#endif
  hOutputFile->cd();
  myTree->Write();
  hOutputFile->Write();   
  hOutputFile->Close() ;
  return ;
}

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 96 of file SimAnalyzerMinbias.cc.

                                                                             {
}

Member Data Documentation

int SimAnalyzerMinbias::cells

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int SimAnalyzerMinbias::depth

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

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

std::string SimAnalyzerMinbias::fOutputFileName

private

Definition at line 52 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and SimAnalyzerMinbias().

TFile* SimAnalyzerMinbias::hOutputFile

private

Definition at line 54 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int SimAnalyzerMinbias::ieta

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int SimAnalyzerMinbias::iphi

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom0_MB

private

Definition at line 59 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom1_MB

private

Definition at line 59 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom2_MB

private

Definition at line 59 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float SimAnalyzerMinbias::mom3_MB

private

Definition at line 59 of file SimAnalyzerMinbias.cc.

Referenced by endJob().

float SimAnalyzerMinbias::mom4_MB

private

Definition at line 59 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int SimAnalyzerMinbias::mydet

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob().

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

private

Definition at line 66 of file SimAnalyzerMinbias.cc.

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

int SimAnalyzerMinbias::mysubd

private

Definition at line 58 of file SimAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

TTree* SimAnalyzerMinbias::myTree

private

Definition at line 55 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 67 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().

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

private

Definition at line 68 of file SimAnalyzerMinbias.cc.

Referenced by analyze(), and SimAnalyzerMinbias().