time plots More...

#include <MuonInsideOutAnalysisPlots.h>

Collaboration diagram for Muon::TimePlots:

Public Member Functions
void	book (TDirectory *dir, const TString &prefix, int type)

void	fill (const MuonValidationTimeBlock &times, int index, float betaTrack)

Public Attributes
int	typeIndex

TH1 *	time

TH1 *	error

TH1 *	timeRes

TH1 *	timePull

TH1 *	beta

TH1 *	betaRes

TH1 *	betaPull

TH1 *	invBetaRes

TH1 *	invBetaPull

TH2 *	timeResScat

TH2 *	timePullScat

TH1 *	error_time

TH1 *	errorRT

TH1 *	errorTrack

TH2 *	errorRT_Track

TH1 *	timeResCor

TH1 *	timePullCor

TH1 *	timeRes0_10

TH1 *	timeRes10_70

TH1 *	timeRes70_150

TH1 *	timeRes150_350

TH1 *	timeRes350_800

TH1 *	timePull0_10

TH1 *	timePull10_70

TH1 *	timePull70_150

TH1 *	timePull150_350

TH1 *	timePull350_800

Detailed Description

time plots

Definition at line 62 of file MuonInsideOutAnalysisPlots.h.

Member Function Documentation

◆ book()

void Muon::TimePlots::book	(	TDirectory *	dir,
		const TString &	prefix,
		int	type
	)

Definition at line 80 of file MuonInsideOutAnalysisPlots.cxx.

                                                                      {
     type = type % 1000; // truncate anything above 1000
     bool isMDT = (type == 1 || type == 10);
     float tmin = -15;
     float tmax =  15;
     float timeRange = 5;
     if( isMDT ){
       tmin = -30;
       tmax =  30;
       if( type == 10 ) timeRange = 800;
     }
     time = new TH1F(prefix+"time","time",100,-20.,100.);
     error = new TH1F(prefix+"error","error",100,0.,15.);
     timeRes = new TH1F(prefix+"timeRes","timeRes",200,tmin,tmax);
     timePull = new TH1F(prefix+"pullRes","pullRes",200,-7,7);
     beta = new TH1F(prefix+"beta","beta",100,.0,2.);
     betaRes = new TH1F(prefix+"betaRes","betaRes",200,-1.2,1.2);
     betaPull = new TH1F(prefix+"betaPull","betaPull",200,-7,7);
     invBetaRes = new TH1F(prefix+"invBetaRes","invBetaRes",200,-1.2,1.2);
     invBetaPull = new TH1F(prefix+"invBetaPull","invBetaPull",200,-7,7);
  
  
     if( type == 10 ){
       timeResScat = new TH2F(prefix+"timeResScat","timeResScat",100,tmin,timeRange,100,tmin,tmax);
       timePullScat = new TH2F(prefix+"pullResScat","pullResScat",100,tmin,timeRange,100,-7.,7.);
  
       error_time = new TH2F(prefix+"error_time","error_time",100,tmin,timeRange,100,0.,15.);
       errorRT = new TH1F(prefix+"errorRT","errorRT",100,-20.,100.);
       errorTrack = new TH1F(prefix+"errorTrack","errorTrack",100,-20.,100.);
       errorRT_Track = new TH2F(prefix+"errorRT_Track","errorRT_Track",100,-20.,100.,100,-20.,100.);
       timeResCor = new TH1F(prefix+"timeResCor","timeResCor",200,tmin,tmax);
       timePullCor  = new TH1F(prefix+"timePullCor","timePullCor",200,-7,7);
       timeRes0_10 = new TH1F(prefix+"timeRes0_10","timeRes0_10",200,tmin,tmax);
       timeRes10_70 = new TH1F(prefix+"timeRes10_70","timeRes10_70",200,tmin,tmax);
       timeRes70_150 = new TH1F(prefix+"timeRes70_150","timeRes70_150",200,tmin,tmax);
       timeRes150_350 = new TH1F(prefix+"timeRes150_350","timeRes150_350",200,tmin,tmax);
       timeRes350_800 = new TH1F(prefix+"timeRes350_800","timeRes350_800",200,tmin,tmax);
  
       timePull0_10 = new TH1F(prefix+"timePull0_10","timePull0_10",200,-7,7);
       timePull10_70 = new TH1F(prefix+"timePull10_70","timePull10_70",200,-7,7);
       timePull70_150 = new TH1F(prefix+"timePull70_150","timePull70_150",200,-7,7);
       timePull150_350 = new TH1F(prefix+"timePull150_350","timePull150_350",200,-7,7);
       timePull350_800 = new TH1F(prefix+"timePull350_800","timePull350_800",200,-7,7);
     }else if( type == 1 ){
       timeResCor = new TH1F(prefix+"timeResLargeError","timeResLargeError",200,tmin,tmax);
       timePullCor  = new TH1F(prefix+"timePullLargError","timePullLargError",200,-7,7);
     } 
   }

◆ fill()

void Muon::TimePlots::fill	(	const MuonValidationTimeBlock &	times,
		int	index,
		float	betaTrack
	)

Definition at line 129 of file MuonInsideOutAnalysisPlots.cxx.

                                                                                          {
     float dist = (*times.d)[index];
     int type = (*times.type)[index];
     float err = (*times.err)[index];
     float t = (*times.time)[index];
     MuonBetaCalculationUtils betaUtils(0.);
     float tof       = betaUtils.calculateTof(1.,dist);
     float timeTrack = betaUtils.calculateTof(betaTrack,dist);
     float t0_       = (*times.time)[index]+tof;
     float betaMeas   = betaUtils.calculateBeta(t0_,dist);
     float betaError  = betaUtils.calculateBetaError(t0_,err,dist);
     float invBetaMeas   = betaUtils.calculateInverseBeta(t0_,dist);
     float invBetaError  = betaUtils.calculateInverseBetaError(err,dist);
     time->Fill(t);
     error->Fill(err);
     beta->Fill(betaMeas);
     float res = t0_-timeTrack;
     float pull = res/err;
     if( type != 1 || err < 7 ){
       timeRes->Fill(res);
       timePull->Fill(pull);
     }
     float betaResidual = betaMeas-betaTrack;
     betaRes->Fill(betaResidual);
     betaPull->Fill(betaResidual/betaError);
     float invBetaResidual = invBetaMeas-1./betaTrack;
     invBetaRes->Fill(invBetaResidual);
     invBetaPull->Fill(invBetaResidual/invBetaError);
     
     float fullTime = (*times.timeProp)[index];
     if( type == 10 ){
       error_time->Fill(fullTime,err);
       timeResScat->Fill(fullTime,res);
       timePullScat->Fill(fullTime,pull);
       float errTrk = (*times.avTimeProp)[index];
       errorTrack->Fill(errTrk);
       float errRT = sqrt(err*err-errTrk*errTrk);
       errorRT->Fill(errRT);
       errorRT_Track->Fill(errRT,errTrk);
       float rescor = res-(*times.timeCor)[index];
       float pullcor = rescor/err;
       timeResCor->Fill(rescor);
       timePullCor->Fill(pullcor);
       if( fullTime < 10 )       timeRes0_10->Fill(res);
       else if( fullTime < 50 )  timeRes10_70->Fill(res);
       else if( fullTime < 150 ) timeRes70_150->Fill(res);
       else if( fullTime < 350 ) timeRes150_350->Fill(res);
       else timeRes350_800->Fill(res);
  
       if( fullTime < 10 )       timePull0_10->Fill(pull);
       else if( fullTime < 50 )  timePull10_70->Fill(pull);
       else if( fullTime < 150 ) timePull70_150->Fill(pull);
       else if( fullTime < 350 ) timePull150_350->Fill(pull);
       else timePull350_800->Fill(pull);
  
     }else if( type == 1 && err > 7 ){
       timeResCor->Fill(res);
       timePullCor->Fill(pull);
     }
   }