TCutG* GetErrorBand(string name, int N, double* Mass, double* Low, double* High)
{
TCutG* cutg = new TCutG(name.c_str(),2*N+2);
cutg->SetFillColor(kGreen-7);
cutg->SetLineStyle(0);
cutg->SetLineColor(0);
int I = 0;
for(int i=0;i<N;i++){cutg->SetPoint(I,Mass[i] , Low[i] );I++; }
cutg->SetPoint(I,Mass[N-1] , Low[N-1] );I++;
cutg->SetPoint(I,Mass[N-1] , High[N-1] );I++;
for(int i=0;i<N;i++){cutg->SetPoint(I,Mass[N-1-i], High[N-1-i]);I++;}
return cutg;
}
TCutG* GetErrorBand(string name, TGraph* Low, TGraph* High)
{
TCutG* cutg = new TCutG(name.c_str(),Low->GetN()+High->GetN()+2);
cutg->SetFillColor(kGreen-7);
cutg->SetLineStyle(0);
cutg->SetLineColor(0);
int I = 0;
for(int i=0;i<Low->GetN();i++){ cutg->SetPoint(I,Low ->GetX()[i] , Low ->GetY()[i] );I++; }
cutg->SetPoint(I,Low ->GetX()[Low ->GetN()-1] , Low ->GetY()[Low ->GetN()-1] );I++;
cutg->SetPoint(I,High->GetX()[High->GetN()-1] , High->GetY()[High->GetN()-1] );I++;
for(int i=0;i<High->GetN() ;i++){cutg->SetPoint(I,High->GetX()[High->GetN()-1-i], High->GetY()[High->GetN()-1-i]);I++;}
return cutg;
}
/*===========================================================================*/
void d2dwin(Char_t *fnam, Char_t *wnam, Int_t n, Char_t *c)
{
c1=(TCanvas *)gROOT->FindObject(c);
c1->cd(n);
TCutG *mycutg;
TFile *fcutg;
fcutg = new TFile(fnam);
mycutg=(TCutG *)fcutg->Get(wnam);
fcutg->Close();
if (mycutg != NULL)
{
mycutg->SetLineColor(2);
mycutg->Draw();
}
f1 = new TFile(rootfilename,"READ");
return;
}
