#include "TCanvas.h"
#include "TFile.h"
#include "TTree.h"
#include "TH1F.h"
#include <sstream>
#include <iostream>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include "curve.h"

Macros
#define	mass_e 0.511e-3

#define	mass_mu 105.658e-3

#define	mass_pi 139.570e-3

#define	mass_k 493.677e-3

#define	mass_p 938.272e-3

#define	P_electron_min 421.97*mass_e

#define	P_electron_max 3030*mass_e

#define	P_electron_bin 80

#define	P_muon_min 10.98001*mass_mu

#define	P_muon_max 115.6*mass_mu

#define	P_muon_bin 20

#define	P_pion_min 0.8

#define	P_pion_max 0.805

#define	P_pion_bin 1

#define	P_kaon_min 0.63001*mass_k

#define	P_kaon_max 0.74*mass_k

#define	P_kaon_bin 5

#define	P_proton_min 0.26*mass_p

#define	P_proton_max 0.63*mass_p

#define	P_proton_bin 320

Functions
void	histgen (string, int, TObjArray &, vector< double > &)

void	histps (string)

bool	event_validate (double, double, double, int)

double	get_hist_range (int, double, int)

void	hist_update ()

const bool	ps_flag (false)

int	main ()

Variables
const int	na = 10

const int	nc = 2

const double	electron_step = (double)(P_electron_max - P_electron_min) / (double)P_electron_bin

const double	muon_step = (double)(P_muon_max - P_muon_min) / (double)P_muon_bin

const double	pion_step = (double)(P_pion_max - P_pion_min) / (double)P_pion_bin

const double	kaon_step = (double)(P_kaon_max - P_kaon_min) / (double)P_kaon_bin

const double	proton_step = (double)(P_proton_max - P_proton_min) / (double)P_proton_bin

const int	run_flag [5] = {0,0,1,0,0}

Macro Definition Documentation

◆ mass_e

#define mass_e 0.511e-3

Definition at line 17 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ mass_k

#define mass_k 493.677e-3

Definition at line 20 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ mass_mu

#define mass_mu 105.658e-3

Definition at line 18 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ mass_p

#define mass_p 938.272e-3

Definition at line 21 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ mass_pi

#define mass_pi 139.570e-3

Definition at line 19 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_electron_bin

#define P_electron_bin 80

Definition at line 25 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_electron_max

#define P_electron_max 3030*mass_e

Definition at line 24 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_electron_min

#define P_electron_min 421.97*mass_e

Definition at line 23 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_kaon_bin

#define P_kaon_bin 5

Definition at line 37 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_kaon_max

#define P_kaon_max 0.74*mass_k

Definition at line 36 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_kaon_min

#define P_kaon_min 0.63001*mass_k

Definition at line 35 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_muon_bin

#define P_muon_bin 20

Definition at line 29 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_muon_max

#define P_muon_max 115.6*mass_mu

Definition at line 28 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_muon_min

#define P_muon_min 10.98001*mass_mu

Definition at line 27 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_pion_bin

#define P_pion_bin 1

Definition at line 33 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_pion_max

#define P_pion_max 0.805

Definition at line 32 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_pion_min

#define P_pion_min 0.8

Definition at line 31 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_proton_bin

#define P_proton_bin 320

Definition at line 41 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_proton_max

#define P_proton_max 0.63*mass_p

Definition at line 40 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ P_proton_min

#define P_proton_min 0.26*mass_p

Definition at line 39 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Function Documentation

◆ event_validate()

bool event_validate	(	double	bitrunc,
		double	ptrk,
		double	costheta,
		int	type
	)

Definition at line 348 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

{
    bool validate = false;
 
    switch (type) {
        case 0:
            //if (bitrunc > 420 / pow(ptrk, 1.2))
      if( ( (ptrk<=0.3 && bitrunc>18842.7-102517*ptrk+157204*ptrk*ptrk) || (ptrk>0.3 && bitrunc>5186.13-14056 *ptrk +11221*ptrk*ptrk) ) && (bitrunc<32695.9-136724*ptrk+217506*ptrk*ptrk-120748*ptrk*ptrk*ptrk) )
            validate = true;
            break;
        case 1:
            /*if (bitrunc > 60 / pow(ptrk, 1.8) + 300 &&
              bitrunc < 60 / pow(ptrk, 2.4) + 600)*/
            validate = true;
            break;
        case 2:
            //if (bitrunc < 15 / ptrk / ptrk + 540)
              validate = true;
            break;
        case 3:
            //if (bitrunc < 2000)
            validate = true;
            break;
        case 4:
            if (bitrunc > 0 && bitrunc < 1200)
            validate = true;
            break;
        case 5:
            if (bitrunc < 700 && bitrunc > 400)
              validate = true;
            break;
        case 6:
            validate = true;
            break;
        default:
            std::cout << "Error: no such sample type!" << std::endl;
            break;
    }
 
    if (ptrk <= 0) validate = false;
    return validate;
}

Referenced by histgen().

◆ get_hist_range()

double get_hist_range	(	int	type,
		double	momentum,
		int	lowhigh
	)

Definition at line 391 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

                                                             {
  double mean_exp[5];
  double sigma_exp[5];
 
  const double theta(1), Nohit(15), t0(1200);
  dedx_pid_exp(momentum, theta, Nohit, t0, mean_exp, sigma_exp);
  
  switch (type) {
  case 0:
    //    cout << "type " << type << "  mean " << mean_exp[4] << "  sigma " << sigma_exp[4] << endl;
    if(lowhigh) return mean_exp[4]+22*(sigma_exp[4]);
    else return mean_exp[4]-10*(sigma_exp[4]) ? mean_exp[4]-10*(sigma_exp[4]) : 0;
    break;
  case 1:
    //    cout << "type " << type << " mean " << mean_exp[3] << "  sigma " << sigma_exp[3] << endl;
    if(lowhigh) return mean_exp[3]+26*(sigma_exp[3]);
    else return mean_exp[3]-10*(sigma_exp[3]) ? mean_exp[3]-10*(sigma_exp[3]) : 0;
    break;
  case 2:
    cout << "type " << type << " mean " << mean_exp[2] << "  sigma " << sigma_exp[2] << endl;
    cout << "mean - 10*sigma" << mean_exp[2]-10*(sigma_exp[2]) << endl;
    cout << "mean + 15*sigma" << mean_exp[2]+15*(sigma_exp[2]) << endl;
    cout << "mean + 25*sigma" << mean_exp[2]+25*(sigma_exp[2]) << endl;
    cout << "mean + 28*sigma" << mean_exp[2]+28*(sigma_exp[2]) << endl;
    if(lowhigh) return mean_exp[2]+28*(sigma_exp[2]);
    else return mean_exp[2]-10*(sigma_exp[2]) ? mean_exp[2]-10*(sigma_exp[2]) : 0;
    break;
  case 3:
    //    cout << "type " << type << " mean " << mean_exp[1] << "  sigma " << sigma_exp[1] << endl;
    if(lowhigh) return mean_exp[1]+25*(sigma_exp[1]);
    else return mean_exp[1]-10*(sigma_exp[1]) ? mean_exp[1]-10*(sigma_exp[1]) : 0;
    break;
  case 4:
    //    cout << "type " << type << " mean " << mean_exp[0] << "  sigma " << sigma_exp[0] << endl;
    if(lowhigh) return mean_exp[0]+25*(sigma_exp[0]);
    else return mean_exp[0]-10*(sigma_exp[0]) ? mean_exp[0]-10*(sigma_exp[0]) : 0;
    break;
  default:
    cout << "no type is determined" << endl;
    return mean_exp[2] * 8;
  }
}

Referenced by histgen().

◆ hist_update()

void hist_update ( )

Definition at line 435 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

                  {
  char hname[200];
  TH1F *h1 = new TH1F();
  TTree *tree = new TTree("TH1F_Col", "TH1F_Col");
  tree->Branch("TH1F_Col", "TH1F", &h1, 3200, 0);
  TFile* f1 = new TFile("DedxSampleHists.root", "UPDATE");
  TTree* bin = (TTree*)f1->Get("bin");
  Int_t totNum;
  bin->SetBranchAddress("totalNum", &totNum);
  bin->GetEntry(0);
 
  for (Int_t i = 0; i < totNum; i++){
    if(i >= 10000) sprintf(hname, "h%05d", i);
    if(i >= 1000)  sprintf(hname, "h%04d", i);
    if(i >= 100)   sprintf(hname, "h%03d", i);
    if(i >= 10)    sprintf(hname, "h%02d", i);
    if(i < 10)     sprintf(hname, "h%01d", i);
    h1 = (TH1F*)f1->Get(hname);
    tree->Fill();
  }
  tree->Write();
  f1->Close();
}

Referenced by main().

◆ histgen()

void histgen	(	string	filename,
		int	type,
		TObjArray &	hist,
		vector< double > &	bg
	)

Definition at line 112 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

                                                                            {
  if(!(run_flag[type])){cout << "skip this type " << type << " particle" << endl; return ;}
  double exraw[100];
  float  ptrk, dEdx_meas, costheta, charge;
  float  chip, chipi, chik, chimu, chie;
  double pmin, pmax, pstep, mass;
  int    nbin, ndedxhit;
 
  TFile* f = new TFile(filename.c_str());
  TTree* trk = (TTree*)f->Get("n103");
  trk->SetBranchAddress("dEdx_meas",  &dEdx_meas);
  trk->SetBranchAddress("ptrk",     &ptrk);
  trk->SetBranchAddress("costheta", &costheta);
  trk->SetBranchAddress("charge",    &charge);
  trk->SetBranchAddress("dEdx_hit",     exraw);
  trk->SetBranchAddress("ndedxhit", &ndedxhit);
  trk->SetBranchAddress("chidedx_p", &chip);
  trk->SetBranchAddress("chidedx_mu", &chimu);
  trk->SetBranchAddress("chidedx_pi", &chipi);
  trk->SetBranchAddress("chidedx_k",  &chik);
  trk->SetBranchAddress("chidedx_e",  &chie);
 
  //initialize variables
  switch (type) {
  case 0:
    pmin  = P_proton_min;
    pmax  = P_proton_max;
    pstep = proton_step;
    nbin  = P_proton_bin;
    mass  = mass_p;
    break;
  case 1:
    pmin  = P_kaon_min;
    pmax  = P_kaon_max;
    pstep = kaon_step;
    nbin  = P_kaon_bin;
    mass  = mass_k;
    break;
  case 2:
    pmin  = P_pion_min;
    pmax  = P_pion_max;
    pstep = pion_step;
    nbin  = P_pion_bin;
    mass  = mass_pi;
    break;
  case 3:
    pmin  = P_muon_min;
    pmax  = P_muon_max;
    pstep = muon_step;
    nbin  = P_muon_bin;
    mass  = mass_mu;
    break;
  case 4:
    pmin  = P_electron_min;
    pmax  = P_electron_max;
    pstep = electron_step;
    nbin  = P_electron_bin;
    mass  = mass_e;
    break;
  default:
    std::cout << "Error: no such sample type!" << std::endl;
  }
 
  //define histograms
  TH1F** h = new TH1F*[nbin * 2 * na];
  stringstream ss;
  double momentum;
  for (int i = 0; i < nbin * 2 * na; i++){ // nbin is determined by pmin, pmax and pstep
    ss.str("");
    ss << "h" << histNum + i;
    momentum = pmin + pstep * ((int)(i / 2 / na) + 0.5);
    h[i] = new TH1F(ss.str().c_str(), ss.str().c_str(), 200, get_hist_range(type, momentum, 0), 
            get_hist_range(type, momentum, 1));
  }
  histNum += nbin * 2 * na;
 
  //loop every hit of n103
  int binidx = -1;
  int angleidx = -1;
  int chargeidx = -1;
  int histidx = -1;
  //  int histidx1 = -1;
    
  for (int i = 0; i < trk->GetEntries(); i++){
    trk->GetEntry(i);
    if (type == 0 && chip > 20) continue;
    //if (type == 1 && (chik < -4 || chik > 7)) continue;
    //if (type == 2 && chipi > 7) continue;
    //        if (type == 3 && ptrk < 10 && chimu > 10 ) continue;
    if (type == 3 && ptrk < 10 && dEdx_meas > 600 ) continue;
    if (!event_validate(dEdx_meas, ptrk, costheta, type)) continue;
 
    // betagamma bin index
    binidx = (int)((ptrk - pmin) / pstep);
        
    // costheta bin index
    if (fabs(costheta) >= 0.93)
      angleidx = na - 1;
    else
      angleidx = (int)(fabs(costheta) * na / 0.93);
        
    // charge bin index
    if (charge == 0) 
      std::cout << "Warning: Neutral track exist!" << std::endl;
    else
      chargeidx = (charge > 0) ? 0 : 1;
        
    // fill histograms
    histidx = -1;
    for (int j = 0; j < ndedxhit; j++) {
      if (binidx >= 0 && binidx < nbin && angleidx >= 0 && angleidx < na){
    bool isLowMomentumProton = (type == 0 && binidx< 100);
    bool isLowMomentumKaon = (type == 1 && binidx <= 40 );
    bool isHighMomentumKaon = false; //(type == 1 && binidx >= 40);
    bool isLowMomentumPion = (type == 2 && binidx < 90);
    bool isHighMomentumPion = (type == 2 && binidx >= 956);
    bool isMiddleMomentumGG4pi = false; //(type == 6 && ptrk > 0.75 && ptrk < 0.9);
    bool isHighMomentumGG4pi = false; //(type == 6 && ptrk > 0.9 && ptrk < 1.2);
    bool isHighMomentumElectron = (type == 4 && ptrk > 3.0/2);
    if (isLowMomentumProton) {
      for(int idx_theta =(int)(angleidx/10)*10; idx_theta<(int)(angleidx/10)*10+10; idx_theta++){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge<(int)(chargeidx/2)*2+2; idx_charge++){
          histidx = binidx*2*na + idx_theta*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        } 
      } 
    }// end of low momentum proton
    else if(isLowMomentumKaon){
      if(angleidx<na-2){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
          histidx = binidx*2*na + angleidx*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        }
      }
      else{
        for(int idx_theta=(int)(angleidx/10)*10; idx_theta < (int)(angleidx/10)*10+10; idx_theta++){ 
          for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
        histidx = binidx*2*na + idx_theta*2 + idx_charge;
        h[histidx]->Fill(exraw[j]);
          } // for
        } // for
      }// else
    }// end of low mometnum kaon
    else if(isHighMomentumKaon){
      for(int idx_bg = (int)(binidx/2)*2; idx_bg < (int)(binidx/2)*2+2; idx_bg++){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
          histidx = idx_bg*2*na + angleidx*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        } // for
        } // for
    } // end of low momentum pion
    else if(isLowMomentumPion){
      for(int idx_theta = (int)(angleidx/10)*10; idx_theta < (int)(angleidx/10)*10+10; idx_theta++){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
          histidx = binidx*2*na + idx_theta*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        } // for
      } // for
    } // end of low momentum pion
    else if(isHighMomentumPion){
      for(int idx_bg = (int)(binidx/2)*2; idx_bg < (int)(binidx/2)*2+2; idx_bg++){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
          histidx = idx_bg*2*na + angleidx*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        } // for
      } // for
    } // end of high momentum pion
    else if(isMiddleMomentumGG4pi){
      for(int idx_bg = (int)(binidx/3)*3; idx_bg < (int)(binidx/3)*3+3; idx_bg++){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
          histidx = idx_bg*2*na + angleidx*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        } // for
      } // for
    } // end of middle momentum gg4pi
    else if(isHighMomentumGG4pi){
      for(int idx_bg = (int)(binidx/6)*6; idx_bg < (int)(binidx/6)*6+6; idx_bg++){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
          histidx = idx_bg*2*na + angleidx*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        } // for
      } // for
    } // end of high momentum gg4pi
    else if(isHighMomentumElectron){
      for(int idx_theta = (int)(angleidx/10)*10; idx_theta < (int)(angleidx/10)*10+10; idx_theta++){
        for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
          histidx = binidx*2*na + idx_theta*2 + idx_charge;
          h[histidx]->Fill(exraw[j]);
        } // for
      } // for
    } // end of high momentum electron
    else{
      for(int idx_charge = (int)(chargeidx/2)*2; idx_charge < (int)(chargeidx/2)*2+2; idx_charge++){
        histidx = binidx*2*na + angleidx*2 + idx_charge;
        h[histidx]->Fill(exraw[j]);
      }
    } // end of else 
      } // end of loop in momentum and angle
    } //loop j, ndedxhit
  } //loop i, entries in trk
 
  //push the histograms and betagammas to the vector
  for (int i = 0; i < nbin * na * 2; i++)
    {
      hist.Add(h[i]);
    }
  for (int i = 0; i < nbin; i++)
    {
      bg.push_back((pmin + pstep * i) / mass);
      bg.push_back((pmin + pstep * (i + 1)) / mass);
    }
}

Referenced by main().

◆ histps()

void histps ( string name )

Definition at line 325 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

                        {
  TFile f(name.c_str());
  TTree* t = (TTree*)f.Get("bin");
  Int_t total, bgb;
  t->SetBranchAddress("totalNum", &total);
  t->SetBranchAddress("betagammaBounds", &bgb);
  t->GetEntry(0);
  TH1F* h;
  stringstream ss;
  TCanvas c1("c1", "canvas", 800, 600);
  c1.Print((name + ".ps[").c_str());
  for (Int_t i = 0; i < total; i++) {
    if(i%100==0) cout << "i: " << i << " ps is written now" << endl;
    ss.str("");
    ss << "h" << i;
    h = (TH1F*)f.Get(ss.str().c_str());
    h->Draw();
    c1.Print((name + ".ps").c_str());
  }
  c1.Print((name + ".ps]").c_str());
}

Referenced by main().

◆ main()

int main ( )

Definition at line 62 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

          {
  std::cout << "in main()" << std::endl;
  //call fuction histgen to generate histograms
  TObjArray histCol;
  vector<double> bg;
  ReadCurPara();
 
  std::cout << "before histgen" << std::endl;
  histgen("../hadron_track/proton/proton.root", 0, histCol, bg);
  std::cout << "proton hists generated" << std::endl;
  histgen("../hadron_track/kaon/kaon.root", 1, histCol, bg);
  std::cout << "kaon hists generated" << std::endl;
  histgen("../hadron_track/pion/pion.root", 2, histCol, bg);
  std::cout << "pion hists generated" << std::endl;
  histgen("../hadron_track/muon/muon.root", 3, histCol, bg);
  std::cout << "muon hists generated" << std::endl;
  histgen("../hadron_track/electron/electron.root", 4, histCol, bg);
  std::cout << "electron hists generated" << std::endl;
 
  //define file and write to it
  string filename = "pion_test.root";
  TFile* f = new TFile(filename.c_str(), "recreate");
 
  int total = histCol.GetLast() + 1;
  int bg_num = bg.size();
  cout << "total: " << total << "  bg num = " << bg_num << endl;
  double betagamma[3000];
  for (unsigned int i = 0; i < bg.size(); i++) betagamma[i] = bg[i];
 
  TTree* t = new TTree("bin", "bin information");
  t->Branch("totalNum", &total, "total/I");
  t->Branch("betagammaBounds", &bg_num, "bg_num/I");
  t->Branch("betagamma", betagamma, "betagamma[bg_num]/D");
  t->Fill();
  t->Write();
    
  std::cout << "begin write hists" << std::endl;
  for (int i=0; i<total; i++){
    if(i%100==0) cout << "i: " << i << " histogram is written now" << endl;
    ((TH1F*)histCol[i])->Write();
  }
  std::cout << "write hists finish" << std::endl;
 
  f->Close();
  if(ps_flag) histps(filename);
  hist_update();
  return 0;
}

◆ ps_flag()

const bool ps_flag ( false )

Referenced by main().

Variable Documentation

◆ electron_step

const double electron_step = (double)(P_electron_max - P_electron_min) / (double)P_electron_bin

Definition at line 46 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Referenced by histgen().

◆ kaon_step

const double kaon_step = (double)(P_kaon_max - P_kaon_min) / (double)P_kaon_bin

Definition at line 49 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Referenced by histgen().

◆ muon_step

const double muon_step = (double)(P_muon_max - P_muon_min) / (double)P_muon_bin

Definition at line 47 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Referenced by histgen().

◆ na

const int na = 10

Definition at line 43 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Referenced by histgen().

◆ nc

const int nc = 2

Definition at line 44 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

◆ pion_step

const double pion_step = (double)(P_pion_max - P_pion_min) / (double)P_pion_bin

Definition at line 48 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Referenced by histgen().

◆ proton_step

const double proton_step = (double)(P_proton_max - P_proton_min) / (double)P_proton_bin

Definition at line 50 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Referenced by histgen().

◆ run_flag

const int run_flag[5] = {0,0,1,0,0}

Definition at line 59 of file Mdc/DedxCalibAlg/DedxCalibAlg-00-01-20/share/template/Simulation/histgen/test.cxx.

Referenced by histgen().

Macros

Functions

Variables

Macro Definition Documentation

◆ mass_e

◆ mass_k

◆ mass_mu

◆ mass_p

◆ mass_pi

◆ P_electron_bin

◆ P_electron_max

◆ P_electron_min

◆ P_kaon_bin

◆ P_kaon_max

◆ P_kaon_min

◆ P_muon_bin

◆ P_muon_max

◆ P_muon_min

◆ P_pion_bin

◆ P_pion_max

◆ P_pion_min

◆ P_proton_bin

◆ P_proton_max

◆ P_proton_min

Function Documentation

◆ event_validate()

◆ get_hist_range()

◆ hist_update()

◆ histgen()

◆ histps()

◆ main()

◆ ps_flag()

Variable Documentation

◆ electron_step

◆ kaon_step

◆ muon_step

◆ na

◆ nc

◆ pion_step

◆ proton_step

◆ run_flag