A set of histograms of residuals and pulls for 1D RecHits. More...

#include <Histograms.h>

Public Member Functions
void	Fill (float distSimHit, float thetaSimHit, float distFESimHit, float distRecHit, float etaSimHit, float phiSimHit, float errRecHit, int station)
	HRes1DHit (std::string name_, DQMStore *dbe_, bool doall=true, bool local=true)
	~HRes1DHit ()
Public Attributes
bool	_doall
MonitorElement *	hDist
MonitorElement *	hPull
MonitorElement *	hPullSt [4]
MonitorElement *	hPullVsAngle
MonitorElement *	hPullVsDistFE
MonitorElement *	hPullVsPos
MonitorElement *	hRes
MonitorElement *	hResSt [4]
MonitorElement *	hResVsAngle
MonitorElement *	hResVsDistFE
MonitorElement *	hResVsEta
MonitorElement *	hResVsPhi
MonitorElement *	hResVsPos
TString	name

Detailed Description

A set of histograms of residuals and pulls for 1D RecHits.

Definition at line 28 of file Histograms.h.

Constructor & Destructor Documentation

HRes1DHit::HRes1DHit	(	std::string	name_,
		DQMStore *	dbe_,
		bool	doall = `true`,
		bool	local = `true`
	)		`[inline]`

Definition at line 30 of file Histograms.h.

References _doall, DQMStore::book1D(), DQMStore::book2D(), hDist, hPull, hPullSt, hPullVsAngle, hPullVsDistFE, hPullVsPos, hRes, hResSt, hResVsAngle, hResVsDistFE, hResVsEta, hResVsPhi, hResVsPos, and DQMStore::setCurrentFolder().

                                                                             {
      std::string pre ="1D_";
      pre += name_;
      _doall = doall;
      // Position, sigma, residual, pull
      //sprintf (histo_n, "1D_%s_hDist",N); sprintf(histo_t, "1D RHit distance from wire");
      //      hDist=0; hDist = dbe_->book1D(histo_n, histo_t, 100, 0,2.5);
      if (local) dbe_->setCurrentFolder("DT/1DRecHits/");
      else dbe_->setCurrentFolder("DT/1DRecHits/Res/");
      
      if(doall){
        hDist=0; hDist = dbe_->book1D(pre + "_hDist" ,"1D RHit distance from wire", 100, 0,2.5);
        //hDist       = new TH1F ("1D_"+N+"_hDist", "1D RHit distance from wire", 100, 0,2.5);
        hResVsAngle = 0; hResVsAngle   = dbe_->book2D(pre+"_hResVsAngle", "1D RHit residual vs impact angle",100, 0.,1.2, 150, -0.5,0.5);    
        hResVsDistFE = 0; hResVsDistFE = dbe_->book2D(pre+"_hResVsDistFE", "1D RHit residual vs FE distance", 100, 0.,400., 150, -0.5,0.5);    
        if(!local) dbe_->setCurrentFolder("DT/1DRecHits/Pull/");
        hPullVsPos= 0; hPullVsPos  = dbe_->book2D (pre+"_hPullVsPos", "1D RHit pull vs position", 100, 0,2.5, 100, -5,5);
        hPullVsAngle = 0; hPullVsAngle  = dbe_->book2D (pre+"_hPullVsAngle", "1D RHit pull vs impact angle",100, 0.,+1.2, 100, -5,5);
        hPullVsDistFE = 0; hPullVsDistFE  = dbe_->book2D (pre+"_hPullVsDistFE", "1D RHit pull vs FE distance", 100, 0., 400., 100, -5,5);
      }
      if(!local) dbe_->setCurrentFolder("DT/1DRecHits/Res/");
      hRes=0; hRes = dbe_->book1D(pre + "_hRes","1D RHit residual", 300, -0.5,0.5);
      hResSt[0] = 0; hResSt[0] = dbe_->book1D(pre + "_hResMB1","1D RHit residual", 300, -0.5,0.5);
      hResSt[1] = 0; hResSt[1] = dbe_->book1D(pre + "_hResMB2","1D RHit residual", 300, -0.5,0.5);
      hResSt[2] = 0; hResSt[2] = dbe_->book1D(pre + "_hResMB3","1D RHit residual", 300, -0.5,0.5);
      hResSt[3] = 0; hResSt[3] = dbe_->book1D(pre + "_hResMB4","1D RHit residual", 300, -0.5,0.5);

      //hRes        = new TH1F ("1D_"+N+"_hRes", "1D RHit residual", 300, -1.5,1.5);
      hResVsEta=0; hResVsEta = dbe_->book2D(pre +"_hResVsEta" , "1D RHit residual vs eta", 50, -1.25,1.25,150,-0.5,0.5);
      //hResVsEta   = new TH2F("1D_"+N+"_hResVsEta", "1D RHit residual vs eta",50, -1.25,1.25, 150, -1.5,1.5);
      hResVsPhi = 0; hResVsPhi   = dbe_->book2D(pre+"_hResVsPhi" , "1D RHit residual vs phi",100, -3.2, 3.2, 150, -0.5,0.5);
      //hResVsPhi   = new TH2F("1D_"+N+"_hResVsPhi", "1D RHit residual vs phi", 100, -3.2, 3.2, 150, -1.5,1.5);
      hResVsPos = 0; hResVsPos   = dbe_->book2D(pre+"_hResVsPos", "1D RHit residual vs position",100, 0, 2.5, 150, -0.5,0.5);    
      //hResVsPos   = new TH2F("1D_"+N+"_hResVsPos", "1D RHit residual vs position",100, 0, 2.5, 150, -1.5,1.5);    
      if (!local) dbe_->setCurrentFolder("DT/1DRecHits/Pull/");
      hPull =0; hPull       = dbe_->book1D (pre+"_hPull", "1D RHit pull", 100, -5,5);
      hPullSt[0] = 0; hPullSt[0] = dbe_->book1D(pre + "_hPullMB1","1D RHit residual", 100, -5,5);
      hPullSt[1] = 0; hPullSt[1] = dbe_->book1D(pre + "_hPullMB2","1D RHit residual", 100, -5,5);
      hPullSt[2] = 0; hPullSt[2] = dbe_->book1D(pre + "_hPullMB3","1D RHit residual", 100, -5,5);
      hPullSt[3] = 0; hPullSt[3] = dbe_->book1D(pre + "_hPullMB4","1D RHit residual", 100, -5,5);
    }

HRes1DHit::~HRes1DHit ( ) [inline]

Definition at line 88 of file Histograms.h.

                {
      //     delete hDist;
      //     delete hRes;
      //     delete hResVsEta;
      //     delete hResVsPhi;
      //     delete hResVsPos;
      //     delete hPull;
    }

Member Function Documentation

void HRes1DHit::Fill	(	float	distSimHit,
		float	thetaSimHit,
		float	distFESimHit,
		float	distRecHit,
		float	etaSimHit,
		float	phiSimHit,
		float	errRecHit,
		int	station
	)		`[inline]`

Definition at line 97 of file Histograms.h.

References _doall, gather_cfg::cout, MonitorElement::Fill(), hDist, hPull, hPullSt, hPullVsAngle, hPullVsDistFE, hPullVsPos, hRes, hResSt, hResVsAngle, hResVsDistFE, hResVsEta, hResVsPhi, and hResVsPos.

Referenced by DTRecHitQuality::compute().

                           {
      // Reso, pull
      float res = distRecHit-distSimHit;
      if(_doall){
        hDist->Fill(distRecHit);
        hResVsAngle->Fill(thetaSimHit,res);
        hResVsDistFE->Fill(distFESimHit,res);
      }
      hRes->Fill(res); hResSt[station-1]->Fill(res);
      hResVsEta->Fill(etaSimHit,res);
      hResVsPhi->Fill(phiSimHit,res);
      hResVsPos->Fill(distSimHit,res);
      if(errRecHit!=0) {
        float pull=res/errRecHit;
        hPull->Fill(pull);hPullSt[station-1]->Fill(pull);
        if(_doall){
          hPullVsPos->Fill(distSimHit,pull);
          hPullVsAngle->Fill(thetaSimHit,pull);
          hPullVsDistFE->Fill(distFESimHit,pull);
        }
      }
      else std::cout<<"Error: RecHit error = 0" << std::endl;
    }