SiPixelTrackResidualSource Class Reference

#include <SiPixelTrackResidualSource.h>

Inheritance diagram for SiPixelTrackResidualSource:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
virtual void	dqmBeginRun (const edm::Run &r, edm::EventSetup const &iSetup)
	SiPixelTrackResidualSource (const edm::ParameterSet &)
void	triplets (double x1, double y1, double z1, double x2, double y2, double z2, double x3, double y3, double z3, double ptsig, double &dc, double &dz, double kap)
	~SiPixelTrackResidualSource ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Public Attributes
std::string	topFolderName_

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
bool	bladeOn
edm::InputTag	clustersrc_
edm::EDGetTokenT < edmNew::DetSetVector < SiPixelCluster > >	clustersrcToken_
bool	debug_
bool	diskOn
bool	firstRun
edm::EDGetTokenT < reco::TrackCollection >	generalTracksToken_
bool	isUpgrade
bool	ladOn
bool	layOn
MonitorElement *	meClChargeNotOnTrack_all
MonitorElement *	meClChargeNotOnTrack_bpix
std::vector< MonitorElement * >	meClChargeNotOnTrack_diskms
std::vector< MonitorElement * >	meClChargeNotOnTrack_diskps
MonitorElement *	meClChargeNotOnTrack_fpix
std::vector< MonitorElement * >	meClChargeNotOnTrack_layers
MonitorElement *	meClChargeOnTrack_all
MonitorElement *	meClChargeOnTrack_bpix
std::vector< MonitorElement * >	meClChargeOnTrack_diskms
std::vector< MonitorElement * >	meClChargeOnTrack_diskps
MonitorElement *	meClChargeOnTrack_fpix
std::vector< MonitorElement * >	meClChargeOnTrack_layers
std::vector< MonitorElement * >	meClPosDisksmzNotOnTrack
std::vector< MonitorElement * >	meClPosDisksmzOnTrack
std::vector< MonitorElement * >	meClPosDiskspzNotOnTrack
std::vector< MonitorElement * >	meClPosDiskspzOnTrack
std::vector< MonitorElement * >	meClPosLayersNotOnTrack
std::vector< MonitorElement * >	meClPosLayersOnTrack
MonitorElement *	meClSizeNotOnTrack_all
MonitorElement *	meClSizeNotOnTrack_bpix
std::vector< MonitorElement * >	meClSizeNotOnTrack_diskms
std::vector< MonitorElement * >	meClSizeNotOnTrack_diskps
MonitorElement *	meClSizeNotOnTrack_fpix
std::vector< MonitorElement * >	meClSizeNotOnTrack_layers
MonitorElement *	meClSizeOnTrack_all
MonitorElement *	meClSizeOnTrack_bpix
std::vector< MonitorElement * >	meClSizeOnTrack_diskms
std::vector< MonitorElement * >	meClSizeOnTrack_diskps
MonitorElement *	meClSizeOnTrack_fpix
std::vector< MonitorElement * >	meClSizeOnTrack_layers
MonitorElement *	meClSizeXNotOnTrack_all
MonitorElement *	meClSizeXNotOnTrack_bpix
std::vector< MonitorElement * >	meClSizeXNotOnTrack_diskms
std::vector< MonitorElement * >	meClSizeXNotOnTrack_diskps
MonitorElement *	meClSizeXNotOnTrack_fpix
std::vector< MonitorElement * >	meClSizeXNotOnTrack_layers
MonitorElement *	meClSizeXOnTrack_all
MonitorElement *	meClSizeXOnTrack_bpix
std::vector< MonitorElement * >	meClSizeXOnTrack_diskms
std::vector< MonitorElement * >	meClSizeXOnTrack_diskps
MonitorElement *	meClSizeXOnTrack_fpix
std::vector< MonitorElement * >	meClSizeXOnTrack_layers
MonitorElement *	meClSizeYNotOnTrack_all
MonitorElement *	meClSizeYNotOnTrack_bpix
std::vector< MonitorElement * >	meClSizeYNotOnTrack_diskms
std::vector< MonitorElement * >	meClSizeYNotOnTrack_diskps
MonitorElement *	meClSizeYNotOnTrack_fpix
std::vector< MonitorElement * >	meClSizeYNotOnTrack_layers
MonitorElement *	meClSizeYOnTrack_all
MonitorElement *	meClSizeYOnTrack_bpix
std::vector< MonitorElement * >	meClSizeYOnTrack_diskms
std::vector< MonitorElement * >	meClSizeYOnTrack_diskps
MonitorElement *	meClSizeYOnTrack_fpix
std::vector< MonitorElement * >	meClSizeYOnTrack_layers
MonitorElement *	meHitProbability
MonitorElement *	meNClustersNotOnTrack_all
MonitorElement *	meNClustersNotOnTrack_bpix
std::vector< MonitorElement * >	meNClustersNotOnTrack_diskms
std::vector< MonitorElement * >	meNClustersNotOnTrack_diskps
MonitorElement *	meNClustersNotOnTrack_fpix
std::vector< MonitorElement * >	meNClustersNotOnTrack_layers
MonitorElement *	meNClustersOnTrack_all
MonitorElement *	meNClustersOnTrack_bpix
std::vector< MonitorElement * >	meNClustersOnTrack_diskms
std::vector< MonitorElement * >	meNClustersOnTrack_diskps
MonitorElement *	meNClustersOnTrack_fpix
std::vector< MonitorElement * >	meNClustersOnTrack_layers
MonitorElement *	meNofClustersNotOnTrack_
MonitorElement *	meNofClustersOnTrack_
MonitorElement *	meNofTracks_
MonitorElement *	meNofTracksInPixVol_
std::vector< MonitorElement * >	meResidualXSummedLay
std::vector< MonitorElement * >	meResidualYSummedLay
MonitorElement *	meSubdetResidualX [3]
MonitorElement *	meSubdetResidualY [3]
bool	modOn
int	NLowProb
int	noOfDisks
int	noOfLayers
int	NTotal
edm::EDGetTokenT < reco::VertexCollection >	offlinePrimaryVerticesToken_
bool	phiOn
edm::ParameterSet	pSet_
double	ptminres_
bool	reducedSet
bool	ringOn
edm::InputTag	src_
std::map< uint32_t, SiPixelTrackResidualModule * >	theSiPixelStructure
edm::EDGetTokenT < TrajTrackAssociationCollection >	trackAssociationToken_
edm::InputTag	tracksrc_
edm::EDGetTokenT< std::vector < Trajectory > >	tracksrcToken_
edm::EDGetTokenT< std::vector < reco::Track > >	trackToken_
std::string	ttrhbuilder_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 48 of file SiPixelTrackResidualSource.h.

Constructor & Destructor Documentation

SiPixelTrackResidualSource::SiPixelTrackResidualSource ( const edm::ParameterSet &  pSet )

explicit

Definition at line 66 of file SiPixelTrackResidualSource.cc.

References beamSpotToken_, bladeOn, clustersrc_, clustersrcToken_, debug_, diskOn, firstRun, generalTracksToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ladOn, layOn, modOn, NLowProb, NTotal, offlinePrimaryVerticesToken_, phiOn, pSet_, ptminres_, ringOn, src_, AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, trackAssociationToken_, tracksrc_, tracksrcToken_, trackToken_, and ttrhbuilder_.

                                                                                  :
  pSet_(pSet),
  modOn( pSet.getUntrackedParameter<bool>("modOn",true) ),
  reducedSet( pSet.getUntrackedParameter<bool>("reducedSet",true) ),
  ladOn( pSet.getUntrackedParameter<bool>("ladOn",false) ), 
  layOn( pSet.getUntrackedParameter<bool>("layOn",false) ), 
  phiOn( pSet.getUntrackedParameter<bool>("phiOn",false) ), 
  ringOn( pSet.getUntrackedParameter<bool>("ringOn",false) ), 
  bladeOn( pSet.getUntrackedParameter<bool>("bladeOn",false) ), 
  diskOn( pSet.getUntrackedParameter<bool>("diskOn",false) ),
  isUpgrade( pSet.getUntrackedParameter<bool>("isUpgrade",false) ),
  noOfLayers(0),
  noOfDisks(0)
{ 
   pSet_ = pSet; 
   debug_ = pSet_.getUntrackedParameter<bool>("debug", false); 
   src_ = pSet_.getParameter<edm::InputTag>("src"); 
   clustersrc_ = pSet_.getParameter<edm::InputTag>("clustersrc");
   tracksrc_ = pSet_.getParameter<edm::InputTag>("trajectoryInput");
   ttrhbuilder_ = pSet_.getParameter<std::string>("TTRHBuilder");
   ptminres_= pSet.getUntrackedParameter<double>("PtMinRes",4.0) ;
   beamSpotToken_ = consumes<reco::BeamSpot>(std::string("offlineBeamSpot"));
   offlinePrimaryVerticesToken_ = consumes<reco::VertexCollection>(std::string("offlinePrimaryVertices"));
   generalTracksToken_ = consumes<reco::TrackCollection>(pSet_.getParameter<edm::InputTag>("tracksrc"));
   tracksrcToken_ = consumes<std::vector<Trajectory> >(pSet_.getParameter<edm::InputTag>("trajectoryInput"));
   trackToken_ = consumes<std::vector<reco::Track> >(pSet_.getParameter<edm::InputTag>("trajectoryInput"));
   trackAssociationToken_ = consumes<TrajTrackAssociationCollection>(pSet_.getParameter<edm::InputTag>("trajectoryInput"));
   clustersrcToken_ = consumes<edmNew::DetSetVector<SiPixelCluster> >(pSet_.getParameter<edm::InputTag>("clustersrc"));

  LogInfo("PixelDQM") << "SiPixelTrackResidualSource constructor" << endl;
  LogInfo ("PixelDQM") << "Mod/Lad/Lay/Phi " << modOn << "/" << ladOn << "/" 
            << layOn << "/" << phiOn << std::endl;
  LogInfo ("PixelDQM") << "Blade/Disk/Ring" << bladeOn << "/" << diskOn << "/" 
            << ringOn << std::endl;

  topFolderName_ = pSet_.getParameter<std::string>("TopFolderName");

  firstRun = true;
  NTotal=0;
  NLowProb=0;
}

SiPixelTrackResidualSource::~SiPixelTrackResidualSource

(

)

Definition at line 109 of file SiPixelTrackResidualSource.cc.

References theSiPixelStructure.

                                                        {
  LogInfo("PixelDQM") << "SiPixelTrackResidualSource destructor" << endl;

  std::map<uint32_t,SiPixelTrackResidualModule*>::iterator struct_iter;
  for (struct_iter = theSiPixelStructure.begin() ; struct_iter != theSiPixelStructure.end() ; struct_iter++){
    delete struct_iter->second;
    struct_iter->second = 0;
  }
}

Member Function Documentation

void SiPixelTrackResidualSource::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 608 of file SiPixelTrackResidualSource.cc.

References funct::abs(), Reference_intrackfit_cff::barrel, beamSpotToken_, edm::AssociationMap< Tag >::begin(), edmNew::DetSetVector< T >::begin(), bladeOn, RecoTauCleanerPlugins::charge, SiPixelRecHit::clusterProbability(), clustersrcToken_, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, gather_cfg::cout, debug_, DetId::det(), PixelEndcapName::diskName(), diskOn, edm::AssociationMap< Tag >::end(), edmNew::DetSetVector< T >::end(), Reference_intrackfit_cff::endcap, edm::HandleBase::failedToGet(), MonitorElement::Fill(), HcalObjRepresent::Fill(), edmNew::DetSetVector< T >::find(), generalTracksToken_, edm::EventSetup::get(), edm::Event::getByToken(), corrVsCorr::gX, corrVsCorr::gY, SiPixelRecHit::hasFilledProb(), TrackingRecHit::hit(), i, TrackerGeometry::idToDet(), reco::TransientTrack::initialFreeState(), reco::TransientTrack::innermostMeasurementState(), edm::Ref< C, T, F >::isNonnull(), isUpgrade, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), edm::Ref< C, T, F >::key(), GetRecoTauVFromDQM_MC_cff::kk, ladOn, PixelBarrelNameUpgrade::layerName(), PixelBarrelName::layerName(), layOn, TrajectoryStateOnSurface::localParameters(), Topology::localPosition(), PV3DBase< T, PVType, FrameType >::mag(), match(), meClChargeNotOnTrack_all, meClChargeNotOnTrack_bpix, meClChargeNotOnTrack_diskms, meClChargeNotOnTrack_diskps, meClChargeNotOnTrack_fpix, meClChargeNotOnTrack_layers, meClChargeOnTrack_all, meClChargeOnTrack_bpix, meClChargeOnTrack_diskms, meClChargeOnTrack_diskps, meClChargeOnTrack_fpix, meClChargeOnTrack_layers, meClPosDisksmzNotOnTrack, meClPosDisksmzOnTrack, meClPosDiskspzNotOnTrack, meClPosDiskspzOnTrack, meClPosLayersNotOnTrack, meClPosLayersOnTrack, meClSizeNotOnTrack_all, meClSizeNotOnTrack_bpix, meClSizeNotOnTrack_diskms, meClSizeNotOnTrack_diskps, meClSizeNotOnTrack_fpix, meClSizeNotOnTrack_layers, meClSizeOnTrack_all, meClSizeOnTrack_bpix, meClSizeOnTrack_diskms, meClSizeOnTrack_diskps, meClSizeOnTrack_fpix, meClSizeOnTrack_layers, meClSizeXNotOnTrack_all, meClSizeXNotOnTrack_bpix, meClSizeXNotOnTrack_diskms, meClSizeXNotOnTrack_diskps, meClSizeXNotOnTrack_fpix, meClSizeXNotOnTrack_layers, meClSizeXOnTrack_all, meClSizeXOnTrack_bpix, meClSizeXOnTrack_diskms, meClSizeXOnTrack_diskps, meClSizeXOnTrack_fpix, meClSizeXOnTrack_layers, meClSizeYNotOnTrack_all, meClSizeYNotOnTrack_bpix, meClSizeYNotOnTrack_diskms, meClSizeYNotOnTrack_diskps, meClSizeYNotOnTrack_fpix, meClSizeYNotOnTrack_layers, meClSizeYOnTrack_all, meClSizeYOnTrack_bpix, meClSizeYOnTrack_diskms, meClSizeYOnTrack_diskps, meClSizeYOnTrack_fpix, meClSizeYOnTrack_layers, meHitProbability, meNClustersNotOnTrack_all, meNClustersNotOnTrack_bpix, meNClustersNotOnTrack_diskms, meNClustersNotOnTrack_diskps, meNClustersNotOnTrack_fpix, meNClustersNotOnTrack_layers, meNClustersOnTrack_all, meNClustersOnTrack_bpix, meNClustersOnTrack_diskms, meNClustersOnTrack_diskps, meNClustersOnTrack_fpix, meNClustersOnTrack_layers, meNofClustersNotOnTrack_, meNofClustersOnTrack_, meNofTracks_, meNofTracksInPixVol_, meResidualXSummedLay, meResidualYSummedLay, modOn, LocalTrajectoryParameters::momentum(), reco::Vertex::ndof(), NLowProb, noOfDisks, noOfLayers, NTotal, offlinePrimaryVerticesToken_, phi, PV3DBase< T, PVType, FrameType >::phi(), phiOn, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, funct::pow(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), EnergyCorrector::pt, ptminres_, TrackerTopology::pxbLayer(), DetId::rawId(), reducedSet, ringOn, edm::AssociationMap< Tag >::size(), edmNew::DetSetVector< T >::size(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), DetId::subdetId(), GeomDet::surface(), funct::tan(), groupFilesInBlocks::temp, theSiPixelStructure, Surface::toGlobal(), trackAssociationToken_, DetId::Tracker, testEve_cfg::tracks, tracksrcToken_, trackToken_, FreeTrajectoryState::transverseCurvature(), triplets(), ttrhbuilder_, HLT_25ns14e33_v1_cff::vertices, x, PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), detailsBasic3DVector::y, PV2DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

                                                                                            {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();
  
  // retrieve TrackerGeometry again and MagneticField for use in transforming 
  // a TrackCandidate's P(ersistent)TrajectoryStateoOnDet (PTSoD) to a TrajectoryStateOnSurface (TSoS)
  ESHandle<TrackerGeometry> TG;
  iSetup.get<TrackerDigiGeometryRecord>().get(TG);
  const TrackerGeometry* theTrackerGeometry = TG.product();
  
  //analytic triplet method to calculate the track residuals in the pixe barrel detector    
  
  //--------------------------------------------------------------------
  // beam spot:
  //
  edm::Handle<reco::BeamSpot> rbs;
  //iEvent.getByLabel( "offlineBeamSpot", rbs );
  iEvent.getByToken( beamSpotToken_, rbs );
  math::XYZPoint bsP = math::XYZPoint(0,0,0);
  if( !rbs.failedToGet() && rbs.isValid() )
    {
      bsP = math::XYZPoint( rbs->x0(), rbs->y0(), rbs->z0() );
    }
  
  //--------------------------------------------------------------------
  // primary vertices:
  //
  edm::Handle<reco::VertexCollection> vertices;
  //iEvent.getByLabel("offlinePrimaryVertices", vertices );
  iEvent.getByToken(  offlinePrimaryVerticesToken_, vertices );
  
  if( vertices.failedToGet() ) return;
  if( !vertices.isValid() ) return;
  
  math::XYZPoint vtxN = math::XYZPoint(0,0,0);
  math::XYZPoint vtxP = math::XYZPoint(0,0,0);
  
  double bestNdof = 0.0;
  double maxSumPt = 0.0;
  reco::Vertex bestPvx;
  for(reco::VertexCollection::const_iterator iVertex = vertices->begin();
      iVertex != vertices->end(); ++iVertex ) {
    if( iVertex->ndof() > bestNdof ) {
      bestNdof = iVertex->ndof();
      vtxN = math::XYZPoint( iVertex->x(), iVertex->y(), iVertex->z() );
    }//ndof
    if( iVertex->p4().pt() > maxSumPt ) {
      maxSumPt = iVertex->p4().pt();
      vtxP = math::XYZPoint( iVertex->x(), iVertex->y(), iVertex->z() );
      bestPvx = *iVertex;
    }//sumpt
    
  }//vertex
  
  if( maxSumPt < 1.0 ) vtxP = vtxN;
  
  //---------------------------------------------
  //get Tracks
  //
  edm:: Handle<reco::TrackCollection> TracksForRes;
  //iEvent.getByLabel( "generalTracks", TracksForRes );
  iEvent.getByToken(  generalTracksToken_, TracksForRes );

  //
  // transient track builder, needs B-field from data base (global tag in .py)
  //
  edm::ESHandle<TransientTrackBuilder> theB;
  iSetup.get<TransientTrackRecord>().get( "TransientTrackBuilder", theB );

  //get the TransienTrackingRecHitBuilder needed for extracting the global position of the hits in the pixel
  edm::ESHandle<TransientTrackingRecHitBuilder> theTrackerRecHitBuilder;
  iSetup.get<TransientRecHitRecord>().get(ttrhbuilder_,theTrackerRecHitBuilder);

  //check that tracks are valid
  if( TracksForRes.failedToGet() ) return;
  if( !TracksForRes.isValid() ) return;

  //get tracker geometry
  edm::ESHandle<TrackerGeometry> pDD;
  iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
  
  if( !pDD.isValid() ) {
    cout << "Unable to find TrackerDigiGeometry. Return\n";
    return;
  }
 
  int kk = -1;
  //----------------------------------------------------------------------------
  // Residuals:
  //
  for( reco::TrackCollection::const_iterator iTrack = TracksForRes->begin();
       iTrack != TracksForRes->end(); ++iTrack ) {
    //count
    kk++;
    //Calculate minimal track pt before curling
    // cpt = cqRB = 0.3*R[m]*B[T] = 1.14*R[m] for B=3.8T
    // D = 2R = 2*pt/1.14
    // calo: D = 1.3 m => pt = 0.74 GeV/c
    double pt = iTrack->pt();
    if( pt < 0.75 ) continue;// curls up
    if( abs( iTrack->dxy(vtxP) ) > 5*iTrack->dxyError() ) continue; // not prompt
    
    double charge = iTrack->charge();
       
    reco::TransientTrack tTrack = theB->build(*iTrack);
    //get curvature of the track, needed for the residuals
    double kap = tTrack.initialFreeState().transverseCurvature();
    //needed for the TransienTrackingRecHitBuilder
    TrajectoryStateOnSurface initialTSOS = tTrack.innermostMeasurementState();
    if( iTrack->extra().isNonnull() &&iTrack->extra().isAvailable() ){
      
      double x1 = 0;
      double y1 = 0;
      double z1 = 0;
      double x2 = 0;
      double y2 = 0;
      double z2 = 0;
      double x3 = 0;
      double y3 = 0;
      double z3 = 0;
      int n1 = 0;
      int n2 = 0;
      int n3 = 0;
      
      //for saving the pixel barrel hits
      vector<TransientTrackingRecHit::RecHitPointer> GoodPixBarrelHits;
      //looping through the RecHits of the track
      for( trackingRecHit_iterator irecHit = iTrack->recHitsBegin();
       irecHit != iTrack->recHitsEnd(); ++irecHit){
    
    if( (*irecHit)->isValid() ){
      DetId detId = (*irecHit)->geographicalId();
      // enum Detector { Tracker=1, Muon=2, Ecal=3, Hcal=4, Calo=5 };
      if( detId.det() != 1 ){
        if(debug_){
          cout << "rec hit ID = " << detId.det() << " not in tracker!?!?\n";
        }
        continue;
      }
      uint32_t subDet = detId.subdetId();
      
      // enum SubDetector{ PixelBarrel=1, PixelEndcap=2 };
      // enum SubDetector{ TIB=3, TID=4, TOB=5, TEC=6 };
      
      TransientTrackingRecHit::RecHitPointer trecHit = theTrackerRecHitBuilder->build(  &*(*irecHit), initialTSOS);
      
      
      double gX = trecHit->globalPosition().x();
      double gY = trecHit->globalPosition().y();
      double gZ = trecHit->globalPosition().z();
          
          
          if( subDet == PixelSubdetector::PixelBarrel ) {

        int ilay = tTopo->pxbLayer(detId);
        
        if( ilay == 1 ){
          n1++;
          x1 = gX;
          y1 = gY;
          z1 = gZ;
          
          GoodPixBarrelHits.push_back((trecHit));
        }//PXB1
        if( ilay == 2 ){
          
          n2++;
          x2 = gX;
          y2 = gY;
          z2 = gZ;
          
          GoodPixBarrelHits.push_back((trecHit));
          
        }//PXB2
        if( ilay == 3 ){
          
          n3++;
          x3 = gX;
          y3 = gY;
          z3 = gZ;
          GoodPixBarrelHits.push_back((trecHit));
        }
          }//PXB
          
          
        }//valid
      }//loop rechits
      
      //CS extra plots    
      
          
      if( n1+n2+n3 == 3 && n1*n2*n3 > 0) {      
        for( unsigned int i = 0; i < GoodPixBarrelHits.size(); i++){
          
          if( GoodPixBarrelHits[i]->isValid() ){
        DetId detId = GoodPixBarrelHits[i]->geographicalId().rawId();
        int ilay = tTopo->pxbLayer(detId);
        if(pt > ptminres_){   
          
          double dca2 = 0.0, dz2=0.0;
          double ptsig = pt;
          if(charge<0.) ptsig = -pt;
          //Filling the histograms in modules
          
          MeasurementPoint Test;
          MeasurementPoint Test2;
          Test=MeasurementPoint(0,0);
          Test2=MeasurementPoint(0,0);
          Measurement2DVector residual;
         
          if( ilay == 1 ){
            
            triplets(x2,y2,z2,x1,y1,z1,x3,y3,z3,ptsig,dca2,dz2, kap);   
                    
            Test=MeasurementPoint(dca2*1E4,dz2*1E4);
            residual=Test-Test2;
          }
          
          if( ilay == 2 ){
            
            triplets(x1,y1,z1,x2,y2,z2,x3,y3,z3,ptsig,dca2,dz2, kap);
            
            Test=MeasurementPoint(dca2*1E4,dz2*1E4);
            residual=Test-Test2;
          
          }
          
          if( ilay == 3 ){
            
            triplets(x1,y1,z1,x3,y3,z3,x2,y2,z2,ptsig,dca2,dz2, kap);
            
            Test=MeasurementPoint(dca2*1E4,dz2*1E4);
            residual=Test-Test2;
        }
          // fill the residual histograms 

          std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find(detId);
          if (pxd!=theSiPixelStructure.end()) (*pxd).second->fill(residual, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);   

          if(ladOn&&pxd!=theSiPixelStructure.end()){  
            meResidualXSummedLay.at(PixelBarrelNameUpgrade((*pxd).first).layerName()-1)->Fill(residual.x());    
            meResidualYSummedLay.at(PixelBarrelNameUpgrade((*pxd).first).layerName()-1)->Fill(residual.y());    
          }
    
        }//three hits
          }//is valid
        }//rechits loop
      }//pt 4
    }
    
    
 }//-----Tracks
  //get trajectories
  edm::Handle<std::vector<Trajectory> > trajCollectionHandle;
  //iEvent.getByLabel(tracksrc_,trajCollectionHandle);
  iEvent.getByToken ( tracksrcToken_, trajCollectionHandle );
  const std::vector<Trajectory> trajColl = *(trajCollectionHandle.product());
   
  //get tracks
  edm::Handle<std::vector<reco::Track> > trackCollectionHandle;
  //iEvent.getByLabel(tracksrc_,trackCollectionHandle);
  iEvent.getByToken( trackToken_, trackCollectionHandle );

  const std::vector<reco::Track> trackColl = *(trackCollectionHandle.product());

  //get the map
  edm::Handle<TrajTrackAssociationCollection> match;
  //iEvent.getByLabel(tracksrc_,match);
  iEvent.getByToken( trackAssociationToken_, match);
  const TrajTrackAssociationCollection ttac = *(match.product());

  // get clusters
  edm::Handle< edmNew::DetSetVector<SiPixelCluster> >  clusterColl;
  //iEvent.getByLabel( clustersrc_, clusterColl );
  iEvent.getByToken( clustersrcToken_, clusterColl );
  const edmNew::DetSetVector<SiPixelCluster> clustColl = *(clusterColl.product());

  if(debug_){
    std::cout << "Trajectories\t : " << trajColl.size() << std::endl;
    std::cout << "recoTracks  \t : " << trackColl.size() << std::endl;
    std::cout << "Map entries \t : " << ttac.size() << std::endl;
  }

  std::set<SiPixelCluster> clusterSet;
  TrajectoryStateCombiner tsoscomb;
  int tracks=0, pixeltracks=0, bpixtracks=0, fpixtracks=0; 
  int trackclusters=0, barreltrackclusters=0, endcaptrackclusters=0;
  int otherclusters=0, barrelotherclusters=0, endcapotherclusters=0;

  //Loop over map entries
  for(TrajTrackAssociationCollection::const_iterator it =  ttac.begin();it !=  ttac.end(); ++it){
    const edm::Ref<std::vector<Trajectory> > traj_iterator = it->key;  
    // Trajectory Map, extract Trajectory for this track
    reco::TrackRef trackref = it->val;
    tracks++;

    bool isBpixtrack = false, isFpixtrack = false, crossesPixVol=false;
    
    //find out whether track crosses pixel fiducial volume (for cosmic tracks)
    
    double d0 = (*trackref).d0(), dz = (*trackref).dz(); 
    
    if(abs(d0)<15 && abs(dz)<50) crossesPixVol = true;

    std::vector<TrajectoryMeasurement> tmeasColl =traj_iterator->measurements();
    std::vector<TrajectoryMeasurement>::const_iterator tmeasIt;
    //loop on measurements to find out whether there are bpix and/or fpix hits
    for(tmeasIt = tmeasColl.begin();tmeasIt!=tmeasColl.end();tmeasIt++){
      if(! tmeasIt->updatedState().isValid()) continue; 
      TransientTrackingRecHit::ConstRecHitPointer testhit = tmeasIt->recHit();
      if(! testhit->isValid() || testhit->geographicalId().det() != DetId::Tracker) continue; 
      uint testSubDetID = (testhit->geographicalId().subdetId()); 
      if(testSubDetID==PixelSubdetector::PixelBarrel) isBpixtrack = true;
      if(testSubDetID==PixelSubdetector::PixelEndcap) isFpixtrack = true;
    }//end loop on measurements
    if(isBpixtrack) {
      bpixtracks++; 
      if(debug_) std::cout << "bpixtrack\n";
    }
    if(isFpixtrack) {
      fpixtracks++; 
      if(debug_) std::cout << "fpixtrack\n";
    }
    if(isBpixtrack || isFpixtrack){
      pixeltracks++;
      
      if(crossesPixVol) meNofTracksInPixVol_->Fill(0,1);

      std::vector<TrajectoryMeasurement> tmeasColl = traj_iterator->measurements();
      for(std::vector<TrajectoryMeasurement>::const_iterator tmeasIt = tmeasColl.begin(); tmeasIt!=tmeasColl.end(); tmeasIt++){   
    if(! tmeasIt->updatedState().isValid()) continue; 
    
    TrajectoryStateOnSurface tsos = tsoscomb( tmeasIt->forwardPredictedState(), tmeasIt->backwardPredictedState() );
    TransientTrackingRecHit::ConstRecHitPointer hit = tmeasIt->recHit();
    if(! hit->isValid() || hit->geographicalId().det() != DetId::Tracker ) {
      continue; 
    } else {
      
//    //residual
          const DetId & hit_detId = hit->geographicalId();
      //uint IntRawDetID = (hit_detId.rawId());
      uint IntSubDetID = (hit_detId.subdetId());
      
      if(IntSubDetID == 0 ) continue; // don't look at SiStrip hits!        
      
      // get the enclosed persistent hit
      const TrackingRecHit *persistentHit = hit->hit();
      // check if it's not null, and if it's a valid pixel hit
      if ((persistentHit != 0) && (typeid(*persistentHit) == typeid(SiPixelRecHit))) {
        // tell the C++ compiler that the hit is a pixel hit
        const SiPixelRecHit* pixhit = static_cast<const SiPixelRecHit*>( hit->hit() );
        //Hit probability:
        float hit_prob = -1.;
        if(pixhit->hasFilledProb()){
          hit_prob = pixhit->clusterProbability(0);
          //std::cout<<"HITPROB= "<<hit_prob<<std::endl;
          if(hit_prob<pow(10.,-15.)) NLowProb++;
          NTotal++;
          if(NTotal>0) meHitProbability->Fill(float(NLowProb/NTotal));
        }
        
        // get the edm::Ref to the cluster
        edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> const& clust = (*pixhit).cluster();
        
        //  check if the ref is not null
        if (clust.isNonnull()) {
          
          //define tracker and pixel geometry and topology
          const TrackerGeometry& theTracker(*theTrackerGeometry);
          const PixelGeomDetUnit* theGeomDet = static_cast<const PixelGeomDetUnit*> (theTracker.idToDet(hit_detId) );
          
          //test if PixelGeomDetUnit exists
          if(theGeomDet == 0) {
        if(debug_) std::cout << "NO THEGEOMDET\n";
        continue;
          }
         
          const PixelTopology * topol = &(theGeomDet->specificTopology());
          //fill histograms for clusters on tracks
          //correct SiPixelTrackResidualModule
          std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find((*hit).geographicalId().rawId());

          //CHARGE CORRECTION (for track impact angle)
          // calculate alpha and beta from cluster position
          LocalTrajectoryParameters ltp = tsos.localParameters();
          LocalVector localDir = ltp.momentum()/ltp.momentum().mag();
         
          float clust_alpha = atan2(localDir.z(), localDir.x());
          float clust_beta = atan2(localDir.z(), localDir.y());
          double corrCharge = clust->charge() * sqrt( 1.0 / ( 1.0/pow( tan(clust_alpha), 2 ) + 
                                  1.0/pow( tan(clust_beta ), 2 ) + 
                                  1.0 )
                              )/1000.;
          
          if (pxd!=theSiPixelStructure.end()) (*pxd).second->fill((*clust), true, corrCharge, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);     


          trackclusters++;
          //CORR CHARGE
          meClChargeOnTrack_all->Fill(corrCharge);
          meClSizeOnTrack_all->Fill((*clust).size());
          meClSizeXOnTrack_all->Fill((*clust).sizeX());
          meClSizeYOnTrack_all->Fill((*clust).sizeY());
          clusterSet.insert(*clust);
          
          //find cluster global position (rphi, z)
          // get cluster center of gravity (of charge)
          float xcenter = clust->x();
          float ycenter = clust->y();
          // get the cluster position in local coordinates (cm) 
          LocalPoint clustlp = topol->localPosition( MeasurementPoint(xcenter, ycenter) );
          // get the cluster position in global coordinates (cm)
          GlobalPoint clustgp = theGeomDet->surface().toGlobal( clustlp );
          
          //find location of hit (barrel or endcap, same for cluster)
          bool barrel = DetId((*hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
          bool endcap = DetId((*hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
          if(barrel) {
        barreltrackclusters++;
        //CORR CHARGE
        meClChargeOnTrack_bpix->Fill(corrCharge);
        meClSizeOnTrack_bpix->Fill((*clust).size());
        meClSizeXOnTrack_bpix->Fill((*clust).sizeX());
        meClSizeYOnTrack_bpix->Fill((*clust).sizeY());
        int DBlayer;
        DBlayer = PixelBarrelName(DetId((*hit).geographicalId()), tTopo, isUpgrade).layerName();
        float phi = clustgp.phi(); 
        float z = clustgp.z();
        for (int i = 0; i < noOfLayers; i++)
          {
          if (DBlayer == i + 1) {
             meClPosLayersOnTrack.at(i)->Fill(z,phi);
             meClChargeOnTrack_layers.at(i)->Fill(corrCharge);
             meClSizeOnTrack_layers.at(i)->Fill((*clust).size());
             meClSizeXOnTrack_layers.at(i)->Fill((*clust).sizeX());
             meClSizeYOnTrack_layers.at(i)->Fill((*clust).sizeY());
          }
          }
          }
          if(endcap) {
        endcaptrackclusters++;
        //CORR CHARGE
        meClChargeOnTrack_fpix->Fill(corrCharge);
        meClSizeOnTrack_fpix->Fill((*clust).size());
        meClSizeXOnTrack_fpix->Fill((*clust).sizeX());
        meClSizeYOnTrack_fpix->Fill((*clust).sizeY());
        int DBdisk = 0;
        DBdisk = PixelEndcapName(DetId((*hit).geographicalId()), tTopo, isUpgrade).diskName();
        float x = clustgp.x(); 
        float y = clustgp.y(); 
        float z = clustgp.z();
        if(z>0){
        for (int i = 0; i < noOfDisks; i++)
        {
            if (DBdisk == i + 1) {
               meClPosDiskspzOnTrack.at(i)->Fill(x,y);
               meClChargeOnTrack_diskps.at(i)->Fill(corrCharge);
               meClSizeOnTrack_diskps.at(i)->Fill((*clust).size());
               meClSizeXOnTrack_diskps.at(i)->Fill((*clust).sizeX());
               meClSizeYOnTrack_diskps.at(i)->Fill((*clust).sizeY());
            }
        }
        }
        else{
        for (int i = 0; i < noOfDisks; i++)
        {
          if (DBdisk == i + 1) {
             meClPosDisksmzOnTrack.at(i)->Fill(x,y);
             meClChargeOnTrack_diskms.at(i)->Fill(corrCharge);
             meClSizeOnTrack_diskms.at(i)->Fill((*clust).size());
             meClSizeXOnTrack_diskms.at(i)->Fill((*clust).sizeX());
             meClSizeYOnTrack_diskms.at(i)->Fill((*clust).sizeY());
          }
        }
        } 
          }
          
        }//end if (cluster exists)

      }//end if (persistent hit exists and is pixel hit)
      
    }//end of else 
    
    
      }//end for (all traj measurements of pixeltrack)
    }//end if (is pixeltrack)
    else {
      if(debug_) std::cout << "no pixeltrack:\n";
      if(crossesPixVol) meNofTracksInPixVol_->Fill(1,1);
    }
    
  }//end loop on map entries

  //find clusters that are NOT on track
  //edmNew::DetSet<SiPixelCluster>::const_iterator  di;
  if(debug_) std::cout << "clusters not on track: (size " << clustColl.size() << ") ";

  for(TrackerGeometry::DetContainer::const_iterator it = TG->dets().begin(); it != TG->dets().end(); it++){
    //if(dynamic_cast<PixelGeomDetUnit const *>((*it))!=0){
    DetId detId = (*it)->geographicalId();
    if(detId>=302055684 && detId<=352477708){ // make sure it's a Pixel module WITHOUT using dynamic_cast!  
      int nofclOnTrack = 0, nofclOffTrack=0; 
      float z=0.; 
      edmNew::DetSetVector<SiPixelCluster>::const_iterator isearch = clustColl.find(detId);
      if( isearch != clustColl.end() ) {  // Not an empty iterator
    edmNew::DetSet<SiPixelCluster>::const_iterator  di;
    for(di=isearch->begin(); di!=isearch->end(); di++){
      unsigned int temp = clusterSet.size();
      clusterSet.insert(*di);
      //check if cluster is off track
      if(clusterSet.size()>temp) {
        otherclusters++;
        nofclOffTrack++; 
        //fill histograms for clusters off tracks
        //correct SiPixelTrackResidualModule
        std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find((*it)->geographicalId().rawId());

        if (pxd!=theSiPixelStructure.end()) (*pxd).second->fill((*di), false, -1., reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn); 
        


        meClSizeNotOnTrack_all->Fill((*di).size());
        meClSizeXNotOnTrack_all->Fill((*di).sizeX());
        meClSizeYNotOnTrack_all->Fill((*di).sizeY());
        meClChargeNotOnTrack_all->Fill((*di).charge()/1000);

        //find cluster global position (rphi, z) get cluster
        //define tracker and pixel geometry and topology
        const TrackerGeometry& theTracker(*theTrackerGeometry);
        const PixelGeomDetUnit* theGeomDet = static_cast<const PixelGeomDetUnit*> (theTracker.idToDet(detId) );
        //test if PixelGeomDetUnit exists
        if(theGeomDet == 0) {
          if(debug_) std::cout << "NO THEGEOMDET\n";
          continue;
        }
        const PixelTopology * topol = &(theGeomDet->specificTopology());
       
        //center of gravity (of charge)
        float xcenter = di->x();
        float ycenter = di->y();
        // get the cluster position in local coordinates (cm) 
        LocalPoint clustlp = topol->localPosition( MeasurementPoint(xcenter, ycenter) );
        // get the cluster position in global coordinates (cm)
        GlobalPoint clustgp = theGeomDet->surface().toGlobal( clustlp );


        //barrel
        if(DetId(detId).subdetId() == 1) {
          meClSizeNotOnTrack_bpix->Fill((*di).size());
          meClSizeXNotOnTrack_bpix->Fill((*di).sizeX());
          meClSizeYNotOnTrack_bpix->Fill((*di).sizeY());
          meClChargeNotOnTrack_bpix->Fill((*di).charge()/1000);
          barrelotherclusters++;
          int DBlayer = PixelBarrelName(DetId(detId), tTopo, isUpgrade).layerName();
          float phi = clustgp.phi(); 
          //float r = clustgp.perp();
          z = clustgp.z();
         for (int i = 0; i < noOfLayers; i++)
         {
           if (DBlayer == i + 1) {
              meClPosLayersNotOnTrack.at(i)->Fill(z,phi);
              meClSizeNotOnTrack_layers.at(i)->Fill((*di).size());
              meClSizeXNotOnTrack_layers.at(i)->Fill((*di).sizeX());
              meClSizeYNotOnTrack_layers.at(i)->Fill((*di).sizeY());
              meClChargeNotOnTrack_layers.at(i)->Fill((*di).charge()/1000);
           }
         }
        }
        //endcap
        if(DetId(detId).subdetId() == 2) {
          meClSizeNotOnTrack_fpix->Fill((*di).size());
          meClSizeXNotOnTrack_fpix->Fill((*di).sizeX());
          meClSizeYNotOnTrack_fpix->Fill((*di).sizeY());
          meClChargeNotOnTrack_fpix->Fill((*di).charge()/1000);
          endcapotherclusters++;
          int DBdisk = PixelEndcapName(DetId(detId), tTopo, isUpgrade).diskName();
          float x = clustgp.x(); 
          float y = clustgp.y(); 
          z = clustgp.z();
          if(z>0){
           for (int i = 0; i < noOfDisks; i++)
           {
             if (DBdisk == i + 1) {
                meClPosDiskspzNotOnTrack.at(i)->Fill(x,y);
                meClSizeNotOnTrack_diskps.at(i)->Fill((*di).size());
                meClSizeXNotOnTrack_diskps.at(i)->Fill((*di).sizeX());
                meClSizeYNotOnTrack_diskps.at(i)->Fill((*di).sizeY());
                meClChargeNotOnTrack_diskps.at(i)->Fill((*di).charge()/1000);
             }
           }
          }
          else{
           for (int i = 0; i < noOfDisks; i++)
           {
             if (DBdisk == i + 1) {
                meClPosDisksmzNotOnTrack.at(i)->Fill(x,y);
                meClSizeNotOnTrack_diskms.at(i)->Fill((*di).size());
                meClSizeXNotOnTrack_diskms.at(i)->Fill((*di).sizeX());
                meClSizeYNotOnTrack_diskms.at(i)->Fill((*di).sizeY());
                meClChargeNotOnTrack_diskms.at(i)->Fill((*di).charge()/1000);
             }
           }
          } 

        }
      }// end "if cluster off track"
      else {
        nofclOnTrack++; 
        if(z == 0){
          //find cluster global position (rphi, z) get cluster
          //define tracker and pixel geometry and topology
          const TrackerGeometry& theTracker(*theTrackerGeometry);
          const PixelGeomDetUnit* theGeomDet = static_cast<const PixelGeomDetUnit*> (theTracker.idToDet(detId) );
          //test if PixelGeomDetUnit exists
          if(theGeomDet == 0) {
        if(debug_) std::cout << "NO THEGEOMDET\n";
        continue;
          }
          const PixelTopology * topol = &(theGeomDet->specificTopology());
          //center of gravity (of charge)
          float xcenter = di->x();
          float ycenter = di->y();
          // get the cluster position in local coordinates (cm) 
          LocalPoint clustlp = topol->localPosition( MeasurementPoint(xcenter, ycenter) );
          // get the cluster position in global coordinates (cm)
          GlobalPoint clustgp = theGeomDet->surface().toGlobal( clustlp );  
          z = clustgp.z();
        }
      }
    }
      }
      //++ fill the number of clusters on a module
      std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find((*it)->geographicalId().rawId());
      if (pxd!=theSiPixelStructure.end()) (*pxd).second->nfill(nofclOnTrack, nofclOffTrack, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn); 
      if(nofclOnTrack!=0) meNClustersOnTrack_all->Fill(nofclOnTrack); 
      if(nofclOffTrack!=0) meNClustersNotOnTrack_all->Fill(nofclOffTrack); 
      //barrel
      if(DetId(detId).subdetId() == 1){
    if(nofclOnTrack!=0) meNClustersOnTrack_bpix->Fill(nofclOnTrack); 
    if(nofclOffTrack!=0) meNClustersNotOnTrack_bpix->Fill(nofclOffTrack); 
    int DBlayer = PixelBarrelName(DetId(detId), tTopo, isUpgrade).layerName();
        for (int i = 0; i < noOfLayers; i++)
        {
          if (DBlayer == i + 1) {
             if(nofclOnTrack!=0) meNClustersOnTrack_layers.at(i)->Fill(nofclOnTrack);
             if(nofclOffTrack!=0) meNClustersNotOnTrack_layers.at(i)->Fill(nofclOffTrack);
          }
        }
      }//end barrel
      //endcap
      if(DetId(detId).subdetId() == 2) {
    int DBdisk = PixelEndcapName(DetId(detId )).diskName();
    //z = clustgp.z();
    if(nofclOnTrack!=0) meNClustersOnTrack_fpix->Fill(nofclOnTrack); 
    if(nofclOffTrack!=0) meNClustersNotOnTrack_fpix->Fill(nofclOffTrack);
    if(z>0){
     for (int i = 0; i < noOfDisks; i++)
     {
       if (DBdisk == i + 1) {
          if(nofclOnTrack!=0) meNClustersOnTrack_diskps.at(i)->Fill(nofclOnTrack);
          if(nofclOffTrack!=0) meNClustersNotOnTrack_diskps.at(i)->Fill(nofclOffTrack);
       }
     }
    }
    if(z<0){
     for (int i = 0; i < noOfDisks; i++)
     {
       if (DBdisk == i + 1) {
         if(nofclOnTrack!=0) meNClustersOnTrack_diskms.at(i)->Fill(nofclOnTrack);
         if(nofclOffTrack!=0) meNClustersNotOnTrack_diskms.at(i)->Fill(nofclOffTrack);
       }
     }
    }
      }

    }//end if it's a Pixel module
  }//end for loop over tracker detector geometry modules

  if(trackclusters>0) (meNofClustersOnTrack_)->Fill(0,trackclusters);
  if(barreltrackclusters>0)(meNofClustersOnTrack_)->Fill(1,barreltrackclusters);
  if(endcaptrackclusters>0)(meNofClustersOnTrack_)->Fill(2,endcaptrackclusters);
  if(otherclusters>0)(meNofClustersNotOnTrack_)->Fill(0,otherclusters);
  if(barrelotherclusters>0)(meNofClustersNotOnTrack_)->Fill(1,barrelotherclusters);
  if(endcapotherclusters>0)(meNofClustersNotOnTrack_)->Fill(2,endcapotherclusters);
  if(tracks>0)(meNofTracks_)->Fill(0,tracks);
  if(pixeltracks>0)(meNofTracks_)->Fill(1,pixeltracks);
  if(bpixtracks>0)(meNofTracks_)->Fill(2,bpixtracks);
  if(fpixtracks>0)(meNofTracks_)->Fill(3,fpixtracks);
}

void SiPixelTrackResidualSource::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  , edm::EventSetup const &  iSetup  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 154 of file SiPixelTrackResidualSource.cc.

References bladeOn, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), clustersrc_, debug_, diskOn, edm::hlt::Exception, firstRun, i, isUpgrade, edm::InputTag::label(), ladOn, layOn, meClChargeNotOnTrack_all, meClChargeNotOnTrack_bpix, meClChargeNotOnTrack_diskms, meClChargeNotOnTrack_diskps, meClChargeNotOnTrack_fpix, meClChargeNotOnTrack_layers, meClChargeOnTrack_all, meClChargeOnTrack_bpix, meClChargeOnTrack_diskms, meClChargeOnTrack_diskps, meClChargeOnTrack_fpix, meClChargeOnTrack_layers, meClPosDisksmzNotOnTrack, meClPosDisksmzOnTrack, meClPosDiskspzNotOnTrack, meClPosDiskspzOnTrack, meClPosLayersNotOnTrack, meClPosLayersOnTrack, meClSizeNotOnTrack_all, meClSizeNotOnTrack_bpix, meClSizeNotOnTrack_diskms, meClSizeNotOnTrack_diskps, meClSizeNotOnTrack_fpix, meClSizeNotOnTrack_layers, meClSizeOnTrack_all, meClSizeOnTrack_bpix, meClSizeOnTrack_diskms, meClSizeOnTrack_diskps, meClSizeOnTrack_fpix, meClSizeOnTrack_layers, meClSizeXNotOnTrack_all, meClSizeXNotOnTrack_bpix, meClSizeXNotOnTrack_diskms, meClSizeXNotOnTrack_diskps, meClSizeXNotOnTrack_fpix, meClSizeXNotOnTrack_layers, meClSizeXOnTrack_all, meClSizeXOnTrack_bpix, meClSizeXOnTrack_diskms, meClSizeXOnTrack_diskps, meClSizeXOnTrack_fpix, meClSizeXOnTrack_layers, meClSizeYNotOnTrack_all, meClSizeYNotOnTrack_bpix, meClSizeYNotOnTrack_diskms, meClSizeYNotOnTrack_diskps, meClSizeYNotOnTrack_fpix, meClSizeYNotOnTrack_layers, meClSizeYOnTrack_all, meClSizeYOnTrack_bpix, meClSizeYOnTrack_diskms, meClSizeYOnTrack_diskps, meClSizeYOnTrack_fpix, meClSizeYOnTrack_layers, meHitProbability, meNClustersNotOnTrack_all, meNClustersNotOnTrack_bpix, meNClustersNotOnTrack_diskms, meNClustersNotOnTrack_diskps, meNClustersNotOnTrack_fpix, meNClustersNotOnTrack_layers, meNClustersOnTrack_all, meNClustersOnTrack_bpix, meNClustersOnTrack_diskms, meNClustersOnTrack_diskps, meNClustersOnTrack_fpix, meNClustersOnTrack_layers, meNofClustersNotOnTrack_, meNofClustersOnTrack_, meNofTracks_, meNofTracksInPixVol_, meResidualXSummedLay, meResidualYSummedLay, meSubdetResidualX, meSubdetResidualY, modOn, noOfDisks, noOfLayers, phiOn, pSet_, reducedSet, ringOn, alignCSCRings::s, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), SiPixelFolderOrganizer::setModuleFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theSiPixelStructure, topFolderName_, and tracksrc_.

                                                                                                                      {
 
  // book residual histograms in theSiPixelFolder - one (x,y) pair of histograms per det
  SiPixelFolderOrganizer theSiPixelFolder(false);
  std::stringstream nameX, titleX, nameY, titleY;

  if(ladOn){
    iBooker.setCurrentFolder(topFolderName_+"/Barrel");
    for (int i = 1; i <= noOfLayers; i++){
      nameX.str(std::string()); nameX <<"siPixelTrackResidualsX_SummedLayer_" << i;
      titleX.str(std::string()); titleX <<"Layer"<< i << "Hit-to-Track Residual in r-phi";
      meResidualXSummedLay.push_back(iBooker.book1D(nameX.str(),titleX.str(),100,-150,150));
      meResidualXSummedLay.at(i-1)->setAxisTitle("hit-to-track residual in r-phi (um)",1);
      nameY.str(std::string()); nameY <<"siPixelTrackResidualsY_SummedLayer_" << i;
      titleY.str(std::string()); titleY <<"Layer"<< i << "Hit-to-Track Residual in Z";
      meResidualYSummedLay.push_back(iBooker.book1D(nameY.str(),titleY.str(),100,-300,300));
      meResidualYSummedLay.at(i-1)->setAxisTitle("hit-to-track residual in z (um)",1);
    }
  }

  for (std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.begin(); 
       pxd!=theSiPixelStructure.end(); pxd++){
   
    if(modOn){
      if (theSiPixelFolder.setModuleFolder(iBooker,(*pxd).first,0,isUpgrade)) (*pxd).second->book(pSet_,iSetup,iBooker,reducedSet,0,isUpgrade);
      else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Folder Creation Failed! "; 
    }
    if(ladOn){
      if (theSiPixelFolder.setModuleFolder(iBooker,(*pxd).first,1,isUpgrade)) {
    
    (*pxd).second->book(pSet_,iSetup,iBooker,reducedSet,1,isUpgrade);
      }
      else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource ladder Folder Creation Failed! "; 
    }
    if(layOn){
      if (theSiPixelFolder.setModuleFolder(iBooker,(*pxd).first,2,isUpgrade)) (*pxd).second->book(pSet_,iSetup,iBooker,reducedSet,2,isUpgrade);
      else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource layer Folder Creation Failed! "; 
    }
    if(phiOn){
      if (theSiPixelFolder.setModuleFolder(iBooker,(*pxd).first,3,isUpgrade)) (*pxd).second->book(pSet_,iSetup,iBooker,reducedSet,3,isUpgrade);
      else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource phi Folder Creation Failed! "; 
    }
    if(bladeOn){
      if (theSiPixelFolder.setModuleFolder(iBooker,(*pxd).first,4,isUpgrade)) (*pxd).second->book(pSet_,iSetup,iBooker,reducedSet,4,isUpgrade);
      else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Blade Folder Creation Failed! "; 
    }
    if(diskOn){
      if (theSiPixelFolder.setModuleFolder(iBooker,(*pxd).first,5,isUpgrade)) (*pxd).second->book(pSet_,iSetup,iBooker,reducedSet,5,isUpgrade);
      else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Disk Folder Creation Failed! "; 
    }
    if(ringOn){
      if (theSiPixelFolder.setModuleFolder(iBooker,(*pxd).first,6,isUpgrade)) (*pxd).second->book(pSet_,iSetup,iBooker,reducedSet,6,isUpgrade);
      else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Ring Folder Creation Failed! "; 
    }
  }

//   edm::InputTag tracksrc = pSet_.getParameter<edm::InputTag>("trajectoryInput");
//   edm::InputTag clustersrc = pSet_.getParameter<edm::InputTag>("clustersrc");

  //number of tracks
  iBooker.setCurrentFolder(topFolderName_+"/Tracks");
  meNofTracks_ = iBooker.book1D("ntracks_" + tracksrc_.label(),"Number of Tracks",4,0,4);
  meNofTracks_->setAxisTitle("Number of Tracks",1);
  meNofTracks_->setBinLabel(1,"All");
  meNofTracks_->setBinLabel(2,"Pixel");
  meNofTracks_->setBinLabel(3,"BPix");
  meNofTracks_->setBinLabel(4,"FPix");

  //number of tracks in pixel fiducial volume
  iBooker.setCurrentFolder(topFolderName_+"/Tracks");
  meNofTracksInPixVol_ = iBooker.book1D("ntracksInPixVol_" + tracksrc_.label(),"Number of Tracks crossing Pixel fiducial Volume",2,0,2);
  meNofTracksInPixVol_->setAxisTitle("Number of Tracks",1);
  meNofTracksInPixVol_->setBinLabel(1,"With Hits");
  meNofTracksInPixVol_->setBinLabel(2,"Without Hits");

  //number of clusters (associated to track / not associated)
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OnTrack");
  meNofClustersOnTrack_ = iBooker.book1D("nclusters_" + clustersrc_.label() + "_tot","Number of Clusters (on track)",3,0,3);
  meNofClustersOnTrack_->setAxisTitle("Number of Clusters on Track",1);
  meNofClustersOnTrack_->setBinLabel(1,"All");
  meNofClustersOnTrack_->setBinLabel(2,"BPix");
  meNofClustersOnTrack_->setBinLabel(3,"FPix");
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OffTrack");
  meNofClustersNotOnTrack_ = iBooker.book1D("nclusters_" + clustersrc_.label() + "_tot","Number of Clusters (off track)",3,0,3);
  meNofClustersNotOnTrack_->setAxisTitle("Number of Clusters off Track",1);
  meNofClustersNotOnTrack_->setBinLabel(1,"All");
  meNofClustersNotOnTrack_->setBinLabel(2,"BPix");
  meNofClustersNotOnTrack_->setBinLabel(3,"FPix");

  //cluster charge and size
  //charge
  //on track
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OnTrack");
  std::stringstream ss1, ss2;
  meClChargeOnTrack_all = iBooker.book1D("charge_" + clustersrc_.label(),"Charge (on track)",500,0.,500.);
  meClChargeOnTrack_all->setAxisTitle("Charge size (in ke)",1);
  meClChargeOnTrack_bpix = iBooker.book1D("charge_" + clustersrc_.label() + "_Barrel","Charge (on track, barrel)",500,0.,500.);
  meClChargeOnTrack_bpix->setAxisTitle("Charge size (in ke)",1);
  meClChargeOnTrack_fpix = iBooker.book1D("charge_" + clustersrc_.label() + "_Endcap","Charge (on track, endcap)",500,0.,500.);
  meClChargeOnTrack_fpix->setAxisTitle("Charge size (in ke)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "charge_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Charge (on track, layer" << i << ")";
    meClChargeOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClChargeOnTrack_layers.at(i-1)->setAxisTitle("Charge size (in ke)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "charge_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "Charge (on track, diskp" << i << ")";
    meClChargeOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClChargeOnTrack_diskps.at(i-1)->setAxisTitle("Charge size (in ke)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "charge_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "Charge (on track, diskm" << i << ")";
    meClChargeOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClChargeOnTrack_diskms.at(i-1)->setAxisTitle("Charge size (in ke)",1);
  }
  //off track
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OffTrack");
  meClChargeNotOnTrack_all = iBooker.book1D("charge_" + clustersrc_.label(),"Charge (off track)",500,0.,500.);
  meClChargeNotOnTrack_all->setAxisTitle("Charge size (in ke)",1);
  meClChargeNotOnTrack_bpix = iBooker.book1D("charge_" + clustersrc_.label() + "_Barrel","Charge (off track, barrel)",500,0.,500.);
  meClChargeNotOnTrack_bpix->setAxisTitle("Charge size (in ke)",1);
  meClChargeNotOnTrack_fpix = iBooker.book1D("charge_" + clustersrc_.label() + "_Endcap","Charge (off track, endcap)",500,0.,500.);
  meClChargeNotOnTrack_fpix->setAxisTitle("Charge size (in ke)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "charge_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Charge (off track, layer" << i << ")";
    meClChargeNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClChargeNotOnTrack_layers.at(i-1)->setAxisTitle("Charge size (in ke)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "charge_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "Charge (off track, diskp" << i << ")";
    meClChargeNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClChargeNotOnTrack_diskps.at(i-1)->setAxisTitle("Charge size (in ke)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "charge_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "Charge (off track, diskm" << i << ")";
    meClChargeNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClChargeNotOnTrack_diskms.at(i-1)->setAxisTitle("Charge size (in ke)",1);
  }

  //size
  //on track
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OnTrack");
  meClSizeOnTrack_all = iBooker.book1D("size_" + clustersrc_.label(),"Size (on track)",100,0.,100.);
  meClSizeOnTrack_all->setAxisTitle("Cluster size (in pixels)",1);
  meClSizeOnTrack_bpix = iBooker.book1D("size_" + clustersrc_.label() + "_Barrel","Size (on track, barrel)",100,0.,100.);
  meClSizeOnTrack_bpix->setAxisTitle("Cluster size (in pixels)",1);
  meClSizeOnTrack_fpix = iBooker.book1D("size_" + clustersrc_.label() + "_Endcap","Size (on track, endcap)",100,0.,100.);
  meClSizeOnTrack_fpix->setAxisTitle("Cluster size (in pixels)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "size_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Size (on track, layer" << i << ")";
    meClSizeOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeOnTrack_layers.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "size_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "Size (on track, diskp" << i << ")";
    meClSizeOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeOnTrack_diskps.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "size_" + clustersrc_.label() + "_Disk_m1" << i;
    ss2.str(std::string()); ss2 << "Size (on track, diskm" << i << ")";
    meClSizeOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeOnTrack_diskms.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  meClSizeXOnTrack_all = iBooker.book1D("sizeX_" + clustersrc_.label(),"SizeX (on track)",100,0.,100.);
  meClSizeXOnTrack_all->setAxisTitle("Cluster sizeX (in pixels)",1);
  meClSizeXOnTrack_bpix = iBooker.book1D("sizeX_" + clustersrc_.label() + "_Barrel","SizeX (on track, barrel)",100,0.,100.);
  meClSizeXOnTrack_bpix->setAxisTitle("Cluster sizeX (in pixels)",1);
  meClSizeXOnTrack_fpix = iBooker.book1D("sizeX_" + clustersrc_.label() + "_Endcap","SizeX (on track, endcap)",100,0.,100.);
  meClSizeXOnTrack_fpix->setAxisTitle("Cluster sizeX (in pixels)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "sizeX_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "SizeX (on track, layer" << i << ")";
    meClSizeXOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeXOnTrack_layers.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeX_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "SizeX (on track, diskp" << i << ")";
    meClSizeXOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeXOnTrack_diskps.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeX_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "SizeX (on track, diskm" << i << ")";
    meClSizeXOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeXOnTrack_diskms.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  meClSizeYOnTrack_all = iBooker.book1D("sizeY_" + clustersrc_.label(),"SizeY (on track)",100,0.,100.);
  meClSizeYOnTrack_all->setAxisTitle("Cluster sizeY (in pixels)",1);
  meClSizeYOnTrack_bpix = iBooker.book1D("sizeY_" + clustersrc_.label() + "_Barrel","SizeY (on track, barrel)",100,0.,100.);
  meClSizeYOnTrack_bpix->setAxisTitle("Cluster sizeY (in pixels)",1);
  meClSizeYOnTrack_fpix = iBooker.book1D("sizeY_" + clustersrc_.label() + "_Endcap","SizeY (on track, endcap)",100,0.,100.);
  meClSizeYOnTrack_fpix->setAxisTitle("Cluster sizeY (in pixels)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "sizeY_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "SizeY (on track, layer" << i << ")";
    meClSizeYOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeYOnTrack_layers.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeY_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "SizeY (on track, diskp" << i << ")";
    meClSizeYOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeYOnTrack_diskps.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeY_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "SizeY (on track, diskm" << i << ")";
    meClSizeYOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeYOnTrack_diskms.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  //off track
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OffTrack");
  meClSizeNotOnTrack_all = iBooker.book1D("size_" + clustersrc_.label(),"Size (off track)",100,0.,100.);
  meClSizeNotOnTrack_all->setAxisTitle("Cluster size (in pixels)",1); 
  meClSizeNotOnTrack_bpix = iBooker.book1D("size_" + clustersrc_.label() + "_Barrel","Size (off track, barrel)",100,0.,100.);
  meClSizeNotOnTrack_bpix->setAxisTitle("Cluster size (in pixels)",1);
  meClSizeNotOnTrack_fpix = iBooker.book1D("size_" + clustersrc_.label() + "_Endcap","Size (off track, endcap)",100,0.,100.);
  meClSizeNotOnTrack_fpix->setAxisTitle("Cluster size (in pixels)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "size_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Size (off track, layer" << i << ")";
    meClSizeNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeNotOnTrack_layers.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "size_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "Size (off track, diskp" << i << ")";
    meClSizeNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeNotOnTrack_diskps.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "size_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "Size (off track, diskm" << i << ")";
    meClSizeNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeNotOnTrack_diskms.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  meClSizeXNotOnTrack_all = iBooker.book1D("sizeX_" + clustersrc_.label(),"SizeX (off track)",100,0.,100.);
  meClSizeXNotOnTrack_all->setAxisTitle("Cluster sizeX (in pixels)",1);
  meClSizeXNotOnTrack_bpix = iBooker.book1D("sizeX_" + clustersrc_.label() + "_Barrel","SizeX (off track, barrel)",100,0.,100.);
  meClSizeXNotOnTrack_bpix->setAxisTitle("Cluster sizeX (in pixels)",1);
  meClSizeXNotOnTrack_fpix = iBooker.book1D("sizeX_" + clustersrc_.label() + "_Endcap","SizeX (off track, endcap)",100,0.,100.);
  meClSizeXNotOnTrack_fpix->setAxisTitle("Cluster sizeX (in pixels)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "sizeX_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "SizeX (off track, layer" << i << ")";
    meClSizeXNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeXNotOnTrack_layers.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeX_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "SizeX (off track, diskp" << i << ")";
    meClSizeXNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeXNotOnTrack_diskps.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeX_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "SizeX (off track, diskm" << i << ")";
    meClSizeXNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeXNotOnTrack_diskms.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  meClSizeYNotOnTrack_all = iBooker.book1D("sizeY_" + clustersrc_.label(),"SizeY (off track)",100,0.,100.);
  meClSizeYNotOnTrack_all->setAxisTitle("Cluster sizeY (in pixels)",1);
  meClSizeYNotOnTrack_bpix = iBooker.book1D("sizeY_" + clustersrc_.label() + "_Barrel","SizeY (off track, barrel)",100,0.,100.);
  meClSizeYNotOnTrack_bpix->setAxisTitle("Cluster sizeY (in pixels)",1);
  meClSizeYNotOnTrack_fpix = iBooker.book1D("sizeY_" + clustersrc_.label() + "_Endcap","SizeY (off track, endcap)",100,0.,100.);
  meClSizeYNotOnTrack_fpix->setAxisTitle("Cluster sizeY (in pixels)",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "sizeY_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "SizeY (off track, layer" << i << ")";
    meClSizeYNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeYNotOnTrack_layers.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeY_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "SizeY (off track, diskp" << i << ")";
    meClSizeYNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeYNotOnTrack_diskps.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "sizeY_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "SizeY (off track, diskm" << i << ")";
    meClSizeYNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meClSizeYNotOnTrack_diskms.at(i-1)->setAxisTitle("Cluster size (in pixels)",1);
  }

  //cluster global position
  //on track
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OnTrack");
  //bpix
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "position_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Clusters Layer" << i << " (on track)";
    meClPosLayersOnTrack.push_back(iBooker.book2D(ss1.str(),ss2.str(),200,-30.,30.,128,-3.2,3.2));
    meClPosLayersOnTrack.at(i-1)->setAxisTitle("Global Z (cm)",1);
    meClPosLayersOnTrack.at(i-1)->setAxisTitle("Global #phi",2);
  }
  //fpix
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "position_" + clustersrc_.label() + "_pz_Disk_" << i;
    ss2.str(std::string()); ss2 << "Clusters +Z Disk" << i << " (on track)";
    meClPosDiskspzOnTrack.push_back(iBooker.book2D(ss1.str(),ss2.str(),80,-20.,20.,80,-20.,20.));
    meClPosDiskspzOnTrack.at(i-1)->setAxisTitle("Global X (cm)",1);
    meClPosDiskspzOnTrack.at(i-1)->setAxisTitle("Global Y (cm)",2);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "position_" + clustersrc_.label() + "_mz_Disk_" << i;
    ss2.str(std::string()); ss2 << "Clusters -Z Disk" << i << " (on track)";
    meClPosDisksmzOnTrack.push_back(iBooker.book2D(ss1.str(),ss2.str(),80,-20.,20.,80,-20.,20.));
    meClPosDisksmzOnTrack.at(i-1)->setAxisTitle("Global X (cm)",1);
    meClPosDisksmzOnTrack.at(i-1)->setAxisTitle("Global Y (cm)",2);
  }
  meNClustersOnTrack_all = iBooker.book1D("nclusters_" + clustersrc_.label(),"Number of Clusters (on Track)",50,0.,50.);
  meNClustersOnTrack_all->setAxisTitle("Number of Clusters",1);
  meNClustersOnTrack_bpix = iBooker.book1D("nclusters_" + clustersrc_.label() + "_Barrel","Number of Clusters (on track, barrel)",50,0.,50.);
  meNClustersOnTrack_bpix->setAxisTitle("Number of Clusters",1);
  meNClustersOnTrack_fpix = iBooker.book1D("nclusters_" + clustersrc_.label() + "_Endcap","Number of Clusters (on track, endcap)",50,0.,50.);
  meNClustersOnTrack_fpix->setAxisTitle("Number of Clusters",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "nclusters_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Number of Clusters (on track, layer" << i << ")";
    meNClustersOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meNClustersOnTrack_layers.at(i-1)->setAxisTitle("Number of Clusters",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "nclusters_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "Number of Clusters (on track, diskp" << i << ")";
    meNClustersOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),50,0.,50.));
    meNClustersOnTrack_diskps.at(i-1)->setAxisTitle("Number of Clusters",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "nclusters_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "Number of Clusters (on track, diskm" << i << ")";
    meNClustersOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meNClustersOnTrack_diskms.at(i-1)->setAxisTitle("Number of Clusters",1);
  }

  //not on track
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OffTrack");
  //bpix
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "position_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Clusters Layer" << i << " (off track)";
    meClPosLayersNotOnTrack.push_back(iBooker.book2D(ss1.str(),ss2.str(),200,-30.,30.,128,-3.2,3.2));
    meClPosLayersNotOnTrack.at(i-1)->setAxisTitle("Global Z (cm)",1);
    meClPosLayersNotOnTrack.at(i-1)->setAxisTitle("Global #phi",2);
  }
  //fpix
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "position_" + clustersrc_.label() + "_pz_Disk_" << i;
    ss2.str(std::string()); ss2 << "Clusters +Z Disk" << i << " (off track)";
    meClPosDiskspzNotOnTrack.push_back(iBooker.book2D(ss1.str(),ss2.str(),80,-20.,20.,80,-20.,20.));
    meClPosDiskspzNotOnTrack.at(i-1)->setAxisTitle("Global X (cm)",1);
    meClPosDiskspzNotOnTrack.at(i-1)->setAxisTitle("Global Y (cm)",2);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "position_" + clustersrc_.label() + "_mz_Disk_" << i;
    ss2.str(std::string()); ss2 << "Clusters -Z Disk" << i << " (off track)";
    meClPosDisksmzNotOnTrack.push_back(iBooker.book2D(ss1.str(),ss2.str(),80,-20.,20.,80,-20.,20.));
    meClPosDisksmzNotOnTrack.at(i-1)->setAxisTitle("Global X (cm)",1);
    meClPosDisksmzNotOnTrack.at(i-1)->setAxisTitle("Global Y (cm)",2);
  }
  meNClustersNotOnTrack_all = iBooker.book1D("nclusters_" + clustersrc_.label(),"Number of Clusters (off Track)",50,0.,50.);
  meNClustersNotOnTrack_all->setAxisTitle("Number of Clusters",1);
  meNClustersNotOnTrack_bpix = iBooker.book1D("nclusters_" + clustersrc_.label() + "_Barrel","Number of Clusters (off track, barrel)",50,0.,50.);
  meNClustersNotOnTrack_bpix->setAxisTitle("Number of Clusters",1);
  meNClustersNotOnTrack_fpix = iBooker.book1D("nclusters_" + clustersrc_.label() + "_Endcap","Number of Clusters (off track, endcap)",50,0.,50.);
  meNClustersNotOnTrack_fpix->setAxisTitle("Number of Clusters",1);
  for (int i = 1; i <= noOfLayers; i++)
  {
    ss1.str(std::string()); ss1 << "nclusters_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string()); ss2 << "Number of Clusters (off track, layer" << i << ")";
    meNClustersNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meNClustersNotOnTrack_layers.at(i-1)->setAxisTitle("Number of Clusters",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "nclusters_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string()); ss2 << "Number of Clusters (off track, diskp" << i << ")";
    meNClustersNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meNClustersNotOnTrack_diskps.at(i-1)->setAxisTitle("Number of Clusters",1);
  }
  for (int i = 1; i <= noOfDisks; i++)
  {
    ss1.str(std::string()); ss1 << "nclusters_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string()); ss2 << "Number of Clusters (off track, diskm" << i << ")";
    meNClustersNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(),ss2.str(),500,0.,500.));
    meNClustersNotOnTrack_diskms.at(i-1)->setAxisTitle("Number of Clusters",1);
  
  }
  //HitProbability
  //on track
  iBooker.setCurrentFolder(topFolderName_+"/Clusters/OnTrack");
  meHitProbability = iBooker.book1D("FractionLowProb","Fraction of hits with low probability;FractionLowProb;#HitsOnTrack",100,0.,1.);

  if (debug_) {
    // book summary residual histograms in a debugging folder - one (x,y) pair of histograms per subdetector 
    iBooker.setCurrentFolder("debugging"); 
    char hisID[80]; 
    for (int s=0; s<3; s++) {
      sprintf(hisID,"residual_x_subdet_%i",s); 
      meSubdetResidualX[s] = iBooker.book1D(hisID,"Pixel Hit-to-Track Residual in X",500,-5.,5.);
      
      sprintf(hisID,"residual_y_subdet_%i",s);
      meSubdetResidualY[s] = iBooker.book1D(hisID,"Pixel Hit-to-Track Residual in Y",500,-5.,5.);  
    }
  }
  
  firstRun = false;
}

void SiPixelTrackResidualSource::dqmBeginRun ( const edm::Run &  r, edm::EventSetup const &  iSetup  )

virtual

Reimplemented from DQMEDAnalyzer.

Definition at line 120 of file SiPixelTrackResidualSource.cc.

References debug_, PixelEndcapName::diskName(), edm::EventSetup::get(), isUpgrade, PixelBarrelName::layerName(), python.rootplot.argparse::module, noOfDisks, noOfLayers, edm::ESHandle< class >::product(), and theSiPixelStructure.

                                                                                         {
  LogInfo("PixelDQM") << "SiPixelTrackResidualSource beginRun()" << endl;
  // retrieve TrackerGeometry for pixel dets
  edm::ESHandle<TrackerGeometry> TG;
  iSetup.get<TrackerDigiGeometryRecord>().get(TG);
  if (debug_) LogVerbatim("PixelDQM") << "TrackerGeometry "<< &(*TG) <<" size is "<< TG->dets().size() << endl;
  edm::ESHandle<TrackerTopology> tTopoHandle;
  iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
  const TrackerTopology *pTT = tTopoHandle.product();

  // build theSiPixelStructure with the pixel barrel and endcap dets from TrackerGeometry
  for (TrackerGeometry::DetContainer::const_iterator pxb = TG->detsPXB().begin();  
       pxb!=TG->detsPXB().end(); pxb++) {
    if (dynamic_cast<PixelGeomDetUnit const *>((*pxb))!=0) {
      SiPixelTrackResidualModule* module = new SiPixelTrackResidualModule((*pxb)->geographicalId().rawId());
      theSiPixelStructure.insert(pair<uint32_t, SiPixelTrackResidualModule*>((*pxb)->geographicalId().rawId(), module));
      //int DBlayer = PixelBarrelNameWrapper(pSet_, DetId((*pxb)->geographicalId())).layerName();
      int DBlayer = PixelBarrelName(DetId((*pxb)->geographicalId()),pTT,isUpgrade).layerName();
      if (noOfLayers < DBlayer) noOfLayers = DBlayer;
    }
  }
  for (TrackerGeometry::DetContainer::const_iterator pxf = TG->detsPXF().begin(); 
       pxf!=TG->detsPXF().end(); pxf++) {
    if (dynamic_cast<PixelGeomDetUnit const *>((*pxf))!=0) {
      SiPixelTrackResidualModule* module = new SiPixelTrackResidualModule((*pxf)->geographicalId().rawId());
      theSiPixelStructure.insert(pair<uint32_t, SiPixelTrackResidualModule*>((*pxf)->geographicalId().rawId(), module));
      int DBdisk;
      DBdisk = PixelEndcapName(DetId((*pxf)->geographicalId()),pTT,isUpgrade).diskName();
      if (noOfDisks < DBdisk) noOfDisks = DBdisk;
    }
  }
  LogInfo("PixelDQM") << "SiPixelStructure size is " << theSiPixelStructure.size() << endl;
}

void SiPixelTrackResidualSource::triplets	(	double	x1,
		double	y1,
		double	z1,
		double	x2,
		double	y2,
		double	z2,
		double	x3,
		double	y3,
		double	z3,
		double	ptsig,
		double &	dc,
		double &	dz,
		double	kap
	)

Definition at line 1302 of file SiPixelTrackResidualSource.cc.

References gather_cfg::cout, dttmaxenums::L, pileupDistInMC::num, diffTwoXMLs::r1, rho, and mathSSE::sqrt().

Referenced by analyze().

                                                                            {
  
  //Define some constants
  using namespace std;

  //Curvature kap from global Track

 //inverse of the curvature is the radius in the transverse plane
  double rho = 1/kap;
  //Check that the hits are in the correct layers
  double r1 = sqrt( x1*x1 + y1*y1 );
  double r3 = sqrt( x3*x3 + y3*y3 );

  if( r3-r1 < 2.0 ) cout << "warn r1 = " << r1 << ", r3 = " << r3 << endl;
  
  // Calculate the centre of the helix in xy-projection with radius rho from the track.
  //start with a line (sekante) connecting the two points (x1,y1) and (x3,y3) vec_L = vec_x3-vec_x1
  //with L being the length of that vector.
  double L=sqrt((x3-x1)*(x3-x1)+(y3-y1)*(y3-y1));
  //lam is the line from the middel point of vec_q towards the center of the circle X0,Y0
  // we already have kap and rho = 1/kap
  double lam = sqrt(rho*rho - L*L/4);
  
  // There are two solutions, the sign of kap gives the information
  // which of them is correct.
  //
  if( kap > 0 ) lam = -lam;

  //
  // ( X0, Y0 ) is the centre of the circle that describes the helix in xy-projection.
  //
  double x0 =  0.5*( x1 + x3 ) + lam/L * ( -y1 + y3 );
  double y0 =  0.5*( y1 + y3 ) + lam/L * (  x1 - x3 );

  // Calculate the dipangle in z direction (needed later for z residual) :
  //Starting from the heliz equation whihc has to hold for both points z1,z3
  double num = ( y3 - y0 ) * ( x1 - x0 ) - ( x3 - x0 ) * ( y1 - y0 );
  double den = ( x1 - x0 ) * ( x3 - x0 ) + ( y1 - y0 ) * ( y3 - y0 );
  double tandip = kap * ( z3 - z1 ) / atan( num / den );

 
  // angle from first hit to dca point:
  //
  double dphi = atan( ( ( x1 - x0 ) * y0 - ( y1 - y0 ) * x0 )
                      / ( ( x1 - x0 ) * x0 + ( y1 - y0 ) * y0 ) );
  //z position of the track based on the middle of the circle
  //track equation for the z component
  double uz0 = z1 + tandip * dphi * rho;

  //RESIDUAL IN R-PHI
  //Calculate distance dca2 from point (x2,y2) to the circle which is given by
  //the distance of the point to the middlepoint dcM = sqrt((x0-x2)^2+(y0-y2)) and rho
  //dca = rho +- dcM
  if(kap>0) dca2=rho-sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2));
  else dca2=rho+sqrt((-x0+x2)*(-x0+x2)+(-y0+y2)*(-y0+y2));

  //RESIDUAL IN Z
  double xx =0 ;
  double yy =0 ;
  //sign of kappa determines the calculation
  //xx and yy are the new coordinates starting from x2, y2 that are on the track itself
  //vec_X2+-dca2*vec(X0-X2)/|(X0-X2)|
  if(kap<0){
    xx =   x2+(dca2*((x0-x2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
    yy = y2+(dca2*((y0-y2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
  }
  else if(kap>=0){
    xx =   x2-(dca2*((x0-x2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
    yy = y2-(dca2*((y0-y2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
  }

  //to get residual in z start with calculating the new uz2 position if one has moved to xx, yy
  //on the track. First calculate the change in phi2 with respect to the center X0, Y0
  double dphi2 = atan( ( ( xx - x0 ) * y0 - ( yy - y0 ) * x0 )
               / ( ( xx - x0 ) * x0 + ( yy - y0 ) * y0 ) );
  //Solve track equation for this new z depending on the dip angle of the track (see above
  //calculated based on X1, X3 and X0, use uz0 as reference point again.
  double  uz2= uz0 - dphi2*tandip*rho;
  
  //subtract new z position from the old one
  dz2=z2-uz2;
  
  //if we are interested in the arclength this is unsigned though
  //  double cosphi2 = (x2*xx+y2*yy)/(sqrt(x2*x2+y2*y2)*sqrt(xx*xx+yy*yy));
  //double arcdca2=sqrt(x2*x2+y2*y2)*acos(cosphi2);
  
}

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> SiPixelTrackResidualSource::beamSpotToken_

private

Definition at line 67 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::bladeOn

private

Definition at line 82 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and SiPixelTrackResidualSource().

edm::InputTag SiPixelTrackResidualSource::clustersrc_

private

Definition at line 64 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > SiPixelTrackResidualSource::clustersrcToken_

private

Definition at line 73 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::debug_

private

Definition at line 76 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), dqmBeginRun(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::diskOn

private

Definition at line 82 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::firstRun

private

Definition at line 85 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

edm::EDGetTokenT<reco::TrackCollection> SiPixelTrackResidualSource::generalTracksToken_

private

Definition at line 69 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::isUpgrade

private

Definition at line 83 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

bool SiPixelTrackResidualSource::ladOn

private

Definition at line 80 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::layOn

private

Definition at line 80 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and SiPixelTrackResidualSource().

MonitorElement* SiPixelTrackResidualSource::meClChargeNotOnTrack_all

private

Definition at line 107 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClChargeNotOnTrack_bpix

private

Definition at line 108 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClChargeNotOnTrack_diskms

private

Definition at line 112 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClChargeNotOnTrack_diskps

private

Definition at line 111 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClChargeNotOnTrack_fpix

private

Definition at line 109 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClChargeNotOnTrack_layers

private

Definition at line 110 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClChargeOnTrack_all

private

Definition at line 101 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClChargeOnTrack_bpix

private

Definition at line 102 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClChargeOnTrack_diskms

private

Definition at line 106 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClChargeOnTrack_diskps

private

Definition at line 105 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClChargeOnTrack_fpix

private

Definition at line 103 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClChargeOnTrack_layers

private

Definition at line 104 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClPosDisksmzNotOnTrack

private

Definition at line 170 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClPosDisksmzOnTrack

private

Definition at line 168 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClPosDiskspzNotOnTrack

private

Definition at line 169 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClPosDiskspzOnTrack

private

Definition at line 167 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClPosLayersNotOnTrack

private

Definition at line 166 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClPosLayersOnTrack

private

Definition at line 165 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeNotOnTrack_all

private

Definition at line 119 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeNotOnTrack_bpix

private

Definition at line 120 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeNotOnTrack_diskms

private

Definition at line 124 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeNotOnTrack_diskps

private

Definition at line 123 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeNotOnTrack_fpix

private

Definition at line 121 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeNotOnTrack_layers

private

Definition at line 122 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeOnTrack_all

private

Definition at line 113 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeOnTrack_bpix

private

Definition at line 114 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeOnTrack_diskms

private

Definition at line 118 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeOnTrack_diskps

private

Definition at line 117 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeOnTrack_fpix

private

Definition at line 115 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeOnTrack_layers

private

Definition at line 116 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeXNotOnTrack_all

private

Definition at line 131 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeXNotOnTrack_bpix

private

Definition at line 132 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeXNotOnTrack_diskms

private

Definition at line 136 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeXNotOnTrack_diskps

private

Definition at line 135 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeXNotOnTrack_fpix

private

Definition at line 133 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeXNotOnTrack_layers

private

Definition at line 134 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeXOnTrack_all

private

Definition at line 125 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeXOnTrack_bpix

private

Definition at line 126 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeXOnTrack_diskms

private

Definition at line 130 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeXOnTrack_diskps

private

Definition at line 129 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeXOnTrack_fpix

private

Definition at line 127 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeXOnTrack_layers

private

Definition at line 128 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeYNotOnTrack_all

private

Definition at line 143 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeYNotOnTrack_bpix

private

Definition at line 144 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeYNotOnTrack_diskms

private

Definition at line 148 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeYNotOnTrack_diskps

private

Definition at line 147 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeYNotOnTrack_fpix

private

Definition at line 145 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeYNotOnTrack_layers

private

Definition at line 146 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeYOnTrack_all

private

Definition at line 137 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeYOnTrack_bpix

private

Definition at line 138 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeYOnTrack_diskms

private

Definition at line 142 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeYOnTrack_diskps

private

Definition at line 141 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meClSizeYOnTrack_fpix

private

Definition at line 139 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meClSizeYOnTrack_layers

private

Definition at line 140 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meHitProbability

private

Definition at line 172 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNClustersNotOnTrack_all

private

Definition at line 157 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNClustersNotOnTrack_bpix

private

Definition at line 158 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meNClustersNotOnTrack_diskms

private

Definition at line 162 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meNClustersNotOnTrack_diskps

private

Definition at line 161 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNClustersNotOnTrack_fpix

private

Definition at line 159 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meNClustersNotOnTrack_layers

private

Definition at line 160 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNClustersOnTrack_all

private

Definition at line 151 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNClustersOnTrack_bpix

private

Definition at line 152 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meNClustersOnTrack_diskms

private

Definition at line 156 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meNClustersOnTrack_diskps

private

Definition at line 155 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNClustersOnTrack_fpix

private

Definition at line 153 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meNClustersOnTrack_layers

private

Definition at line 154 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNofClustersNotOnTrack_

private

Definition at line 100 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNofClustersOnTrack_

private

Definition at line 99 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNofTracks_

private

Definition at line 97 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meNofTracksInPixVol_

private

Definition at line 98 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meResidualXSummedLay

private

Definition at line 94 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> SiPixelTrackResidualSource::meResidualYSummedLay

private

Definition at line 95 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meSubdetResidualX[3]

private

Definition at line 91 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms().

MonitorElement* SiPixelTrackResidualSource::meSubdetResidualY[3]

private

Definition at line 92 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms().

bool SiPixelTrackResidualSource::modOn

private

Definition at line 77 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and SiPixelTrackResidualSource().

int SiPixelTrackResidualSource::NLowProb

private

Definition at line 87 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

int SiPixelTrackResidualSource::noOfDisks

private

Definition at line 175 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

int SiPixelTrackResidualSource::noOfLayers

private

Definition at line 174 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

int SiPixelTrackResidualSource::NTotal

private

Definition at line 86 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

edm::EDGetTokenT<reco::VertexCollection> SiPixelTrackResidualSource::offlinePrimaryVerticesToken_

private

Definition at line 68 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::phiOn

private

Definition at line 80 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and SiPixelTrackResidualSource().

edm::ParameterSet SiPixelTrackResidualSource::pSet_

private

Definition at line 62 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

double SiPixelTrackResidualSource::ptminres_

private

Definition at line 84 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

bool SiPixelTrackResidualSource::reducedSet

private

Definition at line 78 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

bool SiPixelTrackResidualSource::ringOn

private

Definition at line 82 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), and SiPixelTrackResidualSource().

edm::InputTag SiPixelTrackResidualSource::src_

private

Definition at line 63 of file SiPixelTrackResidualSource.h.

Referenced by SiPixelTrackResidualSource().

std::map<uint32_t, SiPixelTrackResidualModule*> SiPixelTrackResidualSource::theSiPixelStructure

private

Definition at line 89 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), dqmBeginRun(), and ~SiPixelTrackResidualSource().

std::string SiPixelTrackResidualSource::topFolderName_

Definition at line 59 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

edm::EDGetTokenT<TrajTrackAssociationCollection> SiPixelTrackResidualSource::trackAssociationToken_

private

Definition at line 72 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

edm::InputTag SiPixelTrackResidualSource::tracksrc_

private

Definition at line 65 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

edm::EDGetTokenT<std::vector<Trajectory> > SiPixelTrackResidualSource::tracksrcToken_

private

Definition at line 70 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

edm::EDGetTokenT<std::vector<reco::Track> > SiPixelTrackResidualSource::trackToken_

private

Definition at line 71 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

std::string SiPixelTrackResidualSource::ttrhbuilder_

private

Definition at line 66 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().