#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
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) 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_

float	trk_alpha

float	trk_beta

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

typedef WorkerT< EDAnalyzer >	WorkerType

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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

CurrentProcessingContext const *	currentContext () const

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 60 of file SiPixelErrorEstimation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 35 of file SiPixelErrorEstimation.cc.

References checkType_, genType_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), include_trk_hits_, outputFile_, src_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                        :tfile_(0), ttree_all_hits_(0), 
                                         ttree_track_hits_(0), ttree_track_hits_strip_(0) 
 {
   //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" );
 }

SiPixelErrorEstimation::~SiPixelErrorEstimation ( )

virtual

Definition at line 50 of file SiPixelErrorEstimation.cc.

51 {}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 364 of file SiPixelErrorEstimation.cc.

 {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<IdealGeometryRecord>().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);
 
   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.getByLabel( src_, 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 std::vector<uint8_t>& stripCharges = cluster->amplitudes();
           uint16_t charge = 0;
           for (unsigned int i = 0; i < stripCharges.size(); ++i) 
         {
           charge += stripCharges.at(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
           const std::vector<uint8_t>& stripCharges = cluster->amplitudes();
           uint16_t charge = 0;
           for (unsigned int i = 0; i < stripCharges.size(); ++i) 
         {
           charge += stripCharges.at(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.getByLabel( "siPixelRecHits", recHitColl);
 
   Handle<edm::SimTrackContainer> simtracks;
   e.getByLabel("g4SimHits", 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.getByLabel(src_, 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 53 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 1665 of file SiPixelErrorEstimation.cc.

References 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 357 of file SiPixelErrorEstimation.cc.

References tfile_.

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

Member Data Documentation

float SiPixelErrorEstimation::all_alpha

private

Definition at line 285 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_beta

private

Definition at line 286 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_bigx

private

Definition at line 282 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_bigy

private

Definition at line 283 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_blade

private

Definition at line 249 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_alpha

private

Definition at line 343 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_beta

private

Definition at line 344 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_col

private

Definition at line 329 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_geoid

private

Definition at line 341 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_maxpixcol

private

Definition at line 336 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_maxpixrow

private

Definition at line 337 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_minpixcol

private

Definition at line 338 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_clust_minpixrow

private

Definition at line 339 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_q

private

Definition at line 334 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_row

private

Definition at line 328 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_x

private

Definition at line 331 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_clust_y

private

Definition at line 332 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_col1

private

Definition at line 315 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_col2

private

Definition at line 316 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_cols

private

Definition at line 256 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_disk

private

Definition at line 248 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_edgex

private

Definition at line 279 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_edgey

private

Definition at line 280 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_eloss

private

Definition at line 302 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_flipped

private

Definition at line 254 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_gx1

private

Definition at line 318 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_gx2

private

Definition at line 319 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_gy1

private

Definition at line 320 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_gy2

private

Definition at line 321 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_gz1

private

Definition at line 322 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_gz2

private

Definition at line 323 of file SiPixelErrorEstimation.h.

Referenced by analyze().

int SiPixelErrorEstimation::all_half

private

Definition at line 253 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_hit_cprob0

private

Definition at line 359 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_hit_cprob1

private

Definition at line 360 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_hit_cprob2

private

Definition at line 361 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_hit_probx

private

Definition at line 357 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_hit_proby

private

Definition at line 358 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_ladder

private

Definition at line 244 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_layer

private

Definition at line 243 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_mod

private

Definition at line 245 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_npix

private

Definition at line 275 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_nsimhit

private

Definition at line 291 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_nxpix

private

Definition at line 276 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_nypix

private

Definition at line 277 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_panel

private

Definition at line 250 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_pidhit

private

Definition at line 292 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixadc[maxpix]

private

Definition at line 349 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixcol[maxpix]

private

Definition at line 348 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixel_clst_err_x

private

Definition at line 364 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixel_clst_err_y

private

Definition at line 365 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_pixel_split

private

Definition at line 363 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixgx[maxpix]

private

Definition at line 353 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixgy[maxpix]

private

Definition at line 354 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixgz[maxpix]

private

Definition at line 355 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixrow[maxpix]

private

Definition at line 347 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixx[maxpix]

private

Definition at line 351 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_pixy[maxpix]

private

Definition at line 352 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_plaq

private

Definition at line 251 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechiterrx

private

Definition at line 266 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechiterry

private

Definition at line 267 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitpullx

private

Definition at line 272 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitpully

private

Definition at line 273 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitresx

private

Definition at line 269 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitresy

private

Definition at line 270 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitx

private

Definition at line 259 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechity

private

Definition at line 260 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_rechitz

private

Definition at line 261 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_row1

private

Definition at line 313 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_row2

private

Definition at line 314 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_rows

private

Definition at line 257 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_side

private

Definition at line 247 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simhitx

private

Definition at line 263 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simhity

private

Definition at line 264 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simphi

private

Definition at line 288 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_simproc

private

Definition at line 293 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simpx

private

Definition at line 298 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simpy

private

Definition at line 299 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simpz

private

Definition at line 300 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simtheta

private

Definition at line 289 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simtrketa

private

Definition at line 325 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_simtrkphi

private

Definition at line 326 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_subdetid

private

Definition at line 241 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::all_trkid

private

Definition at line 304 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_vtxr

private

Definition at line 295 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_vtxz

private

Definition at line 296 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_x1

private

Definition at line 306 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_x2

private

Definition at line 307 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::all_y1

private

Definition at line 308 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_y2

private

Definition at line 309 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_z1

private

Definition at line 310 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::all_z2

private

Definition at line 311 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::alpha

private

Definition at line 195 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::beta

private

Definition at line 196 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::bigx

private

Definition at line 192 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::bigy

private

Definition at line 193 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::blade

private

Definition at line 210 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::charge

private

Definition at line 187 of file SiPixelErrorEstimation.h.

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

bool SiPixelErrorEstimation::checkType_

private

Definition at line 80 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

edm::ParameterSet SiPixelErrorEstimation::conf_

private

Definition at line 77 of file SiPixelErrorEstimation.h.

int SiPixelErrorEstimation::detector_type

private

Definition at line 114 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::disk

private

Definition at line 209 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::edgex

private

Definition at line 189 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::edgey

private

Definition at line 190 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::eta

private

Definition at line 202 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::evt

private

Definition at line 224 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::flipped

private

Definition at line 215 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::genType_

private

Definition at line 81 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

int SiPixelErrorEstimation::half

private

Definition at line 214 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::hit_cprob0

private

Definition at line 229 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_cprob1

private

Definition at line 230 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_cprob2

private

Definition at line 231 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_probx

private

Definition at line 227 of file SiPixelErrorEstimation.h.

Referenced by analyze().

float SiPixelErrorEstimation::hit_proby

private

Definition at line 228 of file SiPixelErrorEstimation.h.

Referenced by analyze().

bool SiPixelErrorEstimation::include_trk_hits_

private

Definition at line 82 of file SiPixelErrorEstimation.h.

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

int SiPixelErrorEstimation::ladder

private

Definition at line 206 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::layer

private

Definition at line 205 of file SiPixelErrorEstimation.h.

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

const int SiPixelErrorEstimation::maxpix = 10000

staticprivate

Definition at line 346 of file SiPixelErrorEstimation.h.

int SiPixelErrorEstimation::mod

private

Definition at line 207 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::npix

private

Definition at line 184 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::nsimhit

private

Definition at line 217 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::nxpix

private

Definition at line 185 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::nypix

private

Definition at line 186 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

std::string SiPixelErrorEstimation::outputFile_

private

Definition at line 78 of file SiPixelErrorEstimation.h.

Referenced by beginJob(), and SiPixelErrorEstimation().

int SiPixelErrorEstimation::panel

private

Definition at line 211 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::phi

private

Definition at line 201 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::pidhit

private

Definition at line 218 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::pixel_clst_err_x

private

Definition at line 235 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::pixel_clst_err_y

private

Definition at line 236 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::pixel_split

private

Definition at line 233 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::plaq

private

Definition at line 212 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechiterrx

private

Definition at line 88 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechiterry

private

Definition at line 89 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitpullx

private

Definition at line 92 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitpully

private

Definition at line 93 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitresx

private

Definition at line 90 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitresy

private

Definition at line 91 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitx

private

Definition at line 85 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechity

private

Definition at line 86 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::rechitz

private

Definition at line 87 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::run

private

Definition at line 225 of file SiPixelErrorEstimation.h.

Referenced by relval_steps.InputInfo::__str__(), generateEDF.LumiInfo::__str__(), analyze(), beginJob(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), relval_steps.InputInfo::dbs(), generateEDF.LumiInfo::fixXingInfo(), relval_steps.InputInfo::lumiRanges(), and dqm_interfaces.DirWalkerFile::walk().

int SiPixelErrorEstimation::side

private

Definition at line 208 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::simhitx

private

Definition at line 221 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::simhity

private

Definition at line 222 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::simproc

private

Definition at line 219 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

std::string SiPixelErrorEstimation::src_

private

Definition at line 79 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

float SiPixelErrorEstimation::strip_charge

private

Definition at line 122 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_clst_err_x

private

Definition at line 181 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_clst_err_y

private

Definition at line 182 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_cotalpha

private

Definition at line 117 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_cotbeta

private

Definition at line 118 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_edge

private

Definition at line 124 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_hit_type

private

Definition at line 113 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_is_stereo

private

Definition at line 112 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_locbx

private

Definition at line 119 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_locby

private

Definition at line 120 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_locbz

private

Definition at line 121 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_nprm

private

Definition at line 163 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_nsimhit

private

Definition at line 125 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit

private

Definition at line 126 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit1

private

Definition at line 165 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit2

private

Definition at line 168 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit3

private

Definition at line 171 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit4

private

Definition at line 174 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_pidhit5

private

Definition at line 177 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_prob

private

Definition at line 160 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_qbin

private

Definition at line 161 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechiterrx

private

Definition at line 102 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechiterry

private

Definition at line 103 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitpullx

private

Definition at line 110 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitpully

private

Definition at line 111 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitresx

private

Definition at line 104 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitresx2

private

Definition at line 107 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitresy

private

Definition at line 109 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitx

private

Definition at line 99 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechity

private

Definition at line 100 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_rechitz

private

Definition at line 101 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc

private

Definition at line 127 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc1

private

Definition at line 166 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc2

private

Definition at line 169 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc3

private

Definition at line 172 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc4

private

Definition at line 175 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_simproc5

private

Definition at line 178 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_size

private

Definition at line 123 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_split

private

Definition at line 180 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_subdet_id

private

Definition at line 129 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_double_side

private

Definition at line 135 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_external_string

private

Definition at line 142 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_internal_string

private

Definition at line 141 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_rphi

private

Definition at line 143 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_stereo

private

Definition at line 144 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_z_minus_side

private

Definition at line 137 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_is_z_plus_side

private

Definition at line 136 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_layer

private

Definition at line 131 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_layer_number

private

Definition at line 138 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_module

private

Definition at line 132 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_module_number

private

Definition at line 140 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_order

private

Definition at line 133 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_side

private

Definition at line 134 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tib_string_number

private

Definition at line 139 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_is_double_side

private

Definition at line 150 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_is_rphi

private

Definition at line 157 of file SiPixelErrorEstimation.h.

Referenced by analyze().

int SiPixelErrorEstimation::strip_tob_is_stereo

private

Definition at line 158 of file SiPixelErrorEstimation.h.

Referenced by analyze().

int SiPixelErrorEstimation::strip_tob_is_z_minus_side

private

Definition at line 152 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_is_z_plus_side

private

Definition at line 151 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_layer

private

Definition at line 146 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_layer_number

private

Definition at line 153 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_module

private

Definition at line 147 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_module_number

private

Definition at line 155 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_rod_number

private

Definition at line 154 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::strip_tob_side

private

Definition at line 149 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::strip_trk_pt

private

Definition at line 116 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

int SiPixelErrorEstimation::subdetId

private

Definition at line 204 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

TFile* SiPixelErrorEstimation::tfile_

private

Definition at line 370 of file SiPixelErrorEstimation.h.

Referenced by beginJob(), and endJob().

float SiPixelErrorEstimation::trk_alpha

private

Definition at line 198 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

float SiPixelErrorEstimation::trk_beta

private

Definition at line 199 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

TTree* SiPixelErrorEstimation::ttree_all_hits_

private

Definition at line 371 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

TTree* SiPixelErrorEstimation::ttree_track_hits_

private

Definition at line 372 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

TTree* SiPixelErrorEstimation::ttree_track_hits_strip_

private

Definition at line 374 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

Public Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation