#include <SiPixelErrorEstimation.h>

Inheritance diagram for SiPixelErrorEstimation:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

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

void	endJob () override

	SiPixelErrorEstimation (const edm::ParameterSet &)

	~SiPixelErrorEstimation () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magneticFieldToken_

int	mod

int	npix

int	nsimhit

int	nxpix

int	nypix

std::string	outputFile_

int	panel

float	phi

int	pidhit

float	pixel_clst_err_x

float	pixel_clst_err_y

int	pixel_split

int	plaq

float	rechiterrx

float	rechiterry

float	rechitpullx

float	rechitpully

float	rechitresx

float	rechitresy

float	rechitx

float	rechity

float	rechitz

int	run

int	side

float	simhitx

float	simhity

int	simproc

std::string	src_

float	strip_charge

float	strip_clst_err_x

float	strip_clst_err_y

float	strip_cotalpha

float	strip_cotbeta

int	strip_edge

int	strip_hit_type

int	strip_is_stereo

float	strip_locbx

float	strip_locby

float	strip_locbz

int	strip_nprm

int	strip_nsimhit

int	strip_pidhit

int	strip_pidhit1

int	strip_pidhit2

int	strip_pidhit3

int	strip_pidhit4

int	strip_pidhit5

float	strip_prob

int	strip_qbin

float	strip_rechiterrx

float	strip_rechiterry

float	strip_rechitpullx

float	strip_rechitpully

float	strip_rechitresx

float	strip_rechitresx2

float	strip_rechitresy

float	strip_rechitx

float	strip_rechity

float	strip_rechitz

int	strip_simproc

int	strip_simproc1

int	strip_simproc2

int	strip_simproc3

int	strip_simproc4

int	strip_simproc5

int	strip_size

int	strip_split

int	strip_subdet_id

int	strip_tib_is_double_side

int	strip_tib_is_external_string

int	strip_tib_is_internal_string

int	strip_tib_is_rphi

int	strip_tib_is_stereo

int	strip_tib_is_z_minus_side

int	strip_tib_is_z_plus_side

int	strip_tib_layer

int	strip_tib_layer_number

int	strip_tib_module

int	strip_tib_module_number

int	strip_tib_order

int	strip_tib_side

int	strip_tib_string_number

int	strip_tob_is_double_side

int	strip_tob_is_rphi

int	strip_tob_is_stereo

int	strip_tob_is_z_minus_side

int	strip_tob_is_z_plus_side

int	strip_tob_layer

int	strip_tob_layer_number

int	strip_tob_module

int	strip_tob_module_number

int	strip_tob_rod_number

int	strip_tob_side

float	strip_trk_pt

int	subdetId

TFile *	tfile_

edm::EDGetTokenT< SiPixelRecHitCollection >	tPixRecHitCollection

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomToken_

TrackerHitAssociator::Config	trackerHitAssociatorConfig_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopoToken_

float	trk_alpha

float	trk_beta

edm::EDGetTokenT< edm::SimTrackContainer >	tSimTrackContainer

edm::EDGetTokenT< reco::TrackCollection >	tTrackCollection

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

TTree *	ttree_all_hits_

TTree *	ttree_track_hits_

TTree *	ttree_track_hits_strip_

Static Private Attributes
static const int	maxpix = 10000

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 62 of file SiPixelErrorEstimation.h.

Constructor & Destructor Documentation

◆ SiPixelErrorEstimation()

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_, ProducerED_cfi::InputTag, magneticFieldToken_, outputFile_, src_, AlCaHLTBitMon_QueryRunRegistry::string, tPixRecHitCollection, trackerGeomToken_, trackerTopoToken_, tSimTrackContainer, tTrackCollection, and tTrajectory.

     : tfile_(nullptr),
       ttree_all_hits_(nullptr),
       ttree_track_hits_(nullptr),
       ttree_track_hits_strip_(nullptr),
       trackerHitAssociatorConfig_(consumesCollector()) {
   //Read config file
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "SiPixelErrorEstimation_Ntuple.root");
 
   // Replace  "ctfWithMaterialTracks" with "generalTracks"
   //src_ = ps.getUntrackedParameter<std::string>( "src", "ctfWithMaterialTracks" );
   src_ = ps.getUntrackedParameter<std::string>("src", "generalTracks");
 
   checkType_ = ps.getParameter<bool>("checkType");
   genType_ = ps.getParameter<int>("genType");
   include_trk_hits_ = ps.getParameter<bool>("include_trk_hits");
 
   tTrajectory = consumes<std::vector<Trajectory>>(src_);
   tPixRecHitCollection = consumes<SiPixelRecHitCollection>(edm::InputTag("siPixelRecHits"));
   tSimTrackContainer = consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
   tTrackCollection = consumes<reco::TrackCollection>(src_);
   trackerTopoToken_ = esConsumes<TrackerTopology, TrackerTopologyRcd>();
   trackerGeomToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
   magneticFieldToken_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
 }

◆ ~SiPixelErrorEstimation()

SiPixelErrorEstimation::~SiPixelErrorEstimation ( )

override

Definition at line 61 of file SiPixelErrorEstimation.cc.

61 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 374 of file SiPixelErrorEstimation.cc.

                                                                                {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle = es.getHandle(trackerTopoToken_);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   using namespace edm;
 
   run = e.id().run();
   evt = e.id().event();
 
   if (evt % 1000 == 0)
     cout << "evt = " << evt << endl;
 
   float math_pi = 3.14159265;
   float radtodeg = 180.0 / math_pi;
 
   LocalPoint position;
   LocalError error;
   float mindist = 999999.9;
 
   std::vector<PSimHit> matched;
   TrackerHitAssociator associate(e, trackerHitAssociatorConfig_);
 
   edm::ESHandle<TrackerGeometry> pDD = es.getHandle(trackerGeomToken_);
   const TrackerGeometry* tracker = &(*pDD);
 
   //cout << "...1..." << endl;
 
   edm::ESHandle<MagneticField> magneticField = es.getHandle(magneticFieldToken_);
   //const MagneticField* magField_ = magFieldHandle.product();
 
   edm::FileInPath FileInPath_;
 
   // Strip hits ==============================================================================================================
 
   edm::Handle<vector<Trajectory>> trajCollectionHandle;
 
   e.getByToken(tTrajectory, trajCollectionHandle);
   //e.getByLabel( "generalTracks", trajCollectionHandle);
 
   for (vector<Trajectory>::const_iterator it = trajCollectionHandle->begin(); it != trajCollectionHandle->end(); ++it) {
     vector<TrajectoryMeasurement> tmColl = it->measurements();
     for (vector<TrajectoryMeasurement>::const_iterator itTraj = tmColl.begin(); itTraj != tmColl.end(); ++itTraj) {
       if (!itTraj->updatedState().isValid())
         continue;
 
       strip_rechitx = -9999999.9;
       strip_rechity = -9999999.9;
       strip_rechitz = -9999999.9;
       strip_rechiterrx = -9999999.9;
       strip_rechiterry = -9999999.9;
       strip_rechitresx = -9999999.9;
       strip_rechitresx2 = -9999999.9;
 
       strip_rechitresy = -9999999.9;
       strip_rechitpullx = -9999999.9;
       strip_rechitpully = -9999999.9;
       strip_is_stereo = -9999999;
       strip_hit_type = -9999999;
       detector_type = -9999999;
 
       strip_trk_pt = -9999999.9;
       strip_cotalpha = -9999999.9;
       strip_cotbeta = -9999999.9;
       strip_locbx = -9999999.9;
       strip_locby = -9999999.9;
       strip_locbz = -9999999.9;
       strip_charge = -9999999.9;
       strip_size = -9999999;
 
       strip_edge = -9999999;
       strip_nsimhit = -9999999;
       strip_pidhit = -9999999;
       strip_simproc = -9999999;
 
       strip_subdet_id = -9999999;
 
       strip_tib_layer = -9999999;
       strip_tib_module = -9999999;
       strip_tib_order = -9999999;
       strip_tib_side = -9999999;
       strip_tib_is_double_side = -9999999;
       strip_tib_is_z_plus_side = -9999999;
       strip_tib_is_z_minus_side = -9999999;
       strip_tib_layer_number = -9999999;
       strip_tib_string_number = -9999999;
       strip_tib_module_number = -9999999;
       strip_tib_is_internal_string = -9999999;
       strip_tib_is_external_string = -9999999;
       strip_tib_is_rphi = -9999999;
       strip_tib_is_stereo = -9999999;
 
       strip_tob_layer = -9999999;
       strip_tob_module = -9999999;
       strip_tob_side = -9999999;
       strip_tob_is_double_side = -9999999;
       strip_tob_is_z_plus_side = -9999999;
       strip_tob_is_z_minus_side = -9999999;
       strip_tob_layer_number = -9999999;
       strip_tob_rod_number = -9999999;
       strip_tob_module_number = -9999999;
 
       strip_tob_is_rphi = -9999999;
       strip_tob_is_stereo = -9999999;
 
       strip_prob = -9999999.9;
       strip_qbin = -9999999;
 
       strip_nprm = -9999999;
 
       strip_pidhit1 = -9999999;
       strip_simproc1 = -9999999;
 
       strip_pidhit2 = -9999999;
       strip_simproc2 = -9999999;
 
       strip_pidhit3 = -9999999;
       strip_simproc3 = -9999999;
 
       strip_pidhit4 = -9999999;
       strip_simproc4 = -9999999;
 
       strip_pidhit5 = -9999999;
       strip_simproc5 = -9999999;
 
       strip_split = -9999999;
       strip_clst_err_x = -9999999.9;
       strip_clst_err_y = -9999999.9;
 
       const TransientTrackingRecHit::ConstRecHitPointer trans_trk_rec_hit_point = itTraj->recHit();
 
       if (trans_trk_rec_hit_point == nullptr)
         continue;
 
       const TrackingRecHit* trk_rec_hit = (*trans_trk_rec_hit_point).hit();
 
       if (trk_rec_hit == nullptr)
         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 != nullptr) {
         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() == SiStripModuleGeometry::IB1) {
         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() == SiStripModuleGeometry::IB2) {
         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() == SiStripModuleGeometry::OB1) {
         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() == SiStripModuleGeometry::OB2) {
         detector_type = 4;
         //cout << "si_strip_det_id.moduleGeometry() = OB2" << endl;
         //cout << "si_strip_det_id.moduleGeometry() = " << si_strip_det_id.moduleGeometry() << endl;
       }
 
       // Store ntuple variables
 
       if ((int)detid.subdetId() == int(StripSubdetector::TIB)) {
         strip_tib_layer = (int)tTopo->tibLayer(detid);
         strip_tib_module = (int)tTopo->tibModule(detid);
         strip_tib_order = (int)tTopo->tibOrder(detid);
         strip_tib_side = (int)tTopo->tibSide(detid);
         strip_tib_is_double_side = (int)tTopo->tibIsDoubleSide(detid);
         strip_tib_is_z_plus_side = (int)tTopo->tibIsZPlusSide(detid);
         strip_tib_is_z_minus_side = (int)tTopo->tibIsZMinusSide(detid);
         strip_tib_layer_number = (int)tTopo->tibLayer(detid);
         strip_tib_string_number = (int)tTopo->tibString(detid);
         strip_tib_module_number = (int)tTopo->tibModule(detid);
         strip_tib_is_internal_string = (int)tTopo->tibIsInternalString(detid);
         strip_tib_is_external_string = (int)tTopo->tibIsExternalString(detid);
         strip_tib_is_rphi = (int)tTopo->tibIsRPhi(detid);
         strip_tib_is_stereo = (int)tTopo->tibIsStereo(detid);
       }
 
       if ((int)detid.subdetId() == int(StripSubdetector::TOB)) {
         strip_tob_layer = (int)tTopo->tobLayer(detid);
         strip_tob_module = (int)tTopo->tobModule(detid);
 
         strip_tob_side = (int)tTopo->tobSide(detid);
         strip_tob_is_double_side = (int)tTopo->tobIsDoubleSide(detid);
         strip_tob_is_z_plus_side = (int)tTopo->tobIsZPlusSide(detid);
         strip_tob_is_z_minus_side = (int)tTopo->tobIsZMinusSide(detid);
         strip_tob_layer_number = (int)tTopo->tobLayer(detid);
         strip_tob_rod_number = (int)tTopo->tobRod(detid);
         strip_tob_module_number = (int)tTopo->tobModule(detid);
 
         strip_tob_is_rphi = (int)tTopo->tobIsRPhi(detid);
         strip_tob_is_stereo = (int)tTopo->tobIsStereo(detid);
       }
 
       if (matchedhit) {
         //cout << endl << endl << endl;
         //cout << "Found a SiStripMatchedRecHit2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! " << endl;
         //cout << endl << endl << endl;
         strip_hit_type = 0;
 
       }  //  if ( matchedhit )
 
       if (hit1d) {
         strip_hit_type = 1;
 
         const SiStripRecHit1D::ClusterRef& cluster = hit1d->cluster();
 
         if (cluster->getSplitClusterError() > 0.0)
           strip_split = 1;
         else
           strip_split = 0;
 
         strip_clst_err_x = cluster->getSplitClusterError();
         //strip_clst_err_y = ...
 
         // Get cluster total charge
         const auto& stripCharges = cluster->amplitudes();
         uint16_t charge = 0;
         for (unsigned int i = 0; i < stripCharges.size(); ++i) {
           charge += stripCharges[i];
         }
 
         strip_charge = (float)charge;
         strip_size = (int)((cluster->amplitudes()).size());
 
         // Association of the rechit to the simhit
         float mindist = 999999.9;
         matched.clear();
 
         matched = associate.associateHit(*hit1d);
 
         strip_nsimhit = (int)matched.size();
 
         if (!matched.empty()) {
           PSimHit closest;
 
           // Get the closest simhit
 
           int strip_nprimaries = 0;
           int current_index = 0;
 
           for (vector<PSimHit>::const_iterator m = matched.begin(); m < matched.end(); ++m) {
             ++current_index;
 
             if ((*m).processType() == 2)
               ++strip_nprimaries;
 
             if (current_index == 1) {
               strip_pidhit1 = (*m).particleType();
               strip_simproc1 = (*m).processType();
             } else if (current_index == 2) {
               strip_pidhit2 = (*m).particleType();
               strip_simproc2 = (*m).processType();
             } else if (current_index == 3) {
               strip_pidhit3 = (*m).particleType();
               strip_simproc3 = (*m).processType();
             } else if (current_index == 4) {
               strip_pidhit4 = (*m).particleType();
               strip_simproc4 = (*m).processType();
             } else if (current_index == 5) {
               strip_pidhit5 = (*m).particleType();
               strip_simproc5 = (*m).processType();
             }
 
             float dist = abs((hit1d)->localPosition().x() - (*m).localPosition().x());
 
             if (dist < mindist) {
               mindist = dist;
               closest = (*m);
             }
 
           }  // for ( vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); ++m)
 
           strip_nprm = strip_nprimaries;
 
           strip_rechitresx = strip_rechitx - closest.localPosition().x();
           strip_rechitresy = strip_rechity - closest.localPosition().y();
 
           strip_rechitresx2 = strip_rechitx - matched[0].localPosition().x();
 
           strip_rechitpullx = strip_rechitresx / strip_rechiterrx;
           strip_rechitpully = strip_rechitresy / strip_rechiterry;
 
           strip_pidhit = closest.particleType();
           strip_simproc = closest.processType();
 
         }  //   if( !matched.empty())
 
         //strip_prob = (hit1d)->getTemplProb();
         //strip_qbin = (hit1d)->getTemplQbin();
 
         //cout << endl;
         //cout << "SiPixelErrorEstimation 1d hit: " << endl;
         //cout << "prob 1d = " << strip_prob << endl;
         //cout << "qbin 1d = " << strip_qbin << endl;
         //cout << endl;
 
         // Is the cluster on edge ?
         /*
                 const auto stripDetInfo = SiStripDetInfoFileReader::read(edm::FileInPath{SiStripDetInfoFileReader::kDefaultFile}.fullPath());
         
         uint16_t firstStrip = cluster->firstStrip();
         uint16_t lastStrip = firstStrip + (cluster->amplitudes()).size() -1;
         unsigned short Nstrips;
         Nstrips = stripDetInfo.getNumberOfApvsAndStripLength(id1).first*128;
         
         if ( firstStrip == 0 || lastStrip == (Nstrips-1) ) 
         strip_edge = 1;
         else
         strip_edge = 0;
           */
 
       }  // if ( hit1d )
 
       if (hit2d) {
         strip_hit_type = 2;
 
         const SiStripRecHit1D::ClusterRef& cluster = hit2d->cluster();
 
         //if ( cluster->getSplitClusterError()  > 0.0 )
         //strip_split = 1;
         //else
         //strip_split = 0;
 
         //strip_clst_err_x = cluster->getSplitClusterError();
         //strip_clst_err_y = ...
 
         // Get cluster total charge
         auto& stripCharges = cluster->amplitudes();
         uint16_t charge = 0;
         for (unsigned int i = 0; i < stripCharges.size(); ++i) {
           charge += stripCharges[i];
         }
 
         strip_charge = (float)charge;
         strip_size = (int)((cluster->amplitudes()).size());
 
         // Association of the rechit to the simhit
         float mindist = 999999.9;
         matched.clear();
 
         matched = associate.associateHit(*hit2d);
 
         strip_nsimhit = (int)matched.size();
 
         if (!matched.empty()) {
           PSimHit closest;
 
           // Get the closest simhit
 
           for (vector<PSimHit>::const_iterator m = matched.begin(); m < matched.end(); ++m) {
             float dist = abs((hit2d)->localPosition().x() - (*m).localPosition().x());
 
             if (dist < mindist) {
               mindist = dist;
               closest = (*m);
             }
           }
 
           strip_rechitresx = strip_rechitx - closest.localPosition().x();
           strip_rechitresy = strip_rechity - closest.localPosition().y();
 
           strip_rechitpullx = strip_rechitresx / strip_rechiterrx;
           strip_rechitpully = strip_rechitresy / strip_rechiterry;
 
           strip_pidhit = closest.particleType();
           strip_simproc = closest.processType();
 
         }  //   if( !matched.empty())
 
         //strip_prob = (hit2d)->getTemplProb();
         //strip_qbin = (hit2d)->getTemplQbin();
 
         //cout << endl;
         //cout << "SiPixelErrorEstimation 2d hit: " << endl;
         //cout << "prob 2d = " << strip_prob << endl;
         //cout << "qbin 2d = " << strip_qbin << endl;
         //cout << endl;
 
       }  // if ( hit2d )
 
       ttree_track_hits_strip_->Fill();
 
     }  // for ( vector<TrajectoryMeasurement>::const_iterator itTraj = tmColl.begin(); itTraj!=tmColl.end(); ++itTraj )
 
   }  // for( vector<Trajectory>::const_iterator it = trajCollectionHandle->begin(); it!=trajCollectionHandle->end(); ++it)
 
   //cout << "...2..." << endl;
 
   // --------------------------------------- all hits -----------------------------------------------------------
   edm::Handle<SiPixelRecHitCollection> recHitColl;
   e.getByToken(tPixRecHitCollection, recHitColl);
 
   Handle<edm::SimTrackContainer> simtracks;
   e.getByToken(tSimTrackContainer, simtracks);
 
   //-----Iterate over detunits
   for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++) {
     DetId detId = ((*it)->geographicalId());
 
     SiPixelRecHitCollection::const_iterator dsmatch = recHitColl->find(detId);
     if (dsmatch == recHitColl->end())
       continue;
 
     SiPixelRecHitCollection::DetSet pixelrechitRange = *dsmatch;
     SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
     SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
     SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
     std::vector<PSimHit> matched;
 
     //----Loop over rechits for this detId
     for (; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter) {
       matched.clear();
       matched = associate.associateHit(*pixeliter);
       // only consider rechits that have associated simhit
       // otherwise cannot determine residiual
       if (matched.empty()) {
         cout << "SiPixelErrorEstimation::analyze: rechits without associated simhit !!!!!!!" << endl;
         continue;
       }
 
       all_subdetid = -9999;
 
       all_layer = -9999;
       all_ladder = -9999;
       all_mod = -9999;
 
       all_side = -9999;
       all_disk = -9999;
       all_blade = -9999;
       all_panel = -9999;
       all_plaq = -9999;
 
       all_half = -9999;
       all_flipped = -9999;
 
       all_cols = -9999;
       all_rows = -9999;
 
       all_rechitx = -9999;
       all_rechity = -9999;
       all_rechitz = -9999;
 
       all_simhitx = -9999;
       all_simhity = -9999;
 
       all_rechiterrx = -9999;
       all_rechiterry = -9999;
 
       all_rechitresx = -9999;
       all_rechitresy = -9999;
 
       all_rechitpullx = -9999;
       all_rechitpully = -9999;
 
       all_npix = -9999;
       all_nxpix = -9999;
       all_nypix = -9999;
 
       all_edgex = -9999;
       all_edgey = -9999;
 
       all_bigx = -9999;
       all_bigy = -9999;
 
       all_alpha = -9999;
       all_beta = -9999;
 
       all_simphi = -9999;
       all_simtheta = -9999;
 
       all_simhitx = -9999;
       all_simhity = -9999;
 
       all_nsimhit = -9999;
       all_pidhit = -9999;
       all_simproc = -9999;
 
       all_vtxr = -9999;
       all_vtxz = -9999;
 
       all_simpx = -9999;
       all_simpy = -9999;
       all_simpz = -9999;
 
       all_eloss = -9999;
 
       all_trkid = -9999;
 
       all_x1 = -9999;
       all_x2 = -9999;
       all_y1 = -9999;
       all_y2 = -9999;
       all_z1 = -9999;
       all_z2 = -9999;
 
       all_row1 = -9999;
       all_row2 = -9999;
       all_col1 = -9999;
       all_col2 = -9999;
 
       all_gx1 = -9999;
       all_gx2 = -9999;
       all_gy1 = -9999;
       all_gy2 = -9999;
       all_gz1 = -9999;
       all_gz2 = -9999;
 
       all_simtrketa = -9999;
       all_simtrkphi = -9999;
 
       all_clust_row = -9999;
       all_clust_col = -9999;
 
       all_clust_x = -9999;
       all_clust_y = -9999;
 
       all_clust_q = -9999;
 
       all_clust_maxpixcol = -9999;
       all_clust_maxpixrow = -9999;
       all_clust_minpixcol = -9999;
       all_clust_minpixrow = -9999;
 
       all_clust_geoid = -9999;
 
       all_clust_alpha = -9999;
       all_clust_beta = -9999;
 
       /*
         for (int i=0; i<all_npix; ++i)
         {
         all_pixrow[i] = -9999;
         all_pixcol[i] = -9999;
         all_pixadc[i] = -9999;
         all_pixx[i] = -9999;
         all_pixy[i] = -9999;
         all_pixgx[i] = -9999;
         all_pixgy[i] = -9999;
         all_pixgz[i] = -9999;
         }
       */
 
       all_hit_probx = -9999;
       all_hit_proby = -9999;
       all_hit_cprob0 = -9999;
       all_hit_cprob1 = -9999;
       all_hit_cprob2 = -9999;
 
       all_pixel_split = -9999;
       all_pixel_clst_err_x = -9999.9;
       all_pixel_clst_err_y = -9999.9;
 
       all_nsimhit = (int)matched.size();
 
       all_subdetid = (int)detId.subdetId();
       // only consider rechits in pixel barrel and pixel forward
       if (!(all_subdetid == 1 || all_subdetid == 2)) {
         cout << "SiPixelErrorEstimation::analyze: Not in a pixel detector !!!!!" << endl;
         continue;
       }
 
       const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(tracker->idToDet(detId));
 
       const PixelTopology* topol = &(theGeomDet->specificTopology());
 
       if (pixeliter->cluster()->getSplitClusterErrorX() > 0.0 && pixeliter->cluster()->getSplitClusterErrorY() > 0.0) {
         all_pixel_split = 1;
       } else {
         all_pixel_split = 0;
       }
 
       all_pixel_clst_err_x = pixeliter->cluster()->getSplitClusterErrorX();
       all_pixel_clst_err_y = pixeliter->cluster()->getSplitClusterErrorY();
 
       const int maxPixelCol = pixeliter->cluster()->maxPixelCol();
       const int maxPixelRow = pixeliter->cluster()->maxPixelRow();
       const int minPixelCol = pixeliter->cluster()->minPixelCol();
       const int minPixelRow = pixeliter->cluster()->minPixelRow();
 
       //all_hit_probx  = (float)pixeliter->probabilityX();
       //all_hit_proby  = (float)pixeliter->probabilityY();
       all_hit_cprob0 = (float)pixeliter->clusterProbability(0);
       all_hit_cprob1 = (float)pixeliter->clusterProbability(1);
       all_hit_cprob2 = (float)pixeliter->clusterProbability(2);
 
       // check whether the cluster is at the module edge
       if (topol->isItEdgePixelInX(minPixelRow) || topol->isItEdgePixelInX(maxPixelRow))
         all_edgex = 1;
       else
         all_edgex = 0;
 
       if (topol->isItEdgePixelInY(minPixelCol) || topol->isItEdgePixelInY(maxPixelCol))
         all_edgey = 1;
       else
         all_edgey = 0;
 
       // check whether this rechit contains big pixels
       if (topol->containsBigPixelInX(minPixelRow, maxPixelRow))
         all_bigx = 1;
       else
         all_bigx = 0;
 
       if (topol->containsBigPixelInY(minPixelCol, maxPixelCol))
         all_bigy = 1;
       else
         all_bigy = 0;
 
       if ((int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel) {
         all_layer = tTopo->pxbLayer(detId);
         all_ladder = tTopo->pxbLadder(detId);
         all_mod = tTopo->pxbModule(detId);
 
         int tmp_nrows = theGeomDet->specificTopology().nrows();
         if (tmp_nrows == 80)
           all_half = 1;
         else if (tmp_nrows == 160)
           all_half = 0;
         else
           cout << "-------------------------------------------------- Wrong module size !!!" << endl;
 
         float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 0.)).perp();
         float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 1.)).perp();
 
         if (tmp2 < tmp1)
           all_flipped = 1;
         else
           all_flipped = 0;
       } else if ((int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap) {
         all_side = tTopo->pxfSide(detId);
         all_disk = tTopo->pxfDisk(detId);
         all_blade = tTopo->pxfBlade(detId);
         all_panel = tTopo->pxfPanel(detId);
         all_plaq = tTopo->pxfModule(detId);  // also known as plaquette
 
       }  // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
       else
         std::cout << "We are not in the pixel detector" << (int)detId.subdetId() << endl;
 
       all_cols = theGeomDet->specificTopology().ncolumns();
       all_rows = theGeomDet->specificTopology().nrows();
 
       LocalPoint lp = pixeliter->localPosition();
       // gavril: change this name
       all_rechitx = lp.x();
       all_rechity = lp.y();
       all_rechitz = lp.z();
 
       LocalError le = pixeliter->localPositionError();
       all_rechiterrx = sqrt(le.xx());
       all_rechiterry = sqrt(le.yy());
 
       bool found_hit_from_generated_particle = false;
 
       //---Loop over sim hits, fill closest
       float closest_dist = 99999.9;
       std::vector<PSimHit>::const_iterator closest_simhit = matched.begin();
 
       for (std::vector<PSimHit>::const_iterator m = matched.begin(); m < matched.end(); m++) {
         if (checkType_) {
           int pid = (*m).particleType();
           if (abs(pid) != genType_)
             continue;
         }
 
         float simhitx = 0.5 * ((*m).entryPoint().x() + (*m).exitPoint().x());
         float simhity = 0.5 * ((*m).entryPoint().y() + (*m).exitPoint().y());
 
         float x_res = simhitx - rechitx;
         float y_res = simhity - rechity;
 
         float dist = sqrt(x_res * x_res + y_res * y_res);
 
         if (dist < closest_dist) {
           closest_dist = dist;
           closest_simhit = m;
           found_hit_from_generated_particle = true;
         }
       }  // end sim hit loop
 
       // If this recHit does not have any simHit with the same particleType as the particles generated
       // ignore it as most probably comes from delta rays.
       if (checkType_ && !found_hit_from_generated_particle)
         continue;
 
       all_x1 = (*closest_simhit).entryPoint().x();  // width (row index, in col direction)
       all_y1 = (*closest_simhit).entryPoint().y();  // length (col index, in row direction)
       all_z1 = (*closest_simhit).entryPoint().z();
       all_x2 = (*closest_simhit).exitPoint().x();
       all_y2 = (*closest_simhit).exitPoint().y();
       all_z2 = (*closest_simhit).exitPoint().z();
       GlobalPoint GP1 = theGeomDet->surface().toGlobal(Local3DPoint(
           (*closest_simhit).entryPoint().x(), (*closest_simhit).entryPoint().y(), (*closest_simhit).entryPoint().z()));
       GlobalPoint GP2 = theGeomDet->surface().toGlobal(Local3DPoint(
           (*closest_simhit).exitPoint().x(), (*closest_simhit).exitPoint().y(), (*closest_simhit).exitPoint().z()));
       all_gx1 = GP1.x();
       all_gx2 = GP2.x();
       all_gy1 = GP1.y();
       all_gy2 = GP2.y();
       all_gz1 = GP1.z();
       all_gz2 = GP2.z();
 
       MeasurementPoint mp1 = topol->measurementPosition(LocalPoint(
           (*closest_simhit).entryPoint().x(), (*closest_simhit).entryPoint().y(), (*closest_simhit).entryPoint().z()));
       MeasurementPoint mp2 = topol->measurementPosition(LocalPoint(
           (*closest_simhit).exitPoint().x(), (*closest_simhit).exitPoint().y(), (*closest_simhit).exitPoint().z()));
       all_row1 = mp1.x();
       all_col1 = mp1.y();
       all_row2 = mp2.x();
       all_col2 = mp2.y();
 
       all_simhitx = 0.5 * (all_x1 + all_x2);
       all_simhity = 0.5 * (all_y1 + all_y2);
 
       all_rechitresx = all_rechitx - all_simhitx;
       all_rechitresy = all_rechity - all_simhity;
 
       all_rechitpullx = all_rechitresx / all_rechiterrx;
       all_rechitpully = all_rechitresy / all_rechiterry;
 
       SiPixelRecHit::ClusterRef const& clust = pixeliter->cluster();
 
       all_npix = clust->size();
       all_nxpix = clust->sizeX();
       all_nypix = clust->sizeY();
 
       all_clust_row = clust->x();
       all_clust_col = clust->y();
 
       LocalPoint lp2 = topol->localPosition(MeasurementPoint(all_clust_row, all_clust_col));
       all_clust_x = lp2.x();
       all_clust_y = lp2.y();
 
       all_clust_q = clust->charge();
 
       all_clust_maxpixcol = clust->maxPixelCol();
       all_clust_maxpixrow = clust->maxPixelRow();
       all_clust_minpixcol = clust->minPixelCol();
       all_clust_minpixrow = clust->minPixelRow();
 
       all_clust_geoid = 0;  // never set!
 
       all_simpx = (*closest_simhit).momentumAtEntry().x();
       all_simpy = (*closest_simhit).momentumAtEntry().y();
       all_simpz = (*closest_simhit).momentumAtEntry().z();
       all_eloss = (*closest_simhit).energyLoss();
       all_simphi = (*closest_simhit).phiAtEntry();
       all_simtheta = (*closest_simhit).thetaAtEntry();
       all_pidhit = (*closest_simhit).particleType();
       all_trkid = (*closest_simhit).trackId();
 
       //--- Fill alpha and beta -- more useful for exploring the residuals...
       all_beta = atan2(all_simpz, all_simpy);
       all_alpha = atan2(all_simpz, all_simpx);
 
       all_simproc = (int)closest_simhit->processType();
 
       const edm::SimTrackContainer& trks = *(simtracks.product());
       SimTrackContainer::const_iterator trksiter;
       for (trksiter = trks.begin(); trksiter != trks.end(); trksiter++)
         if ((int)trksiter->trackId() == all_trkid) {
           all_simtrketa = trksiter->momentum().eta();
           all_simtrkphi = trksiter->momentum().phi();
         }
 
       all_vtxz = theGeomDet->surface().position().z();
       all_vtxr = theGeomDet->surface().position().perp();
 
       //computeAnglesFromDetPosition(clust,
       //               theGeomDet,
       //               all_clust_alpha, all_clust_beta )
 
       const std::vector<SiPixelCluster::Pixel>& pixvector = clust->pixels();
       for (int i = 0; i < (int)pixvector.size(); ++i) {
         SiPixelCluster::Pixel holdpix = pixvector[i];
         all_pixrow[i] = holdpix.x;
         all_pixcol[i] = holdpix.y;
         all_pixadc[i] = holdpix.adc;
         LocalPoint lp = topol->localPosition(MeasurementPoint(holdpix.x, holdpix.y));
         all_pixx[i] = lp.x();
         all_pixy[i] = lp.y();
         GlobalPoint GP = theGeomDet->surface().toGlobal(Local3DPoint(lp.x(), lp.y(), lp.z()));
         all_pixgx[i] = GP.x();
         all_pixgy[i] = GP.y();
         all_pixgz[i] = GP.z();
       }
 
       ttree_all_hits_->Fill();
 
     }  // for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter)
 
   }  // for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++)
 
   // ------------------------------------------------ all hits ---------------------------------------------------------------
 
   //cout << "...3..." << endl;
 
   // ------------------------------------------------ track hits only --------------------------------------------------------
 
   if (include_trk_hits_) {
     // Get tracks
     edm::Handle<reco::TrackCollection> trackCollection;
     e.getByToken(tTrackCollection, trackCollection);
     const reco::TrackCollection* tracks = trackCollection.product();
     reco::TrackCollection::const_iterator tciter;
 
     if (!tracks->empty()) {
       // Loop on tracks
       for (tciter = tracks->begin(); tciter != tracks->end(); ++tciter) {
         // First loop on hits: find matched hits
         for (auto const hit : tciter->recHits()) {
           // Is it a matched hit?
           const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(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 = hit->localPosition();
             error = hit->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;
 
               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) {
                   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;
 
               DetId detId = hit->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 -----------------------------------------------------------
 }

◆ beginJob()

void SiPixelErrorEstimation::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 63 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);
 }

◆ computeAnglesFromDetPosition()

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

Definition at line 1521 of file SiPixelErrorEstimation.cc.

References alpha, cms::cuda::assert(), beta, haddnano::cl, gather_cfg::cout, runTauDisplay::gp, Topology::localPosition(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::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;
 }

◆ endJob()

void SiPixelErrorEstimation::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 369 of file SiPixelErrorEstimation.cc.

References tfile_.

Referenced by o2olib.O2ORunMgr::executeJob().

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

Member Data Documentation

◆ all_alpha

float SiPixelErrorEstimation::all_alpha

private

Definition at line 283 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_beta

float SiPixelErrorEstimation::all_beta

private

Definition at line 284 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_bigx

int SiPixelErrorEstimation::all_bigx

private

Definition at line 280 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_bigy

int SiPixelErrorEstimation::all_bigy

private

Definition at line 281 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_blade

int SiPixelErrorEstimation::all_blade

private

Definition at line 247 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_alpha

float SiPixelErrorEstimation::all_clust_alpha

private

Definition at line 341 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_beta

float SiPixelErrorEstimation::all_clust_beta

private

Definition at line 342 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_col

float SiPixelErrorEstimation::all_clust_col

private

Definition at line 327 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_geoid

int SiPixelErrorEstimation::all_clust_geoid

private

Definition at line 339 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_maxpixcol

int SiPixelErrorEstimation::all_clust_maxpixcol

private

Definition at line 334 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_maxpixrow

int SiPixelErrorEstimation::all_clust_maxpixrow

private

Definition at line 335 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_minpixcol

int SiPixelErrorEstimation::all_clust_minpixcol

private

Definition at line 336 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_minpixrow

int SiPixelErrorEstimation::all_clust_minpixrow

private

Definition at line 337 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_q

float SiPixelErrorEstimation::all_clust_q

private

Definition at line 332 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_row

float SiPixelErrorEstimation::all_clust_row

private

Definition at line 326 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_x

float SiPixelErrorEstimation::all_clust_x

private

Definition at line 329 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_clust_y

float SiPixelErrorEstimation::all_clust_y

private

Definition at line 330 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_col1

float SiPixelErrorEstimation::all_col1

private

Definition at line 313 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_col2

float SiPixelErrorEstimation::all_col2

private

Definition at line 314 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_cols

int SiPixelErrorEstimation::all_cols

private

Definition at line 254 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_disk

int SiPixelErrorEstimation::all_disk

private

Definition at line 246 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_edgex

int SiPixelErrorEstimation::all_edgex

private

Definition at line 277 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_edgey

int SiPixelErrorEstimation::all_edgey

private

Definition at line 278 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_eloss

float SiPixelErrorEstimation::all_eloss

private

Definition at line 300 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_flipped

int SiPixelErrorEstimation::all_flipped

private

Definition at line 252 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_gx1

float SiPixelErrorEstimation::all_gx1

private

Definition at line 316 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_gx2

float SiPixelErrorEstimation::all_gx2

private

Definition at line 317 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_gy1

float SiPixelErrorEstimation::all_gy1

private

Definition at line 318 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_gy2

float SiPixelErrorEstimation::all_gy2

private

Definition at line 319 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_gz1

float SiPixelErrorEstimation::all_gz1

private

Definition at line 320 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_gz2

float SiPixelErrorEstimation::all_gz2

private

Definition at line 321 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_half

int SiPixelErrorEstimation::all_half

private

Definition at line 251 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_hit_cprob0

float SiPixelErrorEstimation::all_hit_cprob0

private

Definition at line 357 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_hit_cprob1

float SiPixelErrorEstimation::all_hit_cprob1

private

Definition at line 358 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_hit_cprob2

float SiPixelErrorEstimation::all_hit_cprob2

private

Definition at line 359 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_hit_probx

float SiPixelErrorEstimation::all_hit_probx

private

Definition at line 355 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_hit_proby

float SiPixelErrorEstimation::all_hit_proby

private

Definition at line 356 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_ladder

int SiPixelErrorEstimation::all_ladder

private

Definition at line 242 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_layer

int SiPixelErrorEstimation::all_layer

private

Definition at line 241 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_mod

int SiPixelErrorEstimation::all_mod

private

Definition at line 243 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_npix

int SiPixelErrorEstimation::all_npix

private

Definition at line 273 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_nsimhit

int SiPixelErrorEstimation::all_nsimhit

private

Definition at line 289 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_nxpix

int SiPixelErrorEstimation::all_nxpix

private

Definition at line 274 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_nypix

int SiPixelErrorEstimation::all_nypix

private

Definition at line 275 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_panel

int SiPixelErrorEstimation::all_panel

private

Definition at line 248 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pidhit

int SiPixelErrorEstimation::all_pidhit

private

Definition at line 290 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixadc

float SiPixelErrorEstimation::all_pixadc[maxpix]

private

Definition at line 347 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixcol

float SiPixelErrorEstimation::all_pixcol[maxpix]

private

Definition at line 346 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixel_clst_err_x

float SiPixelErrorEstimation::all_pixel_clst_err_x

private

Definition at line 362 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixel_clst_err_y

float SiPixelErrorEstimation::all_pixel_clst_err_y

private

Definition at line 363 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixel_split

int SiPixelErrorEstimation::all_pixel_split

private

Definition at line 361 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixgx

float SiPixelErrorEstimation::all_pixgx[maxpix]

private

Definition at line 351 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixgy

float SiPixelErrorEstimation::all_pixgy[maxpix]

private

Definition at line 352 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixgz

float SiPixelErrorEstimation::all_pixgz[maxpix]

private

Definition at line 353 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixrow

float SiPixelErrorEstimation::all_pixrow[maxpix]

private

Definition at line 345 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixx

float SiPixelErrorEstimation::all_pixx[maxpix]

private

Definition at line 349 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_pixy

float SiPixelErrorEstimation::all_pixy[maxpix]

private

Definition at line 350 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_plaq

int SiPixelErrorEstimation::all_plaq

private

Definition at line 249 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechiterrx

float SiPixelErrorEstimation::all_rechiterrx

private

Definition at line 264 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechiterry

float SiPixelErrorEstimation::all_rechiterry

private

Definition at line 265 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechitpullx

float SiPixelErrorEstimation::all_rechitpullx

private

Definition at line 270 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechitpully

float SiPixelErrorEstimation::all_rechitpully

private

Definition at line 271 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechitresx

float SiPixelErrorEstimation::all_rechitresx

private

Definition at line 267 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechitresy

float SiPixelErrorEstimation::all_rechitresy

private

Definition at line 268 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechitx

float SiPixelErrorEstimation::all_rechitx

private

Definition at line 257 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechity

float SiPixelErrorEstimation::all_rechity

private

Definition at line 258 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rechitz

float SiPixelErrorEstimation::all_rechitz

private

Definition at line 259 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_row1

float SiPixelErrorEstimation::all_row1

private

Definition at line 311 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_row2

float SiPixelErrorEstimation::all_row2

private

Definition at line 312 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_rows

int SiPixelErrorEstimation::all_rows

private

Definition at line 255 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_side

int SiPixelErrorEstimation::all_side

private

Definition at line 245 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simhitx

float SiPixelErrorEstimation::all_simhitx

private

Definition at line 261 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simhity

float SiPixelErrorEstimation::all_simhity

private

Definition at line 262 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simphi

float SiPixelErrorEstimation::all_simphi

private

Definition at line 286 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simproc

int SiPixelErrorEstimation::all_simproc

private

Definition at line 291 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simpx

float SiPixelErrorEstimation::all_simpx

private

Definition at line 296 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simpy

float SiPixelErrorEstimation::all_simpy

private

Definition at line 297 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simpz

float SiPixelErrorEstimation::all_simpz

private

Definition at line 298 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simtheta

float SiPixelErrorEstimation::all_simtheta

private

Definition at line 287 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simtrketa

float SiPixelErrorEstimation::all_simtrketa

private

Definition at line 323 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_simtrkphi

float SiPixelErrorEstimation::all_simtrkphi

private

Definition at line 324 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_subdetid

int SiPixelErrorEstimation::all_subdetid

private

Definition at line 239 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_trkid

int SiPixelErrorEstimation::all_trkid

private

Definition at line 302 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_vtxr

float SiPixelErrorEstimation::all_vtxr

private

Definition at line 293 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_vtxz

float SiPixelErrorEstimation::all_vtxz

private

Definition at line 294 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_x1

float SiPixelErrorEstimation::all_x1

private

Definition at line 304 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_x2

float SiPixelErrorEstimation::all_x2

private

Definition at line 305 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ all_y1

float SiPixelErrorEstimation::all_y1

private

Definition at line 306 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_y2

float SiPixelErrorEstimation::all_y2

private

Definition at line 307 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_z1

float SiPixelErrorEstimation::all_z1

private

Definition at line 308 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ all_z2

float SiPixelErrorEstimation::all_z2

private

Definition at line 309 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ alpha

float SiPixelErrorEstimation::alpha

private

Definition at line 194 of file SiPixelErrorEstimation.h.

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

◆ beta

float SiPixelErrorEstimation::beta

private

Definition at line 195 of file SiPixelErrorEstimation.h.

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

◆ bigx

int SiPixelErrorEstimation::bigx

private

Definition at line 191 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ bigy

int SiPixelErrorEstimation::bigy

private

Definition at line 192 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ blade

int SiPixelErrorEstimation::blade

private

Definition at line 209 of file SiPixelErrorEstimation.h.

Referenced by PixelMapPlotter.Forward::__str__(), analyze(), beginJob(), PixelMapPlotter.Forward::convertParts(), and PixelMapPlotter.Forward::GetXYCoords().

◆ charge

float SiPixelErrorEstimation::charge

private

Definition at line 186 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ checkType_

bool SiPixelErrorEstimation::checkType_

private

Definition at line 84 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ conf_

edm::ParameterSet SiPixelErrorEstimation::conf_

private

Definition at line 74 of file SiPixelErrorEstimation.h.

◆ detector_type

int SiPixelErrorEstimation::detector_type

private

Definition at line 113 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ disk

int SiPixelErrorEstimation::disk

private

Definition at line 208 of file SiPixelErrorEstimation.h.

Referenced by PixelMapPlotter.Forward::__str__(), analyze(), beginJob(), PixelMapPlotter.Forward::convertParts(), and PixelMapPlotter.Forward::GetXYCoords().

◆ edgex

int SiPixelErrorEstimation::edgex

private

Definition at line 188 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ edgey

int SiPixelErrorEstimation::edgey

private

Definition at line 189 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ eta

float SiPixelErrorEstimation::eta

private

Definition at line 201 of file SiPixelErrorEstimation.h.

Referenced by Particle.Particle::__str__(), analyze(), beginJob(), datamodel.Object::DeltaR(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

◆ evt

int SiPixelErrorEstimation::evt

private

Definition at line 223 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ flipped

int SiPixelErrorEstimation::flipped

private

Definition at line 214 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ genType_

int SiPixelErrorEstimation::genType_

private

Definition at line 85 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ half

int SiPixelErrorEstimation::half

private

Definition at line 213 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ hit_cprob0

float SiPixelErrorEstimation::hit_cprob0

private

Definition at line 228 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ hit_cprob1

float SiPixelErrorEstimation::hit_cprob1

private

Definition at line 229 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ hit_cprob2

float SiPixelErrorEstimation::hit_cprob2

private

Definition at line 230 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ hit_probx

float SiPixelErrorEstimation::hit_probx

private

Definition at line 226 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ hit_proby

float SiPixelErrorEstimation::hit_proby

private

Definition at line 227 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ include_trk_hits_

bool SiPixelErrorEstimation::include_trk_hits_

private

Definition at line 86 of file SiPixelErrorEstimation.h.

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

◆ ladder

int SiPixelErrorEstimation::ladder

private

Definition at line 205 of file SiPixelErrorEstimation.h.

Referenced by PixelMapPlotter.Barrel::__str__(), analyze(), beginJob(), PixelMapPlotter.Barrel::convertParts(), and PixelMapPlotter.Barrel::GetXYCoords().

◆ layer

int SiPixelErrorEstimation::layer

private

Definition at line 204 of file SiPixelErrorEstimation.h.

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

◆ magneticFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SiPixelErrorEstimation::magneticFieldToken_

private

Definition at line 81 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ maxpix

const int SiPixelErrorEstimation::maxpix = 10000

staticprivate

Definition at line 344 of file SiPixelErrorEstimation.h.

◆ mod

int SiPixelErrorEstimation::mod

private

Definition at line 206 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ npix

int SiPixelErrorEstimation::npix

private

Definition at line 183 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ nsimhit

int SiPixelErrorEstimation::nsimhit

private

Definition at line 216 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ nxpix

int SiPixelErrorEstimation::nxpix

private

Definition at line 184 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ nypix

int SiPixelErrorEstimation::nypix

private

Definition at line 185 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ outputFile_

std::string SiPixelErrorEstimation::outputFile_

private

Definition at line 82 of file SiPixelErrorEstimation.h.

Referenced by beginJob(), and SiPixelErrorEstimation().

◆ panel

int SiPixelErrorEstimation::panel

private

Definition at line 210 of file SiPixelErrorEstimation.h.

Referenced by PixelMapPlotter.Forward::__str__(), analyze(), beginJob(), and PixelMapPlotter.Forward::GetXYCoords().

◆ phi

float SiPixelErrorEstimation::phi

private

Definition at line 200 of file SiPixelErrorEstimation.h.

Referenced by Particle.Particle::__str__(), analyze(), beginJob(), datamodel.Object::DeltaR(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

◆ pidhit

int SiPixelErrorEstimation::pidhit

private

Definition at line 217 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ pixel_clst_err_x

float SiPixelErrorEstimation::pixel_clst_err_x

private

Definition at line 234 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ pixel_clst_err_y

float SiPixelErrorEstimation::pixel_clst_err_y

private

Definition at line 235 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ pixel_split

int SiPixelErrorEstimation::pixel_split

private

Definition at line 232 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ plaq

int SiPixelErrorEstimation::plaq

private

Definition at line 211 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechiterrx

float SiPixelErrorEstimation::rechiterrx

private

Definition at line 92 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechiterry

float SiPixelErrorEstimation::rechiterry

private

Definition at line 93 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechitpullx

float SiPixelErrorEstimation::rechitpullx

private

Definition at line 96 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechitpully

float SiPixelErrorEstimation::rechitpully

private

Definition at line 97 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechitresx

float SiPixelErrorEstimation::rechitresx

private

Definition at line 94 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechitresy

float SiPixelErrorEstimation::rechitresy

private

Definition at line 95 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechitx

float SiPixelErrorEstimation::rechitx

private

Definition at line 89 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechity

float SiPixelErrorEstimation::rechity

private

Definition at line 90 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ rechitz

float SiPixelErrorEstimation::rechitz

private

Definition at line 91 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ run

int SiPixelErrorEstimation::run

private

Definition at line 224 of file SiPixelErrorEstimation.h.

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

◆ side

int SiPixelErrorEstimation::side

private

Definition at line 207 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ simhitx

float SiPixelErrorEstimation::simhitx

private

Definition at line 220 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ simhity

float SiPixelErrorEstimation::simhity

private

Definition at line 221 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ simproc

int SiPixelErrorEstimation::simproc

private

Definition at line 218 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ src_

std::string SiPixelErrorEstimation::src_

private

Definition at line 83 of file SiPixelErrorEstimation.h.

Referenced by SiPixelErrorEstimation().

◆ strip_charge

float SiPixelErrorEstimation::strip_charge

private

Definition at line 121 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_clst_err_x

float SiPixelErrorEstimation::strip_clst_err_x

private

Definition at line 180 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_clst_err_y

float SiPixelErrorEstimation::strip_clst_err_y

private

Definition at line 181 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_cotalpha

float SiPixelErrorEstimation::strip_cotalpha

private

Definition at line 116 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_cotbeta

float SiPixelErrorEstimation::strip_cotbeta

private

Definition at line 117 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_edge

int SiPixelErrorEstimation::strip_edge

private

Definition at line 123 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_hit_type

int SiPixelErrorEstimation::strip_hit_type

private

Definition at line 112 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_is_stereo

int SiPixelErrorEstimation::strip_is_stereo

private

Definition at line 111 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_locbx

float SiPixelErrorEstimation::strip_locbx

private

Definition at line 118 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_locby

float SiPixelErrorEstimation::strip_locby

private

Definition at line 119 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_locbz

float SiPixelErrorEstimation::strip_locbz

private

Definition at line 120 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_nprm

int SiPixelErrorEstimation::strip_nprm

private

Definition at line 162 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_nsimhit

int SiPixelErrorEstimation::strip_nsimhit

private

Definition at line 124 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_pidhit

int SiPixelErrorEstimation::strip_pidhit

private

Definition at line 125 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_pidhit1

int SiPixelErrorEstimation::strip_pidhit1

private

Definition at line 164 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_pidhit2

int SiPixelErrorEstimation::strip_pidhit2

private

Definition at line 167 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_pidhit3

int SiPixelErrorEstimation::strip_pidhit3

private

Definition at line 170 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_pidhit4

int SiPixelErrorEstimation::strip_pidhit4

private

Definition at line 173 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_pidhit5

int SiPixelErrorEstimation::strip_pidhit5

private

Definition at line 176 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_prob

float SiPixelErrorEstimation::strip_prob

private

Definition at line 159 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_qbin

int SiPixelErrorEstimation::strip_qbin

private

Definition at line 160 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechiterrx

float SiPixelErrorEstimation::strip_rechiterrx

private

Definition at line 102 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechiterry

float SiPixelErrorEstimation::strip_rechiterry

private

Definition at line 103 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechitpullx

float SiPixelErrorEstimation::strip_rechitpullx

private

Definition at line 109 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechitpully

float SiPixelErrorEstimation::strip_rechitpully

private

Definition at line 110 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechitresx

float SiPixelErrorEstimation::strip_rechitresx

private

Definition at line 104 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechitresx2

float SiPixelErrorEstimation::strip_rechitresx2

private

Definition at line 106 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechitresy

float SiPixelErrorEstimation::strip_rechitresy

private

Definition at line 108 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechitx

float SiPixelErrorEstimation::strip_rechitx

private

Definition at line 99 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechity

float SiPixelErrorEstimation::strip_rechity

private

Definition at line 100 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_rechitz

float SiPixelErrorEstimation::strip_rechitz

private

Definition at line 101 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_simproc

int SiPixelErrorEstimation::strip_simproc

private

Definition at line 126 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_simproc1

int SiPixelErrorEstimation::strip_simproc1

private

Definition at line 165 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_simproc2

int SiPixelErrorEstimation::strip_simproc2

private

Definition at line 168 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_simproc3

int SiPixelErrorEstimation::strip_simproc3

private

Definition at line 171 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_simproc4

int SiPixelErrorEstimation::strip_simproc4

private

Definition at line 174 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_simproc5

int SiPixelErrorEstimation::strip_simproc5

private

Definition at line 177 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_size

int SiPixelErrorEstimation::strip_size

private

Definition at line 122 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_split

int SiPixelErrorEstimation::strip_split

private

Definition at line 179 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_subdet_id

int SiPixelErrorEstimation::strip_subdet_id

private

Definition at line 128 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_is_double_side

int SiPixelErrorEstimation::strip_tib_is_double_side

private

Definition at line 134 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_is_external_string

int SiPixelErrorEstimation::strip_tib_is_external_string

private

Definition at line 141 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_is_internal_string

int SiPixelErrorEstimation::strip_tib_is_internal_string

private

Definition at line 140 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_is_rphi

int SiPixelErrorEstimation::strip_tib_is_rphi

private

Definition at line 142 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_is_stereo

int SiPixelErrorEstimation::strip_tib_is_stereo

private

Definition at line 143 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_is_z_minus_side

int SiPixelErrorEstimation::strip_tib_is_z_minus_side

private

Definition at line 136 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_is_z_plus_side

int SiPixelErrorEstimation::strip_tib_is_z_plus_side

private

Definition at line 135 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_layer

int SiPixelErrorEstimation::strip_tib_layer

private

Definition at line 130 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_layer_number

int SiPixelErrorEstimation::strip_tib_layer_number

private

Definition at line 137 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_module

int SiPixelErrorEstimation::strip_tib_module

private

Definition at line 131 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_module_number

int SiPixelErrorEstimation::strip_tib_module_number

private

Definition at line 139 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_order

int SiPixelErrorEstimation::strip_tib_order

private

Definition at line 132 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_side

int SiPixelErrorEstimation::strip_tib_side

private

Definition at line 133 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tib_string_number

int SiPixelErrorEstimation::strip_tib_string_number

private

Definition at line 138 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_is_double_side

int SiPixelErrorEstimation::strip_tob_is_double_side

private

Definition at line 149 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_is_rphi

int SiPixelErrorEstimation::strip_tob_is_rphi

private

Definition at line 156 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ strip_tob_is_stereo

int SiPixelErrorEstimation::strip_tob_is_stereo

private

Definition at line 157 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ strip_tob_is_z_minus_side

int SiPixelErrorEstimation::strip_tob_is_z_minus_side

private

Definition at line 151 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_is_z_plus_side

int SiPixelErrorEstimation::strip_tob_is_z_plus_side

private

Definition at line 150 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_layer

int SiPixelErrorEstimation::strip_tob_layer

private

Definition at line 145 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_layer_number

int SiPixelErrorEstimation::strip_tob_layer_number

private

Definition at line 152 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_module

int SiPixelErrorEstimation::strip_tob_module

private

Definition at line 146 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_module_number

int SiPixelErrorEstimation::strip_tob_module_number

private

Definition at line 154 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_rod_number

int SiPixelErrorEstimation::strip_tob_rod_number

private

Definition at line 153 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_tob_side

int SiPixelErrorEstimation::strip_tob_side

private

Definition at line 148 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ strip_trk_pt

float SiPixelErrorEstimation::strip_trk_pt

private

Definition at line 115 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ subdetId

int SiPixelErrorEstimation::subdetId

private

Definition at line 203 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ tfile_

TFile* SiPixelErrorEstimation::tfile_

private

Definition at line 367 of file SiPixelErrorEstimation.h.

Referenced by beginJob(), and endJob().

◆ tPixRecHitCollection

edm::EDGetTokenT<SiPixelRecHitCollection> SiPixelErrorEstimation::tPixRecHitCollection

private

Definition at line 76 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ trackerGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiPixelErrorEstimation::trackerGeomToken_

private

Definition at line 80 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ trackerHitAssociatorConfig_

TrackerHitAssociator::Config SiPixelErrorEstimation::trackerHitAssociatorConfig_

private

Definition at line 373 of file SiPixelErrorEstimation.h.

Referenced by analyze().

◆ trackerTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiPixelErrorEstimation::trackerTopoToken_

private

Definition at line 79 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ trk_alpha

float SiPixelErrorEstimation::trk_alpha

private

Definition at line 197 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ trk_beta

float SiPixelErrorEstimation::trk_beta

private

Definition at line 198 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ tSimTrackContainer

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

private

Definition at line 77 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ tTrackCollection

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

private

Definition at line 78 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ tTrajectory

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

private

Definition at line 75 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and SiPixelErrorEstimation().

◆ ttree_all_hits_

TTree* SiPixelErrorEstimation::ttree_all_hits_

private

Definition at line 368 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ ttree_track_hits_

TTree* SiPixelErrorEstimation::ttree_track_hits_

private

Definition at line 369 of file SiPixelErrorEstimation.h.

Referenced by analyze(), and beginJob().

◆ ttree_track_hits_strip_

TTree* SiPixelErrorEstimation::ttree_track_hits_strip_

private

Definition at line 371 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

◆ SiPixelErrorEstimation()

◆ ~SiPixelErrorEstimation()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ computeAnglesFromDetPosition()

◆ endJob()

Member Data Documentation

◆ all_alpha

◆ all_beta

◆ all_bigx

◆ all_bigy

◆ all_blade

◆ all_clust_alpha

◆ all_clust_beta

◆ all_clust_col

◆ all_clust_geoid

◆ all_clust_maxpixcol

◆ all_clust_maxpixrow

◆ all_clust_minpixcol

◆ all_clust_minpixrow

◆ all_clust_q

◆ all_clust_row

◆ all_clust_x

◆ all_clust_y

◆ all_col1

◆ all_col2

◆ all_cols

◆ all_disk

◆ all_edgex

◆ all_edgey

◆ all_eloss

◆ all_flipped

◆ all_gx1

◆ all_gx2

◆ all_gy1

◆ all_gy2

◆ all_gz1

◆ all_gz2

◆ all_half

◆ all_hit_cprob0

◆ all_hit_cprob1

◆ all_hit_cprob2

◆ all_hit_probx

◆ all_hit_proby

◆ all_ladder

◆ all_layer

◆ all_mod

◆ all_npix

◆ all_nsimhit

◆ all_nxpix

◆ all_nypix

◆ all_panel

◆ all_pidhit

◆ all_pixadc

◆ all_pixcol

◆ all_pixel_clst_err_x

◆ all_pixel_clst_err_y

◆ all_pixel_split

◆ all_pixgx

◆ all_pixgy

◆ all_pixgz

◆ all_pixrow

◆ all_pixx

◆ all_pixy

◆ all_plaq

◆ all_rechiterrx

◆ all_rechiterry

◆ all_rechitpullx

◆ all_rechitpully

◆ all_rechitresx

◆ all_rechitresy

◆ all_rechitx

◆ all_rechity