Limit the bins - to only the first 77 bins - so a range up to ~ 1000 leave the previous selection commented for the time being
the error estimates are always positive ...
Limit the bins - to only the first 77 bins - so a range up to ~ 1000 leave the previous selection commented for the time being
vertex analyses if required ...
  120   std::cout << 
"\nAnalysis_Tier0:: chain specification: " << cname << 
"\t" << cname.
raw() << std::endl;
 
  121   std::cout << 
"\tchain: " << cname.
head()    << std::endl;
 
  122   std::cout << 
"\tkey:   " << cname.
tail()    << std::endl;
 
  123   std::cout << 
"\troi:   " << cname.
roi()     << std::endl;
 
  124   std::cout << 
"\tvtx:   " << cname.
vtx()     << std::endl;
 
  125   std::cout << 
"\tte:    " << cname.
element() << std::endl;
 
  131   m_h_chain->GetXaxis()->SetBinLabel(1, 
"Nrois" );
 
  132   m_h_chain->GetXaxis()->SetBinLabel(2, 
"Nevents" );
 
  133   m_h_chain->GetXaxis()->SetBinLabel(3, 
"N ref tracks" );
 
  134   m_h_chain->GetXaxis()->SetBinLabel(4, 
"N matched tracks" );
 
  135   m_h_chain->GetXaxis()->SetBinLabel(5, 
"N total test tracks" );
 
  146   double vnbins[81] = {
 
  148     0.5,    1.5,    2.5,    3.5,    4.5,    5.5,    6.5,    7.5,    8.5,    9.5,
 
  149     10.5,   11.5,   12.5,   13.5,   14.5,   15.5,   16.5,   17.5,   18.5,   19.5,
 
  150     20.5,   21.5,   22.5,   23.5,   24.5,   25.5,   26.5,   27.5,   28.5,   29.5,
 
  151     31.5,   32.5,   33.5,   34.5,   36.5,   37.5,   39.5,
 
  152     40.5,   42.5,   43.5,   45.5,   47.5,   49.5,
 
  153     50.5,   52.5,   54.5,   57.5,   59.5,
 
  154     61.5,   63.5,   66.5,   69.5,
 
  171   double vnbins[101] = { 
 
  173     0.5,   1.5,   2.5,   3.5,   4.5,   5.5,   6.5,   7.5,   8.5,   9.5,   10.5,   11.5,   12.5,   13.5,   14.5,   15.5,   17.5,   18.5,   19.5,   21.5,
 
  174     23.5,   24.5,   26.5,   28.5,   30.5,   32.5,   35.5,   37.5,   40.5,   43.5,   46.5,   50.5,   53.5,   57.5,   61.5,   66.5,   71.5,   76.5,   81.5,   87.5,
 
  175     93.5,   100.5,   107.5,   114.5,   123.5,   131.5,   141.5,   150.5,   161.5,   172.5,   185.5,   198.5,   211.5,   226.5,   242.5,   259.5,   277.5,   297.5,   317.5,   340.5,
 
  176     363.5,   389.5,   416.5,   445.5,   476.5,   509.5,
 
  177     544.5,   582.5,   623.5,   666.5,   713.5,   762.5,   815.5,   872.5,   933.5,   998.5,   1067.5,
 
  178     1141.5,   1221.5,   1305.5,   1396.5,   1493.5,   1597.5,
 
  179     1708.5,   1827.5,   1953.5,   2089.5,
 
  180     2234.5,   2389.5,   2555.5,
 
  181     2733.5,   2923.5,   3125.5,
 
  192   m_h_ntrk = 
new TH1F( 
"reftrk_N", 
"Reference tracks", 77, vnbins );
 
  200   std::vector<double> ptbins;
 
  202   for ( 
int i=0 ; 
i<=25 ; 
i++ ) {
 
  204     double y = std::log10(1) + 
i*( std::log10(100)-std::log10(1) )/25;
 
  212   m_h_trkpT  = 
new TH1F(
"reftrk_pT" , 
"Reference track pT",  25,   &ptbins[0]   );
 
  214   m_h_trketa = 
new TH1F(
"reftrk_eta", 
"Reference track Eta", 25,   -2.5,     2.5) ;
 
  215   if (
name().
find(
"LRT")!=std::string::npos || 
name().
find(
"lrt")!=std::string::npos) { 
 
  216     m_h_trkd0  = 
new TH1F(
"reftrk_d0" , 
"Reference track d0", 201, -100.0,  100.0 );
 
  218     m_h_trkd0  = 
new TH1F(
"reftrk_d0" , 
"Reference track d0", 101,   -5.0,    5.0 );
 
  220   m_h_trkz0  = 
new TH1F(
"reftrk_z0" , 
"Reference track z0",  50,   -225.,    225.);
 
  223   m_h_trkdd0  = 
new TH1F(
"reftrk_dd0" , 
"Reference track sigma(d0)", 50,  0,     0.5);
 
  224   m_h_trkdz0  = 
new TH1F(
"reftrk_dz0" , 
"Reference track sigma(z0)", 50,  0,     2.5);
 
  226   m_h_trkd0sig = 
new TH1F(
"reftrk_d0sig" , 
"Reference track d0 significance", 101,   -5.,     5.);
 
  253   if (
name().
find(
"LRT")!=std::string::npos || 
name().
find(
"lrt")!=std::string::npos) { 
 
  279   double d0bins_LRT[124] = { -100.5, 
 
  283                  -69.5, -66.5, -63.5, -61.5, 
 
  284                  -59.5, -57.5, -54.5, -52.5, -50.5, 
 
  285                  -49.5, -47.5, -45.5, -43.5, -42.5, -40.5, 
 
  286                  -39.5, -37.5, -36.5, -34.5, -33.5, -32.5, -31.5, 
 
  287                  -29.5, -28.5, -27.5, -26.5, -25.5, -24.5, -23.5, -22.5, -21.5, -20.5, 
 
  288                  -19.5, -18.5, -17.5, -16.5, -15.5, -14.5, -13.5, -12.5, -11.5, -10.5, 
 
  289                  -9.5, -8.5, -7.5, -6.5, -5.5, -4.5, -3.5, -2.5, -1.5, -0.5,
 
  290                  0.5,  1.5,  2.5,  3.5,  4.5,  5.5,  6.5,  7.5,  8.5,  9.5,  
 
  291                  10.5,  11.5,  12.5,  13.5,  14.5,  15.5,  16.5,  17.5,  18.5,  19.5, 
 
  292                  20.5,  21.5,  22.5,  23.5,  24.5,  25.5,  26.5,  27.5,  28.5,  29.5, 
 
  293                  31.5,  32.5,  33.5,  34.5,  36.5,  37.5,  39.5, 
 
  294                  40.5,  42.5,  43.5,  45.5,  47.5,  49.5, 
 
  295                  50.5,  52.5,  54.5,  57.5,  59.5, 
 
  296                  61.5,  63.5,  66.5,  69.5, 
 
  302   double d0bins[40] = { -5.0,  -4.0,  -3.0,  -2.5,   
 
  303             -2.0,  -1.8,  -1.6,  -1.4,  -1.2,  
 
  304             -1.05, -0.95, -0.85, -0.75, -0.65, -0.55, -0.45, -0.35, -0.25, -0.15, -0.05, 
 
  305              0.05,  0.15,  0.25,  0.35,  0.45,  0.55,  0.65,  0.75,  0.85,  0.95,  1.05,  
 
  306              1.2,   1.4,   1.6,   1.8,   2.0,  
 
  307              2.5,   3.0,   4.0,   5.0 };
 
  317   if (
name().
find(
"LRT")!=std::string::npos || 
name().
find(
"lrt")!=std::string::npos) { 
 
  326   m_h_lbeff = 
new TProfile( 
"Eff_lb", 
"efficinecy vs lumiblock", 301, -0.5, 3009.5 );
 
  339   m_h_npixvseta     = 
new TProfile(
"npix_vs_eta",      
"offline npix vs eta;offline #eta;# mean number of offline pixel hits", 30,   -2.5,  2.5);
 
  340   m_h_npixvseta_rec = 
new TProfile(
"npix_vs_eta_rec",  
"trigger npix vs eta;offline #eta;# mean number of trigger pixel hits", 30,   -2.5,  2.5);
 
  342   m_h_nsctvseta     = 
new TProfile(
"nsct_vs_eta",      
"offline nsct vs eta;offline #eta;# mean number of offline sct hits", 30,   -2.5,  2.5);
 
  343   m_h_nsctvseta_rec = 
new TProfile(
"nsct_vs_eta_rec",  
"trigger nsct vs eta;offline #eta;# mean number of trigger sct hits", 30,   -2.5,  2.5);
 
  345   m_h_ntrtvseta     = 
new TProfile(
"ntrt_vs_eta",      
"offline ntrt vs eta;offline #eta;# mean number of offline trt hits", 30,   -2.5, 2.5 );
 
  346   m_h_ntrtvseta_rec = 
new TProfile(
"ntrt_vs_eta_rec",  
"trigger ntrt vs eta;offline #eta;# mean number of trigger trt hits", 30,   -2.5, 2.5 );
 
  357   if (
name().
find(
"LRT")!=std::string::npos || 
name().
find(
"lrt")!=std::string::npos) {
 
  358     m_h_npixvsd0     = 
new TProfile(
"npix_vs_d0",      
"offline npix vs d0;offline #d0;# mean number of offline pixel hits", 123,  d0bins_LRT );
 
  359     m_h_npixvsd0_rec = 
new TProfile(
"npix_vs_d0_rec",  
"trigger npix vs d0;offline #d0;# mean number of trigger pixel hits", 123,  d0bins_LRT );
 
  361     m_h_nsctvsd0     = 
new TProfile(
"nsct_vs_d0",      
"offline nsct vs d0;offline #d0;# mean number of offline sct hits", 123,   d0bins_LRT );
 
  362     m_h_nsctvsd0_rec = 
new TProfile(
"nsct_vs_d0_rec",  
"trigger nsct vs d0;offline #d0;# mean number of trigger sct hits", 123,   d0bins_LRT );
 
  364     m_h_npixvsd0     = 
new TProfile(
"npix_vs_d0",      
"offline npix vs d0;offline #d0;# mean number of offline pixel hits", 39,  d0bins );
 
  365     m_h_npixvsd0_rec = 
new TProfile(
"npix_vs_d0_rec",  
"trigger npix vs d0;offline #d0;# mean number of trigger pixel hits", 39,  d0bins );
 
  367     m_h_nsctvsd0     = 
new TProfile(
"nsct_vs_d0",      
"offline nsct vs d0;offline #d0;# mean number of offline sct hits", 39,   d0bins );
 
  368     m_h_nsctvsd0_rec = 
new TProfile(
"nsct_vs_d0_rec",  
"trigger nsct vs d0;offline #d0;# mean number of trigger sct hits", 39,   d0bins );
 
  371   m_h_npixvspT     = 
new TProfile(
"npix_vs_pT",      
"offline npix vs pT;offline #p_{T};# mean number of offline pixel hits", 25, &ptbins[0] );
 
  372   m_h_npixvspT_rec = 
new TProfile(
"npix_vs_pT_rec",  
"trigger npix vs pT;offline #p_{T};# mean number of trigger pixel hits", 25, &ptbins[0] );
 
  374   m_h_nsctvspT     = 
new TProfile(
"nsct_vs_pT",      
"offline nsct vs pT;offline #p_{T} [GeV];# mean number of offline pixel hits", 25, &ptbins[0] );
 
  375   m_h_nsctvspT_rec = 
new TProfile(
"nsct_vs_pT_rec",  
"trigger nsct vs pT;offline #p_{T} [GeV];# mean number of trigger pixel hits", 25, &ptbins[0] );
 
  385   m_h_layer      = 
new TH1F(
"layer",      
"hit layers", 32,   -0.5,  31.5  );
 
  442   m_h_pTres  = 
new TProfile(
"Res_pT",  
"pT residual;#eta_{ref};#Delta pT [Gev]",  25,   -2.5,  2.5   );
 
  443   m_h_ipTres = 
new TProfile(
"Res_ipT", 
"ipT residual;#eta_{ref};#Delta pT^{-1} [Gev^{-1}]", 25, -2.5, 2.5);
 
  444   m_h_etares = 
new TProfile(
"Res_eta", 
"Eta residual;#eta_{ref};#Delta #eta",     25,   -2.5,  2.5);
 
  445   m_h_phires = 
new TProfile(
"Res_phi", 
"Phi residual;#eta_{ref};#Delta #phi",     25,   -2.5,  2.5);
 
  446   m_h_d0res  = 
new TProfile(
"Res_d0",  
"d0 residual;#eta_{ref};#Delta d0 [mm]",   25,   -2.5,  2.5);
 
  447   m_h_z0res  = 
new TProfile(
"Res_z0",  
"z0 residual;#eta_{ref};#Delta z0 [mm]",   25,   -2.5,  2.5);
 
  526   if ( 
name().
find(
"vtx")!=std::string::npos || 
name().
find(
"Vtx")!=std::string::npos ||
 
  527        name().
find(
"vx")!=std::string::npos || 
name().
find(
"Vx")!=std::string::npos ) {