void
ptdrHistogramStyleAsciiExample(Int_t varIndex=0,Int_t nBins=30,Double_t min=0,Double_t max=10)
{

  if(gSystem->AccessPathName("ptdrHistogramStyle_C.so")!=0){
    cout << "The shared library ptdrHistogramStyle_C.so not found." << endl
	 << "To get it from ptdrHistogramStyle.C, type in ROOT \".L ptdrHistogramStyle.C++\"" << endl;
    gSystem->Exit(0);
  }
  gSystem->Load("ptdrHistogramStyle_C.so");
  gSystem->Load("ptdrHistogramStyleAscii_C.so");

  TString varname;
  if (varIndex==0) varname="btag1";
  if (varIndex==1) varname="btag2";
  if (varIndex==2) varname="ptj1";
  if (varIndex==3) varname="ptj2";
  if (varIndex==4) varname="m12";
  if (varIndex==5) varname="eta2";
  if (varIndex==6) varname="mtop";
  if (varIndex==7) varname="pttop";
  if (varIndex==8) varname="sigmat";
  if (varIndex==9) varname="ht";
  if (varIndex==10) varname="mtb";
  //  Int_t nBins=10;
  //  Double_t min=0;
  //  Double_t max=10;
  TPTDRHistogramStyleAscii* hs = new TPTDRHistogramStyleAscii;
  //  Int_t varIndex=0;

  /** From ptdrHistogramStyle.h :
      TPTDRHistogramStyleAscii::Add(TString variableName, TString fileName, TString lName,
			    Int_t nBins=100, Double_t min=0, Double_t max=0, 
			    Double_t weight=1, Int_t variableIndex);

      reads variableName from the treeName from the fileName and
      puts it to the array of histograms. lName is legend name for the
      histogram. nBins, min and max define the number of bins, 
      maximal and minimal values for the X-axis(range). weight can be
      used when stacking histograms with different weights. */
  hs->AddHistogram(varname,"schannel_4_4_0.dat","s-channel",nBins,min,max,0.117,varIndex);
  hs->AddHistogram(varname,"tchannel_4_4_0.dat","t-channel",nBins,min,max,1.02,varIndex);
  hs->AddHistogram(varname,"ttbar_4_4_0.dat","tt",nBins,min,max,1.82,varIndex);
  hs->AddHistogram(varname,"wt2l_4_4_0.dat","Wt,2l",nBins,min,max,0.66,varIndex);
  hs->AddHistogram(varname,"wt1la_4_4_0.dat","Wt,1l",nBins,min,max,4*1.66,varIndex);
  hs->AddHistogram(varname,"wbb_4_4_0.dat","Wbb",nBins,min,max,1.26,varIndex);
  hs->AddHistogram(varname,"wjj_4_4_0.dat","Wjj",nBins,min,max,54,varIndex);
  
  //  varname="ptj1";
  //  hs->AddHistogram(varname,"schannel.root","sPtj2",nBins,min,max);
  //hs->AddHistogram(varname,"tchannel.root","tPtj2",nBins,min,max);

  cout << "Plotting " << hs->GetHists()->GetSize() << " histograms..." << endl;
  gStyle->SetCanvasColor(kWhite);
  TCanvas * canv = new TCanvas(varname," ",0,0,600,400);
  canv->SetFrameFillColor(kWhite);

  hs->DrawHistograms();
  hs->DrawLegend();
  TString fileps="SCH_";
  fileps.Append(varname);
  fileps.Append(".eps");
  TString filegif=fileps;
  filegif.Append(".gif");
  canv->Print(fileps);
  canv->Print(filegif);
}