SiPixelGenError.h

Go to the documentation of this file.

//
//  SiPixelGenError.h (v1.00)
//
//  Object to contain Lorentz drift and error information for the Generic Algorithm
//
// Created by Morris Swartz on 1/10/2014.
//
//
 
// Build the template storage structure from several pieces 

#ifndef SiPixelGenError_h
#define SiPixelGenError_h 1


#include<vector>
#include<cassert>

#ifndef SI_PIXEL_TEMPLATE_STANDALONE
#include "CondFormats/SiPixelObjects/interface/SiPixelGenErrorDBObject.h"
#include "FWCore/Utilities/interface/Exception.h"
#endif

struct SiPixelGenErrorEntry { 
  int runnum;              
  float cotalpha;          
  float cotbeta;           
  float qavg;              
  float pixmax;            
  float dyone;             
  float syone;             
  float dxone;             
  float sxone;             
  float dytwo;             
  float sytwo;             
  float dxtwo;             
  float sxtwo;             
  float qmin;              
  float qmin2;
  float yavggen[4];        
  float yrmsgen[4];        
  float xavggen[4];        
  float xrmsgen[4];        
} ;





struct SiPixelGenErrorHeader {           
   int ID;                 
   int NTy;                
   int NTyx;               
   int NTxx;               
   int Dtype;              
   float qscale;           
   float lorywidth;        
   float lorxwidth;        
   float lorybias;         
   float lorxbias;         
   float Vbias;            
   float temperature;      
   float fluence;          
   float s50;              
   float ss50;             
   char title[80];         
   int templ_version;      
   float Bfield;           
   float fbin[3];          
   float xsize;            
   float ysize;            
   float zsize;            
} ;



struct SiPixelGenErrorStore { 
  SiPixelGenErrorHeader head;
  float cotbetaY[60];
  float cotbetaX[5];
  float cotalphaX[29];
  SiPixelGenErrorEntry enty[60];     
  SiPixelGenErrorEntry entx[5][29];  
} ;


// ******************************************************************************************
// ******************************************************************************************
class SiPixelGenError {
 public:
  SiPixelGenError(const std::vector< SiPixelGenErrorStore > & thePixelTemp) : thePixelTemp_(thePixelTemp) { id_current_ = -1; index_id_ = -1;} 
  
  static bool pushfile(int filenum, std::vector< SiPixelGenErrorStore > & thePixelTemp_);     // load the private store with info from the 
                                  // file with the index (int) filenum
                                  
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
  static bool pushfile(const SiPixelGenErrorDBObject& dbobject, std::vector< SiPixelGenErrorStore > & thePixelTemp_);     // load the private store with info from db
#endif

// initialize the binary search information;
  static void postInit(std::vector< SiPixelGenErrorStore > & thePixelTemp_);
  
    
// Interpolate input beta angle to estimate the average charge. return qbin flag for input cluster charge, and estimate y/x errors and biases for the Generic Algorithm.
  int qbin(int id, float cotalpha, float cotbeta, float locBz, float qclus, float& pixmx, float& sigmay, float& deltay, float& sigmax, float& deltax,
            float& sy1, float& dy1, float& sy2, float& dy2, float& sx1, float& dx1, float& sx2, float& dx2);
  // Overloaded method to provide only the LA parameters
  int qbin(int id);

  float lorywidth() {return lorywidth_;}                            
  float lorxwidth() {return lorxwidth_;}                            
  float lorybias() {return lorybias_;}                              
  float lorxbias() {return lorxbias_;}                              

  float fbin(int i) {
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
      if(i < 0 || i > 2) {throw cms::Exception("DataCorrupt") << "SiPixelTemplate::chi2xminc2m called with illegal index = " << i << std::endl;}
#else
      assert(i>=0 && i<3);
#endif
  return fbin_[i];} 
  float xsize() {return xsize_;}                                    
  float ysize() {return ysize_;}                                    
  float zsize() {return zsize_;}                                    

 private:
  
  // Keep current template interpolaion parameters  
  
  int id_current_;           
  int index_id_;             
  
  
  // Keep results of last interpolation to return through member functions
  
  float lorxwidth_;         
  float lorywidth_;         
  float lorxbias_;         
  float lorybias_;         
  float fbin_[3];          
  float xsize_;             
  float ysize_;             
  float zsize_;             

  
  // The actual template store is a std::vector container

  const std::vector< SiPixelGenErrorStore > & thePixelTemp_;
} ;


#endif