#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 (std::string name)

	HResolution1DRecHit (const TString &name_, TFile *file)

	~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 295 of file Histograms.h.

Constructor & Destructor Documentation

HResolution1DRecHit::HResolution1DRecHit ( std::string name )

inline

Definition at line 297 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	(	const TString &	name_,
		TFile *	file
	)

inline

Definition at line 334 of file Histograms.h.

334 {}

HResolution1DRecHit::~HResolution1DRecHit ( )

inline

Definition at line 336 of file Histograms.h.

336 {}

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 338 of file Histograms.h.

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

Referenced by MuonTrackResidualAnalyzer::computeResolution().

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