void test(){

   TFile *fout = TFile::Open("sparse.root");

   THnSparse *hsparse = fout->Get("RSN_config1sigmadeltatof_LikePP");   
   
   THnSparse *h_x_to_x1 = (THnSparse*) hsparse->Clone("h1"); //!from ptbin[x] to ptbin[x+1]
   THnSparse *h_0_to_x =  (THnSparse*) hsparse->Clone("h2"); //!from ptbin[0] to ptbin[x]
   THnSparse *h_0_to_x1 = (THnSparse*) hsparse->Clone("h3"); //!from ptbin[0] to ptbin[x+1]
   
   Double_t ptbin[] = {0.1,0.2,0.4,0.6,0.8,1.0,1.2,1.4};
   Int_t nptbin  = (sizeof(ptbin) / sizeof(ptbin[0]))-1; 

   Double_t pmin,  Double_t pmax,  Double_t diff;

   for (Int_t x = 0; x < nptbin; x++) {
    
    pmin= ptbin[x];
    pmax= ptbin[x+1];
 
    h_x_to_x1->GetAxis(1)->SetRangeUser(pmin,pmax); 
    h_0_to_x->GetAxis(1)->SetRangeUser(ptbin[0],pmin);
    h_0_to_x1->GetAxis(1)->SetRangeUser(ptbin[0],pmax);

    TH1D *out_0_to_x = h_0_to_x->Projection(0);
    out_0_to_x->SetName(Form("0tox_%f",x));
    
   
    TH1D *out_0_to_x1 = h_0_to_x1->Projection(0);
    out_0_to_x1->SetName(Form("0tox1_%f",x));
    
    TH1D *out_x_to_x1 = h_x_to_x1->Projection(0);
    out_x_to_x1->SetName(Form("xtox1_%f",x));
    
    diff=(TMath::Abs(out_0_to_x->GetEntries()-out_0_to_x1->GetEntries()))-out_x_to_x1->GetEntries();
    
    cout<<"//////start-->"<<x<<endl;
    cout<<"from " <<ptbin[0]<<" to "<<pmin<<" entries"<<out_0_to_x->GetEntries()<<endl;
    cout<<"from " <<pmin<<" to "<<pmax<<" entries"<<out_x_to_x1->GetEntries()<<endl;   
    cout<<"from " <<ptbin[0]<<" to "<<pmax<<" entries"<<out_0_to_x1->GetEntries()<<endl;  
    cout<<"difference "<<diff<<endl;
    cout<<"//////end-->"<<x<<endl;
    
   }
   
 }