cms::MinBias Class Reference

#include <MinBias.h>

Inheritance diagram for cms::MinBias:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()
	MinBias (const edm::ParameterSet &)
	~MinBias ()
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
bool	allowMissingInputs_
int	depth
float	eta
std::string	fOutputFileName
const CaloGeometry *	geo
std::string	hbheLabel_
edm::EDGetTokenT < HBHERecHitCollection >	hbheToken_
std::string	hfLabel_
edm::EDGetTokenT < HFRecHitCollection >	hfToken_
std::string	hoLabel_
edm::EDGetTokenT < HORecHitCollection >	hoToken_
TFile *	hOutputFile
int	ieta
int	ievent
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 &)
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 42 of file MinBias.h.

Constructor & Destructor Documentation

cms::MinBias::MinBias ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 9 of file MinBias.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_25ns10e33_v2_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

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

}

cms::MinBias::~MinBias

(

)

Definition at line 25 of file MinBias.cc.

{

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

}

Member Function Documentation

void cms::MinBias::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 107 of file MinBias.cc.

References gather_cfg::cout, edm::EventSetup::get(), edm::Event::getByToken(), DetId::Hcal, edm::HandleBase::isValid(), funct::pow(), and edm::ESHandle< class >::product().

{

   using namespace edm;
   if(ievent == 0 ){
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   geo = pG.product();
   std::vector<DetId> did =  geo->getValidDetIds();

   for(std::vector<DetId>::const_iterator id=did.begin(); id != did.end(); id++)
   {
      if( (*id).det() == DetId::Hcal ) {
      theFillDetMap0[*id] = 0.;
      theFillDetMap1[*id] = 0.;
      theFillDetMap2[*id] = 0.;
      theFillDetMap3[*id] = 0.;
      theFillDetMap4[*id] = 0.;
   }
   }
   }



  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(HBHERecHitCollection::const_iterator hbheItr = (*hbhe).begin();
      hbheItr != (*hbhe).end(); ++hbheItr)
    {
      DetId id = (hbheItr)->detid();
      if( (*hbheItr).energy() > 0. ) std::cout<<" Energy = "<<(*hbheItr).energy()<<std::endl;
      theFillDetMap0[id] = theFillDetMap0[id]+ 1.;
      theFillDetMap1[id] = theFillDetMap1[id]+(*hbheItr).energy();
      theFillDetMap2[id] = theFillDetMap2[id]+pow((*hbheItr).energy(),2);
      theFillDetMap3[id] = theFillDetMap3[id]+pow((*hbheItr).energy(),3);
      theFillDetMap4[id] = 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(HORecHitCollection::const_iterator hoItr = (*ho).begin();
      hoItr != (*ho).end(); ++hoItr)
    {
      DetId id = (hoItr)->detid();
      theFillDetMap0[id] = theFillDetMap0[id]+ 1.;
      theFillDetMap1[id] = theFillDetMap1[id]+(*hoItr).energy();
      theFillDetMap2[id] = theFillDetMap2[id]+pow((*hoItr).energy(),2);
      theFillDetMap3[id] = theFillDetMap3[id]+pow((*hoItr).energy(),3);
      theFillDetMap4[id] = 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(HFRecHitCollection::const_iterator hfItr = (*hf).begin();
      hfItr != (*hf).end(); ++hfItr)
    {
      DetId id = (hfItr)->detid();
      theFillDetMap0[id] = theFillDetMap0[id]+ 1.;
      theFillDetMap1[id] = theFillDetMap1[id]+(*hfItr).energy();
      theFillDetMap2[id] = theFillDetMap2[id]+pow((*hfItr).energy(),2);
      theFillDetMap3[id] = theFillDetMap3[id]+pow((*hfItr).energy(),3);
      theFillDetMap4[id] = theFillDetMap4[id]+pow((*hfItr).energy(),4);
    }
    }
  }

}

void cms::MinBias::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 33 of file MinBias.cc.

References HLT_25ns10e33_v2_cff::depth, eta, and phi().

{
   hOutputFile   = new TFile( fOutputFileName.c_str(), "RECREATE" ) ;
   myTree = new 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");

   ievent = 0;

}

void cms::MinBias::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 53 of file MinBias.cc.

References gather_cfg::cout, HcalDetId::depth(), HLT_25ns10e33_v2_cff::depth, eta, PV3DBase< T, PVType, FrameType >::eta(), DetId::Hcal, i, HcalDetId::ieta(), HcalDetId::iphi(), PV3DBase< T, PVType, FrameType >::phi(), phi(), funct::pow(), and HcalDetId::subdet().

{

   const std::vector<DetId>& did =  geo->getSubdetectorGeometry( DetId::Hcal, 1 )->getValidDetIds() ;
   int i=0;
   for(std::vector<DetId>::const_iterator id=did.begin(); id != did.end(); id++)
   {
//      if( (*id).det() == DetId::Hcal ) {
      GlobalPoint pos = geo->getPosition(*id);
      mydet = ((*id).rawId()>>28)&0xF;
      mysubd = ((*id).rawId()>>25)&0x7;
      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.;
      }
      std::cout<<" Detector "<<(*id).rawId()<<" mydet "<<mydet<<" "<<mysubd<<" "<<depth<<" "<<
      HcalDetId(*id).subdet()<<" "<<ieta<<" "<<iphi<<" "<<pos.eta()<<" "<<pos.phi()<<std::endl;
      std::cout<<" Energy "<<mom1<<" "<<mom2<<std::endl;
      myTree->Fill();
      i++;
//      }
   }
   std::cout<<" The number of CaloDet records "<<did.size()<<std::endl;
   std::cout<<" 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 62 of file MinBias.h.

int cms::MinBias::depth

private

Definition at line 68 of file MinBias.h.

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

float cms::MinBias::eta

private

Definition at line 69 of file MinBias.h.

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

std::string cms::MinBias::fOutputFileName

private

Definition at line 61 of file MinBias.h.

const CaloGeometry* cms::MinBias::geo

private

Definition at line 71 of file MinBias.h.

std::string cms::MinBias::hbheLabel_

private

Definition at line 55 of file MinBias.h.

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

private

Definition at line 56 of file MinBias.h.

std::string cms::MinBias::hfLabel_

private

Definition at line 55 of file MinBias.h.

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

private

Definition at line 58 of file MinBias.h.

std::string cms::MinBias::hoLabel_

private

Definition at line 55 of file MinBias.h.

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

private

Definition at line 57 of file MinBias.h.

TFile* cms::MinBias::hOutputFile

private

Definition at line 64 of file MinBias.h.

int cms::MinBias::ieta

private

Definition at line 68 of file MinBias.h.

int cms::MinBias::ievent

private

Definition at line 68 of file MinBias.h.

int cms::MinBias::iphi

private

Definition at line 68 of file MinBias.h.

float cms::MinBias::mom1

private

Definition at line 70 of file MinBias.h.

float cms::MinBias::mom2

private

Definition at line 70 of file MinBias.h.

float cms::MinBias::mom3

private

Definition at line 70 of file MinBias.h.

float cms::MinBias::mom4

private

Definition at line 70 of file MinBias.h.

int cms::MinBias::mydet

private

Definition at line 68 of file MinBias.h.

int cms::MinBias::mysubd

private

Definition at line 68 of file MinBias.h.

TTree* cms::MinBias::myTree

private

Definition at line 66 of file MinBias.h.

float cms::MinBias::occup

private

Definition at line 70 of file MinBias.h.

float cms::MinBias::phi

private

Definition at line 69 of file MinBias.h.

Referenced by Particle.Particle::__str__().

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

private

Definition at line 73 of file MinBias.h.

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

private

Definition at line 74 of file MinBias.h.

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

private

Definition at line 75 of file MinBias.h.

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

private

Definition at line 76 of file MinBias.h.

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

private

Definition at line 77 of file MinBias.h.