#include <DQMHcalPhiSymAlCaReco.h>

Inheritance diagram for DQMHcalPhiSymAlCaReco:

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

	~DQMHcalPhiSymAlCaReco ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

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

void	beginJob ()

void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context)

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

void	endJob ()

void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)

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

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
DQMStore *	dbe_

int	eventCounter_

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

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

edm::InputTag	hbherecoMB
	object to monitor More...

edm::InputTag	hbherecoNoise

MonitorElement *	hFEDsize

edm::InputTag	hfrecoMB

edm::InputTag	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::InputTag	rawInLabel_

bool	saveToFile_
	Write to file. More...

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

DQM Source for phi symmetry stream

Date:
2009/12/14 22:22:19

Revision:
1.2

Author
Stefano Argiro' Andrea Gozzelino - Universita� e INFN Torino
DQM Source for phi symmetry stream

Date:
2009/12/14 22:22:19

Revision:
1.6

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

Definition at line 22 of file DQMHcalPhiSymAlCaReco.h.

Constructor & Destructor Documentation

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

Definition at line 49 of file DQMHcalPhiSymAlCaReco.cc.

References dbe_, 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_, cppFunctionSkipper::operator, period_, rawInLabel_, and saveToFile_.

                                                                         :
 eventCounter_(0)
 {
   dbe_ = Service<DQMStore>().operator->();
   //
   // Input from configurator file 
   //
   folderName_   = ps.getUntrackedParameter<string>("FolderName","ALCAStreamHcalPhiSym");
   
   hbherecoMB    = ps.getParameter<edm::InputTag>("hbheInputMB");
   horecoMB      = ps.getParameter<edm::InputTag>("hoInputMB");
   hfrecoMB      = ps.getParameter<edm::InputTag>("hfInputMB");
   
   hbherecoNoise = ps.getParameter<edm::InputTag>("hbheInputNoise");
   horecoNoise   = ps.getParameter<edm::InputTag>("hoInputNoise");
   hfrecoNoise   = ps.getParameter<edm::InputTag>("hfInputNoise");
 
   rawInLabel_   = 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 ( )

Definition at line 89 of file DQMHcalPhiSymAlCaReco.cc.

90 {}

Member Function Documentation

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

protectedvirtual

Level-1 event number generated by the TTC system

Implements edm::EDAnalyzer.

Definition at line 313 of file DQMHcalPhiSymAlCaReco.cc.

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

                                              {  
  
 
   eventCounter_++;
  
   edm::Handle<FEDRawDataCollection> rawIn;
   iEvent.getByLabel(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 ;
   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) continue;
 
       // 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; }
         }
       } 
     
     } // loop over HcalFEDs
   
   hHcalIsZS->Fill( hcalIsZS );
 
   // get Trigger FED-Id
   const FEDRawData& fedData = rdc->FEDData(FEDNumbering::MINTriggerGTPFEDID) ;
   FEDHeader header(fedData.data()) ;
 
   hL1Id->Fill( (header.lvl1ID())%period_ );
  
   edm::Handle<HBHERecHitCollection> hbheNS;
   iEvent.getByLabel(hbherecoNoise, hbheNS);
   
   if(!hbheNS.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     return ;
   }
   
   edm::Handle<HBHERecHitCollection> hbheMB;
   iEvent.getByLabel(hbherecoMB, hbheMB);
   
   if(!hbheMB.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     return ;
   }
   
   edm::Handle<HFRecHitCollection> hfNS;
   iEvent.getByLabel(hfrecoNoise, hfNS);
   
   if(!hfNS.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     return ;
   }
   
   edm::Handle<HFRecHitCollection> hfMB;
   iEvent.getByLabel(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::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 93 of file DQMHcalPhiSymAlCaReco.cc.

References DQMStore::book1D(), DQMStore::book2D(), dbe_, 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(), and DQMStore::setCurrentFolder().

                                     {
  
   // create and cd into new folder
   dbe_->setCurrentFolder(folderName_);
 
   eventCounter_ = 0;
 
   hFEDsize = dbe_->book1D("hFEDsize","HCAL FED size (kB)",200,-0.5,20.5);
   hFEDsize->setAxisTitle("kB",1);
 
   hHcalIsZS = dbe_->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 = dbe_->book1D("hL1Id", hname,4200,-99.5,4099.5);
   hL1Id->setAxisTitle(hname);
   
 
   // book some histograms 1D
   double xmin = 0.1;
   double xmax = 1.1;
   hiDistrHBHEsize1D_ =
     dbe_->book1D("DistrHBHEsize","Size of HBHE Collection",
                   hiDistr_r_nbin_,
                   xmin,
                   xmax
                 );
   hiDistrHFsize1D_ =
     dbe_->book1D("DistrHFsize","Size of HF Collection",
                   hiDistr_r_nbin_,
                   xmin,
                   xmax
                 );
 
 
   // First moment
   hiDistrMBPl2D_ = 
     dbe_->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_ = 
     dbe_->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_ =
     dbe_->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_ =
     dbe_->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_ =
     dbe_->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_ =
     dbe_->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_ = 
     dbe_->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_ = 
     dbe_->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_ =
     dbe_->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_ =
     dbe_->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_ =
     dbe_->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_ =
     dbe_->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::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	context
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 306 of file DQMHcalPhiSymAlCaReco.cc.

                                 {
   
 }

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 301 of file DQMHcalPhiSymAlCaReco.cc.

                                                                                {
 //   eventCounter_ = 0;
 }

void DQMHcalPhiSymAlCaReco::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 535 of file DQMHcalPhiSymAlCaReco.cc.

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

                                   {
   if (saveToFile_) {
 
   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_;
        hiDistrMBPl2D_->setBinContent(k,j,cc1); 
        float cc2=hiDistrNoisePl2D_->getBinContent(k,j);
        cc2 = cc2 * 1./eventCounter_;
        hiDistrNoisePl2D_->setBinContent(k,j,cc2);
        float cc3=hiDistrMBMin2D_->getBinContent(k,j);
        cc3 = cc3 * 1./eventCounter_;
        hiDistrMBMin2D_->setBinContent(k,j,cc3);
        float cc4=hiDistrNoiseMin2D_->getBinContent(k,j);
        cc4 = cc4 * 1./eventCounter_;
        hiDistrNoiseMin2D_->setBinContent(k,j,cc4);
 // Second moment
        float cc11=hiDistrMB2Pl2D_->getBinContent(k,j);
        cc11 = cc11 * 1./eventCounter_;
        hiDistrMB2Pl2D_->setBinContent(k,j,cc11);
        hiDistrVarMBPl2D_->setBinContent(k,j,cc11-cc1*cc1);
        float cc22=hiDistrNoise2Pl2D_->getBinContent(k,j);
        cc22 = cc22 * 1./eventCounter_;
        hiDistrNoise2Pl2D_->setBinContent(k,j,cc22);
        hiDistrVarNoisePl2D_->setBinContent(k,j,cc22-cc2*cc2);
        float cc33=hiDistrMB2Min2D_->getBinContent(k,j);
        cc33 = cc33 * 1./eventCounter_;
        hiDistrMB2Min2D_->setBinContent(k,j,cc33);
        hiDistrVarMBMin2D_->setBinContent(k,j,cc33-cc3*cc3);
        float cc44=hiDistrNoise2Min2D_->getBinContent(k,j);
        cc44 = cc44 * 1./eventCounter_;
        hiDistrNoise2Min2D_->setBinContent(k,j,cc44);
        hiDistrVarNoiseMin2D_->setBinContent(k,j,cc44-cc4*cc4);
     }
   }
      dbe_->save(fileName_);
   }
 }

void DQMHcalPhiSymAlCaReco::endLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	c
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 497 of file DQMHcalPhiSymAlCaReco.cc.

498 {

499 }

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 501 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_, j, 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);
     }     
   }
   }
 }