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

◆ ~DedxCalib()

DedxCalib::~DedxCalib ( )

inline

Definition at line 15 of file DedxCalib.h.

15{};

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 85 of file DedxCalib.cxx.

85{}

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 81 of file DedxCalib.cxx.

81{}

Referenced by initialize().

◆ cal_dedx()

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

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

◆ execute()

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

◆ FillHists()

void DedxCalib::FillHists ( )

pure virtual

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

Definition at line 83 of file DedxCalib.cxx.

83{}

Referenced by execute().

◆ finalize()

StatusCode DedxCalib::finalize ( )

Definition at line 69 of file DedxCalib.cxx.

{
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "DedxCalib finalize()" << endreq;
 
    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 79 of file DedxCalib.cxx.

79{}

Referenced by execute().

◆ getCurvePar()

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

◆ initialize()

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

◆ initializing()

void DedxCalib::initializing ( )

pure virtual

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

Definition at line 77 of file DedxCalib.cxx.

77{}

Referenced by initialize().

◆ ReadRecFileList()

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

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

◆ WriteHists()

void DedxCalib::WriteHists ( )

pure virtual

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

Definition at line 87 of file DedxCalib.cxx.

87{}

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 30 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 29 of file DedxCalib.h.

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

◆ m_calibflag

int DedxCalib::m_calibflag

protected

Definition at line 51 of file DedxCalib.h.

Referenced by DedxCalib(), and initialize().

◆ m_chi_ex

double DedxCalib::m_chi_ex[5]

protected

Definition at line 46 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 56 of file DedxCalib.h.

Referenced by DedxCalib(), and getCurvePar().

◆ m_dedx_exp

double DedxCalib::m_dedx_exp[5]

protected

Definition at line 43 of file DedxCalib.h.

Referenced by set_dEdx().

◆ m_eventType

std::string DedxCalib::m_eventType

protected

Definition at line 53 of file DedxCalib.h.

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

◆ m_ex_sigma

double DedxCalib::m_ex_sigma[5]

protected

Definition at line 44 of file DedxCalib.h.

Referenced by set_dEdx().

◆ m_phShapeFlag

int DedxCalib::m_phShapeFlag

protected

Definition at line 52 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 45 of file DedxCalib.h.

Referenced by set_dEdx().

◆ m_recFileList

std::string DedxCalib::m_recFileList

protected

Definition at line 54 of file DedxCalib.h.

Referenced by DedxCalib(), and ReadRecFileList().

◆ m_recFileVector

vector<string> DedxCalib::m_recFileVector

protected

Definition at line 48 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 55 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 50 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 33 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