for (j=0;j<muplus.GetLength();++j) { TLorentzVector lv=muplus[j]->P4(); RhoError lvCov=muplus[j]->P4Cov(); hpullpx->Fill( ( lv.Px() - lvmc.Px() ) / sqrt( lvCov[0][0] ) ); hpullpy->Fill( ( lv.Py() - lvmc.Py() ) / sqrt( lvCov[1][1] ) ); hpullpz->Fill( ( lv.Pz() - lvmc.Pz() ) / sqrt( lvCov[2][2] ) ); hpulle->Fill( ( lv.E() - lvmc.E() ) / sqrt( lvCov[3][3] ) ); }