#define tempostrip_cxx
#include "tempostrip.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
#include<string.h>


void tempostrip::Loop(Long64_t HowMany)
{
  
  if (fChain == 0) return;
  TFile *p = new TFile("tempostrip.root","RECREATE");
  
  
  //.............. TEMPOSTRIP Display ....................
  int mynentries=100;
  int counter_good=0;
  TH1D *timerise[mynentries];
  //
  std::string NomeEvt[100];
  for(int ich=0; ich<mynentries; ich++)
    {
      char buffer [100];
      int n = sprintf(buffer,"%d",ich);
      NomeEvt[ich]=buffer;
      std::string hname ="TimeRise"+NomeEvt[ich];
      std::string htitle = "histo time per strip "+NomeEvt[ich];
      timerise[ich]= new TH1D(hname.c_str(), htitle.c_str(),27,0.,27.);
      
    }
  //
  
  Int_t ContaLeStrip =0;
  //................................................  

  Long64_t nentries;
  //  Long64_t nentries = fChain->GetEntriesFast();
  if(HowMany==0){
    nentries = fChain->GetEntriesFast();
  }
  else {
    nentries= HowMany;
  }
  Long64_t nbytes = 0, nb = 0;
  for (Long64_t jentry=0; jentry<nentries;jentry++) {
    Long64_t ientry = LoadTree(jentry);
    if (ientry < 0) break;
    nb = fChain->GetEntry(jentry);   nbytes += nb;
    
    vector<vector<vector<short> > > charge;
    vector<Int_t> chamber;
    vector<short> apvq;
    apvq.clear();
    charge.clear();
    chamber.clear();
    //
    for(int i=0;i<5;i++) 
      {
     	vector<vector<short> > maxq;
     	maxq.clear();
     	charge.push_back(maxq);
      }
    //
    Int_t nstrip =  mm_strip->size();  
    for(int i=0;i<nstrip;++i)
      {
	apvq = apv_q->at(i);
	for(Int_t j=0; j<27; j++) {
	  Float_t q = apvq.at(j);
	  if(ContaLeStrip<100)  timerise[ContaLeStrip]->Fill(j,q);
	}
      
	hist[k]->Write();
	ContaLeStrip++;
      }
    
    
  }
  
  
  p->Write();
  p->Close();
  
}