本文整理汇总了C++中TPaveLabel类的典型用法代码示例。如果您正苦于以下问题:C++ TPaveLabel类的具体用法?C++ TPaveLabel怎么用?C++ TPaveLabel使用的例子?那么恭喜您, 这里精选的类代码示例或许可以为您提供帮助。
在下文中一共展示了TPaveLabel类的20个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的C++代码示例。
示例1: glbox
void glbox()
{
gStyle->SetCanvasPreferGL(kTRUE);
TCanvas *c = new TCanvas("glc","TH3 Drawing", 100, 10, 850, 400);
TPaveLabel *title = new TPaveLabel(0.04, 0.86, 0.96, 0.98,
"\"glbox\" and \"glbox1\" options for TH3.");
title->SetFillColor(32);
title->Draw();
TPad *boxPad = new TPad("box", "box", 0.02, 0.02, 0.48, 0.82);
TPad *box1Pad = new TPad("box1", "box1", 0.52, 0.02, 0.98, 0.82);
boxPad->Draw();
box1Pad->Draw();
TH3F *h31 = new TH3F("h31", "h31", 10, -1, 1, 10, -1, 1, 10, -1, 1);
h31->FillRandom("gaus");
h31->SetFillColor(2);
boxPad->cd();
h31->Draw("glbox");
TH3F *h32 = new TH3F("h32", "h32", 10, -2, 2, 10, -1, 1, 10, -3, 3);
h32->FillRandom("gaus");
h32->SetFillColor(4);
box1Pad->cd();
h32->Draw("glbox1");
}
开发者ID:MycrofD,项目名称:root,代码行数:26,代码来源:glbox.C
示例2: draw_input_labels
void draw_input_labels(Int_t nInputs, Double_t* cy,
Double_t rad, Double_t layerWidth)
{
const Double_t LABEL_HEIGHT = 0.03;
const Double_t LABEL_WIDTH = 0.20;
Double_t width = LABEL_WIDTH + (layerWidth-4*rad);
Double_t margX = 0.01;
Double_t effHeight = 0.8*LABEL_HEIGHT;
TString *varNames = get_var_names(nInputs);
if (varNames == 0) exit(1);
TString input;
for (Int_t i = 0; i < nInputs; i++) {
if (i != nInputs-1) input = varNames[i];
else input = "Bias node";
Double_t x1 = margX;
Double_t x2 = margX + width;
Double_t y1 = cy[i] - effHeight;
Double_t y2 = cy[i] + effHeight;
TPaveLabel *p = new TPaveLabel(x1, y1, x2, y2, input+"", "br");
p->SetFillColor(gStyle->GetTitleFillColor());
p->SetFillStyle(1001);
p->Draw();
if (i == nInputs-1) p->SetTextColor( TMVAGlob::c_NovelBlue );
}
delete[] varNames;
}
开发者ID:TopBrussels,项目名称:TopTreeAnalysisBase,代码行数:31,代码来源:network.C
示例3: help
void help()
{
new TCanvas("chelp","Help on gldemos",200,10,700,600);
TPaveLabel *title = new TPaveLabel(0.04, 0.86, 0.96, 0.98, "These demos show different gl painters.");
title->SetFillColor(32);
title->Draw();
TPaveText *hdemo = new TPaveText(0.04, 0.04, 0.96, 0.8);
hdemo->SetTextAlign(12);
hdemo->SetTextFont(52);
hdemo->SetTextColor(kBlue);
hdemo->AddText("1. Glsurfaces demo shows glsurf4, glsurf1, glsurf3, glsurf1cyl, glsurfpol, gltf3 options.");
hdemo->AddText("2. Glrose demontrates \"glsurf2pol\" drawing option and user-defined palette.");
hdemo->AddText("3. Gltf3 demo shows \"gltf3\" option.");
hdemo->AddText("4. Glbox demo shows \"glbox\" and \"glbox1\" options for TH3.");
hdemo->AddText("5. Glparametric demo shows how to define and display parametric surfaces.");
hdemo->AddText("You can zoom any plot: press 'J', 'K', 'j', 'k' keys, or use mouse wheel.");
hdemo->AddText("Rotate any plot:");
hdemo->AddText(" ---select plot with mouse cursor,");
hdemo->AddText(" ---move mouse cursor, pressing and holding left mouse button ");
hdemo->AddText("Pan plot:");
hdemo->AddText(" ---select with mouse cursor a part of a plot, other than back box planes ,");
hdemo->AddText(" ---move mouse cursor, pressing and holding middle mouse button ");
hdemo->AddText("Selected part of a plot is higlighted (TF3 higlighting is not implemented yet.)");
hdemo->AddText("You can select one of back box planes, press middle mouse button and move cursor-");
hdemo->AddText("this will create \"slice\" (TF does not support yet).");
hdemo->AddText("After the slice was created, you can project it on a back box");
hdemo->AddText(" ---press key 'p' (now implemented only for surf options ).");
hdemo->AddText("Left double click removes all slices/projections.");
hdemo->Draw();
}
开发者ID:digideskio,项目名称:root,代码行数:34,代码来源:gldemos.C
示例4: TCanvas
TCanvas *ellipse(){
TCanvas *c1 = new TCanvas("c1");
c1->Range(0,0,1,1);
TPaveLabel *pel = new TPaveLabel(0.1,0.8,0.9,0.95,"Examples of Ellipses");
pel->SetFillColor(42);
pel->Draw();
TEllipse *el1 = new TEllipse(0.25,0.25,.1,.2);
el1->Draw();
TEllipse *el2 = new TEllipse(0.25,0.6,.2,.1);
el2->SetFillColor(6);
el2->SetFillStyle(3008);
el2->Draw();
TEllipse *el3 = new TEllipse(0.75,0.6,.2,.1,45,315);
el3->SetFillColor(2);
el3->SetFillStyle(1001);
el3->SetLineColor(4);
el3->Draw();
TEllipse *el4 = new TEllipse(0.75,0.25,.2,.15,45,315,62);
el4->SetFillColor(5);
el4->SetFillStyle(1001);
el4->SetLineColor(4);
el4->SetLineWidth(6);
el4->Draw();
return c1;
}
开发者ID:digideskio,项目名称:root,代码行数:25,代码来源:ellipse.C
示例5: TCanvas
//Draw arrows
//Author: Rene Brun
TCanvas *arrow(){
TCanvas *c1 = new TCanvas("c1");
c1->Range(0,0,1,1);
TPaveLabel *par = new TPaveLabel(0.1,0.8,0.9,0.95,"Examples of various arrow formats");
par->SetFillColor(42);
par->Draw();
TArrow *ar1 = new TArrow(0.1,0.1,0.1,0.7);
ar1->Draw();
TArrow *ar2 = new TArrow(0.2,0.1,0.2,0.7,0.05,"|>");
ar2->SetAngle(40);
ar2->SetLineWidth(2);
ar2->Draw();
TArrow *ar3 = new TArrow(0.3,0.1,0.3,0.7,0.05,"<|>");
ar3->SetAngle(40);
ar3->SetLineWidth(2);
ar3->Draw();
TArrow *ar4 = new TArrow(0.46,0.7,0.82,0.42,0.07,"|>");
ar4->SetAngle(60);
ar4->SetLineWidth(2);
ar4->SetFillColor(2);
ar4->Draw();
TArrow *ar5 = new TArrow(0.4,0.25,0.95,0.25,0.15,"<|>");
ar5->SetAngle(60);
ar5->SetLineWidth(4);
ar5->SetLineColor(4);
ar5->SetFillStyle(3008);
ar5->SetFillColor(2);
ar5->Draw();
return c1;
}
开发者ID:ancapopescu13,项目名称:root,代码行数:34,代码来源:drawArrow.C
示例6: testRoot
void testRoot() {
//Fill a 1-D histogram from a parametric function
// To see the output of this macro, click begin_html <a href="gif/fillrandom.gif">here</a>. end_html
//Author: Rene Brun
TCanvas *c1 = new TCanvas("c1","The FillRandom example",200,10,700,900);
c1->SetFillColor(18);
TPad* pad1 = new TPad("pad1","The pad with the function",0.05,0.50,0.95,0.95,21);
TPad* pad2 = new TPad("pad2","The pad with the histogram",0.05,0.05,0.95,0.45,21);
pad1->Draw();
pad2->Draw();
pad1->cd();
gBenchmark->Start("fillrandom");
//
// A function (any dimension) or a formula may reference
// an already defined formula
//
TFormula* form1 = new TFormula("form1","abs(sin(x)/x)");
TF1* sqroot = new TF1("sqroot","x*gaus(0) + [3]*form1",0,10);
sqroot->SetParameters(10,4,1,20);
pad1->SetGridx();
pad1->SetGridy();
pad1->GetFrame()->SetFillColor(42);
pad1->GetFrame()->SetBorderMode(-1);
pad1->GetFrame()->SetBorderSize(5);
sqroot->SetLineColor(4);
sqroot->SetLineWidth(6);
sqroot->Draw();
TPaveLabel* lfunction = new TPaveLabel(5,39,9.8,46,"The sqroot function");
lfunction->SetFillColor(41);
lfunction->Draw();
c1->Update();
//
// Create a one dimensional histogram (one float per bin)
// and fill it following the distribution in function sqroot.
//
pad2->cd();
pad2->GetFrame()->SetFillColor(42);
pad2->GetFrame()->SetBorderMode(-1);
pad2->GetFrame()->SetBorderSize(5);
TH1F* h1f = new TH1F("h1f","Test random numbers",200,0,10);
h1f->SetFillColor(45);
h1f->FillRandom("sqroot",10000);
h1f->Draw();
c1->Update();
//
// Open a ROOT file and save the formula, function and histogram
//
TFile myfile("fillrandom.root","RECREATE");
form1->Write();
sqroot->Write();
h1f->Write();
gBenchmark->Show("fillrandom");
}
开发者ID:hqvigstad,项目名称:alilink,代码行数:58,代码来源:testRoot.C
示例7: showjets
//----------------------------------------------------------------------
/// show the jets contained in filename (as produced by
/// ClusterSequence::print_jets_for_root()), with an optional label
TCanvas * showjets (const char * filename, const char * label = 0) {
// display the various 2-d drawing options
gROOT->Reset();
// set up canvas
TCanvas * lego = new TCanvas("lego","lego options",400,50,800,600);
lego->SetTheta(30.0);
lego->SetPhi(20.0);
// orientation used for plots in subtraction paper
//lego->SetTheta(62.15);
//lego->SetPhi(9.15);
////vector<double> col
set_default_colours(lego);
TPaveLabel pl;
JetHist * jets = new JetHist(filename);
jets->stack.Draw("lego1"); // cyl does not work with 5.16
if (label != 0) {
Float_t x1=0.63, y1=0.875, x2=0.95, y2=0.925;
pl.DrawPaveLabel(x1,y1,x2,y2,label,"brNDC");
} else if (jets->comment() != "") {
Float_t x1=0.15, y1=0.875, x2=0.95, y2=0.925;
pl.DrawPaveLabel(x1,y1,x2,y2,jets->comment().c_str(),"brNDC");
}
// normal histogram labels not working, so draw them by hand
TLatex l;
l.SetTextAlign(22);
l.SetTextSize(0.05);
//l.DrawLatex(0.0,0.85,"anti-k_{t}, R=1");
l.SetTextSize(0.04);
l.DrawLatex(0.20,-0.98,"y");
l.SetTextAlign(32);
l.DrawLatex(-0.7,0.8,"p_{t} [GeV]");
l.DrawLatex(-0.6,-0.78,"#phi");
// do not delete jets -- otherwise you lose everything!;
return lego;
///
lego->Update();
}
开发者ID:Dr15Jones,项目名称:fastjet,代码行数:51,代码来源:jet-plots.C
示例8: DoFit
void DoFit(const char* fitter, TVirtualPad *pad, Int_t npass) {
TStopwatch timer;
TVirtualFitter::SetDefaultFitter(fitter);
pad->SetGrid();
fitFcn->SetParameters(100,0,0,2,7);
fitFcn->Update();
timer.Start();
histo->Fit("fitFcn","0");
timer.Stop();
histo->Draw();
Double_t cputime = timer.CpuTime();
printf("%s, npass=%d : RT=%7.3f s, Cpu=%7.3f s\n",fitter,npass,timer.RealTime(),cputime);
TPaveLabel *p = new TPaveLabel(0.5,0.7,0.85,0.8,Form("%s CPU= %g s",fitter,cputime),"brNDC");
p->Draw();
pad->Update();
}
开发者ID:digideskio,项目名称:root,代码行数:18,代码来源:minuit2FitBench2D.C
示例9: first
void first() {
TCanvas *nut = new TCanvas("nut", "FirstSession",100,10,700,900);
nut->Range(0,0,20,24);
nut->SetFillColor(10);
nut->SetBorderSize(2);
TPaveLabel *pl = new TPaveLabel(3,22,17,23.7,
"My first ROOT interactive session","br");
pl->SetFillColor(18);
pl->Draw();
TText t(0,0,"a");
t.SetTextFont(62);
t.SetTextSize(0.025);
t.SetTextAlign(12);
t.DrawText(2,20.3,"ROOT is based on CINT, a powerful C/C++ interpreter.");
t.DrawText(2,19.3,"Blocks of lines can be entered within {...}.");
t.DrawText(2,18.3,"Previous typed lines can be recalled.");
t.SetTextFont(72);
t.SetTextSize(0.026);
t.DrawText(3,17,"Root > float x=5; float y=7;");
t.DrawText(3,16,"Root > x*sqrt(y)");
t.DrawText(3,14,
"Root > for (int i=2;i<7;i++) printf(\"sqrt(%d) = %f\\n\",i,sqrt(i));");
t.DrawText(3,10,"Root > TF1 f1(\"f1\",\"sin(x)/x\",0,10)");
t.DrawText(3, 9,"Root > f1.Draw()");
t.SetTextFont(81);
t.SetTextSize(0.018);
t.DrawText(4,15,"(float) 13.2288f");
t.DrawText(4,13.3,"sqrt(2) = 1.414214");
t.DrawText(4,12.7,"sqrt(3) = 1.732051");
t.DrawText(4,12.1,"sqrt(4) = 2.000000");
t.DrawText(4,11.5,"sqrt(5) = 2.236068");
t.DrawText(4,10.9,"sqrt(6) = 2.449490");
TPad *pad = new TPad("pad","pad",.2,.05,.8,.35);
pad->Draw();
pad->cd();
pad->SetGrid();
TF1 *f1 = new TF1("f1","sin(x)/x",0,10);
f1->Draw();
}
开发者ID:digideskio,项目名称:root,代码行数:44,代码来源:first.C
示例10: glrose
// Render a TF2 looking like a rose.
// Author: Timur Pocheptsov
void glrose()
{
//Define and set user's palette,
//use polar system.
const Int_t paletteSize = 10;
Float_t rgb[paletteSize * 3] =
{0.80f, 0.55f, 0.40f,
0.85f, 0.60f, 0.45f,
0.90f, 0.65f, 0.50f,
0.95f, 0.70f, 0.55f,
1.f, 0.75f, 0.60f,
1.f, 0.80f, 0.65f,
1.f, 0.85f, 0.70f,
1.f, 0.90f, 0.75f,
1.f, 0.95f, 0.80f,
1.f, 1.f, 0.85f};
Int_t palette[paletteSize] = {0};
for (Int_t i = 0; i < paletteSize; ++i)
palette[i] = TColor::GetColor(rgb[i * 3], rgb[i * 3 + 1], rgb[i * 3 + 2]);
gStyle->SetPalette(paletteSize, palette);
gStyle->SetCanvasPreferGL(1);
TCanvas *cnv = new TCanvas("glc", "Surface sample", 200, 10, 600, 550);
TPaveLabel *title = new TPaveLabel(0.04, 0.86, 0.96, 0.98,
"\"glsurf2pol\" option + user defined palette.");
title->SetFillColor(32);
title->Draw();
TPad *rosePad = new TPad("box", "box", 0.04, 0.04, 0.96, 0.8);
rosePad->Draw();
TF2 *fun = new TF2("a", "cos(y)*sin(x)+cos(x)*sin(y)", -6, 6, -6, 6);
fun->SetContour(paletteSize);
fun->SetNpx(30);
fun->SetNpy(30);
rosePad->cd();
fun->Draw("glsurf2pol");
}
开发者ID:digideskio,项目名称:root,代码行数:45,代码来源:glrose.C
示例11: draw_layer_labels
void draw_layer_labels(Int_t nLayers)
{
const Double_t LABEL_HEIGHT = 0.03;
const Double_t LABEL_WIDTH = 0.20;
Double_t effWidth = 0.8*(1.0-LABEL_WIDTH)/nLayers;
Double_t height = 0.8*LABEL_HEIGHT;
Double_t margY = LABEL_HEIGHT - height;
for (Int_t i = 0; i < nLayers; i++) {
TString label = Form("Layer %i", i);
Double_t cx = i*(1.0-LABEL_WIDTH)/nLayers+1.0/(2.0*nLayers)+LABEL_WIDTH;
Double_t x1 = cx-0.8*effWidth/2.0;
Double_t x2 = cx+0.8*effWidth/2.0;
Double_t y1 = margY;
Double_t y2 = margY + height;
TPaveLabel *p = new TPaveLabel(x1, y1, x2, y2, label+"", "br");
p->SetFillColor(gStyle->GetTitleFillColor());
p->SetFillStyle(1001);
p->Draw();
}
}
开发者ID:TopBrussels,项目名称:TopTreeAnalysisBase,代码行数:22,代码来源:network.C
示例12: KolmogorovTest
////////////////////////////////////////////////////////////
//
// This function performs a compatibility test between two
// histogram based on the Kolmogorov-Smirnof algorithm. It
// also prints the value in a TPaveLabel at the upper-right
// corner.
// The return value contains the result of the test
//
double KolmogorovTest(TH1 *h1, TH1 *h2){
double mya_array[1300], myb_array[1300];
vector<double> mya;
vector<double> myb;
for (int i=0; i<h1->GetNbinsX(); i++){
mya.push_back(h1->GetBinContent(i+1));
myb.push_back(h2->GetBinContent(i+1));
}
sort(mya.begin(),mya.end());
sort(myb.begin(),myb.end());
copy(mya.begin(),mya.end(),mya_array);
copy(myb.begin(),myb.end(),myb_array);
const int nbinsa = h1->GetNbinsX();
const int nbinsb = h2->GetNbinsX();
double kstest = TMath::KolmogorovTest(nbinsa, mya_array,
nbinsb, myb_array,
"UOX");
if (DEBUGP) cout << " + KS value = " << kstest << endl;
// Create text with the value
TString legend = Form("KS=%4.2f", kstest);
// Create a pave text to put the value inside
TPaveLabel* pl = new TPaveLabel(0.79,0.91,0.93,0.96, legend.Data(), "NDC");
// Tune style
//pl->SetTextSize(0.04);
pl->SetLineColor(41);
pl->SetLineWidth(1);
pl->SetLineStyle(1);
pl->SetFillColor(41);
pl->SetBorderSize(3);
if (kstest < 0.7)
pl->SetTextColor(kRed);
pl->Draw();
return kstest;
}
开发者ID:cardinia,项目名称:cmssw,代码行数:55,代码来源:new_PlotHelpers.C
示例13: qa2
void qa2() {
//Fill a 1-D histogram from a parametric function
TCanvas *c1 = new TCanvas("c1","The FillRandom example",0,0,700,500);
c1->SetFillColor(18);
gBenchmark->Start("fillrandom");
//
// A function (any dimension) or a formula may reference
// an already defined formula
//
TFormula *form1 = new TFormula("form1","abs(sin(x)/x)");
TF1 *sqroot = new TF1("sqroot","x*gaus(0) + [3]*form1",0,10);
sqroot->SetParameters(10,4,1,20);
//
// Create a one dimensional histogram (one float per bin)
// and fill it following the distribution in function sqroot.
//
TH1F *h1f = new TH1F("h1f","Test random numbers",200,0,10);
h1f->SetFillColor(45);
h1f->FillRandom("sqroot",100000);
h1f->Draw();
TPaveLabel *lfunction = new TPaveLabel(5,39,9.8,46,"The sqroot function");
lfunction->SetFillColor(41);
c1->SetGridx();
c1->SetGridy();
c1->GetFrame()->SetFillColor(42);
c1->GetFrame()->SetBorderMode(-1);
c1->GetFrame()->SetBorderSize(5);
h1f->SetDirectory(0);
c1->Update();
sqroot->SetParameters(200,4,1,20);
}
开发者ID:digideskio,项目名称:root,代码行数:37,代码来源:qa2.C
示例14: memstat
void memstat(double update=0.01, const char* fname="*") {
// Open the memstat data file, then call TTree::Draw to precompute
// the arrays of positions and nbytes per entry.
// update is the time interval in the data file in seconds after which
// the display is updated. For example is the job producing the memstat.root file
// took 100s to execute, an update of 0.1s will generate 1000 time views of
// the memory use.
// if fname=="*" (default), the most recent file memstat*.root will be taken.
TString s;
if (!fname || strlen(fname) <5 || strstr(fname,"*")) {
//take the most recent file memstat*.root
s = gSystem->GetFromPipe("ls -lrt memstat*.root");
Int_t ns = s.Length();
fname = strstr(s.Data()+ns-25,"memstat");
}
printf("Analyzing file: %s\n",fname);
f = TFile::Open(fname);
if (!f) {
printf("Cannot open file %s\n",fname);
return;
}
T = (TTree*)f->Get("T");
if (!T) {
printf("cannot find the TMemStat TTree named T in file %s\n",fname);
return;
}
if (update <= 0) {
printf("Illegal update value %g, changed to 0.01\n",update);
update = 0.01;
}
if (update < 0.001) printf("Warning update parameter is very small, processing may be slow\n");
Long64_t nentries = T->GetEntries();
T->SetEstimate(nentries+10);
Long64_t nsel = T->Draw("pos:nbytes:time:btid","","goff");
//now we compute the best binning for the histogram
Int_t nbytes;
Double_t pos;
V1 = T->GetV1();
V2 = T->GetV2();
V3 = T->GetV3();
V4 = T->GetV4();
Long64_t imean = (Long64_t)TMath::Mean(nsel,V1);
Long64_t irms = (Long64_t)TMath::RMS(nsel,V1);
//Long64_t bw = 10000;
Long64_t bw = 1000;
imean = imean - imean%bw;
irms = irms -irms%bw;
Int_t nbins = Int_t(4*irms/bw);
Long64_t ivmin = imean -bw*nbins/2;
Long64_t ivmax = ivmin+bw*nbins;
if (ivmax > 2000000000 && ivmin <2000000000) {
//the data set has been likely generated on a 32 bits machine
//we are mostly interested by the small allocations, so we select
//only values below 2 GBytes
printf("memory locations above 2GBytes will be ignored\n");
nsel = T->Draw("pos:nbytes:time:btid","pos <2e9","goff");
V1 = T->GetV1();
V2 = T->GetV2();
V3 = T->GetV3();
V4 = T->GetV4();
imean = (Long64_t)TMath::Mean(nsel,V1);
irms = (Long64_t)TMath::RMS(nsel,V1);
bw = 10000;
imean = imean - imean%bw;
irms = irms -irms%bw;
nbins = Int_t(4*irms/bw);
ivmin = imean -bw*nbins/2;
ivmax = ivmin+bw*nbins;
}
update *= 0.0001*V3[nsel-1]; //convert time per cent in seconds
Long64_t nvm = Long64_t(ivmax-ivmin+1);
Long64_t *nbold = new Long64_t[nvm];
Int_t *ientry = new Int_t[nvm];
memset(nbold,0,nvm*8);
Double_t dv = (ivmax-ivmin)/nbins;
h = new TH1D("h",Form("%s;pos;per cent of pages used",fname),nbins,ivmin,ivmax);
TAxis *axis = h->GetXaxis();
gStyle->SetOptStat("ie");
h->SetFillColor(kRed);
h->SetMinimum(0);
h->SetMaximum(100);
halloc = new TH1D("halloc",Form("%s;pos;number of mallocs",fname),nbins,ivmin,ivmax);
hfree = new TH1D("hfree", Form("%s;pos;number of frees",fname),nbins,ivmin,ivmax);
//open a canvas and draw the empty histogram
TCanvas *c1 = new TCanvas("c1","c1",1200,600);
c1->SetFrameFillColor(kYellow-3);
c1->SetGridx();
c1->SetGridy();
h->Draw();
//create a TPaveText to show the summary results
TPaveText *pvt = new TPaveText(.5,.9,.75,.99,"brNDC");
pvt->Draw();
//create a TPaveLabel to show the time
TPaveLabel *ptime = new TPaveLabel(.905,.7,.995,.76,"time","brNDC");
ptime->SetFillColor(kYellow-3);
ptime->Draw();
//.........这里部分代码省略.........
开发者ID:My-Source,项目名称:root,代码行数:101,代码来源:memstat.C
示例15: advancedNoiseAnalysis
//.........这里部分代码省略.........
tempTitle = "NoiseDist Det. " + toString( iDetector ) + " - Ch. " + toString( iChan ) ;
TH1D * noiseDistCh = new TH1D( tempName.c_str(), tempTitle.c_str(), 50, 0., 10. );
noiseDistCh->SetXTitle("Noise [ADC]");
noiseDistCh->SetLineColor( kColor[iDetector] );
noiseDistCh->SetLineStyle( iChan + 2 );
noiseDistCh->SetLineWidth( 2 );
outputHistoList->Add( noiseDistCh );
// let's start looping on pixels now
for ( size_t yPixel = 1 ; yPixel <= kYPixel ; ++yPixel ) {
for ( size_t xPixel = xLimit[ iChan ] + 1; xPixel <= xLimit[ iChan +1 ] ; ++xPixel ) {
double noise = noiseMap->GetBinContent( xPixel , yPixel );
noiseMapCh->Fill( xPixel - 1 , yPixel - 1, noise );
noiseDistCh->Fill( noise );
}
}
canvas->cd( iChan + 1 ) ;
noiseMapCh->Draw("colz");
canvas->cd( iChan + kNChan + 1 );
noiseDistCh->Draw();
topPad->cd( iDetector + 1 );
if ( iChan == 0 ) {
noiseDistCh->Draw();
} else {
noiseDistCh->Draw("same");
}
middlePad->cd( iChan + 1 );
if ( iDetector == 0 ) {
noiseDistCh->Draw();
} else {
noiseDistCh->Draw("same");
}
noiseMean[ kNChan * iDetector + iChan ] = noiseDistCh->GetMean();
noiseRMS[ kNChan * iDetector + iChan ] = noiseDistCh->GetRMS();
}
canvas->Write();
}
canvasName = "summary";
canvasTitle = "Noise summary";
TCanvas * summaryCanvas = new TCanvas( canvasName.c_str(), canvasTitle.c_str(), 1000, 500 );
summaryCanvas->SetGridx(1);
TLegend * legend = new TLegend(0.5, 4.8, 1.5, 4.3,"","br");;
for ( size_t iDetector = 0 ; iDetector < nDetector ; ++iDetector ) {
TGraphErrors * gr = new TGraphErrors( kNChan, channel, &noiseMean[ iDetector * kNChan ], NULL, &noiseRMS[ iDetector * kNChan ] );
gr->SetName( string( "NoisePerChannel_d" + toString( iDetector )).c_str());
gr->SetTitle(string("Detector " + toString( iDetector )).c_str());
gr->GetXaxis()->SetTitle("Channel #");
gr->GetYaxis()->SetTitle("Noise [ADC]");
gr->GetXaxis()->SetNdivisions( 5 );
gr->GetXaxis()->SetLabelSize( 0 );
gr->SetMarkerStyle( iDetector + 1 );
gr->SetMarkerColor( kColor[iDetector] );
gr->SetLineColor( kColor[iDetector] );
gr->SetLineWidth( 2 );
legend->AddEntry( gr, string("Detector " + toString( iDetector )).c_str(), "LP");
if ( iDetector == 0 ) {
gr->Draw("ALP");
} else {
gr->Draw("LP");
}
}
legend->Draw();
for ( size_t iChan = 0 ; iChan < kNChan ; ++iChan ) {
TPaveLabel * label = new TPaveLabel( iChan - 0.75 , 3.2 , iChan -0.25 , 3, string("Ch " + toString( iChan ) ).c_str());
label->Draw();
}
summaryCanvas->Write();
comparisonCanvas->Write();
outputHistoList->Write();
}
开发者ID:AlexanderMorton,项目名称:eutelescope,代码行数:101,代码来源:advancedNoiseAnalysis.C
示例16: file
void file(){
TCanvas *c1 = new TCanvas("c1","ROOT File description",200,10,700,550);
c1->Range(0,-0.25,21,14);
TPaveLabel *title = new TPaveLabel(5,12,15,13.7,c1->GetTitle());
title->SetFillColor(16);
title->Draw();
// horizonthal file layout
TPave *file = new TPave(1,8.5,20,11);
file->SetFillColor(11);
file->Draw();
TPave *fileh = new TPave(1,8.5,2.5,11);
fileh->SetFillColor(44);
fileh->Draw();
TPave *lrh = new TPave(2.5,8.5,3.3,11,1);
lrh->SetFillColor(33);
lrh->Draw();
lrh->DrawPave(6.9,8.5,7.7,11,1);
lrh->DrawPave(10.5,8.5,11.3,11,1);
lrh->DrawPave(14.5,8.5,15.3,11,1);
TLine *ldot = new TLine(1,8.5,0.5,6.5);
ldot->SetLineStyle(2);
ldot->Draw();
ldot->DrawLine(2.5, 8.5, 9.4, 6.5);
ldot->DrawLine(10.5, 8.5, 10, 6.5);
ldot->DrawLine(11.3, 8.5, 19.5, 6.5);
TLine *line = new TLine(2.6,11,2.6,11.5);
line->Draw();
line->DrawLine(2.6,11.5,7,11.5);
TArrow *arrow = new TArrow(7,11.5,7,11.1,0.01,"|>");
arrow->SetFillStyle(1001);
arrow->Draw();
line->DrawLine( 7, 8.5, 7, 8.0);
line->DrawLine( 7, 8.0, 10.6, 8);
arrow->DrawArrow( 10.6,8, 10.6, 8.4,0.01,"|>");
line->DrawLine( 10.6, 11, 10.6, 11.5);
line->DrawLine( 10.6, 11.5, 14.6, 11.5);
arrow->DrawArrow( 14.6,11.5, 14.6,11.1,0.01,"|>");
line->DrawLine( 14.6, 8.5, 14.6, 8.0);
line->DrawLine( 14.6, 8.0, 16, 8);
ldot->DrawLine(16, 8, 19, 8);
TText *vert = new TText(1.5,9.75,"File");
vert->SetTextAlign(21);
vert->SetTextAngle(90);
vert->SetTextSize(0.025);
vert->Draw();
vert->DrawText(2.0, 9.75,"Header");
vert->DrawText(2.9, 9.75,"Logical Record");
vert->DrawText(3.2, 9.75,"Header");
vert->DrawText(7.3, 9.75,"Logical Record");
vert->DrawText(7.6, 9.75,"Header");
vert->DrawText(10.9,9.75,"Logical Record");
vert->DrawText(11.2,9.75,"Header");
vert->DrawText(14.9,9.75,"Logical Record");
vert->DrawText(15.2,9.75,"Header");
TText *hori = new TText(4.75,10,"Object");
hori->SetTextAlign(22);
hori->SetTextSize(0.035);
hori->Draw();
hori->DrawText(4.75, 9.5,"Data");
hori->DrawText(9.2, 10,"Deleted");
hori->DrawText(9.2, 9.5,"Object");
line->DrawLine( 6.9, 8.5, 10.5, 11);
line->DrawLine( 6.9, 11, 10.5, 8.5);
TText *tbig = new TText(17,9.75,"............");
tbig->SetTextAlign(22);
tbig->SetTextSize(0.03);
tbig->Draw();
tbig->DrawText(2.6, 7, "fBEGIN");
tbig->DrawText(20., 7, "fEND");
arrow->DrawArrow( 2.6,7, 2.6,8.4,0.01,"|>");
arrow->DrawArrow( 20,7, 20,8.4,0.01,"|>");
//file header
TPaveText *header = new TPaveText(0.5,.2,9.4,6.5);
header->SetFillColor(44);
header->Draw();
TText *fh=header->AddText("File Header");
fh->SetTextAlign(22);
fh->SetTextSize(0.04);
header->SetTextSize(0.027);
header->SetTextAlign(12);
header->AddText(" ");
header->AddLine(0,0,0,0);
header->AddText("\"root\": Root File Identifier");
header->AddText("fVersion: File version identifier");
header->AddText("fBEGIN: Pointer to first data record");
header->AddText("fEND: Pointer to first free word at EOF");
header->AddText("fSeekFree: Pointer to FREE data record");
header->AddText("fNbytesFree: Number of bytes in FREE");
header->AddText("fNfree: Number of free data records");
header->AddText("fNbytesName: Number of bytes in name/title");
header->AddText("fUnits: Number of bytes for pointers");
header->AddText("fCompress: Compression level");
//logical record header
TPaveText *lrecord = new TPaveText(10,0.2,19.5,6.5);
lrecord->SetFillColor(33);
//.........这里部分代码省略.........
开发者ID:MycrofD,项目名称:root,代码行数:101,代码来源:file.C
示例17: ST_Monitoring_Eff
// File: ST_Monitoring_Efficiency.C
// Last Modified: 05/27/2015
// Creator: Mahmoud Kamel [email protected]
// Purpose: Displaying histograms for online monitoring purposes
void ST_Monitoring_Eff () {
// Define the directory that contains the histograms
TDirectory *dir = (TDirectory*)gDirectory->FindObjectAny("st_tracking");
if(dir) dir->cd();
// Grab 1D histograms
TH1D *MacropEff = (TH1D*)dir->FindObjectAny("MacropEff");
TH1D *MacropEff_adc = (TH1D*)dir->FindObjectAny("MacropEff_adc");
TH1D *h_phi_sec_pred_hit_cntr = (TH1D*)dir->FindObjectAny("h_phi_sec_pred_hit_cntr");
TH1D *h_phi_sec_hit_cntr = (TH1D*)dir->FindObjectAny("h_phi_sec_hit_cntr");
TH1D *h_phi_sec_adc_cntr = (TH1D*)dir->FindObjectAny("h_phi_sec_adc_cntr");
// get Binomial errors in an efficiency plot
// hit object efficiency
h_phi_sec_hit_cntr->Sumw2();
h_phi_sec_pred_hit_cntr->Sumw2();
MacropEff->Sumw2();
MacropEff->Divide(h_phi_sec_hit_cntr,h_phi_sec_pred_hit_cntr,1,1,"B");
// adc efficiency
h_phi_sec_adc_cntr->Sumw2();
MacropEff_adc->Sumw2();
MacropEff_adc->Divide(h_phi_sec_adc_cntr,h_phi_sec_pred_hit_cntr,1,1,"B");
//Create the canvas
if(gPad == NULL)
{
TCanvas *c1 = new TCanvas("c1","Start Counter 1D Histograms", 200, 10, 600, 480);
c1->cd(0);
c1->Draw();
c1->Update();
}
if(!gPad) return;
TCanvas *c1 = gPad->GetCanvas();
c1->Divide(2,1);
// ST ADC Efficiency histogram
c1->cd(1);
gStyle->SetOptStat(0);
gStyle->SetErrorX(0);
gPad->SetTicks();
gPad->SetGrid();
if (MacropEff_adc) {
MacropEff_adc->Draw("E1");
MacropEff_adc->SetMarkerStyle(21);
MacropEff_adc->SetMarkerSize(1.5);
MacropEff_adc->SetMarkerColor(4.0);
MacropEff_adc->SetAxisRange(0., 1.,"Y");
MacropEff_adc->GetYaxis()->SetTitleOffset(1.26);
Double_t TWA_adc=0;
Double_t sumE_adc=0;
Double_t Final_adc=0;
Double_t error_adc=0.;
Double_t BC_adc=0.;
Double_t WA_adc=0.;
//Double_t Final_adc=0.;
for (int i = 0; i < 30; i++)
{
error_adc = MacropEff_adc->GetBinError(i+2);
sumE_adc = sumE_adc + error_adc;
BC_adc = MacropEff_adc->GetBinContent(i+2);
WA_adc = BC_adc*error_adc;
TWA_adc = TWA_adc + WA_adc ;
Final_adc = TWA_adc/sumE_adc;
// cout << "error_adc = "<< error_adc << endl;
// cout << "BC_adc = "<< BC_adc << endl;
// cout << "Final_adc = "<< Final_adc << endl;
}
TLine *line = new TLine(1,Final_adc,30,Final_adc);
line->SetLineWidth(3);
line->SetLineColor(2);
//Write the eff value on the histogram
char tFinal_adc[40];
char terror_adc[40];
sprintf(tFinal_adc,"ADC Efficiency (%%) = %g",Final_adc*100);
sprintf(terror_adc,"ADC Efficiency error (%%) = %g",error_adc*100);
line->Draw();
TPaveLabel *p = new TPaveLabel(0.3,0.6,0.7,0.7,tFinal_adc,"brNDC");
p->Draw();
TPaveLabel *p1 = new TPaveLabel(0.3,0.5,0.7,0.6,terror_adc,"brNDC");
p1->Draw();
}
// ST Hit Efficiency histogram
c1->cd(2);
gPad->SetTicks();
gPad->SetGrid();
if (MacropEff) {
MacropEff->Draw("E1");
MacropEff->SetMarkerStyle(21);
MacropEff->SetMarkerSize(1.5);
MacropEff->SetMarkerColor(4.0);
MacropEff->SetAxisRange(0., 1.,"Y");
MacropEff->GetYaxis()->SetTitleOffset(1.26);
Double_t TWA=0;
Double_t sumE=0;
Double_t Final=0;
Double_t error=0.;
Double_t BC=0.;
//.........这里部分代码省略.........
开发者ID:JeffersonLab,项目名称:sim-recon,代码行数:101,代码来源:ST_Monitoring_Eff.C
示例18: file
void
show
(
char* inpDir = "", // MuDST directory
char* inpFile = "show.lis", // MuDST file(s);
char* outFile = "show.root",// output tree file
Int_t nFiles = 50, // # of MuDST file(s)
Int_t nEvents = 100 // # of events
)
// remeber to adjust dbase timestamp below !!!!
// what a ... design
{
//gErrorIgnoreLevel=1999;
// load root/root4star libraries
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
// load more libraries :)
gSystem->Load("libmysqlclient");
gSystem->Load("libGeom");
gSystem->Load("StDbLib");
gSystem->Load("StDbBroker");
gSystem->Load("St_db_Maker");
// load even more libraries (EEMC stuff)
gSystem->Load("StEEmcUtil");
gSystem->Load("StEEmcDbMaker");
gSystem->Load("StEEmcPoolTTM");
// create the chain
chain = new StChain("StChain");
//
now = new TDatime;
// for display
TCanvas *c1 = new TCanvas("eemc","eemc",10,10,1000,1000);
TPaveLabel *tlab = new TPaveLabel(-0.99,+0.99,+0.99,+0.90,
"EEMC TOWERS & TPC TRACKS Piotr A Zolnierczuk (IU)");
eventInfo = new TPaveText (-0.99,-0.99,+0.0 ,-0.75);
dateInfo = new TPaveLabel(+0.60,-0.99,+0.99,-0.95,now->AsString());
TGeoManager *gm = new TGeoManager("eemc", "eemc tower display");
TGeoVolume *top = gm->MakeBox("star",0, 200., 200., 350.);
TGeoVolume *smbox = gm->MakeBox("smbox1",0, 2., 2., 2.);
smbox->SetLineColor(kRed);
// eemc
eemc = new EEmcTTDisplay();
eemc->SetMagneticField(0.5); // in Tesla
eemc->SetShowExtrapolatedTracks(true);
TGeoTranslation *etra = new TGeoTranslation(0.0,0.0,0.5*(eemc->getZ1()+eemc->getZ2()));
top->AddNode(smbox, 1,NULL);
top->AddNode(eemc(),1,etra);
gm->SetTopVolume(top);
gm->CloseGeometry();
gm->SetVisLevel(4);
gm->SetVisOption(0);
c1->SetTheta(90);
c1->SetPhi(0);
top->Draw();
tlab->Draw();
gPad->Update();
// now we add Makers to the chain... some of that is black magic :)
muDstMk = new StMuDstMaker(0,0,inpDir,inpFile,"",nFiles); // muDST main chain
StMuDbReader *db = StMuDbReader::instance(); // need the database
St_db_Maker *dbMk = new St_db_Maker("StarDb", "MySQL:StarDb"); // need another db(?)
new StEEmcDbMaker("eemcDb"); // need EEMC database
// now comment in/out/change the below if you want it your way
dbMk->setTimeStampDay(20040331); // format: yyyymmdd
// finally after so many lines we arrive at the good stuff
ttm = new EEmcTTMMaker ("TTM",muDstMk,eemcDbMk);
ttm->Summary(cout); //
StMuDebug::setLevel(0);
chain->Init();
StEEmcDb *eemcDb = (StEEmcDb*)chain->GetDataSet("StEEmcDb");
eemcDb->setSectors(1,12); // request EEMC DB for sectors you need (dafault:1-12)
eemcDb->setPreferedFlavor("onlped","eemcPMTped"); // request alternative flavor
chain->ls(3);
//---------------------------------------------------
next();
}
开发者ID:star-bnl,项目名称:star-emc,代码行数:95,代码来源:show.C
示例19: zlljets_ZpT_recoGenAna
//.........这里部分代码省略.........
Double_t lastEdge;
TH1D *hZpTMuEleRecoRatio;
TCanvas *cZpTMuEleRecoRatio;
if( hZpTreco[0]->GetNbinsX() == hZpTreco[1]->GetNbinsX() ) {
bins = hZpTreco[0]->GetNbinsX();
firstEdge = hZpTreco[0]->GetXaxis()->GetBinLowEdge(1);
lastEdge = hZpTreco[0]->GetXaxis()->GetBinUpEdge(bins);
hZpTMuEleRecoRatio = new TH1D("hZpTMuEleRecoRatio","",bins,firstEdge,lastEdge);
// cout << "hZpTreco[0]->GetNbinsX() : "<< hZpTreco[0]->GetNbinsX() << endl;
// cout << "hZpTreco[1]->GetNbinsX() : "<< hZpTreco[1]->GetNbinsX() << endl;
// cout << "hZpTMuEleGenRatio->GetNbinsX() : "<< hZpTMuEleGenRatio->GetNbinsX() << endl;
if (!hZpTMuEleRecoRatio->Divide(hZpTreco[0],hZpTreco[1])) cout << " Error in hZpTMuEleRecoRatio->Divide(hZpTreco[0],hZpTreco[1])" << endl;
hZpTMuEleRecoRatio->SetStats(kFALSE);
cZpTMuEleRecoRatio = new TCanvas("cZpTMuEleRecoRatio","");
hZpTMuEleRecoRatio->GetXaxis()->SetTitle("Z_{pT}[GeV]");
hZpTMuEleRecoRatio->GetXaxis()->SetTitleSize(0.04);
hZpTMuEleRecoRatio->GetYaxis()->SetTitle("reco pT_{Z(#mu#mu)} / pT_{Z(ee)} ");
hZpTMuEleRecoRatio->GetYaxis()->SetTitleSize(0.045);
hZpTMuEleRecoRatio->GetYaxis()->SetRangeUser(0.5,1.5);
hZpTMuEleRecoRatio->GetYaxis()->CenterTitle();
hZpTMuEleRecoRatio->Draw("E");
gStyle->SetStatStyle(0);
cZpTMuEleRecoRatio->SaveAs( (plotDirectoryPath + cZpTMuEleRecoRatio->GetName() + suffix + plotFileExtension).c_str() );
}
TH1D *hZpTMuEleGenRatio;
TCanvas *cZpTMuEleGenRatio;
if (hZpTgen[0]->GetNbinsX() == hZpTgen[1]->GetNbinsX()) {
bins = hZpTgen[0]->GetNbinsX();
firstEdge = hZpTgen[0]->GetXaxis()->GetBinLowEdge(1);
lastEdge = hZpTgen[0]->GetXaxis()->GetBinUpEdge(bins);
hZpTMuEleGenRatio = new TH1D("hZpTMuEleGenRatio","",bins,firstEdge,lastEdge);
if (!hZpTMuEleGenRatio->Divide(hZpTgen[0],hZpTgen[1])) cout <
|
请发表评论