MrpcCount.cxx

Go to the documentation of this file.

#include "MrpcRec/MrpcTrack.h"
#include "MrpcRec/MrpcCount.h"
#include <iostream>
#include <iomanip>
 
MrpcCount::MrpcCount() {
  extTrkNum = 0;
  barrelTrk = 0;
  endcapTrk = 0;
  inner1  = 0;
  outer1  = 0;
  double1 = 0;
  east1   = 0;
  west1   = 0;
  nohit1  = 0;
  unchanged1 = 0;
  changed1   = 0;
  inner2  = 0;
  outer2  = 0;
  double2 = 0;
  east2   = 0;
  west2   = 0;
  nohit2  = 0;
  unchanged2 = 0;
  changed2   = 0;
  inner3  = 0;
  outer3  = 0;
  double3 = 0;
  east3   = 0;
  west3   = 0;
  nohit3  = 0;
  unchanged3 = 0;
  changed3   = 0;
  inner4  = 0;
  outer4  = 0;
  double4 = 0;
  east4   = 0;
  west4   = 0;
  nohit4  = 0;
  unchanged4 = 0;
  changed4   = 0;
  mrpc_east1_1=0;
  mrpc_east1_2=0; 
  mrpc_east1_3=0;
  mrpc_east1_4=0;
  mrpc_east2_1=0; 
  mrpc_east2_2=0; 
  mrpc_east2_3=0; 
  mrpc_east2_4=0;
  mrpc_west1_1=0;
  mrpc_west1_2=0;
  mrpc_west1_3=0;
  mrpc_west1_4=0;
  mrpc_west2_1=0;
  mrpc_west2_2=0;
  mrpc_west2_3=0;
  mrpc_west2_4=0;
  mrpc_extrapolation_good_3=0;
  mrpc_extrapolation_false_3=0;
  return;
}
 
 
void MrpcCount::init() {
  extTrkNum = 0;
  barrelTrk = 0;
  endcapTrk = 0;
  inner1  = 0;
  outer1  = 0;
  double1 = 0;
  east1   = 0;
  west1   = 0;
  nohit1  = 0;
  unchanged1 = 0;
  changed1   = 0;
  inner2  = 0;
  outer2  = 0;
  double2 = 0;
  east2   = 0;
  west2   = 0;
  nohit2  = 0;
  unchanged2 = 0;
  changed2   = 0;
  inner3  = 0;
  outer3  = 0;
  double3 = 0;
  east3   = 0;
  west3   = 0;
  nohit3  = 0;
  unchanged3 = 0;
  changed3   = 0;
  inner4  = 0;
  outer4  = 0;
  double4 = 0;
  east4   = 0;
  west4   = 0;
  nohit4  = 0;
  unchanged4 = 0;
  changed4   = 0;
  mrpc_east1_1=0;
  mrpc_east1_2=0;
  mrpc_east1_3=0;
  mrpc_east1_4=0;
  mrpc_east2_1=0;
  mrpc_east2_2=0;
  mrpc_east2_3=0;
  mrpc_east2_4=0;
  mrpc_west1_1=0;
  mrpc_west1_2=0;
  mrpc_west1_3=0;
  mrpc_west1_4=0;
  mrpc_west2_1=0;
  mrpc_west2_2=0;
  mrpc_west2_3=0;
  mrpc_west2_4=0;
  mrpc_extrapolation_good_3=0;
  mrpc_extrapolation_false_3=0;
  return;
}
 
 
void MrpcCount::setExtTrackNum( unsigned int ntrk ) {
  extTrkNum = extTrkNum + ntrk;
  return;
}
 
 
void MrpcCount::setTrack1( MrpcTrack*& tof ) {
  if( tof->hitCase() == InnerLayer )       inner1++;
  else if( tof->hitCase() == OuterLayer )  outer1++;
  else if( tof->hitCase() == DoubleLayer ) double1++;
  else if( tof->hitCase() == EastEndcap )  east1++;
  else if( tof->hitCase() == WestEndcap )  west1++;
  else if( tof->hitCase() == eastendcap_1 ) mrpc_east1_1++;
  else if( tof->hitCase() == eastendcap_2 ) mrpc_east2_1++;
  else if( tof->hitCase() == westendcap_1 ) mrpc_west1_1++;
  else if( tof->hitCase() == westendcap_2 ) mrpc_west2_1++;
  else                                     nohit1++;
  return;
}
 
 
void MrpcCount::setTrack1Col( std::vector<MrpcTrack*>*& tofTrackVec ) {
  TofTrackVec::iterator iter_tof = tofTrackVec->begin();
  for( ; iter_tof!=tofTrackVec->end(); iter_tof++ ) {
    setTrack1( *iter_tof );
  }
  return;
}
 
 
void MrpcCount::setTrack2( MrpcTrack*& tof ) {
  if( tof->hitCase() == InnerLayer )       inner2++;
  else if( tof->hitCase() == OuterLayer )  outer2++;
  else if( tof->hitCase() == DoubleLayer ) double2++;
  else if( tof->hitCase() == EastEndcap )  east2++;
  else if( tof->hitCase() == WestEndcap )  west2++;
  else if( tof->hitCase() == eastendcap_1 ) mrpc_east1_2++;
  else if( tof->hitCase() == eastendcap_2 ) mrpc_east2_2++;
  else if( tof->hitCase() == westendcap_1 ) mrpc_west1_2++;
  else if( tof->hitCase() == westendcap_2 ) mrpc_west2_2++;
  else                                     nohit2++;
  //  if( tof->changed() ) changed2++;
  //  else                 unchanged2++;
  return;
}
 
 
void MrpcCount::setTrack3( MrpcTrack*& tof ) {
  if( tof->hitCase() == InnerLayer )       inner3++;
  else if( tof->hitCase() == OuterLayer )  outer3++;
  else if( tof->hitCase() == DoubleLayer ) double3++;
  else if( tof->hitCase() == EastEndcap )  east3++;
  else if( tof->hitCase() == WestEndcap )  west3++;
  else if( tof->hitCase() == eastendcap_1 ) mrpc_east1_3++;
  else if( tof->hitCase() == eastendcap_2 ) mrpc_east2_3++;
  else if( tof->hitCase() == westendcap_1 ) mrpc_west1_3++;
  else if( tof->hitCase() == westendcap_2 ) mrpc_west2_3++;
  else                                     nohit3++;
  //  if( tof->changed() ) changed3++;
  //  else                 unchanged3++;
 
 
 
  if(  ( tof->hitCase() == eastendcap_1 )  ||  ( tof->hitCase() == eastendcap_2 )  ||  ( tof->hitCase() == westendcap_1 )  ||   ( tof->hitCase() == westendcap_2 ))
    {
      if(tof->is_mrpc_extrapolation_good()==true) mrpc_extrapolation_good_3++;
      else mrpc_extrapolation_false_3++;
    }
 
  return;
}
 
 
void MrpcCount::setTrack4( MrpcTrack*& tof ) {
  if( tof->hitCase() == InnerLayer )       inner4++;
  else if( tof->hitCase() == OuterLayer )  outer4++;
  else if( tof->hitCase() == DoubleLayer ) double4++;
  else if( tof->hitCase() == EastEndcap )  east4++;
  else if( tof->hitCase() == WestEndcap )  west4++;
  else if( tof->hitCase() == eastendcap_1 ) mrpc_east1_4++;
  else if( tof->hitCase() == eastendcap_2 ) mrpc_east2_4++;
  else if( tof->hitCase() == westendcap_1 ) mrpc_west1_4++;
  else if( tof->hitCase() == westendcap_2 ) mrpc_west2_4++;
  else                                     nohit4++;
  //  if( tof->changed() ) changed4++;
  //  else                 unchanged4++;
  return;
}
 
/*
void TofCount::final() {
   std::cout<<"############################# TofRec Summary #####################################" << std::endl;
   std::cout<<" ===> The Extend tracks in total:  " << extTrkNum <<std::endl;
   std::cout<<"     ===> The valid Extend tracks in total:  "<< (inner4+outer4+double4+east4+west4) <<std::endl;
   std::cout<<"          ===> The TofRec barrel tracks in total:  " << (inner4+outer4+double4) <<std::endl;
   std::cout<<"              ===> The Inner Layer hit only in total:   " << inner1 << "   " << inner2 << "   " << inner3 << "   " << inner4 << std::endl;
   std::cout<<"              ===> The Outer Layer hit only in total:   " << outer1 << "   " << outer2 << "   " << outer3 << "   " << outer4 << std::endl;
   std::cout<<"              ===> The Double Layer hit in total:       " << double1 << "   " << double2 << "   " << double3 << "   " << double4 << std::endl;
   std::cout<<"          ===> The TofRec endcap tracks in total:  " << (east4+west4) <<std::endl;
   std::cout<<"              ===> The TofRec East endcap tracks in total: " << east1 << "   " << east2 << "   " << east3 << "   " << east4 << std::endl;
   std::cout<<"              ===> The TofRec West endcap tracks in total: " << west1 << "   " << west2 << "   " << west3 << "   " << west4 << std::endl;
   std::cout<<"     ===> The invalid Extend tracks in total:  "<< nohit1 << "   " << nohit2 << "   " << nohit3 << "   " << nohit4 <<std::endl;
   std::cout<<"#############################  End Summary #######################################" <<std::endl;
 
}
*/
 
void MrpcCount::final() {
   std::cout<<"############################# TofRec Summary #####################################" << std::endl;
   std::cout<<" ===> Number of Extrapolated tracks:  " << extTrkNum <<std::endl;
   std::cout<<"   ===> Number of valid Extraploated tracks:  " << (inner1+outer1+double1+east1+west1+mrpc_east1_1+mrpc_east2_1+mrpc_west1_1+mrpc_west2_1) << "  Barrel: " << (inner1+outer1+double1) << "  Endcap: " << (east1+west1+mrpc_east1_1+mrpc_east2_1+mrpc_west1_1+mrpc_west2_1) << std::endl;
   std::cout<<"     ===> TofRec: number of barrel tracks:  " << (inner3+outer3+double3) << "  Efficiency:  " << setprecision(4) << (inner3+outer3+double3)*100.0/((inner1+outer1+double1)*1.0) << "%" << std::endl;
   std::cout<<"       ===> Inner Layer hit only in total:   " << inner1 << " \t" << inner2 << " \t" << inner3 << " \t" << inner4 << std::endl;
   std::cout<<"       ===> Outer Layer hit only in total:   " << outer1 << "\t" << outer2 << " \t" << outer3 << " \t" << outer4 << std::endl;
   std::cout<<"       ===> Double Layer hit in total:       " << double1 << " \t" << double2 << " \t" << double3 << " \t" << double4 << std::endl;
   if((east1+west1+mrpc_east1_1+mrpc_east2_1+mrpc_west1_1+mrpc_west2_1)!=0) std::cout<<"     ===> TofRec: number of endcap tracks:  " << (east3+west3+mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3) << "  Efficiency:  " << setprecision(4) << (east3+west3+mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*100.0/((east1+west1+mrpc_east1_1+mrpc_east2_1+mrpc_west1_1+mrpc_west2_1)*1.0) << "%" << std::endl;
   std::cout<<"       ===> East endcap tracks in total: " << east1 << " \t" << east2 << " \t" << east3 << " \t" << east4 << std::endl;
   std::cout<<"       ===> West endcap tracks in total: " << west1 << " \t" << west2 << " \t" << west3 << " \t" << west4 << std::endl;
   std::cout<<"       ===> East 1 mrpc tracks in total: " << mrpc_east1_1 << " \t" << mrpc_east1_2 << " \t" << mrpc_east1_3 << " \t" << mrpc_east1_4 << std::endl;
   std::cout<<"       ===> West 1 mrpc tracks in total: " << mrpc_west1_1 << " \t" << mrpc_west1_2 << " \t" << mrpc_west1_3 << " \t" << mrpc_west1_4 << std::endl;
   std::cout<<"       ===> East 2 mrpc tracks in total: " << mrpc_east2_1 << " \t" << mrpc_east2_2 << " \t" << mrpc_east2_3 << " \t" << mrpc_east2_4 << std::endl;
   std::cout<<"       ===> West 2 mrpc tracks in total: " << mrpc_west2_1 << " \t" << mrpc_west2_2 << " \t" << mrpc_west2_3 << " \t" << mrpc_west2_4 << std::endl;
 
   std::cout<<"     ===> The invalid Extend tracks in total:  "<< nohit1 << " \t" << nohit2 << " \t" << nohit3 << " \t" << nohit4 <<std::endl;
 
   if((mrpc_extrapolation_good_3+mrpc_extrapolation_false_3)!=0)std::cout<<"  ==> Correct MRPC extrapolation:              " << (mrpc_extrapolation_good_3*100.0)/((mrpc_extrapolation_good_3+mrpc_extrapolation_false_3)*1.0)<< "%"  << std::endl;
 
   std::cout<<"#############################  End Summary #######################################" <<std::endl;
 
}