Functions
	if (limit!=flimit)

TFile	infile ("dedx_out.root","READ")

	switch (ifile)

c	cd ()

h	Draw ()

c2	Divide (cx, cy)

	for (int i=0;i< nsteps;i++)

graphs	Divide (2, 1)

fsum	SetParameters (parcomp)

fsum	Draw ("SAME")

cout<< "ifile "<< ifile<< endl;cout << " momentum sampled over " << nsteps<< " with step width " << 1.5/nsteps<< endl;cout << endl;cout<< "MEAN DEDX PARAMETRIZATION"<< endl;cout << "mom limits "<< mean_inf << " "<< mean_sup<< endl;cout << "mu: ";for(int param=0;param < npardedx;param++) cout << fdedx->	GetParameter (param)<< "

Variables
double	mean_inf = 0.39

double	mean_sup = 1.23

double	sigma_inf = 0.39

double	sigma_sup = 1.23

Double_t	yup = 30.

double	nsteps = 50

int	cx = 5

int	cy = 10

double	limit = 100/nsteps

double	flimit = floor(100/nsteps)

TString	fundedx = "[0] * (1./x)*2 + [1] TMath::Log(x) + [2]"

int	npardedx = 3

TString	funsig = "[0] * (1./x)*2 + [1] TMath::Log(x) + [2]"

TString	funsig_forsum = "([3] * (1./x)*2 + [4] TMath::Log(x) + [5])"

int	nparsig = 3

int	npartot = npardedx + nparsig

TCanvas *	c = new TCanvas("c", "dedx_p")

TH2F *	h

TH1D *	htemp [nsteps]

TF1 *	fgaus [nsteps]

Double_t	x [nsteps]

Double_t	mean [nsteps]

Double_t	sigma [nsteps]

Double_t	meanerr [nsteps]

Double_t	sigmaerr [nsteps]

TCanvas *	c2 = new TCanvas("c2", "slices")

TCanvas *	graphs = new TCanvas("graphs", "graphs")

Double_t	masses [5] = {0.000510999, 0.1056584, 0.139570, 0.49368, 0.9382723 }

Double_t	xerr [nsteps]

TString	funsum = fundedx + " + " + funsig_forsum

TF1 *	fsum = new TF1("fsum", funsum, 0., 1.5)

TString	fundiff = fundedx + " - " + funsig_forsum

TF1 *	fdiff = new TF1("fdiff", fundiff, 0., 1.5)

Function Documentation

c cd ( )

c2 Divide	(	cx	,
		cy
	)

graphs Divide	(	2	,
		1
	)

h Draw ( )

Definition at line 64 of file NHitsPerEvent.C.

            {
 
         std::cout << "mvdHitArray: " << mvdHitArray->GetEntries() << " projectedHitArray " << mvdProjectedArray->GetEntries() << std::endl;
 
 //      gHitPoints->Clear();
         ExtractData(mvdHitArray, gHitPoints);
         ExtractData(mvdProjectedArray, gProjections);
         //if (gHitPoints->GetN() != x.size()) gHitPoints->Expand(x.size());
 
         std::cout << "NHits: " << gHitPoints->GetN() << std::endl;
         gHitPoints->Draw("apl");
         gProjections->Draw("pl");
 }

fsum Draw ( "SAME" )

for ( )

Definition at line 71 of file dedx_bands.C.

References i, name, nsteps, and TString.

                                   {
     double stepsize = 1.5/nsteps;
     x[i] = (i + 0.5) * stepsize;
     
     TString name = "htemp"; name += i;
     TString title = "";
     TString fname = name; 
     TString fgname = name; 
     switch(ifile) {
     case 0: fname.Prepend("e_");  fgname.Prepend("fe_"); break;
     case 1: fname.Prepend("mu_"); fgname.Prepend("fmu_"); break;
     case 2: fname.Prepend("pi_"); fgname.Prepend("fpi_"); break;
     case 3: fname.Prepend("k_");  fgname.Prepend("fk_"); break;
     case 4: fname.Prepend("p_");  fgname.Prepend("fp_"); 
     }
     title = fname; title += "("; title += (x[i] - 0.5 * stepsize); title += ", "; title += (x[i] + 0.5 * stepsize); title += ") GeV"; title += ((i + 1) * limit - 1);
 
     htemp[i] = new TH1D(fname, title, 100, 0, yup);
     h->ProjectionY(fname, limit * i, (i + 1) * limit - 1,"o"); 
     c2->cd(i + 1); 
     htemp[i].Draw();
 
 
     double tmpmu = htemp[i]->GetMean() ;
     double tmpsig = htemp[i]->GetRMS();
     double gaus_inf = tmpmu - 3 * tmpsig;
     double gaus_sup = tmpmu + 3 * tmpsig;
 
     fgaus[i] = new TF1("fgname", "gaus(0)", gaus_inf, gaus_sup);
     fgaus[i]->SetParameters(1, tmpmu, tmpsig);
     htemp[i]->Fit("fgname", "R"); 
     mean[i] = fgaus[i]->GetParameter(1);
     sigma[i] = fgaus[i]->GetParameter(2);
     meanerr[i] = fgaus[i]->GetParError(1);
     sigmaerr[i] = fgaus[i]->GetParError(2);
 }

cout<< "ifile " << ifile << endl; cout << " momentum sampled over " << nsteps << " with step width " << 1.5/nsteps << endl; cout << endl; cout << "MEAN DEDX PARAMETRIZATION" << endl; cout << "mom limits " << mean_inf << " " << mean_sup << endl; cout << "mu: "; for(int param = 0; param < npardedx; param++) cout << fdedx-> GetParameter ( param )

Referenced by PndHypGeSpectrumAnalyser::DrawCalibratedSpectrum(), FitGammaSpectra(), PndHypGeSpectrumAnalyser::FitPeaks(), TtLinFitTask::MyFit(), TtCracowTask::MyFit(), and TtAliTask::MyFit().

if ( limit! = flimit )

Definition at line 20 of file dedx_bands.C.

References flimit.

                       {
     nsteps = 100/flimit;
     limit = flimit;
     cout << "changing nsteps to " << nsteps << endl;
   }

TFile infile	(	"dedx_out.root"	,
		"READ"
	)

Referenced by ana_fast(), createGeoSim(), dbuilder(), MyMainFrame::DoRecreate(), emc_correction_hist(), emc_correction_parametrization(), makeTCands(), PndFtsMapCreator2::PndFtsMapCreator2(), run_HypSimpleAna(), run_HypSimpleAna_Geo2(), simfast(), simfast_cmp(), switch(), writeMicro(), and writeMicro_Geo2().

fdiff SetParameters ( parcomp )

switch ( ifile )

Definition at line 49 of file dedx_bands.C.

References infile().

                 {
   case 0: h = (TH2F*) infile.Get("hdedx_p_e"); break;
   case 1: h = (TH2F*) infile.Get("hdedx_p_mu"); break;
   case 2: h = (TH2F*) infile.Get("hdedx_p_pi"); break;
   case 3: h = (TH2F*) infile.Get("hdedx_p_k"); break;
   case 4: h = (TH2F*) infile.Get("hdedx_p_p");
   }