HadronPrep Class Reference

Class to prepare sample for hadron saturation calibration. More...

#include <HadronPrep.h>

Collaboration diagram for HadronPrep:

Public Member Functions
	HadronPrep ()
	Constructor: Sets the description, the properties and the parameters of the algorithm.
virtual	~HadronPrep ()
	Destructor.
	HadronPrep (int bgbins, double upperbg, double lowerbg, int cosbins, double uppercos, double lowercos, double cut)
	Constructor: set the input variables.
void	prepareSample (std::shared_ptr< TTree > hadron, TFile *&outfile, const std::string &suffix, const std::string &pdg, bool ismakePlots, bool correct)
	function to prepare sample for monitoring plots, bg curve fitting and sigma vs ionz fitting
void	defineHisto (std::vector< TH1F * > &htemp, const std::string &title, const std::string &pdg)
	function to define histograms
void	plotDist (std::map< int, std::vector< TH1F * > > &hist, const std::string &sname, const std::string &pdg)
	function to plot the map of histograms
void	plotDist (std::vector< TH2F * > &hist, const std::string &sname, const std::string &pdg)
	function to plot the 2-D histograms
void	setPars (TFile *&outfile, std::map< int, std::vector< TH1F * > > &hdedx_bgcosth, const std::string &pdg)
	function to fill the parameters like mean and reso in the tree
void	plotGraph (const std::string &sname, const std::string &pdg)
	function to make graph dedx vs cos in different bg bins
void	clear ()
	function to clear the variables

Private Attributes
std::map< int, std::vector< double > >	m_sumcos
	variables to add cos values
std::map< int, std::vector< double > >	m_sumbg
	variables to add bg values
std::map< int, std::vector< double > >	m_means
	mean variable
std::map< int, std::vector< double > >	m_errors
	error variable
std::map< int, std::vector< int > >	m_sumsize
	variables for size
double	m_dedxmax = 0.0
	variables to set maximum dedx mean
double	m_dedxmin = 99999999.0
	variables to set minimum dedx mean
int	m_bgBins
	bins for dedx histogram
double	m_bgMin
	min range of dedx
double	m_bgMax
	max range of dedx
int	m_cosBins
	bins for dedx histogram
double	m_cosMin
	min range of dedx
double	m_cosMax
	max range of dedx
double	m_cut
	cut to clean protons

Detailed Description

Class to prepare sample for hadron saturation calibration.

Definition at line 40 of file HadronPrep.h.

Constructor & Destructor Documentation

HadronPrep() [1/2]

HadronPrep

(

)

Constructor: Sets the description, the properties and the parameters of the algorithm.

Definition at line 12 of file HadronPrep.cc.

{
  m_bgBins = 10;
  m_bgMin = 2.85;
  m_bgMax = 40;
  m_cosBins = 8;
  m_cosMin = 0.0;
  m_cosMax = 0.95;
  m_cut = 0.5;
}

~HadronPrep()

virtual ~HadronPrep ( )

inlinevirtual

Destructor.

Definition at line 52 of file HadronPrep.h.

{};

HadronPrep() [2/2]

HadronPrep	(	int	bgbins,
		double	upperbg,
		double	lowerbg,
		int	cosbins,
		double	uppercos,
		double	lowercos,
		double	cut
	)

Constructor: set the input variables.

Definition at line 23 of file HadronPrep.cc.

{
  m_bgBins = bgbins;
  m_bgMin = lowerbg;
  m_bgMax = upperbg;
  m_cosBins = cosbins;
  m_cosMin = lowercos;
  m_cosMax = uppercos;
  m_cut = cut;
}

Member Function Documentation

clear()

void clear

(

)

function to clear the variables

Definition at line 367 of file HadronPrep.cc.

{
  m_sumcos.clear();
  m_sumbg.clear();
  m_sumsize.clear();
  m_means.clear();
  m_errors.clear();
}

defineHisto()

void defineHisto	(	std::vector< TH1F * > &	htemp,
		const std::string &	title,
		const std::string &	pdg
	)

function to define histograms

Definition at line 157 of file HadronPrep.cc.

{
  for (int j = 0; j < m_cosBins; ++j) {
    // initialize the histograms
    std::string title = Form("%s_%s_%d", pdg.data(), svar.data(), j);
    double dedxMax = 4.0;
    int cosbins = int(50 * dedxMax);
    if (pdg == "pion")  {
      dedxMax = 2.0; cosbins = int(100 * dedxMax);
    }
    if (pdg == "proton")  dedxMax = 20.0;
    if (svar.substr(0, svar.find("_")) == "dedx")
      htemp.push_back(new TH1F(title.data(), title.data(), cosbins, 0, dedxMax));
    else
      htemp.push_back(new TH1F(title.data(), title.data(), cosbins, -30, 30));
 
  }
}

plotDist() [1/2]

void plotDist	(	std::map< int, std::vector< TH1F * > > &	hist,
		const std::string &	sname,
		const std::string &	pdg
	)

function to plot the map of histograms

Definition at line 177 of file HadronPrep.cc.

{
 
  TCanvas ctmp(Form("cdcdedx_%s_%s", sname.data(), pdg.data()), "", 1200, 1200);
  ctmp.Divide(2, 2);
  ctmp.SetBatch(kTRUE);
 
  std::stringstream psname;
  psname << Form("plots/HadronSat/%s_%s.pdf[", sname.data(), pdg.data());
  ctmp.Print(psname.str().c_str());
  psname.str("");
  psname << Form("plots/HadronSat/%s_%s.pdf", sname.data(), pdg.data());
 
  for (int i = 0 ; i < m_bgBins; ++i) {
    for (int j = 0; j < m_cosBins; ++j) {
      hist[i][j]->SetFillColor(kYellow - 9);
 
      ctmp.cd(j % 4 + 1);
      hist[i][j]->Draw();
 
      if (m_cosBins <= 4) {
        if (j + 1 == m_cosBins) ctmp.Print(psname.str().c_str());
      } else if ((j + 1) % 4 == 0 || (j + 1 == m_cosBins)) {
        ctmp.Print(psname.str().c_str());
        ctmp.Clear("D");
      }
    }
  }
  psname.str("");
  psname << Form("plots/HadronSat/%s_%s.pdf]", sname.data(), pdg.data());
  ctmp.Print(psname.str().c_str());
}

plotDist() [2/2]

void plotDist	(	std::vector< TH2F * > &	hist,
		const std::string &	sname,
		const std::string &	pdg
	)

function to plot the 2-D histograms

Definition at line 211 of file HadronPrep.cc.

{
 
  TCanvas ctmp(Form("cdcdedx_%s_%s", sname.data(), pdg.data()), "", 1200, 1200);
  ctmp.Divide(2, 2);
  ctmp.SetBatch(kTRUE);
 
  std::stringstream psname;
  psname << Form("plots/HadronSat/%s_%s.pdf[", sname.data(), pdg.data());
  ctmp.Print(psname.str().c_str());
  psname.str("");
  psname << Form("plots/HadronSat/%s_%s.pdf", sname.data(), pdg.data());
  double bgstep = (m_bgMax - m_bgMin) / m_bgBins;
  double dmax = 2.5;
  if (pdg == "kaon")dmax = 3.5;
  else if (pdg == "proton")dmax = 15.5;
 
  for (int i = 0 ; i < m_bgBins; ++i) {
    ctmp.cd(i % 4 + 1);
    double frombg = m_bgMin + i * bgstep;
    double tobg = frombg + bgstep;
 
    hist[i]->SetTitle(Form("%s, pF=(%0.02f,%0.02f)", hist[i]->GetTitle(), frombg, tobg));
    hist[i]->GetYaxis()->SetRangeUser(-0.5, dmax);
    hist[i]->DrawCopy("colz");
    if ((i + 1) % 4 == 0 || (i + 1) == m_bgBins) {
      ctmp.Print(psname.str().c_str());
      ctmp.Clear("D");
    }
  }
  psname.str("");
  psname << Form("plots/HadronSat/%s_%s.pdf]", sname.data(), pdg.data());
  ctmp.Print(psname.str().c_str());
}

plotGraph()

void plotGraph	(	const std::string &	sname,
		const std::string &	pdg
	)

function to make graph dedx vs cos in different bg bins

Definition at line 304 of file HadronPrep.cc.

{
  std::vector<double> cbcenters(m_cosBins), cberrors(m_cosBins);
  double cosstep = (m_cosMax - m_cosMin) / m_cosBins;
  double bgstep = (m_bgMax - m_bgMin) / m_bgBins;
 
  for (int j = 0; j < m_cosBins; ++j) {
    cbcenters[j] = m_cosMin + 0.5 * cosstep + j * cosstep;
    cberrors[j] = 0;
  }
 
  // Plot the dE/dx means vs. cos(theta) for validation
  TCanvas ctmp6(Form("cdcdedx_costh_%s_%s", sname.data(), pdg.data()), "", 500, 500);
  ctmp6.SetGridy(1);
 
  TH1F base("base", "bla-bla", 25, 0, 1.0);
  base.SetTitle(Form("%s: dE/dx fit means vs. cos(#theta);cos(#theta);dE/dx-fit mean", pdg.data()));
  base.GetXaxis()->SetTitleOffset(1.2);
 
  double margin = 0.1;
  if (pdg == "proton") margin = 0.95;
 
  base.SetMaximum(m_dedxmax + margin);
  base.SetMinimum(m_dedxmin - margin);
  base.SetStats(0);
  base.DrawCopy("");
 
  TLegend legend(0.75, 0.75, 0.95, 0.9);
 
  std::vector<TGraphErrors*> gdedx_costh(m_bgBins);
 
  legend.SetBorderSize(0);
 
  for (int i = 0; i < m_bgBins; ++i) {
    std::vector<double> mean_d(m_cosBins);
    std::vector<double> error_d(m_cosBins);
 
    for (int j = 0; j < m_cosBins; ++j) {
      mean_d [j] = m_means[i][j];
      error_d[j] = m_errors[i][j];
    }
    gdedx_costh[i] = new TGraphErrors(m_cosBins, cbcenters.data(), mean_d.data(), cberrors.data(), error_d.data());
 
    gdedx_costh[i]->SetMarkerSize(0.9);
    gdedx_costh[i]->SetMarkerColor(50 + i * 3);
    gdedx_costh[i]->SetLineColor(i + 1);
    gdedx_costh[i]->SetLineWidth(1);
    gdedx_costh[i]->Draw("same");
    legend.AddEntry(gdedx_costh[i], Form("bg = (%0.03f - %0.03f)", m_bgMin + i * bgstep, m_bgMin + (i + 1) * bgstep), "lep");
  }
 
  legend.Draw("same");
 
  std::stringstream psname;
  psname.str("");
  psname << Form("plots/HadronSat/gr_dedx_vscosth_%s_%s.pdf", sname.data(), pdg.data());
  ctmp6.SaveAs(psname.str().c_str());
 
  for (int i = 0; i < m_bgBins; ++i)
    delete gdedx_costh[i];
}

prepareSample()

void prepareSample	(	std::shared_ptr< TTree >	hadron,
		TFile *&	outfile,
		const std::string &	suffix,
		const std::string &	pdg,
		bool	ismakePlots,
		bool	correct
	)

function to prepare sample for monitoring plots, bg curve fitting and sigma vs ionz fitting

Definition at line 34 of file HadronPrep.cc.

{
 
  std::map<int, std::vector<TH1F*>>  hdedx_bgcosth;
  std::vector<TH2F*> hdedxvscosth_bg;
 
  //define histograms
  for (int i = 0; i < m_bgBins; ++i) {
 
    std::string title = Form("%s_dedxvscos_bg_%d", pdg.data(), i);
    hdedxvscosth_bg.push_back(new TH2F(title.data(), Form("%s;costh;dEdx", title.data()),
                                       440, -1.1, 1.1, 2600, -1.0, 25));
 
    defineHisto(hdedx_bgcosth[i], Form("dedx_bg_%d_cos", i), pdg.data());
  }
 
 
  // Create some containers to calculate averages
  for (int i = 0; i < m_bgBins; ++i) {
    for (int j = 0; j < m_cosBins; ++j) {
      m_sumcos[i].push_back(0.);
      m_sumbg[i].push_back(0.);
      m_sumsize[i].push_back(0.);
      m_means[i].push_back(0.);
      m_errors[i].push_back(0.);
 
    }
  }
 
  // // --------------------------------------------------
  // // LOOP OVER EVENTS AND FILL CONTAINERS
  // // --------------------------------------------------
  // // Fill the histograms to be fitted
 
  double dedxnosat;    // dE/dx w/o HS correction (use to get new parameters)
  double p;       // track momentum
  double costh;   // cosine of track polar angle
  double timereso;
  int nhit;       // number of hits on this track
 
  hadron->SetBranchAddress("dedxnosat", &dedxnosat); //must be without any HS correction
  hadron->SetBranchAddress("p", &p);
  hadron->SetBranchAddress("costh", &costh);
  hadron->SetBranchAddress("timereso", &timereso);
  hadron->SetBranchAddress("nhits", &nhit);
 
  HadronBgPrep prep;
  double mass = prep.getParticleMass(pdg);
  if (mass == 0.0) B2FATAL("Mass of particle " << pdg.data() << " is zero");
 
  CDCDedxHadSat had;
 
  if (correct) {
    std::string infile = "sat-pars.fit.txt";
    std::ifstream parfile(infile);
 
    if (!parfile.fail()) {
      B2INFO("new parameters are using");
      had.setParameters(infile);
    } else {
      B2INFO("default parameters are using");
      had.setParameters();
    }
  }
 
  for (unsigned int index = 0; index < hadron->GetEntries(); ++index) {
 
    hadron->GetEvent(index);
 
    double bg;      // track beta-gamma
    bg = fabs(p) / mass;
 
    double uncosth = costh;
    costh = fabs(costh);
 
    // clean up bad events and restrict the momentum range
    if (nhit <  0 || nhit  > 100)continue;
    if (dedxnosat <= 0)continue;
    if (costh != costh)continue;
    if (bg <= m_bgMin || bg >= m_bgMax)continue;
    if (costh <= m_cosMin || costh >= m_cosMax)continue;
 
    if (pdg == "proton") if ((dedxnosat - 0.45) * abs(p) * abs(p) < m_cut)continue;
 
    int bgBin = (int)((bg - m_bgMin) / (m_bgMax - m_bgMin) * m_bgBins);
    int cosBin = (int)((costh - m_cosMin) / (m_cosMax - m_cosMin) * m_cosBins);
 
    double dedx_new = dedxnosat;
 
    if (correct) dedx_new = had.D2I(costh, had.I2D(costh, 1.0) * dedxnosat);
 
    hdedx_bgcosth[bgBin][cosBin]->Fill(dedx_new);
    hdedxvscosth_bg[bgBin]->Fill(uncosth, dedx_new);
 
    m_sumcos[bgBin][cosBin] += costh;
    m_sumbg[bgBin][cosBin] += bg;
    m_sumsize[bgBin][cosBin] += 1;
 
  }// end of event loop
 
  // // --------------------------------------------------
  // // FIT IN BINS OF BETA-GAMMA AND COS(THETA)
  // // --------------------------------------------------
  // // fit the histograms with Gaussian functions
  // // and extract the means and errors
 
  setPars(outfile, hdedx_bgcosth, pdg.data());
 
  if (ismakePlots) {
 
    plotDist(hdedx_bgcosth, Form("fits_dedx_inbg_%s", suffix.data()), pdg.data());
    plotDist(hdedxvscosth_bg, Form("dist_dedx_vscos_scat_%s", suffix.data()), pdg.data());
    plotGraph(suffix.data(), pdg.data());
  }
 
  hdedxvscosth_bg.clear();
  hdedx_bgcosth.clear();
  clear();
 
}

setPars()

void setPars	(	TFile *&	outfile,
		std::map< int, std::vector< TH1F * > > &	hdedx_bgcosth,
		const std::string &	pdg
	)

function to fill the parameters like mean and reso in the tree

Definition at line 247 of file HadronPrep.cc.

{
  outfile->cd();
  TTree* satTree = new TTree(Form("%s", pdg.data()), "dE/dx means and errors");
  double satbg;          // beta-gamma value for this bin
  double satcosth;       // cos(theta) value for this bin
  double satdedx;        // mean dE/dx value for this bin
  double satdedxerr;     // error on ^
  double satdedxwidth;// width of ^ distribution
  double satbg_avg;      // average beta-gamma value for this sample
  double satcosth_avg;   // average cos(theta) value for this sample
 
  satTree->Branch("bg", &satbg, "bg/D");
  satTree->Branch("costh", &satcosth, "costh/D");
 
  //dEdx related variables
  satTree->Branch("dedx", &satdedx, "dedx/D");
  satTree->Branch("dedxerr", &satdedxerr, "dedxerr/D");
  satTree->Branch("dedxwidth", &satdedxwidth, "dedxwidth/D");
 
  // Other variables
  satTree->Branch("bg_avg", &satbg_avg, "bg_avg/D");
  satTree->Branch("costh_avg", &satcosth_avg, "costh_avg/D");
 
  double bgstep = (m_bgMax - m_bgMin) / m_bgBins;
  double cosstep = (m_cosMax - m_cosMin) / m_cosBins;
 
  for (int i = 0; i < m_bgBins; ++i) {
    for (int j = 0; j < m_cosBins; ++j) {
 
      // fill some details for this bin
      satbg = m_bgMin + 0.5 * bgstep + i * bgstep;
      satcosth = m_cosMin + 0.5 * cosstep + j * cosstep;
 
      satbg_avg = m_sumbg[i][j] / m_sumsize[i][j];
      satcosth_avg = m_sumcos[i][j] / m_sumsize[i][j];
 
      //1. -------------------------
      // fit the dE/dx distribution in bins of beta-gamma and cosine
      HadronBgPrep prep;
      prep.fit(hdedx_bgcosth[i][j],  pdg.data());
      satdedx = m_means[i][j] = hdedx_bgcosth[i][j]->GetFunction("gaus")->GetParameter(1);
      satdedxerr = m_errors[i][j] = hdedx_bgcosth[i][j]->GetFunction("gaus")->GetParError(1);
      satdedxwidth = hdedx_bgcosth[i][j]->GetFunction("gaus")->GetParameter(2);
 
      //Be careful to not set for every particle
      if (satdedx > m_dedxmax) m_dedxmax = satdedx;
      if (satdedx < m_dedxmin) m_dedxmin = satdedx;
 
      // fill the tree for this bin
      satTree->Fill();
    }
  }
  outfile->Write();
}

Member Data Documentation

m_bgBins

int m_bgBins

private

bins for dedx histogram

Definition at line 107 of file HadronPrep.h.

m_bgMax

double m_bgMax

private

max range of dedx

Definition at line 109 of file HadronPrep.h.

m_bgMin

double m_bgMin

private

min range of dedx

Definition at line 108 of file HadronPrep.h.

m_cosBins

int m_cosBins

private

bins for dedx histogram

Definition at line 111 of file HadronPrep.h.

m_cosMax

double m_cosMax

private

max range of dedx

Definition at line 113 of file HadronPrep.h.

m_cosMin

double m_cosMin

private

min range of dedx

Definition at line 112 of file HadronPrep.h.

m_cut

double m_cut

private

cut to clean protons

Definition at line 115 of file HadronPrep.h.

m_dedxmax

double m_dedxmax = 0.0

private

variables to set maximum dedx mean

Definition at line 104 of file HadronPrep.h.

m_dedxmin

double m_dedxmin = 99999999.0

private

variables to set minimum dedx mean

Definition at line 105 of file HadronPrep.h.

m_errors

std::map<int, std::vector<double> > m_errors

private

error variable

Definition at line 100 of file HadronPrep.h.

m_means

std::map<int, std::vector<double> > m_means

private

mean variable

Definition at line 99 of file HadronPrep.h.

m_sumbg

std::map<int, std::vector<double> > m_sumbg

private

variables to add bg values

Definition at line 98 of file HadronPrep.h.

m_sumcos

std::map<int, std::vector<double> > m_sumcos

private

variables to add cos values

Definition at line 97 of file HadronPrep.h.

m_sumsize

std::map<int, std::vector<int> > m_sumsize

private

variables for size

Definition at line 102 of file HadronPrep.h.

---

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

• cdc/calibration/CDCdEdx/include/HadronPrep.h
• cdc/calibration/CDCdEdx/src/HadronPrep.cc