#include <DQMHcalPhiSymAlCaReco.h>

Inheritance diagram for DQMHcalPhiSymAlCaReco:

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

	~DQMHcalPhiSymAlCaReco () override

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

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

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

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

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

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

	~DQMEDAnalyzer () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

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

	ProducerBase ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

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

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

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

void	endRun (const edm::Run &r, const edm::EventSetup &c) override

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Attributes
int	eventCounter_

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

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

edm::EDGetTokenT< HBHERecHitCollection >	hbherecoMB
	object to monitor More...

edm::EDGetTokenT< HBHERecHitCollection >	hbherecoNoise

MonitorElement *	hFEDsize

edm::EDGetTokenT< HFRecHitCollection >	hfrecoMB

edm::EDGetTokenT< HFRecHitCollection >	hfrecoNoise

MonitorElement *	hHcalIsZS

int	hiDistr_r_nbin_

double	hiDistr_x_max_

double	hiDistr_x_min_

int	hiDistr_x_nbin_

double	hiDistr_y_max_

double	hiDistr_y_min_

int	hiDistr_y_nbin_

MonitorElement *	hiDistrHBHEsize1D_

MonitorElement *	hiDistrHFsize1D_

MonitorElement *	hiDistrMB2Min2D_

MonitorElement *	hiDistrMB2Pl2D_

MonitorElement *	hiDistrMBMin2D_

MonitorElement *	hiDistrMBPl2D_

MonitorElement *	hiDistrNoise2Min2D_

MonitorElement *	hiDistrNoise2Pl2D_

MonitorElement *	hiDistrNoiseMin2D_

MonitorElement *	hiDistrNoisePl2D_

MonitorElement *	hiDistrVarMBMin2D_

MonitorElement *	hiDistrVarMBPl2D_

MonitorElement *	hiDistrVarNoiseMin2D_

MonitorElement *	hiDistrVarNoisePl2D_

MonitorElement *	hL1Id

edm::InputTag	horecoMB

edm::InputTag	horecoNoise

double	ihbhe_size_

double	ihf_size_

unsigned int	period_

edm::EDGetTokenT< FEDRawDataCollection >	rawInLabel_

bool	saveToFile_
	Write to file. More...

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

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

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

DQM Source for phi symmetry stream

Author
Stefano Argiro' Andrea Gozzelino - Universita� e INFN Torino

Definition at line 24 of file DQMHcalPhiSymAlCaReco.h.

Constructor & Destructor Documentation

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

Definition at line 45 of file DQMHcalPhiSymAlCaReco.cc.

References fileName_, folderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hbherecoMB, hbherecoNoise, hfrecoMB, hfrecoNoise, hiDistr_r_nbin_, hiDistr_x_max_, hiDistr_x_min_, hiDistr_x_nbin_, hiDistr_y_max_, hiDistr_y_min_, hiDistr_y_nbin_, horecoMB, horecoNoise, ihbhe_size_, ihf_size_, period_, rawInLabel_, and saveToFile_.

                                                                         :
 eventCounter_(0)
 {
   //
   // Input from configurator file 
   //
   folderName_   = ps.getUntrackedParameter<string>("FolderName","ALCAStreamHcalPhiSym");
   
   hbherecoMB    = consumes<HBHERecHitCollection>(ps.getParameter<edm::InputTag>("hbheInputMB"));
   horecoMB      = ps.getParameter<edm::InputTag>("hoInputMB");
   hfrecoMB      = consumes<HFRecHitCollection>(ps.getParameter<edm::InputTag>("hfInputMB"));
   
   hbherecoNoise = consumes<HBHERecHitCollection>(ps.getParameter<edm::InputTag>("hbheInputNoise"));
   horecoNoise   = ps.getParameter<edm::InputTag>("hoInputNoise");
   hfrecoNoise   = consumes<HFRecHitCollection>(ps.getParameter<edm::InputTag>("hfInputNoise"));
 
   rawInLabel_   = consumes<FEDRawDataCollection>(ps.getParameter<edm::InputTag>("rawInputLabel"));
 
   period_       = ps.getParameter<unsigned int>("period") ;
 
   saveToFile_   = ps.getUntrackedParameter<bool>("SaveToFile",false);
   fileName_     = ps.getUntrackedParameter<string>("FileName","MonitorAlCaHcalPhiSym.root");
 
   // histogram parameters
 
   // Distribution of rechits in iPhi, iEta 
   hiDistr_y_nbin_ = ps.getUntrackedParameter<int>("hiDistr_y_nbin",72);
   hiDistr_y_min_  =  ps.getUntrackedParameter<double>("hiDistr_y_min",0.5);
   hiDistr_y_max_  =  ps.getUntrackedParameter<double>("hiDistr_y_max",72.5);
   hiDistr_x_nbin_ = ps.getUntrackedParameter<int>("hiDistr_x_nbin",41);
   hiDistr_x_min_  =  ps.getUntrackedParameter<double>("hiDistr_x_min",0.5);
   hiDistr_x_max_  =  ps.getUntrackedParameter<double>("hiDistr_x_max",41.5);
   // Check for NZS
   hiDistr_r_nbin_ = ps.getUntrackedParameter<int>("hiDistr_r_nbin",100);
   ihbhe_size_     = ps.getUntrackedParameter<double>("ihbhe_size_",5184.);
   ihf_size_       = ps.getUntrackedParameter<double>("ihf_size_",1728.);
     
 }

DQMHcalPhiSymAlCaReco::~DQMHcalPhiSymAlCaReco ( )

override

Definition at line 84 of file DQMHcalPhiSymAlCaReco.cc.

85 {}

Member Function Documentation

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

overrideprotected

Definition at line 304 of file DQMHcalPhiSymAlCaReco.cc.

References edm::SortedCollection< T, SORT >::begin(), FEDRawData::data(), HcalDetId::depth(), edm::SortedCollection< T, SORT >::end(), eventCounter_, FEDRawDataCollection::FEDData(), MonitorElement::Fill(), edm::Event::getByToken(), HcalDCCHeader::getSpigotData(), HcalDCCHeader::getSpigotPresent(), hbherecoMB, hbherecoNoise, hFEDsize, hfrecoMB, hfrecoNoise, hHcalIsZS, hiDistrHBHEsize1D_, hiDistrHFsize1D_, hiDistrMB2Min2D_, hiDistrMB2Pl2D_, hiDistrMBMin2D_, hiDistrMBPl2D_, hiDistrNoise2Min2D_, hiDistrNoise2Pl2D_, hiDistrNoiseMin2D_, hiDistrNoisePl2D_, hL1Id, mps_fire::i, HcalDetId::ieta(), ihbhe_size_, ihf_size_, HcalDetId::iphi(), HcalHTRData::isUnsuppressed(), edm::HandleBase::isValid(), gen::k, LogDebug, FEDHeader::lvl1ID(), FEDNumbering::MAXHCALFEDID, FEDNumbering::MINHCALFEDID, period_, funct::pow(), edm::Handle< T >::product(), rawInLabel_, mathSSE::return(), FEDRawData::size(), edm::SortedCollection< T, SORT >::size(), and HcalDCCHeader::SPIGOT_COUNT.

                                              {  
  
 
   eventCounter_++;
  
   edm::Handle<FEDRawDataCollection> rawIn;
   iEvent.getByToken(rawInLabel_,rawIn);
 
   if(!rawIn.isValid()){
      LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
       return ;
   }
 
   //get HCAL FEDs:
   std::vector<int> selFEDs;
   for (int i=FEDNumbering::MINHCALFEDID; i<=FEDNumbering::MAXHCALFEDID; i++)
     {
       selFEDs.push_back(i);
     }
 
   //  std::cout<<" Size of FED "<<selFEDs.size()<<std::endl;
 
   const FEDRawDataCollection *rdc=rawIn.product(); 
 
   bool hcalIsZS = false;
   int lvl1ID = 0;
   bool lvl1IDFound = false;
   for (unsigned int k=0; k<selFEDs.size(); k++) {
     const FEDRawData & fedData = rdc->FEDData(selFEDs[k]);
     // std::cout<<fedData.size()*std::pow(1024.,-1)<<std::endl;
     hFEDsize->Fill(fedData.size()*std::pow(1024.,-1),1);
     
     // get HCAL DCC Header for each FEDRawData
     const HcalDCCHeader* dccHeader=(const HcalDCCHeader*)(fedData.data());
     if (dccHeader) {
       // walk through the HTR data...
       HcalHTRData htr;
 
       int nspigot =0; 
       for (int spigot=0; spigot<HcalDCCHeader::SPIGOT_COUNT; spigot++) {    
        nspigot++;
 
         if (!dccHeader->getSpigotPresent(spigot)) continue;
 
         // Load the given decoder with the pointer and length from this spigot.
         dccHeader->getSpigotData(spigot,htr, fedData.size()); 
 
         if(k != 20 && nspigot !=14 ) {      
           if ( !htr.isUnsuppressed() ) { hcalIsZS = true; }
         }
       } 
     }
  
     // try to get the lvl1ID from the HCAL fed
     if (!lvl1IDFound && (fedData.size() > 0)) { 
       // get FED Header for FEDRawData
       FEDHeader fedHeader(fedData.data());
       lvl1ID = fedHeader.lvl1ID();
       lvl1IDFound = true;
     }
   } // loop over HcalFEDs
   
   hHcalIsZS->Fill( hcalIsZS );
   hL1Id->Fill( lvl1ID%period_ );
  
   edm::Handle<HBHERecHitCollection> hbheNS;
   iEvent.getByToken(hbherecoNoise, hbheNS);
   
   if(!hbheNS.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     return ;
   }
   
   edm::Handle<HBHERecHitCollection> hbheMB;
   iEvent.getByToken(hbherecoMB, hbheMB);
   
   if(!hbheMB.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     return ;
   }
   
   edm::Handle<HFRecHitCollection> hfNS;
   iEvent.getByToken(hfrecoNoise, hfNS);
   
   if(!hfNS.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     return ;
   }
   
   edm::Handle<HFRecHitCollection> hfMB;
   iEvent.getByToken(hfrecoMB, hfMB);
   
   if(!hfMB.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     return ;
   }
   
   const HBHERecHitCollection HithbheNS = *(hbheNS.product());
   
   hiDistrHBHEsize1D_->Fill(HithbheNS.size()/ihbhe_size_);
   
      
   for(HBHERecHitCollection::const_iterator hbheItr=HithbheNS.begin(); hbheItr!=HithbheNS.end(); hbheItr++)
         {
              DetId id = (*hbheItr).detid(); 
              HcalDetId hid=HcalDetId(id);
                  
          if(hid.depth() == 1) {
                  if( hid.ieta() > 0 ) {
          hiDistrNoisePl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
                  hiDistrNoise2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
                  } else {
          hiDistrNoiseMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
                  hiDistrNoise2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
          }
          }
         }
 
   const HBHERecHitCollection HithbheMB = *(hbheMB.product());
 
   for(HBHERecHitCollection::const_iterator hbheItr=HithbheMB.begin(); hbheItr!=HithbheMB.end(); hbheItr++)
         {
              DetId id = (*hbheItr).detid(); 
              HcalDetId hid=HcalDetId(id);
                  
          if(hid.depth() == 1) {
                  if( hid.ieta() > 0 ) {
          hiDistrMBPl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
                  hiDistrMB2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
                  } else {
          hiDistrMBMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
                  hiDistrMB2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
          }
          }
 
         }
 
   const HFRecHitCollection HithfNS = *(hfNS.product());
 
   hiDistrHFsize1D_->Fill(HithfNS.size()/ihf_size_);
   
   for(HFRecHitCollection::const_iterator hbheItr=HithfNS.begin(); hbheItr!=HithfNS.end(); hbheItr++)
         {
     
              DetId id = (*hbheItr).detid(); 
              HcalDetId hid=HcalDetId(id);
                  
          if(hid.depth() == 1) {
                  if( hid.ieta() > 0 ) {
          hiDistrNoisePl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
                  hiDistrNoise2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
                  } else {
          hiDistrNoiseMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
                  hiDistrNoise2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
          }
          }
     
         }
 
   const HFRecHitCollection HithfMB = *(hfMB.product());
  
   for(HFRecHitCollection::const_iterator hbheItr=HithfMB.begin(); hbheItr!=HithfMB.end(); hbheItr++)
         {
              DetId id = (*hbheItr).detid(); 
              HcalDetId hid=HcalDetId(id);
                  
          if(hid.depth() == 1) {
                  if( hid.ieta() > 0 ) {
          hiDistrMBPl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
                  hiDistrMB2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
                  } else {
          hiDistrMBMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
                  hiDistrMB2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
          }
          }
         }   
     
     
 } //analyze

void DQMHcalPhiSymAlCaReco::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	irun,
		edm::EventSetup const &	isetup
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 88 of file DQMHcalPhiSymAlCaReco.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), eventCounter_, folderName_, hFEDsize, hHcalIsZS, hiDistr_r_nbin_, hiDistr_x_max_, hiDistr_x_min_, hiDistr_x_nbin_, hiDistr_y_max_, hiDistr_y_min_, hiDistr_y_nbin_, hiDistrHBHEsize1D_, hiDistrHFsize1D_, hiDistrMB2Min2D_, hiDistrMB2Pl2D_, hiDistrMBMin2D_, hiDistrMBPl2D_, hiDistrNoise2Min2D_, hiDistrNoise2Pl2D_, hiDistrNoiseMin2D_, hiDistrNoisePl2D_, hiDistrVarMBMin2D_, hiDistrVarMBPl2D_, hiDistrVarNoiseMin2D_, hiDistrVarNoisePl2D_, hL1Id, period_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

                                                                                                                        {
  
   // create and cd into new folder
   ibooker.setCurrentFolder(folderName_);
 
   eventCounter_ = 0;
 
   hFEDsize = ibooker.book1D("hFEDsize","HCAL FED size (kB)",200,-0.5,20.5);
   hFEDsize->setAxisTitle("kB",1);
 
   hHcalIsZS = ibooker.book1D("hHcalIsZS", "Hcal Is ZS", 4, -1.5, 2.5);
   hHcalIsZS->setBinLabel(2, "NZS");
   hHcalIsZS->setBinLabel(3, "ZS");
 
   char hname[50];
   sprintf(hname, "L1 Event Number %% %i", period_);
   hL1Id = ibooker.book1D("hL1Id", hname,4200,-99.5,4099.5);
   hL1Id->setAxisTitle(hname);
   
 
   // book some histograms 1D
   double xmin = 0.1;
   double xmax = 1.1;
   hiDistrHBHEsize1D_ =
     ibooker.book1D("DistrHBHEsize","Size of HBHE Collection",
                   hiDistr_r_nbin_,
                   xmin,
                   xmax
                 );
   hiDistrHFsize1D_ =
     ibooker.book1D("DistrHFsize","Size of HF Collection",
                   hiDistr_r_nbin_,
                   xmin,
                   xmax
                 );
 
 
   // First moment
   hiDistrMBPl2D_ = 
     ibooker.book2D("MBdepthPl1", "iphi- +ieta signal distribution at depth1",
          hiDistr_x_nbin_, 
          hiDistr_x_min_,
          hiDistr_x_max_,
          hiDistr_y_nbin_, 
          hiDistr_y_min_,
          hiDistr_y_max_
          );
 
   hiDistrMBPl2D_->setAxisTitle("i#phi ", 2);
   hiDistrMBPl2D_->setAxisTitle("i#eta ", 1);
 
 
   hiDistrNoisePl2D_ = 
     ibooker.book2D("NoisedepthPl1", "iphi-ieta noise distribution at depth1",
          hiDistr_x_nbin_+1, 
          hiDistr_x_min_-1.,
          hiDistr_x_max_,
          hiDistr_y_nbin_+1, 
          hiDistr_y_min_-1.,
          hiDistr_y_max_
          );
 
   hiDistrNoisePl2D_->setAxisTitle("i#phi ", 2);
   hiDistrNoisePl2D_->setAxisTitle("i#eta ", 1);
 // Second moment
   hiDistrMB2Pl2D_ =
     ibooker.book2D("MB2depthPl1", "iphi- +ieta signal distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrMB2Pl2D_->setAxisTitle("i#phi ", 2);
   hiDistrMB2Pl2D_->setAxisTitle("i#eta ", 1);
 
 
   hiDistrNoise2Pl2D_ =
     ibooker.book2D("Noise2depthPl1", "iphi-ieta noise distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrNoise2Pl2D_->setAxisTitle("i#phi ", 2);
   hiDistrNoise2Pl2D_->setAxisTitle("i#eta ", 1);
 
 // Variance
   hiDistrVarMBPl2D_ =
     ibooker.book2D("VarMBdepthPl1", "iphi- +ieta signal distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrVarMBPl2D_->setAxisTitle("i#phi ", 2);
   hiDistrVarMBPl2D_->setAxisTitle("i#eta ", 1);
 
 
   hiDistrVarNoisePl2D_ =
     ibooker.book2D("VarNoisedepthPl1", "iphi-ieta noise distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrVarNoisePl2D_->setAxisTitle("i#phi ", 2);
   hiDistrVarNoisePl2D_->setAxisTitle("i#eta ", 1);
 
 //==================================================================================
 // First moment
   hiDistrMBMin2D_ = 
     ibooker.book2D("MBdepthMin1", "iphi- +ieta signal distribution at depth1",
          hiDistr_x_nbin_, 
          hiDistr_x_min_,
          hiDistr_x_max_,
          hiDistr_y_nbin_, 
          hiDistr_y_min_,
          hiDistr_y_max_
          );
 
   hiDistrMBMin2D_->setAxisTitle("i#phi ", 2);
   hiDistrMBMin2D_->setAxisTitle("i#eta ", 1);
 
 
   hiDistrNoiseMin2D_ = 
     ibooker.book2D("NoisedepthMin1", "iphi-ieta noise distribution at depth1",
          hiDistr_x_nbin_, 
          hiDistr_x_min_,
          hiDistr_x_max_,
          hiDistr_y_nbin_, 
          hiDistr_y_min_,
          hiDistr_y_max_
          );
 
   hiDistrNoiseMin2D_->setAxisTitle("i#phi ", 2);
   hiDistrNoiseMin2D_->setAxisTitle("i#eta ", 1);
 // Second moment
   hiDistrMB2Min2D_ =
     ibooker.book2D("MB2depthMin1", "iphi- +ieta signal distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrMB2Min2D_->setAxisTitle("i#phi ", 2);
   hiDistrMB2Min2D_->setAxisTitle("i#eta ", 1);
 
 
   hiDistrNoise2Min2D_ =
     ibooker.book2D("Noise2depthMin1", "iphi-ieta noise distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrNoise2Min2D_->setAxisTitle("i#phi ", 2);
   hiDistrNoise2Min2D_->setAxisTitle("i#eta ", 1);
 
 // Variance
   hiDistrVarMBMin2D_ =
     ibooker.book2D("VarMBdepthMin1", "iphi- +ieta signal distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrVarMBMin2D_->setAxisTitle("i#phi ", 2);
   hiDistrVarMBMin2D_->setAxisTitle("i#eta ", 1);
 
 
   hiDistrVarNoiseMin2D_ =
     ibooker.book2D("VarNoisedepthMin1", "iphi-ieta noise distribution at depth1",
                  hiDistr_x_nbin_,
                  hiDistr_x_min_,
                  hiDistr_x_max_,
                  hiDistr_y_nbin_,
                  hiDistr_y_min_,
                  hiDistr_y_max_
                  );
 
   hiDistrVarNoiseMin2D_->setAxisTitle("i#phi ", 2);
   hiDistrVarNoiseMin2D_->setAxisTitle("i#eta ", 1);
 
 
 }

void DQMHcalPhiSymAlCaReco::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overrideprotected

Definition at line 490 of file DQMHcalPhiSymAlCaReco.cc.

References eventCounter_, MonitorElement::getBinContent(), hiDistr_x_nbin_, hiDistr_y_nbin_, hiDistrMB2Min2D_, hiDistrMB2Pl2D_, hiDistrMBMin2D_, hiDistrMBPl2D_, hiDistrNoise2Min2D_, hiDistrNoise2Pl2D_, hiDistrNoiseMin2D_, hiDistrNoisePl2D_, hiDistrVarMBMin2D_, hiDistrVarMBPl2D_, hiDistrVarNoiseMin2D_, hiDistrVarNoisePl2D_, gen::k, and MonitorElement::setBinContent().

                                                                          {
 // Keep Variances
   if(eventCounter_ > 0) {
   for(int k=0; k<=hiDistr_x_nbin_;k++)
   {
     for(int j=0; j<=hiDistr_y_nbin_;j++)
     {
 // First moment
        float cc1=hiDistrMBPl2D_->getBinContent(k,j);
        cc1 = cc1 * 1./eventCounter_;
        float cc2=hiDistrNoisePl2D_->getBinContent(k,j);
        cc2 = cc2 * 1./eventCounter_;
        float cc3=hiDistrMBMin2D_->getBinContent(k,j);
        cc3 = cc3 * 1./eventCounter_;
        float cc4=hiDistrNoiseMin2D_->getBinContent(k,j);
        cc4 = cc4 * 1./eventCounter_;
 // Second moment
        float cc11=hiDistrMB2Pl2D_->getBinContent(k,j);
        cc11 = cc11 * 1./eventCounter_;
        hiDistrVarMBPl2D_->setBinContent(k,j,cc11-cc1*cc1);
        float cc22=hiDistrNoise2Pl2D_->getBinContent(k,j);
        cc22 = cc22 * 1./eventCounter_;
        hiDistrVarNoisePl2D_->setBinContent(k,j,cc22-cc2*cc2);
        float cc33=hiDistrMB2Min2D_->getBinContent(k,j);
        cc33 = cc33 * 1./eventCounter_;
        hiDistrVarMBMin2D_->setBinContent(k,j,cc33-cc3*cc3);
        float cc44=hiDistrNoise2Min2D_->getBinContent(k,j);
        cc44 = cc44 * 1./eventCounter_;
        hiDistrVarNoiseMin2D_->setBinContent(k,j,cc44-cc4*cc4);
     }     
   }
   }
 }