Functions
h_etac_nocut	UseCurrentStyle ()

nc	Draw ()

c2	Divide (1, 2)

c2	cd (1)

c2	cd (2)

h_mphi_final	Fit (f1_phi,"R","", mean_phi-1.6 rms_phi, mean_phi+1.6 rms_phi)

h_mphi_final	Fit (f2_phi,"R","", mean1_phi-1.6 sigma1_phi, mean1_phi+1.6 sigma1_phi)

h_etac_phimass	Fit (f1,"R","", mean-1.6 rms, mean+1.6 rms)

h_etac_phimass	Fit (f2,"R","", mean1-1.6 sigma1, mean1+1.6 sigma1)

	if (vtxfit)

Variables
TString	inFile ="etac_histo_vtx.root"

TFile *	f = TFile::Open(inFile)

TVectorD *	n_etac = (TVectorD*)f->Get("n_etac")

TVectorD *	n_events = (TVectorD*)f->Get("n_events")

TH1F *	h_etac_nocut =(TH1F *)f->Get("h_etac_nocut")

TH1F *	h_etac_pid =(TH1F *)f->Get("h_etac_pid")

TH1F *	h_etac_phimass =(TH1F *)f->Get("h_etac_phimass")

TH1F *	h_etac_vtx =(TH1F *)f->Get("h_etac_vtx")

TH1F *	h_etac_4c =(TH1F *)f->Get("h_etac_4c")

TH1F *	h_mphi_nocuts =(TH1F *)f->Get("h_mphi_nocuts")

TH1F *	h_mphi_pid =(TH1F *)f->Get("h_mphi_pid")

TH1F *	h_mphi_vtx =(TH1F *)f->Get("h_mphi_vtx")

TH1F *	h_mphi_4c =(TH1F *)f->Get("h_mphi_4c")

TH1F *	h_mphi_final =(TH1F *)f->Get("h_mphi_final")

TH1F *	nc =(TH1F *)f->Get("nc")

TH1F *	h_chi2_4c =(TH1F *)f->Get("h_chi2_4c")

TH1F *	h_chi2_vtx =(TH1F *)f->Get("h_chi2_vtx")

TH1F *	hvzpos =(TH1F *)f->Get("hvzpos")

TH2F *	hvpos =(TH2F *)f->Get("hvpos")

TCanvas *	c1 =new TCanvas("c1","c1",600,600)

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

TCanvas *	c3 =new TCanvas("c3","c3",600,600)

TCanvas *	c4 =new TCanvas("c4","c4",600,600)

TCanvas *	c5 =new TCanvas("c5","c5",600,600)

double	mean_phi =1.1

double	rms_phi =0.1

TF1 *	f1_phi = new TF1("f1_phi","gaus",0.9,1.1)

double	sigma1_phi =f1_phi->GetParameter(2)

double	mean1_phi =f1_phi->GetParameter(1)

TF1 *	f2_phi = new TF1("f2_phi","gaus",0.9,1.1)

double	sigma2_phi =f2_phi->GetParameter(2)

double	mean =h_etac_phimass->GetMean()

double	rms =h_etac_phimass->GetRMS()

TF1 *	f1 = new TF1("f1","gaus",2.8,3.2)

double	sigma1 =f1->GetParameter(2)

double	mean1 =f1->GetParameter(1)

TF1 *	f2 = new TF1("f2","gaus",2.8,3.2)

double	sigma2 =f2->GetParameter(2)

Function Documentation

c2 cd ( 1 )

c2 cd ( 2 )

c2 Divide	(	1	,
		2
	)

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

h_mphi_final Fit	(	f1_phi	,
		"R"	,
		""	,
		mean_phi-1.6 *	rms_phi,
		mean_phi+1.6 *	rms_phi
	)

h_mphi_final Fit	(	f2_phi	,
		"R"	,
		""	,
		mean1_phi-1.6 *	sigma1_phi,
		mean1_phi+1.6 *	sigma1_phi
	)

h_etac_phimass Fit	(	f1	,
		"R"	,
		""	,
		mean-1.6 *	rms,
		mean+1.6 *	rms
	)

h_etac_phimass Fit	(	f2	,
		"R"	,
		""	,
		mean1-1.6 *	sigma1,
		mean1+1.6 *	sigma1
	)

if ( vtxfit )

Definition at line 115 of file plot_eta_c.C.

References c10, c7, c8, and c9.

         {
                 TCanvas *c7=new TCanvas("c7","c7",600,600);
                 hvpos->Draw();
 
                 TCanvas *c8=new TCanvas("c8","c8",600,600);
                 hvzpos->Draw();
                 
                 TCanvas *c9=new TCanvas("c9","c9",600,600);
                 h_chi2_vtx->Draw();
                 
                 TCanvas *c10=new TCanvas("c10","c10",600,600);
                 c10->Divide(1,2);
                 c10->cd(1);
                 h_mphi_vtx->Draw();
                 c10->cd(2);
                 h_etac_vtx->Draw();
         }