#include <Histograms.h>

Public Member Functions
void	Fill (double x, double y, double z, double dx, double dy, double dz, double errx, double erry, double errz, double eta, double phi)
	HResolution1DRecHit (TString name_, TFile *file)
	HResolution1DRecHit (std::string name)
	~HResolution1DRecHit ()
Public Attributes
MonitorElement *	hResX
MonitorElement *	hResXVsEta
MonitorElement *	hResY
MonitorElement *	hResYVsEta
MonitorElement *	hResZ
MonitorElement *	hResZVsEta
MonitorElement *	hXPull
MonitorElement *	hXPullVsPos
MonitorElement *	hXResVsPhi
MonitorElement *	hXResVsPos
MonitorElement *	hYPull
MonitorElement *	hYPullVsPos
MonitorElement *	hYResVsPhi
MonitorElement *	hYResVsPos
MonitorElement *	hZPull
MonitorElement *	hZPullVsPos
MonitorElement *	hZResVsPhi
MonitorElement *	hZResVsPos
std::string	theName
Private Attributes
DQMStore *	dbe_

Detailed Description

Definition at line 297 of file Histograms.h.

Constructor & Destructor Documentation

HResolution1DRecHit::HResolution1DRecHit ( std::string name ) [inline]

Definition at line 299 of file Histograms.h.

References DQMStore::book1D(), DQMStore::book2D(), dbe_, hResX, hResXVsEta, hResY, hResYVsEta, hResZ, hResZVsEta, hXPull, hXPullVsPos, hXResVsPhi, hXResVsPos, hYPull, hYPullVsPos, hYResVsPhi, hYResVsPos, hZPull, hZPullVsPos, hZResVsPhi, hZResVsPos, Geom::pi(), and theName.

                                     :theName(name.c_str()){

    // Position, sigma, residual, pull
    hResX        = dbe_->book1D (theName+"_X_Res", "X residual", 5000, -0.5,0.5);
    hResY        = dbe_->book1D (theName+"_Y_Res", "Y residual", 5000, -1.,1.);
    hResZ        = dbe_->book1D (theName+"_Z_Res", "Z residual", 5000, -1.5,1.5);

    hResXVsEta   = dbe_->book2D(theName+"_XResVsEta", "X residual vs eta",
                            200, -2.5,2.5, 5000, -1.5,1.5);
    hResYVsEta   = dbe_->book2D(theName+"_YResVsEta", "Y residual vs eta",
                            200, -2.5,2.5, 5000, -1.,1.);
    hResZVsEta   = dbe_->book2D(theName+"_ZResVsEta", "Z residual vs eta",
                            200, -2.5,2.5, 5000, -1.5,1.5);
    
    hXResVsPhi   = dbe_->book2D(theName+"_XResVsPhi", "X residual vs phi",
                            100,-Geom::pi(),Geom::pi(), 5000, -0.5,0.5);
    hYResVsPhi   = dbe_->book2D(theName+"_YResVsPhi", "Y residual vs phi",
                            100,-Geom::pi(),Geom::pi(), 5000, -1.,1.);
    hZResVsPhi   = dbe_->book2D(theName+"_ZResVsPhi", "Z residual vs phi",
                            100,-Geom::pi(),Geom::pi(), 5000, -1.5,1.5);
    
    hXResVsPos   = dbe_->book2D(theName+"_XResVsPos", "X residual vs position",
                            10000, -750,750, 5000, -0.5,0.5);    
    hYResVsPos   = dbe_->book2D(theName+"_YResVsPos", "Y residual vs position",
                            10000, -740,740, 5000, -1.,1.);    
    hZResVsPos   = dbe_->book2D(theName+"_ZResVsPos", "Z residual vs position",
                            10000, -1100,1100, 5000, -1.5,1.5);   
    
    hXPull       = dbe_->book1D (theName+"_XPull", "X pull", 600, -2,2);
    hYPull       = dbe_->book1D (theName+"_YPull", "Y pull", 600, -2,2);
    hZPull       = dbe_->book1D (theName+"_ZPull", "Z pull", 600, -2,2);

    hXPullVsPos  = dbe_->book2D (theName+"_XPullVsPos", "X pull vs position",10000, -750,750, 600, -2,2);
    hYPullVsPos  = dbe_->book2D (theName+"_YPullVsPos", "Y pull vs position",10000, -740,740, 600, -2,2);
    hZPullVsPos  = dbe_->book2D (theName+"_ZPullVsPos", "Z pull vs position",10000, -1100,1100, 600, -2,2);
  }

HResolution1DRecHit::HResolution1DRecHit	(	TString	name_,
		TFile *	file
	)		`[inline]`

Definition at line 336 of file Histograms.h.

{}

HResolution1DRecHit::~HResolution1DRecHit ( ) [inline]

Definition at line 338 of file Histograms.h.

{}

Member Function Documentation

void HResolution1DRecHit::Fill	(	double	x,
		double	y,
		double	z,
		double	dx,
		double	dy,
		double	dz,
		double	errx,
		double	erry,
		double	errz,
		double	eta,
		double	phi
	)		`[inline]`

Definition at line 340 of file Histograms.h.

References gather_cfg::cout, ExpressReco_HICollisions_FallBack::e, MonitorElement::Fill(), hResX, hResXVsEta, hResY, hResYVsEta, hResZ, hResZVsEta, hXPull, hXPullVsPos, hXResVsPhi, hXResVsPos, hYPull, hYPullVsPos, hYResVsPhi, hYResVsPos, hZPull, hZPullVsPos, hZResVsPhi, hZResVsPos, ExpressReco_HICollisions_FallBack::x, ExpressReco_HICollisions_FallBack::y, and z.

                                    {
    
    double rx = dx/x, ry = dy/y, rz = dz/z;

    hResX->Fill(rx);
    hResY->Fill(ry);
    hResZ->Fill(rz);
    
    hResXVsEta->Fill(eta,rx);
    hResYVsEta->Fill(eta,ry);
    hResZVsEta->Fill(eta,rz);
    
    hXResVsPhi->Fill(phi,rx);
    hYResVsPhi->Fill(phi,ry);
    hZResVsPhi->Fill(phi,rz);
    
    hXResVsPos->Fill(x,rx);
    hYResVsPos->Fill(y,ry);
    hZResVsPos->Fill(z,rz);
    
    if(errx < 1e-6)
      std::cout << "NO proper error set for X: "<<errx<<std::endl;
    else{
      hXPull->Fill(dx/errx);
      hXPullVsPos->Fill(x,dx/errx);
    }
    
    if(erry < 1e-6) 
      std::cout << "NO proper error set for Y: "<<erry<<std::endl;
    else{
      hYPull->Fill(dy/erry);
      hYPullVsPos->Fill(y,dy/erry);
    }
    if(errz < 1e-6)
      std::cout << "NO proper error set for Z: "<<errz<<std::endl; 
    else{
      hZPull->Fill(dz/errz);
      hZPullVsPos->Fill(z,dz/errz);
    }
  }