#include <DQMHcalDiJetsAlCaReco.h>

Inheritance diagram for DQMHcalDiJetsAlCaReco:

Public Member Functions
	DQMHcalDiJetsAlCaReco (const edm::ParameterSet &)

	~DQMHcalDiJetsAlCaReco () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

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)

Private Attributes
bool	allowMissingInputs_

edm::EDGetTokenT< EcalRecHitCollection >	ec_

int	eventCounter_

std::string	fileName_
	Output file name if required. More...

std::string	folderName_
	DQM folder name. More...

edm::EDGetTokenT< HBHERecHitCollection >	hbhe_

edm::EDGetTokenT< HFRecHitCollection >	hf_

MonitorElement *	hiDistrEtThirdJet_

MonitorElement *	hiDistrProbeJetEnergy_

MonitorElement *	hiDistrProbeJetEta_

MonitorElement *	hiDistrProbeJetPhi_

MonitorElement *	hiDistrRecHitEnergyEBEE_

MonitorElement *	hiDistrRecHitEnergyHBHE_

MonitorElement *	hiDistrRecHitEnergyHF_

MonitorElement *	hiDistrRecHitEnergyHO_

MonitorElement *	hiDistrTagJetEnergy_

MonitorElement *	hiDistrTagJetEta_

MonitorElement *	hiDistrTagJetPhi_

edm::EDGetTokenT< HORecHitCollection >	ho_

edm::EDGetTokenT< reco::CaloJetCollection >	jets_
	object to monitor More...

bool	saveToFile_
	Write to file. More...

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 27 of file DQMHcalDiJetsAlCaReco.h.

Constructor & Destructor Documentation

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

Definition at line 42 of file DQMHcalDiJetsAlCaReco.cc.

References ec_, fileName_, folderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hbhe_, hf_, ho_, jets_, and saveToFile_.

                                                                              :
 eventCounter_(0)
 {
 //
 // Input from configurator file 
 //
   folderName_ = iConfig.getUntrackedParameter<string>("FolderName","ALCAStreamHcalDiJets");
 
 
   jets_= consumes<CaloJetCollection>(iConfig.getParameter<edm::InputTag>("jetsInput"));
   ec_= consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("ecInput"));
   hbhe_= consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
   ho_= consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"));
   hf_= consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInput"));
   
   saveToFile_= iConfig.getUntrackedParameter<bool>("SaveToFile",false);
   fileName_=  iConfig.getUntrackedParameter<string>("FileName","MonitorAlCaHcalDiJets.root");
     
 }

DQMHcalDiJetsAlCaReco::~DQMHcalDiJetsAlCaReco ( )

override

Definition at line 62 of file DQMHcalDiJetsAlCaReco.cc.

63 {}

Member Function Documentation

void DQMHcalDiJetsAlCaReco::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotected

Definition at line 123 of file DQMHcalDiJetsAlCaReco.cc.

References edm::SortedCollection< T, SORT >::begin(), RecoEcal_EventContent_cff::ec, ec_, edm::SortedCollection< T, SORT >::end(), reco::LeafCandidate::energy(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), eventCounter_, MonitorElement::Fill(), edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, hbhe_, photonIsolationHIProducer_cfi::hf, hf_, hiDistrEtThirdJet_, hiDistrProbeJetEnergy_, hiDistrProbeJetEta_, hiDistrProbeJetPhi_, hiDistrRecHitEnergyEBEE_, hiDistrRecHitEnergyHBHE_, hiDistrRecHitEnergyHF_, hiDistrRecHitEnergyHO_, hiDistrTagJetEnergy_, hiDistrTagJetEta_, hiDistrTagJetPhi_, photonIsolationHIProducer_cfi::ho, ho_, edm::HandleBase::isValid(), metsig::jet, fwrapper::jets, jets_, LogDebug, reco::LeafCandidate::phi(), and mathSSE::return().

                                              {  
  
 
    eventCounter_++;
 
    CaloJet jet1, jet2, jet3;
    Float_t etVetoJet; 
 
    edm::Handle<CaloJetCollection> jets;
    iEvent.getByToken(jets_,jets);
    
   if(!jets.isValid()){
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't getjet product!" << std::endl;
     return ;
   }
    
    
    
    if(jets->size()>1){
     jet1 = (*jets)[0];
     jet2 = (*jets)[1];
      if(fabs(jet1.eta())>fabs(jet2.eta())){
        CaloJet jet = jet1;
        jet1 = jet2;
        jet2 = jet;
      }
      //     if(fabs(jet1.eta())>eta_1 || (fabs(jet2.eta())-jet_R) < eta_2){ return;}
    } else {return;}
 
    hiDistrTagJetEnergy_->Fill(jet1.energy());
    hiDistrTagJetEta_->Fill(jet1.eta());
    hiDistrTagJetPhi_->Fill(jet1.phi());
 
    hiDistrProbeJetEnergy_->Fill(jet2.energy());
    hiDistrProbeJetEta_->Fill(jet2.eta());
    hiDistrProbeJetPhi_->Fill(jet2.phi());
 
    if(jets->size()>2){
      jet3 = (*jets)[2];
      etVetoJet = jet3.et();
      hiDistrEtThirdJet_->Fill(etVetoJet);
    } else { etVetoJet = 0.; hiDistrEtThirdJet_->Fill(etVetoJet); }
 
 
 
       Handle<EcalRecHitCollection> ec;
       iEvent.getByToken(ec_,ec);
       
   if(!ec.isValid()){
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get ec product!" << std::endl;
     return ;
   }
       
       
       for(EcalRecHitCollection::const_iterator ecItr = (*ec).begin();
                                                 ecItr != (*ec).end(); ++ecItr)
       {
         hiDistrRecHitEnergyEBEE_->Fill(ecItr->energy()); 
       }
 
 
 
 
       Handle<HBHERecHitCollection> hbhe;
       iEvent.getByToken(hbhe_, hbhe);
 
   if(!hbhe.isValid()){
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get hbhe product!" << std::endl;
     return ;
   }
 
 
       for(HBHERecHitCollection::const_iterator hbheItr=hbhe->begin();
                                                  hbheItr!=hbhe->end(); hbheItr++)
       {
     hiDistrRecHitEnergyHBHE_->Fill(hbheItr->energy()); 
       }
 
    
       Handle<HORecHitCollection> ho;
       iEvent.getByToken(ho_, ho);
 
   if(!ho.isValid()){
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get ho product!" << std::endl;
     return ;
   }
 
 
       for(HORecHitCollection::const_iterator hoItr=ho->begin();
                                                hoItr!=ho->end(); hoItr++)
       {
          hiDistrRecHitEnergyHO_->Fill(hoItr->energy());
 
       }
 
 
 
 
       Handle<HFRecHitCollection> hf;
       iEvent.getByToken(hf_, hf);
 
   if(!hf.isValid()){
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get hf product!" << std::endl;
     return ;
   }
 
 
       for(HFRecHitCollection::const_iterator hfItr=hf->begin();
                                                hfItr!=hf->end(); hfItr++)
       {
     hiDistrRecHitEnergyHF_->Fill(hfItr->energy()); 
       }
     
 } //analyze

void DQMHcalDiJetsAlCaReco::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 66 of file DQMHcalDiJetsAlCaReco.cc.

References DQMStore::IBooker::book1D(), folderName_, hiDistrEtThirdJet_, hiDistrProbeJetEnergy_, hiDistrProbeJetEta_, hiDistrProbeJetPhi_, hiDistrRecHitEnergyEBEE_, hiDistrRecHitEnergyHBHE_, hiDistrRecHitEnergyHF_, hiDistrRecHitEnergyHO_, hiDistrTagJetEnergy_, hiDistrTagJetEta_, hiDistrTagJetPhi_, MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

 {  
 
   // create and cd into new folder
   ibooker.setCurrentFolder(folderName_);
 
   // book some histograms 1D
   hiDistrRecHitEnergyEBEE_ = ibooker.book1D("RecHitEnergyEBEE", "the number of hits inside jets", 100, 0, 800);
   hiDistrRecHitEnergyEBEE_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyEBEE_->setAxisTitle("# rechits", 2);
 
   hiDistrRecHitEnergyHBHE_ = ibooker.book1D("RecHitEnergyHBHE", "the number of hits inside jets", 100,0, 800);
   hiDistrRecHitEnergyHBHE_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyHBHE_->setAxisTitle("# rechits", 2);
 
   hiDistrRecHitEnergyHF_ = ibooker.book1D("RecHitEnergyHF", "the number of hits inside jets", 150,0, 1500); 
   hiDistrRecHitEnergyHF_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyHF_->setAxisTitle("# rechits", 2);
 
   hiDistrRecHitEnergyHO_ = ibooker.book1D("RecHitEnergyHO", "the number of hits inside jets", 100,0, 100); 
   hiDistrRecHitEnergyHO_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyHO_->setAxisTitle("# rechits", 2);
 
   hiDistrProbeJetEnergy_ = ibooker.book1D("ProbeJetEnergy", "the energy of probe jets", 250,0, 2500); 
   hiDistrProbeJetEnergy_->setAxisTitle("E, GeV", 1);
   hiDistrProbeJetEnergy_->setAxisTitle("# jets", 2);
 
   hiDistrProbeJetEta_ = ibooker.book1D("ProbeJetEta", "the number of probe jets", 100, -5., 5.); 
   hiDistrProbeJetEta_->setAxisTitle("#eta", 1);
   hiDistrProbeJetEta_->setAxisTitle("# jets", 2);
 
   hiDistrProbeJetPhi_ = ibooker.book1D("ProbeJetPhi", "the number of probe jets", 50, -3.14, 3.14); 
   hiDistrProbeJetPhi_->setAxisTitle("#phi", 1);
   hiDistrProbeJetPhi_->setAxisTitle("# jets", 2);
 
   hiDistrTagJetEnergy_ = ibooker.book1D("TagJetEnergy", "the energy of tsg jets", 250,0, 2500); 
   hiDistrTagJetEnergy_->setAxisTitle("E, GeV", 1);
   hiDistrTagJetEnergy_->setAxisTitle("# jets", 2);
 
   hiDistrTagJetEta_ = ibooker.book1D("TagJetEta", "the number of  tag jets", 100, -5., 5.); 
   hiDistrTagJetEta_->setAxisTitle("#eta", 1);
   hiDistrTagJetEta_->setAxisTitle("# jets", 2);
 
   hiDistrTagJetPhi_ = ibooker.book1D("TagJetPhi", "the number of tag jets", 50, -3.14, 3.14); 
   hiDistrTagJetPhi_->setAxisTitle("#phi", 1);
   hiDistrTagJetPhi_->setAxisTitle("# jets", 2);
 
   hiDistrEtThirdJet_ = ibooker.book1D("EtThirdJet", "Et of the third jet", 90, 0, 90); 
 
 
 //==================================================================================
 
 
 }