RresTbin Class Reference

Inheritance diagram for RresTbin:

Collaboration diagram for RresTbin:

Public Member Functions
	RresTbin (TFile *, const string &, int)

Public Attributes
vector< TH1D * >	reshists
map< double, double >	resmap

Private Attributes
TGraphErrors	m_thegraph
std::vector< double >	m_rdata
std::vector< double >	m_erdata
std::vector< double >	m_tdata
std::vector< double >	m_etdata

Detailed Description

Definition at line 3330 of file TRTCalib_makeplots.cxx.

Constructor & Destructor Documentation

RresTbin()

RresTbin::RresTbin	(	TFile *	file,
		const string &	detname,
		int	det
	)

Definition at line 3345 of file TRTCalib_makeplots.cxx.

{
 
    TF1 ff("fitfunc", "gaus");
    TF1 *ff2 = new TF1("fitfunc2", "pol3");
 
    TDirectory *trt = nullptr;
    TDirectory *detdir = nullptr;
    TDirectory *binhist = nullptr;
    string resname;
    bool isAr = false;
    if (detname.find("Ar") != string::npos)
        isAr = true;
    bool allTRT = false;
    if (detname == "WholeBarrel_1")
    {
        resname = "reshists5";
    }
    else if (detname == "Detector_-1")
    {
        resname = "reshists2";
    }
    else if (detname == "Detector_1")
    {
        resname = "reshists3";
    }
    else if (detname == "Detector_-2")
    {
        resname = "reshists1";
    }
    else if (detname == "Detector_2")
    {
        resname = "reshists4";
    }
    else if (!isAr)
    {
        resname = "reshists";
        allTRT = true;
    }
    else if (detname == "WholeBarrel_Ar_1")
    {
        resname = "reshists_Ar5";
    }
    else if (detname == "Detector_Ar_-1")
    {
        resname = "reshists_Ar2";
    }
    else if (detname == "Detector_Ar_1")
    {
        resname = "reshists_Ar3";
    }
    else if (detname == "Detector_Ar_-2")
    {
        resname = "reshists_Ar1";
    }
    else if (detname == "Detector_Ar_2")
    {
        resname = "reshists_Ar4";
    }
    else
    {
        resname = "reshists_Ar";
        allTRT = true;
    }
 
    if (!isAr)
        trt = (TDirectory *)file->FindKey("TRT_all")->ReadObj();
    else
        trt = (TDirectory *)file->FindKey("TRT_Ar_all")->ReadObj();
 
    if (!isAr)
        cout << "In RresTbin Xenon" << endl;
    if (isAr)
        cout << "In RresTbin Argon" << endl;
    if (trt)
    {
        cout << "   found TRT_all directory " << endl;
        detdir = (TDirectory *)trt->FindKey(detname.c_str())->ReadObj();
        if (detdir)
        {
            cout << "   found binhist directory for detector name " << detname << endl;
            if (!allTRT)
            {
                binhist = (TDirectory *)detdir->FindKey(resname.c_str())->ReadObj();
            }
            else
            {
                binhist = (TDirectory *)trt->FindKey(resname.c_str())->ReadObj();
            }
 
            if (binhist)
            {
                cout << "   found reshists directory for detector number: " << det << endl;
            }
            else
            {
                cout << "   did not find reshists directory for detector number " << det << endl;
            }
        }
        else
        {
            cout << "   did not find binhist directory for detector name " << detname << endl;
        }
    }
    else
    {
        cout << "   did not find TRT_all directory " << endl;
    }
 
    if (binhist)
    {
 
        TVectorD *tbins = (TVectorD *)file->FindKey("tbins")->ReadObj();
 
        if (tbins)
        {
            const size_t nElements = tbins->GetNoElements();
            cout << "   found tbins vector with length " << nElements << endl;
            m_rdata = std::vector<double>(nElements,0.);
            m_erdata = std::vector<double>(nElements,0.);
            m_tdata = std::vector<double>(nElements,0.);
            m_etdata = std::vector<double>(nElements,0.);
            m_thegraph = TGraphErrors(nElements);
 
            int ipoint = 0;
            for (int ihist = 0; ihist < tbins->GetNoElements(); ihist++)
            {
 
                TH1D *reshist = nullptr;
 
                if (det == 4)
                    reshist = (TH1D *)binhist->Get(Form("res_tbin%i_bar", ihist));
                else if (det == 5)
                    reshist = (TH1D *)binhist->Get(Form("res_tbin%i_trt", ihist));
                else
                    reshist = (TH1D *)binhist->Get(Form("res_tbin%i_%i", ihist, det));
 
                if (reshist != nullptr)
                {
                    double lolim = -0.1; // reshist->GetBinCenter(reshist->GetMaximumBin()-10);
                    double hilim = 0.1;  // reshist->GetBinCenter(reshist->GetMaximumBin()+10);
 
                    if (reshist->Fit(&ff, "Q", "", lolim, hilim) + 1)
                    {
                        m_rdata[ipoint] = ff.GetParameter(1);
                        m_erdata[ipoint] = ff.GetParError(1);
                        m_tdata[ipoint] = (*tbins)[ihist];
 
                        m_etdata[ipoint] = 0;
                        // printf("%f %f\n",ff.GetParameter(1),ff.GetParError(1));
                        ipoint++;
                    }
 
                    reshists.push_back(reshist);
                }
            }
            cout << "        succesfully fitted " << ipoint << " histograms " << endl;
            for (int ip = 0; ip < ipoint; ip++)
            {
                // cout << rdata[ip]  << " ";
                resmap[m_tdata[ip]] = m_rdata[ip];
            }
            // cout << endl;
 
            if (ipoint > 10)
            {
                m_thegraph = TGraphErrors(ipoint, m_tdata.data(), m_rdata.data(), m_etdata.data(), m_erdata.data());
 
                cout << "       draw t-binned residual graph " << endl;
                string detstr[8] = {"endcap A", "barrel A", "", "barrel C", "endcap C", "", "whole barrel", "whole TRT"};
 
                m_thegraph.SetMarkerStyle(20);
                m_thegraph.SetTitle(Form("t-binned absolute residual (%s)", detstr[det + 2].c_str()));
                m_thegraph.Draw("ap");
                m_thegraph.GetYaxis()->SetTitle("|r|-|r_{track}| /mm");
                m_thegraph.GetYaxis()->SetRangeUser(-0.150, 0.150);
                m_thegraph.GetXaxis()->SetRangeUser(0, 50);
                m_thegraph.GetXaxis()->SetTitle("t-t0 /ns");
 
                printf("RPOLDIFF %i %e %e %e %e\n", det, ff2->GetParameter(0), ff2->GetParameter(1), ff2->GetParameter(2), ff2->GetParameter(3));
 
                this->SetGrid();
            }
        }
    }
}

Member Data Documentation

m_erdata

std::vector<double> RresTbin::m_erdata

private

Definition at line 3340 of file TRTCalib_makeplots.cxx.

m_etdata

std::vector<double> RresTbin::m_etdata

private

Definition at line 3342 of file TRTCalib_makeplots.cxx.

m_rdata

std::vector<double> RresTbin::m_rdata

private

Definition at line 3339 of file TRTCalib_makeplots.cxx.

m_tdata

std::vector<double> RresTbin::m_tdata

private

Definition at line 3341 of file TRTCalib_makeplots.cxx.

m_thegraph

TGraphErrors RresTbin::m_thegraph

private

Definition at line 3338 of file TRTCalib_makeplots.cxx.

reshists

vector<TH1D *> RresTbin::reshists

Definition at line 3334 of file TRTCalib_makeplots.cxx.

resmap

map<double, double> RresTbin::resmap

Definition at line 3335 of file TRTCalib_makeplots.cxx.

---

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

• TRTCalib_makeplots.cxx