#include "TFile.h"
#include "TTree.h"
#include "TString.h"
#include <assert.h>
#include <iostream>
#include <sstream>
#include "TCanvas.h"
#include "TStyle.h"
#include "TH1D.h"

void clusterres(TString filename){
  
  TFile* file=new TFile(filename,"READ");
  TTree* tr=(TTree*)file->Get("cbmsim");
  assert(tr!=0);

  gStyle->SetOptFit();

  TCanvas* c=new TCanvas("Clusters","Clusters",10,10,500,900); 
  c->Divide(3,4);

  for(int i=1; i<11; ++i){
    c->cd(i);
    std::ostringstream histo;
    histo  <<"hRX"<<i;
    std::ostringstream select;
    select << "TrackFitStat._TpcResX>>"<<histo.str()<<"(100,-1,1)";
    std::ostringstream cut;
    cut<<"TrackFitStat._TpcClusterSize=="<<i;
    tr->Draw(select.str().c_str(),cut.str().c_str());
    TH1F* hX = (TH1F*)gDirectory->Get(histo.str().c_str());
    hX->Fit("gaus");
  }


}