/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/Iro/read_ntuples.C File Reference

#include <iostream>
#include <cmath>
#include <string>
#include <vector>
#include "TChain.h"
#include "TProfile.h"
#include "TFile.h"
#include "TTree.h"
#include "TString.h"

Include dependency graph for read_ntuples.C:

Go to the source code of this file.

Functions
void	read_ntuples ()

---

Function Documentation

void read_ntuples

(

)

Definition at line 12 of file read_ntuples.C.

                   {
  gStyle->SetOptStat(0);

  TFile *f1=new TFile("DSTPionsAllEM.root","READ");
  TTree *t1 = (TTree*)(f1->Get("Pions"));

  int num_entries=t1->GetEntries();
  cout << "entries " << num_entries << endl;
  double energy; double Nhtot; double Nhtot1; double Nhtot2; double Nhtot3; float ShCont1; double NbHitsRear; double NbHitsBoard; double Zstart;
  t1 -> SetBranchAddress("Ebeam",&energy);
  t1 -> SetBranchAddress("Nhtot",&Nhtot);
  t1 -> SetBranchAddress("Nhtot1",&Nhtot1);
  t1 -> SetBranchAddress("Nhtot2",&Nhtot2);
  t1 -> SetBranchAddress("Nhtot3",&Nhtot3);
  t1 -> SetBranchAddress("ShCont1",&ShCont1);
  t1 -> SetBranchAddress("NbHitsRear",&NbHitsRear);
  t1 -> SetBranchAddress("NbHitsBoard",&NbHitsBoard);
  t1 -> SetBranchAddress("Zstart",&Zstart);

  vector<double> energies;
  double energy_step=0;

  for ( int ient = 0; ient < num_entries; ient++ ) {
    t1->GetEntry(ient);

    if (energy!=energy_step) energies.push_back(energy);
    energy_step=energy;

  }

  int size_energies=energies.size();
  cout << "Nb of different energies " << size_energies << endl;

}

---