#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

	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 ()

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 InstallArea/include/DedxCalibAlg/DedxCalibAlg/DedxCalib.h.

Constructor & Destructor Documentation

◆ DedxCalib() [1/2]

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

◆ ~DedxCalib() [1/2]

DedxCalib::~DedxCalib ( )

inline

Definition at line 15 of file InstallArea/include/DedxCalibAlg/DedxCalibAlg/DedxCalib.h.

15{};

◆ DedxCalib() [2/2]

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

◆ ~DedxCalib() [2/2]

DedxCalib::~DedxCalib ( )

inline

Definition at line 15 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-15/DedxCalibAlg/DedxCalib.h.

15{};

Member Function Documentation

◆ AnalyseHists() [1/2]

void DedxCalib::AnalyseHists ( )

pure virtual

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

Definition at line 85 of file DedxCalib.cxx.

85{}

Referenced by execute().

◆ AnalyseHists() [2/2]

virtual void DedxCalib::AnalyseHists ( )

pure virtual

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

◆ BookHists() [1/2]

void DedxCalib::BookHists ( )

pure virtual

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

Definition at line 81 of file DedxCalib.cxx.

81{}

Referenced by initialize().

◆ BookHists() [2/2]

virtual void DedxCalib::BookHists ( )

pure virtual

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

◆ cal_dedx() [1/2]

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

protected

Definition at line 127 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() [2/2]

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

protected

◆ cal_dedx_bitrunc() [1/2]

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

protected

Definition at line 108 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(), and DedxCalibRunByRun::FillHists().

◆ cal_dedx_bitrunc() [2/2]

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

protected

◆ execute() [1/2]

StatusCode DedxCalib::execute ( )

Definition at line 56 of file DedxCalib.cxx.

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

◆ execute() [2/2]

StatusCode DedxCalib::execute ( )

◆ FillHists() [1/2]

void DedxCalib::FillHists ( )

pure virtual

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

Definition at line 83 of file DedxCalib.cxx.

83{}

Referenced by execute().

◆ FillHists() [2/2]

virtual void DedxCalib::FillHists ( )

pure virtual

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

◆ finalize() [1/2]

StatusCode DedxCalib::finalize ( )

Definition at line 69 of file DedxCalib.cxx.

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

◆ finalize() [2/2]

StatusCode DedxCalib::finalize ( )

◆ genNtuple() [1/2]

void DedxCalib::genNtuple ( )

pure virtual

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

Definition at line 79 of file DedxCalib.cxx.

79{}

Referenced by execute().

◆ genNtuple() [2/2]

virtual void DedxCalib::genNtuple ( )

pure virtual

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

◆ getCurvePar() [1/2]

void DedxCalib::getCurvePar ( )

protected

Definition at line 146 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().

◆ getCurvePar() [2/2]

void DedxCalib::getCurvePar ( )

protected

◆ initialize() [1/2]

StatusCode DedxCalib::initialize ( )

Definition at line 24 of file DedxCalib.cxx.

{
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "DedxCalib initialze()" << endreq;
 
    StatusCode gesc = service("MdcGeomSvc", geosvc);
    if (gesc ==  StatusCode::SUCCESS)
    {
        log << MSG::INFO <<"MdcGeomSvc is running"<<endreq;
    }
    else
    {
        log << MSG::ERROR <<"MdcGeomSvc is failed"<<endreq;
        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;
}

◆ initialize() [2/2]

StatusCode DedxCalib::initialize ( )

◆ initializing() [1/2]

void DedxCalib::initializing ( )

pure virtual

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

Definition at line 77 of file DedxCalib.cxx.

77{}

Referenced by initialize().

◆ initializing() [2/2]

virtual void DedxCalib::initializing ( )

pure virtual

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

◆ ReadRecFileList() [1/2]

void DedxCalib::ReadRecFileList ( )

protected

Definition at line 89 of file DedxCalib.cxx.

{   
 
    MsgStream log(msgSvc(), name());
    log<<MSG::INFO<<"DedxCalib::ReadRecFileList()"<<endreq;
 
    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(), and DedxCalibWireGain::BookHists().

◆ ReadRecFileList() [2/2]

void DedxCalib::ReadRecFileList ( )

protected

◆ set_dEdx() [1/2]

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 184 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(), and DedxCalibMomentum::FillHists().

◆ set_dEdx() [2/2]

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
	)