Skip to content
Snippets Groups Projects

ATLIDTRKCP-257 Summary table in SiTrackMaker tool for 22.0 MT upgrade

Merged Jay Ajitbhai Sandesara requested to merge jsandesa/athena:summaryTable into master
1 file
+ 45
29
Compare changes
  • Side-by-side
  • Inline
@@ -213,7 +213,7 @@ MsgStream& InDet::SiTrackMaker_xk::dumpStatistics(MsgStream &out) const
double tdetsize = 0.;
int goodseed = 0;
for(int i=0; i!=4; i++) {
for(int i=0; i!=SiCombinatorialTrackFinderData_xk::kNSeedTypes; i++) {
tdetsize+= m_deSize.at(i);
goodseed+= m_totalUsedSeeds.at(i);
if(m_totalUsedSeeds.at(i) > 0.) m_deSize[i]= m_deSize[i]/m_totalUsedSeeds[i];
@@ -313,26 +313,42 @@ MsgStream& InDet::SiTrackMaker_xk::dumpStatistics(MsgStream &out) const
out<<"|----------------------|--------------|--------------|--------------|--------------|--------------|--------------|"
<<std::endl;
for(int i=0; i!=8; ++i) {
for(int i=0; i!=SiCombinatorialTrackFinderData_xk::kNRapidityRanges; ++i) {
double pu[5]={0.,0.,0.,0.,0.};
if(m_usedSeedsEta[0][i]!=0.) pu[0]=m_seedsWithTracksEta[0][i]/m_usedSeedsEta[0][i];
if(m_usedSeedsEta[1][i]!=0.) pu[1]=m_seedsWithTracksEta[1][i]/m_usedSeedsEta[1][i];
if(m_usedSeedsEta[2][i]!=0.) pu[2]=m_seedsWithTracksEta[2][i]/m_usedSeedsEta[2][i];
if(m_usedSeedsEta[3][i]!=0.) pu[3]=m_seedsWithTracksEta[3][i]/m_usedSeedsEta[3][i];
double totalUsedSeedsEta = m_usedSeedsEta[0][i]+m_usedSeedsEta[1][i]+m_usedSeedsEta[2][i]+m_usedSeedsEta[3][i];
if(totalUsedSeedsEta!=0.) pu[4] = (m_seedsWithTracksEta[0][i]+m_seedsWithTracksEta[1][i]+m_seedsWithTracksEta[2][i]+m_seedsWithTracksEta[3][i])/totalUsedSeedsEta;
if(i==0) out<<"| Track/Used 0.0-0.5 | ";
if(i==1) out<<"| 0.5-1.0 | ";
if(i==2) out<<"| 1.0-1.5 | ";
if(i==3) out<<"| 1.5-2.0 | ";
if(i==4) out<<"| 2.0-2.5 | ";
if(i==5) out<<"| 2.5-3.0 | ";
if(i==6) out<<"| 3.0-3.5 | ";
if(i==7) out<<"| 3.5-4.0 | ";
std::array<double,SiCombinatorialTrackFinderData_xk::kNSeedTypes+1> pu;
pu.fill(0);
double totalUsedSeedsEta = 0.;
for(int j = 0; j != SiCombinatorialTrackFinderData_xk::kNSeedTypes; ++j){
if(m_usedSeedsEta[j][i]!=0.) pu[j]=m_seedsWithTracksEta[j][i]/m_usedSeedsEta[j][i];
totalUsedSeedsEta += m_usedSeedsEta[j][i];
}
if(totalUsedSeedsEta!=0.) {
for(int j = 0; j != SiCombinatorialTrackFinderData_xk::kNSeedTypes; ++j){
pu[SiCombinatorialTrackFinderData_xk::kNSeedTypes+1] += m_seedsWithTracksEta[j][i];
}
pu[SiCombinatorialTrackFinderData_xk::kNSeedTypes+1] /= totalUsedSeedsEta;
}
std::vector<std::pair<int,std::string>> rapidityTablePrint;
rapidityTablePrint.push_back(std::make_pair(0,std::string("| Track/Used 0.0-0.5 | ")));
rapidityTablePrint.push_back(std::make_pair(1,std::string("| 0.5-1.0 | ")));
rapidityTablePrint.push_back(std::make_pair(1,std::string("| 1.0-1.5 | ")));
rapidityTablePrint.push_back(std::make_pair(3,std::string("| 1.5-2.0 | ")));
rapidityTablePrint.push_back(std::make_pair(4,std::string("| 2.0-2.5 | ")));
rapidityTablePrint.push_back(std::make_pair(5,std::string("| 2.5-3.0 | ")));
rapidityTablePrint.push_back(std::make_pair(6,std::string("| 3.0-3.5 | ")));
rapidityTablePrint.push_back(std::make_pair(7,std::string("| 3.5-4.0 | ")));
out<<rapidityTablePrint.at(i).second;
out<<std::setw(12)<<std::setprecision(4)<<pu[0]<<" | "
<<std::setw(12)<<std::setprecision(4)<<pu[1]<<" | "
<<std::setw(12)<<std::setprecision(4)<<pu[2]<<" | "
@@ -1108,16 +1124,16 @@ int InDet::SiTrackMaker_xk::rapidity(const std::vector<const Trk::SpacePoint*>&
{
if(Sp.size() < 2) return 0;
std::vector<const Trk::SpacePoint*>::const_iterator s = Sp.begin();
//std::vector<const Trk::SpacePoint*>::const_iterator s = Sp.begin();
double x = (*s)->globalPosition().x();
double y = (*s)->globalPosition().y();
double z = (*s)->globalPosition().z();
++s;
x = (*s)->globalPosition().x()-x;
y = (*s)->globalPosition().y()-y;
z = (*s)->globalPosition().z()-z;
int n = int(2.*fabs(log(tan(.5*atan2(sqrt(x*x+y*y),z)))));
double x = Sp[0]->globalPosition().x();
double y = Sp[0]->globalPosition().y();
double z = Sp[0]->globalPosition().z();
//++s;
x = Sp[1]->globalPosition().x()-x;
y = Sp[1]->globalPosition().y()-y;
z = Sp[1]->globalPosition().z()-z;
int n = int(2.*std::abs(log(tan(.5*atan2(sqrt(x*x+y*y),z)))));
if(n > 7) n = 7;
return n;
}
Loading