SiStripPI Namespace Reference

Enumerations
enum	estimator { min, max, mean, rms }
enum	palette {   HALFGRAY, GRAY, BLUES, REDS,   ANTIGRAY, FIRE, ANTIFIRE, LOGREDBLUE,   LOGBLUERED, DEFAULT }
enum	TrackerRegion {   TIB1r = 1010, TIB1s = 1011, TIB2r = 1020, TIB2s = 1021,   TIB3r = 1030, TIB4r = 1040, TOB1r = 2010, TOB1s = 2011,   TOB2r = 2020, TOB2s = 2021, TOB3r = 2030, TOB4r = 2040,   TOB5r = 2050, TOB6r = 2060, TEC1r = 3010, TEC1s = 3011,   TEC2r = 3020, TEC2s = 3021, TEC3r = 3030, TEC3s = 3031,   TEC4r = 3040, TEC4s = 3041, TEC5r = 3050, TEC5s = 3051,   TEC6r = 3060, TEC6s = 3061, TEC7r = 3070, TEC7s = 3071,   TEC8r = 3080, TEC8s = 3081, TEC9r = 3090, TEC9s = 3091,   TID1r = 4010, TID1s = 4011, TID2r = 4020, TID2s = 4021,   TID3r = 4030, TID3s = 4031, END_OF_REGIONS }

Functions
void	drawStatBox (std::map< std::string, std::shared_ptr< TH1F >> histos, std::map< std::string, int > colormap, std::vector< std::string > legend, double X=0.15, double Y=0.93, double W=0.15, double H=0.10)
std::string	estimatorType (SiStripPI::estimator e)
std::pair< float, float >	getExtrema (TH1 *h1, TH1 *h2)
std::string	getStringFromSubdet (StripSubdetector::SubDetector sub)
std::pair< float, float >	getTheRange (std::map< uint32_t, float > values, const float nsigma)
void	makeNicePlotStyle (TH1 *hist)
void	printSummary (const std::map< unsigned int, SiStripDetSummary::Values > &map)
std::pair< int, const char * >	regionType (int index)
void	setBadComponents (int i, int component, SiStripQuality::BadComponent &BC, int NBadComponent[4][19][4])
void	setPaletteStyle (SiStripPI::palette palette)

Enumeration Type Documentation

enum SiStripPI::estimator

Enumerator
min
max
mean
rms

Definition at line 17 of file SiStripPayloadInspectorHelper.h.

                 {
    min,
    max,
    mean,
    rms
  };

enum SiStripPI::palette

Enumerator
HALFGRAY
GRAY
BLUES
REDS
ANTIGRAY
FIRE
ANTIFIRE
LOGREDBLUE
LOGBLUERED
DEFAULT

Definition at line 290 of file SiStripPayloadInspectorHelper.h.

{HALFGRAY,GRAY,BLUES,REDS,ANTIGRAY,FIRE,ANTIFIRE,LOGREDBLUE,LOGBLUERED,DEFAULT};

enum SiStripPI::TrackerRegion

Enumerator
TIB1r
TIB1s
TIB2r
TIB2s
TIB3r
TIB4r
TOB1r
TOB1s
TOB2r
TOB2s
TOB3r
TOB4r
TOB5r
TOB6r
TEC1r
TEC1s
TEC2r
TEC2s
TEC3r
TEC3s
TEC4r
TEC4s
TEC5r
TEC5s
TEC6r
TEC6s
TEC7r
TEC7s
TEC8r
TEC8s
TEC9r
TEC9s
TID1r
TID1s
TID2r
TID2s
TID3r
TID3s
END_OF_REGIONS

Definition at line 50 of file SiStripPayloadInspectorHelper.h.

                     { 
    TIB1r = 1010, TIB1s = 1011, 
    TIB2r = 1020, TIB2s = 1021,
    TIB3r = 1030,
    TIB4r = 1040,
    TOB1r = 2010, TOB1s = 2011,
    TOB2r = 2020, TOB2s = 2021,
    TOB3r = 2030,
    TOB4r = 2040,
    TOB5r = 2050,
    TOB6r = 2060,
    TEC1r = 3010, TEC1s = 3011,
    TEC2r = 3020, TEC2s = 3021,
    TEC3r = 3030, TEC3s = 3031,
    TEC4r = 3040, TEC4s = 3041,
    TEC5r = 3050, TEC5s = 3051,
    TEC6r = 3060, TEC6s = 3061,
    TEC7r = 3070, TEC7s = 3071,
    TEC8r = 3080, TEC8s = 3081,
    TEC9r = 3090, TEC9s = 3091,
    TID1r = 4010, TID1s = 4011,
    TID2r = 4020, TID2s = 4021,
    TID3r = 4030, TID3s = 4031,
    END_OF_REGIONS  
  };

Function Documentation

void SiStripPI::drawStatBox	(	std::map< std::string, std::shared_ptr< TH1F >>	histos,
		std::map< std::string, int >	colormap,
		std::vector< std::string >	legend,
		double	X = 0.15,
		double	Y = 0.93,
		double	W = 0.15,
		double	H = 0.10
	)

Definition at line 157 of file SiStripPayloadInspectorHelper.h.

References edmScanValgrind::buffer, class-composition::H, plotFactory::histos, mps_fire::i, X, and DOFs::Y.

  {  
    char   buffer[255];
    
    int i=0;
    for ( const auto &element : legend ){
      TPaveText* stat = new TPaveText(X,Y-(i*H), X+W, Y-(i+1)*H, "NDC");
      i++;
      auto Histo = histos[element];
      sprintf(buffer,"Entries : %i\n",(int)Histo->GetEntries());
      stat->AddText(buffer);
      
      sprintf(buffer,"Mean    : %6.2f\n",Histo->GetMean());
      stat->AddText(buffer);
      
      sprintf(buffer,"RMS     : %6.2f\n",Histo->GetRMS());
      stat->AddText(buffer);
      
      stat->SetFillColor(0);
      stat->SetLineColor(colormap[element]);
      stat->SetTextColor(colormap[element]);
      stat->SetTextSize(0.03);
      stat->SetBorderSize(0);
      stat->SetMargin(0.05);
      stat->SetTextAlign(12);
      stat->Draw();
    }
  }

std::string SiStripPI::estimatorType

(

SiStripPI::estimator

e

)

Definition at line 25 of file SiStripPayloadInspectorHelper.h.

References max, mean, min, and rms.

  {
    switch(e){
    case SiStripPI::min : return "minimum";
    case SiStripPI::max : return "maximum";
    case SiStripPI::mean : return "mean";
    case SiStripPI::rms  : return "RMS";
    default: return "should never be here";
    }
  }

std::pair<float,float> SiStripPI::getExtrema	(	TH1 *	h1,
		TH1 *	h2
	)

Definition at line 188 of file SiStripPayloadInspectorHelper.h.

References mps_fire::result.

  {
    float theMax(-9999.);
    float theMin(9999.);
    theMax =  h1->GetMaximum() > h2->GetMaximum() ?  h1->GetMaximum() :  h2->GetMaximum();
    theMin =  h1->GetMinimum() < h2->GetMaximum() ?  h1->GetMinimum() :  h2->GetMinimum();
    
    float add_min = theMin>0. ? -0.05 :  0.05;
    float add_max = theMax>0. ?  0.05 : -0.05;

    auto result = std::make_pair(theMin*(1+add_min),theMax*(1+add_max));
    return result;
    
  }

std::string SiStripPI::getStringFromSubdet

(

StripSubdetector::SubDetector

sub

)

Definition at line 38 of file SiStripPayloadInspectorHelper.h.

References StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

  {
    switch(sub){
    case StripSubdetector::TIB : return "TIB";
    case StripSubdetector::TOB : return "TOB";
    case StripSubdetector::TID : return "TID";
    case StripSubdetector::TEC : return "TEC";
    default : return "should never be here"; 
    }
  }

std::pair<float,float> SiStripPI::getTheRange	(	std::map< uint32_t, float >	values,
		const float	nsigma
	)

Definition at line 128 of file SiStripPayloadInspectorHelper.h.

References begin, end, funct::m, AlCaHLTBitMon_ParallelJobs::p, mathSSE::sqrt(), and plotscripts::stdev().

  {
    float sum = std::accumulate(std::begin(values), 
                std::end(values), 
                0.0,
                [] (float value, const std::map<uint32_t,float>::value_type& p)
                { return value + p.second; }
                );
    
    float m =  sum / values.size();
    
    float accum = 0.0;
    std::for_each (std::begin(values), 
           std::end(values), 
           [&](const std::map<uint32_t,float>::value_type& p) 
           {accum += (p.second - m) * (p.second - m);}
           );
    
    float stdev = sqrt(accum / (values.size()-1)); 
    
    if(stdev!=0.){
      return std::make_pair(m-nsigma*stdev,m+nsigma*stdev);
    } else {
      return std::make_pair(m>0.? 0.95*m : 1.05*m, m>0? 1.05*m : 0.95*m);
    }
  }

void SiStripPI::makeNicePlotStyle

(

TH1 *

hist

)

Definition at line 206 of file SiStripPayloadInspectorHelper.h.

  { 
    hist->SetStats(kFALSE);  
    hist->SetLineWidth(2);
    hist->GetXaxis()->CenterTitle(true);
    hist->GetYaxis()->CenterTitle(true);
    hist->GetXaxis()->SetTitleFont(42); 
    hist->GetYaxis()->SetTitleFont(42);  
    hist->GetXaxis()->SetTitleSize(0.05);
    hist->GetYaxis()->SetTitleSize(0.05);
    hist->GetXaxis()->SetTitleOffset(0.9);
    hist->GetYaxis()->SetTitleOffset(1.3);
    hist->GetXaxis()->SetLabelFont(42);
    hist->GetYaxis()->SetLabelFont(42);
    hist->GetYaxis()->SetLabelSize(.05);
    hist->GetXaxis()->SetLabelSize(.05);
  }

void SiStripPI::printSummary

(

const std::map< unsigned int, SiStripDetSummary::Values > &

map

)

Definition at line 227 of file SiStripPayloadInspectorHelper.h.

References KineDebug3::count(), gather_cfg::cout, gamEcalExtractorBlocks_cff::detector, mean, regionType(), rms, mathSSE::sqrt(), and AlCaHLTBitMon_QueryRunRegistry::string.

  {
    for (const auto &element : map){
      int count   = element.second.count;
      double mean = count>0 ? (element.second.mean)/count : 0. ;
      double rms  = count>0 ? (element.second.rms)/count - mean*mean : 0.;
      if(rms <= 0)
    rms = 0;
      else
    rms = sqrt(rms);
      
      std::string detector;
      
      switch ((element.first)/1000) 
    {
    case 1:
      detector = "TIB ";
      break;
    case 2:
      detector = "TOB ";
      break;
    case 3:
      detector = "TEC ";
      break;
    case 4:
      detector = "TID ";
      break;
    }
      
      int layer  = (element.first)/10 - (element.first)/1000*100;
      int stereo = (element.first) - (layer*10) -(element.first)/1000*1000;
      
      std::cout<<"key of the map:"<<element.first <<" ( region: "<<regionType(element.first).second <<" ) "  
           << detector<<" layer: "<<layer<<" stereo:"<<stereo
           <<"| count:"<<count<<" mean: "<<mean<<" rms: "<<rms<<std::endl;
      
    }
  }

std::pair<int,const char *> SiStripPI::regionType

(

int

index

)

Definition at line 77 of file SiStripPayloadInspectorHelper.h.

References END_OF_REGIONS, diffTreeTool::index, TEC1r, TEC1s, TEC2r, TEC2s, TEC3r, TEC3s, TEC4r, TEC4s, TEC5r, TEC5s, TEC6r, TEC6s, TEC7r, TEC7s, TEC8r, TEC8s, TEC9r, TEC9s, TIB1r, TIB1s, TIB2r, TIB2s, TIB3r, TIB4r, TID1r, TID1s, TID2r, TID2s, TID3r, TID3s, TOB1r, TOB1s, TOB2r, TOB2s, TOB3r, TOB4r, TOB5r, and TOB6r.

Referenced by printSummary().

  {
    
    auto region = static_cast<std::underlying_type_t<SiStripPI::TrackerRegion> >(index);

    switch(region){
    case SiStripPI::TIB1r: return std::make_pair(1 ,"TIB L1 r-#varphi");
    case SiStripPI::TIB1s: return std::make_pair(2 ,"TIB L1 stereo");
    case SiStripPI::TIB2r: return std::make_pair(3 ,"TIB L2 r-#varphi");
    case SiStripPI::TIB2s: return std::make_pair(4 ,"TIB L2 stereo");
    case SiStripPI::TIB3r: return std::make_pair(5 ,"TIB L3");
    case SiStripPI::TIB4r: return std::make_pair(6 ,"TIB L4");
    case SiStripPI::TOB1r: return std::make_pair(7 ,"TOB L1 r-#varphi");
    case SiStripPI::TOB1s: return std::make_pair(8 ,"TOB L1 stereo");
    case SiStripPI::TOB2r: return std::make_pair(9 ,"TOB L2 r-#varphi");
    case SiStripPI::TOB2s: return std::make_pair(10,"TOB L2 stereo");
    case SiStripPI::TOB3r: return std::make_pair(11,"TOB L3 r-#varphi");
    case SiStripPI::TOB4r: return std::make_pair(12,"TOB L4");
    case SiStripPI::TOB5r: return std::make_pair(13,"TOB L5");
    case SiStripPI::TOB6r: return std::make_pair(14,"TOB L6");
    case SiStripPI::TEC1r: return std::make_pair(15,"TEC D1 r-#varphi");
    case SiStripPI::TEC1s: return std::make_pair(16,"TEC D1 stereo");
    case SiStripPI::TEC2r: return std::make_pair(17,"TEC D2 r-#varphi");
    case SiStripPI::TEC2s: return std::make_pair(18,"TEC D2 stereo");
    case SiStripPI::TEC3r: return std::make_pair(19,"TEC D3 r-#varphi");
    case SiStripPI::TEC3s: return std::make_pair(20,"TEC D3 stereo");
    case SiStripPI::TEC4r: return std::make_pair(21,"TEC D4 r-#varphi");
    case SiStripPI::TEC4s: return std::make_pair(22,"TEC D4 stereo");
    case SiStripPI::TEC5r: return std::make_pair(23,"TEC D5 r-#varphi");
    case SiStripPI::TEC5s: return std::make_pair(24,"TEC D5 stereo");
    case SiStripPI::TEC6r: return std::make_pair(25,"TEC D6 r-#varphi");
    case SiStripPI::TEC6s: return std::make_pair(26,"TEC D6 stereo");
    case SiStripPI::TEC7r: return std::make_pair(27,"TEC D7 r-#varphi");
    case SiStripPI::TEC7s: return std::make_pair(28,"TEC D7 stereo");
    case SiStripPI::TEC8r: return std::make_pair(29,"TEC D8 r-#varphi");
    case SiStripPI::TEC8s: return std::make_pair(30,"TEC D8 stereo");
    case SiStripPI::TEC9r: return std::make_pair(31,"TEC D9 r-#varphi");
    case SiStripPI::TEC9s: return std::make_pair(32,"TEC D9 stereo");
    case SiStripPI::TID1r: return std::make_pair(33,"TID D1 r-#varphi");
    case SiStripPI::TID1s: return std::make_pair(34,"TID D1 stereo");
    case SiStripPI::TID2r: return std::make_pair(35,"TID D2 r-#varphi");
    case SiStripPI::TID2s: return std::make_pair(36,"TID D2 stereo");
    case SiStripPI::TID3r: return std::make_pair(37,"TID D3 r-#varphi"); 
    case SiStripPI::TID3s: return std::make_pair(38,"TID D3 stereo");
    case SiStripPI::END_OF_REGIONS : std::make_pair(-1,"undefined");
    default : return std::make_pair(999,"should never be here");  
    }
  }

void SiStripPI::setBadComponents	(	int	i,
		int	component,
		SiStripQuality::BadComponent &	BC,
		int	NBadComponent[4][19][4]
	)

Definition at line 270 of file SiStripPayloadInspectorHelper.h.

References SiStripQuality::BadComponent::BadApvs, SiStripQuality::BadComponent::BadFibers, SiStripQuality::BadComponent::BadModule, KineDebug3::count(), and mps_fire::i.

  {
   
    if (BC.BadApvs){
      NBadComponent[i][0][2]+= std::bitset<16>(BC.BadApvs&0x3f).count(); 
      NBadComponent[i][component][2]+= std::bitset<16>(BC.BadApvs&0x3f).count(); 
    }

    if (BC.BadFibers){ 
      NBadComponent[i][0][1]+= std::bitset<4>(BC.BadFibers&0x7).count();
      NBadComponent[i][component][1]+= std::bitset<4>(BC.BadFibers&0x7).count();
    }   

    if (BC.BadModule){
      NBadComponent[i][0][0]++;
      NBadComponent[i][component][0]++;
    }
  }

void SiStripPI::setPaletteStyle

(

SiStripPI::palette

palette

)

Definition at line 293 of file SiStripPayloadInspectorHelper.h.

References ANTIFIRE, ANTIGRAY, BLUES, gather_cfg::cout, FIRE, GRAY, cmsRelvalreport::green(), HALFGRAY, LOGBLUERED, LOGREDBLUE, cmsRelvalreport::red(), and REDS.

  {
  
    TStyle *palettestyle = new TStyle("palettestyle","Style for P-TDR");
  
    const int NRGBs = 5;
    const int NCont = 255;
    
    switch(palette){

    case HALFGRAY:
      {
    double stops[NRGBs] = {0.00, 0.34, 0.61, 0.84, 1.00};
    double red[NRGBs]   = {1.00, 0.91, 0.80, 0.67, 1.00};
    double green[NRGBs] = {1.00, 0.91, 0.80, 0.67, 1.00};
    double blue[NRGBs]  = {1.00, 0.91, 0.80, 0.67, 1.00};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
      }
      break;

    case GRAY:
      {
    double stops[NRGBs] = {0.00, 0.01, 0.05, 0.09, 0.1};
    double red[NRGBs]   = {1.00, 0.84, 0.61, 0.34, 0.00};
    double green[NRGBs] = {1.00, 0.84, 0.61, 0.34, 0.00};
    double blue[NRGBs]  = {1.00, 0.84, 0.61, 0.34, 0.00};
        TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
      }
      break;

    case BLUES:
      {
    double stops[NRGBs] = {0.00, 0.34, 0.61, 0.84, 1.00};
    double red[NRGBs]   = {1.00, 0.84, 0.61, 0.34, 0.00};
    double green[NRGBs] = {1.00, 0.84, 0.61, 0.34, 0.00};
    double blue[NRGBs]  = {1.00, 1.00, 1.00, 1.00, 1.00};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
    
      }
      break;
      
    case REDS:
    {
      double stops[NRGBs] = {0.00, 0.34, 0.61, 0.84, 1.00};
      double red[NRGBs]   = {1.00, 1.00, 1.00, 1.00, 1.00};
      double green[NRGBs] = {1.00, 0.84, 0.61, 0.34, 0.00};
      double blue[NRGBs]  = {1.00, 0.84, 0.61, 0.34, 0.00};
      TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);  
    }
    break;
      
    case ANTIGRAY:
      {
    double stops[NRGBs] = {0.00, 0.34, 0.61, 0.84, 1.00};
    double red[NRGBs]   = {0.00, 0.34, 0.61, 0.84, 1.00};
    double green[NRGBs] = {0.00, 0.34, 0.61, 0.84, 1.00};
    double blue[NRGBs]  = {0.00, 0.34, 0.61, 0.84, 1.00};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);        
      }      
      break;
      
    case FIRE:
      {
    double stops[NRGBs] = {0.00, 0.20, 0.80, 1.00};
    double red[NRGBs]   = {1.00, 1.00, 1.00, 0.50};
    double green[NRGBs] = {1.00, 1.00, 0.00, 0.00};
    double blue[NRGBs]  = {0.20, 0.00, 0.00, 0.00};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);    
      }
      break;

    case ANTIFIRE:
      {
    double stops[NRGBs] = {0.00, 0.20, 0.80, 1.00};
    double red[NRGBs]   = {0.50, 1.00, 1.00, 1.00};
    double green[NRGBs] = {0.00, 0.00, 1.00, 1.00};
    double blue[NRGBs]  = {0.00, 0.00, 0.00, 0.20};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);        
      }
      break;

    case LOGREDBLUE:
      {
    double stops[NRGBs] = {0.0001, 0.0010, 0.0100, 0.1000,  1.0000};
    double red[NRGBs]   = {1.00,   0.75,   0.50,   0.25,    0.00};
    double green[NRGBs] = {0.00,   0.00,   0.00,   0.00,    0.00};
    double blue[NRGBs]  = {0.00,   0.25,   0.50,   0.75,    1.00};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);        
      }
      break;
      
    case LOGBLUERED:
      {
    double stops[NRGBs] = {0.0001, 0.0010, 0.0100, 0.1000,  1.0000};
    double red[NRGBs]   = {0.00,   0.25,   0.50,   0.75,    1.00};
    double green[NRGBs] = {0.00,   0.00,   0.00,   0.00,    0.00};
    double blue[NRGBs]  = {1.00,   0.75,   0.50,   0.25,    0.00};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);            
      } 
      break;
    
    case DEFAULT:
      {
    double stops[NRGBs] = {0.00, 0.34, 0.61, 0.84, 1.00};
    double red[NRGBs]   = {0.00, 0.00, 0.87, 1.00, 0.51};
    double green[NRGBs] = {0.00, 0.81, 1.00, 0.20, 0.00};
    double blue[NRGBs]  = {0.51, 1.00, 0.12, 0.00, 0.00};
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);    
      }
      break;
    default:
      std::cout<<"should nevere be here" << std::endl;
      break;
    }
    
    palettestyle->SetNumberContours(NCont);
  }