#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
#include <stdlib.h>
#include <TROOT.h>
#include <TArrayI.h>
#include <TBranch.h>
#include <TCanvas.h>
#include <TClassTable.h>
#include <TClonesArray.h>
#include <TFile.h>
#include <TGeoManager.h>
#include <TH1.h>
#include <TH2.h>
#include <TStopwatch.h>
#include <TMath.h>
#include <TSystem.h>
#include <TRandom.h>
#include <TInterpreter.h>
#include <TObject.h>
#include <TObjArray.h>
#include <TTree.h>
#include <TChain.h>
#include <TCanvas.h>
#include <TView.h>
#include <TNtuple.h>
#include <TPolyMarker3D.h>
#include <TFile.h>
#endif

void calib(char* infile="rivelatore.root",bool output=1, float slope=2.3,float interc=10)
{

  TNtuple *pjmca_calib;
  
  if(output){
    
    char str_in[50]=infile;
    char str_out[50];
    int nb=strlen(str_in)-5;
    strncpy(str_out,str_in,nb);
    str_out[nb]='\0';
    strcat(str_out,"_calib.root");
    //    strcat(str_out,"root");
    TFile *fout=new TFile(str_out,"RECREATE");
    cout<<"output file  "<< str_out<<endl;
    pjmca_calib=new TNtuple("pjmca_calib","","ch0:ch1:ene:ch3"); 
  }
  
  //create a TNtuple with the channel 2 replaced by the calibrated value
  
  TH1F *h0=new TH1F("h0","h0 ",2048,0,2048);
  TH1F *h1 =new TH1F("h1","h1 ",2048,0,2048);
  TH1F *h2 =new TH1F("h2","h2 ",2048,interc,2048*slope+interc);
  TH1F *h3 =new TH1F("h3","h3 ",2048,0,2048);
  
  
  TFile *fin=new TFile(infile);
  
  TTree *pjmca=(TTree*)fin->Get("pjmca");
  float ch1, ch2,ch3,ch0,temp,ene;
  pjmca->SetBranchAddress("ch0",&ch0);
  pjmca->SetBranchAddress("ch1",&ch1);
  pjmca->SetBranchAddress("ch2",&ch2);
  pjmca->SetBranchAddress("ch3",&ch3);
  
  int nentries=pjmca->GetEntries();
  for(int i=0;i<nentries;i++){
    if(i%50000 == 0) printf("event %d \n",i);
    pjmca->GetEntry(i);
    
    //instruction for a cut  if(ch0>10 && ch1<1000){}
    
    h0->Fill(ch0);
    h1->Fill(ch1);
    h3->Fill(ch3);
    
    //the new value for the third channel, a calibrated value!
    
    ene=ch2*slope+interc;
    
    //trick for the binning
    Float_t random = gRandom->Rndm(1);
    ene+=random-0.5;
    
    h2->Fill(ene);
    //fill the new ntuple 
    if(output) pjmca_calib->Fill(ch0,ch1,ene,ch3);
  }
  
  
  if(!output)h2->Draw();
  
  if(output){
    fout->cd();
    fout->Write();
    fout->Close();
  }
  
  
  

}