CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
SiPixelLorentzAnglePCLWorker Class Reference

#include <CalibTracker/SiPixelLorentzAnglePCLWorker/src/SiPixelLorentzAnglePCLWorker.cc>

Inheritance diagram for SiPixelLorentzAnglePCLWorker:
DQMEDAnalyzer edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

Public Member Functions

 SiPixelLorentzAnglePCLWorker (const edm::ParameterSet &)
 
 ~SiPixelLorentzAnglePCLWorker () override=default
 
- Public Member Functions inherited from DQMEDAnalyzer
void accumulate (edm::Event const &event, edm::EventSetup const &setup) final
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
void beginStream (edm::StreamID id) final
 
 DQMEDAnalyzer ()
 
void endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual bool getCanSaveByLumi ()
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static void globalEndJob (DQMEDAnalyzerGlobalCache const *)
 
static void globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
 
static void globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
 
static std::unique_ptr< DQMEDAnalyzerGlobalCache > initializeGlobalCache (edm::ParameterSet const &)
 

Private Member Functions

void analyze (edm::Event const &, edm::EventSetup const &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
LorentzAngleAnalysisTypeEnum convertStringToLorentzAngleAnalysisTypeEnum (std::string type)
 
void dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
 
void dqmEndRun (edm::Run const &, edm::EventSetup const &)
 
const Pixinfo fillPix (const SiPixelCluster &LocPix, const PixelTopology *topol) const
 
const std::pair< LocalPoint, LocalPoint > surface_deformation (const PixelTopology *topol, TrajectoryStateOnSurface &tsos, const SiPixelRecHit *recHitPix) const
 

Private Attributes

LorentzAngleAnalysisTypeEnum analysisType_
 
int bladeF_
 
int bx_
 
double chi2_
 
Clust clust_
 
double clustChargeMaxPerLength_
 
Clust clustF_
 
int clustSizeXMax_
 
std::vector< int > clustSizeYMin_
 
int diskF_
 
float eta_
 
long int event_
 
std::string filename_
 
std::string folder_
 
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomEsToken_
 
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomPerEventEsToken_
 
std::unique_ptr< TFile > hFile_
 
int hist_depth_
 
int hist_drift_
 
SiPixelLorentzAngleCalibrationHistograms iHists
 
int isflipped_
 
int ladder_
 
int layer_
 
int lumiblock_
 
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magneticFieldToken_
 
int module_
 
int moduleF_
 
double ndof_
 
std::vector< std::string > newmodulelist_
 
double normChi2Max_
 
bool notInPCL_
 
int orbit_
 
int panelF_
 
float phi_
 
Pixinfo pixinfo_
 
Pixinfo pixinfoF_
 
float pt_
 
double ptmin_
 
float qScale_
 
float qScaleF_
 
Rechit rechit_
 
Rechit rechitCorr_
 
Rechit rechitCorrF_
 
Rechit rechitF_
 
double residualMax_
 
float rQmQt_
 
float rQmQtF_
 
int run_
 
int sideF_
 
Hit simhit_
 
Hit simhitF_
 
std::unique_ptr< TTree > SiPixelLorentzAngleTreeBarrel_
 
std::unique_ptr< TTree > SiPixelLorentzAngleTreeForward_
 
edm::ESGetToken< SiPixelTemplateDBObject, SiPixelTemplateDBObjectESProducerRcd > siPixelTemplateEsToken_
 
edm::EDGetTokenT< TrajTrackAssociationCollection > t_trajTrack
 
const SiPixelTemplateDBObject * templateDBobject_
 
std::vector< SiPixelTemplateStore > thePixelTemp_
 
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoEsToken_
 
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoPerEventEsToken_
 
Hit trackhit_
 
float trackhitCorrX_
 
float trackhitCorrXF_
 
float trackhitCorrY_
 
float trackhitCorrYF_
 
Hit trackhitF_
 
edm::ESWatcher< SiPixelTemplateDBObjectESProducerRcd > watchSiPixelTemplateRcd_
 

Additional Inherited Members

- Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using CacheTypes = CacheContexts< T... >
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T... >
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 
- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::EDPutTokenT< DQMToken > runToken_
 
unsigned int streamId_
 

Detailed Description

Description: generates the intermediate ALCAPROMPT dataset for the measurement of the SiPixel Lorentz Angle in the Prompt Calibration Loop Implementation: Books and fills 2D histograms of the drift vs depth in bins of pixel module rings to be fed into the SiPixelLorentzAnglePCLHarvester

Definition at line 104 of file SiPixelLorentzAnglePCLWorker.cc.

Constructor & Destructor Documentation

◆ SiPixelLorentzAnglePCLWorker()

SiPixelLorentzAnglePCLWorker::SiPixelLorentzAnglePCLWorker ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 209 of file SiPixelLorentzAnglePCLWorker.cc.

References Pixinfo::adc, bladeF_, bx_, chi2_, clust_, clustF_, Pixinfo::col, diskF_, eta_, event_, filename_, edm::ParameterSet::getParameter(), hFile_, isflipped_, ladder_, layer_, lumiblock_, module_, moduleF_, ndof_, notInPCL_, Pixinfo::npix, orbit_, panelF_, phi_, pixinfo_, pixinfoF_, pt_, qScale_, qScaleF_, rechit_, rechitCorr_, rechitCorrF_, rechitF_, Pixinfo::row, rQmQt_, rQmQtF_, run_, sideF_, SiPixelLorentzAngleTreeBarrel_, SiPixelLorentzAngleTreeForward_, t_trajTrack, trackhit_, trackhitCorrX_, trackhitCorrXF_, trackhitCorrY_, trackhitCorrYF_, trackhitF_, Pixinfo::x, and Pixinfo::y.

210  : analysisType_(convertStringToLorentzAngleAnalysisTypeEnum(iConfig.getParameter<std::string>("analysisType"))),
211  folder_(iConfig.getParameter<std::string>("folder")),
212  notInPCL_(iConfig.getParameter<bool>("notInPCL")),
213  filename_(iConfig.getParameter<std::string>("fileName")),
214  newmodulelist_(iConfig.getParameter<std::vector<std::string>>("newmodulelist")),
215  ptmin_(iConfig.getParameter<double>("ptMin")),
216  normChi2Max_(iConfig.getParameter<double>("normChi2Max")),
217  clustSizeYMin_(iConfig.getParameter<std::vector<int>>("clustSizeYMin")),
218  clustSizeXMax_(iConfig.getParameter<int>("clustSizeXMax")),
219  residualMax_(iConfig.getParameter<double>("residualMax")),
220  clustChargeMaxPerLength_(iConfig.getParameter<double>("clustChargeMaxPerLength")),
221  hist_depth_(iConfig.getParameter<int>("binsDepth")),
222  hist_drift_(iConfig.getParameter<int>("binsDrift")),
223  geomEsToken_(esConsumes<edm::Transition::BeginRun>()),
224  topoEsToken_(esConsumes<edm::Transition::BeginRun>()),
225  siPixelTemplateEsToken_(esConsumes<edm::Transition::BeginRun>()),
226  topoPerEventEsToken_(esConsumes()),
227  geomPerEventEsToken_(esConsumes()),
228  magneticFieldToken_(esConsumes()) {
229  t_trajTrack = consumes<TrajTrackAssociationCollection>(iConfig.getParameter<edm::InputTag>("src"));
230 
231  // now do what ever initialization is needed
232  int bufsize = 64000;
233  // create tree structure
234  // Barrel pixel
235  if (notInPCL_) {
236  hFile_ = std::make_unique<TFile>(filename_.c_str(), "RECREATE");
237  SiPixelLorentzAngleTreeBarrel_ =
238  std::make_unique<TTree>("SiPixelLorentzAngleTreeBarrel_", "SiPixel LorentzAngle tree barrel", bufsize);
239  SiPixelLorentzAngleTreeBarrel_->Branch("run", &run_, "run/I", bufsize);
240  SiPixelLorentzAngleTreeBarrel_->Branch("event", &event_, "event/l", bufsize);
241  SiPixelLorentzAngleTreeBarrel_->Branch("lumiblock", &lumiblock_, "lumiblock/I", bufsize);
242  SiPixelLorentzAngleTreeBarrel_->Branch("bx", &bx_, "bx/I", bufsize);
243  SiPixelLorentzAngleTreeBarrel_->Branch("orbit", &orbit_, "orbit/I", bufsize);
244  SiPixelLorentzAngleTreeBarrel_->Branch("module", &module_, "module/I", bufsize);
245  SiPixelLorentzAngleTreeBarrel_->Branch("ladder", &ladder_, "ladder/I", bufsize);
246  SiPixelLorentzAngleTreeBarrel_->Branch("layer", &layer_, "layer/I", bufsize);
247  SiPixelLorentzAngleTreeBarrel_->Branch("isflipped", &isflipped_, "isflipped/I", bufsize);
248  SiPixelLorentzAngleTreeBarrel_->Branch("pt", &pt_, "pt/F", bufsize);
249  SiPixelLorentzAngleTreeBarrel_->Branch("eta", &eta_, "eta/F", bufsize);
250  SiPixelLorentzAngleTreeBarrel_->Branch("phi", &phi_, "phi/F", bufsize);
251  SiPixelLorentzAngleTreeBarrel_->Branch("chi2", &chi2_, "chi2/D", bufsize);
252  SiPixelLorentzAngleTreeBarrel_->Branch("ndof", &ndof_, "ndof/D", bufsize);
253  SiPixelLorentzAngleTreeBarrel_->Branch("trackhit", &trackhit_, "x/F:y/F:alpha/D:beta/D:gamma_/D", bufsize);
254  SiPixelLorentzAngleTreeBarrel_->Branch("npix", &pixinfo_.npix, "npix/I", bufsize);
255  SiPixelLorentzAngleTreeBarrel_->Branch("rowpix", pixinfo_.row, "row[npix]/F", bufsize);
256  SiPixelLorentzAngleTreeBarrel_->Branch("colpix", pixinfo_.col, "col[npix]/F", bufsize);
257  SiPixelLorentzAngleTreeBarrel_->Branch("adc", pixinfo_.adc, "adc[npix]/F", bufsize);
258  SiPixelLorentzAngleTreeBarrel_->Branch("xpix", pixinfo_.x, "x[npix]/F", bufsize);
259  SiPixelLorentzAngleTreeBarrel_->Branch("ypix", pixinfo_.y, "y[npix]/F", bufsize);
260 
261  SiPixelLorentzAngleTreeBarrel_->Branch(
262  "clust",
263  &clust_,
264  "x/F:y/F:charge/F:size_x/I:size_y/I:maxPixelCol/I:maxPixelRow:minPixelCol/I:minPixelRow/I",
265  bufsize);
266  SiPixelLorentzAngleTreeBarrel_->Branch("rechit", &rechit_, "x/F:y/F", bufsize);
267  SiPixelLorentzAngleTreeBarrel_->Branch("rechit_corr", &rechitCorr_, "x/F:y/F", bufsize);
268  SiPixelLorentzAngleTreeBarrel_->Branch("trackhitcorr_x", &trackhitCorrX_, "trackhitcorr_x/F", bufsize);
269  SiPixelLorentzAngleTreeBarrel_->Branch("trackhitcorr_y", &trackhitCorrY_, "trackhitcorr_y/F", bufsize);
270  SiPixelLorentzAngleTreeBarrel_->Branch("qScale", &qScale_, "qScale/F", bufsize);
271  SiPixelLorentzAngleTreeBarrel_->Branch("rQmQt", &rQmQt_, "rQmQt/F", bufsize);
272  // Forward pixel
273 
274  SiPixelLorentzAngleTreeForward_ =
275  std::make_unique<TTree>("SiPixelLorentzAngleTreeForward_", "SiPixel LorentzAngle tree forward", bufsize);
276  SiPixelLorentzAngleTreeForward_->Branch("run", &run_, "run/I", bufsize);
277  SiPixelLorentzAngleTreeForward_->Branch("event", &event_, "event/l", bufsize);
278  SiPixelLorentzAngleTreeForward_->Branch("lumiblock", &lumiblock_, "lumiblock/I", bufsize);
279  SiPixelLorentzAngleTreeForward_->Branch("bx", &bx_, "bx/I", bufsize);
280  SiPixelLorentzAngleTreeForward_->Branch("orbit", &orbit_, "orbit/I", bufsize);
281  SiPixelLorentzAngleTreeForward_->Branch("side", &sideF_, "side/I", bufsize);
282  SiPixelLorentzAngleTreeForward_->Branch("disk", &diskF_, "disk/I", bufsize);
283  SiPixelLorentzAngleTreeForward_->Branch("blade", &bladeF_, "blade/I", bufsize);
284  SiPixelLorentzAngleTreeForward_->Branch("panel", &panelF_, "panel/I", bufsize);
285  SiPixelLorentzAngleTreeForward_->Branch("module", &moduleF_, "module/I", bufsize);
286  SiPixelLorentzAngleTreeForward_->Branch("pt", &pt_, "pt/F", bufsize);
287  SiPixelLorentzAngleTreeForward_->Branch("eta", &eta_, "eta/F", bufsize);
288  SiPixelLorentzAngleTreeForward_->Branch("phi", &phi_, "phi/F", bufsize);
289  SiPixelLorentzAngleTreeForward_->Branch("chi2", &chi2_, "chi2/D", bufsize);
290  SiPixelLorentzAngleTreeForward_->Branch("ndof", &ndof_, "ndof/D", bufsize);
291  SiPixelLorentzAngleTreeForward_->Branch("trackhit", &trackhitF_, "x/F:y/F:alpha/D:beta/D:gamma_/D", bufsize);
292  SiPixelLorentzAngleTreeForward_->Branch("npix", &pixinfoF_.npix, "npix/I", bufsize);
293  SiPixelLorentzAngleTreeForward_->Branch("rowpix", pixinfoF_.row, "row[npix]/F", bufsize);
294  SiPixelLorentzAngleTreeForward_->Branch("colpix", pixinfoF_.col, "col[npix]/F", bufsize);
295  SiPixelLorentzAngleTreeForward_->Branch("adc", pixinfoF_.adc, "adc[npix]/F", bufsize);
296  SiPixelLorentzAngleTreeForward_->Branch("xpix", pixinfoF_.x, "x[npix]/F", bufsize);
297  SiPixelLorentzAngleTreeForward_->Branch("ypix", pixinfoF_.y, "y[npix]/F", bufsize);
298 
299  SiPixelLorentzAngleTreeForward_->Branch(
300  "clust",
301  &clustF_,
302  "x/F:y/F:charge/F:size_x/I:size_y/I:maxPixelCol/I:maxPixelRow:minPixelCol/I:minPixelRow/I",
303  bufsize);
304  SiPixelLorentzAngleTreeForward_->Branch("rechit", &rechitF_, "x/F:y/F", bufsize);
305  SiPixelLorentzAngleTreeForward_->Branch("rechit_corr", &rechitCorrF_, "x/F:y/F", bufsize);
306  SiPixelLorentzAngleTreeForward_->Branch("trackhitcorr_x", &trackhitCorrXF_, "trackhitcorr_x/F", bufsize);
307  SiPixelLorentzAngleTreeForward_->Branch("trackhitcorr_y", &trackhitCorrYF_, "trackhitcorr_y/F", bufsize);
308  SiPixelLorentzAngleTreeForward_->Branch("qScale", &qScaleF_, "qScale/F", bufsize);
309  SiPixelLorentzAngleTreeForward_->Branch("rQmQt", &rQmQtF_, "rQmQt/F", bufsize);
310  }
311 }
deDxTools::esConsumes
ESGetTokenH3DDVariant esConsumes(std::string const &Record, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
Pixinfo::x
float x[maxpix]
Definition: SiPixelLorentzAngle.h:54
SiPixelLorentzAnglePCLWorker::topoPerEventEsToken_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoPerEventEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:198
Pixinfo::y
float y[maxpix]
Definition: SiPixelLorentzAngle.h:55
SiPixelLorentzAnglePCLWorker::geomPerEventEsToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomPerEventEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:199
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
SiPixelLorentzAnglePCLWorker::topoEsToken_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:196
SiPixelLorentzAnglePCLWorker::SiPixelLorentzAngleTreeBarrel_
std::unique_ptr< TTree > SiPixelLorentzAngleTreeBarrel_
Definition: SiPixelLorentzAnglePCLWorker.cc:191
SiPixelLorentzAnglePCLWorker::magneticFieldToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magneticFieldToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:200
SiPixelLorentzAnglePCLWorker::newmodulelist_
std::vector< std::string > newmodulelist_
Definition: SiPixelLorentzAnglePCLWorker.cc:137
SiPixelLorentzAnglePCLWorker::convertStringToLorentzAngleAnalysisTypeEnum
LorentzAngleAnalysisTypeEnum convertStringToLorentzAngleAnalysisTypeEnum(std::string type)
Definition: SiPixelLorentzAnglePCLWorker.cc:995
Pixinfo::adc
float adc[maxpix]
Definition: SiPixelLorentzAngle.h:53
SiPixelLorentzAnglePCLWorker::geomEsToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:195
Pixinfo::col
float col[maxpix]
Definition: SiPixelLorentzAngle.h:52
SiPixelLorentzAnglePCLWorker::siPixelTemplateEsToken_
edm::ESGetToken< SiPixelTemplateDBObject, SiPixelTemplateDBObjectESProducerRcd > siPixelTemplateEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:197
SiPixelLorentzAnglePCLWorker::t_trajTrack
edm::EDGetTokenT< TrajTrackAssociationCollection > t_trajTrack
Definition: SiPixelLorentzAnglePCLWorker.cc:203
SiPixelLorentzAnglePCLWorker::analysisType_
LorentzAngleAnalysisTypeEnum analysisType_
Definition: SiPixelLorentzAnglePCLWorker.cc:133
Pixinfo::row
float row[maxpix]
Definition: SiPixelLorentzAngle.h:51
SiPixelLorentzAnglePCLWorker::SiPixelLorentzAngleTreeForward_
std::unique_ptr< TTree > SiPixelLorentzAngleTreeForward_
Definition: SiPixelLorentzAnglePCLWorker.cc:192

◆ ~SiPixelLorentzAnglePCLWorker()

SiPixelLorentzAnglePCLWorker::~SiPixelLorentzAnglePCLWorker ( )
overridedefault

Member Function Documentation

◆ analyze()

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

Reimplemented from DQMEDAnalyzer.

Definition at line 319 of file SiPixelLorentzAnglePCLWorker.cc.

References funct::abs(), Pixinfo::adc, Hit::alpha, analysisType_, edm::AssociationMap< Tag >::begin(), Hit::beta, SiPixelLorentzAngleCalibrationHistograms::betaStartIdx_, mkfit::Config::Bfield, bladeF_, Surface::bounds(), SiPixelLorentzAngleCalibrationHistograms::BPixnewDetIds_, bx_, Clust::charge, chi2_, Trajectory::chiSquared(), clust_, clustChargeMaxPerLength_, SiPixelRecHit::cluster(), clustF_, clustSizeXMax_, clustSizeYMin_, siPixelLACalibration::cmToum, Pixinfo::col, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, hcalRecHitTable_cff::depth, hcalRecHitTable_cff::detId, diskF_, HLT_2023v12_cff::distance, shallow::drift(), PVValHelper::dx, PVValHelper::dy, eGrazingAngle, eMinimumClusterSize, edm::AssociationMap< Tag >::empty(), edm::AssociationMap< Tag >::end(), eta_, event_, dqm::impl::MonitorElement::Fill(), fillPix(), spr::find(), Hit::gamma, geomPerEventEsToken_, edm::EventSetup::getData(), SiPixelTemplateDBObject::getTemplateID(), SiPixelLorentzAngleCalibrationHistograms::h_bySectOccupancy_, SiPixelLorentzAngleCalibrationHistograms::h_drift_depth_adc2_, SiPixelLorentzAngleCalibrationHistograms::h_drift_depth_adc_, SiPixelLorentzAngleCalibrationHistograms::h_drift_depth_noadc_, SiPixelLorentzAngleCalibrationHistograms::h_fpixAngleSize_, SiPixelLorentzAngleCalibrationHistograms::h_fpixMagField_, SiPixelLorentzAngleCalibrationHistograms::h_trackChi2_, SiPixelLorentzAngleCalibrationHistograms::h_trackEta_, SiPixelLorentzAngleCalibrationHistograms::h_trackPhi_, SiPixelLorentzAngleCalibrationHistograms::h_trackPt_, SiPixelLorentzAngleCalibrationHistograms::h_tracks_, heavyIonCSV_trainingSettings::idx, iEvent, iHists, SiPixelTemplate::interpolate(), MagneticField::inTesla(), isflipped_, TrajectoryStateOnSurface::isValid(), dqmiolumiharvest::j, ladder_, layer_, TrajectoryStateOnSurface::localDirection(), BaseTrackerRecHit::localPosition(), TrajectoryStateOnSurface::localPosition(), lumiblock_, M_PI, magneticFieldToken_, Clust::maxPixelCol, Clust::maxPixelRow, Trajectory::measurements(), Clust::minPixelCol, Clust::minPixelRow, module_, moduleF_, Trajectory::ndof(), ndof_, SiPixelLorentzAngleCalibrationHistograms::nlay, SiPixelLorentzAngleCalibrationHistograms::nModules_, MagneticField::nominalValue(), normChi2Max_, notInPCL_, SiPixelLorentzAngleCalibrationHistograms::nPanels_, Pixinfo::npix, SiPixelLorentzAngleCalibrationHistograms::nSides_, orbit_, panelF_, PV3DBase< T, PVType, FrameType >::perp(), TrackerGeometry::Ph1PXB, phi_, PixelTopology::pitch(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixinfo_, pixinfoF_, GloballyPositioned< T >::position(), conifer::pow(), pt_, ptmin_, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), TrackerTopology::pxfBlade(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), TrackerTopology::pxfSide(), SiPixelTemplate::qscale(), qScale_, qScaleF_, SiPixelTemplate::r_qMeas_qTrue(), DetId::rawId(), rpcPointValidation_cfi::recHit, rechit_, rechitCorr_, rechitCorrF_, rechitF_, residualMax_, Pixinfo::row, rQmQt_, rQmQtF_, run_, sideF_, SiPixelLorentzAngleTreeBarrel_, SiPixelLorentzAngleTreeForward_, Clust::size_x, Clust::size_y, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), GeomDet::surface(), surface_deformation(), t_trajTrack, funct::tan(), templateDBobject_, thePixelTemp_, Bounds::thickness(), Surface::toGlobal(), GloballyPositioned< T >::toLocal(), topoPerEventEsToken_, HLT_2023v12_cff::track, DetId::Tracker, PbPb_ZMuSkimMuonDPG_cff::tracker, trackhit_, trackhitCorrX_, trackhitCorrXF_, trackhitCorrY_, trackhitCorrYF_, trackhitF_, x, Pixinfo::x, Hit::x, PV3DBase< T, PVType, FrameType >::x(), Clust::x, Rechit::x, y, Pixinfo::y, Hit::y, PV3DBase< T, PVType, FrameType >::y(), Clust::y, Rechit::y, and PV3DBase< T, PVType, FrameType >::z().

319  {
320  // Retrieve tracker topology from geometry
321  const TrackerTopology* const tTopo = &iSetup.getData(topoPerEventEsToken_);
322 
323  // Retrieve track geometry
324  const TrackerGeometry* tracker = &iSetup.getData(geomPerEventEsToken_);
325 
326  // Retrieve magnetic field
327  const MagneticField* magField = &iSetup.getData(magneticFieldToken_);
328 
329  // get the association map between tracks and trajectories
330  edm::Handle<TrajTrackAssociationCollection> trajTrackCollectionHandle;
331  iEvent.getByToken(t_trajTrack, trajTrackCollectionHandle);
332 
333  module_ = -1;
334  layer_ = -1;
335  ladder_ = -1;
336  isflipped_ = -1;
337  pt_ = -999;
338  eta_ = 999;
339  phi_ = 999;
340  pixinfo_.npix = 0;
341 
342  run_ = iEvent.id().run();
343  event_ = iEvent.id().event();
344  lumiblock_ = iEvent.luminosityBlock();
345  bx_ = iEvent.bunchCrossing();
346  orbit_ = iEvent.orbitNumber();
347 
348  if (!trajTrackCollectionHandle->empty()) {
349  for (TrajTrackAssociationCollection::const_iterator it = trajTrackCollectionHandle->begin();
350  it != trajTrackCollectionHandle->end();
351  ++it) {
352  const reco::Track& track = *it->val;
353  const Trajectory& traj = *it->key;
354 
355  // get the trajectory measurements
356  std::vector<TrajectoryMeasurement> tmColl = traj.measurements();
357  pt_ = track.pt();
358  eta_ = track.eta();
359  phi_ = track.phi();
360  chi2_ = traj.chiSquared();
361  ndof_ = traj.ndof();
362 
363  if (pt_ < ptmin_)
364  continue;
365 
366  iHists.h_trackEta_->Fill(eta_);
367  iHists.h_trackPhi_->Fill(phi_);
368  iHists.h_trackPt_->Fill(pt_);
369  iHists.h_trackChi2_->Fill(chi2_ / ndof_);
370  iHists.h_tracks_->Fill(0);
371  bool pixeltrack = false;
372 
373  // iterate over trajectory measurements
374  for (const auto& itTraj : tmColl) {
375  if (!itTraj.updatedState().isValid())
376  continue;
377  const TransientTrackingRecHit::ConstRecHitPointer& recHit = itTraj.recHit();
378  if (!recHit->isValid() || recHit->geographicalId().det() != DetId::Tracker)
379  continue;
380  unsigned int subDetID = (recHit->geographicalId().subdetId());
381  if (subDetID == PixelSubdetector::PixelBarrel || subDetID == PixelSubdetector::PixelEndcap) {
382  if (!pixeltrack) {
383  iHists.h_tracks_->Fill(1);
384  }
385  pixeltrack = true;
386  }
387 
388  if (subDetID == PixelSubdetector::PixelBarrel) {
389  DetId detIdObj = recHit->geographicalId();
390  const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(tracker->idToDet(detIdObj));
391  if (!theGeomDet)
392  continue;
393 
394  const PixelTopology* topol = &(theGeomDet->specificTopology());
395 
396  float ypitch_ = topol->pitch().second;
397  float width_ = theGeomDet->surface().bounds().thickness();
398 
399  if (!topol)
400  continue;
401 
402  layer_ = tTopo->pxbLayer(detIdObj);
403  ladder_ = tTopo->pxbLadder(detIdObj);
404  module_ = tTopo->pxbModule(detIdObj);
405 
406  float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 0.)).perp();
407  float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 1.)).perp();
408 
409  isflipped_ = (tmp2 < tmp1) ? 1 : 0;
410 
411  const SiPixelRecHit* recHitPix = dynamic_cast<const SiPixelRecHit*>((*recHit).hit());
412  if (!recHitPix)
413  continue;
414  rechit_.x = recHitPix->localPosition().x();
415  rechit_.y = recHitPix->localPosition().y();
416  SiPixelRecHit::ClusterRef const& cluster = recHitPix->cluster();
417 
418  pixinfo_ = fillPix(*cluster, topol);
419 
420  // fill entries in clust_
421 
422  clust_.x = (cluster)->x();
423  clust_.y = (cluster)->y();
424  clust_.charge = (cluster->charge()) / 1000.; // clust_.charge: in the unit of 1000e
425  clust_.size_x = cluster->sizeX();
426  clust_.size_y = cluster->sizeY();
427  clust_.maxPixelCol = cluster->maxPixelCol();
428  clust_.maxPixelRow = cluster->maxPixelRow();
429  clust_.minPixelCol = cluster->minPixelCol();
430  clust_.minPixelRow = cluster->minPixelRow();
431 
432  // fill the trackhit info
433  TrajectoryStateOnSurface tsos = itTraj.updatedState();
434  if (!tsos.isValid()) {
435  edm::LogWarning("SiPixelLorentzAnglePCLWorker") << "tsos not valid";
436  continue;
437  }
438  LocalVector trackdirection = tsos.localDirection();
439  LocalPoint trackposition = tsos.localPosition();
440 
441  if (trackdirection.z() == 0)
442  continue;
443  // the local position and direction
444  trackhit_.alpha = atan2(trackdirection.z(), trackdirection.x());
445  trackhit_.beta = atan2(trackdirection.z(), trackdirection.y());
446  trackhit_.gamma = atan2(trackdirection.x(), trackdirection.y());
447  trackhit_.x = trackposition.x();
448  trackhit_.y = trackposition.y();
449 
450  // get qScale_ = templ.qscale() and templ.r_qMeas_qTrue();
451  float cotalpha = trackdirection.x() / trackdirection.z();
452  float cotbeta = trackdirection.y() / trackdirection.z();
453  float cotbeta_min = clustSizeYMin_[layer_ - 1] * ypitch_ / width_;
454  if (std::abs(cotbeta) <= cotbeta_min)
455  continue;
456  double drdz = sqrt(1. + cotalpha * cotalpha + cotbeta * cotbeta);
457  double clusterCharge_cut = clustChargeMaxPerLength_ * drdz;
458 
459  auto detId = detIdObj.rawId();
460  int DetId_index = -1;
461 
462  const auto& newModIt = (std::find(iHists.BPixnewDetIds_.begin(), iHists.BPixnewDetIds_.end(), detId));
463  bool isNewMod = (newModIt != iHists.BPixnewDetIds_.end());
464  if (isNewMod) {
465  DetId_index = std::distance(iHists.BPixnewDetIds_.begin(), newModIt);
466  }
467 
468  if (notInPCL_) {
469  // fill the template from the store (from dqmBeginRun)
470  SiPixelTemplate theTemplate(thePixelTemp_);
471 
472  float locBx = (cotbeta < 0.) ? -1 : 1.;
473  float locBz = (cotalpha < 0.) ? -locBx : locBx;
474 
475  int TemplID = templateDBobject_->getTemplateID(detId);
476  theTemplate.interpolate(TemplID, cotalpha, cotbeta, locBz, locBx);
477  qScale_ = theTemplate.qscale();
478  rQmQt_ = theTemplate.r_qMeas_qTrue();
479  }
480 
481  // Surface deformation
482  const auto& lp_pair = surface_deformation(topol, tsos, recHitPix);
483 
484  LocalPoint lp_track = lp_pair.first;
485  LocalPoint lp_rechit = lp_pair.second;
486 
487  rechitCorr_.x = lp_rechit.x();
488  rechitCorr_.y = lp_rechit.y();
489  trackhitCorrX_ = lp_track.x();
490  trackhitCorrY_ = lp_track.y();
491 
492  if (notInPCL_) {
493  SiPixelLorentzAngleTreeBarrel_->Fill();
494  }
495 
496  if (analysisType_ != eGrazingAngle)
497  continue;
498  // is one pixel in cluster a large pixel ? (hit will be excluded)
499  bool large_pix = false;
500  for (int j = 0; j < pixinfo_.npix; j++) {
501  int colpos = static_cast<int>(pixinfo_.col[j]);
502  if (pixinfo_.row[j] == 0 || pixinfo_.row[j] == 79 || pixinfo_.row[j] == 80 || pixinfo_.row[j] == 159 ||
503  colpos % 52 == 0 || colpos % 52 == 51) {
504  large_pix = true;
505  }
506  }
507 
508  double residualsq = (trackhitCorrX_ - rechitCorr_.x) * (trackhitCorrX_ - rechitCorr_.x) +
509  (trackhitCorrY_ - rechitCorr_.y) * (trackhitCorrY_ - rechitCorr_.y);
510 
511  double xlim1 = trackhitCorrX_ - width_ * cotalpha / 2.;
512  double hypitch_ = ypitch_ / 2.;
513  double ylim1 = trackhitCorrY_ - width_ * cotbeta / 2.;
514  double ylim2 = trackhitCorrY_ + width_ * cotbeta / 2.;
515 
516  int clustSizeY_cut = clustSizeYMin_[layer_ - 1];
517 
518  if (!large_pix && (chi2_ / ndof_) < normChi2Max_ && cluster->sizeY() >= clustSizeY_cut &&
519  residualsq < residualMax_ * residualMax_ && cluster->charge() < clusterCharge_cut &&
520  cluster->sizeX() < clustSizeXMax_) {
521  // iterate over pixels in hit
522  for (int j = 0; j < pixinfo_.npix; j++) {
523  // use trackhits and include bowing correction
524  float ypixlow = pixinfo_.y[j] - hypitch_;
525  float ypixhigh = pixinfo_.y[j] + hypitch_;
526  if (cotbeta > 0.) {
527  if (ylim1 > ypixlow)
528  ypixlow = ylim1;
529  if (ylim2 < ypixhigh)
530  ypixhigh = ylim2;
531  } else {
532  if (ylim2 > ypixlow)
533  ypixlow = ylim2;
534  if (ylim1 < ypixhigh)
535  ypixhigh = ylim1;
536  }
537  float ypixavg = 0.5f * (ypixlow + ypixhigh);
538 
539  float dx = (pixinfo_.x[j] - xlim1) * siPixelLACalibration::cmToum; // dx: in the unit of micrometer
540  float dy = (ypixavg - ylim1) * siPixelLACalibration::cmToum; // dy: in the unit of micrometer
541  float depth = dy * tan(trackhit_.beta);
542  float drift = dx - dy * tan(trackhit_.gamma);
543 
544  if (isNewMod == false) {
545  int i_index = module_ + (layer_ - 1) * iHists.nModules_[layer_ - 1];
546  iHists.h_drift_depth_adc_[i_index]->Fill(drift, depth, pixinfo_.adc[j]);
547  iHists.h_drift_depth_adc2_[i_index]->Fill(drift, depth, pixinfo_.adc[j] * pixinfo_.adc[j]);
548  iHists.h_drift_depth_noadc_[i_index]->Fill(drift, depth, 1.);
549  iHists.h_bySectOccupancy_->Fill(i_index - 1); // histogram starts at 0
550 
551  if (tracker->getDetectorType(subDetID) == TrackerGeometry::ModuleType::Ph1PXB) {
552  if ((module_ == 3 || module_ == 5) && (layer_ == 3 || layer_ == 4)) {
553  int i_index_merge = i_index + 1;
554  iHists.h_drift_depth_adc_[i_index_merge]->Fill(drift, depth, pixinfo_.adc[j]);
555  iHists.h_drift_depth_adc2_[i_index_merge]->Fill(drift, depth, pixinfo_.adc[j] * pixinfo_.adc[j]);
556  iHists.h_drift_depth_noadc_[i_index_merge]->Fill(drift, depth, 1.);
557  iHists.h_bySectOccupancy_->Fill(i_index_merge - 1);
558  }
559  if ((module_ == 4 || module_ == 6) && (layer_ == 3 || layer_ == 4)) {
560  int i_index_merge = i_index - 1;
561  iHists.h_drift_depth_adc_[i_index_merge]->Fill(drift, depth, pixinfo_.adc[j]);
562  iHists.h_drift_depth_adc2_[i_index_merge]->Fill(drift, depth, pixinfo_.adc[j] * pixinfo_.adc[j]);
563  iHists.h_drift_depth_noadc_[i_index_merge]->Fill(drift, depth, 1.);
564  iHists.h_bySectOccupancy_->Fill(i_index_merge - 1);
565  }
566  }
567 
568  } else {
569  int new_index = iHists.nModules_[iHists.nlay - 1] +
570  (iHists.nlay - 1) * iHists.nModules_[iHists.nlay - 1] + 1 + DetId_index;
571 
572  iHists.h_drift_depth_adc_[new_index]->Fill(drift, depth, pixinfo_.adc[j]);
573  iHists.h_drift_depth_adc2_[new_index]->Fill(drift, depth, pixinfo_.adc[j] * pixinfo_.adc[j]);
574  iHists.h_drift_depth_noadc_[new_index]->Fill(drift, depth, 1.);
575  iHists.h_bySectOccupancy_->Fill(new_index - 1); // histogram starts at 0
576  }
577  }
578  }
579  } else if (subDetID == PixelSubdetector::PixelEndcap) {
580  DetId detIdObj = recHit->geographicalId();
581  const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(tracker->idToDet(detIdObj));
582  if (!theGeomDet)
583  continue;
584 
585  const PixelTopology* topol = &(theGeomDet->specificTopology());
586 
587  if (!topol)
588  continue;
589 
590  sideF_ = tTopo->pxfSide(detIdObj);
591  diskF_ = tTopo->pxfDisk(detIdObj);
592  bladeF_ = tTopo->pxfBlade(detIdObj);
593  panelF_ = tTopo->pxfPanel(detIdObj);
594  moduleF_ = tTopo->pxfModule(detIdObj);
595 
596  const SiPixelRecHit* recHitPix = dynamic_cast<const SiPixelRecHit*>((*recHit).hit());
597  if (!recHitPix)
598  continue;
599  rechitF_.x = recHitPix->localPosition().x();
600  rechitF_.y = recHitPix->localPosition().y();
601  SiPixelRecHit::ClusterRef const& cluster = recHitPix->cluster();
602 
603  pixinfoF_ = fillPix(*cluster, topol);
604 
605  // fill entries in clust_
606 
607  clustF_.x = (cluster)->x();
608  clustF_.y = (cluster)->y();
609  clustF_.charge = (cluster->charge()) / 1000.; // clustF_.charge: in the unit of 1000e
610  clustF_.size_x = cluster->sizeX();
611  clustF_.size_y = cluster->sizeY();
612  clustF_.maxPixelCol = cluster->maxPixelCol();
613  clustF_.maxPixelRow = cluster->maxPixelRow();
614  clustF_.minPixelCol = cluster->minPixelCol();
615  clustF_.minPixelRow = cluster->minPixelRow();
616 
617  // fill the trackhit info
618  TrajectoryStateOnSurface tsos = itTraj.updatedState();
619  if (!tsos.isValid()) {
620  edm::LogWarning("SiPixelLorentzAnglePCLWorker") << "tsos not valid";
621  continue;
622  }
623  LocalVector trackdirection = tsos.localDirection();
624  LocalPoint trackposition = tsos.localPosition();
625 
626  if (trackdirection.z() == 0)
627  continue;
628  // the local position and direction
629  trackhitF_.alpha = atan2(trackdirection.z(), trackdirection.x());
630  trackhitF_.beta = atan2(trackdirection.z(), trackdirection.y());
631  trackhitF_.gamma = atan2(trackdirection.x(), trackdirection.y());
632  trackhitF_.x = trackposition.x();
633  trackhitF_.y = trackposition.y();
634 
635  float cotalpha = trackdirection.x() / trackdirection.z();
636  float cotbeta = trackdirection.y() / trackdirection.z();
637 
638  auto detId = detIdObj.rawId();
639 
640  if (notInPCL_) {
641  // fill the template from the store (from dqmBeginRun)
642  SiPixelTemplate theTemplate(thePixelTemp_);
643 
644  float locBx = (cotbeta < 0.) ? -1 : 1.;
645  float locBz = (cotalpha < 0.) ? -locBx : locBx;
646 
647  int TemplID = templateDBobject_->getTemplateID(detId);
648  theTemplate.interpolate(TemplID, cotalpha, cotbeta, locBz, locBx);
649  qScaleF_ = theTemplate.qscale();
650  rQmQtF_ = theTemplate.r_qMeas_qTrue();
651  }
652 
653  // Surface deformation
654  const auto& lp_pair = surface_deformation(topol, tsos, recHitPix);
655 
656  LocalPoint lp_track = lp_pair.first;
657  LocalPoint lp_rechit = lp_pair.second;
658 
659  rechitCorrF_.x = lp_rechit.x();
660  rechitCorrF_.y = lp_rechit.y();
661  trackhitCorrXF_ = lp_track.x();
662  trackhitCorrYF_ = lp_track.y();
663  if (notInPCL_) {
664  SiPixelLorentzAngleTreeForward_->Fill();
665  }
666 
667  if (analysisType_ != eMinimumClusterSize)
668  continue;
669 
670  int theMagField = magField->nominalValue();
671  if (theMagField < 37 || theMagField > 39)
672  continue;
673 
674  double chi2_ndof = chi2_ / ndof_;
675  if (chi2_ndof >= normChi2Max_)
676  continue;
677 
678  //--- large pixel cut
679  bool large_pix = false;
680  for (int j = 0; j < pixinfoF_.npix; j++) {
681  int colpos = static_cast<int>(pixinfoF_.col[j]);
682  if (pixinfoF_.row[j] == 0 || pixinfoF_.row[j] == 79 || pixinfoF_.row[j] == 80 || pixinfoF_.row[j] == 159 ||
683  colpos % 52 == 0 || colpos % 52 == 51) {
684  large_pix = true;
685  }
686  }
687 
688  if (large_pix)
689  continue;
690 
691  //--- residual cut
692  double residual = sqrt(pow(trackhitCorrXF_ - rechitCorrF_.x, 2) + pow(trackhitCorrYF_ - rechitCorrF_.y, 2));
693 
694  if (residual > residualMax_)
695  continue;
696 
697  int ringIdx = bladeF_ <= 22 ? 0 : 1;
698  int panelIdx = panelF_ - 1;
699  int sideIdx = sideF_ - 1;
700  int idx = iHists.nSides_ * iHists.nPanels_ * ringIdx + iHists.nSides_ * panelIdx + sideIdx;
701  int idxBeta = iHists.betaStartIdx_ + idx;
702 
703  double cotanAlpha = std::tan(M_PI / 2. - trackhitF_.alpha);
704  double cotanBeta = std::tan(M_PI / 2. - trackhitF_.beta);
705 
706  LocalVector Bfield = theGeomDet->surface().toLocal(magField->inTesla(theGeomDet->surface().position()));
707  iHists.h_fpixMagField_[0][idx]->Fill(Bfield.x());
708  iHists.h_fpixMagField_[1][idx]->Fill(Bfield.y());
709  iHists.h_fpixMagField_[2][idx]->Fill(Bfield.z());
710 
711  if (clustF_.size_y >= 2) {
712  iHists.h_fpixAngleSize_[idx]->Fill(cotanAlpha, clustF_.size_x);
713  }
714 
715  if (clust_.size_x >= 0) {
716  iHists.h_fpixAngleSize_[idxBeta]->Fill(cotanBeta, clustF_.size_y);
717  }
718  }
719  } //end iteration over trajectory measurements
720  } //end iteration over trajectories
721  }
722 }
SiPixelRecHit::cluster
ClusterRef cluster() const
Definition: SiPixelRecHit.h:47
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:7
TrackerTopology::pxbLayer
unsigned int pxbLayer(const DetId &id) const
Definition: TrackerTopology.h:146
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
SiPixelLorentzAnglePCLWorker::thePixelTemp_
std::vector< SiPixelTemplateStore > thePixelTemp_
Definition: SiPixelLorentzAnglePCLWorker.cc:131
Pixinfo::x
float x[maxpix]
Definition: SiPixelLorentzAngle.h:54
TrackerTopology::pxfBlade
unsigned int pxfBlade(const DetId &id) const
Definition: TrackerTopology.h:483
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
MagneticField::inTesla
virtual GlobalVector inTesla(const GlobalPoint &gp) const =0
Field value ad specified global point, in Tesla.
SiPixelLorentzAnglePCLWorker::topoPerEventEsToken_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoPerEventEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:198
Pixinfo::y
float y[maxpix]
Definition: SiPixelLorentzAngle.h:55
TrackerTopology::pxfModule
unsigned int pxfModule(const DetId &id) const
Definition: TrackerTopology.h:165
shallow::drift
LocalVector drift(const StripGeomDetUnit *, const MagneticField &, const SiStripLorentzAngle &)
Definition: ShallowTools.cc:36
SiPixelLorentzAnglePCLWorker::geomPerEventEsToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomPerEventEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:199
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:24
Trajectory::chiSquared
float chiSquared() const
Definition: Trajectory.h:241
TrackerTopology::pxbLadder
unsigned int pxbLadder(const DetId &id) const
Definition: TrackerTopology.h:157
edm::AssociationMap::empty
bool empty() const
return true if empty
Definition: AssociationMap.h:166
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
GloballyPositioned::toLocal
LocalPoint toLocal(const GlobalPoint &gp) const
Definition: GloballyPositioned.h:98
SiPixelLorentzAnglePCLWorker::fillPix
const Pixinfo fillPix(const SiPixelCluster &LocPix, const PixelTopology *topol) const
Definition: SiPixelLorentzAnglePCLWorker.cc:960
PbPb_ZMuSkimMuonDPG_cff.tracker
tracker
Definition: PbPb_ZMuSkimMuonDPG_cff.py:60
Trajectory::measurements
DataContainer const & measurements() const
Definition: Trajectory.h:178
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:173
SiPixelLorentzAnglePCLWorker::surface_deformation
const std::pair< LocalPoint, LocalPoint > surface_deformation(const PixelTopology *topol, TrajectoryStateOnSurface &tsos, const SiPixelRecHit *recHitPix) const
Definition: SiPixelLorentzAnglePCLWorker.cc:978
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:292
Bounds::thickness
virtual float thickness() const =0
Hit::beta
double beta
Definition: SiPixelLorentzAngle.h:61
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
iEvent
int iEvent
Definition: GenABIO.cc:224
Trajectory::ndof
int ndof(bool bon=true) const
Definition: Trajectory.cc:97
SiPixelLorentzAnglePCLWorker::templateDBobject_
const SiPixelTemplateDBObject * templateDBobject_
Definition: SiPixelLorentzAnglePCLWorker.cc:130
TrackingRecHit::ConstRecHitPointer
std::shared_ptr< TrackingRecHit const > ConstRecHitPointer
Definition: TrackingRecHit.h:25
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
TrajectoryStateOnSurface::localDirection
LocalVector localDirection() const
Definition: TrajectoryStateOnSurface.h:76
SiPixelLorentzAnglePCLWorker::SiPixelLorentzAngleTreeBarrel_
std::unique_ptr< TTree > SiPixelLorentzAngleTreeBarrel_
Definition: SiPixelLorentzAnglePCLWorker.cc:191
mkfit::Config::Bfield
constexpr float Bfield
Definition: Config.h:53
SiPixelLorentzAnglePCLWorker::magneticFieldToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magneticFieldToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:200
TrackerTopology::pxfDisk
unsigned int pxfDisk(const DetId &id) const
Definition: TrackerTopology.h:482
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
Pixinfo::adc
float adc[maxpix]
Definition: SiPixelLorentzAngle.h:53
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
TrackerTopology::pxfPanel
unsigned int pxfPanel(const DetId &id) const
Definition: TrackerTopology.h:486
DetId
Definition: DetId.h:17
TrackerTopology::pxfSide
unsigned int pxfSide(const DetId &id) const
Definition: TrackerTopology.h:194
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
GloballyPositioned::position
const PositionType & position() const
Definition: GloballyPositioned.h:36
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
Hit::alpha
double alpha
Definition: SiPixelLorentzAngle.h:60
SiPixelLorentzAnglePCLWorker::iHists
SiPixelLorentzAngleCalibrationHistograms iHists
Definition: SiPixelLorentzAnglePCLWorker.cc:126
SiPixelTemplateDBObject::getTemplateID
short getTemplateID(const uint32_t &detid) const
Definition: SiPixelTemplateDBObject.h:131
Pixinfo::col
float col[maxpix]
Definition: SiPixelLorentzAngle.h:52
edm::AssociationMap::begin
const_iterator begin() const
first iterator over the map (read only)
Definition: AssociationMap.h:171
hcalRecHitTable_cff.depth
depth
Definition: hcalRecHitTable_cff.py:17
PixelGeomDetUnit::specificTopology
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
Definition: PixelGeomDetUnit.cc:17
BaseTrackerRecHit::localPosition
LocalPoint localPosition() const override
Definition: BaseTrackerRecHit.h:56
SiPixelLorentzAnglePCLWorker::t_trajTrack
edm::EDGetTokenT< TrajTrackAssociationCollection > t_trajTrack
Definition: SiPixelLorentzAnglePCLWorker.cc:203
MagneticField::nominalValue
int nominalValue() const
The nominal field value for this map in kGauss.
Definition: MagneticField.h:49
Hit::gamma
double gamma
Definition: SiPixelLorentzAngle.h:62
PixelTopology::pitch
virtual std::pair< float, float > pitch() const =0
TrackerTopology::pxbModule
unsigned int pxbModule(const DetId &id) const
Definition: TrackerTopology.h:162
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:128
SiPixelLorentzAnglePCLWorker::analysisType_
LorentzAngleAnalysisTypeEnum analysisType_
Definition: SiPixelLorentzAnglePCLWorker.cc:133
Pixinfo::row
float row[maxpix]
Definition: SiPixelLorentzAngle.h:51
HLT_2023v12_cff.track
track
Definition: HLT_2023v12_cff.py:9604
HLT_2023v12_cff.distance
distance
Definition: HLT_2023v12_cff.py:5849
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
SiPixelLorentzAnglePCLWorker::SiPixelLorentzAngleTreeForward_
std::unique_ptr< TTree > SiPixelLorentzAngleTreeForward_
Definition: SiPixelLorentzAnglePCLWorker.cc:192
Surface::bounds
const Bounds & bounds() const
Definition: Surface.h:87
Local3DPoint
Point3DBase< float, LocalTag > Local3DPoint
Definition: LocalPoint.h:9

◆ bookHistograms()

void SiPixelLorentzAnglePCLWorker::bookHistograms ( DQMStore::IBooker &  iBooker,
edm::Run const &  run,
edm::EventSetup const &  iSetup 
)
overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 773 of file SiPixelLorentzAnglePCLWorker.cc.

References analysisType_, SiPixelLorentzAngleCalibrationHistograms::betaStartIdx_, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), SiPixelLorentzAngleCalibrationHistograms::BPixnewDetIds_, SiPixelLorentzAngleCalibrationHistograms::BPixnewmodulename_, eGrazingAngle, folder_, SiPixelLorentzAngleCalibrationHistograms::h_bySectOccupancy_, SiPixelLorentzAngleCalibrationHistograms::h_drift_depth_, SiPixelLorentzAngleCalibrationHistograms::h_drift_depth_adc2_, SiPixelLorentzAngleCalibrationHistograms::h_drift_depth_adc_, SiPixelLorentzAngleCalibrationHistograms::h_drift_depth_noadc_, SiPixelLorentzAngleCalibrationHistograms::h_fpixAngleSize_, SiPixelLorentzAngleCalibrationHistograms::h_fpixMagField_, SiPixelLorentzAngleCalibrationHistograms::h_fpixMean_, SiPixelLorentzAngleCalibrationHistograms::h_mean_, SiPixelLorentzAngleCalibrationHistograms::h_trackChi2_, SiPixelLorentzAngleCalibrationHistograms::h_trackEta_, SiPixelLorentzAngleCalibrationHistograms::h_trackPhi_, SiPixelLorentzAngleCalibrationHistograms::h_trackPt_, SiPixelLorentzAngleCalibrationHistograms::h_tracks_, hist_depth_, hist_drift_, mps_fire::i, heavyIonCSV_trainingSettings::idx, iHists, createfilelist::int, LogDebug, visualization-live-secondInstance_cfg::m, M_PI, Skims_PA_cff::name, SiPixelLorentzAngleCalibrationHistograms::nlay, SiPixelLorentzAngleCalibrationHistograms::nModules_, SiPixelLorentzAngleCalibrationHistograms::nPanels_, SiPixelLorentzAngleCalibrationHistograms::nRings_, SiPixelLorentzAngleCalibrationHistograms::nSides_, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, alignCSCRings::s, dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

775  {
776  std::string name;
777  std::string title;
778  if (analysisType_ == eGrazingAngle) {
779  // book the by partition monitoring
780  const auto maxSect = iHists.nlay * iHists.nModules_[iHists.nlay - 1] + (int)iHists.BPixnewDetIds_.size();
781 
782  iBooker.setCurrentFolder(fmt::sprintf("%s/SectorMonitoring", folder_.data()));
783  iHists.h_bySectOccupancy_ = iBooker.book1D(
784  "h_bySectorOccupancy", "hit occupancy by sector;pixel sector;hits on track", maxSect, -0.5, maxSect - 0.5);
785 
786  iBooker.setCurrentFolder(folder_);
787  static constexpr double min_depth_ = -100.;
788  static constexpr double max_depth_ = 400.;
789  static constexpr double min_drift_ = -500.;
790  static constexpr double max_drift_ = 500.;
791 
792  // book the mean values projections and set the bin names of the by sector monitoring
793  for (int i_layer = 1; i_layer <= iHists.nlay; i_layer++) {
794  for (int i_module = 1; i_module <= iHists.nModules_[i_layer - 1]; i_module++) {
795  unsigned int i_index = i_module + (i_layer - 1) * iHists.nModules_[i_layer - 1];
796  std::string binName = fmt::sprintf("BPix Lay%i Mod%i", i_layer, i_module);
797  LogDebug("SiPixelLorentzAnglePCLWorker") << " i_index: " << i_index << " bin name: " << binName
798  << " (i_layer: " << i_layer << " i_module:" << i_module << ")";
799 
800  iHists.h_bySectOccupancy_->setBinLabel(i_index, binName);
801 
802  name = fmt::sprintf("h_mean_layer%i_module%i", i_layer, i_module);
803  title = fmt::sprintf(
804  "average drift vs depth layer%i module%i; production depth [#mum]; #LTdrift#GT [#mum]", i_layer, i_module);
805  iHists.h_mean_[i_index] = iBooker.book1D(name, title, hist_depth_, min_depth_, max_depth_);
806  }
807  }
808  for (int i = 0; i < (int)iHists.BPixnewDetIds_.size(); i++) {
809  name = fmt::sprintf("h_BPixnew_mean_%s", iHists.BPixnewmodulename_[i].c_str());
810  title = fmt::sprintf("average drift vs depth %s; production depth [#mum]; #LTdrift#GT [#mum]",
811  iHists.BPixnewmodulename_[i].c_str());
812  int new_index = iHists.nModules_[iHists.nlay - 1] + (iHists.nlay - 1) * iHists.nModules_[iHists.nlay - 1] + 1 + i;
813  iHists.h_mean_[new_index] = iBooker.book1D(name, title, hist_depth_, min_depth_, max_depth_);
814 
815  LogDebug("SiPixelLorentzAnglePCLWorker")
816  << "i_index" << new_index << " bin name: " << iHists.BPixnewmodulename_[i];
817 
818  iHists.h_bySectOccupancy_->setBinLabel(new_index, iHists.BPixnewmodulename_[i]);
819  }
820 
821  //book the 2D histograms
822  for (int i_layer = 1; i_layer <= iHists.nlay; i_layer++) {
823  iBooker.setCurrentFolder(fmt::sprintf("%s/BPix/BPixLayer%i", folder_.data(), i_layer));
824  for (int i_module = 1; i_module <= iHists.nModules_[i_layer - 1]; i_module++) {
825  unsigned int i_index = i_module + (i_layer - 1) * iHists.nModules_[i_layer - 1];
826 
827  name = fmt::sprintf("h_drift_depth_adc_layer%i_module%i", i_layer, i_module);
828  title = fmt::sprintf(
829  "depth vs drift (ADC) layer%i module%i; drift [#mum]; production depth [#mum]", i_layer, i_module);
830  iHists.h_drift_depth_adc_[i_index] =
831  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
832 
833  name = fmt::sprintf("h_drift_depth_adc2_layer%i_module%i", i_layer, i_module);
834  title = fmt::sprintf(
835  "depth vs drift (ADC^{2}) layer%i module%i; drift [#mum]; production depth [#mum]", i_layer, i_module);
836  iHists.h_drift_depth_adc2_[i_index] =
837  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
838 
839  name = fmt::sprintf("h_drift_depth_noadc_layer%i_module%i", i_layer, i_module);
840  title = fmt::sprintf(
841  "depth vs drift (no ADC) layer%i module%i; drift [#mum]; production depth [#mum]", i_layer, i_module);
842  iHists.h_drift_depth_noadc_[i_index] =
843  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
844 
845  name = fmt::sprintf("h_drift_depth_layer%i_module%i", i_layer, i_module);
846  title =
847  fmt::sprintf("depth vs drift layer%i module%i; drift [#mum]; production depth [#mum]", i_layer, i_module);
848  iHists.h_drift_depth_[i_index] =
849  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
850  }
851  }
852 
853  // book the "new" modules
854  iBooker.setCurrentFolder(fmt::sprintf("%s/BPix/NewModules", folder_.data()));
855  for (int i = 0; i < (int)iHists.BPixnewDetIds_.size(); i++) {
856  int new_index = iHists.nModules_[iHists.nlay - 1] + (iHists.nlay - 1) * iHists.nModules_[iHists.nlay - 1] + 1 + i;
857 
858  name = fmt::sprintf("h_BPixnew_drift_depth_adc_%s", iHists.BPixnewmodulename_[i].c_str());
859  title = fmt::sprintf("depth vs drift (ADC) %s; drift [#mum]; production depth [#mum]",
860  iHists.BPixnewmodulename_[i].c_str());
861  iHists.h_drift_depth_adc_[new_index] =
862  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
863 
864  name = fmt::sprintf("h_BPixnew_drift_depth_adc2_%s", iHists.BPixnewmodulename_[i].c_str());
865  title = fmt::sprintf("depth vs drift (ADC^{2}) %s; drift [#mum]; production depth [#mum]",
866  iHists.BPixnewmodulename_[i].c_str());
867  iHists.h_drift_depth_adc2_[new_index] =
868  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
869 
870  name = fmt::sprintf("h_BPixnew_drift_depth_noadc_%s", iHists.BPixnewmodulename_[i].c_str());
871  title = fmt::sprintf("depth vs drift (no ADC)%s; drift [#mum]; production depth [#mum]",
872  iHists.BPixnewmodulename_[i].c_str());
873  iHists.h_drift_depth_noadc_[new_index] =
874  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
875 
876  name = fmt::sprintf("h_BPixnew_drift_depth_%s", iHists.BPixnewmodulename_[i].c_str());
877  title = fmt::sprintf("depth vs drift %s; drift [#mum]; production depth [#mum]",
878  iHists.BPixnewmodulename_[i].c_str());
879  iHists.h_drift_depth_[new_index] =
880  iBooker.book2D(name, title, hist_drift_, min_drift_, max_drift_, hist_depth_, min_depth_, max_depth_);
881  }
882  } // end if GrazinAngleAnalysis
883  else {
884  iBooker.setCurrentFolder(folder_);
885  std::string baseName;
886  std::string baseTitle;
887 
888  for (int r = 0; r < iHists.nRings_; ++r) {
889  for (int p = 0; p < iHists.nPanels_; ++p) {
890  for (int s = 0; s < iHists.nSides_; ++s) {
891  baseName = fmt::sprintf("R%d_P%d_z%d", r + 1, p + 1, s + 1);
892  if (s == 0)
893  baseTitle = fmt::sprintf("Ring%d_Panel%d_z-", r + 1, p + 1);
894  else
895  baseTitle = fmt::sprintf("Ring%d_Panel%d_z+", r + 1, p + 1);
896 
897  int idx = iHists.nSides_ * iHists.nPanels_ * r + iHists.nSides_ * p + s;
898  int idxBeta = iHists.betaStartIdx_ + idx;
899 
900  name = fmt::sprintf("%s_alphaMean", baseName);
901  title = fmt::sprintf("%s_alphaMean;cot(#alpha); Average cluster size x (pixel)", baseTitle);
902  iHists.h_fpixMean_[idx] = iBooker.book1D(name, title, 60, -3., 3.);
903  name = fmt::sprintf("%s_betaMean", baseName);
904  title = fmt::sprintf("%s_betaMean;cot(#beta); Average cluster size y (pixel)", baseTitle);
905  iHists.h_fpixMean_[idxBeta] = iBooker.book1D(name, title, 60, -3., 3.);
906 
907  } // loop over sides
908  } // loop over panels
909  } // loop over rings
910  iBooker.setCurrentFolder(fmt::sprintf("%s/FPix", folder_.data()));
911  for (int r = 0; r < iHists.nRings_; ++r) {
912  for (int p = 0; p < iHists.nPanels_; ++p) {
913  for (int s = 0; s < iHists.nSides_; ++s) {
914  baseName = fmt::sprintf("R%d_P%d_z%d", r + 1, p + 1, s + 1);
915  if (s == 0)
916  baseTitle = fmt::sprintf("Ring%d_Panel%d_z-", r + 1, p + 1);
917  else
918  baseTitle = fmt::sprintf("Ring%d_Panel%d_z+", r + 1, p + 1);
919 
920  int idx = iHists.nSides_ * iHists.nPanels_ * r + iHists.nSides_ * p + s;
921  int idxBeta = iHists.betaStartIdx_ + idx;
922 
923  name = fmt::sprintf("%s_alpha", baseName);
924  title = fmt::sprintf("%s_alpha;cot(#alpha); Cluster size x (pixel)", baseTitle);
925  iHists.h_fpixAngleSize_[idx] = iBooker.book2D(name, title, 60, -3., 3., 10, 0.5, 10.5);
926  name = fmt::sprintf("%s_beta", baseName);
927  title = fmt::sprintf("%s_beta;cot(#beta); Cluster size y (pixel) ", baseTitle);
928  iHists.h_fpixAngleSize_[idxBeta] = iBooker.book2D(name, title, 60, -3., 3., 10, 0.5, 10.5);
929  for (int m = 0; m < 3; ++m) {
930  name = fmt::sprintf("%s_B%d", baseName, m);
931  char bComp = m == 0 ? 'x' : (m == 1 ? 'y' : 'z');
932  title = fmt::sprintf("%s_magField%d;B_{%c} [T];Entries", baseTitle, m, bComp);
933  iHists.h_fpixMagField_[m][idx] = iBooker.book1D(name, title, 10000, -5., 5.);
934  } // mag. field comps
935  } // loop over sides
936  } // loop over panels
937  } // loop over rings
938  } // if MinimalClusterSize
939 
940  // book the track monitoring plots
941  iBooker.setCurrentFolder(fmt::sprintf("%s/TrackMonitoring", folder_.data()));
942  iHists.h_tracks_ = iBooker.book1D("h_tracks", ";tracker volume;tracks", 2, -0.5, 1.5);
943  iHists.h_tracks_->setBinLabel(1, "all tracks", 1);
944  iHists.h_tracks_->setBinLabel(2, "has pixel hits", 1);
945  iHists.h_trackEta_ = iBooker.book1D("h_trackEta", ";track #eta; #tracks", 30, -3., 3.);
946  iHists.h_trackPhi_ = iBooker.book1D("h_trackPhi", ";track #phi; #tracks", 48, -M_PI, M_PI);
947  iHists.h_trackPt_ = iBooker.book1D("h_trackPt", ";track p_{T} [GeV]; #tracks", 100, 0., 100.);
948  iHists.h_trackChi2_ = iBooker.book1D("h_trackChi2ndof", ";track #chi^{2}/ndof; #tracks", 100, 0., 10.);
949 }
mps_fire.i
i
Definition: mps_fire.py:429
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
visualization-live-secondInstance_cfg.m
m
Definition: visualization-live-secondInstance_cfg.py:80
dqm::impl::MonitorElement::setBinLabel
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:772
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
SiPixelLorentzAnglePCLWorker::iHists
SiPixelLorentzAngleCalibrationHistograms iHists
Definition: SiPixelLorentzAnglePCLWorker.cc:126
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:212
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
SiPixelLorentzAnglePCLWorker::analysisType_
LorentzAngleAnalysisTypeEnum analysisType_
Definition: SiPixelLorentzAnglePCLWorker.cc:133
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:241

◆ convertStringToLorentzAngleAnalysisTypeEnum()

LorentzAngleAnalysisTypeEnum SiPixelLorentzAnglePCLWorker::convertStringToLorentzAngleAnalysisTypeEnum ( std::string  type)
private

Definition at line 995 of file SiPixelLorentzAnglePCLWorker.cc.

References eGrazingAngle, and eMinimumClusterSize.

996  {
997  return (type == "GrazingAngle") ? eGrazingAngle : eMinimumClusterSize;
998 }
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39

◆ dqmBeginRun()

void SiPixelLorentzAnglePCLWorker::dqmBeginRun ( edm::Run const &  run,
edm::EventSetup const &  iSetup 
)
overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 724 of file SiPixelLorentzAnglePCLWorker.cc.

References PixelEndcapName::bladeName(), SiPixelLorentzAngleCalibrationHistograms::BPixnewDetIds_, SiPixelLorentzAngleCalibrationHistograms::BPixnewLayer_, SiPixelLorentzAngleCalibrationHistograms::BPixnewModule_, SiPixelLorentzAngleCalibrationHistograms::BPixnewmodulename_, edm::ESWatcher< T >::check(), hcalRecHitTable_cff::detId, PixelEndcapName::diskName(), SiPixelLorentzAngleCalibrationHistograms::FPixnewBlade_, SiPixelLorentzAngleCalibrationHistograms::FPixnewDetIds_, SiPixelLorentzAngleCalibrationHistograms::FPixnewDisk_, SiPixelLorentzAngleCalibrationHistograms::FPixnewmodulename_, relativeConstraints::geom, geomEsToken_, edm::EventSetup::getData(), PixelBarrelName::getDetId(), PixelEndcapName::getDetId(), mps_fire::i, iHists, PixelBarrelName::layerName(), genParticles_cff::map, PixelBarrelName::moduleName(), newmodulelist_, SiPixelLorentzAngleCalibrationHistograms::nlay, SiPixelLorentzAngleCalibrationHistograms::nModules_, notInPCL_, PixelSubdetector::PixelBarrel, SiPixelTemplate::pushfile(), siPixelTemplateEsToken_, templateDBobject_, thePixelTemp_, topoEsToken_, and watchSiPixelTemplateRcd_.

724  {
725  // geometry
726  const TrackerGeometry* geom = &iSetup.getData(geomEsToken_);
727  const TrackerTopology* tTopo = &iSetup.getData(topoEsToken_);
728 
729  if (notInPCL_) {
730  // Initialize 1D templates
731  if (watchSiPixelTemplateRcd_.check(iSetup)) {
732  templateDBobject_ = &iSetup.getData(siPixelTemplateEsToken_);
733  if (!SiPixelTemplate::pushfile(*templateDBobject_, thePixelTemp_)) {
734  edm::LogError("SiPixelLorentzAnglePCLWorker")
735  << "Templates not filled correctly. Check the sqlite file. Using SiPixelTemplateDBObject version "
736  << (*templateDBobject_).version();
737  }
738  }
739  }
740 
741  PixelTopologyMap map = PixelTopologyMap(geom, tTopo);
742  iHists.nlay = geom->numberOfLayers(PixelSubdetector::PixelBarrel);
743  iHists.nModules_.resize(iHists.nlay);
744  for (int i = 0; i < iHists.nlay; i++) {
745  iHists.nModules_[i] = map.getPXBModules(i + 1);
746  }
747 
748  // list of modules already filled, then return (we already entered here)
749  if (!iHists.BPixnewDetIds_.empty() || !iHists.FPixnewDetIds_.empty())
750  return;
751 
752  if (!newmodulelist_.empty()) {
753  for (auto const& modulename : newmodulelist_) {
754  if (modulename.find("BPix_") != std::string::npos) {
755  PixelBarrelName bn(modulename, true);
756  const auto& detId = bn.getDetId(tTopo);
757  iHists.BPixnewmodulename_.push_back(modulename);
758  iHists.BPixnewDetIds_.push_back(detId.rawId());
759  iHists.BPixnewModule_.push_back(bn.moduleName());
760  iHists.BPixnewLayer_.push_back(bn.layerName());
761  } else if (modulename.find("FPix_") != std::string::npos) {
762  PixelEndcapName en(modulename, true);
763  const auto& detId = en.getDetId(tTopo);
764  iHists.FPixnewmodulename_.push_back(modulename);
765  iHists.FPixnewDetIds_.push_back(detId.rawId());
766  iHists.FPixnewDisk_.push_back(en.diskName());
767  iHists.FPixnewBlade_.push_back(en.bladeName());
768  }
769  }
770  }
771 }
mps_fire.i
i
Definition: mps_fire.py:429
SiPixelLorentzAnglePCLWorker::thePixelTemp_
std::vector< SiPixelTemplateStore > thePixelTemp_
Definition: SiPixelLorentzAnglePCLWorker.cc:131
SiPixelLorentzAnglePCLWorker::watchSiPixelTemplateRcd_
edm::ESWatcher< SiPixelTemplateDBObjectESProducerRcd > watchSiPixelTemplateRcd_
Definition: SiPixelLorentzAnglePCLWorker.cc:129
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
SiPixelLorentzAnglePCLWorker::templateDBobject_
const SiPixelTemplateDBObject * templateDBobject_
Definition: SiPixelLorentzAnglePCLWorker.cc:130
SiPixelLorentzAnglePCLWorker::topoEsToken_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:196
SiPixelLorentzAnglePCLWorker::newmodulelist_
std::vector< std::string > newmodulelist_
Definition: SiPixelLorentzAnglePCLWorker.cc:137
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
genParticles_cff.map
map
Definition: genParticles_cff.py:11
SiPixelLorentzAnglePCLWorker::iHists
SiPixelLorentzAngleCalibrationHistograms iHists
Definition: SiPixelLorentzAnglePCLWorker.cc:126
SiPixelLorentzAnglePCLWorker::geomEsToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:195
edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:57
SiPixelTemplate::pushfile
static bool pushfile(int filenum, std::vector< SiPixelTemplateStore > &pixelTemp, std::string dir="CalibTracker/SiPixelESProducers/data/")
Definition: SiPixelTemplate.cc:130
SiPixelLorentzAnglePCLWorker::siPixelTemplateEsToken_
edm::ESGetToken< SiPixelTemplateDBObject, SiPixelTemplateDBObjectESProducerRcd > siPixelTemplateEsToken_
Definition: SiPixelLorentzAnglePCLWorker.cc:197

◆ dqmEndRun()

void SiPixelLorentzAnglePCLWorker::dqmEndRun ( edm::Run const &  run,
edm::EventSetup const &  iSetup 
)
private

Definition at line 951 of file SiPixelLorentzAnglePCLWorker.cc.

References hFile_, and notInPCL_.

951  {
952  if (notInPCL_) {
953  hFile_->cd();
954  hFile_->Write();
955  hFile_->Close();
956  }
957 }

◆ fillDescriptions()

void SiPixelLorentzAnglePCLWorker::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 1001 of file SiPixelLorentzAnglePCLWorker.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

1001  {
1002  edm::ParameterSetDescription desc;
1003  desc.setComment("Worker module of the SiPixel Lorentz Angle PCL monitoring workflow");
1004  desc.add<std::string>("analysisType", "GrazingAngle")
1005  ->setComment("analysis type - GrazingAngle (default) or MinimumClusterSize");
1006  desc.add<std::string>("folder", "AlCaReco/SiPixelLorentzAngle")->setComment("directory of PCL Worker output");
1007  desc.add<bool>("notInPCL", false)->setComment("create TTree (true) or not (false)");
1008  desc.add<std::string>("fileName", "testrun.root")->setComment("name of the TTree file if notInPCL = true");
1009  desc.add<std::vector<std::string>>("newmodulelist", {})->setComment("the list of DetIds for new sensors");
1010  desc.add<edm::InputTag>("src", edm::InputTag("TrackRefitter"))->setComment("input track collections");
1011  desc.add<double>("ptMin", 3.)->setComment("minimum pt on tracks");
1012  desc.add<double>("normChi2Max", 2.)->setComment("maximum reduced chi squared");
1013  desc.add<std::vector<int>>("clustSizeYMin", {4, 3, 3, 2})
1014  ->setComment("minimum cluster size on Y axis for all Barrel Layers");
1015  desc.add<int>("clustSizeXMax", 5)->setComment("maximum cluster size on X axis");
1016  desc.add<double>("residualMax", 0.005)->setComment("maximum residual");
1017  desc.add<double>("clustChargeMaxPerLength", 50000)
1018  ->setComment("maximum cluster charge per unit length of pixel depth (z)");
1019  desc.add<int>("binsDepth", 50)->setComment("# bins for electron production depth axis");
1020  desc.add<int>("binsDrift", 100)->setComment("# bins for electron drift axis");
1021  descriptions.addWithDefaultLabel(desc);
1022 }
edm::ConfigurationDescriptions::addWithDefaultLabel
void addWithDefaultLabel(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:87
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256

◆ fillPix()

const Pixinfo SiPixelLorentzAnglePCLWorker::fillPix ( const SiPixelCluster &  LocPix,
const PixelTopology *  topol 
) const
private

Definition at line 960 of file SiPixelLorentzAnglePCLWorker.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().

960  {
961  Pixinfo pixinfo;
962  const std::vector<SiPixelCluster::Pixel>& pixvector = LocPix.pixels();
963  pixinfo.npix = 0;
964  for (std::vector<SiPixelCluster::Pixel>::const_iterator itPix = pixvector.begin(); itPix != pixvector.end();
965  itPix++) {
966  pixinfo.row[pixinfo.npix] = itPix->x;
967  pixinfo.col[pixinfo.npix] = itPix->y;
968  pixinfo.adc[pixinfo.npix] = itPix->adc;
969  LocalPoint lp = topol->localPosition(MeasurementPoint(itPix->x + 0.5, itPix->y + 0.5));
970  pixinfo.x[pixinfo.npix] = lp.x();
971  pixinfo.y[pixinfo.npix] = lp.y();
972  pixinfo.npix++;
973  }
974  return pixinfo;
975 }
Pixinfo::x
float x[maxpix]
Definition: SiPixelLorentzAngle.h:54
Topology::localPosition
virtual LocalPoint localPosition(const MeasurementPoint &) const =0
Pixinfo::y
float y[maxpix]
Definition: SiPixelLorentzAngle.h:55
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
MeasurementPoint
Measurement2DPoint MeasurementPoint
Measurement points are two-dimensional by default.
Definition: MeasurementPoint.h:12
SiPixelCluster::pixels
const std::vector< Pixel > pixels() const
Definition: SiPixelCluster.h:160
Pixinfo::adc
float adc[maxpix]
Definition: SiPixelLorentzAngle.h:53
Pixinfo::col
float col[maxpix]
Definition: SiPixelLorentzAngle.h:52
Pixinfo::row
float row[maxpix]
Definition: SiPixelLorentzAngle.h:51

◆ surface_deformation()

const std::pair< LocalPoint, LocalPoint > SiPixelLorentzAnglePCLWorker::surface_deformation ( const PixelTopology *  topol,
TrajectoryStateOnSurface &  tsos,
const SiPixelRecHit *  recHitPix 
) const
private

Definition at line 978 of file SiPixelLorentzAnglePCLWorker.cc.

References LocalTrajectoryParameters::dxdz(), LocalTrajectoryParameters::dydz(), TrajectoryStateOnSurface::localParameters(), BaseTrackerRecHit::localPosition(), TrajectoryStateOnSurface::localPosition(), Topology::localPosition(), and PixelTopology::pixel().

Referenced by analyze().

979  {
980  LocalPoint trackposition = tsos.localPosition();
981  const LocalTrajectoryParameters& ltp = tsos.localParameters();
982  const Topology::LocalTrackAngles localTrackAngles(ltp.dxdz(), ltp.dydz());
983 
984  std::pair<float, float> pixels_track = topol->pixel(trackposition, localTrackAngles);
985  std::pair<float, float> pixels_rechit = topol->pixel(recHitPix->localPosition(), localTrackAngles);
986 
987  LocalPoint lp_track = topol->localPosition(MeasurementPoint(pixels_track.first, pixels_track.second));
988 
989  LocalPoint lp_rechit = topol->localPosition(MeasurementPoint(pixels_rechit.first, pixels_rechit.second));
990 
991  std::pair<LocalPoint, LocalPoint> lps = std::make_pair(lp_track, lp_rechit);
992  return lps;
993 }
PixelTopology::pixel
virtual std::pair< float, float > pixel(const LocalPoint &p) const =0
Topology::localPosition
virtual LocalPoint localPosition(const MeasurementPoint &) const =0
TrajectoryStateOnSurface::localParameters
const LocalTrajectoryParameters & localParameters() const
Definition: TrajectoryStateOnSurface.h:73
MeasurementPoint
Measurement2DPoint MeasurementPoint
Measurement points are two-dimensional by default.
Definition: MeasurementPoint.h:12
BaseTrackerRecHit::localPosition
LocalPoint localPosition() const override
Definition: BaseTrackerRecHit.h:56

Member Data Documentation

◆ analysisType_

LorentzAngleAnalysisTypeEnum SiPixelLorentzAnglePCLWorker::analysisType_
private

Definition at line 133 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and bookHistograms().

◆ bladeF_

int SiPixelLorentzAnglePCLWorker::bladeF_
private

Definition at line 167 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ bx_

int SiPixelLorentzAnglePCLWorker::bx_
private

Definition at line 143 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ chi2_

double SiPixelLorentzAnglePCLWorker::chi2_
private

Definition at line 152 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ clust_

Clust SiPixelLorentzAnglePCLWorker::clust_
private

Definition at line 156 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ clustChargeMaxPerLength_

double SiPixelLorentzAnglePCLWorker::clustChargeMaxPerLength_
private

Definition at line 186 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ clustF_

Clust SiPixelLorentzAnglePCLWorker::clustF_
private

Definition at line 172 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ clustSizeXMax_

int SiPixelLorentzAnglePCLWorker::clustSizeXMax_
private

Definition at line 184 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ clustSizeYMin_

std::vector<int> SiPixelLorentzAnglePCLWorker::clustSizeYMin_
private

Definition at line 183 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ diskF_

int SiPixelLorentzAnglePCLWorker::diskF_
private

Definition at line 166 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ eta_

float SiPixelLorentzAnglePCLWorker::eta_
private

Definition at line 150 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ event_

long int SiPixelLorentzAnglePCLWorker::event_
private

Definition at line 141 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ filename_

std::string SiPixelLorentzAnglePCLWorker::filename_
private

Definition at line 136 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by SiPixelLorentzAnglePCLWorker().

◆ folder_

std::string SiPixelLorentzAnglePCLWorker::folder_
private

Definition at line 134 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by bookHistograms().

◆ geomEsToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiPixelLorentzAnglePCLWorker::geomEsToken_
private

Definition at line 195 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by dqmBeginRun().

◆ geomPerEventEsToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiPixelLorentzAnglePCLWorker::geomPerEventEsToken_
private

Definition at line 199 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ hFile_

std::unique_ptr<TFile> SiPixelLorentzAnglePCLWorker::hFile_
private

Definition at line 190 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by dqmEndRun(), and SiPixelLorentzAnglePCLWorker().

◆ hist_depth_

int SiPixelLorentzAnglePCLWorker::hist_depth_
private

Definition at line 187 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by bookHistograms().

◆ hist_drift_

int SiPixelLorentzAnglePCLWorker::hist_drift_
private

Definition at line 188 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by bookHistograms().

◆ iHists

SiPixelLorentzAngleCalibrationHistograms SiPixelLorentzAnglePCLWorker::iHists
private

Definition at line 126 of file SiPixelLorentzAnglePCLWorker.cc.

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

◆ isflipped_

int SiPixelLorentzAnglePCLWorker::isflipped_
private

Definition at line 148 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ ladder_

int SiPixelLorentzAnglePCLWorker::ladder_
private

Definition at line 146 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ layer_

int SiPixelLorentzAnglePCLWorker::layer_
private

Definition at line 147 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ lumiblock_

int SiPixelLorentzAnglePCLWorker::lumiblock_
private

Definition at line 142 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ magneticFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SiPixelLorentzAnglePCLWorker::magneticFieldToken_
private

Definition at line 200 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ module_

int SiPixelLorentzAnglePCLWorker::module_
private

Definition at line 145 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ moduleF_

int SiPixelLorentzAnglePCLWorker::moduleF_
private

Definition at line 169 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ ndof_

double SiPixelLorentzAnglePCLWorker::ndof_
private

Definition at line 153 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ newmodulelist_

std::vector<std::string> SiPixelLorentzAnglePCLWorker::newmodulelist_
private

Definition at line 137 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by dqmBeginRun().

◆ normChi2Max_

double SiPixelLorentzAnglePCLWorker::normChi2Max_
private

Definition at line 182 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ notInPCL_

bool SiPixelLorentzAnglePCLWorker::notInPCL_
private

Definition at line 135 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), dqmBeginRun(), dqmEndRun(), and SiPixelLorentzAnglePCLWorker().

◆ orbit_

int SiPixelLorentzAnglePCLWorker::orbit_
private

Definition at line 144 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ panelF_

int SiPixelLorentzAnglePCLWorker::panelF_
private

Definition at line 168 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ phi_

float SiPixelLorentzAnglePCLWorker::phi_
private

Definition at line 151 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ pixinfo_

Pixinfo SiPixelLorentzAnglePCLWorker::pixinfo_
private

Definition at line 154 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ pixinfoF_

Pixinfo SiPixelLorentzAnglePCLWorker::pixinfoF_
private

Definition at line 170 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ pt_

float SiPixelLorentzAnglePCLWorker::pt_
private

Definition at line 149 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ ptmin_

double SiPixelLorentzAnglePCLWorker::ptmin_
private

Definition at line 181 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ qScale_

float SiPixelLorentzAnglePCLWorker::qScale_
private

Definition at line 161 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ qScaleF_

float SiPixelLorentzAnglePCLWorker::qScaleF_
private

Definition at line 177 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ rechit_

Rechit SiPixelLorentzAnglePCLWorker::rechit_
private

Definition at line 157 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ rechitCorr_

Rechit SiPixelLorentzAnglePCLWorker::rechitCorr_
private

Definition at line 158 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ rechitCorrF_

Rechit SiPixelLorentzAnglePCLWorker::rechitCorrF_
private

Definition at line 174 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ rechitF_

Rechit SiPixelLorentzAnglePCLWorker::rechitF_
private

Definition at line 173 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ residualMax_

double SiPixelLorentzAnglePCLWorker::residualMax_
private

Definition at line 185 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ rQmQt_

float SiPixelLorentzAnglePCLWorker::rQmQt_
private

Definition at line 162 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ rQmQtF_

float SiPixelLorentzAnglePCLWorker::rQmQtF_
private

Definition at line 178 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ run_

int SiPixelLorentzAnglePCLWorker::run_
private

Definition at line 140 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ sideF_

int SiPixelLorentzAnglePCLWorker::sideF_
private

Definition at line 165 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ simhit_

Hit SiPixelLorentzAnglePCLWorker::simhit_
private

Definition at line 155 of file SiPixelLorentzAnglePCLWorker.cc.

◆ simhitF_

Hit SiPixelLorentzAnglePCLWorker::simhitF_
private

Definition at line 171 of file SiPixelLorentzAnglePCLWorker.cc.

◆ SiPixelLorentzAngleTreeBarrel_

std::unique_ptr<TTree> SiPixelLorentzAnglePCLWorker::SiPixelLorentzAngleTreeBarrel_
private

Definition at line 191 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ SiPixelLorentzAngleTreeForward_

std::unique_ptr<TTree> SiPixelLorentzAnglePCLWorker::SiPixelLorentzAngleTreeForward_
private

Definition at line 192 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ siPixelTemplateEsToken_

edm::ESGetToken<SiPixelTemplateDBObject, SiPixelTemplateDBObjectESProducerRcd> SiPixelLorentzAnglePCLWorker::siPixelTemplateEsToken_
private

Definition at line 197 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by dqmBeginRun().

◆ t_trajTrack

edm::EDGetTokenT<TrajTrackAssociationCollection> SiPixelLorentzAnglePCLWorker::t_trajTrack
private

Definition at line 203 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ templateDBobject_

const SiPixelTemplateDBObject* SiPixelLorentzAnglePCLWorker::templateDBobject_
private

Definition at line 130 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and dqmBeginRun().

◆ thePixelTemp_

std::vector<SiPixelTemplateStore> SiPixelLorentzAnglePCLWorker::thePixelTemp_
private

Definition at line 131 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and dqmBeginRun().

◆ topoEsToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiPixelLorentzAnglePCLWorker::topoEsToken_
private

Definition at line 196 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by dqmBeginRun().

◆ topoPerEventEsToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiPixelLorentzAnglePCLWorker::topoPerEventEsToken_
private

Definition at line 198 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze().

◆ trackhit_

Hit SiPixelLorentzAnglePCLWorker::trackhit_
private

Definition at line 155 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ trackhitCorrX_

float SiPixelLorentzAnglePCLWorker::trackhitCorrX_
private

Definition at line 159 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ trackhitCorrXF_

float SiPixelLorentzAnglePCLWorker::trackhitCorrXF_
private

Definition at line 175 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ trackhitCorrY_

float SiPixelLorentzAnglePCLWorker::trackhitCorrY_
private

Definition at line 160 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ trackhitCorrYF_

float SiPixelLorentzAnglePCLWorker::trackhitCorrYF_
private

Definition at line 176 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ trackhitF_

Hit SiPixelLorentzAnglePCLWorker::trackhitF_
private

Definition at line 171 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by analyze(), and SiPixelLorentzAnglePCLWorker().

◆ watchSiPixelTemplateRcd_

edm::ESWatcher<SiPixelTemplateDBObjectESProducerRcd> SiPixelLorentzAnglePCLWorker::watchSiPixelTemplateRcd_
private

Definition at line 129 of file SiPixelLorentzAnglePCLWorker.cc.

Referenced by dqmBeginRun().

Generated for CMSSW Reference Manual by doxygen 1.8.14