#include <DedxCalib.h>

Inheritance diagram for DedxCalib:

Public Member Functions
	DedxCalib (const std::string &name, ISvcLocator *pSvcLocator)
	~DedxCalib ()
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
virtual void	genNtuple ()=0
virtual void	initializing ()=0
virtual void	BookHists ()=0
virtual void	FillHists ()=0
virtual void	AnalyseHists ()=0
virtual void	WriteHists ()=0

Protected Member Functions
double	cal_dedx_bitrunc (float truncate, std::vector< double > phlist)
double	cal_dedx (float truncate, std::vector< double > phlist)
void	getCurvePar ()
void	set_dEdx (int landau, float dEdx, int trkalg, int runflag, int dedxhit_use, float ptrk, float theta, float pl_rp, int vflag[3], double t0)
void	ReadRecFileList ()

Protected Attributes
IMdcGeomSvc *	geosvc
IDedxCorrecSvc *	exsvc
float	truncate
vector< double >	Curve_Para
vector< double >	Sigma_Para
int	vFlag [3]
double	m_dedx_exp [5]
double	m_ex_sigma [5]
double	m_pid_prob [5]
double	m_chi_ex [5]
vector< string >	m_recFileVector
int	ParticleFlag
int	m_calibflag
int	m_phShapeFlag
std::string	m_eventType
std::string	m_recFileList
std::string	m_rootfile
std::string	m_curvefile

Detailed Description

Definition at line 11 of file DedxCalib.h.

Constructor & Destructor Documentation

◆ DedxCalib()

DedxCalib::DedxCalib	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 12 of file DedxCalib.cxx.

    : Algorithm( name, pSvcLocator ) {
  declareProperty( "ParticleType", ParticleFlag = -1 );
  declareProperty( "EventType", m_eventType = "isBhabha" );
  declareProperty( "CalibFlag", m_calibflag = 63 );
  declareProperty( "PhShapeFlag", m_phShapeFlag = 0 );
  declareProperty( "TruncRate", truncate = 0.7 );
  declareProperty( "RecFileList", m_recFileList = "no file" );
  declareProperty( "RootFile", m_rootfile = "no file" );
  declareProperty( "CurveFile", m_curvefile = std::string( "no rootfile" ) );
}

Referenced by DedxCalibCostheta::DedxCalibCostheta(), DedxCalibDocaEAng::DedxCalibDocaEAng(), DedxCalibEAng::DedxCalibEAng(), DedxCalibEvent::DedxCalibEvent(), DedxCalibLayerGain::DedxCalibLayerGain(), DedxCalibMomentum::DedxCalibMomentum(), DedxCalibRunByRun::DedxCalibRunByRun(), DedxCalibTzero::DedxCalibTzero(), DedxCalibWireGain::DedxCalibWireGain(), and DedxSaveWireGain::DedxSaveWireGain().

◆ ~DedxCalib()

DedxCalib::~DedxCalib ( )

inline

Definition at line 14 of file DedxCalib.h.

14{};

Member Function Documentation

◆ AnalyseHists()

void DedxCalib::AnalyseHists ( )

pure virtual

Implemented in DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 80 of file DedxCalib.cxx.

80{}

Referenced by execute().

◆ BookHists()

void DedxCalib::BookHists ( )

pure virtual

Implemented in DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 76 of file DedxCalib.cxx.

76{}

Referenced by initialize().

◆ cal_dedx()

double DedxCalib::cal_dedx	(	float	truncate,
		std::vector< double >	phlist )

protected

Definition at line 119 of file DedxCalib.cxx.

                                                                     {
  sort( phlist.begin(), phlist.end() );
  int      smpl    = (int)( phlist.size() * ( truncate + 0.05 ) );
  int      min_cut = 0;
  double   qSum    = 0;
  unsigned i       = 0;
  for ( vector<double>::iterator ql = phlist.begin(); ql != phlist.end(); ql++ )
  {
    i++;
    if ( i <= smpl && i >= min_cut ) qSum += ( *ql );
  }
  double trunc   = -99;
  int    usedhit = smpl - min_cut + 1;
  if ( usedhit > 0 ) trunc = qSum / usedhit;
 
  return trunc;
}

◆ cal_dedx_bitrunc()

double DedxCalib::cal_dedx_bitrunc	(	float	truncate,
		std::vector< double >	phlist )

protected

Definition at line 101 of file DedxCalib.cxx.

                                                                             {
  sort( phlist.begin(), phlist.end() );
  int      smpl    = (int)( phlist.size() * ( truncate + 0.05 ) );
  int      min_cut = (int)( phlist.size() * 0.05 + 0.5 );
  double   qSum    = 0;
  unsigned i       = 0;
  for ( vector<double>::iterator ql = phlist.begin(); ql != phlist.end(); ql++ )
  {
    i++;
    if ( i <= smpl && i >= min_cut ) qSum += ( *ql );
  }
  double trunc   = -99;
  int    usedhit = smpl - min_cut + 1;
  if ( usedhit > 0 ) trunc = qSum / usedhit;
 
  return trunc;
}

Referenced by DedxCalibLayerGain::AnalyseHists(), DedxCalibCostheta::FillHists(), DedxCalibMomentum::FillHists(), DedxCalibRunByRun::FillHists(), and DedxCalibTzero::FillHists().

◆ execute()

StatusCode DedxCalib::execute ( )

Definition at line 53 of file DedxCalib.cxx.

                              {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "DedxCalib execute()" << endmsg;
 
  genNtuple();
  FillHists();
  AnalyseHists();
  WriteHists();
 
  return StatusCode::SUCCESS;
}

◆ FillHists()

void DedxCalib::FillHists ( )

pure virtual

Implemented in DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 78 of file DedxCalib.cxx.

78{}

Referenced by execute().

◆ finalize()

StatusCode DedxCalib::finalize ( )

Definition at line 65 of file DedxCalib.cxx.

                               {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "DedxCalib finalize()" << endmsg;
 
  return StatusCode::SUCCESS;
}

◆ genNtuple()

void DedxCalib::genNtuple ( )

pure virtual

Implemented in DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 74 of file DedxCalib.cxx.

74{}

Referenced by execute().

◆ getCurvePar()

void DedxCalib::getCurvePar ( )

protected

Definition at line 137 of file DedxCalib.cxx.

{
  if ( m_curvefile == "no rootfile" )
  {
    cout << "no curve file! can not calculate chi!!! " << endl;
    return;
  }
 
  double Curve[100];
  double Sigma[100];
  int    CurveSize, SigmaSize;
  TFile* f     = new TFile( m_curvefile.c_str() );
  TTree* curve = (TTree*)f->Get( "curve" );
  TTree* sigma = (TTree*)f->Get( "sigma" );
  curve->SetBranchAddress( "CurveSize", &CurveSize );
  curve->SetBranchAddress( "curve", Curve );
  sigma->SetBranchAddress( "SigmaSize", &SigmaSize );
  sigma->SetBranchAddress( "sigma", Sigma );
  curve->GetEntry( 0 );
  sigma->GetEntry( 0 );
 
  for ( int i = 0; i < CurveSize; i++ ) // get the dedx curve parameters
  {
    cout << Curve[i] << endl;
    if ( i == 0 ) vFlag[0] = (int)Curve[i]; // first element is dedx curve version
    else Curve_Para.push_back( Curve[i] );  // dedx curve parameters
  }
  for ( int i = 0; i < SigmaSize; i++ )
  {
    cout << Sigma[i] << endl;
    if ( i == 0 ) vFlag[1] = (int)Sigma[i]; // dedx sigma version: 2: psip; 3:jpsi
    else Sigma_Para.push_back( Sigma[i] );  // dedx sigma parameters
  }
  // if(runflag>0)   vFlag[2]=0;
  // else            vFlag[2]=1;
}

Referenced by initialize().

◆ initialize()

StatusCode DedxCalib::initialize ( )

Definition at line 24 of file DedxCalib.cxx.

                                 {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "DedxCalib initialze()" << endmsg;
 
  StatusCode gesc = service( "MdcGeomSvc", geosvc );
  if ( gesc == StatusCode::SUCCESS ) { log << MSG::INFO << "MdcGeomSvc is running" << endmsg; }
  else
  {
    log << MSG::ERROR << "MdcGeomSvc is failed" << endmsg;
    return StatusCode::SUCCESS;
  }
 
  StatusCode scex = service( "DedxCorrecSvc", exsvc );
  if ( scex == StatusCode::SUCCESS )
  { log << MSG::INFO << "Hi, DedxCorrectSvc is running" << endl; }
  else
  {
    log << MSG::ERROR << "DedxCorrectSvc is failed" << endl;
    return StatusCode::SUCCESS;
  }
  exsvc->set_flag( m_calibflag );
 
  initializing();
  BookHists();
  getCurvePar();
 
  return StatusCode::SUCCESS;
}

◆ initializing()

void DedxCalib::initializing ( )

pure virtual

Implemented in DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 72 of file DedxCalib.cxx.

72{}

Referenced by initialize().

◆ ReadRecFileList()

void DedxCalib::ReadRecFileList ( )

protected

Definition at line 84 of file DedxCalib.cxx.

                                {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "DedxCalib::ReadRecFileList()" << endmsg;
 
  ifstream filea( m_recFileList.c_str(), ifstream::in );
  cout << m_recFileList.c_str() << endl;
 
  string runlist;
  filea >> runlist;
  do {
    m_recFileVector.push_back( runlist );
    cout << runlist.c_str() << endl;
    filea >> runlist;
  } while ( filea );
}

Referenced by DedxCalibCostheta::BookHists(), DedxCalibDocaEAng::BookHists(), DedxCalibEAng::BookHists(), DedxCalibLayerGain::BookHists(), DedxCalibMomentum::BookHists(), DedxCalibRunByRun::BookHists(), DedxCalibTzero::BookHists(), DedxCalibWireGain::BookHists(), and DedxSaveWireGain::BookHists().

◆ set_dEdx()

void DedxCalib::set_dEdx	(	int	landau,
		float	dEdx,
		int	trkalg,
		int	runflag,
		int	dedxhit_use,
		float	ptrk,
		float	theta,
		float	pl_rp,
		int	vflag[3],
		double	t0 )

protected

Definition at line 175 of file DedxCalib.cxx.

                                                                                          {
  // landau: 1:landau distribution; 0:gaus distribution
  // pl_rp: 1.5, not know what
  // getCurvePar(runflag);
  if ( runflag > 0 ) vFlag[2] = 0;
  else vFlag[2] = 1;
 
  // some old data with old methods
  if ( runflag == 1 || runflag == 2 )
    dedx_pid_exp_old( landau, runflag, (float)dEdx, (int)dedxhit_use, (float)ptrk,
                      (float)theta, (float)t0, (float)pl_rp, m_dedx_exp, m_ex_sigma,
                      m_pid_prob, m_chi_ex );
  // for 2009 psip data and after
  else
    dedx_pid_exp( vflag, (float)dEdx, trkalg, (int)dedxhit_use, (float)ptrk, (float)theta,
                  (float)t0, (float)pl_rp, m_dedx_exp, m_ex_sigma, m_pid_prob, m_chi_ex,
                  Curve_Para, Sigma_Para );
}

Referenced by DedxCalibCostheta::FillHists(), DedxCalibMomentum::FillHists(), and DedxCalibTzero::FillHists().

◆ WriteHists()

void DedxCalib::WriteHists ( )

pure virtual

Implemented in DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 82 of file DedxCalib.cxx.

82{}

Referenced by execute().

Member Data Documentation

◆ Curve_Para

vector<double> DedxCalib::Curve_Para

protected

Definition at line 39 of file DedxCalib.h.

Referenced by getCurvePar(), and set_dEdx().

◆ exsvc

IDedxCorrecSvc* DedxCalib::exsvc

protected

Definition at line 29 of file DedxCalib.h.

Referenced by DedxCalibCostheta::FillHists(), DedxCalibDocaEAng::FillHists(), DedxCalibEAng::FillHists(), DedxCalibLayerGain::FillHists(), DedxCalibMomentum::FillHists(), DedxCalibRunByRun::FillHists(), DedxCalibTzero::FillHists(), DedxCalibWireGain::FillHists(), DedxSaveWireGain::FillHists(), and initialize().

◆ geosvc

IMdcGeomSvc* DedxCalib::geosvc

protected

Definition at line 28 of file DedxCalib.h.

Referenced by DedxCalibEvent::genNtuple(), and initialize().

◆ m_calibflag

int DedxCalib::m_calibflag

protected

Definition at line 53 of file DedxCalib.h.

Referenced by DedxCalib(), and initialize().

◆ m_chi_ex

double DedxCalib::m_chi_ex[5]

protected

Definition at line 48 of file DedxCalib.h.

Referenced by DedxCalibCostheta::FillHists(), DedxCalibMomentum::FillHists(), DedxCalibTzero::FillHists(), and set_dEdx().

◆ m_curvefile

std::string DedxCalib::m_curvefile

protected

Definition at line 58 of file DedxCalib.h.

Referenced by DedxCalib(), and getCurvePar().

◆ m_dedx_exp

double DedxCalib::m_dedx_exp[5]

protected

Definition at line 45 of file DedxCalib.h.

Referenced by set_dEdx().

◆ m_eventType

std::string DedxCalib::m_eventType

protected

Definition at line 55 of file DedxCalib.h.

Referenced by DedxCalib(), and DedxCalibEvent::genNtuple().

◆ m_ex_sigma

double DedxCalib::m_ex_sigma[5]

protected

Definition at line 46 of file DedxCalib.h.

Referenced by set_dEdx().

◆ m_phShapeFlag

int DedxCalib::m_phShapeFlag

protected

Definition at line 54 of file DedxCalib.h.

Referenced by DedxCalibDocaEAng::AnalyseHists(), DedxCalibEAng::AnalyseHists(), DedxCalibLayerGain::AnalyseHists(), DedxCalibWireGain::AnalyseHists(), DedxCalib(), DedxCalibDocaEAng::WriteHists(), DedxCalibEAng::WriteHists(), DedxCalibLayerGain::WriteHists(), and DedxCalibWireGain::WriteHists().

◆ m_pid_prob

double DedxCalib::m_pid_prob[5]

protected

Definition at line 47 of file DedxCalib.h.

Referenced by set_dEdx().

◆ m_recFileList

std::string DedxCalib::m_recFileList

protected

Definition at line 56 of file DedxCalib.h.

Referenced by DedxCalib(), and ReadRecFileList().

◆ m_recFileVector

vector<string> DedxCalib::m_recFileVector

protected

Definition at line 50 of file DedxCalib.h.

Referenced by DedxCalibRunByRun::BookHists(), DedxCalibCostheta::FillHists(), DedxCalibDocaEAng::FillHists(), DedxCalibEAng::FillHists(), DedxCalibLayerGain::FillHists(), DedxCalibMomentum::FillHists(), DedxCalibRunByRun::FillHists(), DedxCalibTzero::FillHists(), DedxCalibWireGain::FillHists(), DedxSaveWireGain::FillHists(), and ReadRecFileList().

◆ m_rootfile

std::string DedxCalib::m_rootfile

protected

Definition at line 57 of file DedxCalib.h.

Referenced by DedxCalib(), DedxCalibCostheta::WriteHists(), DedxCalibDocaEAng::WriteHists(), DedxCalibEAng::WriteHists(), DedxCalibLayerGain::WriteHists(), DedxCalibMomentum::WriteHists(), DedxCalibRunByRun::WriteHists(), DedxCalibTzero::WriteHists(), and DedxCalibWireGain::WriteHists().

◆ ParticleFlag

int DedxCalib::ParticleFlag

protected

Definition at line 52 of file DedxCalib.h.

Referenced by DedxCalib(), DedxCalibCostheta::FillHists(), DedxCalibMomentum::FillHists(), DedxCalibTzero::FillHists(), and DedxCalibEvent::genNtuple().

◆ Sigma_Para

vector<double> DedxCalib::Sigma_Para

protected

Definition at line 40 of file DedxCalib.h.

Referenced by getCurvePar(), and set_dEdx().

◆ truncate

float DedxCalib::truncate

protected

Definition at line 32 of file DedxCalib.h.

Referenced by DedxCalibLayerGain::AnalyseHists(), cal_dedx(), cal_dedx_bitrunc(), DedxCalib(), DedxCalibCostheta::FillHists(), DedxCalibMomentum::FillHists(), DedxCalibRunByRun::FillHists(), and DedxCalibTzero::FillHists().

◆ vFlag

int DedxCalib::vFlag[3]

protected

Definition at line 41 of file DedxCalib.h.

Referenced by DedxCalibCostheta::FillHists(), DedxCalibMomentum::FillHists(), DedxCalibTzero::FillHists(), getCurvePar(), and set_dEdx().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DedxCalib()

◆ ~DedxCalib()

Member Function Documentation

◆ AnalyseHists()

◆ BookHists()

◆ cal_dedx()

◆ cal_dedx_bitrunc()

◆ execute()

◆ FillHists()

◆ finalize()

◆ genNtuple()

◆ getCurvePar()

◆ initialize()

◆ initializing()

◆ ReadRecFileList()

◆ set_dEdx()

◆ WriteHists()

Member Data Documentation

◆ Curve_Para

◆ exsvc

◆ geosvc

◆ m_calibflag

◆ m_chi_ex

◆ m_curvefile

◆ m_dedx_exp

◆ m_eventType

◆ m_ex_sigma

◆ m_phShapeFlag

◆ m_pid_prob

◆ m_recFileList

◆ m_recFileVector

◆ m_rootfile

◆ ParticleFlag

◆ Sigma_Para

◆ truncate

◆ vFlag