#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 25 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, 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");
 
   // 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 74 of file HcalTB06Histo.cc.

74 {}

Member Function Documentation

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

Definition at line 93 of file HcalTB06Histo.cc.

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

Referenced by HcalTB06Analysis::analyze().

                                                                        { 
 
   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 80 of file HcalTB06Histo.cc.

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

Referenced by HcalTB06Analysis::analyze().

                                                                      {
 
   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 109 of file HcalTB06Histo.cc.

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

Referenced by HcalTB06Analysis::analyze().

                                            {
 
   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_;
   }
 }