SiPixelVCalReader Class Reference

Inheritance diagram for SiPixelVCalReader:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	SiPixelVCalReader (const edm::ParameterSet &)
	~SiPixelVCalReader () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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
bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
TH1F *	offsetBPix_
TH1F *	offsetFPix_
const bool	printdebug_
const edm::ESGetToken< SiPixelVCal, SiPixelVCalSimRcd >	siPixelVCalSimToken_
const edm::ESGetToken< SiPixelVCal, SiPixelVCalRcd >	siPixelVCalToken_
TH1F *	slopeBPix_
TH1F *	slopeFPix_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tkTopoToken_
const bool	useSimRcd_

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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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 36 of file SiPixelVCalReader.cc.

Constructor & Destructor Documentation

SiPixelVCalReader()

SiPixelVCalReader::SiPixelVCalReader ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 60 of file SiPixelVCalReader.cc.

References TFileService::kSharedResource.

    : siPixelVCalSimToken_(esConsumes()),
      siPixelVCalToken_(esConsumes()),
      tkGeomToken_(esConsumes()),
      tkTopoToken_(esConsumes()),
      printdebug_(iConfig.getUntrackedParameter<bool>("printDebug", false)),
      useSimRcd_(iConfig.getParameter<bool>("useSimRcd")) {
  usesResource(TFileService::kSharedResource);
}

~SiPixelVCalReader()

SiPixelVCalReader::~SiPixelVCalReader ( )

overridedefault

Member Function Documentation

analyze()

void SiPixelVCalReader::analyze ( const edm::Event &  e, const edm::EventSetup &  iSetup  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 72 of file SiPixelVCalReader.cc.

References TrackerGeometry::detUnits(), compareTotals::fs, edm::EventSetup::getData(), SiPixelVCal::getSlopeAndOffset(), PVValHelper::ladder, phase1PixelTopology::layer, hltrates_dqm_sourceclient-live_cfg::offset, offsetBPix_, offsetFPix_, PVValHelper::phase1, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), relativeConstraints::ring, PixelEndcapName::ringName(), siPixelVCalSimToken_, siPixelVCalToken_, slope, slopeBPix_, slopeFPix_, DetId::subdetId(), tkGeomToken_, tkTopoToken_, and useSimRcd_.

                                                                              {
  const SiPixelVCal* siPixelVCal;

  // Get record & file service
  if (useSimRcd_) {
    siPixelVCal = &iSetup.getData(siPixelVCalSimToken_);
  } else {
    siPixelVCal = &iSetup.getData(siPixelVCalToken_);
  }
  edm::LogInfo("SiPixelVCalReader") << "[SiPixelVCalReader::analyze] End Reading SiPixelVCal" << std::endl;
  edm::Service<TFileService> fs;

  // Retrieve tracker topology from geometry
  const TrackerTopology* const tTopo = &iSetup.getData(tkTopoToken_);

  // Retrieve old style tracker geometry from geometry
  const TrackerGeometry* pDD = &iSetup.getData(tkGeomToken_);
  edm::LogPrint("SiPixelVCalReader") << " There are " << pDD->detUnits().size() << " modules" << std::endl;

  // Phase
  bool phase1 = true;

  // Prepare tree
  TTree* tree = new TTree("tree", "tree");
  uint32_t detid, subdet, layer, ladder, side, disk, ring;
  double slope, offset;
  tree->Branch("detid", &detid, "detid/I");
  tree->Branch("subdet", &subdet, "subdet/I");
  tree->Branch("layer", &layer, "layer/I");
  tree->Branch("ladder", &ladder, "ladder/I");
  tree->Branch("side", &side, "side/I");
  tree->Branch("disk", &disk, "disk/I");
  tree->Branch("ring", &ring, "ring/I");
  tree->Branch("slope", &slope, "slope/D");
  tree->Branch("offset", &offset, "offset/D");

  // Prepare histograms
  slopeBPix_ = fs->make<TH1F>("VCalSlopeBarrelPixel", "VCalSlopeBarrelPixel", 150, 0, 100);
  slopeFPix_ = fs->make<TH1F>("VCalSlopeForwardPixel", "VCalSlopeForwardPixel", 150, 0, 100);
  offsetBPix_ = fs->make<TH1F>("VCalOffsetBarrelPixel", "VCalOffsetBarrelPixel", 200, -900, 100);
  offsetFPix_ = fs->make<TH1F>("VCalOffsetForwardPixel", "VCalOffsetForwardPixel", 200, -900, 100);
  std::map<unsigned int, SiPixelVCal::VCal> vcal = siPixelVCal->getSlopeAndOffset();
  std::map<unsigned int, SiPixelVCal::VCal>::const_iterator it;

  // Fill histograms
  edm::LogPrint("SiPixelVCalReader") << std::setw(12) << "detid" << std::setw(8) << "subdet" << std::setw(8) << "layer"
                                     << std::setw(8) << "disk" << std::setw(14) << "VCal slope" << std::setw(8)
                                     << "offset" << std::endl;
  for (it = vcal.begin(); it != vcal.end(); it++) {
    detid = it->first;
    slope = it->second.slope;
    offset = it->second.offset;
    const DetId detIdObj(detid);
    PixelEndcapName fpix(detid, tTopo, phase1);
    subdet = detIdObj.subdetId();
    layer = tTopo->pxbLayer(detIdObj);    // 1, 2, 3, 4
    ladder = tTopo->pxbLadder(detIdObj);  // 1-12/28/44/64
    side = tTopo->pxfSide(detIdObj);      // 1, 2
    disk = tTopo->pxfDisk(detIdObj);      // 1, 2, 3
    ring = fpix.ringName();               // 1 (lower), 2 (upper)
    edm::LogPrint("SiPixelVCalReader") << std::setw(12) << detid << std::setw(8) << subdet << std::setw(8) << layer
                                       << std::setw(8) << disk << std::setw(14) << slope << std::setw(8) << offset
                                       << std::endl;
    // edm::LogPrint("SiPixelVCalReader") << "detid " << detid << ", subdet " << subdet << ", layer " <<
    // layer << ", disk " << disk
    //          << ", VCal slope " << slope << ", offset " << offset <<
    //          std::endl;
    // edm::LogInfo("SiPixelVCalReader") << "detid " << detid << ", subdet " <<
    // subdet << ", layer " << layer
    //                                  << ", VCal slope " << slope << ", offset
    //                                  " << offset;
    if (subdet == static_cast<int>(PixelSubdetector::PixelBarrel)) {
      slopeBPix_->Fill(slope);
      offsetBPix_->Fill(offset);
    } else if (subdet == static_cast<int>(PixelSubdetector::PixelEndcap)) {
      slopeFPix_->Fill(slope);
      offsetFPix_->Fill(offset);
    }
    tree->Fill();
  }
}

Member Data Documentation

offsetBPix_

TH1F* SiPixelVCalReader::offsetBPix_

private

Definition at line 54 of file SiPixelVCalReader.cc.

Referenced by analyze().

offsetFPix_

TH1F* SiPixelVCalReader::offsetFPix_

private

Definition at line 55 of file SiPixelVCalReader.cc.

Referenced by analyze().

printdebug_

const bool SiPixelVCalReader::printdebug_

private

Definition at line 49 of file SiPixelVCalReader.cc.

siPixelVCalSimToken_

const edm::ESGetToken<SiPixelVCal, SiPixelVCalSimRcd> SiPixelVCalReader::siPixelVCalSimToken_

private

Definition at line 44 of file SiPixelVCalReader.cc.

Referenced by analyze().

siPixelVCalToken_

const edm::ESGetToken<SiPixelVCal, SiPixelVCalRcd> SiPixelVCalReader::siPixelVCalToken_

private

Definition at line 45 of file SiPixelVCalReader.cc.

Referenced by analyze().

slopeBPix_

TH1F* SiPixelVCalReader::slopeBPix_

private

Definition at line 52 of file SiPixelVCalReader.cc.

Referenced by analyze().

slopeFPix_

TH1F* SiPixelVCalReader::slopeFPix_

private

Definition at line 53 of file SiPixelVCalReader.cc.

Referenced by analyze().

tkGeomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiPixelVCalReader::tkGeomToken_

private

Definition at line 46 of file SiPixelVCalReader.cc.

Referenced by analyze().

tkTopoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiPixelVCalReader::tkTopoToken_

private

Definition at line 47 of file SiPixelVCalReader.cc.

Referenced by analyze().

useSimRcd_

const bool SiPixelVCalReader::useSimRcd_

private

Definition at line 50 of file SiPixelVCalReader.cc.

Referenced by analyze().