#include <SimG4CMS/HcalTestBeam/interface/HcalTB02Histo.h>

Public Member Functions
void	fillAllTime (float v)

void	fillProfile (int ilayer, float value)

void	fillTransProf (float u, float v)

float	getMean (int ilayer)

float	getRMS (int ilayer)

	HcalTB02Histo (const edm::ParameterSet &ps)

virtual	~HcalTB02Histo ()

Private Attributes
std::string	fileName

std::vector< TH1D * >	rt_histoProf

TH1D *	rt_tbTimes

TH2D *	rt_TransProf

bool	verbose

Detailed Description

Description: Histogram handling for Hcal Test Beam 2002 studies

Usage: Sets up histograms and stores in a file

Definition at line 30 of file HcalTB02Histo.h.

Constructor & Destructor Documentation

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

Definition at line 26 of file HcalTB02Histo.cc.

References fileName, edm::ParameterSet::getUntrackedParameter(), h, edm::Service< T >::isAvailable(), TFileService::make(), rt_histoProf, rt_tbTimes, rt_TransProf, compare::tfile, runGCPTkAlMap::title, and verbose.

                                                       : rt_tbTimes(nullptr), rt_TransProf(nullptr) {
   verbose = ps.getUntrackedParameter<bool>("Verbose", false);
   edm::LogVerbatim("HcalTBSim") << "HcalTB02Histo:: initialised with o/p file " << fileName << " verbosity " << verbose;
 
   // Book histograms
   edm::Service<TFileService> tfile;
 
   if (!tfile.isAvailable())
     throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                       << "please add it to config file";
 
   char title[80];
   for (int ilayer = 0; ilayer < 19; ilayer++) {
     sprintf(title, "Scint. Energy in Layer %d ", ilayer);
     TH1D* h = tfile->make<TH1D>(title, title, 500, 0., 1.5);
     rt_histoProf.push_back(h);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalTBSim") << "HcalTB02Histo:: Initialise Histo " << title;
 #endif
   }
   sprintf(title, "All Hit Time slices");
   rt_tbTimes = tfile->make<TH1D>(title, title, 100, 0., 100.);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB02Histo:: Initialise Histo " << title;
 #endif
   sprintf(title, "Transv. Shower Profile");
   rt_TransProf = tfile->make<TH2D>(title, title, 100, 0., 1., 1000, 0., 0.5);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB02Histo:: Initialise Histo " << title;
 #endif
 }

HcalTB02Histo::~HcalTB02Histo ( )

virtual

Definition at line 58 of file HcalTB02Histo.cc.

                               {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << " Deleting HcalTB02Histo";
 #endif
 }

Member Function Documentation

void HcalTB02Histo::fillAllTime ( float v )

Definition at line 68 of file HcalTB02Histo.cc.

References LogDebug, rt_tbTimes, and findQualityFiles::v.

                                        {
   LogDebug("HcalTBSim") << "HcalTB02Histo::Fill Time histo with " << v;
   if (rt_tbTimes) {
     rt_tbTimes->Fill(v);
   }
 }

void HcalTB02Histo::fillProfile	(	int	ilayer,
		float	value
	)

Definition at line 83 of file HcalTB02Histo.cc.

References LogDebug, rt_histoProf, and relativeConstraints::value.

                                                        {
   if (ilayer < (int)(rt_histoProf.size())) {
     rt_histoProf[ilayer]->Fill(value);
     LogDebug("HcalTBSim") << "HcalTB02Histo::Fill profile " << ilayer << " with " << value;
   }
 }

void HcalTB02Histo::fillTransProf	(	float	u,
		float	v
	)

Definition at line 75 of file HcalTB02Histo.cc.

References LogDebug, rt_TransProf, and findQualityFiles::v.

                                                   {
   LogDebug("HcalTBSim") << "HcalTB02Histo:::Fill Shower Transv. Profile histo"
                         << " with " << u << " and " << v;
   if (rt_TransProf) {
     rt_TransProf->Fill(u, v);
   }
 }

float HcalTB02Histo::getMean ( int ilayer )

Definition at line 90 of file HcalTB02Histo.cc.

References rt_histoProf.

                                        {
   if (ilayer < (int)(rt_histoProf.size())) {
     return rt_histoProf[ilayer]->GetMean();
   } else {
     return 0;
   }
 }

float HcalTB02Histo::getRMS ( int ilayer )

Definition at line 98 of file HcalTB02Histo.cc.

References rt_histoProf.

                                       {
   if (ilayer < (int)(rt_histoProf.size())) {
     return rt_histoProf[ilayer]->GetRMS();
   } else {
     return 0;
   }
 }