MrpcDBSCount.cxx

Go to the documentation of this file.

#include "MrpcRecDBS/MrpcDBSTrack.h"
#include "MrpcRecDBS/MrpcDBSCount.h"
#include <iostream>
#include <iomanip>
 
MrpcDBSCount::MrpcDBSCount() {
  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_trans_with_ext_true_2=0;
  mrpc_trans_with_ext_false_2=0;
  mrpc_trans_with_dbs_true_2=0;
  mrpc_trans_with_dbs_false_2=0;
  mrpc_trans_with_ext_true_3=0;
  mrpc_trans_with_ext_false_3=0;
  mrpc_trans_with_dbs_true_3=0;
  mrpc_trans_with_dbs_false_3=0;
  mrpc_trans_with_ext_true_4=0;
  mrpc_trans_with_ext_false_4=0;
  mrpc_trans_with_dbs_true_4=0;
  mrpc_trans_with_dbs_false_4=0;
 
  return;
}
 
 
void MrpcDBSCount::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_trans_with_ext_true_2=0;
  mrpc_trans_with_ext_false_2=0;
  mrpc_trans_with_dbs_true_2=0;
  mrpc_trans_with_dbs_false_2=0;
 
  mrpc_trans_with_ext_true_3=0;
  mrpc_trans_with_ext_false_3=0;
  mrpc_trans_with_dbs_true_3=0;
  mrpc_trans_with_dbs_false_3=0;
 
  mrpc_trans_with_ext_true_4=0;
  mrpc_trans_with_ext_false_4=0;
  mrpc_trans_with_dbs_true_4=0;
  mrpc_trans_with_dbs_false_4=0;
  return;
}
 
 
void MrpcDBSCount::setExtTrackNum( unsigned int ntrk ) {
  extTrkNum = extTrkNum + ntrk;
  return;
}
 
 
void MrpcDBSCount::setTrack1( MrpcDBSTrack*& 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 MrpcDBSCount::setTrack1Col( std::vector<MrpcDBSTrack*>*& tofTrackVec ) {
  TofTrackVec::iterator iter_tof = tofTrackVec->begin();
  for( ; iter_tof!=tofTrackVec->end(); iter_tof++ ) {
    setTrack1( *iter_tof );
  }
  return;
}
 
 
void MrpcDBSCount::setTrack2( MrpcDBSTrack*& 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 MrpcDBSCount::setTrack3( MrpcDBSTrack*& 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->trans_with_extrap_track()) mrpc_trans_with_ext_true_3++;
      else if(tof->trans_with_dbs_info() && !(tof->trans_with_extrap_track())) mrpc_trans_with_dbs_true_3++;
      else mrpc_trans_with_dbs_false_3++;
      
 
 
    }
 
  return;
}
 
 
void MrpcDBSCount::setTrack4( MrpcDBSTrack*& 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->hitCase() == eastendcap_1 )  ||  ( tof->hitCase() == eastendcap_2 )  ||  ( tof->hitCase() == westendcap_1 )  ||   ( tof->hitCase() == westendcap_2 ))
    {
 
      if(tof->trans_with_extrap_track()) mrpc_trans_with_ext_true_4++;
      else if(tof->trans_with_dbs_info() && !(tof->trans_with_extrap_track())) mrpc_trans_with_dbs_true_4++;
      else mrpc_trans_with_dbs_false_4++;
      
 
 
    }
  
  return;
}
 
void MrpcDBSCount::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;
   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_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3) >0)
     {
   std::cout<<"  ==> MRPC extrapolation (before calibration):              " << endl;
   std::cout<<"      Determine transmision time without extrapolation :    "<< setprecision(3)  << (double)mrpc_trans_with_dbs_true_3*100.0/((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*1.0)<< " % "<<endl ;
   std::cout<<"      Determine transmision time with extrapolation :       "<< setprecision(3)  << (double)mrpc_trans_with_ext_true_3*100.0/((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*1.0)<< " % "<<endl ;
   std::cout<<"      Determine transmision time with average value :       "<< setprecision(3)  << 
     ((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)-mrpc_trans_with_dbs_true_3-mrpc_trans_with_ext_true_3)*100.0/((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*1.0)<< " % "<<endl ;
     }
   if((mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4) >0)
     {
   std::cout<<"  ==> MRPC extrapolation (after calibration):               " << endl;
   std::cout<<"      Determine transmision time without extrapolation :    "<< setprecision(3)  << (double)mrpc_trans_with_dbs_true_4*100.0/((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)*1.0)<< " % "<<endl ;
   std::cout<<"      Determine transmision time with extrapolation :       "<< setprecision(3)  << (double)mrpc_trans_with_ext_true_4*100.0/((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)*1.0)<< " % "<<endl ;
   std::cout<<"      Determine transmision time with average value :       "<< setprecision(3)  << 
     ((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)-mrpc_trans_with_dbs_true_4-mrpc_trans_with_ext_true_4)*100.0/((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)*1.0)<< " % "<<endl ;
     }
 
 
 
   std::cout<<"#############################  End Summary #######################################" <<std::endl;
 
}