#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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)

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_

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 45 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 116 of file EcalRecHitsValidation.cc.

116 {}

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 276 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(), 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 &
	)

overrideprotected

Definition at line 118 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 787 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 817 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 771 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 140 of file EcalRecHitsValidation.h.

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

private

Definition at line 83 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EBrechitCollection_

private

Definition at line 71 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EBRecHitCollection> EcalRecHitsValidation::EBrechitCollection_Token_

private

Definition at line 78 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EBuncalibrechitCollection_

private

Definition at line 74 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EBUncalibratedRecHitCollection> EcalRecHitsValidation::EBuncalibrechitCollection_Token_

private

Definition at line 81 of file EcalRecHitsValidation.h.

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

private

Definition at line 84 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EErechitCollection_

private

Definition at line 72 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EERecHitCollection> EcalRecHitsValidation::EErechitCollection_Token_

private

Definition at line 79 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EEuncalibrechitCollection_

private

Definition at line 75 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EEUncalibratedRecHitCollection> EcalRecHitsValidation::EEuncalibrechitCollection_Token_

private

Definition at line 82 of file EcalRecHitsValidation.h.

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

private

Definition at line 85 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::ESrechitCollection_

private

Definition at line 73 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<ESRecHitCollection> EcalRecHitsValidation::ESrechitCollection_Token_

private

Definition at line 80 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::HepMCLabel

private

Definition at line 64 of file EcalRecHitsValidation.h.

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

private

Definition at line 77 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::hitsProducer_

private

Definition at line 65 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEB5x5RecHitSimHitvsSimHitFlag8_

private

Definition at line 133 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5_

private

Definition at line 105 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5OverGun_

private

Definition at line 107 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5OverSimHits_

private

Definition at line 106 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitFlags_

private

Definition at line 129 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy5x5Contr_

private

Definition at line 118 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy_

private

Definition at line 112 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10EnergyContr_

private

Definition at line 115 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag6_

private

Definition at line 131 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag7_

private

Definition at line 132 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio1011_

private

Definition at line 93 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio12_

private

Definition at line 95 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio13_

private

Definition at line 97 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio_

private

Definition at line 90 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatioGt35_

private

Definition at line 99 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitvsSimHitFlag5_6_

private

Definition at line 130 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag5_6_

private

Definition at line 121 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag8_9_

private

Definition at line 122 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatio_

private

Definition at line 101 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatioGt35_

private

Definition at line 103 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5_

private

Definition at line 108 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5OverGun_

private

Definition at line 110 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5OverSimHits_

private

Definition at line 109 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitFlags_

private

Definition at line 135 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy5x5Contr_

private

Definition at line 119 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy_

private

Definition at line 113 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10EnergyContr_

private

Definition at line 116 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag6_

private

Definition at line 137 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag7_

private

Definition at line 138 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio1011_

private

Definition at line 94 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio12_

private

Definition at line 96 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio13_

private

Definition at line 98 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio_

private

Definition at line 91 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatioGt35_

private

Definition at line 100 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitvsSimHitFlag5_6_

private

Definition at line 136 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag5_6_

private

Definition at line 125 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag8_9_

private

Definition at line 127 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag5_6_

private

Definition at line 124 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag8_9_

private

Definition at line 126 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatio_

private

Definition at line 102 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatioGt35_

private

Definition at line 104 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitLog10Energy_

private

Definition at line 114 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitLog10EnergyContr_

private

Definition at line 117 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitSimHitRatio_

private

Definition at line 92 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunEnergy_

private

Definition at line 87 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunEta_

private

Definition at line 88 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunPhi_

private

Definition at line 89 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::outputFile_

private

Definition at line 69 of file EcalRecHitsValidation.h.

bool EcalRecHitsValidation::verbose_

private

Definition at line 67 of file EcalRecHitsValidation.h.

Public Member Functions

Protected Member Functions

Private Types

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation