Calibrator Class Reference

An instance of this class is created for each sub-level of the TRT by the TRTCalibrator. More...

#include <Calibrator.h>

Collaboration diagram for Calibrator:

Public Member Functions
	Calibrator ()
	Calibrator (int, const std::string &, int, int, const std::string &, const std::string &, float)
	Constructor.
	Calibrator (const Calibrator &)=delete
Calibrator &	operator= (const Calibrator &)=delete
	~Calibrator ()
	Destructor.
int	Simple1dHist (float, float, int, float)
	A 1D histograming function.
int	Simple2dHist (float, float, int, float, float, int, float, float)
	A 2D histograming function.
float	AccumulativeMean (float, float, float)
	Increments a single bin in all three histograms in a sub-module.
int	AddHit (const std::string &, const databundle &, int *, bool)
	Adds hits to a sub-module either in the form of a single straw hit or a histogram containing all the hits in a single straw.
TDirectory *Calibrate	ATLAS_NOT_THREAD_SAFE (TDirectory *, const std::string &, const std::string &, caldata *)
	Creates root histograms, performs the t0 and R-t calibration for a given sub-module and writes the resulting histograms to a root directory tree (if not told otherwise).
float FitRt	ATLAS_NOT_THREAD_SAFE (const std::string &, const std::string &, TH2F *, TDirectory *)
	Makes the R-t fit.
float FitTimeResidual	ATLAS_NOT_THREAD_SAFE (const std::string &, TH1F *)
	Makes the time residual fit.
float FitResidual	ATLAS_NOT_THREAD_SAFE (const std::string &, TH1F *)
	Makes the residual fit.
void	WriteStat (TDirectory *)
	Creates an ntuple with entries containing data associated with the sub-modules in a sub level.
void	DumpConstants ()
std::string	GetBinnedRt (const std::string &)
std::string	PrintStat ()
	Prints some sub-level statistics.
std::string	PrintInfo ()
bool	HasKey (const std::string &) const
	...
int	UpdateOldConstants ()
	...
bool	CheckSelection (int)
	Checks if a given sub-module is selected for calibration.
bool	Skip ()
	...
std::string	GetOptString () const
	Creates a string summarizing what is being done at this sub-level.
int	GetNDirs ()
	...
std::string	GetSelString ()
	Creates a string summarizing which modules are being calibrated at this sub-level.
float	oldt0 (const std::string &key)
float	xmean (const std::string &key)
float	ymean (const std::string &key)
float	zmean (const std::string &key)

Public Attributes
std::map< std::string, caldata >	data
	A map between the sub-module identifier string and the calibration data structure (caldata).
bool	dort {}
	if true an rt fit is made, if false the value from the level above is used
bool	dot0 {}
	if true a time residual fit is made, if false the value from the level above is used
bool	dores {}
	if true a residual fit is made
bool	nort {}
	if true the old rt parameters are copied to the new ones
bool	not0 {}
	if true the old t0 valus is copied to the new one
bool	usebref {}
	if true chip reference t0 values (offset from board mean) are used
bool	bequiet {}
	if true no histograms are written to the root file
bool	printlog {}
	if true a log entry is prined for each sub-modile in this sub-level
bool	printt0 {}
	if true a t0 entry in the calibration output file is prined for each sub-module in this sub-level
bool	printrt {}
	if true an rt entry in the calibration output file is prined for each sub-module in this sub-level
bool	usep0 {}
	if true the 0th order coeficcient of the rt fit function is not set to 0 in the calibration output file
bool	floatp3 {}
	if true the 3rd order coeficcient of the rt fit function is not fixed to 0
bool	useshortstraws {}
	if true a shift of -0.75 ns is applied for straws in layer ==0, strawlayer < 9, and when doing calibration in a granuralyty different from chip or straw
std::set< int >	selection
	a set containing the sub-modules to be calibrated
int	level {}
	the sub-level of the Calibrator instance

Private Attributes
std::string	m_name
	The name of the Calibrator instance.
std::string	m_rtbinning
	The direction to do the R-t binning.
int	m_minrtstat {}
	minimum number of hits in a sub-module required to do an R-t calibration
int	m_mint0stat {}
	minimum number of hits in a sub-module required to do a t0 calibration
float	m_t0shift {}
	the t0 shift
int	m_nbinsr {}
	number of r-bins in the 2D rt histogram
int	m_nbinst {}
	number of t-bins in the 2D rt histogram
int	m_nbinstres {}
	number of bins in the 1D time residual histogram
int	m_nbinsres {}
	number of bins in the 1D residual histogram
float	m_minr {}
	lower limit of r in 2D rt histogram
float	m_maxr {}
	upper limit of r in 2D rt histogram
float	m_mint {}
	lower limit of t in 2D rt histogram
float	m_maxt {}
	upper limit of t in 2D rt histogram
float	m_mintres {}
	lower limit of 1D time residual histogram
float	m_maxtres {}
	upper limit of 1D time residual histogram
float	m_minres {}
	lower limit of 1D residual histogram
float	m_maxres {}
	upper limit of 1D residual histogram
std::map< std::string, TDirectory * >	m_hdirs
	...
std::map< std::string, TH1F * >	m_resHists
	...
std::map< std::string, TH1F * >	m_tresHists
	...
std::map< std::string, TH2F * >	m_rtHists
	...
std::map< std::string, RtGraph * >	m_rtgraphs {}
	...
int	m_ntreshits {}
	...
int	m_nreshits {}
	...
int	m_nrthits {}
	...
int	m_nhits {}
	...
bool	m_isdines
	true if the rt relation is Dines'

Detailed Description

An instance of this class is created for each sub-level of the TRT by the TRTCalibrator.

It contains three histograms for each sub-module in the level:

• a 2D R-t histogram (caldata::rthist)
• a 1D time residual histogram (caldata::treshist)
• a 1D residual histogram (caldata::reshist)

The histograms are inside a caldata structure which also contains other variables with information about the calibration of a sub-module. There is one caldata instance for each sub-module in the level. They are contained in a map (Calibrator::data) with a key string which identifies the sub-module.

Methods exists for adding hits to the histograms (Calibrator::AddHits) and for fitting them (Calibrator::Calibrate).

Definition at line 178 of file Calibrator.h.

Constructor & Destructor Documentation

Calibrator() [1/3]

Calibrator::Calibrator

(

)

Definition at line 339 of file Calibrator.cxx.

  : dort(false),
    dot0(false),
    dores(false),
    nort(false),
    not0(false),
    usebref(false),
    bequiet(false),
    printlog(false),
    printt0(false),
    printrt(false),
    usep0(false),
    floatp3(false),
    useshortstraws(true),
 
    m_name ("None"),
    m_rtbinning ("None"),
    m_minrtstat (-10),
    m_mint0stat (-10),
    m_t0shift (-100.),
 
    m_nbinsr(100),
    m_nbinst(55),
    m_nbinstres(100),
    m_nbinsres(100),
    m_minr(0),
    m_maxr(2),
    m_mint(-5),
    m_maxt(50),
    m_mintres(-10),
    m_maxtres(10),
    m_minres(-0.6),
    m_maxres(0.6),
 
    m_ntreshits(0),
    m_nreshits(0),
    m_nrthits(0),
    m_nhits(0),
 
    m_isdines ( false)
{
  selection.insert(-3);
  level =-10;
}

Calibrator() [2/3]

Calibrator::Calibrator	(	int	lev,
		const std::string &	nme,
		int	mint0,
		int	minrt,
		const std::string &	rtr,
		const std::string &	rtb,
		float	t0sft )

Constructor.

Parameters

[in]	lev	sub level number 0=TRT 1=barrel/ec 2=layer 3=phi-module 4=board 5=chip 6=straw
[in]	nme	the name of the calibrator
[in]	mint0	minimum number of hits in a sub-module required to do a t0 calibration
[in]	minrt	minimum number of hits in a sub-module required to do an rt calibration
[in]	rtr	which rt-relation to use
[in]	rtb	which rt binning to use

Definition at line 384 of file Calibrator.cxx.

  :  dort(false),
     dot0(false),
     dores(false),
     nort(false),
     not0(false),
     usebref(false),
     bequiet(false),
     printlog(false),
     printt0(false),
     printrt(false),
     usep0(false),
     floatp3(false),
     useshortstraws(true),
 
     m_name(nme), 
     m_rtbinning(rtb),
     m_minrtstat(minrt),
     m_mint0stat(mint0),
     m_t0shift(t0sft),
 
     m_nbinsr(100),
     m_nbinst(55),
     m_nbinstres(100),
     m_nbinsres(100),
     m_minr(0),
     m_maxr(2),
     m_mint(-5),
     m_maxt(50),
     m_mintres(-10),
     m_maxtres(10),
     m_minres(-0.6),
     m_maxres(0.6),
 
     m_ntreshits(0),
     m_nreshits(0),
     m_nrthits(0),
     m_nhits(0),
 
     m_isdines ( rtr.find("dines")!=std::string::npos)
{
  level=lev;
  selection.insert(-3);
}

Calibrator() [3/3]

Calibrator::Calibrator ( const Calibrator & )

delete

~Calibrator()

Calibrator::~Calibrator

(

)

Destructor.

Definition at line 429 of file Calibrator.cxx.

                       {
}

Member Function Documentation

AccumulativeMean()

float Calibrator::AccumulativeMean	(	float	n,
		float	oldmean,
		float	newvalue )

Increments a single bin in all three histograms in a sub-module.

Parameters

[in]	key	the identifier of the histogram
[in]	d	the databundle containing single straw hit data

Returns

1

Definition at line 447 of file Calibrator.cxx.

                                                                        {
  if((int)n==0) return oldmean;
  return oldmean*((n-1)/n)+newvalue/n;
}

AddHit()

int Calibrator::AddHit	(	const std::string &	key,
		const databundle &	d,
		int *	binhist,
		bool	makehist )

Adds hits to a sub-module either in the form of a single straw hit or a histogram containing all the hits in a single straw.

Parameters

[in]	key	the identifier of the sub-module
[in]	d	the databundle containing single straw hit data
[in]	binhist	an integer arry containing the histogram data for a single straw

Returns

1 if the sub-module is seen for the first time and 0 if it has been seen before

Definition at line 962 of file Calibrator.cxx.

                                                                                             {
  
  int tresbin=Simple1dHist(m_mintres,m_maxtres,m_nbinstres,d.tres);
  int resbin=Simple1dHist(m_minres,m_maxres,m_nbinsres,d.res);
  int rtbin=Simple2dHist(m_mint,m_maxt,m_nbinst,m_minr,m_maxr,m_nbinsr,d.t,d.r);
  int npop;
  int ibin;
  int value;
  
  //if the it is the first hit or histogram
  if (data.find(key) == data.end()){
    
    //create a new histogram object
    caldata* hist=new caldata(makehist,m_nbinst,m_nbinsr);    
 
    //out of memory?
    if (hist == nullptr){
      std::cout << "OUT OF MEMORY!" << std::endl;
      return -1;
    }
 
    hist->ntres=0;
    hist->nres=0;
    hist->nrt=0;
    hist->sumt0=0;
 
    //zero histogram bins
    for (int i=0;i<m_nbinsr*m_nbinst+200;i++) {
      if (makehist) hist->rthist[i]=0;
      if (i<100) {
        if (makehist) hist->m_treshist[i]=0;
        if (makehist) hist->reshist[i]=0;
      }
    }
 
    if (binhist==nullptr){ //if it is a hit
      
      if (tresbin>=0){ 
        if (makehist) hist->m_treshist[tresbin]=d.weight; //add value to bin
        hist->ntres=1; //set bin content to 1
        m_nhits=1;
      }
      if (resbin>=0){ 
        if (makehist) hist->reshist[resbin]=1;
        hist->nres=1;
      }
      if (rtbin>=0){ 
        if (makehist) hist->rthist[rtbin]=1;
        hist->nrt=1;
      }
 
      hist->sumt0=d.t0;
    }
    else { //if it is a histogram
 
      npop=binhist[0]; //the number of populated (non-zero) bins
      if(level<2) std::cout << " Calibrator::AddHit called for key " << key << " add histogram with "<< npop << " non-zero bins " << std::endl;
      
      for(int ipop=2;ipop<2*npop+2;ipop=ipop+2) { //loop over the data
 
        ibin=binhist[ipop]; //bin number
        value=binhist[ipop+1]; //value
 
        if (ibin<100) { //timeresidual histogram
          m_ntreshits+=value;
          hist->ntres+=value;
          m_nhits+=value;
          if (makehist) hist->m_treshist[ibin]=(float)value;
          hist->sumt0+=d.t0*value;
        }
        else if (ibin>=100 && ibin<200) { //residual histogram
          m_nreshits+=value;
          hist->nres+=value;
          if (makehist) hist->reshist[ibin-100]=(float)value;
        }
        else { //rt histogram
          m_nrthits+=value;
          hist->nrt+=value;
          if (makehist) hist->rthist[ibin-200]=(float)value;
        }
 
      }    
    }
 
    //initialize parameters
    if (level>0) hist->det=d.det; else hist->det=-3; //
    if (level>1) hist->lay=d.lay; else hist->lay=-3;
    if (level>2) hist->mod=d.mod; else hist->mod=-3;
    if (level>3) hist->brd=d.brd; else hist->brd=-3;
    if (level>4) hist->chp=d.chp; else hist->chp=-3;
    if (level>5) hist->stw=d.stw; else hist->stw=-3;
    if (level>2) hist->stl=d.stl; else hist->stl=-3;
    if (level>5) hist->sid=d.sid; else hist->sid=-3;
    hist->res=0;
    hist->resMean=0;
    hist->reserr=0;
    hist->tres=0;
    hist->tresMean=0;
    hist->t0=0;
    if (level==5) hist->reft0=d.rt0; else hist->reft0=0;
    hist->t0off=0;
    hist->t0err=0;
    hist->nhits=1;
    hist->t0fittype=0;
    
    hist->x=d.x; //straw x coordinate
    hist->y=d.y; //straw y coordinate
    hist->z=d.z; //straw z coordinate
 
    hist->oldt02=d.t0; //straw old t0 value
    hist->oldrtpar[0]=d.rtpar[0];
    hist->oldrtpar[1]=d.rtpar[1];
    hist->oldrtpar[2]=d.rtpar[2];
    hist->oldrtpar[3]=d.rtpar[3];
 
    hist->sumx=0;
    hist->sumy=0;
    hist->sumz=0;
 
    hist->rtflag=false;
    hist->t0flag=false;
    hist->calflag=false;
 
    
 
    for (unsigned int i =0; i < 4; i++){
      hist->rtpar[i]=-10.0;
    }
 
    data[key]=*hist; //save the histogram in the map 
 
 
 
    data[key].oldt02 = AccumulativeMean(data[key].nhits, data[key].oldt02, d.t0); //update old t0 m_mean value
 
    data[key].nhits++; // increment m_nhits
 
    if(level<2) std::cout << "   First hit for det " << data[key].det << " lay " << data[key].lay << " t0 " << data[key].t0 << " upd t0 " << data[key].oldt02 << " nhits: " << data[key].nhits << std::endl;    
 
    delete hist;
 
    return 1;
}
  else { //if not the first hit or histogram
 
    //increment histogram bins
    if (binhist==nullptr){
      
      if (tresbin>=0){
        if (makehist) data[key].m_treshist[tresbin]=data[key].m_treshist[tresbin]+d.weight;
        data[key].ntres++;
        m_nhits=1;
      }
      if (resbin>=0){
        if (makehist) data[key].reshist[resbin]++;
        data[key].nres++;
      }
      if (rtbin>=0){
        if (makehist) data[key].rthist[rtbin]++;
        data[key].nrt++;
      }
 
      data[key].sumt0+=d.t0;
 
    }
 
    else {
 
      npop=binhist[0];
 
      for(int ipop=2;ipop<2*npop+2;ipop=ipop+2) {
 
        ibin=binhist[ipop];
        value=binhist[ipop+1];
 
        if (ibin<100) {
          m_ntreshits+=value;
          m_nhits+=value;
          data[key].ntres+=value;
          if (makehist) data[key].m_treshist[ibin]+=(float)value;
          data[key].sumt0+=d.t0*value;
        }
        else if (ibin>=100 && ibin<200) {
          m_nreshits+=value;
          data[key].nres+=value;
          if (makehist) data[key].reshist[ibin-100]+=(float)value;
        }
        else {
          m_nrthits+=value;
          data[key].nrt+=value;
          if (makehist) data[key].rthist[ibin-200]+=(float)value;
        }
 
      }    
      
    }
 
        
    data[key].sumx+=d.x;
    data[key].sumy+=d.y;
    data[key].sumz+=d.z;
    data[key].oldt02 = AccumulativeMean(data[key].nhits, data[key].oldt02, d.t0);
    data[key].oldrtpar[0] = AccumulativeMean(data[key].nhits, data[key].oldrtpar[0], d.rtpar[0]);
    data[key].oldrtpar[1] = AccumulativeMean(data[key].nhits, data[key].oldrtpar[1], d.rtpar[1]);
    data[key].oldrtpar[2] = AccumulativeMean(data[key].nhits, data[key].oldrtpar[2], d.rtpar[2]);
    data[key].oldrtpar[3] = AccumulativeMean(data[key].nhits, data[key].oldrtpar[3], d.rtpar[3]);
    data[key].x = AccumulativeMean(data[key].nhits, data[key].x, d.x);
    data[key].y = AccumulativeMean(data[key].nhits, data[key].y, d.y);
    data[key].z = AccumulativeMean(data[key].nhits, data[key].z, d.z);
 
    data[key].nhits++; //increment hit counts
    
    return 0;
  }    
 
 
}

ATLAS_NOT_THREAD_SAFE() [1/4]

float FitRt Calibrator::ATLAS_NOT_THREAD_SAFE	(	const std::string &	,
		const std::string &	,
		TH2F *	,
		TDirectory *	)

Makes the R-t fit.

Parameters

[in]	key	the identifier of the sub-module to calibrate
[in]	rtHist	the 2D root histogram with the R-t data

Returns

the t0 from the R-t fit

ATLAS_NOT_THREAD_SAFE() [2/4]

float FitResidual Calibrator::ATLAS_NOT_THREAD_SAFE	(	const std::string &	,
		TH1F *	)

Makes the residual fit.

Parameters

[in]	key	the identifier of the sub-module to calibrate
[in]	resHist	the 1D root histogram with the residuals

Returns

the mean value of the residual fit

ATLAS_NOT_THREAD_SAFE() [3/4]

float FitTimeResidual Calibrator::ATLAS_NOT_THREAD_SAFE	(	const std::string &	,
		TH1F *	)

Makes the time residual fit.

Parameters

[in]	key	the identifier of the sub-module to calibrate
[in]	tresHist	the 1D root histogram with the time residuals

Returns

the mean value of the time residual fit

ATLAS_NOT_THREAD_SAFE() [4/4]

TDirectory *Calibrate Calibrator::ATLAS_NOT_THREAD_SAFE	(	TDirectory *	,
		const std::string &	,
		const std::string &	,
		caldata *	)

Creates root histograms, performs the t0 and R-t calibration for a given sub-module and writes the resulting histograms to a root directory tree (if not told otherwise).

The new t0 values are calculated here.

Parameters

[in]	dir	the root directory of the sub level above
[in]	key	the identifier of the sub-module to calibrate
[in]	caldata_above	the caldata object from the sub-module above the one to be calibrated

Returns

the root directory where the histgrams were written

CheckSelection()

bool Calibrator::CheckSelection

(

int

level

)

Checks if a given sub-module is selected for calibration.

Parameters

[in]

level

the number of the sub-module to check

Returns

returns true if the module is found

Definition at line 457 of file Calibrator.cxx.

                                        {
  if (selection.find(level)!=selection.end() || selection.find(-3)!=selection.end()) return true;
  else return false;
}

DumpConstants()

void Calibrator::DumpConstants

(

)

Definition at line 1225 of file Calibrator.cxx.

                              {
 
  std::ofstream dumpfile( "calib_constants_out.txt", std::ios::out | std::ios::app);
 
  for (const std::pair<const std::string, caldata>& p : data) {
    dumpfile << p.first << " " << p.second.t0 << " " << std::endl;
    
  }
  
  dumpfile.close();
  
}

GetBinnedRt()

std::string Calibrator::GetBinnedRt

(

const std::string &

key

)

Definition at line 514 of file Calibrator.cxx.

                                                     {
 
  int nbins=data[key].rtgraph->tval.size();
  double maxt=data[key].rtgraph->tval[nbins-1];
  double mint=data[key].rtgraph->tval[0];
 
  std::string brt = "";
  brt += Form("%f %f %i ",mint,maxt,nbins);
 
  for (int ip=0; ip<nbins;ip++){
    brt += Form("%f ",data[key].rtgraph->rval[ip]);
  }
  return brt;
}

GetNDirs()

int Calibrator::GetNDirs

(

)

...

Definition at line 502 of file Calibrator.cxx.

                        {
  return data.size();
}

GetOptString()

std::string Calibrator::GetOptString

(

)

const

Creates a string summarizing what is being done at this sub-level.

Returns

the string

Definition at line 485 of file Calibrator.cxx.

                                         {
  std::string optstr="";
  if (dort) optstr += 'R';
  if (dot0) optstr += 'T';
  if (printlog) optstr +='P';
  if (printt0) optstr += 'F';
  if (printrt) optstr += 'G';
  if (bequiet) optstr += 'Q';
  if (usebref) optstr += 'B';
  if (useshortstraws) optstr += 'S';
  if (usep0 && dort) optstr += '0';
  if (floatp3 && dort) optstr += '3';
  if (!(dort || dot0 || usebref || bequiet || printlog || printt0 || printrt)) optstr += 'N';
  
  return optstr;
}

GetSelString()

std::string Calibrator::GetSelString

(

)

Creates a string summarizing which modules are being calibrated at this sub-level.

Returns

the string

Definition at line 506 of file Calibrator.cxx.

                                  {
  std::string selstr="";
  if (selection.find(-3)!=selection.end()) selstr="*";
  else if (selection.find(-4)!=selection.end()) selstr="-";
  else for (int isel : selection)  selstr += std::string(Form("%i",isel));
  return selstr;
}

HasKey()

bool Calibrator::HasKey

(

const std::string &

key

)

const

...

Definition at line 452 of file Calibrator.cxx.

                                                  {
 
  return data.find(key) != data.end();
}

oldt0()

float Calibrator::oldt0 ( const std::string & key )

inline

Definition at line 333 of file Calibrator.h.

{return data[key].sumt0/float(data[key].ntres);}

operator=()

Calibrator & Calibrator::operator= ( const Calibrator & )

delete

PrintInfo()

std::string Calibrator::PrintInfo

(

)

Definition at line 467 of file Calibrator.cxx.

                               {
  std::string yn[2]={"NO","YES"};
  std::string info = std::string(Form("CONFIGURATION %-16s: dort=%-3s, dot0=%-3s, dores=%-3s, selection=",m_name.data(),yn[dort].data(),yn[dot0].data(),yn[dores].data()));
  for (int isel : selection) {
    if (isel==-3) info += std::string("all"); 
    else if (isel==-4) info += std::string("none"); 
    else info += std::string(Form("%2i,",isel)); 
  } 
  return info;
}

PrintStat()

std::string Calibrator::PrintStat

(

)

Prints some sub-level statistics.

Definition at line 478 of file Calibrator.cxx.

                               {
  int units = (int)data.size();
  if(units==0) units++;
  std::string info = std::string(Form("STATISTICS %16s: nunits=%8i, nhits=%9i, hits/unit=%11.1f", m_name.data(), units, m_nhits, (float)m_nhits/(float)units ));
  return info;
}

Simple1dHist()

int Calibrator::Simple1dHist	(	float	min,
		float	max,
		int	nbins,
		float	value )

A 1D histograming function.

Parameters

[in]	min	histogram lower limit
[in]	max	histogram upper limit
[in]	nbins	nymber of bins
[in]	value	the value

Returns

the histogram bin

Definition at line 432 of file Calibrator.cxx.

                                                                        {
  if ( (value<min) || (value>max) ) return -1;
  if ( max-min == 0 ) return -1;
  int binno=(int)(nbins*((value-min)/(max-min)));
  return binno;
}

Simple2dHist()

int Calibrator::Simple2dHist	(	float	minx,
		float	maxx,
		int	nbinsx,
		float	miny,
		float	maxy,
		int	nbinsy,
		float	valuex,
		float	valuey )

A 2D histograming function.

Parameters

[in]	minx	histogram lower x-limit
[in]	maxx	histogram upper x-limit
[in]	nbinsx	nymber of x bins
[in]	valuex	the x value
[in]	miny	histogram lower y-limit
[in]	maxy	histogram upper y-limit
[in]	nbinsy	nymber of y bins
[in]	valuey	the y value

Returns

the histogram bin in the form (biny*nbinsx+binx)

Definition at line 439 of file Calibrator.cxx.

                                                                                                                              {
  if ( (valuex<minx) || (valuex>maxx) ||  (valuey<miny) || (valuey>maxy) ) return -1;
  if(maxx-minx == 0 || maxy-miny==0) return -1;
  int binnox=(int)(nbinsx*((valuex-minx)/(maxx-minx)));
  int binnoy=(int)(nbinsy*((valuey-miny)/(maxy-miny)));
  return binnoy*nbinsx+binnox;
}

Skip()

bool Calibrator::Skip

(

)

...

Definition at line 462 of file Calibrator.cxx.

                     {
  if (selection.find(-4)!=selection.end()) return true;
  else return false;
}

UpdateOldConstants()

int Calibrator::UpdateOldConstants

(

)

...

Definition at line 529 of file Calibrator.cxx.

                                  {
 
  std::string line;
  std::ifstream oldconstfile("calib_constants_in.txt");
  std::string key;
  float t0;
  
  if (oldconstfile.is_open())
    {
      while (!oldconstfile.eof() )
      {
        oldconstfile >> key >> t0;
        if (data.find(key) != data.end()){
          std::cout << "UPDATED OLD T0: " << key << " " << data[key].oldt02 << " -> " << t0 << std::endl; 
          data[key].oldt02=t0;
        }
      }
      oldconstfile.close();
      return 0;
    }  
  else {
    std::cout << "NO OLD T0S FILE FOUND. USING AVERAGE VALUES!" << std::endl;
    return -1;
  } 
}

WriteStat()

void Calibrator::WriteStat

(

TDirectory *

dir

)

Creates an ntuple with entries containing data associated with the sub-modules in a sub level.

Parameters

[in]

dir

the root directory where to write the ntuple

Definition at line 1180 of file Calibrator.cxx.

                                         {
  dir->cd();
  TNtuple* stattup = new TNtuple(Form("%stuple",m_name.data()),"statistics","det:lay:mod:brd:chp:sid:stl:stw:rtflag:t0flag:t0:oldt0:rt0:dt0:t0offset:ftype:nhits:nt0:nrt:res:resMean:dres:tres:tresMean:x:y:z");
  for(std::map<std::string,caldata>::iterator ihist=data.begin(); ihist!=data.end(); ++ihist){
    if ((ihist->second).calflag) {
      //std::cout << " Writing TNtuple for name: " << m_name << " det: " << (ihist->second).det << " lay " << (ihist->second).lay << " mod: " << (ihist->second).mod << " brd " << (ihist->second).brd << " res width " << (ihist->second).res << " tres mean " << (ihist->second).tresMean << std::endl;
      float const ntvar[27]={
        float((ihist->second).det),
        float((ihist->second).lay),
        float((ihist->second).mod),
        float((ihist->second).brd),
        float((ihist->second).chp),
        float((ihist->second).sid),
        float((ihist->second).stl),
        float((ihist->second).stw),
        float((int)(ihist->second).rtflag),
        float((int)(ihist->second).t0flag),
        (ihist->second).t0,
        float((ihist->second).oldt02),
        //oldt0(ihist->first),
        (ihist->second).reft0,
        (ihist->second).t0err,
        (ihist->second).t0off,
        float((ihist->second).t0fittype),
        (ihist->second).nhits,
        (float)(ihist->second).ntres,
        (float)(ihist->second).nrt,
        (float)(ihist->second).res,
        (float)(ihist->second).resMean,
        (ihist->second).reserr,
        (ihist->second).tres,
        (ihist->second).tresMean,
        (ihist->second).x,
        (ihist->second).y,
        (ihist->second).z
      };      
      stattup->Fill(ntvar);
    }
      
  }
 
  stattup->Write();
 
}

xmean()

float Calibrator::xmean ( const std::string & key )

inline

Definition at line 334 of file Calibrator.h.

{return data[key].sumx/float(data.size());}

ymean()

float Calibrator::ymean ( const std::string & key )

inline

Definition at line 335 of file Calibrator.h.

{return data[key].sumy/float(data.size());}

zmean()

float Calibrator::zmean ( const std::string & key )

inline

Definition at line 336 of file Calibrator.h.

{return data[key].sumz/float(data.size());}

Member Data Documentation

bequiet

bool Calibrator::bequiet {}

if true no histograms are written to the root file

Definition at line 347 of file Calibrator.h.

{};

data

std::map<std::string,caldata> Calibrator::data

A map between the sub-module identifier string and the calibration data structure (caldata).

Definition at line 339 of file Calibrator.h.

dores

bool Calibrator::dores {}

if true a residual fit is made

Definition at line 343 of file Calibrator.h.

{};

dort

bool Calibrator::dort {}

if true an rt fit is made, if false the value from the level above is used

Definition at line 341 of file Calibrator.h.

{};

dot0

bool Calibrator::dot0 {}

if true a time residual fit is made, if false the value from the level above is used

Definition at line 342 of file Calibrator.h.

{};

floatp3

bool Calibrator::floatp3 {}

if true the 3rd order coeficcient of the rt fit function is not fixed to 0

Definition at line 352 of file Calibrator.h.

{};

level

int Calibrator::level {}

the sub-level of the Calibrator instance

Definition at line 355 of file Calibrator.h.

{};

m_hdirs

std::map<std::string,TDirectory*> Calibrator::m_hdirs

private

...

Definition at line 376 of file Calibrator.h.

m_isdines

bool Calibrator::m_isdines

private

true if the rt relation is Dines'

Definition at line 388 of file Calibrator.h.

m_maxr

float Calibrator::m_maxr {}

private

upper limit of r in 2D rt histogram

Definition at line 369 of file Calibrator.h.

{};

m_maxres

float Calibrator::m_maxres {}

private

upper limit of 1D residual histogram

Definition at line 375 of file Calibrator.h.

{};

m_maxt

float Calibrator::m_maxt {}

private

upper limit of t in 2D rt histogram

Definition at line 371 of file Calibrator.h.

{};

m_maxtres

float Calibrator::m_maxtres {}

private

upper limit of 1D time residual histogram

Definition at line 373 of file Calibrator.h.

{};

m_minr

float Calibrator::m_minr {}

private

lower limit of r in 2D rt histogram

Definition at line 368 of file Calibrator.h.

{};

m_minres

float Calibrator::m_minres {}

private

lower limit of 1D residual histogram

Definition at line 374 of file Calibrator.h.

{};

m_minrtstat

int Calibrator::m_minrtstat {}

private

minimum number of hits in a sub-module required to do an R-t calibration

Definition at line 361 of file Calibrator.h.

{};

m_mint

float Calibrator::m_mint {}

private

lower limit of t in 2D rt histogram

Definition at line 370 of file Calibrator.h.

{};

m_mint0stat

int Calibrator::m_mint0stat {}

private

minimum number of hits in a sub-module required to do a t0 calibration

Definition at line 362 of file Calibrator.h.

{};

m_mintres

float Calibrator::m_mintres {}

private

lower limit of 1D time residual histogram

Definition at line 372 of file Calibrator.h.

{};

m_name

std::string Calibrator::m_name

private

The name of the Calibrator instance.

Definition at line 359 of file Calibrator.h.

m_nbinsr

int Calibrator::m_nbinsr {}

private

number of r-bins in the 2D rt histogram

Definition at line 364 of file Calibrator.h.

{};

m_nbinsres

int Calibrator::m_nbinsres {}

private

number of bins in the 1D residual histogram

Definition at line 367 of file Calibrator.h.

{};

m_nbinst

int Calibrator::m_nbinst {}

private

number of t-bins in the 2D rt histogram

Definition at line 365 of file Calibrator.h.

{};

m_nbinstres

int Calibrator::m_nbinstres {}

private

number of bins in the 1D time residual histogram

Definition at line 366 of file Calibrator.h.

{};

m_nhits

int Calibrator::m_nhits {}

private

...

Definition at line 384 of file Calibrator.h.

{};

m_nreshits

int Calibrator::m_nreshits {}

private

...

Definition at line 382 of file Calibrator.h.

{};

m_nrthits

int Calibrator::m_nrthits {}

private

...

Definition at line 383 of file Calibrator.h.

{};

m_ntreshits

int Calibrator::m_ntreshits {}

private

...

Definition at line 381 of file Calibrator.h.

{};

m_resHists

std::map<std::string,TH1F*> Calibrator::m_resHists

private

...

Definition at line 377 of file Calibrator.h.

m_rtbinning

std::string Calibrator::m_rtbinning

private

The direction to do the R-t binning.

Definition at line 360 of file Calibrator.h.

m_rtgraphs

std::map<std::string,RtGraph*> Calibrator::m_rtgraphs {}

private

...

Definition at line 380 of file Calibrator.h.

m_rtHists

std::map<std::string,TH2F*> Calibrator::m_rtHists

private

...

Definition at line 379 of file Calibrator.h.

m_t0shift

float Calibrator::m_t0shift {}

private

the t0 shift

Definition at line 363 of file Calibrator.h.

{};

m_tresHists

std::map<std::string,TH1F*> Calibrator::m_tresHists

private

...

Definition at line 378 of file Calibrator.h.

nort

bool Calibrator::nort {}

if true the old rt parameters are copied to the new ones

Definition at line 344 of file Calibrator.h.

{};

not0

bool Calibrator::not0 {}

if true the old t0 valus is copied to the new one

Definition at line 345 of file Calibrator.h.

{};

printlog

bool Calibrator::printlog {}

if true a log entry is prined for each sub-modile in this sub-level

Definition at line 348 of file Calibrator.h.

{};

printrt

bool Calibrator::printrt {}

if true an rt entry in the calibration output file is prined for each sub-module in this sub-level

Definition at line 350 of file Calibrator.h.

{};

printt0

bool Calibrator::printt0 {}

if true a t0 entry in the calibration output file is prined for each sub-module in this sub-level

Definition at line 349 of file Calibrator.h.

{};

selection

std::set<int> Calibrator::selection

a set containing the sub-modules to be calibrated

Definition at line 354 of file Calibrator.h.

usebref

bool Calibrator::usebref {}

if true chip reference t0 values (offset from board mean) are used

Definition at line 346 of file Calibrator.h.

{};

usep0

bool Calibrator::usep0 {}

if true the 0th order coeficcient of the rt fit function is not set to 0 in the calibration output file

Definition at line 351 of file Calibrator.h.

{};

useshortstraws

bool Calibrator::useshortstraws {}

if true a shift of -0.75 ns is applied for straws in layer ==0, strawlayer < 9, and when doing calibration in a granuralyty different from chip or straw

Definition at line 353 of file Calibrator.h.

{};

---

The documentation for this class was generated from the following files:

• Calibrator.h
• Calibrator.cxx