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 *, 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 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() [3/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() [4/4]

TDirectory *Calibrate Calibrator::ATLAS_NOT_THREAD_SAFE	(	TDirectory *	,
		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