/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/RecoLocalTracker/SiPixelRecHits/interface/PixelCPEGeneric.h

Go to the documentation of this file.

#ifndef RecoLocalTracker_SiPixelRecHits_PixelCPEGeneric_H
#define RecoLocalTracker_SiPixelRecHits_PixelCPEGeneric_H

// \class PixelCPEGeneric  -- a generalized CPE reco for the idealized detector
//
// The basic idea of this class is to use generic formulae in order
// to achieve clean and minimal code.  It should work for
// - both normal and big pixels
// - both barrel and forward
// - both "FromDetPosition" and "FromTrackAngles" (i.e. by the track fit)
//
// This is possible since, in its nature, the original "ORCA" algorithm by 
// Danek and Susana is the same in both X and Y directions, provided that
// one correctly computes angles alpha_ and beta_ up front.  Thus, all
// geometrical and special corrections are dropped, since the presumption
// is that alpha_ and beta_ are determined as best as possible.  That means
// that they either come from the track, or, if they come from the 
// position of the DetUnit, they include all geometrical information 
// possible for this DetUnit:
// - for both the barrel and the forward, we use the cluster position 
//   instead of the center of the module/plaquette
// - for the forward, the tilt of the blades is included too
//
// In addtion, anything which is special for the computation of the lorentz
// angle is done in setTheDet() method.  So the algorithm per se does not
// need to worry about it.  This includes extra E*B term (a.k.a. "alpha2Order")
// and extra tilt in the forward.
//
// Thus, the formula for the computation of the hit position is very
// simple, and is described in Morris's note (IN ???) on the generalizaton
// of the pixel algorithm.

#include "RecoLocalTracker/SiPixelRecHits/interface/PixelCPEBase.h"
#include "CalibTracker/SiPixelESProducers/interface/SiPixelCPEGenericDBErrorParametrization.h"

// Already defined in the base class
//#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
//#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
//#include "Geometry/TrackerGeometryBuilder/interface/RectangularPixelTopology.h"
//#include "Geometry/CommonDetAlgo/interface/MeasurementPoint.h"
//#include "Geometry/CommonDetAlgo/interface/MeasurementError.h"
//#include "Geometry/Surface/interface/GloballyPositioned.h"
//#include "FWCore/ParameterSet/interface/ParameterSet.h"

// The template header files
#include "RecoLocalTracker/SiPixelRecHits/interface/SiPixelTemplateReco.h"
#include "RecoLocalTracker/SiPixelRecHits/interface/SiPixelTemplate.h"


#include <utility>
#include <vector>


#if 0

#endif

class MagneticField;
class PixelCPEGeneric : public PixelCPEBase
{
 public:
  // PixelCPEGeneric( const DetUnit& det );
  PixelCPEGeneric(edm::ParameterSet const& conf, const MagneticField *, const SiPixelLorentzAngle *, const SiPixelCPEGenericErrorParm *, const SiPixelTemplateDBObject *);
  ~PixelCPEGeneric() {;}

  LocalPoint localPosition (const SiPixelCluster& cluster, const GeomDetUnit & det) const; 
  
  // However, we do need to implement localError().
  LocalError localError   (const SiPixelCluster& cl, const GeomDetUnit & det) const;
  
  MeasurementPoint measurementPosition ( const SiPixelCluster&, 
                                         const GeomDetUnit & det) const;
/*   MeasurementError measurementError    ( const SiPixelCluster&,  */
/*                                        const GeomDetUnit & det) const; */


 private:
  //--------------------------------------------------------------------
  //  Methods.
  //------------------------------------------------------------------
  double
    generic_position_formula( int size,                
                              double Q_f,              
                              double Q_l,              
                              double upper_edge_first_pix, 
                              double lower_edge_last_pix,  
                              double half_lorentz_shift,   
                              double cot_angle,            
                              double pitch,            
                              bool first_is_big,       
                              bool last_is_big,        
                              double eff_charge_cut_low, 
                              double eff_charge_cut_high,
                              double size_cut           
                              ) const;

  void
    collect_edge_charges(const SiPixelCluster& cluster,  
                         float & Q_f_X,              
                         float & Q_l_X,              
                         float & Q_m_X,              
                         float & Q_f_Y,              
                         float & Q_l_Y,              
                         float & Q_m_Y               
                         ) const;
  
  
  //--- Errors squared in x and y.  &&& Need to be revisited.
  float err2X(bool&, int&) const;
  float err2Y(bool&, int&) const;

  //--- Cuts made externally settable
  double the_eff_charge_cut_lowX;
  double the_eff_charge_cut_lowY;
  double the_eff_charge_cut_highX;
  double the_eff_charge_cut_highY;
  double the_size_cutX;
  double the_size_cutY;

  bool inflate_errors;
  bool inflate_all_errors_no_trk_angle;

  bool UseErrorsFromTemplates_;
  bool DoCosmics_;
  bool LoadTemplatesFromDB_;
  bool TruncatePixelCharge_;
  bool IrradiationBiasCorrection_;

  double EdgeClusterErrorX_;
  double EdgeClusterErrorY_;

  //--- DB Error Parametrization object
  SiPixelCPEGenericDBErrorParametrization * genErrorsFromDB_;

  mutable SiPixelTemplate templ_;
  mutable int templID_; 

        // The truncation value pix_maximum is an angle-dependent cutoff on the
        // individual pixel signals. It should be applied to all pixels in the
        // cluster [signal_i = fminf(signal_i, pixmax)] before the column and row
        // sums are made. Morris
        mutable float pixmx;

        // These are errors predicted by PIXELAV
        mutable float sigmay; // CPE Generic y-error for multi-pixel cluster
        mutable float sigmax; // CPE Generic x-error for multi-pixel cluster
        mutable float sy1   ; // CPE Generic y-error for single single-pixel
        mutable float sy2   ; // CPE Generic y-error for single double-pixel cluster
        mutable float sx1   ; // CPE Generic x-error for single single-pixel cluster
        mutable float sx2   ; // CPE Generic x-error for single double-pixel cluster

        // These are irradiation bias corrections
        mutable float deltay; // CPE Generic y-bias for multi-pixel cluster
        mutable float deltax; // CPE Generic x-bias for multi-pixel cluster
        mutable float dy1   ; // CPE Generic y-bias for single single-pixel cluster
        mutable float dy2   ; // CPE Generic y-bias for single double-pixel cluster
        mutable float dx1   ; // CPE Generic x-bias for single single-pixel cluster
        mutable float dx2   ; // CPE Generic x-bias for single double-pixel cluster

         
 protected:
  //--- These functions are no longer needed, yet they are declared 
  //--- pure virtual in the base class.
  float xpos( const SiPixelCluster& ) const { return -999000.0; }  // &&& should abort
  float ypos( const SiPixelCluster& ) const { return -999000.0; }  // &&& should abort

};

#endif