#include <DQMHcalPhiSymAlCaReco.h>

Inheritance diagram for DQMHcalPhiSymAlCaReco:

Public Member Functions
	DQMHcalPhiSymAlCaReco (const edm::ParameterSet &)
	~DQMHcalPhiSymAlCaReco ()
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)
Private Attributes
DQMStore *	dbe_
int	eventCounter_
std::string	fileName_
	Output file name if required.
std::string	folderName_
	DQM folder name.
edm::InputTag	hbherecoMB
	object to monitor
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.

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_, cmsCodeRules::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.

{}

Member Function Documentation

void DQMHcalPhiSymAlCaReco::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

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, MultipleCompare::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_, 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 ) [protected, virtual]

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
	)		`[protected, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 306 of file DQMHcalPhiSymAlCaReco.cc.

                                {
  
}

void DQMHcalPhiSymAlCaReco::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 301 of file DQMHcalPhiSymAlCaReco.cc.

                                                                               {
//   eventCounter_ = 0;
}

void DQMHcalPhiSymAlCaReco::endJob ( void ) [protected, virtual]

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
	)		`[protected, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 497 of file DQMHcalPhiSymAlCaReco.cc.

{
}

void DQMHcalPhiSymAlCaReco::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

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);
    }     
  }
  }
}