[u/mrichter/AliRoot.git] / PWGLF / FORWARD / analysis2 / qa / DrawMCResult.C

/**
 * @file   DrawMCResult.C
 * @author Christian Holm Christensen <cholm@dalsgaard.hehi.nbi.dk>
 * @date   Thu Jul  7 10:57:01 2011
 * 
 * @brief  Script to draw steps (deprecated version - use DrawSteps.C)
 * 
 * 
 * @deprecated Use DrawSteps.C instead
 * @ingroup pwglf_forward_scripts_qa
 */
/** 
 * 
 * 
 * @param forward Forward folder
 * @param sub     Sub-folder to get
 * @param name    Name of stack
 * 
 * @return Stack, or null
 *
 * @deprecated Use DrawSteps.C instead
 * @ingroup pwglf_forward_scripts_qa
 */
THStack*
GetStack(const TList& forward,  const char* sub, const char* name)
{
  TList* lsub = static_cast<TList*>(forward.FindObject(sub));
  if (!lsub) { 
    Warning("GetStack", "Sub list %s not found in %s", sub, forward.GetName());
    return 0;
  }
  THStack* ret = static_cast<THStack*>(lsub->FindObject(name));
  if (!ret) 
    Warning("GetStack" "Stack %s not found in %s", name, sub);
  return ret;
}

/** 
 * Rebin a histogram
 * 
 * @param h      Histogram
 * @param rebin  Rebinning factor
 * 
 * @return Pointer to histogram
 *
 * @deprecated Use DrawSteps.C instead
 * @ingroup pwglf_forward_scripts_qa
 */
TH1* 
Rebin(TH1* h, Int_t rebin)
{
  if (rebin <= 1) return h;
  h->Rebin(rebin);
  h->Scale(1. / rebin);
  return h;
}

/** 
 * Take ratio of two histograms
 * 
 * @param h1  Numerator
 * @param h2  Denominator
 * 
 * @return Newly allocated histogram containg ratio
 *
 * @deprecated Use DrawSteps.C instead
 * @ingroup pwglf_forward_scripts_qa
 */
TH1*
Ratio(const TH1* h1, const TH1* h2)
{
  if (!h1) return;
  if (!h2) return;
  
  TH1* copy = static_cast<TH1*>(h2->Clone("tmp"));
  copy->SetName(Form("%s_%s", h2->GetName(), h1->GetName()));
  copy->SetTitle(Form("%s/%s", h2->GetTitle(), h1->GetTitle()));
  copy->SetDirectory(0);
  copy->Divide(h1);

  return copy;
}
/** 
 * Take ratio of histograms in stacks
 * 
 * @param r   Return stack of ratios
 * @param h1  Numerators
 * @param h2  Denominators
 * 
 * @return How many histograms in the return stack
 *
 * @deprecated Use DrawSteps.C instead
 * @ingroup pwglf_forward_scripts_qa
 */
Int_t 
Ratio(THStack* r, const THStack* h1, const THStack* h2)
{
  if (!h1) return 0;
  if (!h2) return 0;

  int n1 = h1->GetHists()->GetEntries();
  int n2 = h2->GetHists()->GetEntries();
  int nH = 0;
  for (int i = 0; i < n1 && i < n2; i++) { 
    TH1* hh1 = static_cast<TH1*>(h1->GetHists()->At(i));
    TH1* hh2 = static_cast<TH1*>(h2->GetHists()->At(i));
    TH1* h   = Ratio(hh1, hh2);
    if (!h) continue;
    nH++;
    r->Add(h);
  }
  return nH;
}
/** 
 * Draw MC results
 * 
 * @param filename  Input file name
 * @param rebin     Rebinning factor
 * @param ratios    Whether to show ratios
 *
 * @deprecated Use DrawSteps.C instead
 * @ingroup pwglf_forward_scripts_qa
 */
void
DrawMCResult(const char* filename="forward.root", Int_t rebin=1,
	     Bool_t ratios=true)
{
  gStyle->SetPalette(1);
  gStyle->SetOptFit(0);
  gStyle->SetOptStat(0);

  TFile* file = TFile::Open(filename, "READ");
  if (!file) { 
    Error("DrawMCResult", "failed to open %s", filename);
    return;
  }

  TList* forward = static_cast<TList*>(file->Get("Forward"));
  if (!forward) { 
    Error("DrawMCResult", "List Forward not found in %s", filename);
    return;
  }
  THStack* res    = GetStack(*forward, "ringResults", "all");
  THStack* mcRes  = GetStack(*forward, "mcRingResults", "all");
  THStack* deltas = GetStack(*forward, "fmdSharingFilter", "sums");
  THStack* nchs   = GetStack(*forward, "fmdDensityCalculator", "sums");
  THStack* prims  = GetStack(*forward, "fmdCorrector", "sums");
  
  TH1* sumEta = static_cast<TH1*>(forward->FindObject("mcSumEta"));
  if (!sumEta) { 
    Warning("DrawMCResults", "mcSumEta not found in Forward");
  }


  gStyle->SetTitleBorderSize(0);
  gStyle->SetTitleFillColor(0);
  gStyle->SetTitleStyle(0);
  // gStyle->SetTitleColor(kBlack);


  TCanvas* c = new TCanvas("c", "C", 900, 700);
  c->SetFillColor(0);
  c->SetBorderSize(0);
  c->SetTopMargin(0.05);
  c->SetRightMargin(0.05);

  Double_t y1 = (ratios ? .3 : 0);
  TPad* p1 = new TPad("p1", "p1", 0, y1, 1, 1, 0, 0, 0); 
  p1->SetBottomMargin(ratios ? 0.01 : .10);
  p1->SetFillColor(0);
  p1->SetBorderSize(0);
  p1->SetTopMargin(0.05);
  p1->SetRightMargin(0.05);
  p1->Draw();
  p1->cd();

  THStack* all = new THStack("all", "Analysis steps");
  if (res) {
    res->SetTitle("dN_{ch}/d#eta");
    TIter next(res->GetHists());
    TH1* h = 0;
    while ((h = static_cast<TH1*>(next()))) all->Add(Rebin(h,rebin));
  }
  if (mcRes) {
    mcRes->SetTitle("Track-Refs");
    TIter next(mcRes->GetHists());
    TH1* h = 0;
    while ((h = static_cast<TH1*>(next()))) all->Add(Rebin(h,rebin));
  }
  if (deltas) {
    deltas->SetTitle("#sum_{} #Delta/#Delta_{mip}");
    TIter next(deltas->GetHists());
    TH1* h = 0;
    while ((h = static_cast<TH1*>(next()))) { 
      h->SetMarkerStyle(25);
      all->Add(Rebin(h,rebin));
    }
  }
  if (nchs) {
    nchs->SetTitle("#sum_{} N_{ch,incl}");
    TIter next(nchs->GetHists());
    TH1* h = 0;
    while ((h = static_cast<TH1*>(next()))) { 
      h->SetMarkerStyle(21);
      all->Add(Rebin(h,rebin));
    }
  }
  if (prims) {
    prims->SetTitle("#sum_{} N_{ch,primary}");
    TIter next(prims->GetHists());
    TH1* h = 0;
    while ((h = static_cast<TH1*>(next()))) { 
      h->SetMarkerStyle(22);
      all->Add(Rebin(h,rebin));
    }
  }
  if (sumEta) all->Add(sumEta);
  all->Draw("nostack");
  all->GetHistogram()->SetXTitle("#eta");
  all->GetHistogram()->SetYTitle("signal");
  all->GetHistogram()->GetXaxis()->SetLabelFont(132);
  all->GetHistogram()->GetXaxis()->SetTitleFont(132);
  all->GetHistogram()->GetYaxis()->SetLabelFont(132);
  all->GetHistogram()->GetYaxis()->SetTitleFont(132);
  c->SetGridx();

  TLegend* l = new TLegend(.33, .2, .53, .9);
  TLegendEntry* e = 0;
  l->SetFillColor(0);
  l->SetFillStyle(0);
  l->SetBorderSize(0);
  l->SetNColumns(1);
  l->SetTextFont(132);
  Int_t i = 0;
  if (res) { 
    TIter next(res->GetHists());
    TH1* h = 0;
    while ((h = static_cast<TH1*>(next()))) {
      e = l->AddEntry(Form("dummy%02d", i++),h->GetTitle(),"pl");
      e->SetMarkerStyle(20);
      e->SetMarkerColor(h->GetMarkerColor());
    }
    e = l->AddEntry(Form("dummy%02d", i++),res->GetTitle(),"pl");
    e->SetMarkerStyle(20);
    e->SetMarkerColor(kBlack);
  }
  if (deltas) { 
    e = l->AddEntry(Form("dummy%02d", i++), deltas->GetTitle(),"pl");
    TH1* h = static_cast<TH1*>(deltas->GetHists()->At(0));
    e->SetMarkerStyle(h->GetMarkerStyle());
    e->SetMarkerColor(kBlack);
  }
  if (nchs) { 
    e = l->AddEntry(Form("dummy%02d",i++),nchs->GetTitle(),"pl");
    TH1* h = static_cast<TH1*>(nchs->GetHists()->At(0));
    e->SetMarkerStyle(h->GetMarkerStyle());
    e->SetMarkerColor(kBlack);
  }
  if (prims) { 
    e = l->AddEntry(Form("dummy%02d", i++), prims->GetTitle(),"pl");
    TH1* h = static_cast<TH1*>(prims->GetHists()->At(0));
    e->SetMarkerStyle(h->GetMarkerStyle());
    e->SetMarkerColor(kBlack);
  }

  if (mcRes) { 
    e = l->AddEntry(Form("dummy%02d", i++), mcRes->GetTitle(), "pl");
    TH1* h = static_cast<TH1*>(mcRes->GetHists()->At(0));
    e->SetMarkerStyle(h->GetMarkerStyle());
    e->SetMarkerColor(kBlack);
  }

  if (sumEta) l->AddEntry(sumEta);
  l->Draw();


  if (!ratios) return;

  c->cd();
  TPad* p2 = new TPad("p2", "p2", 0, 0, 1, y1, 0, 0, 0); 
  p2->SetTopMargin(0);
  p2->SetFillColor(0);
  p2->SetBorderSize(0);
  p2->SetRightMargin(0.05);
  p2->Draw();
  p2->cd();

  THStack* rs = new THStack("ratios", "Ratios");
  Int_t nDN = Ratio(rs, deltas, nchs);
  Int_t nNR = Ratio(rs, nchs, res);
  Int_t nRP = Ratio(rs, res, prims);
  rs->Draw("nostack");

  TLegend* ll = new TLegend(.38, .2, .48, .9);
  ll->SetFillColor(0);
  ll->SetFillStyle(0);
  ll->SetBorderSize(0);
  ll->SetNColumns(1);
  ll->SetTextFont(132);
  if (nDN) {
    e = ll->AddEntry("d1", Form("#frac{%s}{%s}", 
					      nchs->GetTitle(), 
					      deltas->GetTitle()), "lp");
    e->SetMarkerStyle(21);
  }
  if (nNR) {
    e = ll->AddEntry("d2", Form("#frac{%s}{%s}", res->GetTitle(), 
				nchs->GetTitle()), "lp");
    e->SetMarkerStyle(20);
  }
  if (nRP) {
    e = ll->AddEntry("d3", Form("#frac{%s}{%s}", prims->GetTitle(), 
			    res->GetTitle()), "lp");
    e->SetMarkerStyle(22);
  }
  ll->Draw();
  

}
 
//
// EOF
//
Commit	Line	Data
56199f2b	1	/**
	2	* @file DrawMCResult.C
	3	* @author Christian Holm Christensen <cholm@dalsgaard.hehi.nbi.dk>
	4	* @date Thu Jul 7 10:57:01 2011
	5	*
	6	* @brief Script to draw steps (deprecated version - use DrawSteps.C)
	7	*
	8	*
2dbde04b	9	* @deprecated Use DrawSteps.C instead
bd6f5206	10	* @ingroup pwglf_forward_scripts_qa
56199f2b	11	*/
	12	/**
	13	*
	14	*
2dbde04b	15	* @param forward Forward folder
	16	* @param sub Sub-folder to get
	17	* @param name Name of stack
56199f2b	18	*
2dbde04b	19	* @return Stack, or null
56199f2b	20	*
2dbde04b	21	* @deprecated Use DrawSteps.C instead
bd6f5206	22	* @ingroup pwglf_forward_scripts_qa
56199f2b	23	*/
0919aa56	24	THStack*
	25	GetStack(const TList& forward, const char* sub, const char* name)
	26	{
	27	TList* lsub = static_cast<TList*>(forward.FindObject(sub));
	28	if (!lsub) {
	29	Warning("GetStack", "Sub list %s not found in %s", sub, forward.GetName());
	30	return 0;
	31	}
	32	THStack* ret = static_cast<THStack*>(lsub->FindObject(name));
	33	if (!ret)
	34	Warning("GetStack" "Stack %s not found in %s", name, sub);
	35	return ret;
	36	}
	37
56199f2b	38	/**
2dbde04b	39	* Rebin a histogram
56199f2b	40	*
2dbde04b	41	* @param h Histogram
2dbde04b	42	* @param rebin Rebinning factor
56199f2b	43	*
2dbde04b	44	* @return Pointer to histogram
56199f2b	45	*
2dbde04b	46	* @deprecated Use DrawSteps.C instead
bd6f5206	47	* @ingroup pwglf_forward_scripts_qa
56199f2b	48	*/
0919aa56	49	TH1*
	50	Rebin(TH1* h, Int_t rebin)
	51	{
	52	if (rebin <= 1) return h;
	53	h->Rebin(rebin);
	54	h->Scale(1. / rebin);
	55	return h;
	56	}
	57
56199f2b	58	/**
2dbde04b	59	* Take ratio of two histograms
56199f2b	60	*
2dbde04b	61	* @param h1 Numerator
2dbde04b	62	* @param h2 Denominator
56199f2b	63	*
2dbde04b	64	* @return Newly allocated histogram containg ratio
56199f2b	65	*
2dbde04b	66	* @deprecated Use DrawSteps.C instead
bd6f5206	67	* @ingroup pwglf_forward_scripts_qa
56199f2b	68	*/
8b1c306c	69	TH1*
	70	Ratio(const TH1* h1, const TH1* h2)
	71	{
	72	if (!h1) return;
	73	if (!h2) return;
	74
	75	TH1* copy = static_cast<TH1*>(h2->Clone("tmp"));
	76	copy->SetName(Form("%s_%s", h2->GetName(), h1->GetName()));
	77	copy->SetTitle(Form("%s/%s", h2->GetTitle(), h1->GetTitle()));
	78	copy->SetDirectory(0);
	79	copy->Divide(h1);
	80
	81	return copy;
	82	}
56199f2b	83	/**
2dbde04b	84	* Take ratio of histograms in stacks
56199f2b	85	*
2dbde04b	86	* @param r Return stack of ratios
	87	* @param h1 Numerators
	88	* @param h2 Denominators
56199f2b	89	*
2dbde04b	90	* @return How many histograms in the return stack
56199f2b	91	*
2dbde04b	92	* @deprecated Use DrawSteps.C instead
bd6f5206	93	* @ingroup pwglf_forward_scripts_qa
56199f2b	94	*/
8b1c306c	95	Int_t
	96	Ratio(THStack* r, const THStack* h1, const THStack* h2)
	97	{
	98	if (!h1) return 0;
	99	if (!h2) return 0;
	100
	101	int n1 = h1->GetHists()->GetEntries();
	102	int n2 = h2->GetHists()->GetEntries();
	103	int nH = 0;
	104	for (int i = 0; i < n1 && i < n2; i++) {
	105	TH1* hh1 = static_cast<TH1*>(h1->GetHists()->At(i));
	106	TH1* hh2 = static_cast<TH1*>(h2->GetHists()->At(i));
	107	TH1* h = Ratio(hh1, hh2);
	108	if (!h) continue;
	109	nH++;
	110	r->Add(h);
	111	}
	112	return nH;
	113	}
56199f2b	114	/**
2dbde04b	115	* Draw MC results
56199f2b	116	*
2dbde04b	117	* @param filename Input file name
	118	* @param rebin Rebinning factor
	119	* @param ratios Whether to show ratios
56199f2b	120	*
2dbde04b	121	* @deprecated Use DrawSteps.C instead
bd6f5206	122	* @ingroup pwglf_forward_scripts_qa
56199f2b	123	*/
0919aa56	124	void
8b1c306c	125	DrawMCResult(const char* filename="forward.root", Int_t rebin=1,
8b1c306c	126	Bool_t ratios=true)
0919aa56	127	{
	128	gStyle->SetPalette(1);
	129	gStyle->SetOptFit(0);
	130	gStyle->SetOptStat(0);
	131
	132	TFile* file = TFile::Open(filename, "READ");
	133	if (!file) {
	134	Error("DrawMCResult", "failed to open %s", filename);
	135	return;
	136	}
	137
	138	TList* forward = static_cast<TList*>(file->Get("Forward"));
	139	if (!forward) {
	140	Error("DrawMCResult", "List Forward not found in %s", filename);
	141	return;
	142	}
	143	THStack* res = GetStack(*forward, "ringResults", "all");
	144	THStack* mcRes = GetStack(*forward, "mcRingResults", "all");
	145	THStack* deltas = GetStack(*forward, "fmdSharingFilter", "sums");
	146	THStack* nchs = GetStack(*forward, "fmdDensityCalculator", "sums");
	147	THStack* prims = GetStack(*forward, "fmdCorrector", "sums");
	148
	149	TH1* sumEta = static_cast<TH1*>(forward->FindObject("mcSumEta"));
	150	if (!sumEta) {
	151	Warning("DrawMCResults", "mcSumEta not found in Forward");
	152	}
	153
	154
8b1c306c	155	gStyle->SetTitleBorderSize(0);
	156	gStyle->SetTitleFillColor(0);
	157	gStyle->SetTitleStyle(0);
	158	// gStyle->SetTitleColor(kBlack);
	159
0919aa56	160
	161	TCanvas* c = new TCanvas("c", "C", 900, 700);
	162	c->SetFillColor(0);
	163	c->SetBorderSize(0);
	164	c->SetTopMargin(0.05);
	165	c->SetRightMargin(0.05);
8b1c306c	166
	167	Double_t y1 = (ratios ? .3 : 0);
	168	TPad* p1 = new TPad("p1", "p1", 0, y1, 1, 1, 0, 0, 0);
	169	p1->SetBottomMargin(ratios ? 0.01 : .10);
	170	p1->SetFillColor(0);
	171	p1->SetBorderSize(0);
	172	p1->SetTopMargin(0.05);
	173	p1->SetRightMargin(0.05);
	174	p1->Draw();
	175	p1->cd();
0919aa56	176
	177	THStack* all = new THStack("all", "Analysis steps");
	178	if (res) {
8b1c306c	179	res->SetTitle("dN_{ch}/d#eta");
0919aa56	180	TIter next(res->GetHists());
	181	TH1* h = 0;
	182	while ((h = static_cast<TH1*>(next()))) all->Add(Rebin(h,rebin));
	183	}
	184	if (mcRes) {
8b1c306c	185	mcRes->SetTitle("Track-Refs");
0919aa56	186	TIter next(mcRes->GetHists());
	187	TH1* h = 0;
	188	while ((h = static_cast<TH1*>(next()))) all->Add(Rebin(h,rebin));
	189	}
	190	if (deltas) {
8b1c306c	191	deltas->SetTitle("#sum_{} #Delta/#Delta_{mip}");
0919aa56	192	TIter next(deltas->GetHists());
	193	TH1* h = 0;
	194	while ((h = static_cast<TH1*>(next()))) {
	195	h->SetMarkerStyle(25);
	196	all->Add(Rebin(h,rebin));
	197	}
	198	}
	199	if (nchs) {
8b1c306c	200	nchs->SetTitle("#sum_{} N_{ch,incl}");
0919aa56	201	TIter next(nchs->GetHists());
	202	TH1* h = 0;
	203	while ((h = static_cast<TH1*>(next()))) {
	204	h->SetMarkerStyle(21);
	205	all->Add(Rebin(h,rebin));
	206	}
	207	}
	208	if (prims) {
8b1c306c	209	prims->SetTitle("#sum_{} N_{ch,primary}");
0919aa56	210	TIter next(prims->GetHists());
	211	TH1* h = 0;
	212	while ((h = static_cast<TH1*>(next()))) {
	213	h->SetMarkerStyle(22);
	214	all->Add(Rebin(h,rebin));
	215	}
	216	}
	217	if (sumEta) all->Add(sumEta);
	218	all->Draw("nostack");
	219	all->GetHistogram()->SetXTitle("#eta");
	220	all->GetHistogram()->SetYTitle("signal");
	221	all->GetHistogram()->GetXaxis()->SetLabelFont(132);
	222	all->GetHistogram()->GetXaxis()->SetTitleFont(132);
	223	all->GetHistogram()->GetYaxis()->SetLabelFont(132);
	224	all->GetHistogram()->GetYaxis()->SetTitleFont(132);
	225	c->SetGridx();
	226
	227	TLegend* l = new TLegend(.33, .2, .53, .9);
8b1c306c	228	TLegendEntry* e = 0;
0919aa56	229	l->SetFillColor(0);
	230	l->SetFillStyle(0);
	231	l->SetBorderSize(0);
	232	l->SetNColumns(1);
	233	l->SetTextFont(132);
	234	Int_t i = 0;
	235	if (res) {
	236	TIter next(res->GetHists());
	237	TH1* h = 0;
	238	while ((h = static_cast<TH1*>(next()))) {
8b1c306c	239	e = l->AddEntry(Form("dummy%02d", i++),h->GetTitle(),"pl");
0919aa56	240	e->SetMarkerStyle(20);
	241	e->SetMarkerColor(h->GetMarkerColor());
	242	}
8b1c306c	243	e = l->AddEntry(Form("dummy%02d", i++),res->GetTitle(),"pl");
	244	e->SetMarkerStyle(20);
	245	e->SetMarkerColor(kBlack);
0919aa56	246	}
0919aa56	247	if (deltas) {
8b1c306c	248	e = l->AddEntry(Form("dummy%02d", i++), deltas->GetTitle(),"pl");
0919aa56	249	TH1* h = static_cast<TH1*>(deltas->GetHists()->At(0));
8b1c306c	250	e->SetMarkerStyle(h->GetMarkerStyle());
8b1c306c	251	e->SetMarkerColor(kBlack);
0919aa56	252	}
0919aa56	253	if (nchs) {
8b1c306c	254	e = l->AddEntry(Form("dummy%02d",i++),nchs->GetTitle(),"pl");
0919aa56	255	TH1* h = static_cast<TH1*>(nchs->GetHists()->At(0));
8b1c306c	256	e->SetMarkerStyle(h->GetMarkerStyle());
8b1c306c	257	e->SetMarkerColor(kBlack);
0919aa56	258	}
0919aa56	259	if (prims) {
8b1c306c	260	e = l->AddEntry(Form("dummy%02d", i++), prims->GetTitle(),"pl");
0919aa56	261	TH1* h = static_cast<TH1*>(prims->GetHists()->At(0));
8b1c306c	262	e->SetMarkerStyle(h->GetMarkerStyle());
8b1c306c	263	e->SetMarkerColor(kBlack);
0919aa56	264	}
	265
	266	if (mcRes) {
8b1c306c	267	e = l->AddEntry(Form("dummy%02d", i++), mcRes->GetTitle(), "pl");
0919aa56	268	TH1* h = static_cast<TH1*>(mcRes->GetHists()->At(0));
8b1c306c	269	e->SetMarkerStyle(h->GetMarkerStyle());
8b1c306c	270	e->SetMarkerColor(kBlack);
0919aa56	271	}
	272
	273	if (sumEta) l->AddEntry(sumEta);
	274	l->Draw();
	275
8b1c306c	276
	277	if (!ratios) return;
	278
	279	c->cd();
	280	TPad* p2 = new TPad("p2", "p2", 0, 0, 1, y1, 0, 0, 0);
	281	p2->SetTopMargin(0);
	282	p2->SetFillColor(0);
	283	p2->SetBorderSize(0);
	284	p2->SetRightMargin(0.05);
	285	p2->Draw();
	286	p2->cd();
	287
	288	THStack* rs = new THStack("ratios", "Ratios");
	289	Int_t nDN = Ratio(rs, deltas, nchs);
	290	Int_t nNR = Ratio(rs, nchs, res);
	291	Int_t nRP = Ratio(rs, res, prims);
	292	rs->Draw("nostack");
	293
	294	TLegend* ll = new TLegend(.38, .2, .48, .9);
	295	ll->SetFillColor(0);
	296	ll->SetFillStyle(0);
	297	ll->SetBorderSize(0);
	298	ll->SetNColumns(1);
	299	ll->SetTextFont(132);
	300	if (nDN) {
	301	e = ll->AddEntry("d1", Form("#frac{%s}{%s}",
	302	nchs->GetTitle(),
	303	deltas->GetTitle()), "lp");
	304	e->SetMarkerStyle(21);
	305	}
	306	if (nNR) {
	307	e = ll->AddEntry("d2", Form("#frac{%s}{%s}", res->GetTitle(),
	308	nchs->GetTitle()), "lp");
	309	e->SetMarkerStyle(20);
	310	}
	311	if (nRP) {
	312	e = ll->AddEntry("d3", Form("#frac{%s}{%s}", prims->GetTitle(),
	313	res->GetTitle()), "lp");
	314	e->SetMarkerStyle(22);
	315	}
	316	ll->Draw();
	317
	318
0919aa56	319	}
0919aa56	320
56199f2b	321	//
	322	// EOF
	323	//