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