{ 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; i19&&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(); }