{

int nrun;
cout<<"Enter run number: ";
cin>>nrun;
char fname[50];
sprintf(fname,"./rootfiles/cere_%d.root",nrun);
char command[50];
sprintf(command,"test -e %s",fname);
int is_exist=system(command);
if(is_exist!=0)  { 
  cout<<"The file you want to analyze does not exist,"<<endl
      <<"please enter the same run number here to get the histograms you need"<<endl;
  return;
}

TFile *f=new TFile(fname);

TCanvas* c1 = new TCanvas("c1","",600,800);
TLine* a0[40];
TLine* a1[40];

TH1F* tmph;

//c1->SetFillColor(0);
//a1
c1->Print("tmpps.ps[");
for(int i=0; i<34; i++){
  int a=i%8;
  if(!a){
    c1->Clear();
    c1->Divide(2,4);
  }
  c1->cd(a+1);
  T->Draw(Form("L.prl1.a_c[%d]>>Lprl1_%d(50,100,500)",i,i),Form("abs(L.prl1.x[%d]-%f)<0.01",i,1.2-(i%17)*0.15),"",500000);
  a0[i] = new TLine(200,0.,200,((TH1F*)gROOT->FindObject(Form("Lprl1_%d",i)))->GetMaximum());
  a0[i]->SetLineColor(2);
  a0[i]->Draw();
  c1->Update();

  if(a==7||i==33){
    c1->Print("tmpps.ps");
  }
}

//a2
for(int i=0; i<34; i++){
  int a=i%8;
  if(!a){
    c1->Clear();
    c1->Divide(2,4);
  }
  c1->cd(a+1);
  T->Draw(Form("L.prl2.a_c[%d]>>Lprl2_%d(50,100,500)",i,i),Form("abs(L.prl2.x[%d]-%f)<0.01",i,1.15-(i%17)*0.15),"",500000);
  a1[i] = new TLine(200,0.,200,((TH1F*)gROOT->FindObject(Form("Lprl2_%d",i)))->GetMaximum());
  a1[i]->SetLineColor(4);
  a1[i]->Draw();
  c1->Update();

  if(a==7||i==33){
    c1->Print("tmpps.ps");
  }
}

c1->Print("tmpps.ps]");

char command[256];
sprintf(command,"mv tmpps.ps prl_%d.ps",nrun);
system(command);
}