Class for Time walk calibration. More...

#include <TimeWalkCalibration.h>

Public Member Functions
	TimeWalkCalibration ()
	Constructor.

virtual	~TimeWalkCalibration ()
	Destructor.

virtual void	setUseDB (bool useDB)
	setting for use database or text mode

virtual void	setStoreHisto (bool storeHist)
	Store Histogram or not.

virtual void	setMinimumNDF (double minndf)
	minimum number of degree freedom requirement for track

virtual void	setMinimumPval (double minPval)
	minimum chi2 prob requirement for track

virtual void	InputFileNames (std::string inputname)
	input root file names, output of collector

virtual void	InputTWFileName (std::string tw_name)
	Input tw file name, incase use text mode.

virtual void	OutputTWFileName (std::string tw_out_name)
	output tw file name, for text mode

virtual void	OutputT0FileName (std::string t0_out_name)
	output t0 file name, for text mode

virtual void	setMode (unsigned short mode)
	Set time walk mode.

virtual void	execute ()
	function to run algorithm

Protected Member Functions
virtual bool	calibrate ()
	Run algorithm.

virtual void	CreateHisto ()
	Create histo for calibrate.

virtual void	Write ()
	save calibration

virtual void	updateT0 ()
	update constant term to t0 database.

virtual void	readTW ()
	read tw from database

virtual void	fitToExponentialFunc (TH1D *h1)
	fit tw histogram

Private Attributes
TH1D *	m_h1 [300]
	Mean of residual as function of ADC of each board.

TH2D *	m_h2 [300]
	2D histogram of residual vs ADC for each board

double	m_xmin = 0.07
	minimum value cut of drift length.

double	m_ndfmin = 20
	minimum number of degree of freedom required for track.

double	m_Pvalmin = 0.
	minimum number of Prob(chi2) of fitted track.

double	m_constTerm [300] = {0.}
	const term in fitting, it will be added to T0 instead tw

std::vector< float >	m_tw_old [300]
	tw list old.

std::vector< float >	m_tw_new [300]
	tw list new.

bool	m_storeHisto = true
	Store all Histogram or not.

bool	m_useDB = false
	flag to switch btw text mode and database.

unsigned short	m_flag [300]
	flag for fit status

std::string	m_InputTWFileName = "tw.dat"
	Old tw file name.

std::string	m_InputT0FileName = "t0.dat"
	Old t0 file name.

std::string	m_outputTWFileName = "tw_new.dat"
	Output tw file name for time walk.

std::string	m_outputT0FileName = "t0.dat"
	Output tw file name for time walk.

std::string	m_InputRootFileName = "rootfile/output*"
	root input file name.

int	m_firstExperiment = 0
	First experiment.

int	m_firstRun = 0
	First run.

int	m_lastExperiment = -1
	Last experiment.

int	m_lastRun = -1
	Last run.

unsigned short	m_twParamMode_old
	=0 for P0/Sqrt(ADC); =1 for P0Exp(-P1ADC).

unsigned short	m_twParamMode_new = 1
	=0 for P0/Sqrt(ADC); =1 for P0Exp(-P1ADC).

unsigned short	m_nTwParams_new = 2
	No.

unsigned short	m_nTwParams_old
	No.

Detailed Description

Class for Time walk calibration.

Time walk is calibrated for each board as function of ADC

Definition at line 27 of file TimeWalkCalibration.h.

Constructor & Destructor Documentation

◆ TimeWalkCalibration()

TimeWalkCalibration ( )

Constructor.

Definition at line 28 of file TimeWalkCalibration.cc.

{
  //  setDescription("CDC Time Walk Calibration");
}

◆ ~TimeWalkCalibration()

virtual ~TimeWalkCalibration ( )

inlinevirtual

Destructor.

Definition at line 32 of file TimeWalkCalibration.h.

32{}

Member Function Documentation

◆ calibrate()

bool calibrate ( )

protectedvirtual

Run algorithm.

Definition at line 95 of file TimeWalkCalibration.cc.

{
  gROOT->SetBatch(1);
  readTW();
  CreateHisto();
  TF1* fold = nullptr;
  if (m_twParamMode_old == 0)
    fold = new TF1("fold", "[0]/sqrt(x)", 0, 500);
  else if (m_twParamMode_old == 1)
    fold = new TF1("fold", "[0]*exp(-1*[1]*x)", 0, 500);
  else
    B2FATAL("Mode " << m_twParamMode_new << " haven't implemented yet.");
 
  B2INFO("Old time walk formular: ");
  fold->Print();
  B2INFO("Time walk mode new: " << m_twParamMode_new << " with " << m_nTwParams_new << " parameters");
  //  double p3, p1;
  m_tw_new[0].resize(m_nTwParams_new, 0.0); //for board 0, no available
  for (int ib = 1; ib < 300; ++ib) {
    m_flag[ib] = 1;
    B2DEBUG(199, "Board ID:" << ib);
    m_h2[ib]->SetDirectory(0);
 
    // ignore if histogram is low statistic
    if (m_h2[ib]->GetEntries() < 500) { m_flag[ib] = 0; continue;}
    // do slide fit
    m_h2[ib]->FitSlicesY(0, 1, -1, 10);
    TString name = m_h2[ib]->GetName();
    TString hm_name = name + "_1";
    m_h1[ib] = (TH1D*)gDirectory->Get(hm_name);
    if (!m_h1[ib]) {m_flag[ib] = 0; continue;}
    m_h1[ib]->SetDirectory(0);
 
    //if number of poist after slide fit quite low, ignore
    if (m_h1[ib]->GetEntries() < 5) {
      m_flag[ib] = 0;
      B2WARNING("Low statistic, number of points after slice fit: " << m_h1[ib]->GetEntries());
      continue;
    }
 
    // Add previous correction to this
    for (int p = 0; p < m_nTwParams_old; ++p) {
      fold->SetParameter(p, m_tw_old[ib][p]);
    }
    m_h1[ib]->Add(fold);
    if (m_twParamMode_new == 0) {
      TF1* func = new TF1("func", "[0]+[1]/sqrt(x)", 0, 500);
      func->SetParameters(-4, 28);
      m_h1[ib]->Fit("func", "MQ", "", 20, 150);
    }//fit for mode 0
    else if (m_twParamMode_new == 1) {
      fitToExponentialFunc(m_h1[ib]);//do fitting for mode 1
    } else {
      B2FATAL("Mode " << m_twParamMode_new << " is not available, please check again");
    }
 
    // read fitted parameters
    TF1* f1 = m_h1[ib]->GetFunction("func");
    if (!f1) { m_flag[ib] = 0; continue;}
    m_constTerm[ib] = f1->GetParameter(0);
    m_tw_new[ib].resize(m_nTwParams_new, 0.0);
    for (int i = 1; i <= m_nTwParams_new; ++i) {
      m_tw_new[ib][i - 1] = f1->GetParameter(i);
    }
    B2DEBUG(199, "Prob of fitting:" << f1->GetProb());
    B2DEBUG(199, "Fitting Param 0-1:" << f1->GetParameter(0) << " - " << f1->GetParameter(1));
  }
  //Write histogram to file
  if (m_storeHisto) {
    B2INFO("Storing histogram");
    TFile* fhist = new TFile("tw_histo.root", "recreate");
    TDirectory* old = gDirectory;
    TDirectory* h1D = old->mkdir("h1D");
    TDirectory* h2D = old->mkdir("h2D");
    h1D->cd();
    for (int ib = 1; ib < 300; ++ib) {
      if (!m_h1[ib] || m_flag[ib] != 1) continue;
      m_h1[ib]->SetDirectory(0);
      if (m_h1[ib]->GetEntries() < 5) continue;
 
      m_h1[ib]->SetMinimum(-5);
      m_h1[ib]->SetMaximum(15);
      m_h1[ib]->Write();
    }
    B2DEBUG(199, "Store 2D histogram");
    h2D->cd();
    for (int ib = 1; ib < 300; ++ib) {
      if (m_h2[ib]) {
        m_h2[ib]->SetDirectory(0);
        m_h2[ib]->Write();
      }
    }
    fhist->Close();
    B2INFO("All hitograms were stored");
  }
  Write();
  updateT0();
  return true;
}

◆ CreateHisto()

void CreateHisto ( )

protectedvirtual

Create histo for calibrate.

Definition at line 32 of file TimeWalkCalibration.cc.

{
  double x;
  double t_mea;
  double w;
  double t_fit;
  double ndf;
  double Pval;
  unsigned short adc;
  //  int board;
  int IWire;
  int lay;
 
  //  auto& tree = getObject<TTree>("cdcCalib");
  TChain* tree = new TChain("tree");
  tree->Add(m_InputRootFileName.c_str());
 
  tree->SetBranchAddress("lay", &lay);
  tree->SetBranchAddress("IWire", &IWire);
  tree->SetBranchAddress("x_u", &x);
  tree->SetBranchAddress("t", &t_mea);
  tree->SetBranchAddress("t_fit", &t_fit);
  tree->SetBranchAddress("weight", &w);
  tree->SetBranchAddress("ndf", &ndf);
  tree->SetBranchAddress("Pval", &Pval);
  tree->SetBranchAddress("adc", &adc);
 
  /* Disable unused branch */
  std::vector<TString> list_vars = {"lay", "IWire", "x_u", "t", "t_fit",  "weight", "Pval", "ndf", "adc"};
  tree->SetBranchStatus("*", 0);
 
  for (TString brname : list_vars) {
    tree->SetBranchStatus(brname, 1);
  }
 
 
  static CDCGeometryPar& cdcgeo = CDCGeometryPar::Instance();
  double halfCSize[56];
  for (int i = 0; i < 56; ++i) {
    double R = cdcgeo.senseWireR(i);
    double nW = cdcgeo.nWiresInLayer(i);
    halfCSize[i] = M_PI * R / nW;
  }
  int nEntries = tree->GetEntries();
  B2INFO("Number of entry: " << nEntries);
  //Histogram of each board
  for (int i = 0; i < 300; ++i) {
    m_h2[i] = new TH2D(Form("board_%d", i), Form("board_%d", i), 50, 0., 500, 60, -30, 30);
  }
 
  //read data
  for (int i = 0; i < nEntries; ++i) {
    tree->GetEntry(i);
    double xmax = halfCSize[lay] - 0.12;
    if ((fabs(x) < m_xmin) || (fabs(x) > xmax)
        || (ndf < m_ndfmin)
        || (Pval < m_Pvalmin)) continue; /*select good region*/
 
    m_h2[cdcgeo.getBoardID(WireID(lay, IWire))]->Fill(adc, fabs(t_mea) - fabs(t_fit));
  }
  B2INFO("Reading data and filling histogram are done");
}

◆ execute()

virtual void execute ( )

inlinevirtual

function to run algorithm

Definition at line 61 of file TimeWalkCalibration.h.

      {
        calibrate();
      }

◆ fitToExponentialFunc()

void fitToExponentialFunc ( TH1D * h1 )

protectedvirtual

fit tw histogram

Definition at line 263 of file TimeWalkCalibration.cc.

{
  h1->SetDirectory(0);
  int max = h1->GetMaximumBin();
  double maxX = h1->GetBinCenter(max);
  double maxY = h1->GetBinContent(max);
  B2DEBUG(199, "Max: id - x - y : " << max << "  " << maxX << "  " << maxY);
 
  //search for p0
  double p0 = -1;
  h1->Fit("pol0", "MQE", "", maxX + 125, 500);
  if (h1->GetFunction("pol0")) {
    p0 = h1->GetFunction("pol0")->GetParameter(0);
  }
  //create histo = old-p0;
  // fit with expo function to find intial parameters
  TH1D* hshift  = new TH1D("hshift", "shift", h1->GetNbinsX(), 0, 500);
  hshift->SetDirectory(0);
  for (int i = 0; i <= h1->GetNbinsX(); ++i) {
    hshift->SetBinContent(i, h1->GetBinContent(i) - p0);
  }
  hshift->Fit("expo", "MQE", "", 0, 400);
  double p1 = maxY + p0;
  double p2 = -0.04;
  if (hshift->GetFunction("expo")) {
    p1 =  exp(hshift->GetFunction("expo")->GetParameter(0));
    p2 =  hshift->GetFunction("expo")->GetParameter(1);
  }
 
  // fit with function
  TF1* func = new TF1("func", "[0]+ [1]*exp(-1*[2]*x)", 0, 500);
  func->SetParameters(p0, p1, -1 * p2);
  func->SetParLimits(0, -5, 5);
  func->SetParLimits(1, -5, 500);
  func->SetParLimits(2, -0.01, 0.1);
  h1->Fit("func", "MQE", "", 10, 400);
}

◆ InputFileNames()

virtual void InputFileNames ( std::string inputname )

inlinevirtual

input root file names, output of collector

Definition at line 42 of file TimeWalkCalibration.h.

42{m_InputRootFileName.assign(inputname);}

◆ InputTWFileName()

virtual void InputTWFileName ( std::string tw_name )

inlinevirtual

Input tw file name, incase use text mode.

Definition at line 44 of file TimeWalkCalibration.h.

44{m_InputTWFileName.assign(tw_name);}

Belle2::CDC::TimeWalkCalibration::m_InputTWFileName

std::string m_InputTWFileName

Old tw file name.

Definition: TimeWalkCalibration.h:95

◆ OutputT0FileName()

virtual void OutputT0FileName ( std::string t0_out_name )

inlinevirtual

output t0 file name, for text mode

Definition at line 48 of file TimeWalkCalibration.h.

48{m_outputT0FileName.assign(t0_out_name);}

Belle2::CDC::TimeWalkCalibration::m_outputT0FileName

std::string m_outputT0FileName

Output tw file name for time walk.

Definition: TimeWalkCalibration.h:98

◆ OutputTWFileName()

virtual void OutputTWFileName ( std::string tw_out_name )

inlinevirtual

output tw file name, for text mode

Definition at line 46 of file TimeWalkCalibration.h.

46{m_outputTWFileName.assign(tw_out_name);}

Belle2::CDC::TimeWalkCalibration::m_outputTWFileName

std::string m_outputTWFileName

Output tw file name for time walk.

Definition: TimeWalkCalibration.h:97

◆ readTW()

void readTW ( )

protectedvirtual

read tw from database

Definition at line 300 of file TimeWalkCalibration.cc.

{
  if (m_useDB) {
    DBObjPtr<CDCTimeWalks> oldDB;
    m_twParamMode_old = oldDB->getTwParamMode();
    for (int ib = 0; ib < 300; ++ib) {
      m_tw_old[ib] = oldDB->getTimeWalkParams(ib);
    }
  } else {
    //For text mode
    std::string fileName1 = "/data/cdc" + m_InputTWFileName;
    std::string fileName = FileSystem::findFile(fileName1);
    ifstream ifs;
    if (fileName == "") {
      fileName = FileSystem::findFile(m_InputTWFileName);
    }
    if (fileName == "") {
      B2FATAL("CDCGeometryPar: " << fileName1 << " not exist!");
    } else {
      B2INFO("Time Walk Calibration: open " << fileName);
      ifs.open(fileName.c_str());
      if (!ifs) B2FATAL("CDCGeometryPar: cannot open " << fileName1 << " !");
    }
 
    unsigned short iBoard(0);
    int nRead(0);
    double dumy(0.);// dumy2(0.), dumy3(0.);
    ifs >> m_twParamMode_old >> m_nTwParams_old;
    while (ifs >> iBoard) {
      m_tw_old[iBoard].resize(m_nTwParams_old);
      for (int i = 0; i < m_nTwParams_old; ++i) {
        ifs >> dumy;
        m_tw_old[iBoard][i] = dumy;
      }
      ++nRead;
    }
    if (nRead - 1 != 299) B2FATAL("#lines read-in (=" << nRead << ") is not equal #boards 299 !");
    ifs.close();
    B2INFO("TW from database: mode = " << m_twParamMode_old << "  with " << m_nTwParams_old << " parameters");
  }
}

◆ setMinimumNDF()

virtual void setMinimumNDF ( double minndf )

inlinevirtual

minimum number of degree freedom requirement for track

Definition at line 38 of file TimeWalkCalibration.h.

38{m_ndfmin = minndf;}

◆ setMinimumPval()

virtual void setMinimumPval ( double minPval )

inlinevirtual

minimum chi2 prob requirement for track

Definition at line 40 of file TimeWalkCalibration.h.

40{m_Pvalmin = minPval;}

◆ setMode()

virtual void setMode ( unsigned short mode )

inlinevirtual

Set time walk mode.

Definition at line 50 of file TimeWalkCalibration.h.

      {
        m_twParamMode_new = mode;
        if (mode == 0)
          m_nTwParams_new = 1;
        else if (mode == 1)
          m_nTwParams_new = 2;
        else
          B2FATAL("Mode hasn't implemented yet");
      }

◆ setStoreHisto()

virtual void setStoreHisto ( bool storeHist )

inlinevirtual

Store Histogram or not.

Definition at line 36 of file TimeWalkCalibration.h.

36{m_storeHisto = storeHist;}

◆ setUseDB()

virtual void setUseDB ( bool useDB )

inlinevirtual

setting for use database or text mode

Definition at line 34 of file TimeWalkCalibration.h.

34{m_useDB = useDB; }

◆ updateT0()

void updateT0 ( )

protectedvirtual

update constant term to t0 database.

Definition at line 230 of file TimeWalkCalibration.cc.

{
  B2INFO("Add constant term into T0 database");
  static CDCGeometryPar& cdcgeo = CDCGeometryPar::Instance();
  ofstream ofs(m_outputT0FileName.c_str());
  DBImportObjPtr<CDCTimeZeros> tz;// = new CDCTimeZeros();
  tz.construct();
  double T0;
  //correct T0 and write
  for (int ilay = 0; ilay < 56; ++ilay) {
    for (unsigned int iwire = 0; iwire < cdcgeo.nWiresInLayer(ilay); ++iwire) {
      WireID wireid(ilay, iwire);
      int bID = cdcgeo.getBoardID(wireid);
      T0 = cdcgeo.getT0(wireid);
      if (m_flag[bID] == 1) {
        ofs <<  ilay << "\t" << iwire << "\t" << T0 - m_constTerm[bID] << std::endl;
        if (m_useDB)
          tz->setT0(wireid, T0 - m_constTerm[bID]);
      } else {
        ofs <<  ilay << "\t" << iwire << "\t" << T0  << std::endl;
        if (m_useDB)
          tz->setT0(wireid, T0 - m_constTerm[bID]);
      }
    }
  }
  ofs.close();
  if (m_useDB) {
    IntervalOfValidity iov(m_firstExperiment, m_firstRun,
                           m_lastExperiment, m_lastRun);
    tz.import(iov);
  }
  B2RESULT("updated constant term of tw correction to T0 constant.");
}

◆ Write()

void Write ( )

protectedvirtual

save calibration

Definition at line 194 of file TimeWalkCalibration.cc.

{
  B2INFO("Store calibrated time walk parameters");
  int nfailure(-1);//for no exist board; board 0;
  //  CDCTimeWalks* dbTw = new CDCTimeWalks();
  DBImportObjPtr<CDCTimeWalks> dbTw;
  dbTw.construct();
 
  for (int ib = 0; ib < 300; ++ib) {
    // If fitting fail and new mode is same as previous input mode. use old param
    if (m_flag[ib] != 1) {
      nfailure += 1;
      if (m_twParamMode_old == m_twParamMode_new) {
        dbTw->setTimeWalkParams(ib, m_tw_old[ib]);
      } else {
        //and even calibrated fail but mode is different from previous.
        //in this case, param is zero
        m_tw_new[ib].resize(m_nTwParams_new, 0.0);
        dbTw->setTimeWalkParams(ib, m_tw_new[ib]);
      }
    } else {
      //Use new param if board is successfuly calibrated
      dbTw->setTimeWalkParams(ib, m_tw_new[ib]);
    }
  }
  dbTw->setTwParamMode(m_twParamMode_new);
  dbTw->outputToFile(m_outputTWFileName);
  if (m_useDB) {
    IntervalOfValidity iov(m_firstExperiment, m_firstRun,
                           m_lastExperiment, m_lastRun);
    dbTw.import(iov);
  }
 
  B2RESULT("Failure to calibrate time-walk for " << nfailure << " board");
  B2RESULT("Time-walk coeff. table wrote to " << m_outputTWFileName.c_str());
}

Member Data Documentation

◆ m_constTerm

double m_constTerm[300] = {0.}

private

const term in fitting, it will be added to T0 instead tw

Definition at line 89 of file TimeWalkCalibration.h.

◆ m_firstExperiment

int m_firstExperiment = 0

private

First experiment.

Definition at line 100 of file TimeWalkCalibration.h.

◆ m_firstRun

int m_firstRun = 0

private

First run.

Definition at line 101 of file TimeWalkCalibration.h.

◆ m_flag

unsigned short m_flag[300]

private

flag for fit status

Definition at line 94 of file TimeWalkCalibration.h.

◆ m_h1

TH1D* m_h1[300]

private

Mean of residual as function of ADC of each board.

Definition at line 81 of file TimeWalkCalibration.h.

◆ m_h2

TH2D* m_h2[300]

private

2D histogram of residual vs ADC for each board

Definition at line 82 of file TimeWalkCalibration.h.

◆ m_InputRootFileName

std::string m_InputRootFileName = "rootfile/output*"

private

root input file name.

Definition at line 99 of file TimeWalkCalibration.h.

◆ m_InputT0FileName

std::string m_InputT0FileName = "t0.dat"

private

Old t0 file name.

Definition at line 96 of file TimeWalkCalibration.h.

◆ m_InputTWFileName

std::string m_InputTWFileName = "tw.dat"

private

Old tw file name.

Definition at line 95 of file TimeWalkCalibration.h.

◆ m_lastExperiment

int m_lastExperiment = -1

private

Last experiment.

Definition at line 102 of file TimeWalkCalibration.h.

◆ m_lastRun

int m_lastRun = -1

private

Last run.

Definition at line 103 of file TimeWalkCalibration.h.

◆ m_ndfmin

double m_ndfmin = 20

private

minimum number of degree of freedom required for track.

Definition at line 85 of file TimeWalkCalibration.h.

◆ m_nTwParams_new

unsigned short m_nTwParams_new = 2

private

No.

of tw parameters. for new database

Definition at line 106 of file TimeWalkCalibration.h.

◆ m_nTwParams_old

unsigned short m_nTwParams_old

private

No.

of tw parameters. in old database

Definition at line 107 of file TimeWalkCalibration.h.

◆ m_outputT0FileName

std::string m_outputT0FileName = "t0.dat"

private

Output tw file name for time walk.

Definition at line 98 of file TimeWalkCalibration.h.

◆ m_outputTWFileName

std::string m_outputTWFileName = "tw_new.dat"

private

Output tw file name for time walk.

Definition at line 97 of file TimeWalkCalibration.h.

◆ m_Pvalmin

double m_Pvalmin = 0.

private

minimum number of Prob(chi2) of fitted track.

Definition at line 86 of file TimeWalkCalibration.h.

◆ m_storeHisto

bool m_storeHisto = true

private

Store all Histogram or not.

Definition at line 92 of file TimeWalkCalibration.h.

◆ m_tw_new

std::vector<float> m_tw_new[300]

private

tw list new.

Definition at line 91 of file TimeWalkCalibration.h.

◆ m_tw_old

std::vector<float> m_tw_old[300]

private

tw list old.

Definition at line 90 of file TimeWalkCalibration.h.

◆ m_twParamMode_new

unsigned short m_twParamMode_new = 1

private

=0 for P0/Sqrt(ADC); =1 for P0*Exp(-P1*ADC).

Definition at line 105 of file TimeWalkCalibration.h.

◆ m_twParamMode_old

unsigned short m_twParamMode_old

private

=0 for P0/Sqrt(ADC); =1 for P0*Exp(-P1*ADC).

Definition at line 104 of file TimeWalkCalibration.h.

◆ m_useDB

bool m_useDB = false

private

flag to switch btw text mode and database.

Definition at line 93 of file TimeWalkCalibration.h.

◆ m_xmin

double m_xmin = 0.07

private

minimum value cut of drift length.

Definition at line 84 of file TimeWalkCalibration.h.

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

cdc/calibration/include/TimeWalkCalibration.h
cdc/calibration/src/TimeWalkCalibration.cc

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TimeWalkCalibration()

◆ ~TimeWalkCalibration()

Member Function Documentation

◆ calibrate()

◆ CreateHisto()

◆ execute()

◆ fitToExponentialFunc()

◆ InputFileNames()

◆ InputTWFileName()

◆ OutputT0FileName()

◆ OutputTWFileName()

◆ readTW()

◆ setMinimumNDF()

◆ setMinimumPval()

◆ setMode()

◆ setStoreHisto()

◆ setUseDB()

◆ updateT0()

◆ Write()

Member Data Documentation

◆ m_constTerm

◆ m_firstExperiment

◆ m_firstRun

◆ m_flag

◆ m_h1

◆ m_h2

◆ m_InputRootFileName

◆ m_InputT0FileName

◆ m_InputTWFileName

◆ m_lastExperiment

◆ m_lastRun

◆ m_ndfmin

◆ m_nTwParams_new

◆ m_nTwParams_old

◆ m_outputT0FileName

◆ m_outputTWFileName

◆ m_Pvalmin

◆ m_storeHisto

◆ m_tw_new

◆ m_tw_old

◆ m_twParamMode_new

◆ m_twParamMode_old

◆ m_useDB

◆ m_xmin