{
  EdbShowerP* dd;
  EdbDataProc     *dset;
  dset = new EdbDataProc("lnk_all.def");
  dset->InitVolume(0, "");        // Read in (all) BT from cp.root (always).
  gAli= new EdbPVRec();
  gAli = dset->PVR();
  
  int npat = gAli->Npatterns();      //number of plates

  
  dd = new EdbShowerP(12);

  for (Int_t nr=0; nr<npat; ++nr) {
    for (Int_t n=0; n<gAli->GetPattern(nr)->GetN(); ++n) {
      e1=gAli->GetPattern(nr)->GetSegment(n);
      e1->PrintNice();
      dd->AddSegment(e1);
    }
  }
 cout << "_________________"<< endl << endl;
  
  
  dd->Update();
  dd->PrintBasics();

   TFile* file2= new TFile("TESTFILE2.root","RECREATE");
   cout << "file created"<<endl;
   dd->Write();
   cout << "write done...."<<endl;
   file2->Close();
   cout << "file closed"<<endl;
  
   return;
}