{

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[30];
TLine* a1[30];
TLine* a2[20];

TH1F* tmph;

//c1->SetFillColor(0);
//a1
c1->Print("tmpps.ps[");
for(int i=0; i<22; i++){
  int a=i%8;
  if(!a){
    c1->Clear();
    c1->Divide(2,4);
  }
  c1->cd(a+1);
  T->Draw(Form("L.a1.a_c[%d]>>La1_%d(50,10,300)",i,i),Form("abs(L.a1.t[%d]-1500)<300",i),"",50000);
  a0[i] = new TLine(100,0.,100,((TH1F*)gROOT->FindObject(Form("La1_%d",i)))->GetMaximum());
  a0[i]->SetLineColor(2);
  a0[i]->Draw();
  c1->Update();

  if(a==7||i==21){
    c1->Print("tmpps.ps");
  }
}
//c1->Print("La1.ps]");
//a2
//c1->Print("La2.ps[");
for(int i=0; i<26; i++){
  int a=i%8;
  if(!a){
    c1->Clear();
    c1->Divide(2,4);
  }
  c1->cd(a+1);
  T->Draw(Form("L.a2.a_c[%d]>>La2_%d(50,10,200)",i,i),Form("abs(L.a2.t[%d]-1500)<300",i),"",50000);
  a1[i] = new TLine(60,0.,60,((TH1F*)gROOT->FindObject(Form("La2_%d",i)))->GetMaximum());
  a1[i]->SetLineColor(4);
  a1[i]->Draw();
  c1->Update();

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

// gas cherenkov

for(int i=0; i<10; i++){
  int a=i%10;
  if(!a){
    c1->Clear();
    c1->Divide(2,5);
  }
  c1->cd(a+1);
  T->Draw(Form("R.cer.a_c[%d]>>Rcer_%d(50,10,500)",i,i),Form("abs(R.cer.t[%d]-1450)<50",i),"",50000);
  a2[i] = new TLine(100,0.,100,((TH1F*)gROOT->FindObject(Form("Rcer_%d",i)))->GetMaximum());
  a2[i]->SetLineColor(3);
  a2[i]->Draw();
  c1->Update();

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

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

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