SiPixelErrorEstimation Class Reference

#include <SiPixelErrorEstimation.h>

Inheritance diagram for SiPixelErrorEstimation:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
void	computeAnglesFromDetPosition (const SiPixelCluster &cl, const GeomDetUnit &det, float &alpha, float &beta)
virtual void	endJob ()
	SiPixelErrorEstimation (const edm::ParameterSet &)
virtual	~SiPixelErrorEstimation ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
float	all_alpha
float	all_beta
int	all_bigx
int	all_bigy
int	all_blade
float	all_clust_alpha
float	all_clust_beta
float	all_clust_col
int	all_clust_geoid
int	all_clust_maxpixcol
int	all_clust_maxpixrow
int	all_clust_minpixcol
int	all_clust_minpixrow
float	all_clust_q
float	all_clust_row
float	all_clust_x
float	all_clust_y
float	all_col1
float	all_col2
int	all_cols
int	all_disk
int	all_edgex
int	all_edgey
float	all_eloss
int	all_flipped
float	all_gx1
float	all_gx2
float	all_gy1
float	all_gy2
float	all_gz1
float	all_gz2
int	all_half
float	all_hit_cprob0
float	all_hit_cprob1
float	all_hit_cprob2
float	all_hit_probx
float	all_hit_proby
int	all_ladder
int	all_layer
int	all_mod
int	all_npix
int	all_nsimhit
int	all_nxpix
int	all_nypix
int	all_panel
int	all_pidhit
float	all_pixadc [maxpix]
float	all_pixcol [maxpix]
float	all_pixel_clst_err_x
float	all_pixel_clst_err_y
int	all_pixel_split
float	all_pixgx [maxpix]
float	all_pixgy [maxpix]
float	all_pixgz [maxpix]
float	all_pixrow [maxpix]
float	all_pixx [maxpix]
float	all_pixy [maxpix]
int	all_plaq
float	all_rechiterrx
float	all_rechiterry
float	all_rechitpullx
float	all_rechitpully
float	all_rechitresx
float	all_rechitresy
float	all_rechitx
float	all_rechity
float	all_rechitz
float	all_row1
float	all_row2
int	all_rows
int	all_side
float	all_simhitx
float	all_simhity
float	all_simphi
int	all_simproc
float	all_simpx
float	all_simpy
float	all_simpz
float	all_simtheta
float	all_simtrketa
float	all_simtrkphi
int	all_subdetid
int	all_trkid
float	all_vtxr
float	all_vtxz
float	all_x1
float	all_x2
float	all_y1
float	all_y2
float	all_z1
float	all_z2
float	alpha
float	beta
int	bigx
int	bigy
int	blade
float	charge
bool	checkType_
edm::ParameterSet	conf_
int	detector_type
int	disk
int	edgex
int	edgey
float	eta
int	evt
int	flipped
int	genType_
int	half
float	hit_cprob0
float	hit_cprob1
float	hit_cprob2
float	hit_probx
float	hit_proby
bool	include_trk_hits_
int	ladder
int	layer
int	mod
int	npix
int	nsimhit
int	nxpix
int	nypix
std::string	outputFile_
int	panel
float	phi
int	pidhit
float	pixel_clst_err_x
float	pixel_clst_err_y
int	pixel_split
int	plaq
float	rechiterrx
float	rechiterry
float	rechitpullx
float	rechitpully
float	rechitresx
float	rechitresy
float	rechitx
float	rechity
float	rechitz
int	run
int	side
float	simhitx
float	simhity
int	simproc
std::string	src_
float	strip_charge
float	strip_clst_err_x
float	strip_clst_err_y
float	strip_cotalpha
float	strip_cotbeta
int	strip_edge
int	strip_hit_type
int	strip_is_stereo
float	strip_locbx
float	strip_locby
float	strip_locbz
int	strip_nprm
int	strip_nsimhit
int	strip_pidhit
int	strip_pidhit1
int	strip_pidhit2
int	strip_pidhit3
int	strip_pidhit4
int	strip_pidhit5
float	strip_prob
int	strip_qbin
float	strip_rechiterrx
float	strip_rechiterry
float	strip_rechitpullx
float	strip_rechitpully
float	strip_rechitresx
float	strip_rechitresx2
float	strip_rechitresy
float	strip_rechitx
float	strip_rechity
float	strip_rechitz
int	strip_simproc
int	strip_simproc1
int	strip_simproc2
int	strip_simproc3
int	strip_simproc4
int	strip_simproc5
int	strip_size
int	strip_split
int	strip_subdet_id
int	strip_tib_is_double_side
int	strip_tib_is_external_string
int	strip_tib_is_internal_string
int	strip_tib_is_rphi
int	strip_tib_is_stereo
int	strip_tib_is_z_minus_side
int	strip_tib_is_z_plus_side
int	strip_tib_layer
int	strip_tib_layer_number
int	strip_tib_module
int	strip_tib_module_number
int	strip_tib_order
int	strip_tib_side
int	strip_tib_string_number
int	strip_tob_is_double_side
int	strip_tob_is_rphi
int	strip_tob_is_stereo
int	strip_tob_is_z_minus_side
int	strip_tob_is_z_plus_side
int	strip_tob_layer
int	strip_tob_layer_number
int	strip_tob_module
int	strip_tob_module_number
int	strip_tob_rod_number
int	strip_tob_side
float	strip_trk_pt
int	subdetId
TFile *	tfile_
edm::EDGetTokenT < SiPixelRecHitCollection >	tPixRecHitCollection
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
float	trk_alpha
float	trk_beta
edm::EDGetTokenT < edm::SimTrackContainer >	tSimTrackContainer
edm::EDGetTokenT < reco::TrackCollection >	tTrackCollection
edm::EDGetTokenT< std::vector < Trajectory > >	tTrajectory
TTree *	ttree_all_hits_
TTree *	ttree_track_hits_
TTree *	ttree_track_hits_strip_

Static Private Attributes
static const int	maxpix = 10000

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 63 of file SiPixelErrorEstimation.h.

Constructor & Destructor Documentation

SiPixelErrorEstimation::SiPixelErrorEstimation ( const edm::ParameterSet &  ps )

explicit

Definition at line 36 of file SiPixelErrorEstimation.cc.

References checkType_, genType_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), include_trk_hits_, HLT_25ns14e33_v1_cff::InputTag, outputFile_, src_, AlCaHLTBitMon_QueryRunRegistry::string, tPixRecHitCollection, tSimTrackContainer, tTrackCollection, and tTrajectory.

                                                                       :tfile_(0), ttree_all_hits_(0), 
                                        ttree_track_hits_(0), ttree_track_hits_strip_(0),
                                        trackerHitAssociatorConfig_(consumesCollector())
{
  //Read config file
  outputFile_ = ps.getUntrackedParameter<string>( "outputFile", "SiPixelErrorEstimation_Ntuple.root" );
  
  // Replace  "ctfWithMaterialTracks" with "generalTracks"
  //src_ = ps.getUntrackedParameter<std::string>( "src", "ctfWithMaterialTracks" );
  src_ = ps.getUntrackedParameter<std::string>( "src", "generalTracks" );

  checkType_ = ps.getParameter<bool>( "checkType" );
  genType_ = ps.getParameter<int>( "genType" );
  include_trk_hits_ = ps.getParameter<bool>( "include_trk_hits" );

  tTrajectory = consumes<std::vector<Trajectory>> (src_);
  tPixRecHitCollection = consumes<SiPixelRecHitCollection>(edm::InputTag( "siPixelRecHits"));
  tSimTrackContainer = consumes <edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
  tTrackCollection = consumes<reco::TrackCollection>(src_);

}

SiPixelErrorEstimation::~SiPixelErrorEstimation ( )

virtual

Definition at line 58 of file SiPixelErrorEstimation.cc.

{}

Member Function Documentation

void SiPixelErrorEstimation::analyze ( const edm::Event &  e, const edm::EventSetup &  es  )

virtual

Implements edm::EDAnalyzer.

Definition at line 372 of file SiPixelErrorEstimation.cc.

References funct::abs(), SiPixelCluster::Pixel::adc, all_alpha, all_beta, all_bigx, all_bigy, all_blade, all_clust_alpha, all_clust_beta, all_clust_col, all_clust_geoid, all_clust_maxpixcol, all_clust_maxpixrow, all_clust_minpixcol, all_clust_minpixrow, all_clust_q, all_clust_row, all_clust_x, all_clust_y, all_col1, all_col2, all_cols, all_disk, all_edgex, all_edgey, all_eloss, all_flipped, all_gx1, all_gx2, all_gy1, all_gy2, all_gz1, all_gz2, all_half, all_hit_cprob0, all_hit_cprob1, all_hit_cprob2, all_hit_probx, all_hit_proby, all_ladder, all_layer, all_mod, all_npix, all_nsimhit, all_nxpix, all_nypix, all_panel, all_pidhit, all_pixadc, all_pixcol, all_pixel_clst_err_x, all_pixel_clst_err_y, all_pixel_split, all_pixgx, all_pixgy, all_pixgz, all_pixrow, all_pixx, all_pixy, all_plaq, all_rechiterrx, all_rechiterry, all_rechitpullx, all_rechitpully, all_rechitresx, all_rechitresy, all_rechitx, all_rechity, all_rechitz, all_row1, all_row2, all_rows, all_side, all_simhitx, all_simhity, all_simphi, all_simproc, all_simpx, all_simpy, all_simpz, all_simtheta, all_simtrketa, all_simtrkphi, all_subdetid, all_trkid, all_vtxr, all_vtxz, all_x1, all_x2, all_y1, all_y2, all_z1, all_z2, alpha, TrackerHitAssociator::associateHit(), edmNew::DetSet< T >::begin(), beta, bigx, bigy, blade, charge, checkType_, SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), SiPixelRecHit::cluster(), SiPixelRecHit::clusterProbability(), PixelTopology::containsBigPixelInX(), PixelTopology::containsBigPixelInY(), gather_cfg::cout, detector_type, cond::rpcobgas::detid, disk, edgex, edgey, edmNew::DetSet< T >::end(), relativeConstraints::error, eta, edm::EventID::event(), evt, flipped, genType_, edm::EventSetup::get(), edm::Event::getByToken(), half, TrackingRecHit::hit(), hit_cprob0, hit_cprob1, hit_cprob2, hit_probx, hit_proby, i, edm::EventBase::id(), include_trk_hits_, PixelTopology::isItEdgePixelInX(), PixelTopology::isItEdgePixelInY(), ladder, layer, PSimHit::localPosition(), Topology::localPosition(), visualization-live-secondInstance_cfg::m, PV3DBase< T, PVType, FrameType >::mag(), HLT_25ns14e33_v1_cff::magneticField, math_pi, Topology::measurementPosition(), mod, SiStripDetId::moduleGeometry(), LocalTrajectoryParameters::momentum(), PixelTopology::ncolumns(), npix, PixelTopology::nrows(), nsimhit, NULL, nxpix, nypix, panel, PSimHit::particleType(), PV3DBase< T, PVType, FrameType >::perp(), phi, sysUtil::pid, pidhit, pixel_clst_err_x, pixel_clst_err_y, pixel_split, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, plaq, GloballyPositioned< T >::position(), position, PSimHit::processType(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), TrackerTopology::pxfBlade(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), TrackerTopology::pxfSide(), rechiterrx, rechiterry, rechitpullx, rechitpully, rechitresx, rechitresy, rechitx, rechity, rechitz, edm::EventID::run(), run, side, simhitx, simhity, simproc, findQualityFiles::size, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), StripSubdetector::stereo(), strip_charge, strip_clst_err_x, strip_clst_err_y, strip_cotalpha, strip_cotbeta, strip_edge, strip_hit_type, strip_is_stereo, strip_locbx, strip_locby, strip_locbz, strip_nprm, strip_nsimhit, strip_pidhit, strip_pidhit1, strip_pidhit2, strip_pidhit3, strip_pidhit4, strip_pidhit5, strip_prob, strip_qbin, strip_rechiterrx, strip_rechiterry, strip_rechitpullx, strip_rechitpully, strip_rechitresx, strip_rechitresx2, strip_rechitresy, strip_rechitx, strip_rechity, strip_rechitz, strip_simproc, strip_simproc1, strip_simproc2, strip_simproc3, strip_simproc4, strip_simproc5, strip_size, strip_split, strip_subdet_id, strip_tib_is_double_side, strip_tib_is_external_string, strip_tib_is_internal_string, strip_tib_is_rphi, strip_tib_is_stereo, strip_tib_is_z_minus_side, strip_tib_is_z_plus_side, strip_tib_layer, strip_tib_layer_number, strip_tib_module, strip_tib_module_number, strip_tib_order, strip_tib_side, strip_tib_string_number, strip_tob_is_double_side, strip_tob_is_rphi, strip_tob_is_stereo, strip_tob_is_z_minus_side, strip_tob_is_z_plus_side, strip_tob_layer, strip_tob_layer_number, strip_tob_module, strip_tob_module_number, strip_tob_rod_number, strip_tob_side, strip_trk_pt, DetId::subdetId(), subdetId, GeomDet::surface(), StripSubdetector::TIB, TrackerTopology::tibIsDoubleSide(), TrackerTopology::tibIsExternalString(), TrackerTopology::tibIsInternalString(), TrackerTopology::tibIsRPhi(), TrackerTopology::tibIsStereo(), TrackerTopology::tibIsZMinusSide(), TrackerTopology::tibIsZPlusSide(), TrackerTopology::tibLayer(), TrackerTopology::tibModule(), TrackerTopology::tibOrder(), TrackerTopology::tibSide(), TrackerTopology::tibString(), StripSubdetector::TOB, TrackerTopology::tobIsDoubleSide(), TrackerTopology::tobIsRPhi(), TrackerTopology::tobIsStereo(), TrackerTopology::tobIsZMinusSide(), TrackerTopology::tobIsZPlusSide(), TrackerTopology::tobLayer(), TrackerTopology::tobModule(), TrackerTopology::tobRod(), TrackerTopology::tobSide(), Surface::toGlobal(), toLocal(), tPixRecHitCollection, HLT_25ns14e33_v3_cff::trackCollection, patCandidatesForDimuonsSequences_cff::tracker, trackerHitAssociatorConfig_, testEve_cfg::tracks, trk_alpha, trk_beta, tSimTrackContainer, tTrackCollection, tTrajectory, ttree_all_hits_, ttree_track_hits_, ttree_track_hits_strip_, x, SiPixelCluster::Pixel::x, PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), SiPixelCluster::Pixel::y, PV2DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

{
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  es.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();

  using namespace edm;
  
  run = e.id().run();
  evt = e.id().event();
  
  if ( evt%1000 == 0 ) 
    cout << "evt = " << evt << endl;

  float math_pi = 3.14159265;
  float radtodeg = 180.0 / math_pi;
    
  LocalPoint position;
  LocalError error;
  float mindist = 999999.9;

  std::vector<PSimHit> matched;
  TrackerHitAssociator associate(e, trackerHitAssociatorConfig_);

  edm::ESHandle<TrackerGeometry> pDD;
  es.get<TrackerDigiGeometryRecord> ().get (pDD);
  const TrackerGeometry* tracker = &(* pDD);
  

  //cout << "...1..." << endl;


  edm::ESHandle<MagneticField> magneticField;
  es.get<IdealMagneticFieldRecord>().get( magneticField );
  //const MagneticField* magField_ = magFieldHandle.product();

  edm::FileInPath FileInPath_;
    


  // Strip hits ==============================================================================================================


  edm::Handle<vector<Trajectory> > trajCollectionHandle;
  
  e.getByToken( tTrajectory, trajCollectionHandle);
  //e.getByLabel( "generalTracks", trajCollectionHandle);

  for ( vector<Trajectory>::const_iterator it = trajCollectionHandle->begin(); it!=trajCollectionHandle->end(); ++it )
    {
      
      vector<TrajectoryMeasurement> tmColl = it->measurements();
      for ( vector<TrajectoryMeasurement>::const_iterator itTraj = tmColl.begin(); itTraj!=tmColl.end(); ++itTraj )
    {
      
      if ( !itTraj->updatedState().isValid() ) 
        continue;
       
      strip_rechitx     = -9999999.9;
      strip_rechity     = -9999999.9;
      strip_rechitz     = -9999999.9;
      strip_rechiterrx  = -9999999.9;
      strip_rechiterry  = -9999999.9;
      strip_rechitresx  = -9999999.9;
      strip_rechitresx2  = -9999999.9;
      

      strip_rechitresy  = -9999999.9;
      strip_rechitpullx = -9999999.9;
      strip_rechitpully = -9999999.9;
      strip_is_stereo   = -9999999  ;
      strip_hit_type    = -9999999  ;
      detector_type     = -9999999  ;
      
      strip_trk_pt    = -9999999.9;
      strip_cotalpha  = -9999999.9;
      strip_cotbeta   = -9999999.9;
      strip_locbx     = -9999999.9;
      strip_locby     = -9999999.9;
      strip_locbz     = -9999999.9;
      strip_charge    = -9999999.9;
      strip_size        = -9999999  ;

      strip_edge        = -9999999  ;
      strip_nsimhit     = -9999999  ;
      strip_pidhit      = -9999999  ;
      strip_simproc     = -9999999  ;


      strip_subdet_id = -9999999;
 
      strip_tib_layer             = -9999999;
      strip_tib_module            = -9999999;
      strip_tib_order             = -9999999;
      strip_tib_side              = -9999999;
      strip_tib_is_double_side    = -9999999;
      strip_tib_is_z_plus_side    = -9999999;
      strip_tib_is_z_minus_side   = -9999999;
      strip_tib_layer_number      = -9999999;
      strip_tib_string_number     = -9999999;
      strip_tib_module_number     = -9999999;
      strip_tib_is_internal_string= -9999999;
      strip_tib_is_external_string= -9999999;
      strip_tib_is_rphi           = -9999999;
      strip_tib_is_stereo         = -9999999;          
      
      strip_tob_layer             = -9999999;
      strip_tob_module            = -9999999;
      strip_tob_side              = -9999999;
      strip_tob_is_double_side    = -9999999;
      strip_tob_is_z_plus_side    = -9999999;
      strip_tob_is_z_minus_side   = -9999999;
      strip_tob_layer_number      = -9999999;
      strip_tob_rod_number        = -9999999;
      strip_tob_module_number     = -9999999;
          
      strip_tob_is_rphi           = -9999999;
      strip_tob_is_stereo         = -9999999;     

      strip_prob                  = -9999999.9;
      strip_qbin                  = -9999999;

      strip_nprm                  = -9999999;
      
      strip_pidhit1      = -9999999  ;
      strip_simproc1     = -9999999  ;
      
      strip_pidhit2      = -9999999  ;
      strip_simproc2     = -9999999  ;
      
      strip_pidhit3      = -9999999  ;
      strip_simproc3     = -9999999  ;
      
      strip_pidhit4      = -9999999  ;
      strip_simproc4     = -9999999  ;

      strip_pidhit5      = -9999999  ;
      strip_simproc5     = -9999999  ;

      strip_split        = -9999999;   
      strip_clst_err_x   = -9999999.9;
      strip_clst_err_y   = -9999999.9;

      const TransientTrackingRecHit::ConstRecHitPointer trans_trk_rec_hit_point = itTraj->recHit();
     
      if ( trans_trk_rec_hit_point == NULL )
        continue;

      const TrackingRecHit *trk_rec_hit = (*trans_trk_rec_hit_point).hit();

      if ( trk_rec_hit == NULL )
        continue;

      DetId detid = (trk_rec_hit)->geographicalId();

      strip_subdet_id = (int)detid.subdetId();
      
      if ( (int)detid.subdetId() != (int)(StripSubdetector::TIB) && (int)detid.subdetId() != (int)(StripSubdetector::TOB) )
        continue;

      const SiStripMatchedRecHit2D* matchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>( (*trans_trk_rec_hit_point).hit() );
      const SiStripRecHit2D       * hit2d      = dynamic_cast<const SiStripRecHit2D       *>( (*trans_trk_rec_hit_point).hit() );
      const SiStripRecHit1D       * hit1d      = dynamic_cast<const SiStripRecHit1D       *>( (*trans_trk_rec_hit_point).hit() );
      
      if ( !matchedhit && !hit2d && !hit1d )
        continue;
      
      position = (trk_rec_hit)->localPosition(); 
      error = (trk_rec_hit)->localPositionError();
      
      strip_rechitx = position.x();
      strip_rechity = position.y();
      strip_rechitz = position.z();
      strip_rechiterrx = sqrt( error.xx() );
      strip_rechiterry = sqrt( error.yy() );
      
      
      //cout << "strip_rechitx = " << strip_rechitx << endl;          
      //cout << "strip_rechity = " << strip_rechity << endl;
      //cout << "strip_rechitz = " << strip_rechitz << endl;
      
      //cout << "strip_rechiterrx = " << strip_rechiterrx << endl;
      //cout << "strip_rechiterry = " << strip_rechiterry << endl;
      
      TrajectoryStateOnSurface tsos = itTraj->updatedState(); 
      
      strip_trk_pt = tsos.globalMomentum().perp();
      
      LocalTrajectoryParameters ltp = tsos.localParameters();
      
      LocalVector localDir = ltp.momentum()/ltp.momentum().mag();
      
      float locx = localDir.x();
      float locy = localDir.y();
      float locz = localDir.z();
      
      //alpha_ = atan2( locz, locx );
      //beta_  = atan2( locz, locy );
      
      strip_cotalpha = locx/locz;
      strip_cotbeta  = locy/locz;
          

      StripSubdetector StripSubdet = (StripSubdetector)detid;

      if ( StripSubdet.stereo() == 0 )
        strip_is_stereo = 0;
      else
        strip_is_stereo = 1;


      SiStripDetId si_strip_det_id = SiStripDetId( detid );
      
      //const StripGeomDetUnit* strip_geom_det_unit = dynamic_cast<const StripGeomDetUnit*> ( tracker->idToDet( detid ) );
      const StripGeomDetUnit* strip_geom_det_unit = (const StripGeomDetUnit*)tracker->idToDetUnit( detid );
      
      if ( strip_geom_det_unit != NULL )
        {
          LocalVector lbfield 
        = (strip_geom_det_unit->surface()).toLocal( (*magneticField).inTesla( strip_geom_det_unit->surface().position() ) ); 
          
          strip_locbx = lbfield.x();
          strip_locby = lbfield.y();
          strip_locbz = lbfield.z();
        }
      

      //enum ModuleGeometry {UNKNOWNGEOMETRY, IB1, IB2, OB1, OB2, W1A, W2A, W3A, W1B, W2B, W3B, W4, W5, W6, W7};
      
      if ( si_strip_det_id.moduleGeometry() == 1 )
        {
          detector_type = 1;
          //cout << "si_strip_det_id.moduleGeometry() = IB1" << endl;
          //cout << "si_strip_det_id.moduleGeometry() = " << si_strip_det_id.moduleGeometry() << endl;
        }
      else if ( si_strip_det_id.moduleGeometry() == 2 )
        {
          detector_type = 2;
          //cout << "si_strip_det_id.moduleGeometry() = IB2" << endl;
          //cout << "si_strip_det_id.moduleGeometry() = " << si_strip_det_id.moduleGeometry() << endl;
        }
      else if ( si_strip_det_id.moduleGeometry() == 3 )
        {
          detector_type = 3;
          //cout << "si_strip_det_id.moduleGeometry() = OB1" << endl;
          //cout << "si_strip_det_id.moduleGeometry() = " << si_strip_det_id.moduleGeometry() << endl;
        }
      else if ( si_strip_det_id.moduleGeometry() == 4 )
        {
          detector_type = 4;
          //cout << "si_strip_det_id.moduleGeometry() = OB2" << endl;
          //cout << "si_strip_det_id.moduleGeometry() = " << si_strip_det_id.moduleGeometry() << endl;
        }
      

      // Store ntuple variables 
      
      if ( (int)detid.subdetId() == int(StripSubdetector::TIB) )
        {
          
          
          
          strip_tib_layer              = (int)tTopo->tibLayer( detid );
          strip_tib_module             = (int)tTopo->tibModule( detid ); 
          strip_tib_order              = (int)tTopo->tibOrder( detid );
          strip_tib_side               = (int)tTopo->tibSide( detid );
          strip_tib_is_double_side     = (int)tTopo->tibIsDoubleSide( detid );
          strip_tib_is_z_plus_side     = (int)tTopo->tibIsZPlusSide( detid );
          strip_tib_is_z_minus_side    = (int)tTopo->tibIsZMinusSide( detid );
          strip_tib_layer_number       = (int)tTopo->tibLayer( detid );
          strip_tib_string_number      = (int)tTopo->tibString( detid ) ;
          strip_tib_module_number      = (int)tTopo->tibModule( detid );
          strip_tib_is_internal_string = (int)tTopo->tibIsInternalString( detid );
          strip_tib_is_external_string = (int)tTopo->tibIsExternalString( detid );
          strip_tib_is_rphi            = (int)tTopo->tibIsRPhi( detid );
          strip_tib_is_stereo          = (int)tTopo->tibIsStereo( detid );
        }
      
      
      if ( (int)detid.subdetId() == int(StripSubdetector::TOB) )
        {
          
          
          
          strip_tob_layer              = (int)tTopo->tobLayer( detid );
          strip_tob_module             = (int)tTopo->tobModule( detid );
          
          strip_tob_side               = (int)tTopo->tobSide( detid );
          strip_tob_is_double_side     = (int)tTopo->tobIsDoubleSide( detid );
          strip_tob_is_z_plus_side     = (int)tTopo->tobIsZPlusSide( detid );
          strip_tob_is_z_minus_side    = (int)tTopo->tobIsZMinusSide( detid );
          strip_tob_layer_number       = (int)tTopo->tobLayer( detid );
          strip_tob_rod_number         = (int)tTopo->tobRod( detid );
          strip_tob_module_number      = (int)tTopo->tobModule( detid );
          
          
          strip_tob_is_rphi            = (int)tTopo->tobIsRPhi( detid );
          strip_tob_is_stereo          = (int)tTopo->tobIsStereo( detid );
              
        }
     
     
      if ( matchedhit ) 
        {
          //cout << endl << endl << endl;
          //cout << "Found a SiStripMatchedRecHit2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! " << endl;
          //cout << endl << endl << endl;
          strip_hit_type = 0;

        } //  if ( matchedhit )

   
      if ( hit1d )
        {
          strip_hit_type = 1;

          const SiStripRecHit1D::ClusterRef & cluster = hit1d->cluster();

          if ( cluster->getSplitClusterError()  > 0.0 )
        strip_split = 1;
          else 
        strip_split = 0;
          
          strip_clst_err_x = cluster->getSplitClusterError();
          //strip_clst_err_y = ... 

          // Get cluster total charge
          const auto & stripCharges = cluster->amplitudes();
          uint16_t charge = 0;
          for (unsigned int i = 0; i < stripCharges.size(); ++i) 
        {
          charge += stripCharges[i];
        }
          
          strip_charge = (float)charge;
          strip_size = (int)( (cluster->amplitudes()).size() );


          // Association of the rechit to the simhit
          float mindist = 999999.9;
          matched.clear();  
          
          matched = associate.associateHit(*hit1d);

          strip_nsimhit = (int)matched.size();
          
          if ( !matched.empty()) 
        {
          PSimHit closest;
          
          // Get the closest simhit
          
          int strip_nprimaries = 0; 
          int current_index = 0;
         
          for ( vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); ++m)
            {
              ++current_index;
              
              if ( (*m).processType() == 2 )
            ++strip_nprimaries;

              if ( current_index == 1 )
            {
              strip_pidhit1 = (*m).particleType();
              strip_simproc1 = (*m).processType();
            }
              else if ( current_index == 2 )
            {
              strip_pidhit2 = (*m).particleType();
              strip_simproc2 = (*m).processType();
            }
              else if ( current_index == 3 )
            {
              strip_pidhit3 = (*m).particleType();
              strip_simproc3 = (*m).processType();
            }
              else if ( current_index == 4 )
            {
              strip_pidhit4 = (*m).particleType();
              strip_simproc4 = (*m).processType();
            }
              else if ( current_index == 5 )
            {
              strip_pidhit5 = (*m).particleType();
              strip_simproc5 = (*m).processType();
            }


              float dist = abs( (hit1d)->localPosition().x() - (*m).localPosition().x() );
              
              if ( dist<mindist )
            {
              mindist = dist;
              closest = (*m);
            }
            
            } // for ( vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); ++m)
          
          strip_nprm = strip_nprimaries;

          strip_rechitresx = strip_rechitx - closest.localPosition().x();
          strip_rechitresy = strip_rechity - closest.localPosition().y();
 
          strip_rechitresx2 = strip_rechitx - matched[0].localPosition().x();

          strip_rechitpullx = strip_rechitresx / strip_rechiterrx;
          strip_rechitpully = strip_rechitresy / strip_rechiterry;

          strip_pidhit = closest.particleType();
          strip_simproc = closest.processType();

          
        } //   if( !matched.empty()) 
          

          //strip_prob = (hit1d)->getTemplProb();
          //strip_qbin = (hit1d)->getTemplQbin();

          //cout << endl;
          //cout << "SiPixelErrorEstimation 1d hit: " << endl;
          //cout << "prob 1d = " << strip_prob << endl;
          //cout << "qbin 1d = " << strip_qbin << endl;
          //cout << endl;

          
          // Is the cluster on edge ?
          /*
        SiStripDetInfoFileReader* reader;
        
        FileInPath_("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat");
        
        reader = new SiStripDetInfoFileReader( FileInPath_.fullPath() );
        
        uint16_t firstStrip = cluster->firstStrip();
        uint16_t lastStrip = firstStrip + (cluster->amplitudes()).size() -1;
        unsigned short Nstrips;
        Nstrips = reader->getNumberOfApvsAndStripLength(id1).first*128;
        
        if ( firstStrip == 0 || lastStrip == (Nstrips-1) ) 
        strip_edge = 1;
        else
        strip_edge = 0;
          */

        } // if ( hit1d )


      if ( hit2d )
        {
          strip_hit_type = 2;
          
          const SiStripRecHit1D::ClusterRef & cluster = hit2d->cluster();
          
          //if ( cluster->getSplitClusterError()  > 0.0 )
          //strip_split = 1;
          //else 
          //strip_split = 0;

          //strip_clst_err_x = cluster->getSplitClusterError();
          //strip_clst_err_y = ... 

          // Get cluster total charge
          auto & stripCharges = cluster->amplitudes();
          uint16_t charge = 0;
          for (unsigned int i = 0; i < stripCharges.size(); ++i) 
        {
          charge += stripCharges[i];
        }
          
          strip_charge = (float)charge;
          strip_size = (int)( (cluster->amplitudes()).size() );
         
          // Association of the rechit to the simhit
          float mindist = 999999.9;
          matched.clear();  
          
          matched = associate.associateHit(*hit2d);

          strip_nsimhit = (int)matched.size();
          
          if ( !matched.empty()) 
        {
          PSimHit closest;
          
          // Get the closest simhit
          
          for ( vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); ++m)
            {
              float dist = abs( (hit2d)->localPosition().x() - (*m).localPosition().x() );
              
              if ( dist<mindist )
            {
              mindist = dist;
              closest = (*m);
            }
            }
          
          strip_rechitresx = strip_rechitx - closest.localPosition().x();
          strip_rechitresy = strip_rechity - closest.localPosition().y();
 
          strip_rechitpullx = strip_rechitresx / strip_rechiterrx;
          strip_rechitpully = strip_rechitresy / strip_rechiterry;

          strip_pidhit = closest.particleType();
          strip_simproc = closest.processType();

          
        } //   if( !matched.empty()) 
          

          //strip_prob = (hit2d)->getTemplProb();
          //strip_qbin = (hit2d)->getTemplQbin();

          //cout << endl;
          //cout << "SiPixelErrorEstimation 2d hit: " << endl;
          //cout << "prob 2d = " << strip_prob << endl;
          //cout << "qbin 2d = " << strip_qbin << endl;
          //cout << endl;
          
          
        } // if ( hit2d )


      ttree_track_hits_strip_->Fill();


    } // for ( vector<TrajectoryMeasurement>::const_iterator itTraj = tmColl.begin(); itTraj!=tmColl.end(); ++itTraj )

    } // for( vector<Trajectory>::const_iterator it = trajCollectionHandle->begin(); it!=trajCollectionHandle->end(); ++it)








  //cout << "...2..." << endl;






  // --------------------------------------- all hits -----------------------------------------------------------
  edm::Handle<SiPixelRecHitCollection> recHitColl;
  e.getByToken(tPixRecHitCollection, recHitColl);

  Handle<edm::SimTrackContainer> simtracks;
  e.getByToken(tSimTrackContainer, simtracks);

  //-----Iterate over detunits
  for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++) 
    {
      DetId detId = ((*it)->geographicalId());
      
      SiPixelRecHitCollection::const_iterator dsmatch = recHitColl->find(detId);
      if (dsmatch == recHitColl->end()) continue;

      SiPixelRecHitCollection::DetSet pixelrechitRange = *dsmatch;
      SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
      SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
      SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
      std::vector<PSimHit> matched;
      
      //----Loop over rechits for this detId
      for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter) 
    {
      matched.clear();
      matched = associate.associateHit(*pixeliter);
      // only consider rechits that have associated simhit
      // otherwise cannot determine residiual
      if ( matched.empty() )
        {
          cout << "SiPixelErrorEstimation::analyze: rechits without associated simhit !!!!!!!" 
           << endl;
          continue;
        }
        
      all_subdetid = -9999;
      
      all_layer = -9999;
      all_ladder = -9999;
      all_mod = -9999;
      
      all_side = -9999;
      all_disk = -9999;
      all_blade = -9999;
      all_panel = -9999;
      all_plaq = -9999;
      
      all_half = -9999;
      all_flipped = -9999;
      
      all_cols = -9999;
      all_rows = -9999;
      
      all_rechitx = -9999;
      all_rechity = -9999;
      all_rechitz = -9999;
      
      all_simhitx = -9999;
      all_simhity = -9999;

      all_rechiterrx = -9999;
      all_rechiterry = -9999;
      
      all_rechitresx = -9999;
      all_rechitresy = -9999;
      
      all_rechitpullx = -9999;
      all_rechitpully = -9999;
      
      all_npix = -9999;
      all_nxpix = -9999;
      all_nypix = -9999;
          
      all_edgex = -9999;
      all_edgey = -9999;
      
      all_bigx = -9999;
      all_bigy = -9999;
      
      all_alpha = -9999;
      all_beta = -9999;
      
      all_simphi = -9999;
      all_simtheta = -9999;
      
      all_simhitx = -9999;
      all_simhity = -9999;
      
      all_nsimhit = -9999;
      all_pidhit = -9999;
      all_simproc = -9999;
      
      all_vtxr = -9999;
      all_vtxz = -9999;
      
      all_simpx = -9999;
      all_simpy = -9999;
      all_simpz = -9999;
      
      all_eloss = -9999;
          
      all_trkid = -9999;
      
      all_x1 = -9999;
      all_x2 = -9999;
      all_y1 = -9999;
      all_y2 = -9999;
      all_z1 = -9999;
      all_z2 = -9999;
      
      all_row1 = -9999;
      all_row2 = -9999;
      all_col1 = -9999;
      all_col2 = -9999;
      
      all_gx1 = -9999;
      all_gx2 = -9999;
      all_gy1 = -9999;
      all_gy2 = -9999;
      all_gz1 = -9999;
      all_gz2 = -9999;
      
      all_simtrketa = -9999;
      all_simtrkphi = -9999;
      
      all_clust_row = -9999;
      all_clust_col = -9999;
      
      all_clust_x = -9999;
      all_clust_y = -9999;
      
      all_clust_q = -9999;
      
      all_clust_maxpixcol = -9999;
      all_clust_maxpixrow = -9999;
      all_clust_minpixcol = -9999;
      all_clust_minpixrow = -9999;
      
      all_clust_geoid = -9999;
      
      all_clust_alpha = -9999;
      all_clust_beta = -9999;
      
      /*
        for (int i=0; i<all_npix; ++i)
        {
        all_pixrow[i] = -9999;
        all_pixcol[i] = -9999;
        all_pixadc[i] = -9999;
        all_pixx[i] = -9999;
        all_pixy[i] = -9999;
        all_pixgx[i] = -9999;
        all_pixgy[i] = -9999;
        all_pixgz[i] = -9999;
        }
      */

      all_hit_probx = -9999;
      all_hit_proby = -9999;
      all_hit_cprob0 = -9999;
      all_hit_cprob1 = -9999;
      all_hit_cprob2 = -9999;

      all_pixel_split = -9999;
      all_pixel_clst_err_x = -9999.9;
      all_pixel_clst_err_y = -9999.9;

      all_nsimhit = (int)matched.size();
      
      all_subdetid = (int)detId.subdetId();
      // only consider rechits in pixel barrel and pixel forward 
      if ( !(all_subdetid==1 || all_subdetid==2) )
        {
          cout << "SiPixelErrorEstimation::analyze: Not in a pixel detector !!!!!" << endl; 
          continue;
        }

      const PixelGeomDetUnit* theGeomDet 
        = dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(detId) );
      
      const PixelTopology* topol = &(theGeomDet->specificTopology());

      if ( pixeliter->cluster()->getSplitClusterErrorX()  > 0.0 && 
           pixeliter->cluster()->getSplitClusterErrorY()  > 0.0 )
        {
          all_pixel_split = 1;
        }
      else
        {
          all_pixel_split = 0;
        }
    
      all_pixel_clst_err_x = pixeliter->cluster()->getSplitClusterErrorX();
      all_pixel_clst_err_y = pixeliter->cluster()->getSplitClusterErrorY();


      const int maxPixelCol = pixeliter->cluster()->maxPixelCol();
      const int maxPixelRow = pixeliter->cluster()->maxPixelRow();
      const int minPixelCol = pixeliter->cluster()->minPixelCol();
      const int minPixelRow = pixeliter->cluster()->minPixelRow();
      
      //all_hit_probx  = (float)pixeliter->probabilityX();
      //all_hit_proby  = (float)pixeliter->probabilityY();
      all_hit_cprob0 = (float)pixeliter->clusterProbability(0);
      all_hit_cprob1 = (float)pixeliter->clusterProbability(1);
      all_hit_cprob2 = (float)pixeliter->clusterProbability(2);
      
      // check whether the cluster is at the module edge 
      if ( topol->isItEdgePixelInX( minPixelRow ) || 
           topol->isItEdgePixelInX( maxPixelRow ) )
        all_edgex = 1;
      else 
        all_edgex = 0;
      
      if ( topol->isItEdgePixelInY( minPixelCol ) || 
           topol->isItEdgePixelInY( maxPixelCol ) )
        all_edgey = 1;
      else 
        all_edgey = 0;
      
      // check whether this rechit contains big pixels
      if ( topol->containsBigPixelInX(minPixelRow, maxPixelRow) )
        all_bigx = 1;
      else 
        all_bigx = 0;
      
      if ( topol->containsBigPixelInY(minPixelCol, maxPixelCol) )
        all_bigy = 1;
      else 
        all_bigy = 0;
      
      if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel ) 
        {
          
          all_layer = tTopo->pxbLayer(detId);
          all_ladder = tTopo->pxbLadder(detId);
          all_mod = tTopo->pxbModule(detId);
          
          int tmp_nrows = theGeomDet->specificTopology().nrows();
          if ( tmp_nrows == 80 ) 
        all_half = 1;
          else if ( tmp_nrows == 160 ) 
        all_half = 0;
          else 
        cout << "-------------------------------------------------- Wrong module size !!!" << endl;
          
          float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
          float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
          
          if ( tmp2<tmp1 ) 
        all_flipped = 1;
          else 
        all_flipped = 0;
        }
      else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
        {
          
          all_side  = tTopo->pxfSide(detId);
          all_disk  = tTopo->pxfDisk(detId);
          all_blade = tTopo->pxfBlade(detId);
          all_panel = tTopo->pxfPanel(detId);
          all_plaq  = tTopo->pxfModule(detId); // also known as plaquette
          
        } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
      else std::cout << "We are not in the pixel detector" << (int)detId.subdetId() << endl;
      
      all_cols = theGeomDet->specificTopology().ncolumns();
      all_rows = theGeomDet->specificTopology().nrows();
          
      LocalPoint lp = pixeliter->localPosition();
      // gavril: change this name
      all_rechitx = lp.x();
      all_rechity = lp.y();
      all_rechitz = lp.z();
      
      LocalError le = pixeliter->localPositionError();
      all_rechiterrx = sqrt( le.xx() );
      all_rechiterry = sqrt( le.yy() );

      bool found_hit_from_generated_particle = false;
      
      //---Loop over sim hits, fill closest
      float closest_dist = 99999.9;
      std::vector<PSimHit>::const_iterator closest_simhit = matched.begin();
      
      for (std::vector<PSimHit>::const_iterator m = matched.begin(); m < matched.end(); m++) 
        {
          if ( checkType_ )
        {
          int pid = (*m).particleType();
          if ( abs(pid) != genType_ )
            continue;
        } 
          
          float simhitx = 0.5 * ( (*m).entryPoint().x() + (*m).exitPoint().x() );
          float simhity = 0.5 * ( (*m).entryPoint().y() + (*m).exitPoint().y() );
          
          float x_res = simhitx - rechitx;
          float y_res = simhity - rechity;
          
          float dist = sqrt( x_res*x_res + y_res*y_res );         
          
          if ( dist < closest_dist ) 
        {
          closest_dist = dist;
          closest_simhit = m;
          found_hit_from_generated_particle = true;
        } 
        } // end sim hit loop
      
      // If this recHit does not have any simHit with the same particleType as the particles generated
      // ignore it as most probably comes from delta rays.
      if ( checkType_ && !found_hit_from_generated_particle )
        continue; 

      all_x1 = (*closest_simhit).entryPoint().x(); // width (row index, in col direction)
      all_y1 = (*closest_simhit).entryPoint().y(); // length (col index, in row direction) 
      all_z1 = (*closest_simhit).entryPoint().z(); 
      all_x2 = (*closest_simhit).exitPoint().x();
      all_y2 = (*closest_simhit).exitPoint().y();
      all_z2 = (*closest_simhit).exitPoint().z();
      GlobalPoint GP1 = 
        theGeomDet->surface().toGlobal( Local3DPoint( (*closest_simhit).entryPoint().x(),
                              (*closest_simhit).entryPoint().y(),
                              (*closest_simhit).entryPoint().z() ) );
      GlobalPoint GP2 = 
        theGeomDet->surface().toGlobal (Local3DPoint( (*closest_simhit).exitPoint().x(),
                              (*closest_simhit).exitPoint().y(),
                              (*closest_simhit).exitPoint().z() ) );
      all_gx1 = GP1.x();
      all_gx2 = GP2.x();
      all_gy1 = GP1.y();
      all_gy2 = GP2.y();
      all_gz1 = GP1.z();
      all_gz2 = GP2.z();
      
      MeasurementPoint mp1 =
        topol->measurementPosition( LocalPoint( (*closest_simhit).entryPoint().x(),
                            (*closest_simhit).entryPoint().y(),
                            (*closest_simhit).entryPoint().z() ) );
      MeasurementPoint mp2 =
        topol->measurementPosition( LocalPoint( (*closest_simhit).exitPoint().x(),
                            (*closest_simhit).exitPoint().y(), 
                            (*closest_simhit).exitPoint().z() ) );
      all_row1 = mp1.x();
      all_col1 = mp1.y();
      all_row2 = mp2.x();
      all_col2 = mp2.y();
      
      all_simhitx = 0.5*(all_x1+all_x2);  
      all_simhity = 0.5*(all_y1+all_y2);  
      
      all_rechitresx = all_rechitx - all_simhitx;
      all_rechitresy = all_rechity - all_simhity;

      all_rechitpullx = all_rechitresx / all_rechiterrx;
      all_rechitpully = all_rechitresy / all_rechiterry;
      
      SiPixelRecHit::ClusterRef const& clust = pixeliter->cluster();
      
      all_npix = clust->size();
      all_nxpix = clust->sizeX();
      all_nypix = clust->sizeY();

      all_clust_row = clust->x();
      all_clust_col = clust->y();
      
      LocalPoint lp2 = topol->localPosition( MeasurementPoint( all_clust_row, all_clust_col ) );
      all_clust_x = lp2.x();
      all_clust_y = lp2.y();

      all_clust_q = clust->charge();

      all_clust_maxpixcol = clust->maxPixelCol();
      all_clust_maxpixrow = clust->maxPixelRow();
      all_clust_minpixcol = clust->minPixelCol();
      all_clust_minpixrow = clust->minPixelRow();
      
      all_clust_geoid = 0;  // never set!
  
      all_simpx  = (*closest_simhit).momentumAtEntry().x();
      all_simpy  = (*closest_simhit).momentumAtEntry().y();
      all_simpz  = (*closest_simhit).momentumAtEntry().z();
      all_eloss = (*closest_simhit).energyLoss();
      all_simphi   = (*closest_simhit).phiAtEntry();
      all_simtheta = (*closest_simhit).thetaAtEntry();
      all_pidhit = (*closest_simhit).particleType();
      all_trkid = (*closest_simhit).trackId();
      
      //--- Fill alpha and beta -- more useful for exploring the residuals...
      all_beta  = atan2(all_simpz, all_simpy);
      all_alpha = atan2(all_simpz, all_simpx);
      
      all_simproc = (int)closest_simhit->processType();
      
      const edm::SimTrackContainer& trks = *(simtracks.product());
      SimTrackContainer::const_iterator trksiter;
      for (trksiter = trks.begin(); trksiter != trks.end(); trksiter++) 
        if ( (int)trksiter->trackId() == all_trkid ) 
          {
        all_simtrketa = trksiter->momentum().eta();
        all_simtrkphi = trksiter->momentum().phi();
          }
      
      all_vtxz = theGeomDet->surface().position().z();
      all_vtxr = theGeomDet->surface().position().perp();
      
      //computeAnglesFromDetPosition(clust, 
      //               theGeomDet, 
      //               all_clust_alpha, all_clust_beta )

      const std::vector<SiPixelCluster::Pixel>& pixvector = clust->pixels();
      for ( int i=0;  i<(int)pixvector.size(); ++i)
        {
          SiPixelCluster::Pixel holdpix = pixvector[i];
          all_pixrow[i] = holdpix.x;
          all_pixcol[i] = holdpix.y;
          all_pixadc[i] = holdpix.adc;
          LocalPoint lp = topol->localPosition(MeasurementPoint(holdpix.x, holdpix.y));
          all_pixx[i] = lp.x();
          all_pixy[i]= lp.y();
          GlobalPoint GP =  theGeomDet->surface().toGlobal(Local3DPoint(lp.x(),lp.y(),lp.z()));
          all_pixgx[i] = GP.x();    
          all_pixgy[i]= GP.y();
          all_pixgz[i]= GP.z();
        }

      ttree_all_hits_->Fill();
      
    } // for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter)

    } // for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++) 
 
  // ------------------------------------------------ all hits ---------------------------------------------------------------


  //cout << "...3..." << endl;


  // ------------------------------------------------ track hits only -------------------------------------------------------- 

  
  
  if ( include_trk_hits_ )
    {
      // Get tracks
      edm::Handle<reco::TrackCollection> trackCollection;
      e.getByToken(tTrackCollection, trackCollection);
      const reco::TrackCollection *tracks = trackCollection.product();
      reco::TrackCollection::const_iterator tciter;
      
      if ( tracks->size() > 0 )
    {
      // Loop on tracks
      for ( tciter=tracks->begin(); tciter!=tracks->end(); ++tciter)
        {
          // First loop on hits: find matched hits
          for ( trackingRecHit_iterator it = tciter->recHitsBegin(); it != tciter->recHitsEnd(); ++it) 
        {
          const TrackingRecHit &trk_rec_hit = **it;
          // Is it a matched hit?
          const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(&trk_rec_hit);
          
          if ( matchedhit ) 
            {
              rechitx = -9999.9;
              rechity = -9999.9;
              rechitz = -9999.9;
              rechiterrx = -9999.9;
              rechiterry = -9999.9;           
              rechitresx = -9999.9;
              rechitresy = -9999.9;
              rechitpullx = -9999.9;
              rechitpully = -9999.9;
              
              npix = -9999;
              nxpix = -9999;
              nypix = -9999;
              charge = -9999.9;
              
              edgex = -9999;
              edgey = -9999;
              
              bigx = -9999;
              bigy = -9999;
              
              alpha = -9999.9;
              beta  = -9999.9;
              
              phi = -9999.9;
              eta = -9999.9;
              
              subdetId = -9999;
              
              layer  = -9999; 
              ladder = -9999; 
              mod    = -9999; 
              side   = -9999;  
              disk   = -9999;  
              blade  = -9999; 
              panel  = -9999; 
              plaq   = -9999; 
              
              half = -9999;
              flipped = -9999;
              
              nsimhit = -9999;
              pidhit  = -9999;
              simproc = -9999;
              
              simhitx = -9999.9;
              simhity = -9999.9;

              hit_probx = -9999.9;
              hit_proby = -9999.9;
              hit_cprob0 = -9999.9;
              hit_cprob1 = -9999.9;
              hit_cprob2 = -9999.9;
  
              pixel_split = -9999;

              pixel_clst_err_x = -9999.9;
              pixel_clst_err_y = -9999.9;
              
              position = (*it)->localPosition();
              error = (*it)->localPositionError();
              
              rechitx = position.x();
              rechity = position.y();
              rechitz = position.z();
              rechiterrx = sqrt(error.xx());
              rechiterry = sqrt(error.yy());
              
              npix = matchedhit->cluster()->size();
              nxpix = matchedhit->cluster()->sizeX();
              nypix = matchedhit->cluster()->sizeY();
              charge = matchedhit->cluster()->charge();

              if ( matchedhit->cluster()->getSplitClusterErrorX() > 0.0 && 
               matchedhit->cluster()->getSplitClusterErrorY() > 0.0 )
            pixel_split = 1;
              else
            pixel_split = 0;
              
              pixel_clst_err_x = matchedhit->cluster()->getSplitClusterErrorX();
              pixel_clst_err_y = matchedhit->cluster()->getSplitClusterErrorY();

              //hit_probx  = (float)matchedhit->probabilityX();
              //hit_proby  = (float)matchedhit->probabilityY();
              hit_cprob0 = (float)matchedhit->clusterProbability(0);
              hit_cprob1 = (float)matchedhit->clusterProbability(1);
              hit_cprob2 = (float)matchedhit->clusterProbability(2);
              
              
              //Association of the rechit to the simhit
              matched.clear();
              matched = associate.associateHit(*matchedhit);
              
              nsimhit = (int)matched.size();
              
              if ( !matched.empty() ) 
            {
              mindist = 999999.9;
              float distx, disty, dist;
              bool found_hit_from_generated_particle = false;
              
              int n_assoc_muon = 0;
              
              vector<PSimHit>::const_iterator closestit = matched.begin();
              for (vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); ++m)
                {
                  if ( checkType_ )
                { // only consider associated simhits with the generated pid (muons)
                  int pid = (*m).particleType();
                  if ( abs(pid) != genType_ )
                    continue;
                }
                  
                  float simhitx = 0.5 * ( (*m).entryPoint().x() + (*m).exitPoint().x() );
                  float simhity = 0.5 * ( (*m).entryPoint().y() + (*m).exitPoint().y() );
                  
                  distx = fabs(rechitx - simhitx);
                  disty = fabs(rechity - simhity);
                  dist = sqrt( distx*distx + disty*disty );
                  
                  if ( dist < mindist )
                {
                  n_assoc_muon++;
                  
                  mindist = dist;
                  closestit = m;
                  found_hit_from_generated_particle = true;
                }
                } // for (vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++)
              
              // This recHit does not have any simHit with the same particleType as the particles generated
              // Ignore it as most probably come from delta rays.
              if ( checkType_ && !found_hit_from_generated_particle )
                continue; 
              
              //if ( n_assoc_muon > 1 )
              //{
              //  // cout << " ----- This is not good: n_assoc_muon = " << n_assoc_muon << endl;
              //  // cout << "evt = " << evt << endl;
              //}
              
              DetId detId = (*it)->geographicalId();

              const PixelGeomDetUnit* theGeomDet =
                dynamic_cast<const PixelGeomDetUnit*> ((*tracker).idToDet(detId) );
              
              const PixelTopology* theTopol = &(theGeomDet->specificTopology());

              pidhit = (*closestit).particleType();
              simproc = (int)(*closestit).processType();
              
              simhitx = 0.5*( (*closestit).entryPoint().x() + (*closestit).exitPoint().x() );
              simhity = 0.5*( (*closestit).entryPoint().y() + (*closestit).exitPoint().y() );
              
              rechitresx = rechitx - simhitx;
              rechitresy = rechity - simhity;
              rechitpullx = ( rechitx - simhitx ) / sqrt(error.xx());
              rechitpully = ( rechity - simhity ) / sqrt(error.yy());
              
              float simhitpx = (*closestit).momentumAtEntry().x();
              float simhitpy = (*closestit).momentumAtEntry().y();
              float simhitpz = (*closestit).momentumAtEntry().z();
              
              beta = atan2(simhitpz, simhitpy) * radtodeg;
              alpha = atan2(simhitpz, simhitpx) * radtodeg;
              
              //beta  = fabs(atan2(simhitpz, simhitpy)) * radtodeg;
              //alpha = fabs(atan2(simhitpz, simhitpx)) * radtodeg;

              // calculate alpha and beta exactly as in PixelCPEBase.cc
              float locx = simhitpx;
              float locy = simhitpy;
              float locz = simhitpz;
              
              bool isFlipped = false;
              float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
              float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
              if ( tmp2<tmp1 ) 
                isFlipped = true;
              else 
                isFlipped = false;    

              trk_alpha = acos(locx/sqrt(locx*locx+locz*locz)) * radtodeg;
              if ( isFlipped )                    // &&& check for FPIX !!!
                trk_alpha = 180.0 - trk_alpha ;
              
              trk_beta = acos(locy/sqrt(locy*locy+locz*locz)) * radtodeg;
              

              phi = tciter->momentum().phi() / math_pi*180.0;
              eta = tciter->momentum().eta();
              
              const int maxPixelCol = (*matchedhit).cluster()->maxPixelCol();
              const int maxPixelRow = (*matchedhit).cluster()->maxPixelRow();
              const int minPixelCol = (*matchedhit).cluster()->minPixelCol();
              const int minPixelRow = (*matchedhit).cluster()->minPixelRow();
                      
              // check whether the cluster is at the module edge 
              if ( theTopol->isItEdgePixelInX( minPixelRow ) || 
                   theTopol->isItEdgePixelInX( maxPixelRow ) )
                edgex = 1;
              else 
                edgex = 0;
              
              if ( theTopol->isItEdgePixelInY( minPixelCol ) || 
                   theTopol->isItEdgePixelInY( maxPixelCol ) )
                edgey = 1;
              else 
                edgey = 0;
              
              // check whether this rechit contains big pixels
              if ( theTopol->containsBigPixelInX(minPixelRow, maxPixelRow) )
                bigx = 1;
              else 
                bigx = 0;
              
              if ( theTopol->containsBigPixelInY(minPixelCol, maxPixelCol) )
                bigy = 1;
              else 
                bigy = 0;
              
              subdetId = (int)detId.subdetId();
              
              if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel ) 
                { 
      
                  int tmp_nrows = theGeomDet->specificTopology().nrows();
                  if ( tmp_nrows == 80 ) 
                half = 1;
                  else if ( tmp_nrows == 160 ) 
                half = 0;
                  else 
                cout << "-------------------------------------------------- Wrong module size !!!" << endl;
                  
                  float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
                  float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
                  
                  if ( tmp2<tmp1 ) 
                flipped = 1;
                  else 
                flipped = 0;
                  
                  
                  layer  = tTopo->pxbLayer(detId);   // Layer: 1,2,3.
                  ladder = tTopo->pxbLadder(detId);  // Ladder: 1-20, 32, 44. 
                  mod   = tTopo->pxbModule(detId);  // Mod: 1-8.
                }             
              else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
                {
                  
                  side  = tTopo->pxfSide(detId);
                  disk  = tTopo->pxfDisk(detId);
                  blade = tTopo->pxfBlade(detId);
                  panel = tTopo->pxfPanel(detId);
                  plaq  = tTopo->pxfModule(detId); // also known as plaquette
                  
                  float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
                  float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
                  
                  if ( tmp2<tmp1 ) 
                flipped = 1;
                  else 
                flipped = 0;
                  
                } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
              //else std::// cout << "We are not in the pixel detector. detId.subdetId() = " << (int)detId.subdetId() << endl;
              
              ttree_track_hits_->Fill();
              
            } // if ( !matched.empty() )
              else
            cout << "---------------- RecHit with no associated SimHit !!! -------------------------- " << endl;
              
            } // if ( matchedhit )
          
        } // end of loop on hits
          
        } //end of loop on track 
      
    } // tracks > 0.
     
    } // if ( include_trk_hits_ )

 

  // ----------------------------------------------- track hits only -----------------------------------------------------------
  
}

void SiPixelErrorEstimation::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file SiPixelErrorEstimation.cc.

References all_alpha, all_beta, all_bigx, all_bigy, all_blade, all_clust_alpha, all_clust_beta, all_clust_col, all_clust_geoid, all_clust_maxpixcol, all_clust_maxpixrow, all_clust_minpixcol, all_clust_minpixrow, all_clust_q, all_clust_row, all_clust_x, all_clust_y, all_col1, all_col2, all_cols, all_disk, all_edgex, all_edgey, all_eloss, all_flipped, all_gx1, all_gx2, all_half, all_hit_probx, all_hit_proby, all_ladder, all_layer, all_mod, all_npix, all_nsimhit, all_nxpix, all_nypix, all_panel, all_pidhit, all_pixadc, all_pixcol, all_pixel_clst_err_x, all_pixel_clst_err_y, all_pixel_split, all_pixgx, all_pixgy, all_pixgz, all_pixrow, all_pixx, all_pixy, all_plaq, all_rechiterrx, all_rechiterry, all_rechitpullx, all_rechitpully, all_rechitresx, all_rechitresy, all_rechitx, all_rechity, all_rechitz, all_row1, all_row2, all_rows, all_side, all_simhitx, all_simhity, all_simphi, all_simproc, all_simpx, all_simpy, all_simpz, all_simtheta, all_simtrketa, all_simtrkphi, all_subdetid, all_trkid, all_vtxr, all_vtxz, all_x1, all_x2, alpha, beta, bigx, bigy, blade, charge, detector_type, disk, edgex, edgey, eta, evt, flipped, half, include_trk_hits_, ladder, layer, mod, npix, nsimhit, nxpix, nypix, outputFile_, panel, phi, pidhit, pixel_clst_err_x, pixel_clst_err_y, pixel_split, plaq, rechiterrx, rechiterry, rechitpullx, rechitpully, rechitresx, rechitresy, rechitx, rechity, rechitz, run, side, simhitx, simhity, simproc, strip_charge, strip_clst_err_x, strip_clst_err_y, strip_cotalpha, strip_cotbeta, strip_edge, strip_hit_type, strip_is_stereo, strip_locbx, strip_locby, strip_locbz, strip_nprm, strip_nsimhit, strip_pidhit, strip_pidhit1, strip_pidhit2, strip_pidhit3, strip_pidhit4, strip_pidhit5, strip_prob, strip_qbin, strip_rechiterrx, strip_rechiterry, strip_rechitpullx, strip_rechitpully, strip_rechitresx, strip_rechitresx2, strip_rechitresy, strip_rechitx, strip_rechity, strip_rechitz, strip_simproc, strip_simproc1, strip_simproc2, strip_simproc3, strip_simproc4, strip_simproc5, strip_size, strip_split, strip_subdet_id, strip_tib_is_double_side, strip_tib_is_external_string, strip_tib_is_internal_string, strip_tib_is_rphi, strip_tib_is_stereo, strip_tib_is_z_minus_side, strip_tib_is_z_plus_side, strip_tib_layer, strip_tib_layer_number, strip_tib_module, strip_tib_module_number, strip_tib_order, strip_tib_side, strip_tib_string_number, strip_tob_is_double_side, strip_tob_is_z_minus_side, strip_tob_is_z_plus_side, strip_tob_layer, strip_tob_layer_number, strip_tob_module, strip_tob_module_number, strip_tob_rod_number, strip_tob_side, strip_trk_pt, subdetId, tfile_, trk_alpha, trk_beta, ttree_all_hits_, ttree_track_hits_, and ttree_track_hits_strip_.

{
   int bufsize = 64000;


  tfile_ = new TFile ( outputFile_.c_str() , "RECREATE");

  
  ttree_track_hits_strip_ = new TTree("TrackHitNtupleStrip", "TrackHitNtupleStrip");
  
  ttree_track_hits_strip_->Branch("strip_rechitx", &strip_rechitx, "strip_rechitx/F"    , bufsize);
  ttree_track_hits_strip_->Branch("strip_rechity", &strip_rechity, "strip_rechity/F"    , bufsize);
  ttree_track_hits_strip_->Branch("strip_rechitz", &strip_rechitz, "strip_rechitz/F"    , bufsize);
  
  ttree_track_hits_strip_->Branch("strip_rechiterrx", &strip_rechiterrx, "strip_rechiterrx/F" , bufsize);
  ttree_track_hits_strip_->Branch("strip_rechiterry", &strip_rechiterry, "strip_rechiterry/F" , bufsize);

  ttree_track_hits_strip_->Branch("strip_rechitresx", &strip_rechitresx, "strip_rechitresx/F" , bufsize);
  
  ttree_track_hits_strip_->Branch("strip_rechitresx2", &strip_rechitresx2, "strip_rechitresx2/F" , bufsize);


  ttree_track_hits_strip_->Branch("strip_rechitresy", &strip_rechitresy, "strip_rechitresy/F" , bufsize);
  
  ttree_track_hits_strip_->Branch("strip_rechitpullx", &strip_rechitpullx, "strip_rechitpullx/F", bufsize);
  ttree_track_hits_strip_->Branch("strip_rechitpully", &strip_rechitpully, "strip_rechitpully/F", bufsize);

  ttree_track_hits_strip_->Branch("strip_is_stereo", &strip_is_stereo, "strip_is_stereo/I", bufsize);
  ttree_track_hits_strip_->Branch("strip_hit_type" , &strip_hit_type , "strip_hit_type/I" , bufsize);
  ttree_track_hits_strip_->Branch("detector_type"  , &detector_type  , "detector_type/I"  , bufsize);

  ttree_track_hits_strip_->Branch("strip_trk_pt"   , &strip_trk_pt   , "strip_trk_pt/F"   , bufsize);
  ttree_track_hits_strip_->Branch("strip_cotalpha" , &strip_cotalpha , "strip_cotalpha/F" , bufsize);
  ttree_track_hits_strip_->Branch("strip_cotbeta"  , &strip_cotbeta  , "strip_cotbeta/F"  , bufsize);
  ttree_track_hits_strip_->Branch("strip_locbx"    , &strip_locbx    , "strip_locbx/F"    , bufsize);
  ttree_track_hits_strip_->Branch("strip_locby"    , &strip_locby    , "strip_locby/F"    , bufsize);
  ttree_track_hits_strip_->Branch("strip_locbz"    , &strip_locbz    , "strip_locbz/F"    , bufsize);
  ttree_track_hits_strip_->Branch("strip_charge"   , &strip_charge   , "strip_charge/F"   , bufsize);
  ttree_track_hits_strip_->Branch("strip_size"     , &strip_size     , "strip_size/I"     , bufsize);

  
  ttree_track_hits_strip_->Branch("strip_edge"   , &strip_edge   , "strip_edge/I"    , bufsize);
  ttree_track_hits_strip_->Branch("strip_nsimhit", &strip_nsimhit, "strip_nsimhit/I" , bufsize);
  ttree_track_hits_strip_->Branch("strip_pidhit" , &strip_pidhit , "strip_pidhit/I"  , bufsize);
  ttree_track_hits_strip_->Branch("strip_simproc", &strip_simproc, "strip_simproc/I" , bufsize);


ttree_track_hits_strip_->Branch("strip_subdet_id"             , &strip_subdet_id             , "strip_subdet_id/I"            , bufsize);
 
ttree_track_hits_strip_->Branch("strip_tib_layer"             , &strip_tib_layer             , "strip_tib_layer/I"            , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_module"            , &strip_tib_module            , "strip_tib_module/I"           , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_order"             , &strip_tib_order             , "strip_tib_order/I"            , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_side"              , &strip_tib_side              , "strip_tib_side/I"             , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_is_double_side"    , &strip_tib_is_double_side    , "strip_tib_is_double_side/I"   , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_is_z_plus_side"    , &strip_tib_is_z_plus_side    , "strip_tib_is_z_plus_side/I"   , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_is_z_minus_side"   , &strip_tib_is_z_minus_side   , "strip_tib_is_z_minus_side/I"  , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_layer_number"      , &strip_tib_layer_number      , "strip_tib_layer_number/I"     , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_string_number"     , &strip_tib_string_number     , "strip_tib_string_number/I"    , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_module_number"     , &strip_tib_module_number     ,"strip_tib_module_number/I"     , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_is_internal_string", &strip_tib_is_internal_string,"strip_tib_is_internal_string/I", bufsize);
ttree_track_hits_strip_->Branch("strip_tib_is_external_string", &strip_tib_is_external_string,"strip_tib_is_external_string/I", bufsize);
ttree_track_hits_strip_->Branch("strip_tib_is_rphi"           , &strip_tib_is_rphi           , "strip_tib_is_rphi/I"          , bufsize);
ttree_track_hits_strip_->Branch("strip_tib_is_stereo"         , &strip_tib_is_stereo         , "strip_tib_is_stereo/I"        , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_layer"             , &strip_tob_layer             , "strip_tob_layer/I"            , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_module"            , &strip_tob_module            , "strip_tob_module/I"           , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_side"              , &strip_tob_side              , "strip_tob_side/I"             , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_is_double_side"    , &strip_tob_is_double_side    , "strip_tob_is_double_side/I"   , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_is_z_plus_side"    , &strip_tob_is_z_plus_side    , "strip_tob_is_z_plus_side/I"   , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_is_z_minus_side"   , &strip_tob_is_z_minus_side   , "strip_tob_is_z_minus_side/I"  , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_layer_number"      , &strip_tob_layer_number      , "strip_tob_layer_number/I"     , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_rod_number"        , &strip_tob_rod_number        , "strip_tob_rod_number/I"       , bufsize);
ttree_track_hits_strip_->Branch("strip_tob_module_number"     , &strip_tob_module_number     , "strip_tob_module_number/I"    , bufsize);


ttree_track_hits_strip_->Branch("strip_prob", &strip_prob, "strip_prob/F"   , bufsize);
ttree_track_hits_strip_->Branch("strip_qbin", &strip_qbin, "strip_qbin/I", bufsize);

ttree_track_hits_strip_->Branch("strip_nprm", &strip_nprm, "strip_nprm/I", bufsize);

ttree_track_hits_strip_->Branch("strip_pidhit1" , &strip_pidhit1 , "strip_pidhit1/I"  , bufsize);
ttree_track_hits_strip_->Branch("strip_simproc1", &strip_simproc1, "strip_simproc1/I" , bufsize);

ttree_track_hits_strip_->Branch("strip_pidhit2" , &strip_pidhit2 , "strip_pidhit2/I"  , bufsize);
ttree_track_hits_strip_->Branch("strip_simproc2", &strip_simproc2, "strip_simproc2/I" , bufsize);

ttree_track_hits_strip_->Branch("strip_pidhit3" , &strip_pidhit3 , "strip_pidhit3/I"  , bufsize);
ttree_track_hits_strip_->Branch("strip_simproc3", &strip_simproc3, "strip_simproc3/I" , bufsize);

ttree_track_hits_strip_->Branch("strip_pidhit4" , &strip_pidhit4 , "strip_pidhit4/I"  , bufsize);
ttree_track_hits_strip_->Branch("strip_simproc4", &strip_simproc4, "strip_simproc4/I" , bufsize);

ttree_track_hits_strip_->Branch("strip_pidhit5" , &strip_pidhit5 , "strip_pidhit5/I"  , bufsize);
ttree_track_hits_strip_->Branch("strip_simproc5", &strip_simproc5, "strip_simproc5/I" , bufsize);

ttree_track_hits_strip_->Branch("strip_split", &strip_split, "strip_split/I" , bufsize);

ttree_track_hits_strip_->Branch("strip_clst_err_x", &strip_clst_err_x, "strip_clst_err_x/F"   , bufsize);
ttree_track_hits_strip_->Branch("strip_clst_err_y", &strip_clst_err_y, "strip_clst_err_y/F"   , bufsize);

 if ( include_trk_hits_ )
   {
     //tfile_ = new TFile ("SiPixelErrorEstimation_Ntuple.root" , "RECREATE");
     //const char* tmp_name = outputFile_.c_str();
     
     
     ttree_track_hits_ = new TTree("TrackHitNtuple", "TrackHitNtuple");
     
      ttree_track_hits_->Branch("evt", &evt, "evt/I", bufsize);
      ttree_track_hits_->Branch("run", &run, "run/I", bufsize);
      
      ttree_track_hits_->Branch("subdetId", &subdetId, "subdetId/I", bufsize);
      
      ttree_track_hits_->Branch("layer" , &layer , "layer/I" , bufsize);
      ttree_track_hits_->Branch("ladder", &ladder, "ladder/I", bufsize);
      ttree_track_hits_->Branch("mod"   , &mod   , "mod/I"   , bufsize);
      
      ttree_track_hits_->Branch("side"  , &side  , "side/I"  , bufsize);
      ttree_track_hits_->Branch("disk"  , &disk  , "disk/I"  , bufsize);
      ttree_track_hits_->Branch("blade" , &blade , "blade/I" , bufsize);
      ttree_track_hits_->Branch("panel" , &panel , "panel/I" , bufsize);
      ttree_track_hits_->Branch("plaq"  , &plaq  , "plaq/I"  , bufsize);
      
      ttree_track_hits_->Branch("half"   , &half   , "half/I"   , bufsize);
      ttree_track_hits_->Branch("flipped", &flipped, "flipped/I", bufsize);
      
      ttree_track_hits_->Branch("rechitx", &rechitx, "rechitx/F"    , bufsize);
      ttree_track_hits_->Branch("rechity", &rechity, "rechity/F"    , bufsize);
      ttree_track_hits_->Branch("rechitz", &rechitz, "rechitz/F"    , bufsize);
      
      ttree_track_hits_->Branch("rechiterrx", &rechiterrx, "rechiterrx/F" , bufsize);
      ttree_track_hits_->Branch("rechiterry", &rechiterry, "rechiterry/F" , bufsize);
      
      ttree_track_hits_->Branch("rechitresx", &rechitresx, "rechitresx/F" , bufsize);
      ttree_track_hits_->Branch("rechitresy", &rechitresy, "rechitresy/F" , bufsize);
      
      ttree_track_hits_->Branch("rechitpullx", &rechitpullx, "rechitpullx/F", bufsize);
      ttree_track_hits_->Branch("rechitpully", &rechitpully, "rechitpully/F", bufsize);

      
      ttree_track_hits_->Branch("npix"  , &npix  , "npix/I"  , bufsize);
      ttree_track_hits_->Branch("nxpix" , &nxpix , "nxpix/I" , bufsize);
      ttree_track_hits_->Branch("nypix" , &nypix , "nypix/I" , bufsize);
      ttree_track_hits_->Branch("charge", &charge, "charge/F", bufsize);
      
      ttree_track_hits_->Branch("edgex", &edgex, "edgex/I", bufsize);
      ttree_track_hits_->Branch("edgey", &edgey, "edgey/I", bufsize);
      
      ttree_track_hits_->Branch("bigx", &bigx, "bigx/I", bufsize);
      ttree_track_hits_->Branch("bigy", &bigy, "bigy/I", bufsize);
      
      ttree_track_hits_->Branch("alpha", &alpha, "alpha/F", bufsize);
      ttree_track_hits_->Branch("beta" , &beta , "beta/F" , bufsize);
      
      ttree_track_hits_->Branch("trk_alpha", &trk_alpha, "trk_alpha/F", bufsize);
      ttree_track_hits_->Branch("trk_beta" , &trk_beta , "trk_beta/F" , bufsize);

      ttree_track_hits_->Branch("phi", &phi, "phi/F", bufsize);
      ttree_track_hits_->Branch("eta", &eta, "eta/F", bufsize);
      
      ttree_track_hits_->Branch("simhitx", &simhitx, "simhitx/F", bufsize);
      ttree_track_hits_->Branch("simhity", &simhity, "simhity/F", bufsize);
      
      ttree_track_hits_->Branch("nsimhit", &nsimhit, "nsimhit/I", bufsize);
      ttree_track_hits_->Branch("pidhit" , &pidhit , "pidhit/I" , bufsize);
      ttree_track_hits_->Branch("simproc", &simproc, "simproc/I", bufsize);
      
      ttree_track_hits_->Branch("pixel_split", &pixel_split, "pixel_split/I", bufsize);

      ttree_track_hits_->Branch("pixel_clst_err_x", &pixel_clst_err_x, "pixel_clst_err_x/F"   , bufsize);
      ttree_track_hits_->Branch("pixel_clst_err_y", &pixel_clst_err_y, "pixel_clst_err_y/F"   , bufsize);

    } // if ( include_trk_hits_ )

  // ----------------------------------------------------------------------
  
  ttree_all_hits_ = new TTree("AllHitNtuple", "AllHitNtuple");

  ttree_all_hits_->Branch("evt", &evt, "evt/I", bufsize);
  ttree_all_hits_->Branch("run", &run, "run/I", bufsize);

  ttree_all_hits_->Branch("subdetid", &all_subdetid, "subdetid/I", bufsize);
  
  ttree_all_hits_->Branch("layer" , &all_layer , "layer/I" , bufsize);
  ttree_all_hits_->Branch("ladder", &all_ladder, "ladder/I", bufsize);
  ttree_all_hits_->Branch("mod"   , &all_mod   , "mod/I"   , bufsize);
  
  ttree_all_hits_->Branch("side"  , &all_side  , "side/I"  , bufsize);
  ttree_all_hits_->Branch("disk"  , &all_disk  , "disk/I"  , bufsize);
  ttree_all_hits_->Branch("blade" , &all_blade , "blade/I" , bufsize);
  ttree_all_hits_->Branch("panel" , &all_panel , "panel/I" , bufsize);
  ttree_all_hits_->Branch("plaq"  , &all_plaq  , "plaq/I"  , bufsize);

  ttree_all_hits_->Branch("half"   , &all_half   , "half/I"   , bufsize);
  ttree_all_hits_->Branch("flipped", &all_flipped, "flipped/I", bufsize);

  ttree_all_hits_->Branch("cols", &all_cols, "cols/I", bufsize);
  ttree_all_hits_->Branch("rows", &all_rows, "rows/I", bufsize);

  ttree_all_hits_->Branch("rechitx"    , &all_rechitx    , "rechitx/F"    , bufsize);
  ttree_all_hits_->Branch("rechity"    , &all_rechity    , "rechity/F"    , bufsize);
  ttree_all_hits_->Branch("rechitz"    , &all_rechitz    , "rechitz/F"    , bufsize);
 
  ttree_all_hits_->Branch("rechiterrx" , &all_rechiterrx , "rechiterrx/F" , bufsize);
  ttree_all_hits_->Branch("rechiterry" , &all_rechiterry , "rechiterry/F" , bufsize);
  
  ttree_all_hits_->Branch("rechitresx" , &all_rechitresx , "rechitresx/F" , bufsize);
  ttree_all_hits_->Branch("rechitresy" , &all_rechitresy , "rechitresy/F" , bufsize);
  
  ttree_all_hits_->Branch("rechitpullx", &all_rechitpullx, "rechitpullx/F", bufsize);
  ttree_all_hits_->Branch("rechitpully", &all_rechitpully, "rechitpully/F", bufsize);

  ttree_all_hits_->Branch("npix"  , &all_npix  , "npix/I"  , bufsize);
  ttree_all_hits_->Branch("nxpix" , &all_nxpix , "nxpix/I" , bufsize);
  ttree_all_hits_->Branch("nypix" , &all_nypix , "nypix/I" , bufsize);

  ttree_all_hits_->Branch("edgex", &all_edgex, "edgex/I", bufsize);
  ttree_all_hits_->Branch("edgey", &all_edgey, "edgey/I", bufsize);
 
  ttree_all_hits_->Branch("bigx", &all_bigx, "bigx/I", bufsize);
  ttree_all_hits_->Branch("bigy", &all_bigy, "bigy/I", bufsize);
  
  ttree_all_hits_->Branch("alpha", &all_alpha, "alpha/F", bufsize);
  ttree_all_hits_->Branch("beta" , &all_beta , "beta/F" , bufsize);

  ttree_all_hits_->Branch("simhitx", &all_simhitx, "simhitx/F", bufsize);
  ttree_all_hits_->Branch("simhity", &all_simhity, "simhity/F", bufsize);

  ttree_all_hits_->Branch("nsimhit", &all_nsimhit, "nsimhit/I", bufsize);
  ttree_all_hits_->Branch("pidhit" , &all_pidhit , "pidhit/I" , bufsize);
  ttree_all_hits_->Branch("simproc", &all_simproc, "simproc/I", bufsize);

  ttree_all_hits_->Branch("vtxr", &all_vtxr, "vtxr/F", bufsize);
  ttree_all_hits_->Branch("vtxz", &all_vtxz, "vtxz/F", bufsize);

  ttree_all_hits_->Branch("simpx", &all_simpx, "simpx/F", bufsize);
  ttree_all_hits_->Branch("simpy", &all_simpy, "simpy/F", bufsize);
  ttree_all_hits_->Branch("simpz", &all_simpz, "simpz/F", bufsize);

  ttree_all_hits_->Branch("eloss", &all_eloss, "eloss/F", bufsize);
  
  ttree_all_hits_->Branch("simphi", &all_simphi, "simphi/F", bufsize);
  ttree_all_hits_->Branch("simtheta", &all_simtheta, "simtheta/F", bufsize);
  
  ttree_all_hits_->Branch("trkid", &all_trkid, "trkid/I", bufsize);
  
  ttree_all_hits_->Branch("x1", &all_x1, "x1/F", bufsize);
  ttree_all_hits_->Branch("x2", &all_x2, "x2/F", bufsize);
  ttree_all_hits_->Branch("y1", &all_x1, "y1/F", bufsize);
  ttree_all_hits_->Branch("y2", &all_x2, "y2/F", bufsize);
  ttree_all_hits_->Branch("z1", &all_x1, "z1/F", bufsize);
  ttree_all_hits_->Branch("z2", &all_x2, "z2/F", bufsize);

  ttree_all_hits_->Branch("row1", &all_row1, "row1/F", bufsize);
  ttree_all_hits_->Branch("row2", &all_row2, "row2/F", bufsize);
  ttree_all_hits_->Branch("col1", &all_col1, "col1/F", bufsize);
  ttree_all_hits_->Branch("col2", &all_col2, "col2/F", bufsize);

  ttree_all_hits_->Branch("gx1", &all_gx1, "gx1/F", bufsize);
  ttree_all_hits_->Branch("gx2", &all_gx2, "gx2/F", bufsize);
  ttree_all_hits_->Branch("gy1", &all_gx1, "gy1/F", bufsize);
  ttree_all_hits_->Branch("gy2", &all_gx2, "gy2/F", bufsize);
  ttree_all_hits_->Branch("gz1", &all_gx1, "gz1/F", bufsize);
  ttree_all_hits_->Branch("gz2", &all_gx2, "gz2/F", bufsize);
  
  ttree_all_hits_->Branch("simtrketa", &all_simtrketa, "simtrketa/F", bufsize);
  ttree_all_hits_->Branch("simtrkphi", &all_simtrkphi, "simtrkphi/F", bufsize);

  ttree_all_hits_->Branch("clust_row", &all_clust_row, "clust_row/F", bufsize);
  ttree_all_hits_->Branch("clust_col", &all_clust_col, "clust_col/F", bufsize);
  
  ttree_all_hits_->Branch("clust_x", &all_clust_x, "clust_x/F", bufsize);
  ttree_all_hits_->Branch("clust_y", &all_clust_y, "clust_y/F", bufsize);

  ttree_all_hits_->Branch("clust_q", &all_clust_q, "clust_q/F", bufsize);

  ttree_all_hits_->Branch("clust_maxpixcol", &all_clust_maxpixcol, "clust_maxpixcol/I", bufsize);
  ttree_all_hits_->Branch("clust_maxpixrow", &all_clust_maxpixrow, "clust_maxpixrow/I", bufsize);
  ttree_all_hits_->Branch("clust_minpixcol", &all_clust_minpixcol, "clust_minpixcol/I", bufsize);
  ttree_all_hits_->Branch("clust_minpixrow", &all_clust_minpixrow, "clust_minpixrow/I", bufsize);
  
  ttree_all_hits_->Branch("clust_geoid", &all_clust_geoid, "clust_geoid/I", bufsize);
  
  ttree_all_hits_->Branch("clust_alpha", &all_clust_alpha, "clust_alpha/F", bufsize);
  ttree_all_hits_->Branch("clust_beta" , &all_clust_beta , "clust_beta/F" , bufsize);

  ttree_all_hits_->Branch("rowpix", all_pixrow, "row[npix]/F", bufsize);
  ttree_all_hits_->Branch("colpix", all_pixcol, "col[npix]/F", bufsize);
  ttree_all_hits_->Branch("adc", all_pixadc, "adc[npix]/F", bufsize);
  ttree_all_hits_->Branch("xpix", all_pixx, "x[npix]/F", bufsize);
  ttree_all_hits_->Branch("ypix", all_pixy, "y[npix]/F", bufsize);
  ttree_all_hits_->Branch("gxpix", all_pixgx, "gx[npix]/F", bufsize);
  ttree_all_hits_->Branch("gypix", all_pixgy, "gy[npix]/F", bufsize);
  ttree_all_hits_->Branch("gzpix", all_pixgz, "gz[npix]/F", bufsize);
  
  ttree_all_hits_->Branch("hit_probx", &all_hit_probx, "hit_probx/F" , bufsize);
  ttree_all_hits_->Branch("hit_proby", &all_hit_proby, "hit_proby/F" , bufsize);

  ttree_all_hits_->Branch("all_pixel_split", &all_pixel_split, "all_pixel_split/I" , bufsize);

  ttree_all_hits_->Branch("all_pixel_clst_err_x", &all_pixel_clst_err_x, "all_pixel_clst_err_x/F"   , bufsize);
  ttree_all_hits_->Branch("all_pixel_clst_err_y", &all_pixel_clst_err_y, "all_pixel_clst_err_y/F"   , bufsize);

}

void SiPixelErrorEstimation::computeAnglesFromDetPosition	(	const SiPixelCluster &	cl,
		const GeomDetUnit &	det,
		float &	alpha,
		float &	beta
	)

Definition at line 1673 of file SiPixelErrorEstimation.cc.

References assert(), gather_cfg::cout, Topology::localPosition(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), SiPixelCluster::x(), PV3DBase< T, PVType, FrameType >::y(), SiPixelCluster::y(), and PV3DBase< T, PVType, FrameType >::z().

{
  //--- This is a new det unit, so cache it
  const PixelGeomDetUnit* theDet = dynamic_cast<const PixelGeomDetUnit*>( &det );
  if (! theDet) 
    {
      cout << "---------------------------------------------- Not a pixel detector !!!!!!!!!!!!!!" << endl;
      assert(0);
    }

  const PixelTopology* theTopol = &(theDet->specificTopology());

  // get cluster center of gravity (of charge)
  float xcenter = cl.x();
  float ycenter = cl.y();
  
  // get the cluster position in local coordinates (cm) 
  LocalPoint lp = theTopol->localPosition( MeasurementPoint(xcenter, ycenter) );
  //float lp_mod = sqrt( lp.x()*lp.x() + lp.y()*lp.y() + lp.z()*lp.z() );

  // get the cluster position in global coordinates (cm)
  GlobalPoint gp = theDet->surface().toGlobal( lp );
  float gp_mod = sqrt( gp.x()*gp.x() + gp.y()*gp.y() + gp.z()*gp.z() );

  // normalize
  float gpx = gp.x()/gp_mod;
  float gpy = gp.y()/gp_mod;
  float gpz = gp.z()/gp_mod;

  // make a global vector out of the global point; this vector will point from the 
  // origin of the detector to the cluster
  GlobalVector gv(gpx, gpy, gpz);

  // make local unit vector along local X axis
  const Local3DVector lvx(1.0, 0.0, 0.0);

  // get the unit X vector in global coordinates/
  GlobalVector gvx = theDet->surface().toGlobal( lvx );

  // make local unit vector along local Y axis
  const Local3DVector lvy(0.0, 1.0, 0.0);

  // get the unit Y vector in global coordinates
  GlobalVector gvy = theDet->surface().toGlobal( lvy );
   
  // make local unit vector along local Z axis
  const Local3DVector lvz(0.0, 0.0, 1.0);

  // get the unit Z vector in global coordinates
  GlobalVector gvz = theDet->surface().toGlobal( lvz );
    
  // calculate the components of gv (the unit vector pointing to the cluster) 
  // in the local coordinate system given by the basis {gvx, gvy, gvz}
  // note that both gv and the basis {gvx, gvy, gvz} are given in global coordinates
  float gv_dot_gvx = gv.x()*gvx.x() + gv.y()*gvx.y() + gv.z()*gvx.z();
  float gv_dot_gvy = gv.x()*gvy.x() + gv.y()*gvy.y() + gv.z()*gvy.z();
  float gv_dot_gvz = gv.x()*gvz.x() + gv.y()*gvz.y() + gv.z()*gvz.z();

  // calculate angles
  alpha = atan2( gv_dot_gvz, gv_dot_gvx );
  beta  = atan2( gv_dot_gvz, gv_dot_gvy );

  // calculate cotalpha and cotbeta
  //   cotalpha_ = 1.0/tan(alpha_);
  //   cotbeta_  = 1.0/tan(beta_ );
  // or like this
  //cotalpha_ = gv_dot_gvx / gv_dot_gvz;
  //cotbeta_  = gv_dot_gvy / gv_dot_gvz;
}

void SiPixelErrorEstimation::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 365 of file SiPixelErrorEstimation.cc.

References tfile_.

{
  tfile_->Write();
  tfile_->Close();
}

Member Data Documentation

float SiPixelErrorEstimation::all_alpha

private

Definition at line 292 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_beta

private

Definition at line 293 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_bigx

private

Definition at line 289 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_bigy

private

Definition at line 290 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_blade

private

Definition at line 256 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_alpha

private

Definition at line 350 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_beta

private

Definition at line 351 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_col

private

Definition at line 336 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_geoid

private

Definition at line 348 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_maxpixcol

private

Definition at line 343 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_maxpixrow

private

Definition at line 344 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_minpixcol

private

Definition at line 345 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_minpixrow

private

Definition at line 346 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_q

private

Definition at line 341 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_row

private

Definition at line 335 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_x

private

Definition at line 338 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_y

private

Definition at line 339 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_col1

private

Definition at line 322 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_col2

private

Definition at line 323 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_cols

private

Definition at line 263 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_disk

private

Definition at line 255 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_edgex

private

Definition at line 286 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_edgey

private

Definition at line 287 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_eloss

private

Definition at line 309 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_flipped

private

Definition at line 261 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_gx1

private

Definition at line 325 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_gx2

private

Definition at line 326 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_gy1

private

Definition at line 327 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_gy2

private

Definition at line 328 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_gz1

private

Definition at line 329 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_gz2

private

Definition at line 330 of file SiPixelErrorEstimation.h.

Referenced by analyze().

int SiPixelErrorEstimation::all_half

private

Definition at line 260 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_hit_cprob0

private

Definition at line 366 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_hit_cprob1

private

Definition at line 367 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_hit_cprob2

private

Definition at line 368 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_hit_probx

private

Definition at line 364 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_hit_proby

private

Definition at line 365 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_ladder

private

Definition at line 251 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_layer

private

Definition at line 250 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_mod

private

Definition at line 252 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_npix

private

Definition at line 282 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_nsimhit

private

Definition at line 298 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_nxpix

private

Definition at line 283 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_nypix

private

Definition at line 284 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_panel

private

Definition at line 257 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_pidhit

private

Definition at line 299 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixadc[maxpix]

private

Definition at line 356 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixcol[maxpix]

private

Definition at line 355 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixel_clst_err_x

private

Definition at line 371 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixel_clst_err_y

private

Definition at line 372 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_pixel_split

private

Definition at line 370 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixgx[maxpix]

private

Definition at line 360 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixgy[maxpix]

private

Definition at line 361 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixgz[maxpix]

private

Definition at line 362 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixrow[maxpix]

private

Definition at line 354 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixx[maxpix]

private

Definition at line 358 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixy[maxpix]

private

Definition at line 359 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_plaq

private

Definition at line 258 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechiterrx

private

Definition at line 273 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechiterry

private

Definition at line 274 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitpullx

private

Definition at line 279 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitpully

private

Definition at line 280 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitresx

private

Definition at line 276 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitresy

private

Definition at line 277 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitx

private

Definition at line 266 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechity

private

Definition at line 267 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitz

private

Definition at line 268 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_row1

private

Definition at line 320 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_row2

private

Definition at line 321 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_rows

private

Definition at line 264 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_side

private

Definition at line 254 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simhitx

private

Definition at line 270 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simhity

private

Definition at line 271 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simphi

private

Definition at line 295 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_simproc

private

Definition at line 300 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simpx

private

Definition at line 305 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simpy

private

Definition at line 306 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simpz

private

Definition at line 307 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simtheta

private

Definition at line 296 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simtrketa

private

Definition at line 332 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simtrkphi

private

Definition at line 333 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_subdetid

private

Definition at line 248 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_trkid

private

Definition at line 311 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_vtxr

private

Definition at line 302 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_vtxz

private

Definition at line 303 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_x1

private

Definition at line 313 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_x2

private

Definition at line 314 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_y1

private

Definition at line 315 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_y2

private

Definition at line 316 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_z1

private

Definition at line 317 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_z2

private

Definition at line 318 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::alpha

private

Definition at line 202 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::beta

private

Definition at line 203 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::bigx

private

Definition at line 199 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::bigy

private

Definition at line 200 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::blade

private

Definition at line 217 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::charge

private

Definition at line 194 of file SiPixelErrorEstimation.h.

Referenced by Vispa.Plugins.EdmBrowser.ParticleDataList.ParticleData::__repr__(), analyze(), and beginJob().

bool SiPixelErrorEstimation::checkType_

private

Definition at line 87 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

edm::ParameterSet SiPixelErrorEstimation::conf_

private

Definition at line 80 of file SiPixelErrorEstimation.h.

int SiPixelErrorEstimation::detector_type

private

Definition at line 121 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::disk

private

Definition at line 216 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::edgex

private

Definition at line 196 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::edgey

private

Definition at line 197 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::eta

private

Definition at line 209 of file SiPixelErrorEstimation.h.

Referenced by Particle.Particle::__str__(), analyze(), beginJob(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

int SiPixelErrorEstimation::evt

private

Definition at line 231 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::flipped

private

Definition at line 222 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::genType_

private

Definition at line 88 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

int SiPixelErrorEstimation::half

private

Definition at line 221 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::hit_cprob0

private

Definition at line 236 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_cprob1

private

Definition at line 237 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_cprob2

private

Definition at line 238 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_probx

private

Definition at line 234 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_proby

private

Definition at line 235 of file SiPixelErrorEstimation.h.

Referenced by analyze().

bool SiPixelErrorEstimation::include_trk_hits_

private

Definition at line 89 of file SiPixelErrorEstimation.h.

Referenced by analyze(), beginJob(), and SiPixelErrorEstimation().

int SiPixelErrorEstimation::ladder

private

Definition at line 213 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::layer

private

Definition at line 212 of file SiPixelErrorEstimation.h.

Referenced by analyze(), beginJob(), geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

const int SiPixelErrorEstimation::maxpix = 10000

staticprivate

Definition at line 353 of file SiPixelErrorEstimation.h.

int SiPixelErrorEstimation::mod

private

Definition at line 214 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::npix

private

Definition at line 191 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::nsimhit

private

Definition at line 224 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::nxpix

private

Definition at line 192 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::nypix

private

Definition at line 193 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

std::string SiPixelErrorEstimation::outputFile_

private

Definition at line 85 of file SiPixelErrorEstimation.h.

Referenced by beginJob(), and SiPixelErrorEstimation().

int SiPixelErrorEstimation::panel

private

Definition at line 218 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::phi

private

Definition at line 208 of file SiPixelErrorEstimation.h.

Referenced by Particle.Particle::__str__(), analyze(), and beginJob().

int SiPixelErrorEstimation::pidhit

private

Definition at line 225 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::pixel_clst_err_x

private

Definition at line 242 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::pixel_clst_err_y

private

Definition at line 243 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::pixel_split

private

Definition at line 240 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::plaq

private

Definition at line 219 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechiterrx

private

Definition at line 95 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechiterry

private

Definition at line 96 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitpullx

private

Definition at line 99 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitpully

private

Definition at line 100 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitresx

private

Definition at line 97 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitresy

private

Definition at line 98 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitx

private

Definition at line 92 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechity

private

Definition at line 93 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitz

private

Definition at line 94 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::run

private

Definition at line 232 of file SiPixelErrorEstimation.h.

Referenced by generateEDF.LumiInfo::__str__(), MatrixUtil.InputInfo::__str__(), analyze(), beginJob(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), MatrixUtil.InputInfo::das(), personalPlayback.Playback::do_create_lumi(), personalPlayback.Playback::do_exec(), generateEDF.LumiInfo::fixXingInfo(), MatrixUtil.InputInfo::lumiRanges(), MatrixUtil.InputInfo::queries(), and dqm_interfaces.DirWalkerFile::walk().

int SiPixelErrorEstimation::side

private

Definition at line 215 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::simhitx

private

Definition at line 228 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::simhity

private

Definition at line 229 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::simproc

private

Definition at line 226 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

std::string SiPixelErrorEstimation::src_

private

Definition at line 86 of file SiPixelErrorEstimation.h.

Referenced by SiPixelErrorEstimation().

float SiPixelErrorEstimation::strip_charge

private

Definition at line 129 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_clst_err_x

private

Definition at line 188 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_clst_err_y

private

Definition at line 189 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_cotalpha

private

Definition at line 124 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_cotbeta

private

Definition at line 125 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_edge

private

Definition at line 131 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_hit_type

private

Definition at line 120 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_is_stereo

private

Definition at line 119 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_locbx

private

Definition at line 126 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_locby

private

Definition at line 127 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_locbz

private

Definition at line 128 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_nprm

private

Definition at line 170 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_nsimhit

private

Definition at line 132 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit

private

Definition at line 133 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit1

private

Definition at line 172 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit2

private

Definition at line 175 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit3

private

Definition at line 178 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit4

private

Definition at line 181 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit5

private

Definition at line 184 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_prob

private

Definition at line 167 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_qbin

private

Definition at line 168 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechiterrx

private

Definition at line 109 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechiterry

private

Definition at line 110 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitpullx

private

Definition at line 117 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitpully

private

Definition at line 118 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitresx

private

Definition at line 111 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitresx2

private

Definition at line 114 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitresy

private

Definition at line 116 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitx

private

Definition at line 106 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechity

private

Definition at line 107 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitz

private

Definition at line 108 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc

private

Definition at line 134 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc1

private

Definition at line 173 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc2

private

Definition at line 176 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc3

private

Definition at line 179 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc4

private

Definition at line 182 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc5

private

Definition at line 185 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_size

private

Definition at line 130 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_split

private

Definition at line 187 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_subdet_id

private

Definition at line 136 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_double_side

private

Definition at line 142 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_external_string

private

Definition at line 149 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_internal_string

private

Definition at line 148 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_rphi

private

Definition at line 150 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_stereo

private

Definition at line 151 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_z_minus_side

private

Definition at line 144 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_z_plus_side

private

Definition at line 143 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_layer

private

Definition at line 138 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_layer_number

private

Definition at line 145 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_module

private

Definition at line 139 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_module_number

private

Definition at line 147 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_order

private

Definition at line 140 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_side

private

Definition at line 141 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_string_number

private

Definition at line 146 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_is_double_side

private

Definition at line 157 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_is_rphi

private

Definition at line 164 of file SiPixelErrorEstimation.h.

Referenced by analyze().

int SiPixelErrorEstimation::strip_tob_is_stereo

private

Definition at line 165 of file SiPixelErrorEstimation.h.

Referenced by analyze().

int SiPixelErrorEstimation::strip_tob_is_z_minus_side

private

Definition at line 159 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_is_z_plus_side

private

Definition at line 158 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_layer

private

Definition at line 153 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_layer_number

private

Definition at line 160 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_module

private

Definition at line 154 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_module_number

private

Definition at line 162 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_rod_number

private

Definition at line 161 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_side

private

Definition at line 156 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_trk_pt

private

Definition at line 123 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::subdetId

private

Definition at line 211 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

TFile* SiPixelErrorEstimation::tfile_

private

Definition at line 377 of file SiPixelErrorEstimation.h.

Referenced by beginJob(), and endJob().

edm::EDGetTokenT<SiPixelRecHitCollection> SiPixelErrorEstimation::tPixRecHitCollection

private

Definition at line 82 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

TrackerHitAssociator::Config SiPixelErrorEstimation::trackerHitAssociatorConfig_

private

Definition at line 383 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::trk_alpha

private

Definition at line 205 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::trk_beta

private

Definition at line 206 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

edm::EDGetTokenT<edm::SimTrackContainer> SiPixelErrorEstimation::tSimTrackContainer

private

Definition at line 83 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

edm::EDGetTokenT<reco::TrackCollection> SiPixelErrorEstimation::tTrackCollection

private

Definition at line 84 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

edm::EDGetTokenT<std::vector<Trajectory> > SiPixelErrorEstimation::tTrajectory

private

Definition at line 81 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

TTree* SiPixelErrorEstimation::ttree_all_hits_

private

Definition at line 378 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

TTree* SiPixelErrorEstimation::ttree_track_hits_

private

Definition at line 379 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

TTree* SiPixelErrorEstimation::ttree_track_hits_strip_

private

Definition at line 381 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().