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

Public Member Functions
void	fillEdep (double etots, double eecals, double ehcals)

void	fillPrimary (double energy, double eta, double phi)

void	fillTree (std::vector< double > &ecalo, std::vector< double > &etrig)

	HcalTB06Histo (const edm::ParameterSet &ps)

virtual	~HcalTB06Histo ()

Private Attributes
double	eBeam_

TH1D *	edecN

TH1D *	edecS

TH2D *	edehS

double	edepEC_

double	edepHB_

double	edepHO_

TH1D *	edepN

TH1D *	edepS

double	edepS1_

double	edepS2_

double	edepS3_

double	edepS4_

double	edepS7_

double	edepS8_

double	edepVC_

TH1D *	edhcN

TH1D *	edhcS

TH1D *	emhcN

double	etaBeam_

TH1D *	iEta

TH1D *	iniE

TH1D *	iPhi

double	mip_

bool	mkTree_

double	noiseEC_

double	noiseHB_

double	noiseHO_

double	phiBeam_

TTree *	tree_

bool	verbose_

Detailed Description

Description: Histogram handling for Hcal Test Beam 2006 studies

Usage: Sets up histograms and stores in a file

Definition at line 31 of file HcalTB06Histo.h.

Constructor & Destructor Documentation

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

Definition at line 26 of file HcalTB06Histo.cc.

References eBeam_, edecN, edecS, edehS, edepEC_, edepHB_, edepHO_, edepN, edepS, edepS1_, edepS2_, edepS3_, edepS4_, edepS7_, edepS8_, edepVC_, edhcN, edhcS, emhcN, etaBeam_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), iEta, iniE, iPhi, edm::Service< T >::isAvailable(), TFileService::make(), mip_, mkTree_, noiseEC_, noiseHB_, noiseHO_, phiBeam_, compare::tfile, tree_, and verbose_.

                                                       {
   verbose_ = ps.getUntrackedParameter<bool>("Verbose", false);
   double em1 = ps.getUntrackedParameter<double>("ETtotMax", 400.);
   double em2 = ps.getUntrackedParameter<double>("EHCalMax", 4.0);
   mkTree_ = ps.getUntrackedParameter<bool>("MakeTree", false);
   eBeam_ = 50.;
   mip_ = ps.getParameter<double>("MIP");
   edm::LogInfo("HcalTBSim") << "Verbose :" << verbose_ << " MakeTree: " << mkTree_ << " EMax: " << em1 << ":" << em2
                             << " MIP " << mip_;
 
   // Book histograms
   edm::Service<TFileService> tfile;
 
   if (!tfile.isAvailable())
     throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                       << "please add it to config file";
   iniE = tfile->make<TH1D>("iniE", "Incident Energy (GeV)", 4000, 0., em1);
   iEta = tfile->make<TH1D>("iEta", "Eta at incidence     ", 300, 0., 3.);
   iPhi = tfile->make<TH1D>("iPhi", "Phi at incidence     ", 300, -1., 1.);
   edepS = tfile->make<TH1D>("edepS", "Energy deposit == Total", 4000, 0., em1);
   edecS = tfile->make<TH1D>("edecS", "Energy deposit == ECal ", 300, -2., 28.);
   edhcS = tfile->make<TH1D>("edhcS", "Energy deposit == HCal ", 4000, 0., em2);
   edepN = tfile->make<TH1D>("edepN", "Etot/Ebeam   ", 200, -2.5, 2.5);
   edecN = tfile->make<TH1D>("edecN", "Eecal/Ebeam  ", 200, -2.5, 2.5);
   edhcN = tfile->make<TH1D>("edhcN", "Ehcal/Ebeam  ", 200, -2.5, 2.5);
   emhcN = tfile->make<TH1D>("emhcN", "Ehcal/Ebeam MIP in Ecal", 200, -2.5, 2.5);
   edehS = tfile->make<TH2D>("edehS", "Hcal vs Ecal", 100, 0., em1, 100, 0., em2);
 
   if (mkTree_) {
     tree_ = tfile->make<TTree>("TB06Sim", "TB06Sim");
     tree_->Branch("eBeam_", &eBeam_, "eBeam_/D");
     tree_->Branch("etaBeam_", &etaBeam_, "etaBeam_/D");
     tree_->Branch("phiBeam_", &phiBeam_, "phiBeam_/D");
     tree_->Branch("edepEC_", &edepEC_, "edepEC_/D");
     tree_->Branch("edepHB_", &edepHB_, "edepHB_/D");
     tree_->Branch("edepHO_", &edepHO_, "edepHO_/D");
     tree_->Branch("noiseEC_", &noiseEC_, "noiseEC_/D");
     tree_->Branch("noiseHB_", &noiseHB_, "noiseHB_/D");
     tree_->Branch("noiseHO_", &noiseHO_, "noiseHO_/D");
     tree_->Branch("edepS1_", &edepS1_, "edepS1_/D");
     tree_->Branch("edepS2_", &edepS2_, "edepS2_/D");
     tree_->Branch("edepS3_", &edepS3_, "edepS3_/D");
     tree_->Branch("edepS4_", &edepS4_, "edepS4_/D");
     tree_->Branch("edepVC_", &edepVC_, "edepVC_/D");
     tree_->Branch("edepS7_", &edepS7_, "edepS7_/D");
     tree_->Branch("edepS8_", &edepS8_, "edepS8_/D");
   }
 }

HcalTB06Histo::~HcalTB06Histo ( )

virtual

Definition at line 75 of file HcalTB06Histo.cc.

75 {}

Member Function Documentation

void HcalTB06Histo::fillEdep	(	double	etots,
		double	eecals,
		double	ehcals
	)

Definition at line 92 of file HcalTB06Histo.cc.

References eBeam_, edecN, edecS, edehS, edepN, edepS, edhcN, edhcS, emhcN, mip_, and verbose_.

                                                                        {
   if (verbose_)
     edm::LogInfo("HcalTBSim") << "HcalTB06Histo:::fillEdep: Simulated Total " << etots << " ECal " << eecals << " HCal "
                               << ehcals;
   edepS->Fill(etots);
   edecS->Fill(eecals);
   edhcS->Fill(ehcals);
   edepN->Fill(etots / eBeam_);
   edecN->Fill(eecals / eBeam_);
   edhcN->Fill(ehcals / eBeam_);
   if (eecals <= mip_) {
     emhcN->Fill(etots / eBeam_);
   }
   edehS->Fill(eecals, ehcals);
 }

void HcalTB06Histo::fillPrimary	(	double	energy,
		double	eta,
		double	phi
	)

Definition at line 81 of file HcalTB06Histo.cc.

References eBeam_, relval_parameters_module::energy, PVValHelper::eta, etaBeam_, iEta, iniE, iPhi, phi, phiBeam_, and verbose_.

                                                                      {
   if (verbose_)
     edm::LogInfo("HcalTBSim") << "HcalTB06Histo::fillPrimary: Energy " << energy << " Eta " << eta << " Phi " << phi;
   eBeam_ = energy;
   etaBeam_ = eta;
   phiBeam_ = phi;
   iniE->Fill(energy);
   iEta->Fill(eta);
   iPhi->Fill(phi);
 }

void HcalTB06Histo::fillTree	(	std::vector< double > &	ecalo,
		std::vector< double > &	etrig
	)

Definition at line 108 of file HcalTB06Histo.cc.

References edepEC_, edepHB_, edepHO_, edepS1_, edepS2_, edepS3_, edepS4_, edepS7_, edepS8_, edepVC_, mkTree_, noiseEC_, noiseHB_, noiseHO_, tree_, and verbose_.

                                                                                {
   if (mkTree_) {
     edepEC_ = ecalo[0];
     noiseEC_ = ecalo[1];
     edepHB_ = ecalo[2];
     noiseHB_ = ecalo[3];
     edepHO_ = ecalo[4];
     noiseHO_ = ecalo[5];
     edepS1_ = etrig[0];
     edepS2_ = etrig[1];
     edepS3_ = etrig[2];
     edepS4_ = etrig[3];
     edepVC_ = etrig[4];
     edepS7_ = etrig[5];
     edepS8_ = etrig[6];
     tree_->Fill();
     if (verbose_)
       edm::LogInfo("HcalTBSim") << "HcalTB06Histo:::fillTree: Energies " << edepEC_ << ":" << noiseEC_ << ":" << edepHB_
                                 << ":" << noiseHB_ << ":" << edepHO_ << ":" << noiseHO_ << " Trigger counters "
                                 << edepS1_ << ":" << edepS2_ << ":" << edepS3_ << ":" << edepS4_ << ":" << edepVC_
                                 << ":" << edepS7_ << ":" << edepS8_;
   }
 }