analyzer::SiPixelLorentzAngle Class Reference

#include <SiPixelLorentzAngle.h>

Inheritance diagram for analyzer::SiPixelLorentzAngle:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginJob () override
void	endJob () override
	SiPixelLorentzAngle (const edm::ParameterSet &conf)
	~SiPixelLorentzAngle () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	fillPix (const SiPixelCluster &LocPix, const PixelTopology *topol, Pixinfo &pixinfo)
void	findMean (int i, int i_ring)

Private Attributes
std::map< int, TH2F * >	_h_drift_depth_
std::map< int, TH2F * >	_h_drift_depth_adc2_
std::map< int, TH2F * >	_h_drift_depth_adc_
std::map< int, TH2F * >	_h_drift_depth_noadc_
std::map< int, TH1F * >	_h_mean_
int	bladeF_
double	chi2_
double	chi2F_
Clust	clust_
double	clustChargeMax_
Clust	clustF_
int	clustSizeYMin_
int	diskF_
float	eta_
float	etaF_
int	event_
int	event_counter_
int	eventF_
std::string	filename_
std::string	filenameFit_
TH2F *	h_cluster_shape_
TH2F *	h_cluster_shape_adc_
TH2F *	h_cluster_shape_adc_rot_
TH2F *	h_cluster_shape_noadc_
TH2F *	h_cluster_shape_noadc_rot_
TH2F *	h_cluster_shape_rot_
TH1F *	h_drift_depth_adc_slice_
TH1F *	h_tracks_
TFile *	hFile_
int	hist_depth_
int	hist_drift_
int	hist_x_
int	hist_y_
int	hitCounter_
int	isflipped_
int	ladder_
int	layer_
double	max_depth_
double	max_drift_
double	max_x_
double	max_y_
double	min_depth_
double	min_drift_
double	min_x_
double	min_y_
int	module_
int	moduleF_
double	ndof_
double	ndofF_
double	normChi2Max_
int	panelF_
float	phi_
float	phiF_
int	pixelTracksCounter_
Pixinfo	pixinfo_
Pixinfo	pixinfoF_
float	pt_
float	ptF_
double	ptmin_
Rechit	rechit_
Rechit	rechitF_
double	residualMax_
const TransientTrackingRecHitBuilder *	RHBuilder
int	run_
int	runF_
int	sideF_
bool	simData_
Hit	simhit_
Hit	simhitF_
TTree *	SiPixelLorentzAngleTree_
TTree *	SiPixelLorentzAngleTreeForward_
edm::EDGetTokenT< TrajTrackAssociationCollection >	t_trajTrack
Chi2MeasurementEstimator *	theEstimator
const KFTrajectoryFitter *	theFitter
PropagatorWithMaterial *	thePropagator
PropagatorWithMaterial *	thePropagatorOp
const KFTrajectorySmoother *	theSmoother
KFUpdator *	theUpdator
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomToken_
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopoToken_
int	trackEventsCounter_
Hit	trackhit_
Hit	trackhitF_
TrajectoryStateTransform	tsTransform
int	usedHitCounter_
double	width_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 83 of file SiPixelLorentzAngle.h.

Constructor & Destructor Documentation

SiPixelLorentzAngle()

SiPixelLorentzAngle::SiPixelLorentzAngle ( const edm::ParameterSet &  conf )

explicit

Definition at line 38 of file SiPixelLorentzAngle.cc.

References edm::ParameterSet::getParameter(), hist_x_, hist_y_, max_depth_, max_drift_, max_x_, max_y_, min_depth_, min_drift_, min_x_, min_y_, t_trajTrack, trackerGeomToken_, trackerTopoToken_, and width_.

    : filename_(conf.getParameter<std::string>("fileName")),
      filenameFit_(conf.getParameter<std::string>("fileNameFit")),
      ptmin_(conf.getParameter<double>("ptMin")),
      simData_(conf.getParameter<bool>("simData")),
      normChi2Max_(conf.getParameter<double>("normChi2Max")),
      clustSizeYMin_(conf.getParameter<int>("clustSizeYMin")),
      residualMax_(conf.getParameter<double>("residualMax")),
      clustChargeMax_(conf.getParameter<double>("clustChargeMax")),
      hist_depth_(conf.getParameter<int>("binsDepth")),
      hist_drift_(conf.getParameter<int>("binsDrift")),
      trackerHitAssociatorConfig_(consumesCollector()) {
  //    anglefinder_=new  TrackLocalAngle(conf);
  hist_x_ = 50;
  hist_y_ = 100;
  min_x_ = -500.;
  max_x_ = 500.;
  min_y_ = -1500.;
  max_y_ = 500.;
  width_ = 0.0285;
  min_depth_ = -100.;
  max_depth_ = 400.;
  min_drift_ = -1000.;  //-200.;(conf.getParameter<double>("residualMax"))
  max_drift_ = 1000.;   //400.;

  t_trajTrack = consumes<TrajTrackAssociationCollection>(conf.getParameter<edm::InputTag>("src"));
  trackerTopoToken_ = esConsumes<TrackerTopology, TrackerTopologyRcd>();
  trackerGeomToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
}

~SiPixelLorentzAngle()

SiPixelLorentzAngle::~SiPixelLorentzAngle ( )

overridedefault

Member Function Documentation

analyze()

void SiPixelLorentzAngle::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 181 of file SiPixelLorentzAngle.cc.

References _h_drift_depth_adc2_, _h_drift_depth_adc_, _h_drift_depth_noadc_, Pixinfo::adc, Hit::alpha, edm::AssociationMap< Tag >::begin(), Hit::beta, bladeF_, Clust::charge, chi2_, Trajectory::chiSquared(), clust_, clustChargeMax_, SiPixelRecHit::cluster(), clustF_, clustSizeYMin_, Pixinfo::col, gather_cfg::cout, hcalRecHitTable_cff::depth, diskF_, shallow::drift(), PVValHelper::dx, PVValHelper::dy, MillePedeFileConverter_cfg::e, edm::AssociationMap< Tag >::empty(), edm::AssociationMap< Tag >::end(), eta_, event_, event_counter_, HcalObjRepresent::Fill(), fillPix(), Hit::gamma, edm::EventSetup::getHandle(), h_cluster_shape_adc_, h_cluster_shape_adc_rot_, h_cluster_shape_noadc_, h_cluster_shape_noadc_rot_, h_tracks_, hitCounter_, isflipped_, TrajectoryStateOnSurface::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, ladder_, layer_, TrajectoryStateOnSurface::localDirection(), BaseTrackerRecHit::localPosition(), TrajectoryStateOnSurface::localPosition(), muonTagProbeFilters_cff::matched, Clust::maxPixelCol, Clust::maxPixelRow, Trajectory::measurements(), Clust::minPixelCol, Clust::minPixelRow, module_, moduleF_, Trajectory::ndof(), ndof_, normChi2Max_, Pixinfo::npix, panelF_, PV3DBase< T, PVType, FrameType >::perp(), phi_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixelTracksCounter_, pixinfo_, pixinfoF_, edm::ESHandle< T >::product(), pt_, ptmin_, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), TrackerTopology::pxfBlade(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), TrackerTopology::pxfSide(), rpcPointValidation_cfi::recHit, rechit_, rechitF_, residualMax_, Pixinfo::row, run_, sideF_, simData_, simhit_, simhitF_, SiPixelLorentzAngleTree_, SiPixelLorentzAngleTreeForward_, Clust::size_x, Clust::size_y, PixelGeomDetUnit::specificTopology(), GeomDet::surface(), t_trajTrack, funct::tan(), Surface::toGlobal(), HLT_2024v14_cff::track, DetId::Tracker, PbPb_ZMuSkimMuonDPG_cff::tracker, trackerGeomToken_, trackerHitAssociatorConfig_, trackerTopoToken_, trackEventsCounter_, trackhit_, trackhitF_, usedHitCounter_, width_, Pixinfo::x, x, Hit::x, PV3DBase< T, PVType, FrameType >::x(), Clust::x, Rechit::x, Pixinfo::y, Hit::y, PV3DBase< T, PVType, FrameType >::y(), Clust::y, Rechit::y, and PV3DBase< T, PVType, FrameType >::z().

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

  event_counter_++;
  //    if(event_counter_ % 500 == 0) cout << "event number " << event_counter_ << endl;
  cout << "event number " << event_counter_ << endl;

  edm::ESHandle<TrackerGeometry> estracker = es.getHandle(trackerGeomToken_);
  const TrackerGeometry* tracker = &(*estracker);

  std::unique_ptr<TrackerHitAssociator> associate;
  if (simData_)
    associate = std::make_unique<TrackerHitAssociator>(e, trackerHitAssociatorConfig_);
  // restet values
  module_ = -1;
  layer_ = -1;
  ladder_ = -1;
  isflipped_ = -1;
  pt_ = -999;
  eta_ = 999;
  phi_ = 999;
  pixinfo_.npix = 0;

  run_ = e.id().run();
  event_ = e.id().event();

  // get the association map between tracks and trajectories
  edm::Handle<TrajTrackAssociationCollection> trajTrackCollectionHandle;
  e.getByToken(t_trajTrack, trajTrackCollectionHandle);
  if (!trajTrackCollectionHandle->empty()) {
    trackEventsCounter_++;
    for (TrajTrackAssociationCollection::const_iterator it = trajTrackCollectionHandle->begin();
         it != trajTrackCollectionHandle->end();
         ++it) {
      const Track& track = *it->val;
      const Trajectory& traj = *it->key;

      // get the trajectory measurements
      std::vector<TrajectoryMeasurement> tmColl = traj.measurements();
      //            TrajectoryStateOnSurface tsos = tsoscomb( itTraj->forwardPredictedState(), itTraj->backwardPredictedState() );
      pt_ = track.pt();
      eta_ = track.eta();
      phi_ = track.phi();
      chi2_ = traj.chiSquared();
      ndof_ = traj.ndof();
      if (pt_ < ptmin_)
        continue;
      // iterate over trajectory measurements
      std::vector<PSimHit> matched;
      h_tracks_->Fill(0);
      bool pixeltrack = false;
      for (std::vector<TrajectoryMeasurement>::const_iterator itTraj = tmColl.begin(); itTraj != tmColl.end();
           itTraj++) {
        if (!itTraj->updatedState().isValid())
          continue;
        TransientTrackingRecHit::ConstRecHitPointer recHit = itTraj->recHit();
        if (!recHit->isValid() || recHit->geographicalId().det() != DetId::Tracker)
          continue;
        unsigned int subDetID = (recHit->geographicalId().subdetId());
        if (subDetID == PixelSubdetector::PixelBarrel || subDetID == PixelSubdetector::PixelEndcap) {
          if (!pixeltrack) {
            h_tracks_->Fill(1);
            pixelTracksCounter_++;
          }
          pixeltrack = true;
        }
        if (subDetID == PixelSubdetector::PixelBarrel) {
          hitCounter_++;

          DetId detIdObj = recHit->geographicalId();
          const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(tracker->idToDet(detIdObj));
          if (!theGeomDet)
            continue;

          const PixelTopology* topol = &(theGeomDet->specificTopology());

          if (!topol)
            continue;

          layer_ = tTopo->pxbLayer(detIdObj);
          ladder_ = tTopo->pxbLadder(detIdObj);
          module_ = tTopo->pxbModule(detIdObj);
          float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 0.)).perp();
          float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 1.)).perp();
          if (tmp2 < tmp1)
            isflipped_ = 1;
          else
            isflipped_ = 0;
          const SiPixelRecHit* recHitPix = dynamic_cast<const SiPixelRecHit*>((*recHit).hit());
          if (!recHitPix)
            continue;
          rechit_.x = recHitPix->localPosition().x();
          rechit_.y = recHitPix->localPosition().y();
          SiPixelRecHit::ClusterRef const& cluster = recHitPix->cluster();

          // fill entries in clust_
          clust_.x = (cluster)->x();
          clust_.y = (cluster)->y();
          clust_.charge = (cluster->charge()) / 1000.;
          clust_.size_x = cluster->sizeX();
          clust_.size_y = cluster->sizeY();
          clust_.maxPixelCol = cluster->maxPixelCol();
          clust_.maxPixelRow = cluster->maxPixelRow();
          clust_.minPixelCol = cluster->minPixelCol();
          clust_.minPixelRow = cluster->minPixelRow();
          // fill entries in pixinfo_:
          fillPix(*cluster, topol, pixinfo_);
          // fill the trackhit info
          TrajectoryStateOnSurface tsos = itTraj->updatedState();
          if (!tsos.isValid()) {
            cout << "tsos not valid" << endl;
            continue;
          }
          LocalVector trackdirection = tsos.localDirection();
          LocalPoint trackposition = tsos.localPosition();

          if (trackdirection.z() == 0)
            continue;
          // the local position and direction
          trackhit_.alpha = atan2(trackdirection.z(), trackdirection.x());
          trackhit_.beta = atan2(trackdirection.z(), trackdirection.y());
          trackhit_.gamma = atan2(trackdirection.x(), trackdirection.y());
          trackhit_.x = trackposition.x();
          trackhit_.y = trackposition.y();

          // fill entries in simhit_:
          if (simData_) {
            matched.clear();
            matched = associate->associateHit((*recHitPix));
            float dr_start = 9999.;
            for (std::vector<PSimHit>::iterator isim = matched.begin(); isim != matched.end(); ++isim) {
              DetId simdetIdObj((*isim).detUnitId());
              if (simdetIdObj == detIdObj) {
                float sim_x1 = (*isim).entryPoint().x();  // width (row index, in col direction)
                float sim_y1 = (*isim).entryPoint().y();  // length (col index, in row direction)
                float sim_x2 = (*isim).exitPoint().x();
                float sim_y2 = (*isim).exitPoint().y();
                float sim_xpos = 0.5 * (sim_x1 + sim_x2);
                float sim_ypos = 0.5 * (sim_y1 + sim_y2);
                float sim_px = (*isim).momentumAtEntry().x();
                float sim_py = (*isim).momentumAtEntry().y();
                float sim_pz = (*isim).momentumAtEntry().z();

                float dr = (sim_xpos - (recHitPix->localPosition().x())) * (sim_xpos - recHitPix->localPosition().x()) +
                           (sim_ypos - recHitPix->localPosition().y()) * (sim_ypos - recHitPix->localPosition().y());
                if (dr < dr_start) {
                  simhit_.x = sim_xpos;
                  simhit_.y = sim_ypos;
                  simhit_.alpha = atan2(sim_pz, sim_px);
                  simhit_.beta = atan2(sim_pz, sim_py);
                  simhit_.gamma = atan2(sim_px, sim_py);
                  dr_start = dr;
                }
              }
            }  // end of filling simhit_
          }
          // is one pixel in cluster a large pixel ? (hit will be excluded)
          bool large_pix = false;
          for (int j = 0; j < pixinfo_.npix; j++) {
            int colpos = static_cast<int>(pixinfo_.col[j]);
            if (pixinfo_.row[j] == 0 || pixinfo_.row[j] == 79 || pixinfo_.row[j] == 80 || pixinfo_.row[j] == 159 ||
                colpos % 52 == 0 || colpos % 52 == 51) {
              large_pix = true;
            }
          }

          double residual = TMath::Sqrt((trackhit_.x - rechit_.x) * (trackhit_.x - rechit_.x) +
                                        (trackhit_.y - rechit_.y) * (trackhit_.y - rechit_.y));

          SiPixelLorentzAngleTree_->Fill();
          if (!large_pix && (chi2_ / ndof_) < normChi2Max_ && cluster->sizeY() >= clustSizeYMin_ &&
              residual < residualMax_ && (cluster->charge() < clustChargeMax_)) {
            usedHitCounter_++;
            // iterate over pixels in hit
            for (int j = 0; j < pixinfo_.npix; j++) {
              // use trackhits
              float dx = (pixinfo_.x[j] - (trackhit_.x - width_ / 2. / TMath::Tan(trackhit_.alpha))) * 10000.;
              float dy = (pixinfo_.y[j] - (trackhit_.y - width_ / 2. / TMath::Tan(trackhit_.beta))) * 10000.;
              float depth = dy * tan(trackhit_.beta);
              float drift = dx - dy * tan(trackhit_.gamma);
              _h_drift_depth_adc_[module_ + (layer_ - 1) * 8]->Fill(drift, depth, pixinfo_.adc[j]);
              _h_drift_depth_adc2_[module_ + (layer_ - 1) * 8]->Fill(drift, depth, pixinfo_.adc[j] * pixinfo_.adc[j]);
              _h_drift_depth_noadc_[module_ + (layer_ - 1) * 8]->Fill(drift, depth);
              if (layer_ == 3 && module_ == 1 && isflipped_) {
                float dx_rot = dx * TMath::Cos(trackhit_.gamma) + dy * TMath::Sin(trackhit_.gamma);
                float dy_rot = dy * TMath::Cos(trackhit_.gamma) - dx * TMath::Sin(trackhit_.gamma);
                h_cluster_shape_adc_->Fill(dx, dy, pixinfo_.adc[j]);
                h_cluster_shape_noadc_->Fill(dx, dy);
                h_cluster_shape_adc_rot_->Fill(dx_rot, dy_rot, pixinfo_.adc[j]);
                h_cluster_shape_noadc_rot_->Fill(dx_rot, dy_rot);
              }
            }
          }
        } else if (subDetID == PixelSubdetector::PixelEndcap) {
          DetId detIdObj = recHit->geographicalId();
          const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(tracker->idToDet(detIdObj));
          if (!theGeomDet)
            continue;

          const PixelTopology* topol = &(theGeomDet->specificTopology());

          if (!topol)
            continue;

          sideF_ = tTopo->pxfSide(detIdObj);
          diskF_ = tTopo->pxfDisk(detIdObj);
          bladeF_ = tTopo->pxfBlade(detIdObj);
          panelF_ = tTopo->pxfPanel(detIdObj);
          moduleF_ = tTopo->pxfModule(detIdObj);
          //float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
          //float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
          //if ( tmp2<tmp1 ) isflipped_ = 1;
          //else isflipped_ = 0;
          const SiPixelRecHit* recHitPix = dynamic_cast<const SiPixelRecHit*>((*recHit).hit());
          if (!recHitPix)
            continue;
          rechitF_.x = recHitPix->localPosition().x();
          rechitF_.y = recHitPix->localPosition().y();
          SiPixelRecHit::ClusterRef const& cluster = recHitPix->cluster();

          // fill entries in clust_
          clustF_.x = (cluster)->x();
          clustF_.y = (cluster)->y();
          clustF_.charge = (cluster->charge()) / 1000.;
          clustF_.size_x = cluster->sizeX();
          clustF_.size_y = cluster->sizeY();
          clustF_.maxPixelCol = cluster->maxPixelCol();
          clustF_.maxPixelRow = cluster->maxPixelRow();
          clustF_.minPixelCol = cluster->minPixelCol();
          clustF_.minPixelRow = cluster->minPixelRow();
          // fill entries in pixinfo_:
          fillPix(*cluster, topol, pixinfoF_);
          // fill the trackhit info
          TrajectoryStateOnSurface tsos = itTraj->updatedState();
          if (!tsos.isValid()) {
            cout << "tsos not valid" << endl;
            continue;
          }
          LocalVector trackdirection = tsos.localDirection();
          LocalPoint trackposition = tsos.localPosition();

          if (trackdirection.z() == 0)
            continue;
          // the local position and direction
          trackhitF_.alpha = atan2(trackdirection.z(), trackdirection.x());
          trackhitF_.beta = atan2(trackdirection.z(), trackdirection.y());
          trackhitF_.gamma = atan2(trackdirection.x(), trackdirection.y());
          trackhitF_.x = trackposition.x();
          trackhitF_.y = trackposition.y();

          // fill entries in simhit_:
          if (simData_) {
            matched.clear();
            matched = associate->associateHit((*recHitPix));
            float dr_start = 9999.;
            for (std::vector<PSimHit>::iterator isim = matched.begin(); isim != matched.end(); ++isim) {
              DetId simdetIdObj((*isim).detUnitId());
              if (simdetIdObj == detIdObj) {
                float sim_x1 = (*isim).entryPoint().x();  // width (row index, in col direction)
                float sim_y1 = (*isim).entryPoint().y();  // length (col index, in row direction)
                float sim_x2 = (*isim).exitPoint().x();
                float sim_y2 = (*isim).exitPoint().y();
                float sim_xpos = 0.5 * (sim_x1 + sim_x2);
                float sim_ypos = 0.5 * (sim_y1 + sim_y2);
                float sim_px = (*isim).momentumAtEntry().x();
                float sim_py = (*isim).momentumAtEntry().y();
                float sim_pz = (*isim).momentumAtEntry().z();

                float dr = (sim_xpos - (recHitPix->localPosition().x())) * (sim_xpos - recHitPix->localPosition().x()) +
                           (sim_ypos - recHitPix->localPosition().y()) * (sim_ypos - recHitPix->localPosition().y());
                if (dr < dr_start) {
                  simhitF_.x = sim_xpos;
                  simhitF_.y = sim_ypos;
                  simhitF_.alpha = atan2(sim_pz, sim_px);
                  simhitF_.beta = atan2(sim_pz, sim_py);
                  simhitF_.gamma = atan2(sim_px, sim_py);
                  dr_start = dr;
                }
              }
            }  // end of filling simhit_
          }
          SiPixelLorentzAngleTreeForward_->Fill();
        }
      }  //end iteration over trajectory measurements
    }    //end iteration over trajectories
  }
}

beginJob()

void SiPixelLorentzAngle::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 71 of file SiPixelLorentzAngle.cc.

References _h_drift_depth_, _h_drift_depth_adc2_, _h_drift_depth_adc_, _h_drift_depth_noadc_, _h_mean_, Pixinfo::adc, bladeF_, chi2_, clust_, clustF_, Pixinfo::col, diskF_, eta_, event_, event_counter_, filename_, h_cluster_shape_, h_cluster_shape_adc_, h_cluster_shape_adc_rot_, h_cluster_shape_noadc_, h_cluster_shape_noadc_rot_, h_cluster_shape_rot_, h_tracks_, hFile_, hist_depth_, hist_drift_, hist_x_, hist_y_, hitCounter_, isflipped_, ladder_, layer_, max_depth_, max_drift_, max_x_, max_y_, min_depth_, min_drift_, min_x_, min_y_, module_, moduleF_, Skims_PA_cff::name, ndof_, Pixinfo::npix, panelF_, phi_, pixelTracksCounter_, pixinfo_, pixinfoF_, pt_, rechit_, rechitF_, Pixinfo::row, run_, sideF_, simhit_, simhitF_, SiPixelLorentzAngleTree_, SiPixelLorentzAngleTreeForward_, trackEventsCounter_, trackhit_, trackhitF_, usedHitCounter_, Pixinfo::x, and Pixinfo::y.

                                   {
  //    cout << "started SiPixelLorentzAngle" << endl;
  hFile_ = new TFile(filename_.c_str(), "RECREATE");
  int bufsize = 64000;
  // create tree structure
  SiPixelLorentzAngleTree_ = new TTree("SiPixelLorentzAngleTree_", "SiPixel LorentzAngle tree", bufsize);
  SiPixelLorentzAngleTree_->Branch("run", &run_, "run/I", bufsize);
  SiPixelLorentzAngleTree_->Branch("event", &event_, "event/I", bufsize);
  SiPixelLorentzAngleTree_->Branch("module", &module_, "module/I", bufsize);
  SiPixelLorentzAngleTree_->Branch("ladder", &ladder_, "ladder/I", bufsize);
  SiPixelLorentzAngleTree_->Branch("layer", &layer_, "layer/I", bufsize);
  SiPixelLorentzAngleTree_->Branch("isflipped", &isflipped_, "isflipped/I", bufsize);
  SiPixelLorentzAngleTree_->Branch("pt", &pt_, "pt/F", bufsize);
  SiPixelLorentzAngleTree_->Branch("eta", &eta_, "eta/F", bufsize);
  SiPixelLorentzAngleTree_->Branch("phi", &phi_, "phi/F", bufsize);
  SiPixelLorentzAngleTree_->Branch("chi2", &chi2_, "chi2/D", bufsize);
  SiPixelLorentzAngleTree_->Branch("ndof", &ndof_, "ndof/D", bufsize);
  SiPixelLorentzAngleTree_->Branch("trackhit", &trackhit_, "x/F:y/F:alpha/D:beta/D:gamma_/D", bufsize);
  SiPixelLorentzAngleTree_->Branch("simhit", &simhit_, "x/F:y/F:alpha/D:beta/D:gamma_/D", bufsize);
  SiPixelLorentzAngleTree_->Branch("npix", &pixinfo_.npix, "npix/I", bufsize);
  SiPixelLorentzAngleTree_->Branch("rowpix", pixinfo_.row, "row[npix]/F", bufsize);
  SiPixelLorentzAngleTree_->Branch("colpix", pixinfo_.col, "col[npix]/F", bufsize);
  SiPixelLorentzAngleTree_->Branch("adc", pixinfo_.adc, "adc[npix]/F", bufsize);
  SiPixelLorentzAngleTree_->Branch("xpix", pixinfo_.x, "x[npix]/F", bufsize);
  SiPixelLorentzAngleTree_->Branch("ypix", pixinfo_.y, "y[npix]/F", bufsize);
  SiPixelLorentzAngleTree_->Branch(
      "clust",
      &clust_,
      "x/F:y/F:charge/F:size_x/I:size_y/I:maxPixelCol/I:maxPixelRow:minPixelCol/I:minPixelRow/I",
      bufsize);
  SiPixelLorentzAngleTree_->Branch("rechit", &rechit_, "x/F:y/F", bufsize);

  SiPixelLorentzAngleTreeForward_ =
      new TTree("SiPixelLorentzAngleTreeForward_", "SiPixel LorentzAngle tree forward", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("run", &run_, "run/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("event", &event_, "event/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("side", &sideF_, "side/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("disk", &diskF_, "disk/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("blade", &bladeF_, "blade/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("panel", &panelF_, "panel/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("module", &moduleF_, "module/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("pt", &pt_, "pt/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("eta", &eta_, "eta/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("phi", &phi_, "phi/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("chi2", &chi2_, "chi2/D", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("ndof", &ndof_, "ndof/D", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("trackhit", &trackhitF_, "x/F:y/F:alpha/D:beta/D:gamma_/D", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("simhit", &simhitF_, "x/F:y/F:alpha/D:beta/D:gamma_/D", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("npix", &pixinfoF_.npix, "npix/I", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("rowpix", pixinfoF_.row, "row[npix]/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("colpix", pixinfoF_.col, "col[npix]/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("adc", pixinfoF_.adc, "adc[npix]/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("xpix", pixinfoF_.x, "x[npix]/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("ypix", pixinfoF_.y, "y[npix]/F", bufsize);
  SiPixelLorentzAngleTreeForward_->Branch(
      "clust",
      &clustF_,
      "x/F:y/F:charge/F:size_x/I:size_y/I:maxPixelCol/I:maxPixelRow:minPixelCol/I:minPixelRow/I",
      bufsize);
  SiPixelLorentzAngleTreeForward_->Branch("rechit", &rechitF_, "x/F:y/F", bufsize);

  //book histograms
  char name[128];
  for (int i_module = 1; i_module <= 8; i_module++) {
    for (int i_layer = 1; i_layer <= 3; i_layer++) {
      sprintf(name, "h_drift_depth_adc_layer%i_module%i", i_layer, i_module);
      _h_drift_depth_adc_[i_module + (i_layer - 1) * 8] =
          new TH2F(name, name, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
      sprintf(name, "h_drift_depth_adc2_layer%i_module%i", i_layer, i_module);
      _h_drift_depth_adc2_[i_module + (i_layer - 1) * 8] =
          new TH2F(name, name, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
      sprintf(name, "h_drift_depth_noadc_layer%i_module%i", i_layer, i_module);
      _h_drift_depth_noadc_[i_module + (i_layer - 1) * 8] =
          new TH2F(name, name, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
      sprintf(name, "h_drift_depth_layer%i_module%i", i_layer, i_module);
      _h_drift_depth_[i_module + (i_layer - 1) * 8] =
          new TH2F(name, name, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
      sprintf(name, "h_mean_layer%i_module%i", i_layer, i_module);
      _h_mean_[i_module + (i_layer - 1) * 8] = new TH1F(name, name, hist_depth_, min_depth_, max_depth_);
    }
  }

  // just for some expaining plots
  h_cluster_shape_adc_ = new TH2F(
      "h_cluster_shape_adc", "cluster shape with adc weight", hist_x_, min_x_, max_x_, hist_y_, min_y_, max_y_);
  h_cluster_shape_noadc_ = new TH2F(
      "h_cluster_shape_noadc", "cluster shape without adc weight", hist_x_, min_x_, max_x_, hist_y_, min_y_, max_y_);
  h_cluster_shape_ = new TH2F("h_cluster_shape", "cluster shape", hist_x_, min_x_, max_x_, hist_y_, min_y_, max_y_);
  h_cluster_shape_adc_rot_ = new TH2F(
      "h_cluster_shape_adc_rot", "cluster shape with adc weight", hist_x_, min_x_, max_x_, hist_y_, -max_y_, -min_y_);
  h_cluster_shape_noadc_rot_ = new TH2F("h_cluster_shape_noadc_rot",
                                        "cluster shape without adc weight",
                                        hist_x_,
                                        min_x_,
                                        max_x_,
                                        hist_y_,
                                        -max_y_,
                                        -min_y_);
  h_cluster_shape_rot_ =
      new TH2F("h_cluster_shape_rot", "cluster shape", hist_x_, min_x_, max_x_, hist_y_, -max_y_, -min_y_);
  h_tracks_ = new TH1F("h_tracks", "h_tracks", 2, 0., 2.);
  event_counter_ = 0;
  trackEventsCounter_ = 0;
  //    trackcounter_ = 0;
  hitCounter_ = 0;
  usedHitCounter_ = 0;
  pixelTracksCounter_ = 0;
}

endJob()

void SiPixelLorentzAngle::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 471 of file SiPixelLorentzAngle.cc.

References _h_drift_depth_, _h_drift_depth_adc2_, _h_drift_depth_adc_, _h_drift_depth_noadc_, _h_mean_, nano_mu_local_reco_cff::chi2, gather_cfg::cout, StorageManager_cfg::e1, event_counter_, DeadROC_duringRun::f1, filenameFit_, findMean(), h_cluster_shape_adc_, h_cluster_shape_adc_rot_, h_cluster_shape_noadc_, h_cluster_shape_noadc_rot_, h_drift_depth_adc_slice_, h_tracks_, hFile_, hist_depth_, hist_drift_, hitCounter_, mps_fire::i, max_drift_, min_drift_, LaserDQM_cfg::p1, pixelTracksCounter_, TtFullHadEvtBuilder_cfi::prob, trackEventsCounter_, pileupReCalc_HLTpaths::trunc, and usedHitCounter_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                 {
  // produce histograms with the average adc counts
  for (int i_ring = 1; i_ring <= 24; i_ring++) {
    _h_drift_depth_[i_ring]->Divide(_h_drift_depth_adc_[i_ring], _h_drift_depth_noadc_[i_ring]);
  }

  h_drift_depth_adc_slice_ =
      new TH1F("h_drift_depth_adc_slice", "slice of adc histogram", hist_drift_, min_drift_, max_drift_);

  TF1* f1 = new TF1("f1", "[0] + [1]*x", 50., 235.);
  f1->SetParName(0, "p0");
  f1->SetParName(1, "p1");
  f1->SetParameter(0, 0);
  f1->SetParameter(1, 0.4);
  ofstream fLorentzFit(filenameFit_.c_str(), ios::trunc);
  cout.precision(4);
  fLorentzFit << "module"
              << "\t"
              << "layer"
              << "\t"
              << "offset"
              << "\t"
              << "error"
              << "\t"
              << "slope"
              << "\t"
              << "error"
              << "\t"
                 "rel.err"
              << "\t"
                 "pull"
              << "\t"
              << "chi2"
              << "\t"
              << "prob" << endl;
  //loop over modlues and layers to fit the lorentz angle
  for (int i_layer = 1; i_layer <= 3; i_layer++) {
    for (int i_module = 1; i_module <= 8; i_module++) {
      //loop over bins in depth (z-local-coordinate) (in order to fit slices)
      for (int i = 1; i <= hist_depth_; i++) {
        findMean(i, (i_module + (i_layer - 1) * 8));
      }  // end loop over bins in depth
      _h_mean_[i_module + (i_layer - 1) * 8]->Fit(f1, "ERQ");
      double p0 = f1->GetParameter(0);
      double e0 = f1->GetParError(0);
      double p1 = f1->GetParameter(1);
      double e1 = f1->GetParError(1);
      double chi2 = f1->GetChisquare();
      double prob = f1->GetProb();
      fLorentzFit << setprecision(4) << i_module << "\t" << i_layer << "\t" << p0 << "\t" << e0 << "\t" << p1
                  << setprecision(3) << "\t" << e1 << "\t" << e1 / p1 * 100. << "\t" << (p1 - 0.424) / e1 << "\t"
                  << chi2 << "\t" << prob << endl;
    }
  }  // end loop over modules and layers
  fLorentzFit.close();
  hFile_->cd();
  for (int i_module = 1; i_module <= 8; i_module++) {
    for (int i_layer = 1; i_layer <= 3; i_layer++) {
      _h_drift_depth_adc_[i_module + (i_layer - 1) * 8]->Write();
      _h_drift_depth_adc2_[i_module + (i_layer - 1) * 8]->Write();
      _h_drift_depth_noadc_[i_module + (i_layer - 1) * 8]->Write();
      _h_drift_depth_[i_module + (i_layer - 1) * 8]->Write();
      _h_mean_[i_module + (i_layer - 1) * 8]->Write();
    }
  }
  h_cluster_shape_adc_->Write();
  h_cluster_shape_noadc_->Write();
  h_cluster_shape_adc_rot_->Write();
  h_cluster_shape_noadc_rot_->Write();
  h_tracks_->Write();

  hFile_->Write();
  hFile_->Close();
  cout << "events: " << event_counter_ << endl;
  cout << "events with tracks: " << trackEventsCounter_ << endl;
  cout << "events with pixeltracks: " << pixelTracksCounter_ << endl;
  cout << "hits in the pixel: " << hitCounter_ << endl;
  cout << "number of used Hits: " << usedHitCounter_ << endl;
}

fillPix()

void SiPixelLorentzAngle::fillPix ( const SiPixelCluster &  LocPix, const PixelTopology *  topol, Pixinfo &  pixinfo  )

inlineprivate

Definition at line 551 of file SiPixelLorentzAngle.cc.

References Pixinfo::adc, Pixinfo::col, Topology::localPosition(), Pixinfo::npix, SiPixelCluster::pixels(), Pixinfo::row, Pixinfo::x, PV3DBase< T, PVType, FrameType >::x(), Pixinfo::y, and PV3DBase< T, PVType, FrameType >::y().

Referenced by analyze().

{
  const std::vector<SiPixelCluster::Pixel>& pixvector = LocPix.pixels();
  pixinfo.npix = 0;
  for (std::vector<SiPixelCluster::Pixel>::const_iterator itPix = pixvector.begin(); itPix != pixvector.end();
       itPix++) {
    //  for(pixinfo.npix = 0; pixinfo.npix < static_cast<int>(pixvector.size()); ++pixinfo.npix) {
    pixinfo.row[pixinfo.npix] = itPix->x;
    pixinfo.col[pixinfo.npix] = itPix->y;
    pixinfo.adc[pixinfo.npix] = itPix->adc;
    LocalPoint lp = topol->localPosition(MeasurementPoint(itPix->x + 0.5, itPix->y + 0.5));
    pixinfo.x[pixinfo.npix] = lp.x();
    pixinfo.y[pixinfo.npix] = lp.y();
    pixinfo.npix++;
  }
}

findMean()

void SiPixelLorentzAngle::findMean ( int  i, int  i_ring  )

private

Definition at line 569 of file SiPixelLorentzAngle.cc.

References _h_drift_depth_, _h_drift_depth_adc2_, _h_drift_depth_adc_, _h_drift_depth_noadc_, _h_mean_, relativeConstraints::error, h_drift_depth_adc_slice_, hist_drift_, mps_fire::i, dqmiolumiharvest::j, SiStripPI::mean, and mathSSE::sqrt().

Referenced by endJob().

                                                    {
  double nentries = 0;

  h_drift_depth_adc_slice_->Reset("ICE");

  // determine sigma and sigma^2 of the adc counts and average adc counts
  //loop over bins in drift width
  for (int j = 1; j <= hist_drift_; j++) {
    if (_h_drift_depth_noadc_[i_ring]->GetBinContent(j, i) >= 1) {
      double adc_error2 =
          (_h_drift_depth_adc2_[i_ring]->GetBinContent(j, i) - _h_drift_depth_adc_[i_ring]->GetBinContent(j, i) *
                                                                   _h_drift_depth_adc_[i_ring]->GetBinContent(j, i) /
                                                                   _h_drift_depth_noadc_[i_ring]->GetBinContent(j, i)) /
          _h_drift_depth_noadc_[i_ring]->GetBinContent(j, i);
      _h_drift_depth_adc_[i_ring]->SetBinError(j, i, sqrt(adc_error2));
      double error2 = adc_error2 / (_h_drift_depth_noadc_[i_ring]->GetBinContent(j, i) - 1.);
      _h_drift_depth_[i_ring]->SetBinError(j, i, sqrt(error2));
    } else {
      _h_drift_depth_[i_ring]->SetBinError(j, i, 0);
      _h_drift_depth_adc_[i_ring]->SetBinError(j, i, 0);
    }
    h_drift_depth_adc_slice_->SetBinContent(j, _h_drift_depth_adc_[i_ring]->GetBinContent(j, i));
    h_drift_depth_adc_slice_->SetBinError(j, _h_drift_depth_adc_[i_ring]->GetBinError(j, i));
    nentries += _h_drift_depth_noadc_[i_ring]->GetBinContent(j, i);
  }  // end loop over bins in drift width

  double mean = h_drift_depth_adc_slice_->GetMean(1);
  double error = 0;
  if (nentries != 0) {
    error = h_drift_depth_adc_slice_->GetRMS(1) / sqrt(nentries);
  }

  _h_mean_[i_ring]->SetBinContent(i, mean);
  _h_mean_[i_ring]->SetBinError(i, error);
}

Member Data Documentation

_h_drift_depth_

std::map<int, TH2F *> analyzer::SiPixelLorentzAngle::_h_drift_depth_

private

Definition at line 167 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), endJob(), and findMean().

_h_drift_depth_adc2_

std::map<int, TH2F *> analyzer::SiPixelLorentzAngle::_h_drift_depth_adc2_

private

Definition at line 165 of file SiPixelLorentzAngle.h.

Referenced by analyze(), beginJob(), endJob(), and findMean().

_h_drift_depth_adc_

std::map<int, TH2F *> analyzer::SiPixelLorentzAngle::_h_drift_depth_adc_

private

Definition at line 164 of file SiPixelLorentzAngle.h.

Referenced by analyze(), beginJob(), endJob(), and findMean().

_h_drift_depth_noadc_

std::map<int, TH2F *> analyzer::SiPixelLorentzAngle::_h_drift_depth_noadc_

private

Definition at line 166 of file SiPixelLorentzAngle.h.

Referenced by analyze(), beginJob(), endJob(), and findMean().

_h_mean_

std::map<int, TH1F *> analyzer::SiPixelLorentzAngle::_h_mean_

private

Definition at line 169 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), endJob(), and findMean().

bladeF_

int analyzer::SiPixelLorentzAngle::bladeF_

private

Definition at line 123 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

chi2_

double analyzer::SiPixelLorentzAngle::chi2_

private

Definition at line 111 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

chi2F_

double analyzer::SiPixelLorentzAngle::chi2F_

private

Definition at line 129 of file SiPixelLorentzAngle.h.

clust_

Clust analyzer::SiPixelLorentzAngle::clust_

private

Definition at line 115 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

clustChargeMax_

double analyzer::SiPixelLorentzAngle::clustChargeMax_

private

Definition at line 144 of file SiPixelLorentzAngle.h.

Referenced by analyze().

clustF_

Clust analyzer::SiPixelLorentzAngle::clustF_

private

Definition at line 133 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

clustSizeYMin_

int analyzer::SiPixelLorentzAngle::clustSizeYMin_

private

Definition at line 142 of file SiPixelLorentzAngle.h.

Referenced by analyze().

diskF_

int analyzer::SiPixelLorentzAngle::diskF_

private

Definition at line 122 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

eta_

float analyzer::SiPixelLorentzAngle::eta_

private

Definition at line 109 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

etaF_

float analyzer::SiPixelLorentzAngle::etaF_

private

Definition at line 127 of file SiPixelLorentzAngle.h.

event_

int analyzer::SiPixelLorentzAngle::event_

private

Definition at line 103 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

event_counter_

int analyzer::SiPixelLorentzAngle::event_counter_

private

Definition at line 178 of file SiPixelLorentzAngle.h.

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

eventF_

int analyzer::SiPixelLorentzAngle::eventF_

private

Definition at line 120 of file SiPixelLorentzAngle.h.

filename_

std::string analyzer::SiPixelLorentzAngle::filename_

private

Definition at line 137 of file SiPixelLorentzAngle.h.

Referenced by beginJob().

filenameFit_

std::string analyzer::SiPixelLorentzAngle::filenameFit_

private

Definition at line 138 of file SiPixelLorentzAngle.h.

Referenced by endJob().

h_cluster_shape_

TH2F* analyzer::SiPixelLorentzAngle::h_cluster_shape_

private

Definition at line 172 of file SiPixelLorentzAngle.h.

Referenced by beginJob().

h_cluster_shape_adc_

TH2F* analyzer::SiPixelLorentzAngle::h_cluster_shape_adc_

private

Definition at line 170 of file SiPixelLorentzAngle.h.

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

h_cluster_shape_adc_rot_

TH2F* analyzer::SiPixelLorentzAngle::h_cluster_shape_adc_rot_

private

Definition at line 173 of file SiPixelLorentzAngle.h.

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

h_cluster_shape_noadc_

TH2F* analyzer::SiPixelLorentzAngle::h_cluster_shape_noadc_

private

Definition at line 171 of file SiPixelLorentzAngle.h.

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

h_cluster_shape_noadc_rot_

TH2F* analyzer::SiPixelLorentzAngle::h_cluster_shape_noadc_rot_

private

Definition at line 174 of file SiPixelLorentzAngle.h.

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

h_cluster_shape_rot_

TH2F* analyzer::SiPixelLorentzAngle::h_cluster_shape_rot_

private

Definition at line 175 of file SiPixelLorentzAngle.h.

Referenced by beginJob().

h_drift_depth_adc_slice_

TH1F* analyzer::SiPixelLorentzAngle::h_drift_depth_adc_slice_

private

Definition at line 168 of file SiPixelLorentzAngle.h.

Referenced by endJob(), and findMean().

h_tracks_

TH1F* analyzer::SiPixelLorentzAngle::h_tracks_

private

Definition at line 176 of file SiPixelLorentzAngle.h.

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

hFile_

TFile* analyzer::SiPixelLorentzAngle::hFile_

private

Definition at line 97 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and endJob().

hist_depth_

int analyzer::SiPixelLorentzAngle::hist_depth_

private

Definition at line 145 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and endJob().

hist_drift_

int analyzer::SiPixelLorentzAngle::hist_drift_

private

Definition at line 146 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), endJob(), and findMean().

hist_x_

int analyzer::SiPixelLorentzAngle::hist_x_

private

Definition at line 152 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

hist_y_

int analyzer::SiPixelLorentzAngle::hist_y_

private

Definition at line 153 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

hitCounter_

int analyzer::SiPixelLorentzAngle::hitCounter_

private

Definition at line 178 of file SiPixelLorentzAngle.h.

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

isflipped_

int analyzer::SiPixelLorentzAngle::isflipped_

private

Definition at line 107 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

ladder_

int analyzer::SiPixelLorentzAngle::ladder_

private

Definition at line 105 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

layer_

int analyzer::SiPixelLorentzAngle::layer_

private

Definition at line 106 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

max_depth_

double analyzer::SiPixelLorentzAngle::max_depth_

private

Definition at line 160 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

max_drift_

double analyzer::SiPixelLorentzAngle::max_drift_

private

Definition at line 162 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), endJob(), and SiPixelLorentzAngle().

max_x_

double analyzer::SiPixelLorentzAngle::max_x_

private

Definition at line 155 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

max_y_

double analyzer::SiPixelLorentzAngle::max_y_

private

Definition at line 157 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

min_depth_

double analyzer::SiPixelLorentzAngle::min_depth_

private

Definition at line 159 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

min_drift_

double analyzer::SiPixelLorentzAngle::min_drift_

private

Definition at line 161 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), endJob(), and SiPixelLorentzAngle().

min_x_

double analyzer::SiPixelLorentzAngle::min_x_

private

Definition at line 154 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

min_y_

double analyzer::SiPixelLorentzAngle::min_y_

private

Definition at line 156 of file SiPixelLorentzAngle.h.

Referenced by beginJob(), and SiPixelLorentzAngle().

module_

int analyzer::SiPixelLorentzAngle::module_

private

Definition at line 104 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

moduleF_

int analyzer::SiPixelLorentzAngle::moduleF_

private

Definition at line 125 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

ndof_

double analyzer::SiPixelLorentzAngle::ndof_

private

Definition at line 112 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

ndofF_

double analyzer::SiPixelLorentzAngle::ndofF_

private

Definition at line 130 of file SiPixelLorentzAngle.h.

normChi2Max_

double analyzer::SiPixelLorentzAngle::normChi2Max_

private

Definition at line 141 of file SiPixelLorentzAngle.h.

Referenced by analyze().

panelF_

int analyzer::SiPixelLorentzAngle::panelF_

private

Definition at line 124 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

phi_

float analyzer::SiPixelLorentzAngle::phi_

private

Definition at line 110 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

phiF_

float analyzer::SiPixelLorentzAngle::phiF_

private

Definition at line 128 of file SiPixelLorentzAngle.h.

pixelTracksCounter_

int analyzer::SiPixelLorentzAngle::pixelTracksCounter_

private

Definition at line 178 of file SiPixelLorentzAngle.h.

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

pixinfo_

Pixinfo analyzer::SiPixelLorentzAngle::pixinfo_

private

Definition at line 113 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

pixinfoF_

Pixinfo analyzer::SiPixelLorentzAngle::pixinfoF_

private

Definition at line 131 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

pt_

float analyzer::SiPixelLorentzAngle::pt_

private

Definition at line 108 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

ptF_

float analyzer::SiPixelLorentzAngle::ptF_

private

Definition at line 126 of file SiPixelLorentzAngle.h.

ptmin_

double analyzer::SiPixelLorentzAngle::ptmin_

private

Definition at line 139 of file SiPixelLorentzAngle.h.

Referenced by analyze().

rechit_

Rechit analyzer::SiPixelLorentzAngle::rechit_

private

Definition at line 116 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

rechitF_

Rechit analyzer::SiPixelLorentzAngle::rechitF_

private

Definition at line 134 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

residualMax_

double analyzer::SiPixelLorentzAngle::residualMax_

private

Definition at line 143 of file SiPixelLorentzAngle.h.

Referenced by analyze().

RHBuilder

const TransientTrackingRecHitBuilder* analyzer::SiPixelLorentzAngle::RHBuilder

private

Definition at line 185 of file SiPixelLorentzAngle.h.

run_

int analyzer::SiPixelLorentzAngle::run_

private

Definition at line 102 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

runF_

int analyzer::SiPixelLorentzAngle::runF_

private

Definition at line 119 of file SiPixelLorentzAngle.h.

sideF_

int analyzer::SiPixelLorentzAngle::sideF_

private

Definition at line 121 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

simData_

bool analyzer::SiPixelLorentzAngle::simData_

private

Definition at line 140 of file SiPixelLorentzAngle.h.

Referenced by analyze().

simhit_

Hit analyzer::SiPixelLorentzAngle::simhit_

private

Definition at line 114 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

simhitF_

Hit analyzer::SiPixelLorentzAngle::simhitF_

private

Definition at line 132 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

SiPixelLorentzAngleTree_

TTree* analyzer::SiPixelLorentzAngle::SiPixelLorentzAngleTree_

private

Definition at line 98 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

SiPixelLorentzAngleTreeForward_

TTree* analyzer::SiPixelLorentzAngle::SiPixelLorentzAngleTreeForward_

private

Definition at line 99 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

t_trajTrack

edm::EDGetTokenT<TrajTrackAssociationCollection> analyzer::SiPixelLorentzAngle::t_trajTrack

private

Definition at line 189 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and SiPixelLorentzAngle().

theEstimator

Chi2MeasurementEstimator* analyzer::SiPixelLorentzAngle::theEstimator

private

Definition at line 184 of file SiPixelLorentzAngle.h.

theFitter

const KFTrajectoryFitter* analyzer::SiPixelLorentzAngle::theFitter

private

Definition at line 187 of file SiPixelLorentzAngle.h.

thePropagator

PropagatorWithMaterial* analyzer::SiPixelLorentzAngle::thePropagator

private

Definition at line 181 of file SiPixelLorentzAngle.h.

thePropagatorOp

PropagatorWithMaterial* analyzer::SiPixelLorentzAngle::thePropagatorOp

private

Definition at line 182 of file SiPixelLorentzAngle.h.

theSmoother

const KFTrajectorySmoother* analyzer::SiPixelLorentzAngle::theSmoother

private

Definition at line 186 of file SiPixelLorentzAngle.h.

theUpdator

KFUpdator* analyzer::SiPixelLorentzAngle::theUpdator

private

Definition at line 183 of file SiPixelLorentzAngle.h.

trackerGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> analyzer::SiPixelLorentzAngle::trackerGeomToken_

private

Definition at line 191 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and SiPixelLorentzAngle().

trackerHitAssociatorConfig_

TrackerHitAssociator::Config analyzer::SiPixelLorentzAngle::trackerHitAssociatorConfig_

private

Definition at line 149 of file SiPixelLorentzAngle.h.

Referenced by analyze().

trackerTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> analyzer::SiPixelLorentzAngle::trackerTopoToken_

private

Definition at line 190 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and SiPixelLorentzAngle().

trackEventsCounter_

int analyzer::SiPixelLorentzAngle::trackEventsCounter_

private

Definition at line 178 of file SiPixelLorentzAngle.h.

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

trackhit_

Hit analyzer::SiPixelLorentzAngle::trackhit_

private

Definition at line 114 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

trackhitF_

Hit analyzer::SiPixelLorentzAngle::trackhitF_

private

Definition at line 132 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and beginJob().

tsTransform

TrajectoryStateTransform analyzer::SiPixelLorentzAngle::tsTransform

private

Definition at line 188 of file SiPixelLorentzAngle.h.

usedHitCounter_

int analyzer::SiPixelLorentzAngle::usedHitCounter_

private

Definition at line 178 of file SiPixelLorentzAngle.h.

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

width_

double analyzer::SiPixelLorentzAngle::width_

private

Definition at line 158 of file SiPixelLorentzAngle.h.

Referenced by analyze(), and SiPixelLorentzAngle().