#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-20/DedxCalibAlg/DedxCalib.h.

15{};

Member Function Documentation

◆ AnalyseHists() [1/2]

void DedxCalib::AnalyseHists ( )

pure virtual

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

Definition at line 85 of file DedxCalib.cxx.

85{}

Referenced by execute().

◆ AnalyseHists() [2/2]

virtual void DedxCalib::AnalyseHists ( )

pure virtual

◆ BookHists() [1/2]

void DedxCalib::BookHists ( )

pure virtual

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

Definition at line 81 of file DedxCalib.cxx.

81{}

Referenced by initialize().

◆ BookHists() [2/2]

virtual void DedxCalib::BookHists ( )

pure virtual

◆ 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(), DedxCalibRunByRun::FillHists(), and DedxCalibTzero::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, DedxCalibTzero, DedxCalibWireGain, DedxSaveWireGain, DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 83 of file DedxCalib.cxx.

83{}

Referenced by execute().

◆ FillHists() [2/2]

virtual void DedxCalib::FillHists ( )

pure virtual

◆ 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, DedxCalibTzero, DedxCalibWireGain, DedxSaveWireGain, DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 79 of file DedxCalib.cxx.

79{}

Referenced by execute().

◆ genNtuple() [2/2]

virtual void DedxCalib::genNtuple ( )

pure virtual

◆ 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, DedxCalibTzero, DedxCalibWireGain, DedxSaveWireGain, DedxCalibCostheta, DedxCalibDocaEAng, DedxCalibEAng, DedxCalibEvent, DedxCalibLayerGain, DedxCalibMomentum, DedxCalibRunByRun, DedxCalibTzero, DedxCalibWireGain, and DedxSaveWireGain.

Definition at line 77 of file DedxCalib.cxx.

77{}

Referenced by initialize().

◆ initializing() [2/2]

virtual void DedxCalib::initializing ( )

pure virtual

◆ 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(), DedxCalibTzero::BookHists(), DedxCalibWireGain::BookHists(), and DedxSaveWireGain::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(), DedxCalibMomentum::FillHists(), and DedxCalibTzero::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
	)