#include <HcalPedestalAnalysis.h>

Public Member Functions
int	done (const HcalPedestals fInputPedestals, const HcalPedestalWidths fInputWidths, HcalPedestals fOutputPedestals, HcalPedestalWidths fOutputWidths)

	HcalPedestalAnalysis (const edm::ParameterSet &ps)
	Constructor. More...

void	processEvent (const HBHEDigiCollection &hbhe, const HODigiCollection &ho, const HFDigiCollection &hf, const HcalDbService &cond)

void	SampleAnalysis ()

void	setTopology (const HcalTopology *htopo)

void	setup (const std::string &m_outputFileROOT)

	~HcalPedestalAnalysis ()
	Destructor. More...

Static Public Member Functions
static int	HcalPedVal (int nstat[4], const HcalPedestals fRefPedestals, const HcalPedestalWidths fRefPedestalWidths, HcalPedestals fRawPedestals, HcalPedestalWidths fRawPedestalWidths, HcalPedestals fValPedestals, HcalPedestalWidths fValPedestalWidths)

Private Types
typedef std::pair< TH1F , std::pair< std::map< int, std::vector< double > >, std::vector< TH1F > > >	PEDBUNCH

Private Member Functions
void	AllChanHists (const HcalDetId detid, const HcalQIESample &qie0, const HcalQIESample &qie1, const HcalQIESample &qie2, const HcalQIESample &qie3, const HcalQIESample &qie4, const HcalQIESample &qie5, std::map< HcalDetId, std::map< int, PEDBUNCH > > &toolT)

void	GetPedConst (std::map< HcalDetId, std::map< int, PEDBUNCH > > &toolT, TH1F PedMeans, TH1F PedWidths)

void	per2CapsHists (int flag, int id, const HcalDetId detid, const HcalQIESample &qie1, const HcalQIESample &qie2, std::map< HcalDetId, std::map< int, PEDBUNCH > > &toolT, const HcalDbService &cond)

void	Trendings (std::map< HcalDetId, std::map< int, PEDBUNCH > > &toolT, TH1F Chi2, TH1F CapidAverage, TH1F *CapidChi2)

Private Attributes
std::map< HcalDetId, std::map< int, PEDBUNCH > >::iterator	_meot

int	evt

int	evt_curr

HcalPedestals *	fRawPedestals

HcalPedestalWidths *	fRawPedestalWidths

const HcalPedestals *	fRefPedestals

const HcalPedestalWidths *	fRefPedestalWidths

const HcalTopology *	fTopology

HcalPedestals *	fValPedestals

HcalPedestalWidths *	fValPedestalWidths

struct {
TH1F * ALLPEDS

TH1F * CAPID_AVERAGE

TH1F * CAPID_CHI2

TH1F * CHI2

TH1F * PEDMEAN

TH1F * PEDRMS

std::map< HcalDetId, std::map< int, PEDBUNCH > > PEDTRENDS

}	hbHists

struct {
TH1F * ALLPEDS

TH1F * CAPID_AVERAGE

TH1F * CAPID_CHI2

TH1F * CHI2

TH1F * PEDMEAN

TH1F * PEDRMS

std::map< HcalDetId, std::map< int, PEDBUNCH > > PEDTRENDS

}	hfHists

struct {
TH1F * ALLPEDS

TH1F * CAPID_AVERAGE

TH1F * CAPID_CHI2

TH1F * CHI2

TH1F * PEDMEAN

TH1F * PEDRMS

std::map< HcalDetId, std::map< int, PEDBUNCH > > PEDTRENDS

}	hoHists

int	m_AllPedsOK

const HcalQIECoder *	m_coder

int	m_endTS

TFile *	m_file

int	m_hiSaveflag

std::ofstream	m_logFile

int	m_nevtsample

std::string	m_outputFileMean

std::string	m_outputFileROOT

std::string	m_outputFileWidth

int	m_pedsinADC

int	m_pedValflag

const HcalQIEShape *	m_shape

int	m_startTS

float	m_stat [4]

int	sample

std::vector< bool >	state

Static Private Attributes
static const int	fitflag = 0

Detailed Description

Definition at line 31 of file HcalPedestalAnalysis.h.

Member Typedef Documentation

◆ PEDBUNCH

typedef std::pair<TH1F*, std::pair<std::map<int, std::vector<double> >, std::vector<TH1F*> > > HcalPedestalAnalysis::PEDBUNCH

private

Definition at line 82 of file HcalPedestalAnalysis.h.

Constructor & Destructor Documentation

◆ HcalPedestalAnalysis()

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

Constructor.

Definition at line 16 of file HcalPedestalAnalysis.cc.

References gather_cfg::cout, evt, edm::ParameterSet::getUntrackedParameter(), hbHists, hfHists, hoHists, mps_fire::i, dqmdumpme::k, m_AllPedsOK, m_coder, m_endTS, m_file, m_hiSaveflag, m_nevtsample, m_outputFileMean, m_outputFileROOT, m_outputFileWidth, m_pedsinADC, m_pedValflag, m_shape, m_startTS, m_stat, sample, and state.

     : fRefPedestals(nullptr),
       fRefPedestalWidths(nullptr),
       fRawPedestals(nullptr),
       fRawPedestalWidths(nullptr),
       fValPedestals(nullptr),
       fValPedestalWidths(nullptr),
       fTopology(nullptr) {
   m_coder = nullptr;
   m_shape = nullptr;
   evt = 0;
   sample = 0;
   m_file = nullptr;
   m_AllPedsOK = 0;
   for (int i = 0; i < 4; i++)
     m_stat[i] = 0;
   for (int k = 0; k < 4; k++)
     state.push_back(true);
 
   // user cfg parameters
   m_outputFileMean = ps.getUntrackedParameter<string>("outputFileMeans", "");
   if (!m_outputFileMean.empty()) {
     cout << "Hcal pedestal means will be saved to " << m_outputFileMean.c_str() << endl;
   }
   m_outputFileWidth = ps.getUntrackedParameter<string>("outputFileWidths", "");
   if (!m_outputFileWidth.empty()) {
     cout << "Hcal pedestal widths will be saved to " << m_outputFileWidth.c_str() << endl;
   }
   m_outputFileROOT = ps.getUntrackedParameter<string>("outputFileHist", "");
   if (!m_outputFileROOT.empty()) {
     cout << "Hcal pedestal histograms will be saved to " << m_outputFileROOT.c_str() << endl;
   }
   m_nevtsample = ps.getUntrackedParameter<int>("nevtsample", 0);
   // for compatibility with previous versions
   if (m_nevtsample == 9999999)
     m_nevtsample = 0;
   m_pedsinADC = ps.getUntrackedParameter<int>("pedsinADC", 0);
   m_hiSaveflag = ps.getUntrackedParameter<int>("hiSaveflag", 0);
   m_pedValflag = ps.getUntrackedParameter<int>("pedValflag", 0);
   if (m_pedValflag < 0)
     m_pedValflag = 0;
   if (m_nevtsample > 0 && m_pedValflag > 0) {
     cout << "WARNING - incompatible cfg options: nevtsample = " << m_nevtsample << ", pedValflag = " << m_pedValflag
          << endl;
     cout << "Setting pedValflag = 0" << endl;
     m_pedValflag = 0;
   }
   if (m_pedValflag > 1)
     m_pedValflag = 1;
   m_startTS = ps.getUntrackedParameter<int>("firstTS", 0);
   if (m_startTS < 0)
     m_startTS = 0;
   m_endTS = ps.getUntrackedParameter<int>("lastTS", 9);
 
   //  m_logFile.open("HcalPedestalAnalysis.log");
 
   hbHists.ALLPEDS = new TH1F("HBHE All Pedestals", "HBHE All Peds", 10, 0, 9);
   hbHists.PEDRMS = new TH1F("HBHE All Pedestal Widths", "HBHE All Pedestal RMS", 100, 0, 3);
   hbHists.PEDMEAN = new TH1F("HBHE All Pedestal Means", "HBHE All Pedestal Means", 100, 0, 9);
   hbHists.CHI2 = new TH1F("HBHE Chi2/ndf for whole range Gauss fit", "HBHE Chi2/ndf Gauss", 200, 0., 50.);
 
   hoHists.ALLPEDS = new TH1F("HO All Pedestals", "HO All Peds", 10, 0, 9);
   hoHists.PEDRMS = new TH1F("HO All Pedestal Widths", "HO All Pedestal RMS", 100, 0, 3);
   hoHists.PEDMEAN = new TH1F("HO All Pedestal Means", "HO All Pedestal Means", 100, 0, 9);
   hoHists.CHI2 = new TH1F("HO Chi2/ndf for whole range Gauss fit", "HO Chi2/ndf Gauss", 200, 0., 50.);
 
   hfHists.ALLPEDS = new TH1F("HF All Pedestals", "HF All Peds", 10, 0, 9);
   hfHists.PEDRMS = new TH1F("HF All Pedestal Widths", "HF All Pedestal RMS", 100, 0, 3);
   hfHists.PEDMEAN = new TH1F("HF All Pedestal Means", "HF All Pedestal Means", 100, 0, 9);
   hfHists.CHI2 = new TH1F("HF Chi2/ndf for whole range Gauss fit", "HF Chi2/ndf Gauss", 200, 0., 50.);
 }

◆ ~HcalPedestalAnalysis()

HcalPedestalAnalysis::~HcalPedestalAnalysis ( )

Destructor.

Definition at line 89 of file HcalPedestalAnalysis.cc.

References _meot, hbHists, hfHists, hoHists, and mps_fire::i.

                                             {
   for (_meot = hbHists.PEDTRENDS.begin(); _meot != hbHists.PEDTRENDS.end(); _meot++) {
     for (int i = 0; i < 16; i++)
       _meot->second[i].first->Delete();
   }
   for (_meot = hoHists.PEDTRENDS.begin(); _meot != hoHists.PEDTRENDS.end(); _meot++) {
     for (int i = 0; i < 16; i++)
       _meot->second[i].first->Delete();
   }
   for (_meot = hfHists.PEDTRENDS.begin(); _meot != hfHists.PEDTRENDS.end(); _meot++) {
     for (int i = 0; i < 16; i++)
       _meot->second[i].first->Delete();
   }
   hbHists.ALLPEDS->Delete();
   hbHists.PEDRMS->Delete();
   hbHists.PEDMEAN->Delete();
   hbHists.CHI2->Delete();
 
   hoHists.ALLPEDS->Delete();
   hoHists.PEDRMS->Delete();
   hoHists.PEDMEAN->Delete();
   hoHists.CHI2->Delete();
 
   hfHists.ALLPEDS->Delete();
   hfHists.PEDRMS->Delete();
   hfHists.PEDMEAN->Delete();
   hfHists.CHI2->Delete();
 }

Member Function Documentation

◆ AllChanHists()

void HcalPedestalAnalysis::AllChanHists	(	const HcalDetId	detid,
		const HcalQIESample &	qie0,
		const HcalQIESample &	qie1,
		const HcalQIESample &	qie2,
		const HcalQIESample &	qie3,
		const HcalQIESample &	qie4,
		const HcalQIESample &	qie5,
		std::map< HcalDetId, std::map< int, PEDBUNCH > > &	toolT
	)

private

Definition at line 392 of file HcalPedestalAnalysis.cc.

References _meot, and HcalQIESample::adc().

Referenced by processEvent().

                                                                                     {
   // this function is due to be called for every channel
 
   _meot = toolT.find(detid);
   map<int, PEDBUNCH> _mei = _meot->second;
   _mei[16].first->Fill(qie4.adc() + qie5.adc() - 1.);
   _mei[17].first->Fill(qie4.adc() + qie5.adc() - qie2.adc() - qie3.adc());
   _mei[18].first->Fill(qie4.adc() + qie5.adc() - (qie0.adc() + qie1.adc() + qie2.adc() + qie3.adc()) / 2.);
 }

◆ done()

int HcalPedestalAnalysis::done	(	const HcalPedestals *	fInputPedestals,
		const HcalPedestalWidths *	fInputWidths,
		HcalPedestals *	fOutputPedestals,
		HcalPedestalWidths *	fOutputWidths
	)

Definition at line 591 of file HcalPedestalAnalysis.cc.

References gather_cfg::cout, evt, fRawPedestals, fRawPedestalWidths, fRefPedestals, fRefPedestalWidths, fTopology, fValPedestals, fValPedestalWidths, hbHists, HcalPedVal(), hfHists, hoHists, mps_fire::i, m_AllPedsOK, m_file, m_nevtsample, m_outputFileROOT, m_pedsinADC, m_pedValflag, m_stat, SampleAnalysis(), and Trendings().

                                                                           {
   int nstat[4];
 
   // Pedestal objects
   // inputs...
   fRefPedestals = fInputPedestals;
   fRefPedestalWidths = fInputPedestalWidths;
 
   // outputs...
   if (m_pedValflag > 0) {
     fValPedestals = fOutputPedestals;
     fValPedestalWidths = fOutputPedestalWidths;
     fRawPedestals = new HcalPedestals(fTopology, m_pedsinADC);
     fRawPedestalWidths = new HcalPedestalWidths(fTopology, m_pedsinADC);
   } else {
     fRawPedestals = fOutputPedestals;
     fRawPedestalWidths = fOutputPedestalWidths;
     fValPedestals = new HcalPedestals(fTopology, m_pedsinADC);
     fValPedestalWidths = new HcalPedestalWidths(fTopology, m_pedsinADC);
   }
 
   // compute pedestal constants
   if (m_nevtsample < 1)
     SampleAnalysis();
   if (m_nevtsample > 0) {
     if (evt % m_nevtsample != 0)
       SampleAnalysis();
   }
 
   // trending histos
   if (m_nevtsample > 0) {
     m_file->cd();
     m_file->cd("HB");
     Trendings(hbHists.PEDTRENDS, hbHists.CHI2, hbHists.CAPID_AVERAGE, hbHists.CAPID_CHI2);
     m_file->cd();
     m_file->cd("HO");
     Trendings(hoHists.PEDTRENDS, hoHists.CHI2, hoHists.CAPID_AVERAGE, hoHists.CAPID_CHI2);
     m_file->cd();
     m_file->cd("HF");
     Trendings(hfHists.PEDTRENDS, hfHists.CHI2, hfHists.CAPID_AVERAGE, hfHists.CAPID_CHI2);
   }
 
   if (m_nevtsample < 1) {
     // pedestal validation: m_AllPedsOK=-1 means not validated,
     //                                   0 everything OK,
     //                                   N>0 : mod(N,100000) drifts + width changes
     //                                         int(N/100000) missing channels
     m_AllPedsOK = -1;
     if (m_pedValflag > 0) {
       for (int i = 0; i < 4; i++)
         nstat[i] = (int)m_stat[i];
       int NPedErrors = HcalPedVal(nstat,
                                   fRefPedestals,
                                   fRefPedestalWidths,
                                   fRawPedestals,
                                   fRawPedestalWidths,
                                   fValPedestals,
                                   fValPedestalWidths);
       m_AllPedsOK = NPedErrors;
     }
     // setting m_AllPedsOK=-2 will inhibit writing pedestals out
     //    if(m_pedValflag==1){
     //      if(evt<100)m_AllPedsOK=-2;
     //    }
   }
 
   // Write other histograms.
   // HB
   m_file->cd();
   m_file->cd("HB");
   hbHists.ALLPEDS->Write();
   hbHists.PEDRMS->Write();
   hbHists.PEDMEAN->Write();
   // HO
   m_file->cd();
   m_file->cd("HO");
   hoHists.ALLPEDS->Write();
   hoHists.PEDRMS->Write();
   hoHists.PEDMEAN->Write();
   // HF
   m_file->cd();
   m_file->cd("HF");
   hfHists.ALLPEDS->Write();
   hfHists.PEDRMS->Write();
   hfHists.PEDMEAN->Write();
 
   m_file->Close();
   cout << "Hcal histograms written to " << m_outputFileROOT.c_str() << endl;
   return (int)m_AllPedsOK;
 }

◆ GetPedConst()

void HcalPedestalAnalysis::GetPedConst	(	std::map< HcalDetId, std::map< int, PEDBUNCH > > &	toolT,
		TH1F *	PedMeans,
		TH1F *	PedWidths
	)

private

Definition at line 433 of file HcalPedestalAnalysis.cc.

References _meot, HcalCondObjectContainer< Item >::addValues(), HLT_2022v15_cff::cap, hltPixelTracks_cff::chi2, fitflag, fRawPedestals, fRawPedestalWidths, mps_fire::i, B2GTnPMonitor_cfi::item, dqmiolumiharvest::j, m_endTS, m_hiSaveflag, m_nevtsample, m_pedsinADC, m_startTS, m_stat, HcalPedestalWidth::setSigma(), and mathSSE::sqrt().

Referenced by SampleAnalysis().

                                                                                                                   {
   // Completely rewritten version oct 2006
   // Compute pedestal constants and fill into HcalPedestals and HcalPedestalWidths objects
   float cap[4];
   float sig[4][4];
   float dcap[4];
   float dsig[4][4];
   float chi2[4];
 
   for (_meot = toolT.begin(); _meot != toolT.end(); _meot++) {
     HcalDetId detid = _meot->first;
 
     // take mean and width from a Gaussian fit or directly from the histo
     if (fitflag > 0) {
       for (int i = 0; i < 4; i++) {
         TF1* fit = _meot->second[i].first->GetFunction("gaus");
         chi2[i] = 0;
         if (fit->GetNDF() != 0)
           chi2[i] = fit->GetChisquare() / fit->GetNDF();
         cap[i] = fit->GetParameter(1);
         sig[i][i] = fit->GetParameter(2);
         dcap[i] = fit->GetParError(1);
         dsig[i][i] = fit->GetParError(2);
       }
     } else {
       for (int i = 0; i < 4; i++) {
         cap[i] = _meot->second[i].first->GetMean();
         sig[i][i] = _meot->second[i].first->GetRMS();
         m_stat[i] = 0;
 
         for (int j = m_startTS; j < m_endTS + 1; j++) {
           m_stat[i] += _meot->second[i].first->GetBinContent(j + 1);
         }
         dcap[i] = sig[i][i] / sqrt(m_stat[i]);
         //        dsig[i][i] = dcap[i]*sig[i][i]/cap[i];
         dsig[i][i] = sig[i][i] / sqrt(2. * m_stat[i]);
         chi2[i] = 0.;
       }
     }
 
     for (int i = 0; i < 4; i++) {
       if (m_hiSaveflag > 0) {
         if (m_pedsinADC)
           _meot->second[i].first->GetXaxis()->SetTitle("ADC");
         else
           _meot->second[i].first->GetXaxis()->SetTitle("Charge, fC");
         _meot->second[i].first->GetYaxis()->SetTitle("CapID samplings");
         _meot->second[i].first->Write();
       }
       if (m_nevtsample > 0) {
         _meot->second[i].second.first[0].push_back(cap[i]);
         _meot->second[i].second.first[1].push_back(dcap[i]);
         _meot->second[i].second.first[2].push_back(sig[i][i]);
         _meot->second[i].second.first[3].push_back(dsig[i][i]);
         _meot->second[i].second.first[4].push_back(chi2[i]);
       }
       PedMeans->Fill(cap[i]);
       PedWidths->Fill(sig[i][i]);
     }
 
     // special histos for Shuichi
     if (m_hiSaveflag == -100) {
       for (int i = 16; i < 19; i++) {
         if (m_pedsinADC)
           _meot->second[i].first->GetXaxis()->SetTitle("ADC");
         else
           _meot->second[i].first->GetXaxis()->SetTitle("Charge, fC");
         _meot->second[i].first->GetYaxis()->SetTitle("Events");
         _meot->second[i].first->Write();
       }
     }
 
     // diagonal sigma is width squared
     sig[0][0] = sig[0][0] * sig[0][0];
     sig[1][1] = sig[1][1] * sig[1][1];
     sig[2][2] = sig[2][2] * sig[2][2];
     sig[3][3] = sig[3][3] * sig[3][3];
 
     // off diagonal sigmas (correlations) are computed from 3 histograms
     // here we still have all 4*3=12 combinations
     sig[0][1] = _meot->second[4].first->GetMean() - cap[0] * cap[1];
     sig[0][2] = _meot->second[8].first->GetMean() - cap[0] * cap[2];
     sig[1][2] = _meot->second[5].first->GetMean() - cap[1] * cap[2];
     sig[1][3] = _meot->second[9].first->GetMean() - cap[1] * cap[3];
     sig[2][3] = _meot->second[6].first->GetMean() - cap[2] * cap[3];
     sig[0][3] = _meot->second[12].first->GetMean() - cap[0] * cap[3];
     sig[1][0] = _meot->second[13].first->GetMean() - cap[1] * cap[0];
     sig[2][0] = _meot->second[10].first->GetMean() - cap[2] * cap[0];
     sig[2][1] = _meot->second[14].first->GetMean() - cap[2] * cap[1];
     sig[3][1] = _meot->second[11].first->GetMean() - cap[3] * cap[1];
     sig[3][2] = _meot->second[15].first->GetMean() - cap[3] * cap[2];
     sig[3][0] = _meot->second[7].first->GetMean() - cap[3] * cap[0];
 
     // there is no proper error calculation for the correlation coefficients
     for (int i = 0; i < 4; i++) {
       if (m_nevtsample > 0) {
         _meot->second[i].second.first[5].push_back(sig[i][(i + 1) % 4]);
         _meot->second[i].second.first[6].push_back(2 * sig[i][i] * dsig[i][i]);
         _meot->second[i].second.first[7].push_back(sig[i][(i + 2) % 4]);
         _meot->second[i].second.first[8].push_back(2 * sig[i][i] * dsig[i][i]);
         _meot->second[i].second.first[9].push_back(sig[i][(i + 3) % 4]);
         _meot->second[i].second.first[10].push_back(2 * sig[i][i] * dsig[i][i]);
       }
       // save product histos if desired
       if (m_hiSaveflag > 10) {
         if (m_pedsinADC)
           _meot->second[i + 4].first->GetXaxis()->SetTitle("ADC^2");
         else
           _meot->second[i + 4].first->GetXaxis()->SetTitle("Charge^2, fC^2");
         _meot->second[i + 4].first->GetYaxis()->SetTitle("2-CapID samplings");
         _meot->second[i + 4].first->Write();
         if (m_pedsinADC)
           _meot->second[i + 8].first->GetXaxis()->SetTitle("ADC^2");
         else
           _meot->second[i + 8].first->GetXaxis()->SetTitle("Charge^2, fC^2");
         _meot->second[i + 8].first->GetYaxis()->SetTitle("2-CapID samplings");
         _meot->second[i + 8].first->Write();
         if (m_pedsinADC)
           _meot->second[i + 12].first->GetXaxis()->SetTitle("ADC^2");
         else
           _meot->second[i + 12].first->GetXaxis()->SetTitle("Charge^2, fC^2");
         _meot->second[i + 12].first->GetYaxis()->SetTitle("2-CapID samplings");
         _meot->second[i + 12].first->Write();
       }
     }
 
     // fill the objects - at this point only close and medium correlations are stored
     // and the matrix is assumed symmetric
     if (m_nevtsample < 1) {
       sig[1][0] = sig[0][1];
       sig[2][0] = sig[0][2];
       sig[2][1] = sig[1][2];
       sig[3][1] = sig[1][3];
       sig[3][2] = sig[2][3];
       sig[0][3] = sig[3][0];
       if (fRawPedestals) {
         HcalPedestal item(detid, cap[0], cap[1], cap[2], cap[3]);
         fRawPedestals->addValues(item);
       }
       if (fRawPedestalWidths) {
         HcalPedestalWidth widthsp(detid);
         widthsp.setSigma(0, 0, sig[0][0]);
         widthsp.setSigma(0, 1, sig[0][1]);
         widthsp.setSigma(0, 2, sig[0][2]);
         widthsp.setSigma(1, 1, sig[1][1]);
         widthsp.setSigma(1, 2, sig[1][2]);
         widthsp.setSigma(1, 3, sig[1][3]);
         widthsp.setSigma(2, 2, sig[2][2]);
         widthsp.setSigma(2, 3, sig[2][3]);
         widthsp.setSigma(3, 3, sig[3][3]);
         widthsp.setSigma(3, 0, sig[0][3]);
         fRawPedestalWidths->addValues(widthsp);
       }
     }
   }
 }

◆ HcalPedVal()

int HcalPedestalAnalysis::HcalPedVal	(	int	nstat[4],
		const HcalPedestals *	fRefPedestals,
		const HcalPedestalWidths *	fRefPedestalWidths,
		HcalPedestals *	fRawPedestals,
		HcalPedestalWidths *	fRawPedestalWidths,
		HcalPedestals *	fValPedestals,
		HcalPedestalWidths *	fValPedestalWidths
	)

static

Definition at line 845 of file HcalPedestalAnalysis.cc.

References HcalCondObjectContainer< Item >::addValues(), DMR_cfg::cerr, fRawPedestals, fRawPedestalWidths, fRefPedestals, fRefPedestalWidths, fValPedestals, fValPedestalWidths, HcalCondObjectContainer< Item >::getAllChannels(), HcalPedestalWidth::getSigma(), HcalPedestal::getValue(), HcalCondObjectContainer< Item >::getValues(), HcalForward, mps_fire::i, createfilelist::int, B2GTnPMonitor_cfi::item, HcalPedestalWidth::setSigma(), mathSSE::sqrt(), and ApeEstimator_cff::width.

Referenced by done().

                                                                              {
   // new version of pedestal validation - it is designed to be as independent of
   // all the rest as possible - you only need to provide valid pedestal objects
   // and a vector of statistics per capID to use this as standalone code
   HcalDetId detid;
   float RefPedVals[4];
   float RefPedSigs[4][4];
   float RawPedVals[4];
   float RawPedSigs[4][4];
   map<HcalDetId, bool> isinRaw;
   map<HcalDetId, bool> isinRef;
   std::vector<DetId> RefChanns = fRefPedestals->getAllChannels();
   std::vector<DetId> RawChanns = fRawPedestals->getAllChannels();
   std::ofstream PedValLog;
   PedValLog.open("HcalPedVal.log");
 
   if (nstat[0] + nstat[1] + nstat[2] + nstat[3] < 2500)
     PedValLog << "HcalPedVal: warning - low statistics" << std::endl;
   // find complete list of channels in current data and reference
   for (int i = 0; i < (int)RawChanns.size(); i++) {
     isinRef[HcalDetId(RawChanns[i])] = false;
   }
   for (int i = 0; i < (int)RefChanns.size(); i++) {
     detid = HcalDetId(RefChanns[i]);
     isinRaw[detid] = false;
     isinRef[detid] = true;
   }
   for (int i = 0; i < (int)RawChanns.size(); i++) {
     detid = HcalDetId(RawChanns[i]);
     isinRaw[detid] = true;
     if (isinRef[detid] == false) {
       PedValLog << "HcalPedVal: channel " << detid << " not found in reference set" << std::endl;
       std::cerr << "HcalPedVal: channel " << detid << " not found in reference set" << std::endl;
     }
   }
 
   // main loop over channels
   int erflag = 0;
   for (int i = 0; i < (int)RefChanns.size(); i++) {
     detid = HcalDetId(RefChanns[i]);
     for (int icap = 0; icap < 4; icap++) {
       RefPedVals[icap] = fRefPedestals->getValues(detid)->getValue(icap);
       for (int icap2 = icap; icap2 < 4; icap2++) {
         RefPedSigs[icap][icap2] = fRefPedestalWidths->getValues(detid)->getSigma(icap, icap2);
         if (icap2 != icap)
           RefPedSigs[icap2][icap] = RefPedSigs[icap][icap2];
       }
     }
 
     // read new raw values
     if (isinRaw[detid]) {
       for (int icap = 0; icap < 4; icap++) {
         RawPedVals[icap] = fRawPedestals->getValues(detid)->getValue(icap);
         for (int icap2 = icap; icap2 < 4; icap2++) {
           RawPedSigs[icap][icap2] = fRawPedestalWidths->getValues(detid)->getSigma(icap, icap2);
           if (icap2 != icap)
             RawPedSigs[icap2][icap] = RawPedSigs[icap][icap2];
         }
       }
 
       // first quick check if raw values make sense: if not, the channel is treated like absent
       for (int icap = 0; icap < 4; icap++) {
         if (RawPedVals[icap] < 1. || RawPedSigs[icap][icap] < 0.01)
           isinRaw[detid] = false;
         for (int icap2 = icap; icap2 < 4; icap2++) {
           if (fabs(RawPedSigs[icap][icap2] / sqrt(RawPedSigs[icap][icap] * RawPedSigs[icap2][icap2])) > 1.)
             isinRaw[detid] = false;
         }
       }
     }
 
     // check raw values against reference
     if (isinRaw[detid]) {
       for (int icap = 0; icap < 4; icap++) {
         int icap2 = (icap + 1) % 4;
         float width = sqrt(RawPedSigs[icap][icap]);
         float erof1 = width / sqrt((float)nstat[icap]);
         float erof2 = sqrt(erof1 * erof1 + RawPedSigs[icap][icap] / (float)nstat[icap]);
         float erofwidth = width / sqrt(2. * nstat[icap]);
         float diffof1 = RawPedVals[icap] - RefPedVals[icap];
         float diffof2 = RawPedVals[icap] + RawPedVals[icap2] - RefPedVals[icap] - RefPedVals[icap2];
         float diffofw = width - sqrt(RefPedSigs[icap][icap]);
 
         // validation in 2 TS for HB, HE, HO, in 1 TS for HF
         int nTS = 2;
         if (detid.subdet() == HcalForward)
           nTS = 1;
         if (nTS == 1 && fabs(diffof1) > 0.5 + erof1) {
           erflag += 1;
           PedValLog << "HcalPedVal: drift in channel " << detid << " cap " << icap << ": " << RawPedVals[icap] << " - "
                     << RefPedVals[icap] << " = " << diffof1 << std::endl;
         }
         if (nTS == 2 && fabs(diffof2) > 0.5 + erof2) {
           erflag += 1;
           PedValLog << "HcalPedVal: drift in channel " << detid << " caps " << icap << "+" << icap2 << ": "
                     << RawPedVals[icap] << "+" << RawPedVals[icap2] << " - " << RefPedVals[icap] << "+"
                     << RefPedVals[icap2] << " = " << diffof2 << std::endl;
         }
         if (fabs(diffofw) > 0.15 * width + erofwidth) {
           erflag += 1;
           PedValLog << "HcalPedVal: width changed in channel " << detid << " cap " << icap << ": " << width << " - "
                     << sqrt(RefPedSigs[icap][icap]) << " = " << diffofw << std::endl;
         }
       }
     }
 
     // for disconnected/bad channels restore reference values
     else {
       PedValLog << "HcalPedVal: no valid data from channel " << detid << std::endl;
       erflag += 100000;
       HcalPedestal item(detid, RefPedVals[0], RefPedVals[1], RefPedVals[2], RefPedVals[3]);
       fValPedestals->addValues(item);
       HcalPedestalWidth widthsp(detid);
       for (int icap = 0; icap < 4; icap++) {
         for (int icap2 = icap; icap2 < 4; icap2++)
           widthsp.setSigma(icap2, icap, RefPedSigs[icap2][icap]);
       }
       fValPedestalWidths->addValues(widthsp);
     }
 
     // end of channel loop
   }
 
   if (erflag == 0)
     PedValLog << "HcalPedVal: all pedestals checked OK" << std::endl;
 
   // now construct the remaining part of the validated objects
   // if nothing changed outside tolerance, validated set = reference set
   if (erflag % 100000 == 0) {
     for (int i = 0; i < (int)RefChanns.size(); i++) {
       detid = HcalDetId(RefChanns[i]);
       if (isinRaw[detid]) {
         HcalPedestalWidth widthsp(detid);
         for (int icap = 0; icap < 4; icap++) {
           RefPedVals[icap] = fRefPedestals->getValues(detid)->getValue(icap);
           for (int icap2 = icap; icap2 < 4; icap2++) {
             RefPedSigs[icap][icap2] = fRefPedestalWidths->getValues(detid)->getSigma(icap, icap2);
             if (icap2 != icap)
               RefPedSigs[icap2][icap] = RefPedSigs[icap][icap2];
             widthsp.setSigma(icap2, icap, RefPedSigs[icap2][icap]);
           }
         }
         fValPedestalWidths->addValues(widthsp);
         HcalPedestal item(detid, RefPedVals[0], RefPedVals[1], RefPedVals[2], RefPedVals[3]);
         fValPedestals->addValues(item);
       }
     }
   }
 
   // if anything changed, validated set = raw set + reference for missing/bad channels
   else {
     for (int i = 0; i < (int)RawChanns.size(); i++) {
       detid = HcalDetId(RawChanns[i]);
       if (isinRaw[detid]) {
         HcalPedestalWidth widthsp(detid);
         for (int icap = 0; icap < 4; icap++) {
           RawPedVals[icap] = fRawPedestals->getValues(detid)->getValue(icap);
           for (int icap2 = icap; icap2 < 4; icap2++) {
             RawPedSigs[icap][icap2] = fRawPedestalWidths->getValues(detid)->getSigma(icap, icap2);
             if (icap2 != icap)
               RawPedSigs[icap2][icap] = RawPedSigs[icap][icap2];
             widthsp.setSigma(icap2, icap, RawPedSigs[icap2][icap]);
           }
         }
         fValPedestalWidths->addValues(widthsp);
         HcalPedestal item(detid, RawPedVals[0], RawPedVals[1], RawPedVals[2], RawPedVals[3]);
         fValPedestals->addValues(item);
       }
     }
   }
   return erflag;
 }

◆ per2CapsHists()

void HcalPedestalAnalysis::per2CapsHists	(	int	flag,
		int	id,
		const HcalDetId	detid,
		const HcalQIESample &	qie1,
		const HcalQIESample &	qie2,
		std::map< HcalDetId, std::map< int, PEDBUNCH > > &	toolT,
		const HcalDbService &	cond
	)

private

Definition at line 243 of file HcalPedestalAnalysis.cc.

References _meot, HcalQIESample::adc(), trigObjTnPSource_cfi::bins, HcalQIESample::capid(), HcalQIECoder::charge(), HcalDetId::depth(), evt, RemoveAddSevLevel::flag, hbHists, hfHists, hoHists, mps_fire::i, HcalDetId::ieta(), edm::eventsetup::heterocontainer::insert(), HcalDetId::iphi(), m_nevtsample, m_pedsinADC, Skims_PA_cff::name, and state.

Referenced by processEvent().

                                                                     {
   // this function is due to be called for every time slice, it fills either a charge
   // histo for a single capID (flag=0) or a product histo for two capIDs (flag>0)
 
   static const int bins = 10;
   static const int bins2 = 100;
   float lo = -0.5;
   float hi = 9.5;
   map<int, PEDBUNCH> _mei;
   static map<HcalDetId, map<int, float> > QieCalibMap;
   string type = "HBHE";
 
   if (id == 0) {
     if (detid.ieta() < 16)
       type = "HB";
     if (detid.ieta() > 16)
       type = "HE";
     if (detid.ieta() == 16) {
       if (detid.depth() < 3)
         type = "HB";
       if (detid.depth() == 3)
         type = "HE";
     }
   } else if (id == 1)
     type = "HO";
   else if (id == 2)
     type = "HF";
 
   _meot = toolT.find(detid);
 
   // if histos for the current channel do not exist, first create them,
   if (_meot == toolT.end()) {
     map<int, PEDBUNCH> insert;
     map<int, float> qiecalib;
     char name[1024];
     for (int i = 0; i < 4; i++) {
       lo = -0.5;
       // fix from Andy: if you convert to fC and then bin in units of 1, you may 'skip' a bin while
       // filling, since the ADCs are quantized
       if (m_pedsinADC)
         hi = 9.5;
       else
         hi = 11.5;
       sprintf(
           name, "%s Pedestal, eta=%d phi=%d d=%d cap=%d", type.c_str(), detid.ieta(), detid.iphi(), detid.depth(), i);
       insert[i].first = new TH1F(name, name, bins, lo, hi);
       sprintf(name,
               "%s Product, eta=%d phi=%d d=%d caps=%d*%d",
               type.c_str(),
               detid.ieta(),
               detid.iphi(),
               detid.depth(),
               i,
               (i + 1) % 4);
       insert[4 + i].first = new TH1F(name, name, bins2, 0., 100.);
       sprintf(name,
               "%s Product, eta=%d phi=%d d=%d caps=%d*%d",
               type.c_str(),
               detid.ieta(),
               detid.iphi(),
               detid.depth(),
               i,
               (i + 2) % 4);
       insert[8 + i].first = new TH1F(name, name, bins2, 0., 100.);
       sprintf(name,
               "%s Product, eta=%d phi=%d d=%d caps=%d*%d",
               type.c_str(),
               detid.ieta(),
               detid.iphi(),
               detid.depth(),
               i,
               (i + 3) % 4);
       insert[12 + i].first = new TH1F(name, name, bins2, 0., 100.);
     }
     sprintf(name, "%s Signal in TS 4+5, eta=%d phi=%d d=%d", type.c_str(), detid.ieta(), detid.iphi(), detid.depth());
     insert[16].first = new TH1F(name, name, 21, -0.5, 20.5);
     sprintf(
         name, "%s Signal in TS 4+5-2-3, eta=%d phi=%d d=%d", type.c_str(), detid.ieta(), detid.iphi(), detid.depth());
     insert[17].first = new TH1F(name, name, 21, -10.5, 10.5);
     sprintf(name,
             "%s Signal in TS 4+5-(0+1+2+3)/2., eta=%d phi=%d d=%d",
             type.c_str(),
             detid.ieta(),
             detid.iphi(),
             detid.depth());
     insert[18].first = new TH1F(name, name, 21, -10.5, 10.5);
     toolT[detid] = insert;
     _meot = toolT.find(detid);
     // store QIE calibrations in a map for later reuse
     QieCalibMap[detid] = qiecalib;
   }
 
   _mei = _meot->second;
 
   const HcalQIECoder* coder = cond.getHcalCoder(detid);
   const HcalQIEShape* shape = cond.getHcalShape(coder);
   float charge1 = coder->charge(*shape, qie1.adc(), qie1.capid());
   float charge2 = coder->charge(*shape, qie2.adc(), qie2.capid());
 
   // fill single capID histo
   if (flag == 0) {
     if (m_nevtsample > 0) {
       if ((evt - 1) % m_nevtsample == 0 && state[qie1.capid()]) {
         state[qie1.capid()] = false;
         _mei[qie1.capid()].first->Reset();
         _mei[qie1.capid() + 4].first->Reset();
         _mei[qie1.capid() + 8].first->Reset();
         _mei[qie1.capid() + 12].first->Reset();
       }
     }
     if (qie1.adc() < bins) {
       if (m_pedsinADC)
         _mei[qie1.capid()].first->Fill(qie1.adc());
       else
         _mei[qie1.capid()].first->Fill(charge1);
     } else if (qie1.adc() >= bins) {
       _mei[qie1.capid()].first->AddBinContent(bins + 1, 1);
     }
   }
 
   // fill 2 capID histo
   if (flag > 0) {
     map<int, float> qiecalib = QieCalibMap[detid];
     //float charge1=(qie1.adc()-qiecalib[qie1.capid()+4])/qiecalib[qie1.capid()];
     //float charge2=(qie2.adc()-qiecalib[qie2.capid()+4])/qiecalib[qie2.capid()];
     if (charge1 * charge2 < bins2) {
       _mei[qie1.capid() + 4 * flag].first->Fill(charge1 * charge2);
     } else {
       _mei[qie1.capid() + 4 * flag].first->Fill(bins2);
     }
   }
 
   if (flag == 0) {
     if (id == 0)
       hbHists.ALLPEDS->Fill(qie1.adc());
     else if (id == 1)
       hoHists.ALLPEDS->Fill(qie1.adc());
     else if (id == 2)
       hfHists.ALLPEDS->Fill(qie1.adc());
   }
 }

◆ processEvent()

void HcalPedestalAnalysis::processEvent	(	const HBHEDigiCollection &	hbhe,
		const HODigiCollection &	ho,
		const HFDigiCollection &	hf,
		const HcalDbService &	cond
	)

Definition at line 131 of file HcalPedestalAnalysis.cc.

References AllChanHists(), evt, evt_curr, RemoveAddSevLevel::flag, photonIsolationHIProducer_cfi::hbhe, hbHists, photonIsolationHIProducer_cfi::hf, hfHists, photonIsolationHIProducer_cfi::ho, hoHists, mps_fire::i, HBHEDataFrame::id(), HODataFrame::id(), HFDataFrame::id(), createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, m_coder, m_endTS, m_nevtsample, m_shape, m_startTS, per2CapsHists(), HFDataFrame::sample(), HBHEDataFrame::sample(), HODataFrame::sample(), sample, SampleAnalysis(), HODataFrame::size(), HBHEDataFrame::size(), HFDataFrame::size(), and state.

                                                                    {
   evt++;
   sample = 1;
   evt_curr = evt;
   if (m_nevtsample > 0) {
     sample = (evt - 1) / m_nevtsample + 1;
     evt_curr = evt % m_nevtsample;
     if (evt_curr == 0)
       evt_curr = m_nevtsample;
   }
 
   // Get data for every CAPID.
   // HBHE
   try {
     if (hbhe.empty())
       throw(int) hbhe.size();
     for (HBHEDigiCollection::const_iterator j = hbhe.begin(); j != hbhe.end(); ++j) {
       const HBHEDataFrame digi = (const HBHEDataFrame)(*j);
       m_coder = cond.getHcalCoder(digi.id());
       m_shape = cond.getHcalShape(m_coder);
       for (int k = 0; k < (int)state.size(); k++)
         state[k] = true;
       // here we loop over pairs of time slices, it is more convenient
       // in order to extract the correlation matrix
       for (int i = m_startTS; i < digi.size() && i <= m_endTS; i++) {
         for (int flag = 0; flag < 4; flag++) {
           if (i + flag < digi.size() && i + flag <= m_endTS) {
             per2CapsHists(flag, 0, digi.id(), digi.sample(i), digi.sample(i + flag), hbHists.PEDTRENDS, cond);
           }
         }
       }
       if (m_startTS == 0 && m_endTS > 4) {
         AllChanHists(digi.id(),
                      digi.sample(0),
                      digi.sample(1),
                      digi.sample(2),
                      digi.sample(3),
                      digi.sample(4),
                      digi.sample(5),
                      hbHists.PEDTRENDS);
       }
     }
   } catch (int i) {
     //    m_logFile<< "Event with " << i<<" HBHE Digis passed." << std::endl;
   }
   // HO
   try {
     if (ho.empty())
       throw(int) ho.size();
     for (HODigiCollection::const_iterator j = ho.begin(); j != ho.end(); ++j) {
       const HODataFrame digi = (const HODataFrame)(*j);
       m_coder = cond.getHcalCoder(digi.id());
       for (int i = m_startTS; i < digi.size() && i <= m_endTS; i++) {
         for (int flag = 0; flag < 4; flag++) {
           if (i + flag < digi.size() && i + flag <= m_endTS) {
             per2CapsHists(flag, 1, digi.id(), digi.sample(i), digi.sample(i + flag), hoHists.PEDTRENDS, cond);
           }
         }
       }
       if (m_startTS == 0 && m_endTS > 4) {
         AllChanHists(digi.id(),
                      digi.sample(0),
                      digi.sample(1),
                      digi.sample(2),
                      digi.sample(3),
                      digi.sample(4),
                      digi.sample(5),
                      hoHists.PEDTRENDS);
       }
     }
   } catch (int i) {
     //    m_logFile << "Event with " << i<<" HO Digis passed." << std::endl;
   }
   // HF
   try {
     if (hf.empty())
       throw(int) hf.size();
     for (HFDigiCollection::const_iterator j = hf.begin(); j != hf.end(); ++j) {
       const HFDataFrame digi = (const HFDataFrame)(*j);
       m_coder = cond.getHcalCoder(digi.id());
       for (int i = m_startTS; i < digi.size() && i <= m_endTS; i++) {
         for (int flag = 0; flag < 4; flag++) {
           if (i + flag < digi.size() && i + flag <= m_endTS) {
             per2CapsHists(flag, 2, digi.id(), digi.sample(i), digi.sample(i + flag), hfHists.PEDTRENDS, cond);
           }
         }
       }
       if (m_startTS == 0 && m_endTS > 4) {
         AllChanHists(digi.id(),
                      digi.sample(0),
                      digi.sample(1),
                      digi.sample(2),
                      digi.sample(3),
                      digi.sample(4),
                      digi.sample(5),
                      hfHists.PEDTRENDS);
       }
     }
   } catch (int i) {
     //    m_logFile << "Event with " << i<<" HF Digis passed." << std::endl;
   }
   // Call the function every m_nevtsample events
   if (m_nevtsample > 0) {
     if (evt % m_nevtsample == 0)
       SampleAnalysis();
   }
 }

◆ SampleAnalysis()

void HcalPedestalAnalysis::SampleAnalysis ( )

Definition at line 410 of file HcalPedestalAnalysis.cc.

References GetPedConst(), hbHists, hfHists, hoHists, m_file, and sample.

Referenced by done(), and processEvent().

                                           {
   // it is called every m_nevtsample events (a sample) and the end of run
   char PedSampleNum[20];
 
   // Compute pedestal constants for each HBHE, HO, HF
   sprintf(PedSampleNum, "HB_Sample%d", sample);
   m_file->cd();
   m_file->mkdir(PedSampleNum);
   m_file->cd(PedSampleNum);
   GetPedConst(hbHists.PEDTRENDS, hbHists.PEDMEAN, hbHists.PEDRMS);
   sprintf(PedSampleNum, "HO_Sample%d", sample);
   m_file->cd();
   m_file->mkdir(PedSampleNum);
   m_file->cd(PedSampleNum);
   GetPedConst(hoHists.PEDTRENDS, hoHists.PEDMEAN, hoHists.PEDRMS);
   sprintf(PedSampleNum, "HF_Sample%d", sample);
   m_file->cd();
   m_file->mkdir(PedSampleNum);
   m_file->cd(PedSampleNum);
   GetPedConst(hfHists.PEDTRENDS, hfHists.PEDMEAN, hfHists.PEDRMS);
 }

◆ setTopology()

void HcalPedestalAnalysis::setTopology ( const HcalTopology * htopo )

inline

Definition at line 42 of file HcalPedestalAnalysis.h.

References fTopology.

42 { fTopology = htopo; }

HcalPedestalAnalysis::fTopology

const HcalTopology * fTopology

Definition: HcalPedestalAnalysis.h:142

◆ setup()

void HcalPedestalAnalysis::setup ( const std::string & m_outputFileROOT )

Definition at line 119 of file HcalPedestalAnalysis.cc.

References m_file, and m_outputFileROOT.

Referenced by o2olib.O2OTool::execute().

                                                                   {
   // open the histogram file, create directories within
   m_file = new TFile(m_outputFileROOT.c_str(), "RECREATE");
   m_file->mkdir("HB");
   m_file->cd();
   m_file->mkdir("HO");
   m_file->cd();
   m_file->mkdir("HF");
   m_file->cd();
 }

◆ Trendings()

void HcalPedestalAnalysis::Trendings	(	std::map< HcalDetId, std::map< int, PEDBUNCH > > &	toolT,
		TH1F *	Chi2,
		TH1F *	CapidAverage,
		TH1F *	CapidChi2
	)

private

Definition at line 686 of file HcalPedestalAnalysis.cc.

References _meot, trigObjTnPSource_cfi::bins, HcalDetId::depth(), dqmMemoryStats::float, mps_fire::i, HcalDetId::ieta(), HcalDetId::iphi(), dqmiolumiharvest::j, m_nevtsample, Skims_PA_cff::name, sample, and findQualityFiles::size.

Referenced by done().

                                                       {
   // check stability of pedestal constants in a single long run
 
   map<int, std::vector<double> > AverageValues;
 
   for (_meot = toolT.begin(); _meot != toolT.end(); _meot++) {
     for (int i = 0; i < 4; i++) {
       char name[1024];
       HcalDetId detid = _meot->first;
       sprintf(name, "Pedestal trend, eta=%d phi=%d d=%d cap=%d", detid.ieta(), detid.iphi(), detid.depth(), i);
       int bins = _meot->second[i].second.first[0].size();
       float lo = 0.5;
       float hi = (float)bins + 0.5;
       _meot->second[i].second.second.push_back(new TH1F(name, name, bins, lo, hi));
       sprintf(name, "Width trend, eta=%d phi=%d d=%d cap=%d", detid.ieta(), detid.iphi(), detid.depth(), i);
       bins = _meot->second[i].second.first[2].size();
       hi = (float)bins + 0.5;
       _meot->second[i].second.second.push_back(new TH1F(name, name, bins, lo, hi));
       sprintf(name,
               "Correlation trend, eta=%d phi=%d d=%d caps=%d*%d",
               detid.ieta(),
               detid.iphi(),
               detid.depth(),
               i,
               (i + 1) % 4);
       bins = _meot->second[i].second.first[5].size();
       hi = (float)bins + 0.5;
       _meot->second[i].second.second.push_back(new TH1F(name, name, bins, lo, hi));
       /*      sprintf(name,"Correlation trend, eta=%d phi=%d d=%d caps=%d*%d",detid.ieta(),detid.iphi(),detid.depth(),i,(i+2)%4);
       bins = _meot->second[i].second.first[7].size();
       hi = (float)bins+0.5;
       _meot->second[i].second.second.push_back(new TH1F(name,name,bins,lo,hi));
       sprintf(name,"Correlation trend, eta=%d phi=%d d=%d caps=%d*%d",detid.ieta(),detid.iphi(),detid.depth(),i,(i+3)%4);
       bins = _meot->second[i].second.first[9].size();
       hi = (float)bins+0.5;
       _meot->second[i].second.second.push_back(new TH1F(name,name,bins,lo,hi)); */
 
       std::vector<double>::iterator sample_it;
       // Pedestal mean - put content and errors
       int j = 0;
       for (sample_it = _meot->second[i].second.first[0].begin(); sample_it != _meot->second[i].second.first[0].end();
            ++sample_it) {
         _meot->second[i].second.second[0]->SetBinContent(++j, *sample_it);
       }
       j = 0;
       for (sample_it = _meot->second[i].second.first[1].begin(); sample_it != _meot->second[i].second.first[1].end();
            ++sample_it) {
         _meot->second[i].second.second[0]->SetBinError(++j, *sample_it);
       }
       // fit with a constant - extract parameters
       _meot->second[i].second.second[0]->Fit("pol0", "Q");
       TF1* fit = _meot->second[i].second.second[0]->GetFunction("pol0");
       AverageValues[0].push_back(fit->GetParameter(0));
       AverageValues[1].push_back(fit->GetParError(0));
       if (sample > 1)
         AverageValues[2].push_back(fit->GetChisquare() / fit->GetNDF());
       else
         AverageValues[2].push_back(fit->GetChisquare());
       sprintf(name, "Sample (%d events)", m_nevtsample);
       _meot->second[i].second.second[0]->GetXaxis()->SetTitle(name);
       _meot->second[i].second.second[0]->GetYaxis()->SetTitle("Pedestal value");
       _meot->second[i].second.second[0]->Write();
       // Pedestal width - put content and errors
       j = 0;
       for (sample_it = _meot->second[i].second.first[2].begin(); sample_it != _meot->second[i].second.first[2].end();
            ++sample_it) {
         _meot->second[i].second.second[1]->SetBinContent(++j, *sample_it);
       }
       j = 0;
       for (sample_it = _meot->second[i].second.first[3].begin(); sample_it != _meot->second[i].second.first[3].end();
            ++sample_it) {
         _meot->second[i].second.second[1]->SetBinError(++j, *sample_it);
       }
       _meot->second[i].second.second[1]->GetXaxis()->SetTitle(name);
       _meot->second[i].second.second[1]->GetYaxis()->SetTitle("Pedestal width");
       _meot->second[i].second.second[1]->Write();
       // Correlation coeffs - put contents and errors
       j = 0;
       for (sample_it = _meot->second[i].second.first[5].begin(); sample_it != _meot->second[i].second.first[5].end();
            ++sample_it) {
         _meot->second[i].second.second[2]->SetBinContent(++j, *sample_it);
       }
       j = 0;
       for (sample_it = _meot->second[i].second.first[6].begin(); sample_it != _meot->second[i].second.first[6].end();
            ++sample_it) {
         _meot->second[i].second.second[2]->SetBinError(++j, *sample_it);
       }
       _meot->second[i].second.second[2]->GetXaxis()->SetTitle(name);
       _meot->second[i].second.second[2]->GetYaxis()->SetTitle("Close correlation");
       _meot->second[i].second.second[2]->Write();
       /*     j=0;
       for(sample_it=_meot->second[i].second.first[7].begin();
           sample_it!=_meot->second[i].second.first[7].end();sample_it++){
         _meot->second[i].second.second[3]->SetBinContent(++j,*sample_it);
       }
       j=0;
       for(sample_it=_meot->second[i].second.first[8].begin();
           sample_it!=_meot->second[i].second.first[8].end();sample_it++){
         _meot->second[i].second.second[3]->SetBinError(++j,*sample_it);
       }
       _meot->second[i].second.second[3]->GetXaxis()->SetTitle(name);
       _meot->second[i].second.second[3]->GetYaxis()->SetTitle("Intermediate correlation");
       _meot->second[i].second.second[3]->Write();
       j=0;
       for(sample_it=_meot->second[i].second.first[9].begin();
           sample_it!=_meot->second[i].second.first[9].end();sample_it++){
         _meot->second[i].second.second[4]->SetBinContent(++j,*sample_it);
       }
       j=0;
       for(sample_it=_meot->second[i].second.first[10].begin();
           sample_it!=_meot->second[i].second.first[10].end();sample_it++){
         _meot->second[i].second.second[4]->SetBinError(++j,*sample_it);
       }
       _meot->second[i].second.second[4]->GetXaxis()->SetTitle(name);
       _meot->second[i].second.second[4]->GetYaxis()->SetTitle("Distant correlation");
       _meot->second[i].second.second[4]->Write(); */
       // chi2
       for (sample_it = _meot->second[i].second.first[4].begin(); sample_it != _meot->second[i].second.first[4].end();
            ++sample_it) {
         Chi2->Fill(*sample_it);
       }
     }
   }
   CapidAverage = new TH1F("Constant fit: Pedestal Values",
                           "Constant fit: Pedestal Values",
                           AverageValues[0].size(),
                           0.,
                           AverageValues[0].size());
   std::vector<double>::iterator sample_it;
   int j = 0;
   for (sample_it = AverageValues[0].begin(); sample_it != AverageValues[0].end(); ++sample_it) {
     CapidAverage->SetBinContent(++j, *sample_it);
   }
   j = 0;
   for (sample_it = AverageValues[1].begin(); sample_it != AverageValues[1].end(); ++sample_it) {
     CapidAverage->SetBinError(++j, *sample_it);
   }
   CapidChi2 = new TH1F(
       "Constant fit: Chi2/ndf", "Constant fit: Chi2/ndf", AverageValues[2].size(), 0., AverageValues[2].size());
   j = 0;
   for (sample_it = AverageValues[2].begin(); sample_it != AverageValues[2].end(); ++sample_it) {
     CapidChi2->SetBinContent(++j, *sample_it);
     //CapidChi2->SetBinError(++j,0);
   }
   Chi2->GetXaxis()->SetTitle("Chi2/ndf");
   Chi2->GetYaxis()->SetTitle("50 x [(16+2) x 4 x 4] `events`");
   Chi2->Write();
   CapidAverage->GetYaxis()->SetTitle("Pedestal value");
   CapidAverage->GetXaxis()->SetTitle("(16+2) x 4 x 4 `events`");
   CapidAverage->Write();
   CapidChi2->GetYaxis()->SetTitle("Chi2/ndf");
   CapidChi2->GetXaxis()->SetTitle("(16+2) x 4 x 4 `events`");
   CapidChi2->Write();
 }

Member Data Documentation

◆ _meot

std::map<HcalDetId, std::map<int, PEDBUNCH> >::iterator HcalPedestalAnalysis::_meot

private

Definition at line 130 of file HcalPedestalAnalysis.h.

Referenced by AllChanHists(), GetPedConst(), per2CapsHists(), Trendings(), and ~HcalPedestalAnalysis().

◆ ALLPEDS

TH1F* HcalPedestalAnalysis::ALLPEDS

Definition at line 123 of file HcalPedestalAnalysis.h.

◆ CAPID_AVERAGE

TH1F* HcalPedestalAnalysis::CAPID_AVERAGE

Definition at line 127 of file HcalPedestalAnalysis.h.

◆ CAPID_CHI2

TH1F* HcalPedestalAnalysis::CAPID_CHI2

Definition at line 128 of file HcalPedestalAnalysis.h.

◆ CHI2

TH1F* HcalPedestalAnalysis::CHI2

Definition at line 126 of file HcalPedestalAnalysis.h.

◆ evt

int HcalPedestalAnalysis::evt

private

Definition at line 137 of file HcalPedestalAnalysis.h.

Referenced by done(), HcalPedestalAnalysis(), per2CapsHists(), and processEvent().

◆ evt_curr

int HcalPedestalAnalysis::evt_curr

private

Definition at line 139 of file HcalPedestalAnalysis.h.

Referenced by processEvent().

◆ fitflag

const int HcalPedestalAnalysis::fitflag = 0

staticprivate

Definition at line 144 of file HcalPedestalAnalysis.h.

Referenced by GetPedConst().

◆ fRawPedestals

HcalPedestals* HcalPedestalAnalysis::fRawPedestals

private

Definition at line 133 of file HcalPedestalAnalysis.h.

Referenced by done(), GetPedConst(), and HcalPedVal().

◆ fRawPedestalWidths

HcalPedestalWidths* HcalPedestalAnalysis::fRawPedestalWidths

private

Definition at line 134 of file HcalPedestalAnalysis.h.

Referenced by done(), GetPedConst(), and HcalPedVal().

◆ fRefPedestals

const HcalPedestals* HcalPedestalAnalysis::fRefPedestals

private

Definition at line 131 of file HcalPedestalAnalysis.h.

Referenced by done(), and HcalPedVal().

◆ fRefPedestalWidths

const HcalPedestalWidths* HcalPedestalAnalysis::fRefPedestalWidths

private

Definition at line 132 of file HcalPedestalAnalysis.h.

Referenced by done(), and HcalPedVal().

◆ fTopology

const HcalTopology* HcalPedestalAnalysis::fTopology

private

Definition at line 142 of file HcalPedestalAnalysis.h.

Referenced by done(), and setTopology().

◆ fValPedestals

HcalPedestals* HcalPedestalAnalysis::fValPedestals

private

Definition at line 135 of file HcalPedestalAnalysis.h.

Referenced by done(), and HcalPedVal().

◆ fValPedestalWidths

HcalPedestalWidths* HcalPedestalAnalysis::fValPedestalWidths

private

Definition at line 136 of file HcalPedestalAnalysis.h.

Referenced by done(), and HcalPedVal().

◆ hbHists

struct { ... } HcalPedestalAnalysis::hbHists

Referenced by done(), HcalPedestalAnalysis(), per2CapsHists(), processEvent(), SampleAnalysis(), and ~HcalPedestalAnalysis().

◆ hfHists

struct { ... } HcalPedestalAnalysis::hfHists

Referenced by done(), HcalPedestalAnalysis(), per2CapsHists(), processEvent(), SampleAnalysis(), and ~HcalPedestalAnalysis().

◆ hoHists

struct { ... } HcalPedestalAnalysis::hoHists

Referenced by done(), HcalPedestalAnalysis(), per2CapsHists(), processEvent(), SampleAnalysis(), and ~HcalPedestalAnalysis().

◆ m_AllPedsOK

int HcalPedestalAnalysis::m_AllPedsOK

private

Definition at line 117 of file HcalPedestalAnalysis.h.

Referenced by done(), and HcalPedestalAnalysis().

◆ m_coder

const HcalQIECoder* HcalPedestalAnalysis::m_coder

private

Definition at line 120 of file HcalPedestalAnalysis.h.

Referenced by HcalPedestalAnalysis(), and processEvent().

◆ m_endTS

int HcalPedestalAnalysis::m_endTS

private

Definition at line 112 of file HcalPedestalAnalysis.h.

Referenced by GetPedConst(), HcalPedestalAnalysis(), and processEvent().

◆ m_file

TFile* HcalPedestalAnalysis::m_file

private

Definition at line 105 of file HcalPedestalAnalysis.h.

Referenced by done(), HcalPedestalAnalysis(), SampleAnalysis(), and setup().

◆ m_hiSaveflag

int HcalPedestalAnalysis::m_hiSaveflag

private

Definition at line 115 of file HcalPedestalAnalysis.h.

Referenced by GetPedConst(), and HcalPedestalAnalysis().

◆ m_logFile

std::ofstream HcalPedestalAnalysis::m_logFile

private

Definition at line 110 of file HcalPedestalAnalysis.h.

◆ m_nevtsample

int HcalPedestalAnalysis::m_nevtsample

private

Definition at line 113 of file HcalPedestalAnalysis.h.

Referenced by done(), GetPedConst(), HcalPedestalAnalysis(), per2CapsHists(), processEvent(), and Trendings().

◆ m_outputFileMean

std::string HcalPedestalAnalysis::m_outputFileMean

private

Definition at line 108 of file HcalPedestalAnalysis.h.

Referenced by HcalPedestalAnalysis().

◆ m_outputFileROOT

std::string HcalPedestalAnalysis::m_outputFileROOT

private

Definition at line 107 of file HcalPedestalAnalysis.h.

Referenced by done(), HcalPedestalAnalysis(), and setup().

◆ m_outputFileWidth

std::string HcalPedestalAnalysis::m_outputFileWidth

private

Definition at line 109 of file HcalPedestalAnalysis.h.

Referenced by HcalPedestalAnalysis().

◆ m_pedsinADC

int HcalPedestalAnalysis::m_pedsinADC

private

Definition at line 114 of file HcalPedestalAnalysis.h.

Referenced by done(), GetPedConst(), HcalPedestalAnalysis(), and per2CapsHists().

◆ m_pedValflag

int HcalPedestalAnalysis::m_pedValflag

private

Definition at line 116 of file HcalPedestalAnalysis.h.

Referenced by done(), and HcalPedestalAnalysis().

◆ m_shape

const HcalQIEShape* HcalPedestalAnalysis::m_shape

private

Definition at line 119 of file HcalPedestalAnalysis.h.

Referenced by HcalPedestalAnalysis(), and processEvent().

◆ m_startTS

int HcalPedestalAnalysis::m_startTS

private

Definition at line 111 of file HcalPedestalAnalysis.h.

Referenced by GetPedConst(), HcalPedestalAnalysis(), and processEvent().

◆ m_stat

float HcalPedestalAnalysis::m_stat[4]

private

Definition at line 140 of file HcalPedestalAnalysis.h.

Referenced by done(), GetPedConst(), and HcalPedestalAnalysis().

◆ PEDMEAN

TH1F* HcalPedestalAnalysis::PEDMEAN

Definition at line 125 of file HcalPedestalAnalysis.h.

◆ PEDRMS

TH1F* HcalPedestalAnalysis::PEDRMS

Definition at line 124 of file HcalPedestalAnalysis.h.

◆ PEDTRENDS

std::map<HcalDetId, std::map<int, PEDBUNCH> > HcalPedestalAnalysis::PEDTRENDS

Definition at line 122 of file HcalPedestalAnalysis.h.

◆ sample

int HcalPedestalAnalysis::sample

private

Definition at line 138 of file HcalPedestalAnalysis.h.

Referenced by HcalPedestalAnalysis(), svgfig.Curve::Path(), processEvent(), SampleAnalysis(), and Trendings().

◆ state

std::vector<bool> HcalPedestalAnalysis::state

private

Definition at line 141 of file HcalPedestalAnalysis.h.

Referenced by crabFunctions.CrabTask::handleNoState(), HcalPedestalAnalysis(), per2CapsHists(), processEvent(), and crabFunctions.CrabTask::update().

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ PEDBUNCH

Constructor & Destructor Documentation

◆ HcalPedestalAnalysis()

◆ ~HcalPedestalAnalysis()

Member Function Documentation

◆ AllChanHists()

◆ done()

◆ GetPedConst()

◆ HcalPedVal()

◆ per2CapsHists()

◆ processEvent()

◆ SampleAnalysis()

◆ setTopology()

◆ setup()

◆ Trendings()

Member Data Documentation

◆ _meot

◆ ALLPEDS

◆ CAPID_AVERAGE

◆ CAPID_CHI2

◆ CHI2

◆ evt

◆ evt_curr

◆ fitflag

◆ fRawPedestals

◆ fRawPedestalWidths

◆ fRefPedestals

◆ fRefPedestalWidths

◆ fTopology

◆ fValPedestals

◆ fValPedestalWidths

◆ hbHists

◆ hfHists

◆ hoHists

◆ m_AllPedsOK

◆ m_coder

◆ m_endTS

◆ m_file

◆ m_hiSaveflag

◆ m_logFile

◆ m_nevtsample

◆ m_outputFileMean

◆ m_outputFileROOT

◆ m_outputFileWidth

◆ m_pedsinADC

◆ m_pedValflag

◆ m_shape

◆ m_startTS

◆ m_stat

◆ PEDMEAN

◆ PEDRMS

◆ PEDTRENDS

◆ sample

◆ state