{
  gROOT->Reset();

  TChain chain("h3");
  chain.Add("copy_elena_1.root"); 
  chain.Add("copy_elena_2.root");
  chain.Add("copy_elena_3.root");
  chain.Add("copy_elena_4.root");
  chain.Add("copy_elena_5.root");
  chain.Add("copy_elena_6.root");
  chain.Add("copy_elena_7.root");
  chain.Add("copy_elena_8.root");
  chain.Add("copy_elena_9.root");
  chain.Add("copy_elena_10.root");
  chain.Add("copy_elena_11.root");
  chain.Add("copy_elena_12.root");

//Reset ROOT and connect tree file
//   gROOT->Reset();
//   TFile *f = (TFile*)gROOT->GetListOfFiles()->FindObject("copy_elena_9.root");
//   if (!f) {
//      f = new TFile("copy_elena_9.root");
//   }
//   TTree *h3 = (TTree*)gDirectory->Get("h3");

//Declaration of leaves types
   UInt_t          Gt;
   UInt_t          Ymd;
   UInt_t          Hms;
   UInt_t          Mmn;
   UShort_t        Irun;
   UChar_t         Ifil;
   UInt_t          Ieve;
   Float_t         T200;
   Float_t         P200;
   UShort_t        Iact;
   UShort_t        Idat;
   UShort_t        Itrg;
   Float_t         Azm;
   Float_t         Zem;
   Float_t         Azs;
   Float_t         Zes;
   Float_t         Ra;
   Float_t         Dec;
   Float_t         Gall;
   Float_t         Galb;
   Float_t         Rag;
   Float_t         Decg;
   Float_t         Gallg;
   Float_t         Galbg;
   Int_t           Idev;
   UShort_t        Msch;
   UShort_t        Mscd;
   Float_t         Az;
   Float_t         Ze;
   Float_t         Xc;
   Float_t         Yc;
   Float_t         Size;
   Float_t         Age;
   Float_t         Nmu;
   Float_t         T0;
   Float_t         Lmuo;
   Float_t         Lnel;
   Float_t         Size2;
   Float_t         Age2;
   UInt_t          Hit7;
   UInt_t          Hit4;
   Int_t           Ghp1;
   UChar_t         Ghp2;
   Char_t          Ngrs;
   UChar_t         Novg;
   Char_t          Idmx;
   Float_t         Edmx;
   Float_t         Azg;
   Float_t         Zeg;
   Float_t         T0g;
   Float_t         Xcg;
   Float_t         Ycg;
   Float_t         Sizeg;
   Float_t         Ageg;
   Float_t         Sizmg;
   UInt_t          Nctot;
   UInt_t          Nmtot;
   UChar_t         Ilev;
   Float_t         Dnext;
   Char_t          Nflg;
   Float_t         Gred[37];
   Float_t         Grat1[37];
   UShort_t        Ibgr;
   Int_t           Nbgr;
   UShort_t        Npgr[500];
   Float_t         Argr[500];
   UShort_t        Ibpd;
   Int_t           Nbpd;
   UShort_t        Nppd[500];
   Float_t         Arpd[500];
   UShort_t        Ibed;
   Int_t           Nbed;
   UShort_t        Nped[500];
   Float_t         Ared[500];
   UChar_t         Fanka;

   //   TH1F *hangle = new TH1F("hangle","Zenit Angle",100,0,1);
   //TH2F *hscatt = new TH2F("hscatt","Scatter plot",100,3.2,7.5,100,3,9);
   TH1F *hdistn = new TH1F("hdistr","",100,0.,0.25);

   // Set branch addresses.
   chain.SetBranchAddress("Gt",&Gt);
   chain.SetBranchAddress("Ymd",&Ymd);
   chain.SetBranchAddress("Hms",&Hms);
   chain.SetBranchAddress("Mmn",&Mmn);
   chain.SetBranchAddress("Irun",&Irun);
   chain.SetBranchAddress("Ifil",&Ifil);
   chain.SetBranchAddress("Ieve",&Ieve);
   chain.SetBranchAddress("T200",&T200);
   chain.SetBranchAddress("P200",&P200);
   chain.SetBranchAddress("Iact",&Iact);
   chain.SetBranchAddress("Idat",&Idat);
   chain.SetBranchAddress("Itrg",&Itrg);
   chain.SetBranchAddress("Azm",&Azm);
   chain.SetBranchAddress("Zem",&Zem);
   chain.SetBranchAddress("Azs",&Azs);
   chain.SetBranchAddress("Zes",&Zes);
   chain.SetBranchAddress("Ra",&Ra);
   chain.SetBranchAddress("Dec",&Dec);
   chain.SetBranchAddress("Gall",&Gall);
   chain.SetBranchAddress("Galb",&Galb);
   chain.SetBranchAddress("Rag",&Rag);
   chain.SetBranchAddress("Decg",&Decg);
   chain.SetBranchAddress("Gallg",&Gallg);
   chain.SetBranchAddress("Galbg",&Galbg);
   chain.SetBranchAddress("Idev",&Idev);
   chain.SetBranchAddress("Msch",&Msch);
   chain.SetBranchAddress("Mscd",&Mscd);
   chain.SetBranchAddress("Az",&Az);
   chain.SetBranchAddress("Ze",&Ze);
   chain.SetBranchAddress("Xc",&Xc);
   chain.SetBranchAddress("Yc",&Yc);
   chain.SetBranchAddress("Size",&Size);
   chain.SetBranchAddress("Age",&Age);
   chain.SetBranchAddress("Nmu",&Nmu);
   chain.SetBranchAddress("T0",&T0);
   chain.SetBranchAddress("Lmuo",&Lmuo);
   chain.SetBranchAddress("Lnel",&Lnel);
   chain.SetBranchAddress("Size2",&Size2);
   chain.SetBranchAddress("Age2",&Age2);
   chain.SetBranchAddress("Hit7",&Hit7);
   chain.SetBranchAddress("Hit4",&Hit4);
   chain.SetBranchAddress("Ghp1",&Ghp1);
   chain.SetBranchAddress("Ghp2",&Ghp2);
   chain.SetBranchAddress("Ngrs",&Ngrs);
   chain.SetBranchAddress("Novg",&Novg);
   chain.SetBranchAddress("Idmx",&Idmx);
   chain.SetBranchAddress("Edmx",&Edmx);
   chain.SetBranchAddress("Azg",&Azg);
   chain.SetBranchAddress("Zeg",&Zeg);
   chain.SetBranchAddress("T0g",&T0g);
   chain.SetBranchAddress("Xcg",&Xcg);
   chain.SetBranchAddress("Ycg",&Ycg);
   chain.SetBranchAddress("Sizeg",&Sizeg);
   chain.SetBranchAddress("Ageg",&Ageg);
   chain.SetBranchAddress("Sizmg",&Sizmg);
   chain.SetBranchAddress("Nctot",&Nctot);
   chain.SetBranchAddress("Nmtot",&Nmtot);
   chain.SetBranchAddress("Ilev",&Ilev);
   chain.SetBranchAddress("Dnext",&Dnext);
   chain.SetBranchAddress("Nflg",&Nflg);
   chain.SetBranchAddress("Gred",Gred);
   chain.SetBranchAddress("Grat1",Grat1);
   chain.SetBranchAddress("Ibgr",&Ibgr);
   chain.SetBranchAddress("Nbgr",&Nbgr);
   chain.SetBranchAddress("Npgr",Npgr);
   chain.SetBranchAddress("Argr",Argr);
   chain.SetBranchAddress("Ibpd",&Ibpd);
   chain.SetBranchAddress("Nbpd",&Nbpd);
   chain.SetBranchAddress("Nppd",Nppd);
   chain.SetBranchAddress("Arpd",Arpd);
   chain.SetBranchAddress("Ibed",&Ibed);
   chain.SetBranchAddress("Nbed",&Nbed);
   chain.SetBranchAddress("Nped",Nped);
   chain.SetBranchAddress("Ared",Ared);
   chain.SetBranchAddress("Fanka",&Fanka);

//     This is the loop skeleton
//       To read only selected branches, Insert statements like:
// h3->SetBranchStatus("*",0);  // disable all branches
// TTreePlayer->SetBranchStatus("branchname",1);  // activate branchname

   Long64_t nentries = chain.GetEntries();

   Long64_t nbytes = 0;
   for (Long64_t i=0; i<nentries;i++) {
     
     nbytes += chain.GetEntry(i);
   
     if (Fanka<4&&Ngrs>19&&Nflg>0&&Ageg>=0.4&&Ageg<=1.4&&Xcg>=-460&&Xcg<=-45&&Ycg>=-550&&Ycg<=20&&Zeg>0&&Zeg<0.43&&Sizeg>pow(10,7)&&Sizmg>pow(10,3.1))
       {
         
	 //         hangle->Fill(Zeg);
         //hscatt->Fill(log10(Sizmg),log10(Sizeg));
         hdistn->Fill(Sizmg/Sizeg);
       }
     
   }

   //hscatt->Draw();
   hdistn->Draw();
}