#include <iostream>

#include "grandedq.h"
#include "timecalculator2.h"

#include <TChain.h>
#include <TTree.h>
#include <TROOT.h>
#include <TBranch.h>
#include <TClonesArray.h>
#include <TFile.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TMath.h>

using namespace std;
using namespace kg;

int main(){

//reset ROOT and connect tree file
   gROOT->Reset();

   //   TFile *f = (TFile*)gROOT->GetListOfFiles()->FindObject("SEL8.06470f04-06481f03.root");
   //if (!f) {
   //   f = new TFile("SEL8.06470f04-06481f03.root");
   //}
   //TTree *h3 = (TTree*)gDirectory->Get("h3");

   TChain chain("h3");
   chain.Add("SEL8.04685f01-04691f43.root");
   chain.Add("SEL8.04691f44-04703f04.root");
   chain.Add("SEL8.04703f05-04726f02.root");
   chain.Add("SEL8.04726f03-04736f21.root");
   chain.Add("SEL8.04736f22-04748f07.root");
   chain.Add("SEL8.04748f08-04768f04.root");
   chain.Add("SEL8.04768f05-04792f08.root");
   chain.Add("SEL8.04792f09-04816f09.root");
   chain.Add("SEL8.04816f10-04835f05.root");
   chain.Add("SEL8.04835f06-04872f18.root");
   //chain.Add("SEL8.04872f19-04953f04.root");
   //chain.Add("SEL8.04953f05-05025f08.root");
   //chain.Add("SEL8.05025f09-05140f40.root");
   //chain.Add("SEL8.05140f41-05163f07.root");
   //chain.Add("SEL8.05163f08-05183f46.root");
   //chain.Add("SEL8.05183f47-05189f20.root");
   //chain.Add("SEL8.05189f21-05210f01.root");
   //chain.Add("SEL8.05210f02-05219f10.root");
   //chain.Add("SEL8.05219f11-05231f20.root");
   //chain.Add("SEL8.05231f21-05249f15.root");
   //chain.Add("SEL8.05249f16-05266f36.root");
   //chain.Add("SEL8.05266f37-05283f23.root");
   //chain.Add("SEL8.05283f24-05311f05.root");
   //chain.Add("SEL8.05312f01-05326f07.root");
   //chain.Add("SEL8.05326f08-05349f38.root");
   //chain.Add("SEL8.05349f39-05376f01.root");
   //chain.Add("SEL8.05376f02-05391f22.root");
   //chain.Add("SEL8.05391f23-05403f15.root");
   //chain.Add("SEL8.05403f16-05422f21.root");
   //chain.Add("SEL8.05422f22-05433f40.root");
   //chain.Add("SEL8.05433f41-05449f22.root");
   //chain.Add("SEL8.05449f23-05490f25.root");
   //chain.Add("SEL8.05490f26-05499f11.root");
   //chain.Add("SEL8.05499f12-05514f28.root");
   //chain.Add("SEL8.05514f29-05521f02.root");
   //chain.Add("SEL8.05555f01-05560f07.root");
   //chain.Add("SEL8.05560f08-05569f06.root");
   //chain.Add("SEL8.05569f07-05578f19.root");
   //chain.Add("SEL8.05578f20-05600f04.root");
   //chain.Add("SEL8.05600f05-05634f01.root");
   //chain.Add("SEL8.05634f02-05643f05.root");
   //chain.Add("SEL8.05643f06-05662f06.root");
   //chain.Add("SEL8.05662f07-05675f03.root");
   //chain.Add("SEL8.05675f04-05693f05.root");
   //chain.Add("SEL8.05693f06-05700f02.root");
   //chain.Add("SEL8.05700f03-05705f14.root");
   //chain.Add("SEL8.05705f15-05722f11.root");
   //chain.Add("SEL8.05722f12-05750f18.root");
   //chain.Add("SEL8.05750f19-05780f06.root");
   //chain.Add("SEL8.05780f07-05800f15.root");
   //chain.Add("SEL8.05800f16-05816f12.root");
   //chain.Add("SEL8.05816f13-05853f14.root");
   //chain.Add("SEL8.05853f15-05876f28.root");
   //chain.Add("SEL8.05876f29-05885f23.root");
   //chain.Add("SEL8.05897f01-05904f16.root");
   //chain.Add("SEL8.05904f17-05910f02.root");
   //chain.Add("SEL8.05910f03-05923f31.root");
   //chain.Add("SEL8.05923f32-05936f06.root");
   //chain.Add("SEL8.05936f07-05940f25.root");
   //chain.Add("SEL8.05940f26-05955f07.root");
   //chain.Add("SEL8.05955f08-05972f07.root");
   //chain.Add("SEL8.05972f08-05985f21.root");
   //chain.Add("SEL8.05985f22-05997f09.root");
   //chain.Add("SEL8.05997f10-06005f04.root");
   //chain.Add("SEL8.06005f05-06005f37.root");
   //chain.Add("SEL8.06031f01-06052f07.root");
   //chain.Add("SEL8.06052f08-06069f07.root");
   //chain.Add("SEL8.06069f08-06108f07.root");
   //chain.Add("SEL8.06108f08-06132f04.root");
   //chain.Add("SEL8.06132f05-06176f03.root");
   //chain.Add("SEL8.06176f04-06205f12.root");
   //chain.Add("SEL8.06205f13-06236f05.root");
   //chain.Add("SEL8.06236f06-06255f16.root");
   //chain.Add("SEL8.06255f17-06272f01.root");
   //chain.Add("SEL8.06272f02-06285f12.root");
   //chain.Add("SEL8.06285f13-06308f08.root");
   //chain.Add("SEL8.06308f09-06324f05.root");
   //chain.Add("SEL8.06324f06-06335f17.root");
   //chain.Add("SEL8.06335f18-06345f07.root");
   //chain.Add("SEL8.06345f08-06351f27.root");
   //chain.Add("SEL8.06352f01-06356f09.root");
   //chain.Add("SEL8.06356f10-06363f11.root");
   //chain.Add("SEL8.06363f12-06368f08.root");
   //chain.Add("SEL8.06368f09-06373f16.root");
   //chain.Add("SEL8.06373f17-06380f08.root");
   //chain.Add("SEL8.06380f09-06382f11.root");
   //chain.Add("SEL8.06382f12-06385f18.root");
   //chain.Add("SEL8.06385f19-06388f04.root");
   //chain.Add("SEL8.06388f05-06389f21.root");
   //chain.Add("SEL8.06389f22-06394f02.root");
   //chain.Add("SEL8.06394f03-06411f05.root");
   //chain.Add("SEL8.06411f06-06415f39.root");
   //chain.Add("SEL8.06415f40-06421f05.root");
   //chain.Add("SEL8.06421f06-06424f35.root");
   //chain.Add("SEL8.06424f36-06434f04.root");
   //chain.Add("SEL8.06434f05-06445f03.root");
   //chain.Add("SEL8.06446f01-06456f08.root");
   //chain.Add("SEL8.06456f09-06462f02.root");
   //chain.Add("SEL8.06462f03-06470f03.root");
   //chain.Add("SEL8.06470f04-06481f03.root");

   //TH1F *hangle = new TH1F("hangle","Zenit Angle",100,0,1);   
   //TH2F *hscatt = new TH2F("hscatt","Scatter plot",100,2.,9.,100,2.,9.);
   //TH2F *hxy = new TH2F("hxy","Array",1000,-740.,200.,1000,-800.,300.);

   //TH1F *hdistn1 = new TH1F("hdistn1","n1b1",20,0.,0.15);
   //TH1F *hdistn2 = new TH1F("hdistn2","n2b1",20,0.,0.15);
   //TH1F *hdistn3 = new TH1F("hdistn3","n3b1",20,0.,0.15);

   //TH1F *hdistn1b2 = new TH1F("hdistn1b2","n1b2",28,0.,0.21);
   //TH1F *hdistn2b2 = new TH1F("hdistn2b2","n2b2",28,0.,0.21);
   //TH1F *hdistn3b2 = new TH1F("hdistn3b2","n3b2",28,0.,0.21);

//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;
   UChar_t         Trgi;
   UChar_t         Trga;
   UShort_t        Trgs;
   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         Szl1;
   Float_t         Szl2;
   Float_t         Nmu;
   Float_t         Agemu;
   Short_t         Nedh;
   Short_t         Npdh;
   Short_t         Neds;
   Short_t         Npds;
   UChar_t         Nove;
   UChar_t         Novp;
   Float_t         Slop;
   Float_t         T0;
   Int_t           Imec;
   Float_t         Vdig;
   Float_t         Rmean;
   Float_t         Xcmx;
   Float_t         Ycmx;
   Short_t         Dt0a;
   UChar_t         Clu0;
   Float_t         Lmuo;
   Float_t         Lnel;
   Float_t         Size2;
   Float_t         Age2;
   Int_t           Mchflg;
   Int_t           Mchtrk;
   Int_t           Mchhad;
   Float_t         Mchsar;
   Float_t         Mchcod;
   Float_t         Mchdp6;
   Float_t         Mchdn6;
   Int_t           Mchdma;
   Float_t         Mchdaa;
   Float_t         Mchdca;
   Int_t           Mchdmg;
   Float_t         Mchdag;
   Float_t         Mchdcg;
   Int_t           ntr3;
   Char_t          ifitst[200];
   Float_t         chi2t3[200];
   Float_t         tqt3[200];
   Float_t         xt3[200];
   Float_t         yt3[200];
   Float_t         azt3[200];
   Float_t         zet3[200];
   Float_t         rat3[200];
   Float_t         tat3[200];
   Float_t         hit3[200][16][3];
   UShort_t        patt3[200];
   Float_t         az3m;
   Float_t         ze3m;
   Float_t         Rat3g[200];
   Float_t         Tat3g[200];
   Int_t           ntr2;
   Float_t         xt2[200];
   Float_t         yt2[200];
   Float_t         dxyt2[200];
   Float_t         azt2[200];
   Float_t         zet2[200];
   Float_t         rat2[200];
   Float_t         tat2[200];
   Float_t         hit2[200][16][2];
   UShort_t       patt2[200];
   Float_t         az2m;
   Float_t         ze2m;
   Float_t         Rat2g[200];
   Float_t         Tat2g[200];
   Int_t           Lstflg;
   Int_t           Lsthit;
   Float_t         Lstcd9n;
   Float_t         Lstcd9p;
   Float_t         Lstsd2n;
   Float_t         Lstsd2p;
   UInt_t          Idgr;
   UShort_t        Irng;
   UChar_t         Trgg;
   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         Erat;
   Short_t         Dt0g;
   Float_t         Tdc38;
   Char_t          Ixgr;
   Char_t          Iygr;
   Float_t         Azg;
   Float_t         Zeg;
   Float_t         T0g;
   Float_t         Xcg;
   Float_t         Ycg;
   Float_t         Sizeg;
   Float_t         Ageg;
   Float_t         Sizmg;
   Float_t         R0mg;
   UInt_t          Nctot;
   UInt_t          Nmtot;
   UChar_t         Ilev;
   Float_t         Dnext;
   Float_t         Szold;
   Char_t          Nflg;
   Float_t         Fmeg;
   Float_t         Fmtg;
   Float_t         Fmpg;
   UChar_t         Ndeg;
   UChar_t         Ndtg;
   UChar_t         Ndpg;
   Float_t         Azj;
   Float_t         Zej;
   Float_t         T0j;
   Float_t         Xcj;
   Float_t         Ycj;
   Float_t         Sizej;
   Float_t         Agej;
   Float_t         Sizmj;
   Float_t         R0mj;
   Float_t         Gred[37];
   Float_t         Grat1[37];
   UShort_t        Ibgr;
   Int_t           Nbgr;
   UShort_t        Npgr[500];
   Float_t         Argr[500];
   Float_t         Tigr[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;
   Int_t           Naic;
   Char_t          Trnn;
   Char_t          Naib;
   Char_t          Nesh;
   UShort_t        N16h;
   Char_t          Namh;
   Float_t         T16e;

   //Float_t SGred;
   //SGred = 0.;

   // 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("Trgi",&Trgi);
   chain.SetBranchAddress("Trga",&Trga);
   chain.SetBranchAddress("Trgs",&Trgs);
   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("Szl1",&Szl1);
   chain.SetBranchAddress("Szl2",&Szl2);
   chain.SetBranchAddress("Nmu",&Nmu);
   chain.SetBranchAddress("Agemu",&Agemu);
   chain.SetBranchAddress("Nedh",&Nedh);
   chain.SetBranchAddress("Npdh",&Npdh);
   chain.SetBranchAddress("Neds",&Neds);
   chain.SetBranchAddress("Npds",&Npds);
   chain.SetBranchAddress("Nove",&Nove);
   chain.SetBranchAddress("Novp",&Novp);
   chain.SetBranchAddress("Slop",&Slop);
   chain.SetBranchAddress("T0",&T0);
   chain.SetBranchAddress("Imec",&Imec);
   chain.SetBranchAddress("Vdig",&Vdig);
   chain.SetBranchAddress("Rmean",&Rmean);
   chain.SetBranchAddress("Xcmx",&Xcmx);
   chain.SetBranchAddress("Ycmx",&Ycmx);
   chain.SetBranchAddress("Dt0a",&Dt0a);
   chain.SetBranchAddress("Clu0",&Clu0);
   chain.SetBranchAddress("Lmuo",&Lmuo);
   chain.SetBranchAddress("Lnel",&Lnel);
   chain.SetBranchAddress("Size2",&Size2);
   chain.SetBranchAddress("Age2",&Age2);
   chain.SetBranchAddress("Mchflg",&Mchflg);
   chain.SetBranchAddress("Mchtrk",&Mchtrk);
   chain.SetBranchAddress("Mchhad",&Mchhad);
   chain.SetBranchAddress("Mchsar",&Mchsar);
   chain.SetBranchAddress("Mchcod",&Mchcod);
   chain.SetBranchAddress("Mchdp6",&Mchdp6);
   chain.SetBranchAddress("Mchdn6",&Mchdn6);
   chain.SetBranchAddress("Mchdma",&Mchdma);
   chain.SetBranchAddress("Mchdaa",&Mchdaa);
   chain.SetBranchAddress("Mchdca",&Mchdca);
   chain.SetBranchAddress("Mchdmg",&Mchdmg);
   chain.SetBranchAddress("Mchdag",&Mchdag);
   chain.SetBranchAddress("Mchdcg",&Mchdcg);
   chain.SetBranchAddress("ntr3",&ntr3);
   chain.SetBranchAddress("ifitst",ifitst);
   chain.SetBranchAddress("chi2t3",chi2t3);
   chain.SetBranchAddress("tqt3",tqt3);
   chain.SetBranchAddress("xt3",xt3);
   chain.SetBranchAddress("yt3",yt3);
   chain.SetBranchAddress("azt3",azt3);
   chain.SetBranchAddress("zet3",zet3);
   chain.SetBranchAddress("rat3",rat3);
   chain.SetBranchAddress("tat3",tat3);
   chain.SetBranchAddress("hit3",hit3);
   chain.SetBranchAddress("patt3",patt3);
   chain.SetBranchAddress("az3m",&az3m);
   chain.SetBranchAddress("ze3m",&ze3m);
   chain.SetBranchAddress("Rat3g",Rat3g);
   chain.SetBranchAddress("Tat3g",Tat3g);
   chain.SetBranchAddress("ntr2",&ntr2);
   chain.SetBranchAddress("xt2",xt2);
   chain.SetBranchAddress("yt2",yt2);
   chain.SetBranchAddress("dxyt2",dxyt2);
   chain.SetBranchAddress("azt2",azt2);
   chain.SetBranchAddress("zet2",zet2);
   chain.SetBranchAddress("rat2",rat2);
   chain.SetBranchAddress("tat2",tat2);
   chain.SetBranchAddress("hit2",hit2);
   chain.SetBranchAddress("patt2",patt2);
   chain.SetBranchAddress("az2m",&az2m);
   chain.SetBranchAddress("ze2m",&ze2m);
   chain.SetBranchAddress("Rat2g",Rat2g);
   chain.SetBranchAddress("Tat2g",Tat2g);
   chain.SetBranchAddress("Lstflg",&Lstflg);
   chain.SetBranchAddress("Lsthit",&Lsthit);
   chain.SetBranchAddress("Lstcd9n",&Lstcd9n);
   chain.SetBranchAddress("Lstcd9p",&Lstcd9p);
   chain.SetBranchAddress("Lstsd2n",&Lstsd2n);
   chain.SetBranchAddress("Lstsd2p",&Lstsd2p);
   chain.SetBranchAddress("Idgr",&Idgr);
   chain.SetBranchAddress("Irng",&Irng);
   chain.SetBranchAddress("Trgg",&Trgg);
   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("Erat",&Erat);
   chain.SetBranchAddress("Dt0g",&Dt0g);
   chain.SetBranchAddress("Tdc38",&Tdc38);
   chain.SetBranchAddress("Ixgr",&Ixgr);
   chain.SetBranchAddress("Iygr",&Iygr);
   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("R0mg",&R0mg);
   chain.SetBranchAddress("Nctot",&Nctot);
   chain.SetBranchAddress("Nmtot",&Nmtot);
   chain.SetBranchAddress("Ilev",&Ilev);
   chain.SetBranchAddress("Dnext",&Dnext);
   chain.SetBranchAddress("Szold",&Szold);
   chain.SetBranchAddress("Nflg",&Nflg);
   chain.SetBranchAddress("Fmeg",&Fmeg);
   chain.SetBranchAddress("Fmtg",&Fmtg);
   chain.SetBranchAddress("Fmpg",&Fmpg);
   chain.SetBranchAddress("Ndeg",&Ndeg);
   chain.SetBranchAddress("Ndtg",&Ndtg);
   chain.SetBranchAddress("Ndpg",&Ndpg);
   chain.SetBranchAddress("Azj",&Azj);
   chain.SetBranchAddress("Zej",&Zej);
   chain.SetBranchAddress("T0j",&T0j);
   chain.SetBranchAddress("Xcj",&Xcj);
   chain.SetBranchAddress("Ycj",&Ycj);
   chain.SetBranchAddress("Sizej",&Sizej);
   chain.SetBranchAddress("Agej",&Agej);
   chain.SetBranchAddress("Sizmj",&Sizmj);
   chain.SetBranchAddress("R0mj",&R0mj);
   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("Tigr",Tigr);
   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);
   chain.SetBranchAddress("Naic",&Naic);
   chain.SetBranchAddress("Trnn",&Trnn);
   chain.SetBranchAddress("Naib",&Naib);
   chain.SetBranchAddress("Nesh",&Nesh);
   chain.SetBranchAddress("N16h",&N16h);
   chain.SetBranchAddress("Namh",&Namh);
   chain.SetBranchAddress("T16e",&T16e);

   TTree *qt = chain.CloneTree(0);
   qt->SetMaxTreeSize(1000000000);

   const grandedq DQ("grandedq.db");

   grandedq::cut mycut;
   mycut.standard();//chiedo le stazioni interne
   //mycut.require_station(7).require_station(8).require_station(9).require_station(10).require_station(13).require_station(14).require_station(15).require_station(16).require_station(18).require_station(19).require_station(20).require_station(21).require_station(24).require_station(25).require_station(26).require_station(29).require_station(30).require_station(31).require_station(32); 
                    //chiedo tutte le stazioni
   //   mycut.require_station(1).require_station(2).require_station(3).require_station(4).require_station(5).require_station(6).require_station(7).require_station(8).require_station(9).require_station(10).require_station(11).require_station(12).require_station(13).require_station(14).require_station(15).require_station(16).require_station(17).require_station(18).require_station(19).require_station(20).require_station(21).require_station(22).require_station(23).require_station(24).require_station(25).require_station(26).require_station(27).require_station(28).require_station(29).require_station(30).require_station(31).require_station(32).require_station(33).require_station(34).require_station(35).require_station(36).require_station(37); 
                    //chiedo 18 clusters
   mycut.require_clusters(18);

   timecalculator2 Gtc(120.);//object per calcolo del tempo effettivo

   //ofstream myfile1("muexp.dat",ios::out);
   //ofstream myfile2("neexp.dat",ios::out);

//     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);//prendo ogni singolo evento
     //SGred = 0.;

     //for (int i=0; i<37; i++){
     //  if (Gred[i]>=0){SGred = SGred + Gred[i];}
     //}

     if (DQ.ankaevent(Irun,Ieve)){continue;}//se e` un anka event lo salto
     
     if (not (DQ.getstatus(Irun,Ieve).flags() & grandedq::rawdata_available))
      {continue;}//se non ha rawdata available lo salto
     
     if (Irun<4775){continue;} 
     if(Irun>=5341&&Irun<=5347){continue;}//se ha run bacato lo salto

     //se soddisfa i tagli di sven lo tengo 
     if (DQ.getstatus(Irun,Ieve) & mycut){
      //nuovo tree
      qt->Fill();

      Gtc.put(Gt,Mmn,Ieve,Irun,Ifil);//inizio il conteggio del tempo

      //analisi 
      //tagli fabi + area mia e 10 in formula AC
      //if (Fanka<4&&Nflg>0&&Idmx>0&&Ndtg>19&&round(Ageg*100)/100>0.49&&round(Ageg*100)/100<1.49&&log10(SGred/8.5)>2.9*log10(Sizeg)/3.5-10/3.5&&Zeg>0&&Zeg<0.43&&Sizmg>1111&&Sizeg>11111&&Xcg>=-420&&Xcg<=-60&&Ycg>=-550&&Ycg<=-20){
  
       //if (Sizeg>=pow(10,7.)&&Sizeg<pow(10,7.33)){
       // hdistn1->Fill(Sizmg/Sizeg);
       // myfile1 << Sizmg << endl;
       // myfile2 << Sizeg << endl; 
       // }
       //if (Sizeg>=pow(10,7.33)&&Sizeg<pow(10,7.66)){
       // hdistn2->Fill(Sizmg/Sizeg); 
       // }
       //if (Sizeg>=pow(10,7.66)&&Sizeg<=pow(10,7.99)){
       // hdistn3->Fill(Sizmg/Sizeg);
       // }
        
      //}

     }  

   }

   Gtc.getresult();//in secondi!!!
   cout << Gtc.getresult() << endl;
   chain.Reset();
   TFile *f = new TFile("sel8-03_dqcuts_1.root","recreate");
   qt->Print();
   //qt->SetAutoSave(100000000);
   qt->AutoSave();
   qt->Write();
   f->Close();
   //myfile1.close();
   //myfile2.close();
   //TFile *ana = new TFile("ana_histos_third.root","recreate");
   //ana->cd();
   //hangle->Write();
   //hscatt->Write();
   //hxy->Write();   
   //hdistn1->Write();
   //hdistn2->Write();
   //hdistn3->Write();   
   //ana->Close();

   return 0;  
  
}//end main