#include <EcalRecHitsValidation.h>

Inheritance diagram for EcalRecHitsValidation:

Public Member Functions
	EcalRecHitsValidation (const edm::ParameterSet &ps)
	Constructor. More...

	~EcalRecHitsValidation () override
	Destructor. More...

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

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

void	findBarrelMatrix (int nCellInEta, int nCellInPhi, int CentralEta, int CentralPhi, int CentralZ, MapType &themap)

void	findEndcapMatrix (int nCellInX, int nCellInY, int CentralX, int CentralY, int CentralZ, MapType &themap)

uint32_t	getUnitWithMaxEnergy (MapType &themap)

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 Types
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType

Private Attributes
std::vector< uint32_t >	crystalMatrix

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	EBHits_Token_

edm::InputTag	EBrechitCollection_

edm::EDGetTokenT< EBRecHitCollection >	EBrechitCollection_Token_

edm::InputTag	EBuncalibrechitCollection_

edm::EDGetTokenT< EBUncalibratedRecHitCollection >	EBuncalibrechitCollection_Token_

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	EEHits_Token_

edm::InputTag	EErechitCollection_

edm::EDGetTokenT< EERecHitCollection >	EErechitCollection_Token_

edm::InputTag	EEuncalibrechitCollection_

edm::EDGetTokenT< EEUncalibratedRecHitCollection >	EEuncalibrechitCollection_Token_

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	ESHits_Token_

edm::InputTag	ESrechitCollection_

edm::EDGetTokenT< ESRecHitCollection >	ESrechitCollection_Token_

std::string	HepMCLabel

edm::EDGetTokenT< edm::HepMCProduct >	HepMCLabel_Token_

std::string	hitsProducer_

MonitorElement *	meEB5x5RecHitSimHitvsSimHitFlag8_

MonitorElement *	meEBe5x5_

MonitorElement *	meEBe5x5OverGun_

MonitorElement *	meEBe5x5OverSimHits_

MonitorElement *	meEBRecHitFlags_

MonitorElement *	meEBRecHitLog10Energy5x5Contr_

MonitorElement *	meEBRecHitLog10Energy_

MonitorElement *	meEBRecHitLog10EnergyContr_

MonitorElement *	meEBRecHitSimHitFlag6_

MonitorElement *	meEBRecHitSimHitFlag7_

MonitorElement *	meEBRecHitSimHitRatio1011_

MonitorElement *	meEBRecHitSimHitRatio12_

MonitorElement *	meEBRecHitSimHitRatio13_

MonitorElement *	meEBRecHitSimHitRatio_

MonitorElement *	meEBRecHitSimHitRatioGt35_

MonitorElement *	meEBRecHitSimHitvsSimHitFlag5_6_

MonitorElement *	meEBRecHitsOccupancyFlag5_6_

MonitorElement *	meEBRecHitsOccupancyFlag8_9_

MonitorElement *	meEBUnRecHitSimHitRatio_

MonitorElement *	meEBUnRecHitSimHitRatioGt35_

MonitorElement *	meEEe5x5_

MonitorElement *	meEEe5x5OverGun_

MonitorElement *	meEEe5x5OverSimHits_

MonitorElement *	meEERecHitFlags_

MonitorElement *	meEERecHitLog10Energy5x5Contr_

MonitorElement *	meEERecHitLog10Energy_

MonitorElement *	meEERecHitLog10EnergyContr_

MonitorElement *	meEERecHitSimHitFlag6_

MonitorElement *	meEERecHitSimHitFlag7_

MonitorElement *	meEERecHitSimHitRatio1011_

MonitorElement *	meEERecHitSimHitRatio12_

MonitorElement *	meEERecHitSimHitRatio13_

MonitorElement *	meEERecHitSimHitRatio_

MonitorElement *	meEERecHitSimHitRatioGt35_

MonitorElement *	meEERecHitSimHitvsSimHitFlag5_6_

MonitorElement *	meEERecHitsOccupancyMinusFlag5_6_

MonitorElement *	meEERecHitsOccupancyMinusFlag8_9_

MonitorElement *	meEERecHitsOccupancyPlusFlag5_6_

MonitorElement *	meEERecHitsOccupancyPlusFlag8_9_

MonitorElement *	meEEUnRecHitSimHitRatio_

MonitorElement *	meEEUnRecHitSimHitRatioGt35_

MonitorElement *	meESRecHitLog10Energy_

MonitorElement *	meESRecHitLog10EnergyContr_

MonitorElement *	meESRecHitSimHitRatio_

MonitorElement *	meGunEnergy_

MonitorElement *	meGunEta_

MonitorElement *	meGunPhi_

std::string	outputFile_

bool	verbose_

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 44 of file EcalRecHitsValidation.h.

Member Typedef Documentation

typedef std::map<uint32_t,float,std::less<uint32_t> > EcalRecHitsValidation::MapType

private

Definition at line 46 of file EcalRecHitsValidation.h.

Constructor & Destructor Documentation

EcalRecHitsValidation::EcalRecHitsValidation ( const edm::ParameterSet & ps )

Constructor.

Definition at line 23 of file EcalRecHitsValidation.cc.

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

                                                                   {
 
   // ---------------------- 
   HepMCLabel                 = ps.getParameter<std::string>("moduleLabelMC"); 
   hitsProducer_              = ps.getParameter<std::string>("hitsProducer");
   EBrechitCollection_        = ps.getParameter<edm::InputTag>("EBrechitCollection");
   EErechitCollection_        = ps.getParameter<edm::InputTag>("EErechitCollection");
   ESrechitCollection_        = ps.getParameter<edm::InputTag>("ESrechitCollection");
   EBuncalibrechitCollection_ = ps.getParameter<edm::InputTag>("EBuncalibrechitCollection");
   EEuncalibrechitCollection_ = ps.getParameter<edm::InputTag>("EEuncalibrechitCollection");
   // fix for consumes
   HepMCLabel_Token_         = consumes<HepMCProduct>(ps.getParameter<std::string>("moduleLabelMC"));
   EBrechitCollection_Token_ = consumes<EBRecHitCollection>(ps.getParameter<edm::InputTag>("EBrechitCollection"));
   EErechitCollection_Token_ = consumes<EERecHitCollection>(ps.getParameter<edm::InputTag>("EErechitCollection"));
   ESrechitCollection_Token_ = consumes<ESRecHitCollection>(ps.getParameter<edm::InputTag>("ESrechitCollection"));
   EBuncalibrechitCollection_Token_ = consumes<EBUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("EBuncalibrechitCollection"));
   EEuncalibrechitCollection_Token_ = consumes<EEUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("EEuncalibrechitCollection"));  
   EBHits_Token_ = consumes<CrossingFrame<PCaloHit> >(edm::InputTag(std::string("mix"), ps.getParameter<std::string>("hitsProducer") + std::string("EcalHitsEB")));
   EEHits_Token_ = consumes<CrossingFrame<PCaloHit> >(edm::InputTag(std::string("mix"), ps.getParameter<std::string>("hitsProducer") + std::string("EcalHitsEE")));
   ESHits_Token_ = consumes<CrossingFrame<PCaloHit> >(edm::InputTag(std::string("mix"), ps.getParameter<std::string>("hitsProducer") + std::string("EcalHitsES")));
 
   // ---------------------- 
   // DQM ROOT output 
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
 
   if ( !outputFile_.empty() ) {
     LogInfo("OutputInfo") << " Ecal RecHits Task histograms will be saved to '" << outputFile_.c_str() << "'";
   } else {
     LogInfo("OutputInfo") << " Ecal RecHits Task histograms will NOT be saved";
   }
  
   // ---------------------- 
   // verbosity switch 
   verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
 
   // ----------------------   
   meGunEnergy_                 = nullptr;
   meGunEta_                    = nullptr;   
   meGunPhi_                    = nullptr;   
   meEBRecHitSimHitRatio_       = nullptr;
   meEERecHitSimHitRatio_       = nullptr;
   meESRecHitSimHitRatio_       = nullptr;
   meEBRecHitSimHitRatio1011_   = nullptr;
   meEERecHitSimHitRatio1011_   = nullptr;
   meEBRecHitSimHitRatio12_     = nullptr;
   meEERecHitSimHitRatio12_     = nullptr;
   meEBRecHitSimHitRatio13_     = nullptr;
   meEERecHitSimHitRatio13_     = nullptr;
   meEBRecHitSimHitRatioGt35_   = nullptr;
   meEERecHitSimHitRatioGt35_   = nullptr;
   meEBUnRecHitSimHitRatio_     = nullptr;
   meEEUnRecHitSimHitRatio_     = nullptr;
   meEBUnRecHitSimHitRatioGt35_ = nullptr;
   meEEUnRecHitSimHitRatioGt35_ = nullptr;
   meEBe5x5_                    = nullptr;
   meEBe5x5OverSimHits_         = nullptr;
   meEBe5x5OverGun_             = nullptr;
   meEEe5x5_                    = nullptr;
   meEEe5x5OverSimHits_         = nullptr;
   meEEe5x5OverGun_             = nullptr;
   
   meEBRecHitLog10Energy_         = nullptr;
   meEERecHitLog10Energy_         = nullptr;
   meESRecHitLog10Energy_         = nullptr;
   meEBRecHitLog10EnergyContr_    = nullptr;
   meEERecHitLog10EnergyContr_    = nullptr;
   meESRecHitLog10EnergyContr_    = nullptr;
   meEBRecHitLog10Energy5x5Contr_ = nullptr;
   meEERecHitLog10Energy5x5Contr_ = nullptr;
 
   meEBRecHitsOccupancyFlag5_6_      = nullptr;
   meEBRecHitsOccupancyFlag8_9_      = nullptr;
   meEERecHitsOccupancyPlusFlag5_6_  = nullptr;
   meEERecHitsOccupancyMinusFlag5_6_ = nullptr;
   meEERecHitsOccupancyPlusFlag8_9_  = nullptr;
   meEERecHitsOccupancyMinusFlag8_9_ = nullptr;
   
   meEBRecHitFlags_                   = nullptr;
   meEBRecHitSimHitvsSimHitFlag5_6_   = nullptr;
   meEBRecHitSimHitFlag6_             = nullptr;
   meEBRecHitSimHitFlag7_             = nullptr;
   meEB5x5RecHitSimHitvsSimHitFlag8_  = nullptr;
 
   meEERecHitFlags_                   = nullptr;
   meEERecHitSimHitvsSimHitFlag5_6_   = nullptr;
   meEERecHitSimHitFlag6_             = nullptr;
   meEERecHitSimHitFlag7_             = nullptr;
 
 }

EcalRecHitsValidation::~EcalRecHitsValidation ( )

override

Destructor.

Definition at line 113 of file EcalRecHitsValidation.cc.

                                              {   
 
 }

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 242 of file EcalRecHitsValidation.cc.

References funct::abs(), cosmicBasicClusters_cfi::barrelHits, edm::SortedCollection< T, SORT >::begin(), EcalTrigTowerConstituentsMap::constituentsOf(), EcalCondObjectContainer< T >::end(), edm::SortedCollection< T, SORT >::end(), cosmicBasicClusters_cfi::endcapHits, edm::EventID::event(), HcalObjRepresent::Fill(), EcalCondObjectContainer< T >::find(), edm::SortedCollection< T, SORT >::find(), RemoveAddSevLevel::flag, edm::EventSetup::get(), edm::Event::getByToken(), EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), edm::HepMCProduct::GetEvent(), EcalChannelStatusCode::getStatusCode(), mps_fire::i, edm::EventBase::id(), EBDetId::ieta(), EBDetId::ietaAbs(), createfilelist::int, EBDetId::iphi(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), EcalRecHit::kDead, EcalRecHit::kLeadingEdgeRecovered, EcalRecHit::kNeighboursRecovered, EcalRecHit::kSaturated, EcalRecHit::kTowerRecovered, cmsBatch::log, LogDebug, M_PI, AlCaHLTBitMon_ParallelJobs::p, gedPhotons_cfi::preshowerHits, edm::Handle< T >::product(), edm::ESHandle< T >::product(), DetId::rawId(), rpcPointValidation_cfi::recHit, edm::EventID::run(), SimDataFormats::CaloAnalysis::sc, MixCollection< T >::size(), funct::tan(), EBDetId::tower(), EBDetId::zside(), and EEDetId::zside().

                                                                       {
 
   //Temporary stuff
 
 
   
   LogInfo("EcalRecHitsTask, EventInfo: ") << " Run = " << e.id().run() << " Event = " << e.id().event();
 
   // ADC -> GeV Scale
   edm::ESHandle<EcalADCToGeVConstant> pAgc;
   c.get<EcalADCToGeVConstantRcd>().get(pAgc);
   const EcalADCToGeVConstant* agc = pAgc.product();
   const double barrelADCtoGeV_  = agc->getEBValue();
   const double endcapADCtoGeV_ = agc->getEEValue();
   
   Handle<HepMCProduct> MCEvt;                   
   bool skipMC = false;
   e.getByToken(HepMCLabel_Token_, MCEvt);  
   if (!MCEvt.isValid()) { skipMC = true; }
 
   edm::Handle<CrossingFrame<PCaloHit> > crossingFrame;
 
   bool skipBarrel = false;
   const EBUncalibratedRecHitCollection *EBUncalibRecHit =nullptr;
   Handle< EBUncalibratedRecHitCollection > EcalUncalibRecHitEB;
   e.getByToken( EBuncalibrechitCollection_Token_, EcalUncalibRecHitEB);
   if (EcalUncalibRecHitEB.isValid()) {
     EBUncalibRecHit = EcalUncalibRecHitEB.product() ;    
   } else {
     skipBarrel = true;
   }
 
   bool skipEndcap = false;
   const EEUncalibratedRecHitCollection *EEUncalibRecHit = nullptr;
   Handle< EEUncalibratedRecHitCollection > EcalUncalibRecHitEE;
   e.getByToken( EEuncalibrechitCollection_Token_, EcalUncalibRecHitEE);
   if (EcalUncalibRecHitEE.isValid()){ 
     EEUncalibRecHit = EcalUncalibRecHitEE.product () ;
   } else {
     skipEndcap = true;
   }
 
   const EBRecHitCollection *EBRecHit = nullptr;
   Handle<EBRecHitCollection> EcalRecHitEB;
   e.getByToken( EBrechitCollection_Token_, EcalRecHitEB);
   if (EcalRecHitEB.isValid()){ 
     EBRecHit = EcalRecHitEB.product();
   } else {
     skipBarrel = true;
   }
 
   const EERecHitCollection *EERecHit = nullptr;
   Handle<EERecHitCollection> EcalRecHitEE;
   e.getByToken( EErechitCollection_Token_, EcalRecHitEE);
   if (EcalRecHitEE.isValid()){
     EERecHit = EcalRecHitEE.product ();
   } else {
     skipEndcap = true;
   }
 
   bool skipPreshower = false;
   const ESRecHitCollection *ESRecHit = nullptr;
   Handle<ESRecHitCollection> EcalRecHitES;
   e.getByToken( ESrechitCollection_Token_, EcalRecHitES);
   if (EcalRecHitES.isValid()) {
     ESRecHit = EcalRecHitES.product ();      
   } else {
     skipPreshower = true;
   }
 
 
   // ---------------------- 
   // gun
   double eGun = 0.;
   if ( ! skipMC ) {
     for ( HepMC::GenEvent::particle_const_iterator p = MCEvt->GetEvent()->particles_begin(); p != MCEvt->GetEvent()->particles_end(); ++p )  {      
       double htheta = (*p)->momentum().theta();
       double heta = -99999.;
       if( tan(htheta * 0.5) > 0 ) {
     heta = -log(tan(htheta * 0.5));
       }
       double hphi = (*p)->momentum().phi();
       hphi = (hphi>=0) ? hphi : hphi+2*M_PI;
       hphi = hphi / M_PI * 180.;
       
       LogDebug("EventInfo") << "EcalRecHitsTask: Particle gun type form MC = " << abs((*p)->pdg_id()) 
                 << "\n" << "Energy = "<< (*p)->momentum().e() 
                 << "\n" << "Eta = "   << heta 
                 << "\n" << "Phi = "   << hphi;
       
       if ( (*p)->momentum().e() > eGun ) eGun = (*p)->momentum().e();
 
       if (meGunEnergy_) meGunEnergy_->Fill((*p)->momentum().e());
       if (meGunEta_)    meGunEta_   ->Fill(heta);
       if (meGunPhi_)    meGunPhi_   ->Fill(hphi); 
     }
   }
 
   // -------------------------------------------------------------------
   // BARREL
 
   if ( ! skipBarrel) {
 
     // 1) loop over simHits  
     e.getByToken(EBHits_Token_,crossingFrame);
     const MixCollection<PCaloHit> barrelHits(crossingFrame.product());
     
     MapType ebSimMap;
     MapType ebRecMap;
     const int ebcSize = 90;
     double ebcontr[ebcSize]; 
     double ebcontr25[ebcSize];
     for( int i=0; i<ebcSize; i++ ) { ebcontr[i] = 0.0; ebcontr25[i] = 0.0; } 
     double ebtotal = 0.;
 
     for ( auto const & iHit : barrelHits ) {
       EBDetId ebid = EBDetId(iHit.id());
       
       LogDebug("SimHitInfo, barrel") 
     << "CaloHit "   << iHit.getName() << " DetID = " << iHit.id()   << "\n" 
     << "Energy = "  << iHit.energy()  << " Time = "  << iHit.time() << "\n"
     << "EBDetId = " << ebid.ieta()       << " "         << ebid.iphi();
       
       uint32_t crystid = ebid.rawId();
       ebSimMap[crystid] += iHit.energy();
     }
     
     
     
     // 2) loop over RecHits 
     for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EBUncalibRecHit->begin(); uncalibRecHit != EBUncalibRecHit->end() ; ++uncalibRecHit) {
       EBDetId EBid = EBDetId(uncalibRecHit->id());
       
       // Find corresponding recHit
       EcalRecHitCollection::const_iterator myRecHit = EBRecHit->find(EBid);
       if( myRecHit == EBRecHit->end() ) continue; 
       ebRecMap[EBid.rawId()] += myRecHit->energy();
       
       // Fill log10(Energy) stuff...   
       ebtotal += myRecHit->energy();
       if( myRecHit->energy() > 0 ) {
     if( meEBRecHitLog10Energy_ ) meEBRecHitLog10Energy_->Fill( log10( myRecHit->energy() ) );
     int log10i = int( ( log10( myRecHit->energy() ) + 5. ) * 10. );
     if( log10i >=0 and log10i < ebcSize ) ebcontr[ log10i ] += myRecHit->energy();
       }
 
       // comparison Rec/Sim hit
       if ( ebSimMap[EBid.rawId()] != 0. ) {
     double uncEnergy = uncalibRecHit->amplitude()*barrelADCtoGeV_;
     if (meEBUnRecHitSimHitRatio_)                                {meEBUnRecHitSimHitRatio_    ->Fill(uncEnergy/ebSimMap[EBid.rawId()]);}
     if (meEBUnRecHitSimHitRatioGt35_ && (myRecHit->energy()>3.5)){meEBUnRecHitSimHitRatioGt35_->Fill(uncEnergy/ebSimMap[EBid.rawId()]);}
       }
       
       if (myRecHit != EBRecHit->end()) {
     if ( ebSimMap[EBid.rawId()] != 0. ) {
       if (meEBRecHitSimHitRatio_)                                {meEBRecHitSimHitRatio_    ->Fill(myRecHit->energy()/ebSimMap[EBid.rawId()]);}
       if (meEBRecHitSimHitRatioGt35_ && (myRecHit->energy()>3.5)){meEBRecHitSimHitRatioGt35_->Fill(myRecHit->energy()/ebSimMap[EBid.rawId()]);}
       uint16_t sc = 0;
       edm::ESHandle<EcalChannelStatus> pEcs;
       c.get<EcalChannelStatusRcd>().get(pEcs); 
       const EcalChannelStatus* ecs = nullptr;
       if( pEcs.isValid() ) ecs = pEcs.product();
       if( ecs != nullptr ) {
         EcalChannelStatusMap::const_iterator csmi = ecs->find(EBid.rawId());
         EcalChannelStatusCode csc = 0;
         if( csmi != ecs->end() ) csc = *csmi;
         sc = csc.getStatusCode();
       }
 
       if( meEBRecHitSimHitRatio1011_ != nullptr && 
           ( sc == 10 || sc == 11 ) ) { meEBRecHitSimHitRatio1011_->Fill(myRecHit->energy()/ebSimMap[EBid.rawId()]); }
       if( meEBRecHitSimHitRatio12_ != nullptr && sc == 12 ) { meEBRecHitSimHitRatio12_->Fill(myRecHit->energy()/ebSimMap[EBid.rawId()]); }
 
       edm::ESHandle<EcalTrigTowerConstituentsMap> pttMap;
       c.get<IdealGeometryRecord>().get(pttMap);
       const EcalTrigTowerConstituentsMap* ttMap = nullptr;
       if( pttMap.isValid() ) ttMap = pttMap.product();
       double ttSimEnergy = 0;
       if( ttMap != nullptr ) {
         EcalTrigTowerDetId ttDetId = EBid.tower();
         std::vector<DetId> vid = ttMap->constituentsOf( ttDetId );
         for( std::vector<DetId>::const_iterator dit = vid.begin(); dit != vid.end(); dit++ ) {
           EBDetId ttEBid = EBDetId(*dit);
           ttSimEnergy += ebSimMap[ttEBid.rawId()];
         }
         if( !vid.empty() ) ttSimEnergy = ttSimEnergy / vid.size();
       }
       if( meEBRecHitSimHitRatio13_ != nullptr && sc == 13 && ttSimEnergy != 0 ) 
         meEBRecHitSimHitRatio13_->Fill(myRecHit->energy()/ttSimEnergy); 
 
       int flag = myRecHit->recoFlag();
       if( meEBRecHitFlags_ != nullptr ) meEBRecHitFlags_->Fill( flag );
       if( meEBRecHitSimHitvsSimHitFlag5_6_  && ( flag == EcalRecHit::kSaturated || flag == EcalRecHit::kLeadingEdgeRecovered ))
         meEBRecHitSimHitvsSimHitFlag5_6_->Fill( myRecHit->energy()/ebSimMap[EBid.rawId()], ebSimMap[EBid.rawId()] );
       if( meEBRecHitSimHitFlag6_  && ( flag == EcalRecHit::kLeadingEdgeRecovered ))
         meEBRecHitSimHitFlag6_->Fill( myRecHit->energy()/ebSimMap[EBid.rawId()] );
       if( meEBRecHitSimHitFlag7_  && ( flag == EcalRecHit::kNeighboursRecovered ))
         meEBRecHitSimHitFlag6_->Fill( myRecHit->energy()/ebSimMap[EBid.rawId()] );
       if( meEB5x5RecHitSimHitvsSimHitFlag8_  && ( flag == EcalRecHit::kTowerRecovered ) && ttSimEnergy != 0 )
         meEB5x5RecHitSimHitvsSimHitFlag8_->Fill( myRecHit->energy()/ttSimEnergy, ttSimEnergy );
 
       if (meEBRecHitsOccupancyFlag5_6_ && ( (flag==EcalRecHit::kSaturated) || (flag==EcalRecHit::kLeadingEdgeRecovered) ) ) 
         meEBRecHitsOccupancyFlag5_6_  -> Fill(EBid.ieta(), EBid.iphi());
       if (meEBRecHitsOccupancyFlag8_9_ && ( (flag==EcalRecHit::kTowerRecovered) || (flag==EcalRecHit::kDead) ) ) 
         meEBRecHitsOccupancyFlag8_9_  -> Fill(EBid.ieta(), EBid.iphi());
 
     }
       }
       else
     continue;
     }  // loop over the UncalibratedRecHitCollection
     
     // RecHits matrix
     uint32_t  ebcenterid = getUnitWithMaxEnergy(ebRecMap);
     EBDetId myEBid(ebcenterid);
     int bx = myEBid.ietaAbs();
     int by = myEBid.iphi();
     int bz = myEBid.zside();
     findBarrelMatrix(5,5,bx,by,bz,ebRecMap);
     double e5x5rec = 0.;
     double e5x5sim = 0.;
     for ( unsigned int i = 0; i < crystalMatrix.size(); i++ ) {
       e5x5rec += ebRecMap[crystalMatrix[i]];
       e5x5sim += ebSimMap[crystalMatrix[i]];
       if( ebRecMap[crystalMatrix[i]] > 0 ) {
     int log10i25 = int( ( log10( ebRecMap[crystalMatrix[i]] ) + 5. ) * 10. );
     if( log10i25 >=0 && log10i25 < ebcSize ) ebcontr25[ log10i25 ] += ebRecMap[crystalMatrix[i]];
       }
     }
     
     if( meEBe5x5_ ) meEBe5x5_->Fill(e5x5rec);
     if ( e5x5sim > 0. && meEBe5x5OverSimHits_ ) meEBe5x5OverSimHits_->Fill(e5x5rec/e5x5sim);
     if ( eGun > 0. && meEBe5x5OverGun_ ) meEBe5x5OverGun_->Fill(e5x5rec/eGun);
     
     if( meEBRecHitLog10EnergyContr_  && ebtotal != 0 ) {
       for( int i=0; i<ebcSize; i++ ) {
     meEBRecHitLog10EnergyContr_->Fill( -5.+(float(i)+0.5)/10., ebcontr[i]/ebtotal );
       }
     }
     
     if( meEBRecHitLog10Energy5x5Contr_  && e5x5rec != 0 ) {
       for( int i=0; i<ebcSize; i++ ) {
     meEBRecHitLog10Energy5x5Contr_->Fill( -5.+(float(i)+0.5)/10., ebcontr25[i]/e5x5rec );
       }
     }
   }
 
   // -------------------------------------------------------------------
   // ENDCAP
 
   if ( ! skipEndcap ) {
 
     // 1) loop over simHits
     e.getByToken(EEHits_Token_,crossingFrame);
     const MixCollection<PCaloHit> endcapHits(crossingFrame.product());
   
     MapType eeSimMap;
     MapType eeRecMap;
     const int eecSize = 90;
     double eecontr[eecSize];
     double eecontr25[eecSize];
     for( int i=0; i<eecSize; i++ ) { eecontr[i] = 0.0; eecontr25[i] = 0.0; } 
     double eetotal = 0.;
  
     for ( auto const & iHit : endcapHits ) {
       EEDetId eeid(iHit.id());
 
       LogDebug("Endcap, HitInfo")
     <<" CaloHit "      << iHit.getName() << " DetID = "        << iHit.id()   << "\n"
     << "Energy = "     << iHit.energy()  << " Time = "         << iHit.time() << "\n"
     << "EEDetId side " << eeid.zside()      << " = " << eeid.ix() << " " << eeid.iy();
       
       uint32_t crystid = eeid.rawId();
       eeSimMap[crystid] += iHit.energy();
     }
 
 
 
     // 2) loop over RecHits
     for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EEUncalibRecHit->begin(); uncalibRecHit != EEUncalibRecHit->end(); ++uncalibRecHit) {
       EEDetId EEid = EEDetId(uncalibRecHit->id());
       
       // Find corresponding recHit
       EcalRecHitCollection::const_iterator myRecHit = EERecHit->find(EEid);
       if( myRecHit == EERecHit->end() ) continue; 
       eeRecMap[EEid.rawId()] += myRecHit->energy();
 
       // Fill log10(Energy) stuff...
       eetotal += myRecHit->energy();   
       if( myRecHit->energy() > 0 ) {
     if( meEERecHitLog10Energy_ ) meEERecHitLog10Energy_->Fill( log10( myRecHit->energy() ) );
     int log10i = int( ( log10( myRecHit->energy() ) + 5. ) * 10. );
     if( log10i >=0 and log10i < eecSize ) eecontr[ log10i ] += myRecHit->energy();
       }
 
       // comparison Rec/Sim hit
       if ( eeSimMap[EEid.rawId()] != 0. ) {
     double uncEnergy = uncalibRecHit->amplitude()*endcapADCtoGeV_;
     if (meEEUnRecHitSimHitRatio_)                                {meEEUnRecHitSimHitRatio_    ->Fill(uncEnergy/eeSimMap[EEid.rawId()]);}
     if (meEEUnRecHitSimHitRatioGt35_ && (myRecHit->energy()>3.5)){meEEUnRecHitSimHitRatioGt35_->Fill(uncEnergy/eeSimMap[EEid.rawId()]);}
       }
 
       if (myRecHit != EERecHit->end()) {
     if ( eeSimMap[EEid.rawId()] != 0. ) {
       if (meEERecHitSimHitRatio_)                                {meEERecHitSimHitRatio_    ->Fill(myRecHit->energy()/eeSimMap[EEid.rawId()]); }
       if (meEERecHitSimHitRatioGt35_ && (myRecHit->energy()>3.5)){meEERecHitSimHitRatioGt35_->Fill(myRecHit->energy()/eeSimMap[EEid.rawId()]); }
 
       edm::ESHandle<EcalChannelStatus> pEcs;
       c.get<EcalChannelStatusRcd>().get(pEcs);
       const EcalChannelStatus* ecs = nullptr;
       if( pEcs.isValid() ) ecs = pEcs.product();
       if( ecs != nullptr ) {
         EcalChannelStatusMap::const_iterator csmi = ecs->find(EEid.rawId());
         EcalChannelStatusCode csc = 0;
         if( csmi != ecs->end() ) csc = *csmi;
         uint16_t sc = csc.getStatusCode();
         if( meEERecHitSimHitRatio1011_ != nullptr && 
         ( sc == 10 || sc == 11 ) ) { meEERecHitSimHitRatio1011_->Fill(myRecHit->energy()/eeSimMap[EEid.rawId()]); }
         if( meEERecHitSimHitRatio12_ != nullptr && sc == 12 ) { meEERecHitSimHitRatio12_->Fill(myRecHit->energy()/eeSimMap[EEid.rawId()]); }
         if( meEERecHitSimHitRatio13_ != nullptr && sc == 13 ) { meEERecHitSimHitRatio13_->Fill(myRecHit->energy()/eeSimMap[EEid.rawId()]); }
       }
 
       int flag = myRecHit->recoFlag();
       if( meEERecHitFlags_ != nullptr ) meEERecHitFlags_->Fill( flag );
       if( meEERecHitSimHitvsSimHitFlag5_6_  && ( flag == EcalRecHit::kSaturated || flag == EcalRecHit::kLeadingEdgeRecovered ))
         meEERecHitSimHitvsSimHitFlag5_6_->Fill( myRecHit->energy()/eeSimMap[EEid.rawId()], eeSimMap[EEid.rawId()] );
       if( meEERecHitSimHitFlag6_  && ( flag == EcalRecHit::kLeadingEdgeRecovered ))
         meEERecHitSimHitFlag6_->Fill( myRecHit->energy()/eeSimMap[EEid.rawId()] );
       if( meEERecHitSimHitFlag7_  && ( flag == EcalRecHit::kNeighboursRecovered ))
         meEERecHitSimHitFlag6_->Fill( myRecHit->energy()/eeSimMap[EEid.rawId()] );
 
       if (EEid.zside() > 0) { 
         if (meEERecHitsOccupancyPlusFlag5_6_ && (( flag == EcalRecHit::kSaturated ) || ( flag == EcalRecHit::kLeadingEdgeRecovered ) ))  
           meEERecHitsOccupancyPlusFlag5_6_  ->Fill(EEid.ix(), EEid.iy()); 
         if (meEERecHitsOccupancyPlusFlag8_9_ && (( flag == EcalRecHit::kTowerRecovered ) || ( flag == EcalRecHit::kDead ) ))  
           meEERecHitsOccupancyPlusFlag8_9_  ->Fill(EEid.ix(), EEid.iy()); 
       }
       if (EEid.zside() < 0) { 
         if (meEERecHitsOccupancyMinusFlag5_6_ && (( flag == EcalRecHit::kSaturated ) || ( flag == EcalRecHit::kLeadingEdgeRecovered ) ))  
           meEERecHitsOccupancyMinusFlag5_6_  ->Fill(EEid.ix(), EEid.iy()); 
         if (meEERecHitsOccupancyMinusFlag8_9_ && (( flag == EcalRecHit::kTowerRecovered ) || ( flag == EcalRecHit::kDead ) ))  
           meEERecHitsOccupancyMinusFlag8_9_  ->Fill(EEid.ix(), EEid.iy()); 
       }
     }
       }
       else
     continue;
     }  // loop over the UncalibratedechitCollection
   
     // RecHits matrix
     uint32_t  eecenterid = getUnitWithMaxEnergy(eeRecMap);
     EEDetId myEEid(eecenterid);
     int bx = myEEid.ix();
     int by = myEEid.iy();
     int bz = myEEid.zside();
     findEndcapMatrix(5,5,bx,by,bz,eeRecMap);
     double e5x5rec = 0.;
     double e5x5sim = 0.;
     for ( unsigned int i = 0; i < crystalMatrix.size(); i++ ) {
       e5x5rec += eeRecMap[crystalMatrix[i]];
       e5x5sim += eeSimMap[crystalMatrix[i]];
       if( eeRecMap[crystalMatrix[i]] > 0 ) {
     int log10i25 = int( ( log10( eeRecMap[crystalMatrix[i]] ) + 5. ) * 10. );
     if( log10i25 >=0 && log10i25 < eecSize ) eecontr25[ log10i25 ] += eeRecMap[crystalMatrix[i]];
       }
     }
 
     if( meEEe5x5_ ) meEEe5x5_->Fill(e5x5rec);
     if ( e5x5sim > 0. && meEEe5x5OverSimHits_ ) meEEe5x5OverSimHits_->Fill(e5x5rec/e5x5sim);
     if ( eGun > 0. && meEEe5x5OverGun_ ) meEEe5x5OverGun_->Fill(e5x5rec/eGun);
     
 
     if( meEERecHitLog10EnergyContr_  && eetotal != 0 ) {
       for( int i=0; i<eecSize; i++ ) {
     meEERecHitLog10EnergyContr_->Fill( -5.+(float(i)+0.5)/10., eecontr[i]/eetotal );
       }
     }
     
     if( meEERecHitLog10Energy5x5Contr_  && e5x5rec != 0 ) {
       for( int i=0; i<eecSize; i++ ) {
     meEERecHitLog10Energy5x5Contr_->Fill( -5.+(float(i)+0.5)/10., eecontr25[i]/e5x5rec );
       }
     }
   }
 
   // -------------------------------------------------------------------
   // PRESHOWER
 
   if ( ! skipPreshower ) {
 
     // 1) loop over simHits
     e.getByToken(ESHits_Token_,crossingFrame);
     const MixCollection<PCaloHit> preshowerHits(crossingFrame.product());
 
     MapType esSimMap;
     const int escSize = 90;
     double escontr[escSize];
     for( int i=0; i<escSize; i++ ) { escontr[i] = 0.0; }
     double estotal = 0.;
 
     for ( auto const & iHit : preshowerHits ) {
       ESDetId esid(iHit.id());
       LogDebug("Preshower, HitInfo")
     <<" CaloHit "       << iHit.getName() << " DetID = "         << iHit.id()   << "\n"
     << "Energy = "      << iHit.energy()  << " Time = "          << iHit.time() << "\n"
     << "ESDetId strip " << esid.strip()      << " = " << esid.six() << " " << esid.siy();
       
       uint32_t crystid = esid.rawId();
       esSimMap[crystid] += iHit.energy();
     }
 
 
     // 2) loop over RecHits
     for (EcalRecHitCollection::const_iterator recHit = ESRecHit->begin(); recHit != ESRecHit->end(); ++recHit) {
       ESDetId ESid = ESDetId(recHit->id());
       if ( esSimMap[ESid.rawId()] != 0. ) { 
     
     // Fill log10(Energy) stuff...
     estotal += recHit->energy();   
     if( recHit->energy() > 0 ) {
       if( meESRecHitLog10Energy_ ) meESRecHitLog10Energy_->Fill( log10( recHit->energy() ) );
       int log10i = int( ( log10( recHit->energy() ) + 5. ) * 10. );
       if( log10i >=0 and log10i < escSize ) escontr[ log10i ] += recHit->energy();
     }
 
     if (meESRecHitSimHitRatio_) { 
       meESRecHitSimHitRatio_ ->Fill(recHit->energy()/esSimMap[ESid.rawId()]); 
     }
       }
       else
     continue;
     }  // loop over the RechitCollection
 
     if( meESRecHitLog10EnergyContr_  && estotal != 0 ) {
       for( int i=0; i<escSize; i++ ) {
     meESRecHitLog10EnergyContr_->Fill( -5.+(float(i)+0.5)/10., escontr[i]/estotal );
       }
     }
 
     
   }
   
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 117 of file EcalRecHitsValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), trackerHits::histo, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                        {
 
   std::string histo;
    
   ibooker.setCurrentFolder("EcalRecHitsV/EcalRecHitsTask");
   
   histo = "EcalRecHitsTask Gun Momentum";
   meGunEnergy_ = ibooker.book1D(histo.c_str(), histo.c_str(), 100, 0., 1000.);
 
   histo = "EcalRecHitsTask Gun Eta";      
   meGunEta_ = ibooker.book1D(histo.c_str(), histo.c_str(), 700, -3.5, 3.5);
    
   histo = "EcalRecHitsTask Gun Phi";  
   meGunPhi_ = ibooker.book1D(histo.c_str(), histo.c_str(), 360, 0., 360.);    
   
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio";  
   meEBRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);   
 
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio"; 
   meEERecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Preshower RecSimHit Ratio"; 
   meESRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio gt 3p5 GeV";  
   meEBRecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);   
 
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio gt 3p5 GeV"; 
   meEERecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Barrel Unc RecSimHit Ratio";  
   meEBUnRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);   
 
   histo = "EcalRecHitsTask Endcap Unc RecSimHit Ratio"; 
   meEEUnRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Channel Status=10 11";  
   meEBRecHitSimHitRatio1011_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);   
 
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio Channel Status=10 11"; 
   meEERecHitSimHitRatio1011_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Channel Status=12";  
   meEBRecHitSimHitRatio12_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);   
 
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio Channel Status=12"; 
   meEERecHitSimHitRatio12_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Channel Status=13";  
   meEBRecHitSimHitRatio13_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);   
 
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio Channel Status=13"; 
   meEERecHitSimHitRatio13_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel Unc RecSimHit Ratio gt 3p5 GeV";  
   meEBUnRecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);   
 
   histo = "EcalRecHitsTask Endcap Unc RecSimHit Ratio gt 3p5 GeV"; 
   meEEUnRecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Barrel Rec E5x5";
   meEBe5x5_ = ibooker.book1D(histo.c_str(), histo.c_str(), 4000, 0., 400.);
 
   histo = "EcalRecHitsTask Barrel Rec E5x5 over Sim E5x5";
   meEBe5x5OverSimHits_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Barrel Rec E5x5 over gun energy";
   meEBe5x5OverGun_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Endcap Rec E5x5";
   meEEe5x5_ = ibooker.book1D(histo.c_str(), histo.c_str(), 4000, 0., 400.);
 
   histo = "EcalRecHitsTask Endcap Rec E5x5 over Sim E5x5";
   meEEe5x5OverSimHits_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Endcap Rec E5x5 over gun energy";
   meEEe5x5OverGun_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   meEBRecHitLog10Energy_ = ibooker.book1D( "EcalRecHitsTask Barrel Log10 Energy", "EcalRecHitsTask Barrel Log10 Energy", 90, -5., 4. ); 
   meEERecHitLog10Energy_ = ibooker.book1D( "EcalRecHitsTask Endcap Log10 Energy", "EcalRecHitsTask Endcap Log10 Energy", 90, -5., 4. ); 
   meESRecHitLog10Energy_ = ibooker.book1D( "EcalRecHitsTask Preshower Log10 Energy", "EcalRecHitsTask Preshower Log10 Energy", 90, -5., 4. ); 
   meEBRecHitLog10EnergyContr_ = ibooker.bookProfile( "EcalRecHits Barrel Log10En vs Hit Contribution", "EcalRecHits Barrel Log10En vs Hit Contribution", 90, -5., 4., 100, 0., 1. ); 
   meEERecHitLog10EnergyContr_ = ibooker.bookProfile( "EcalRecHits Endcap Log10En vs Hit Contribution", "EcalRecHits Endcap Log10En vs Hit Contribution", 90, -5., 4., 100, 0., 1. ); 
   meESRecHitLog10EnergyContr_ = ibooker.bookProfile( "EcalRecHits Preshower Log10En vs Hit Contribution", "EcalRecHits Preshower Log10En vs Hit Contribution", 90, -5., 4., 100, 0., 1. ); 
   meEBRecHitLog10Energy5x5Contr_ = ibooker.bookProfile( "EcalRecHits Barrel Log10En5x5 vs Hit Contribution", "EcalRecHits Barrel Log10En5x5 vs Hit Contribution", 90, -5., 4., 100, 0., 1. ); 
   meEERecHitLog10Energy5x5Contr_ = ibooker.bookProfile( "EcalRecHits Endcap Log10En5x5 vs Hit Contribution", "EcalRecHits Endcap Log10En5x5 vs Hit Contribution", 90, -5., 4., 100, 0., 1. ); 
   
 
   histo = "EB Occupancy Flag=5 6";  
   meEBRecHitsOccupancyFlag5_6_ = ibooker.book2D(histo, histo, 170, -85., 85., 360, 0., 360.);
   histo = "EB Occupancy Flag=8 9";  
   meEBRecHitsOccupancyFlag8_9_ = ibooker.book2D(histo, histo, 170, -85., 85., 360, 0., 360.);
 
   histo = "EE+ Occupancy Flag=5 6";  
   meEERecHitsOccupancyPlusFlag5_6_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
   histo = "EE- Occupancy Flag=5 6";  
   meEERecHitsOccupancyMinusFlag5_6_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
   histo = "EE+ Occupancy Flag=8 9";  
   meEERecHitsOccupancyPlusFlag8_9_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
   histo = "EE- Occupancy Flag=8 9";  
   meEERecHitsOccupancyMinusFlag8_9_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
 
 
   histo = "EcalRecHitsTask Barrel Reco Flags";  
   meEBRecHitFlags_ = ibooker.book1D(histo.c_str(), histo.c_str(), 10, 0., 10.);   
   histo = "EcalRecHitsTask Endcap Reco Flags";  
   meEERecHitFlags_ = ibooker.book1D(histo.c_str(), histo.c_str(), 10, 0., 10.);   
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio vs SimHit Flag=5 6";  
   meEBRecHitSimHitvsSimHitFlag5_6_ = ibooker.book2D(histo.c_str(), histo.c_str(), 80, 0., 2., 4000, 0., 400. );   
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio vs SimHit Flag=5 6"; 
   meEERecHitSimHitvsSimHitFlag5_6_ = ibooker.book2D(histo.c_str(), histo.c_str(), 80, 0., 2., 4000, 0., 400. );
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Flag=6";  
   meEBRecHitSimHitFlag6_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);   
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio Flag=6"; 
   meEERecHitSimHitFlag6_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Flag=7";  
   meEBRecHitSimHitFlag7_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);   
   histo = "EcalRecHitsTask Endcap RecSimHit Ratio Flag=7"; 
   meEERecHitSimHitFlag7_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
   histo = "EcalRecHitsTask Barrel 5x5 RecSimHit Ratio vs SimHit Flag=8";  
   meEB5x5RecHitSimHitvsSimHitFlag8_ = ibooker.book2D(histo.c_str(), histo.c_str(), 80, 0., 2., 4000, 0., 400. );   
 
 }

void EcalRecHitsValidation::findBarrelMatrix	(	int	nCellInEta,
		int	nCellInPhi,
		int	CentralEta,
		int	CentralPhi,
		int	CentralZ,
		MapType &	themap
	)

protected

Definition at line 705 of file EcalRecHitsValidation.cc.

References funct::abs(), mps_fire::i, and DetId::rawId().

                                                               {
   
   int goBackInEta = nCellInEta/2;
   int goBackInPhi = nCellInPhi/2;
   int matrixSize = nCellInEta*nCellInPhi; 
   crystalMatrix.clear();
   crystalMatrix.resize(matrixSize);
 
   int startEta  =  CentralZ*CentralEta - goBackInEta;
   int startPhi  =  CentralPhi - goBackInPhi;
   
   int i = 0 ;
   for ( int ieta = startEta; ieta < startEta+nCellInEta; ieta ++ ) {
     for( int iphi = startPhi; iphi < startPhi + nCellInPhi; iphi++ ) {
       uint32_t  index;
       if (abs(ieta) > 85 || abs(ieta)<1 ) { continue; }
       if (iphi< 1)      { index = EBDetId(ieta,iphi+360).rawId(); }
       else if(iphi>360) { index = EBDetId(ieta,iphi-360).rawId(); }
       else              { index = EBDetId(ieta,iphi).rawId();     }
       crystalMatrix[i++] = index;
     }
   }
   
 }

void EcalRecHitsValidation::findEndcapMatrix	(	int	nCellInX,
		int	nCellInY,
		int	CentralX,
		int	CentralY,
		int	CentralZ,
		MapType &	themap
	)

protected

Definition at line 732 of file EcalRecHitsValidation.cc.

References DetId::rawId(), and EEDetId::validDetId().

                                                                {
   int goBackInX = nCellInX/2;
   int goBackInY = nCellInY/2;
   crystalMatrix.clear();
 
    int startX  =  CentralX - goBackInX;
    int startY  =  CentralY - goBackInY;
 
    for ( int ix = startX; ix < startX+nCellInX; ix ++ ) {
 
       for( int iy = startY; iy < startY + nCellInY; iy++ ) {
 
         uint32_t index ;
 
     if(EEDetId::validDetId(ix,iy,CentralZ)) {
           index = EEDetId(ix,iy,CentralZ).rawId();
     }
     else { continue; }
         crystalMatrix.push_back(index);
       }
    }
 }

uint32_t EcalRecHitsValidation::getUnitWithMaxEnergy ( MapType & themap )

protected

Definition at line 687 of file EcalRecHitsValidation.cc.

References particleFlowClusterECALTimeSelected_cfi::maxEnergy.

                                                                     {
   
   //look for max
   uint32_t unitWithMaxEnergy = 0;
   float    maxEnergy = 0.;
   
   MapType::iterator iter;
   for (iter = themap.begin(); iter != themap.end(); iter++) {
     
     if (maxEnergy < (*iter).second) {
       maxEnergy = (*iter).second;       
       unitWithMaxEnergy = (*iter).first;
     }                           
   }
   
   return unitWithMaxEnergy;
 }

Member Data Documentation

std::vector<uint32_t> EcalRecHitsValidation::crystalMatrix

private

Definition at line 148 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalRecHitsValidation::EBHits_Token_

private

Definition at line 91 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EBrechitCollection_

private

Definition at line 79 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EBRecHitCollection> EcalRecHitsValidation::EBrechitCollection_Token_

private

Definition at line 86 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EBuncalibrechitCollection_

private

Definition at line 82 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EBUncalibratedRecHitCollection> EcalRecHitsValidation::EBuncalibrechitCollection_Token_

private

Definition at line 89 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalRecHitsValidation::EEHits_Token_

private

Definition at line 92 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EErechitCollection_

private

Definition at line 80 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EERecHitCollection> EcalRecHitsValidation::EErechitCollection_Token_

private

Definition at line 87 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EEuncalibrechitCollection_

private

Definition at line 83 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EEUncalibratedRecHitCollection> EcalRecHitsValidation::EEuncalibrechitCollection_Token_

private

Definition at line 90 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalRecHitsValidation::ESHits_Token_

private

Definition at line 93 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::ESrechitCollection_

private

Definition at line 81 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<ESRecHitCollection> EcalRecHitsValidation::ESrechitCollection_Token_

private

Definition at line 88 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::HepMCLabel

private

Definition at line 72 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<edm::HepMCProduct> EcalRecHitsValidation::HepMCLabel_Token_

private

Definition at line 85 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::hitsProducer_

private

Definition at line 73 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEB5x5RecHitSimHitvsSimHitFlag8_

private

Definition at line 141 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5_

private

Definition at line 113 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5OverGun_

private

Definition at line 115 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5OverSimHits_

private

Definition at line 114 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitFlags_

private

Definition at line 137 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy5x5Contr_

private

Definition at line 126 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy_

private

Definition at line 120 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10EnergyContr_

private

Definition at line 123 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag6_

private

Definition at line 139 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag7_

private

Definition at line 140 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio1011_

private

Definition at line 101 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio12_

private

Definition at line 103 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio13_

private

Definition at line 105 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio_

private

Definition at line 98 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatioGt35_

private

Definition at line 107 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitvsSimHitFlag5_6_

private

Definition at line 138 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag5_6_

private

Definition at line 129 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag8_9_

private

Definition at line 130 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatio_

private

Definition at line 109 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatioGt35_

private

Definition at line 111 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5_

private

Definition at line 116 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5OverGun_

private

Definition at line 118 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5OverSimHits_

private

Definition at line 117 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitFlags_

private

Definition at line 143 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy5x5Contr_

private

Definition at line 127 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy_

private

Definition at line 121 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10EnergyContr_

private

Definition at line 124 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag6_

private

Definition at line 145 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag7_

private

Definition at line 146 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio1011_

private

Definition at line 102 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio12_

private

Definition at line 104 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio13_

private

Definition at line 106 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio_

private

Definition at line 99 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatioGt35_

private

Definition at line 108 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitvsSimHitFlag5_6_

private

Definition at line 144 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag5_6_

private

Definition at line 133 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag8_9_

private

Definition at line 135 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag5_6_

private

Definition at line 132 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag8_9_

private

Definition at line 134 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatio_

private

Definition at line 110 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatioGt35_

private

Definition at line 112 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitLog10Energy_

private

Definition at line 122 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitLog10EnergyContr_

private

Definition at line 125 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitSimHitRatio_

private

Definition at line 100 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunEnergy_

private

Definition at line 95 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunEta_

private

Definition at line 96 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunPhi_

private

Definition at line 97 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::outputFile_

private

Definition at line 77 of file EcalRecHitsValidation.h.

bool EcalRecHitsValidation::verbose_

private

Definition at line 75 of file EcalRecHitsValidation.h.

Public Member Functions

Protected Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation