TrackerOfflineValidationSummary Class Reference

#include <Alignment/TrackerOfflineValidationSummary/src/TrackerOfflineValidationSummary.cc>

Inheritance diagram for TrackerOfflineValidationSummary:

Classes
struct	HarvestingHierarchy
struct	HarvestingHistos
	Put here the histograms created during harvesting. More...
struct	ModuleHistos

Public Types
typedef dqm::legacy::DQMStore	DQMStore
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (const edm::Event &evt, const edm::EventSetup &) override
void	applyHarvestingHierarchy (TTree &treeMem)
const std::string	associateModuleHistsWithTree (const TkOffTreeVariables &treeMem, TrackerOfflineValidationSummary::ModuleHistos &moduleHists, std::map< std::string, std::string > &substructureName)
void	beginRun (const edm::Run &, const edm::EventSetup &iSetup) override
void	bookHarvestingHists ()
void	collateHarvestingHists (TTree &tree)
void	endJob () override
void	endRun (const edm::Run &, const edm::EventSetup &iSetup) override
void	fillHarvestingHists (TTree &tree)
void	fillTree (TTree &tree, std::map< int, TrackerOfflineValidationSummary::ModuleHistos > &moduleHist, TkOffTreeVariables &treeMem, const TrackerGeometry &tkgeom, std::map< std::string, std::string > &substructureName, const TrackerTopology *tTopo)
std::pair< float, float >	fitResiduals (TH1 *hist) const
void	getBinning (const std::string &binningPSetName, int &nBinsX, double &lowerBoundX, double &upperBoundX) const
float	getMedian (const TH1 *hist) const

Private Attributes
DQMStore *	dbe_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_
const std::string	moduleDirectory_
bool	moduleMapsInitialized
std::map< int, TrackerOfflineValidationSummary::ModuleHistos >	mPxbResiduals_
std::map< int, TrackerOfflineValidationSummary::ModuleHistos >	mPxeResiduals_
std::map< int, TrackerOfflineValidationSummary::ModuleHistos >	mTecResiduals_
std::map< int, TrackerOfflineValidationSummary::ModuleHistos >	mTibResiduals_
std::map< int, TrackerOfflineValidationSummary::ModuleHistos >	mTidResiduals_
std::map< int, TrackerOfflineValidationSummary::ModuleHistos >	mTobResiduals_
const edm::ParameterSet	parSet_
edm::ESHandle< TrackerGeometry >	tkGeom_
std::unique_ptr< TrackerTopology >	tTopo_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_
const bool	useFit_
std::vector< HarvestingHierarchy >	vHarvestingHierarchy_

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

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 58 of file TrackerOfflineValidationSummary.cc.

Member Typedef Documentation

DQMStore

typedef dqm::legacy::DQMStore TrackerOfflineValidationSummary::DQMStore

Definition at line 60 of file TrackerOfflineValidationSummary.cc.

Constructor & Destructor Documentation

TrackerOfflineValidationSummary()

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

explicit

Definition at line 153 of file TrackerOfflineValidationSummary.cc.

References dbe_, and Utilities::operator.

    : geomToken_(esConsumes()),
      tTopoToken_(esConsumes<edm::Transition::EndRun>()),
      parSet_(iConfig),
      tTopo_(nullptr),
      moduleDirectory_(parSet_.getParameter<std::string>("moduleDirectoryInOutput")),
      useFit_(parSet_.getParameter<bool>("useFit")),
      dbe_(nullptr),
      moduleMapsInitialized(false) {
  //now do what ever initialization is needed
  dbe_ = edm::Service<DQMStore>().operator->();
}

~TrackerOfflineValidationSummary()

TrackerOfflineValidationSummary::~TrackerOfflineValidationSummary ( )

overridedefault

Member Function Documentation

analyze()

void TrackerOfflineValidationSummary::analyze ( const edm::Event &  evt, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 193 of file TrackerOfflineValidationSummary.cc.

References hcalRecHitTable_cff::detId, TrackerGeometry::detIds(), Exception, geomToken_, edm::EventSetup::getHandle(), moduleMapsInitialized, mPxbResiduals_, mPxeResiduals_, mTecResiduals_, mTibResiduals_, mTidResiduals_, mTobResiduals_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, nano_mu_digi_cff::rawId, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, tkGeom_, and StripSubdetector::TOB.

                                                                                                 {
  // Access of EventSetup is needed to get the list of silicon-modules and their IDs
  // Since they do not change, it is accessed only once
  if (moduleMapsInitialized)
    return;
  tkGeom_ = iSetup.getHandle(geomToken_);
  const TrackerGeometry* bareTkGeomPtr = &(*tkGeom_);

  const TrackingGeometry::DetIdContainer& detIdContainer = bareTkGeomPtr->detIds();
  std::vector<DetId>::const_iterator iDet;
  for (iDet = detIdContainer.begin(); iDet != detIdContainer.end(); ++iDet) {
    const DetId& detId = *iDet;
    const uint32_t rawId = detId.rawId();
    const unsigned int subdetId = detId.subdetId();
    if (subdetId == PixelSubdetector::PixelBarrel)
      mPxbResiduals_[rawId];
    else if (subdetId == PixelSubdetector::PixelEndcap)
      mPxeResiduals_[rawId];
    else if (subdetId == StripSubdetector::TIB)
      mTibResiduals_[rawId];
    else if (subdetId == StripSubdetector::TID)
      mTidResiduals_[rawId];
    else if (subdetId == StripSubdetector::TOB)
      mTobResiduals_[rawId];
    else if (subdetId == StripSubdetector::TEC)
      mTecResiduals_[rawId];
    else {
      throw cms::Exception("Geometry Error")
          << "[TrackerOfflineValidationSummary] Error, tried to get reference for non-tracker subdet " << subdetId
          << " from detector " << detId.det();
    }
  }
  moduleMapsInitialized = true;
}

applyHarvestingHierarchy()

void TrackerOfflineValidationSummary::applyHarvestingHierarchy ( TTree &  treeMem )

private

Definition at line 747 of file TrackerOfflineValidationSummary.cc.

References bookHarvestingHists(), TkOffTreeVariables::isDoubleSide, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, TkOffTreeVariables::side, AlCaHLTBitMon_QueryRunRegistry::string, TkOffTreeVariables::subDetId, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, and vHarvestingHierarchy_.

Referenced by collateHarvestingHists().

                                                                          {
  TkOffTreeVariables* treeMemPtr = nullptr;
  std::map<std::string, std::string>* substructureName = nullptr;
  tree.SetBranchAddress("TkOffTreeVariables", &treeMemPtr);
  tree.SetBranchAddress("SubstructureName", &substructureName);

  // Loop over modules to select accumulation criteria for harvesting plots
  for (unsigned int iSubdet = 1; iSubdet < 7; ++iSubdet) {
    std::string hierarchyName("");
    std::string componentName("");
    std::vector<unsigned int> treeEntries;
    for (unsigned int iSide = 1; iSide < 3; ++iSide) {
      // Set up only one collection for Barrels, not separated for side
      if (iSide == 1 && (iSubdet == PixelSubdetector::PixelBarrel || iSubdet == StripSubdetector::TIB ||
                         iSubdet == StripSubdetector::TOB))
        continue;
      for (int iTree = 0; iTree < tree.GetEntries(); ++iTree) {
        tree.GetEntry(iTree);
        // Do not use glued Dets
        if (treeMemPtr->isDoubleSide)
          continue;
        if (treeMemPtr->subDetId == iSubdet) {
          if (iSide != treeMemPtr->side && (iSubdet == PixelSubdetector::PixelEndcap ||
                                            iSubdet == StripSubdetector::TID || iSubdet == StripSubdetector::TEC))
            continue;
          treeEntries.push_back(iTree);
          if (hierarchyName.length() == 0) {
            hierarchyName = (*substructureName)["subdet"];
            componentName = (*substructureName)["component"];
          }
        }
      }
      HarvestingHierarchy harvestingHierarchy(hierarchyName, componentName, treeEntries);
      vHarvestingHierarchy_.push_back(harvestingHierarchy);
      hierarchyName = "";
      componentName = "";
      treeEntries.clear();
    }
  }
  // Here could be a further separation of the HarvestingHierarchy.
  // E.g. separate the existing ones by layer and add them to the vector without deleting any element from the vector.
  // The existing hists will stay and the new ones are added

  // Now, book the corresponding histos
  this->bookHarvestingHists();
}

associateModuleHistsWithTree()

const std::string TrackerOfflineValidationSummary::associateModuleHistsWithTree ( const TkOffTreeVariables &  treeMem, TrackerOfflineValidationSummary::ModuleHistos &  moduleHists, std::map< std::string, std::string > &  substructureName  )

private

Definition at line 501 of file TrackerOfflineValidationSummary.cc.

References TkOffTreeVariables::blade, dbe_, contentValuesFiles::fullPath, dqm::implementation::IGetter::get(), dqm::legacy::MonitorElement::getTH1(), TkOffTreeVariables::half, L1TStage2EMTFEventInfoClient_cfi::histDir, TkOffTreeVariables::isDoubleSide, PVValHelper::ladder, nano_mu_digi_cff::layer, TkOffTreeVariables::layer, TkOffTreeVariables::module, moduleDirectory_, TkOffTreeVariables::moduleId, TrackerOfflineValidationSummary::ModuleHistos::NormResHisto, TrackerOfflineValidationSummary::ModuleHistos::NormResXprimeHisto, TrackerOfflineValidationSummary::ModuleHistos::NormResYprimeHisto, TkOffTreeVariables::outerInner, TkOffTreeVariables::panel, TkOffTreeVariables::petal, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, TrackerOfflineValidationSummary::ModuleHistos::ResHisto, TrackerOfflineValidationSummary::ModuleHistos::ResXprimeHisto, TrackerOfflineValidationSummary::ModuleHistos::ResYprimeHisto, TkOffTreeVariables::ring, relativeConstraints::ring, TkOffTreeVariables::rod, TkOffTreeVariables::side, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), AlCaHLTBitMon_QueryRunRegistry::string, TkOffTreeVariables::subDetId, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

Referenced by fillTree().

                                                    {
  std::stringstream histDir, sSubdetName;
  std::string componentName;
  if (moduleDirectory_.length() != 0)
    histDir << moduleDirectory_ << "/";
  std::string wheelOrLayer("_layer_");
  if (treeMem.subDetId == PixelSubdetector::PixelBarrel) {
    unsigned int half(treeMem.half), layer(treeMem.layer), ladder(0);
    if (layer == 1) {
      if (half == 2)
        ladder = treeMem.rod - 5;
      else if (treeMem.rod > 15)
        ladder = treeMem.rod - 10;
      else
        ladder = treeMem.rod;
    } else if (layer == 2) {
      if (half == 2)
        ladder = treeMem.rod - 8;
      else if (treeMem.rod > 24)
        ladder = treeMem.rod - 16;
      else
        ladder = treeMem.rod;
    } else if (layer == 3) {
      if (half == 2)
        ladder = treeMem.rod - 11;
      else if (treeMem.rod > 33)
        ladder = treeMem.rod - 22;
      else
        ladder = treeMem.rod;
    }
    componentName = "Pixel";
    sSubdetName << "TPBBarrel_1";
    histDir << componentName << "/" << sSubdetName.str() << "/TPBHalfBarrel_" << treeMem.half << "/TPBLayer_"
            << treeMem.layer << "/TPBLadder_" << ladder;
  } else if (treeMem.subDetId == PixelSubdetector::PixelEndcap) {
    unsigned int side(treeMem.side), half(treeMem.half), blade(0);
    if (side == 1)
      side = 3;
    if (half == 2)
      blade = treeMem.blade - 6;
    else if (treeMem.blade > 18)
      blade = treeMem.blade - 12;
    else
      blade = treeMem.blade;
    componentName = "Pixel";
    sSubdetName << "TPEEndcap_" << side;
    histDir << componentName << "/" << sSubdetName.str() << "/TPEHalfCylinder_" << treeMem.half << "/TPEHalfDisk_"
            << treeMem.layer << "/TPEBlade_" << blade << "/TPEPanel_" << treeMem.panel;
    wheelOrLayer = "_wheel_";
  } else if (treeMem.subDetId == StripSubdetector::TIB) {
    unsigned int half(treeMem.half), layer(treeMem.layer), surface(treeMem.outerInner), string(0);
    if (half == 2) {
      if (layer == 1) {
        if (surface == 1)
          string = treeMem.rod - 13;
        else if (surface == 2)
          string = treeMem.rod - 15;
      }
      if (layer == 2) {
        if (surface == 1)
          string = treeMem.rod - 17;
        else if (surface == 2)
          string = treeMem.rod - 19;
      }
      if (layer == 3) {
        if (surface == 1)
          string = treeMem.rod - 22;
        else if (surface == 2)
          string = treeMem.rod - 23;
      }
      if (layer == 4) {
        if (surface == 1)
          string = treeMem.rod - 26;
        else if (surface == 2)
          string = treeMem.rod - 28;
      }
    } else
      string = treeMem.rod;
    std::stringstream detString;
    if (treeMem.layer < 3 && !treeMem.isDoubleSide)
      detString << "/Det_" << treeMem.module;
    else
      detString << "";
    componentName = "Strip";
    sSubdetName << "TIBBarrel_1";
    histDir << componentName << "/" << sSubdetName.str() << "/TIBHalfBarrel_" << treeMem.side << "/TIBLayer_"
            << treeMem.layer << "/TIBHalfShell_" << treeMem.half << "/TIBSurface_" << treeMem.outerInner
            << "/TIBString_" << string << detString.str();
  } else if (treeMem.subDetId == StripSubdetector::TID) {
    unsigned int side(treeMem.side), outerInner(0);
    if (side == 1)
      side = 3;
    if (treeMem.outerInner == 1)
      outerInner = 2;
    else if (treeMem.outerInner == 2)
      outerInner = 1;
    std::stringstream detString;
    if (treeMem.ring < 3 && !treeMem.isDoubleSide)
      detString << "/Det_" << treeMem.module;
    else
      detString << "";
    componentName = "Strip";
    sSubdetName << "TIDEndcap_" << side;
    histDir << componentName << "/" << sSubdetName.str() << "/TIDDisk_" << treeMem.layer << "/TIDRing_" << treeMem.ring
            << "/TIDSide_" << outerInner << detString.str();
    wheelOrLayer = "_wheel_";
  } else if (treeMem.subDetId == StripSubdetector::TOB) {
    std::stringstream detString;
    if (treeMem.layer < 3 && !treeMem.isDoubleSide)
      detString << "/Det_" << treeMem.module;
    else
      detString << "";
    componentName = "Strip";
    sSubdetName << "TOBBarrel_4";
    histDir << componentName << "/" << sSubdetName.str() << "/TOBHalfBarrel_" << treeMem.side << "/TOBLayer_"
            << treeMem.layer << "/TOBRod_" << treeMem.rod << detString.str();
  } else if (treeMem.subDetId == StripSubdetector::TEC) {
    unsigned int side(0), outerInner(0), ring(0);
    if (treeMem.side == 1)
      side = 6;
    else if (treeMem.side == 2)
      side = 5;
    if (treeMem.outerInner == 1)
      outerInner = 2;
    else if (treeMem.outerInner == 2)
      outerInner = 1;
    if (treeMem.layer > 3 && treeMem.layer < 7)
      ring = treeMem.ring - 1;
    else if (treeMem.layer == 7 || treeMem.layer == 8)
      ring = treeMem.ring - 2;
    else if (treeMem.layer == 9)
      ring = treeMem.ring - 3;
    else
      ring = treeMem.ring;
    std::stringstream detString;
    if ((treeMem.ring < 3 || treeMem.ring == 5) && !treeMem.isDoubleSide)
      detString << "/Det_" << treeMem.module;
    else
      detString << "";
    componentName = "Strip";
    sSubdetName << "TECEndcap_" << side;
    histDir << componentName << "/" << sSubdetName.str() << "/TECDisk_" << treeMem.layer << "/TECSide_" << outerInner
            << "/TECPetal_" << treeMem.petal << "/TECRing_" << ring << detString.str();
    wheelOrLayer = "_wheel_";
  }

  substructureName["component"] = componentName;
  substructureName["subdet"] = sSubdetName.str();

  std::stringstream histName;
  histName << "residuals_subdet_" << treeMem.subDetId << wheelOrLayer << treeMem.layer << "_module_"
           << treeMem.moduleId;

  std::string fullPath;
  fullPath = histDir.str() + "/h_xprime_" + histName.str();
  if (dbe_->get(fullPath))
    moduleHists.ResXprimeHisto = dbe_->get(fullPath)->getTH1();
  else {
    edm::LogError("TrackerOfflineValidationSummary")
        << "Problem with names in input file produced in TrackerOfflineValidation ...\n"
        << "This histogram should exist in every configuration, "
        << "but no histogram with name " << fullPath << " is found!";
    return "";
  }
  fullPath = histDir.str() + "/h_normxprime" + histName.str();
  if (dbe_->get(fullPath))
    moduleHists.NormResXprimeHisto = dbe_->get(fullPath)->getTH1();
  fullPath = histDir.str() + "/h_yprime_" + histName.str();
  if (dbe_->get(fullPath))
    moduleHists.ResYprimeHisto = dbe_->get(fullPath)->getTH1();
  fullPath = histDir.str() + "/h_normyprime" + histName.str();
  if (dbe_->get(fullPath))
    moduleHists.NormResYprimeHisto = dbe_->get(fullPath)->getTH1();
  fullPath = histDir.str() + "/h_" + histName.str();
  if (dbe_->get(fullPath))
    moduleHists.ResHisto = dbe_->get(fullPath)->getTH1();
  fullPath = histDir.str() + "/h_norm" + histName.str();
  if (dbe_->get(fullPath))
    moduleHists.NormResHisto = dbe_->get(fullPath)->getTH1();

  return histDir.str();
}

beginRun()

void TrackerOfflineValidationSummary::beginRun ( const edm::Run &  , const edm::EventSetup &  iSetup  )

inlineoverrideprivate

Definition at line 101 of file TrackerOfflineValidationSummary.cc.

{};

bookHarvestingHists()

void TrackerOfflineValidationSummary::bookHarvestingHists ( )

private

Definition at line 794 of file TrackerOfflineValidationSummary.cc.

References dqm::implementation::IBooker::book1D(), dbe_, getBinning(), edm::ParameterSet::getParameter(), moduleDirectory_, parSet_, dqm::implementation::DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, vHarvestingHierarchy_, multiplicitycorr_cfi::xMax, and photonAnalyzer_cfi::xMin.

Referenced by applyHarvestingHierarchy().

                                                          {
  edm::LogInfo("TrackerOfflineValidationSummary") << "Harvesting histograms will be booked for "
                                                  << vHarvestingHierarchy_.size() << " different hierarchy selections";
  for (std::vector<HarvestingHierarchy>::iterator iHier = vHarvestingHierarchy_.begin();
       iHier != vHarvestingHierarchy_.end();
       ++iHier) {
    std::stringstream dmrXprimeHistoName, dmrYprimeHistoName, dmrXprimeHistoTitle, dmrYprimeHistoTitle;
    dmrXprimeHistoName << "h_DmrXprime_" << iHier->hierarchyName;
    dmrYprimeHistoName << "h_DmrYprime_" << iHier->hierarchyName;
    dmrXprimeHistoTitle << "DMR for " << iHier->hierarchyName << ";<#DeltaX> [cm];# modules";
    dmrYprimeHistoTitle << "DMR for " << iHier->hierarchyName << ";<#DeltaY> [cm];# modules";

    std::string directoryString(moduleDirectory_);
    if (directoryString.length() != 0)
      directoryString += "/";
    directoryString += iHier->componentName;
    dbe_->setCurrentFolder(directoryString);

    int nBinsX(0);
    double xMin(0.), xMax(0.);
    if (iHier->componentName == "Pixel") {
      this->getBinning("TH1DmrXprimePixelModules", nBinsX, xMin, xMax);
      iHier->harvestingHistos.DmrXprime =
          dbe_->book1D(dmrXprimeHistoName.str(), dmrXprimeHistoTitle.str(), nBinsX, xMin, xMax)->getTH1();
      this->getBinning("TH1DmrYprimePixelModules", nBinsX, xMin, xMax);
      iHier->harvestingHistos.DmrYprime =
          dbe_->book1D(dmrYprimeHistoName.str(), dmrYprimeHistoTitle.str(), nBinsX, xMin, xMax)->getTH1();
    } else if (iHier->componentName == "Strip") {
      this->getBinning("TH1DmrXprimeStripModules", nBinsX, xMin, xMax);
      iHier->harvestingHistos.DmrXprime =
          dbe_->book1D(dmrXprimeHistoName.str(), dmrXprimeHistoTitle.str(), nBinsX, xMin, xMax)->getTH1();
      if (!parSet_.getParameter<bool>("stripYDmrs"))
        continue;
      this->getBinning("TH1DmrYprimeStripModules", nBinsX, xMin, xMax);
      iHier->harvestingHistos.DmrYprime =
          dbe_->book1D(dmrYprimeHistoName.str(), dmrYprimeHistoTitle.str(), nBinsX, xMin, xMax)->getTH1();
    }
  }
}

collateHarvestingHists()

void TrackerOfflineValidationSummary::collateHarvestingHists ( TTree &  tree )

private

Definition at line 742 of file TrackerOfflineValidationSummary.cc.

References applyHarvestingHierarchy(), and fillHarvestingHists().

Referenced by endJob().

                                                                        {
  this->applyHarvestingHierarchy(tree);
  this->fillHarvestingHists(tree);
}

endJob()

void TrackerOfflineValidationSummary::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 236 of file TrackerOfflineValidationSummary.cc.

References collateHarvestingHists(), fillTree(), mPxbResiduals_, mPxeResiduals_, mTecResiduals_, mTibResiduals_, mTidResiduals_, mTobResiduals_, tkGeom_, compare::tree, and tTopo_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                             {
  AlignableTracker aliTracker(&(*tkGeom_), tTopo_.get());

  TTree* tree = new TTree("TkOffVal", "TkOffVal");

  TkOffTreeVariables* treeMemPtr = new TkOffTreeVariables;
  // We create branches for all members of 'TkOffTreeVariables' (even if not needed).
  // This works because we have a dictionary for 'TkOffTreeVariables'
  // (see src/classes_def.xml and src/classes.h):
  tree->Branch("TkOffTreeVariables", &treeMemPtr);  // address of pointer!
  // second branch needed for assigning names and titles to harvesting histograms consistent to others
  std::map<std::string, std::string>* substructureName = new std::map<std::string, std::string>;
  tree->Branch("SubstructureName", &substructureName, 32000, 00);  // SplitLevel must be set to zero

  this->fillTree(*tree, mPxbResiduals_, *treeMemPtr, *tkGeom_, *substructureName, tTopo_.get());
  this->fillTree(*tree, mPxeResiduals_, *treeMemPtr, *tkGeom_, *substructureName, tTopo_.get());
  this->fillTree(*tree, mTibResiduals_, *treeMemPtr, *tkGeom_, *substructureName, tTopo_.get());
  this->fillTree(*tree, mTidResiduals_, *treeMemPtr, *tkGeom_, *substructureName, tTopo_.get());
  this->fillTree(*tree, mTobResiduals_, *treeMemPtr, *tkGeom_, *substructureName, tTopo_.get());
  this->fillTree(*tree, mTecResiduals_, *treeMemPtr, *tkGeom_, *substructureName, tTopo_.get());

  //dbe_->save("dqmOut.root");

  // Method for filling histograms which show summarized values (mean, rms, median ...)
  // of the module-based histograms from TrackerOfflineValidation
  this->collateHarvestingHists(*tree);

  delete tree;
  tree = nullptr;
  delete treeMemPtr;
  treeMemPtr = nullptr;
  delete substructureName;
  substructureName = nullptr;
}

endRun()

void TrackerOfflineValidationSummary::endRun ( const edm::Run &  , const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 229 of file TrackerOfflineValidationSummary.cc.

References edm::EventSetup::getData(), tTopo_, and tTopoToken_.

                                                                                       {
  if (!tTopo_) {
    tTopo_ = std::make_unique<TrackerTopology>(iSetup.getData(tTopoToken_));
  }
}

fillDescriptions()

void TrackerOfflineValidationSummary::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 166 of file TrackerOfflineValidationSummary.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                 {
  edm::ParameterSetDescription desc;
  desc.setComment("Summary of alignment payloads validation by evaluating unbiased track-hit resisuals");
  desc.add<std::string>("moduleDirectoryInOutput", "Alignment/Tracker");
  desc.add<bool>("useFit", false);
  desc.add<bool>("stripYDmrs", false);
  desc.add<unsigned int>("minEntriesPerModuleForDmr", 100);

  // fill in the residuals details
  const std::vector<std::string> listOfDMRPSets = {
      "TH1DmrXprimeStripModules", "TH1DmrYprimeStripModules", "TH1DmrXprimePixelModules", "TH1DmrYprimePixelModules"};

  for (const auto& myPSetName : listOfDMRPSets) {
    edm::ParameterSetDescription myPSet;
    myPSet.add<int>("Nbinx", 100);
    myPSet.add<double>("xmin", -5.f);
    myPSet.add<double>("xmax", 5.f);
    desc.add<edm::ParameterSetDescription>(myPSetName, myPSet);
  }
  descriptions.addWithDefaultLabel(desc);
}

fillHarvestingHists()

void TrackerOfflineValidationSummary::fillHarvestingHists ( TTree &  tree )

private

Definition at line 844 of file TrackerOfflineValidationSummary.cc.

References TkOffTreeVariables::entries, edm::ParameterSet::getParameter(), TkOffTreeVariables::medianX, TkOffTreeVariables::medianY, parSet_, and vHarvestingHierarchy_.

Referenced by collateHarvestingHists().

                                                                     {
  TkOffTreeVariables* treeMemPtr = nullptr;
  std::map<std::string, std::string>* substructureName = nullptr;
  tree.SetBranchAddress("TkOffTreeVariables", &treeMemPtr);
  tree.SetBranchAddress("SubstructureName", &substructureName);

  const unsigned int minEntriesPerModule(parSet_.getParameter<unsigned int>("minEntriesPerModuleForDmr"));
  edm::LogInfo("TrackerOfflineValidationSummary")
      << "Median of a module is added to DMR plots if it contains at least " << minEntriesPerModule << " hits";

  for (std::vector<HarvestingHierarchy>::iterator iHier = vHarvestingHierarchy_.begin();
       iHier != vHarvestingHierarchy_.end();
       ++iHier) {
    for (std::vector<unsigned int>::const_iterator iTreeEntries = iHier->treeEntries.begin();
         iTreeEntries != iHier->treeEntries.end();
         ++iTreeEntries) {
      tree.GetEntry(*iTreeEntries);
      if (treeMemPtr->entries < minEntriesPerModule)
        continue;
      iHier->harvestingHistos.DmrXprime->Fill(treeMemPtr->medianX);
      if (iHier->harvestingHistos.DmrYprime)
        iHier->harvestingHistos.DmrYprime->Fill(treeMemPtr->medianY);
    }
  }
}

fillTree()

void TrackerOfflineValidationSummary::fillTree ( TTree &  tree, std::map< int, TrackerOfflineValidationSummary::ModuleHistos > &  moduleHist, TkOffTreeVariables &  treeMem, const TrackerGeometry &  tkgeom, std::map< std::string, std::string > &  substructureName, const TrackerTopology *  tTopo  )

private

Definition at line 271 of file TrackerOfflineValidationSummary.cc.

References associateModuleHistsWithTree(), PV3DBase< T, PVType, FrameType >::barePhi(), TkOffTreeVariables::blade, TkOffTreeVariables::chi2PerDofX, TkOffTreeVariables::chi2PerDofY, TkOffTreeVariables::clear(), SiPixelRawToDigiRegional_cfi::deltaPhi, hcalRecHitTable_cff::detId, HLT_2024v10_cff::dPhi, HGC3DClusterGenMatchSelector_cfi::dR, l1ctLayer1_cff::dZ, TkOffTreeVariables::entries, PV3DBase< T, PVType, FrameType >::eta(), TkOffTreeVariables::fitMeanNormX, TkOffTreeVariables::fitMeanNormY, TkOffTreeVariables::fitMeanX, TkOffTreeVariables::fitMeanY, fitResiduals(), TkOffTreeVariables::fitSigmaNormX, TkOffTreeVariables::fitSigmaNormY, TkOffTreeVariables::fitSigmaX, TkOffTreeVariables::fitSigmaY, getMedian(), h, TkOffTreeVariables::half, L1TStage2EMTFEventInfoClient_cfi::histDir, TkOffTreeVariables::histNameLocalX, TkOffTreeVariables::histNameNormLocalX, TkOffTreeVariables::histNameNormX, TkOffTreeVariables::histNameNormY, TkOffTreeVariables::histNameX, TkOffTreeVariables::histNameY, TrackerGeometry::idToDet(), TkOffTreeVariables::isDoubleSide, TkOffTreeVariables::isStereo, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, TkOffTreeVariables::layer, TkOffTreeVariables::meanLocalX, TkOffTreeVariables::meanNormLocalX, TkOffTreeVariables::meanNormX, TkOffTreeVariables::meanNormY, TkOffTreeVariables::meanX, TkOffTreeVariables::meanY, TkOffTreeVariables::medianX, TkOffTreeVariables::medianY, TkOffTreeVariables::module, TkOffTreeVariables::moduleId, PV_cfg::numberOfBins, TkOffTreeVariables::numberOfOutliers, TkOffTreeVariables::numberOfOverflows, TkOffTreeVariables::numberOfUnderflows, TkOffTreeVariables::outerInner, TkOffTreeVariables::panel, PV3DBase< T, PVType, FrameType >::perp(), TkOffTreeVariables::petal, PV3DBase< T, PVType, FrameType >::phi(), TkOffTreeVariables::phiDirection, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, TkOffTreeVariables::posEta, GeomDet::position(), TkOffTreeVariables::posPhi, TkOffTreeVariables::posR, TkOffTreeVariables::posX, TkOffTreeVariables::posY, TkOffTreeVariables::posZ, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), TrackerTopology::pxfBlade(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), TrackerTopology::pxfSide(), nano_mu_digi_cff::rawId, TkOffTreeVariables::rDirection, TkOffTreeVariables::ring, TkOffTreeVariables::rmsLocalX, TkOffTreeVariables::rmsNormLocalX, TkOffTreeVariables::rmsNormX, TkOffTreeVariables::rmsNormY, TkOffTreeVariables::rmsX, TkOffTreeVariables::rmsY, TkOffTreeVariables::rod, TkOffTreeVariables::rOrZDirection, TkOffTreeVariables::side, dqmMemoryStats::stats, AlCaHLTBitMon_QueryRunRegistry::string, TkOffTreeVariables::subDetId, GeomDet::surface(), StripSubdetector::TEC, TrackerTopology::tecIsDoubleSide(), TrackerTopology::tecModule(), TrackerTopology::tecPetalInfo(), TrackerTopology::tecRing(), TrackerTopology::tecSide(), TrackerTopology::tecStereo(), TrackerTopology::tecWheel(), StripSubdetector::TIB, TrackerTopology::tibIsDoubleSide(), TrackerTopology::tibLayer(), TrackerTopology::tibModule(), TrackerTopology::tibStereo(), TrackerTopology::tibStringInfo(), StripSubdetector::TID, TrackerTopology::tidIsDoubleSide(), TrackerTopology::tidModuleInfo(), TrackerTopology::tidRing(), TrackerTopology::tidSide(), TrackerTopology::tidStereo(), TrackerTopology::tidWheel(), StripSubdetector::TOB, TrackerTopology::tobIsDoubleSide(), TrackerTopology::tobLayer(), TrackerTopology::tobModule(), TrackerTopology::tobRodInfo(), TrackerTopology::tobStereo(), Surface::toGlobal(), useFit_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and TkOffTreeVariables::zDirection.

Referenced by endJob(), and core.AutoFillTreeProducer.AutoFillTreeProducer::process().

                                                                             {
  for (std::map<int, TrackerOfflineValidationSummary::ModuleHistos>::iterator it = moduleHist.begin(),
                                                                              itEnd = moduleHist.end();
       it != itEnd;
       ++it) {
    treeMem.clear();  // make empty/default

    //variables concerning the tracker components/hierarchy levels
    const DetId detId = it->first;
    treeMem.moduleId = detId;
    treeMem.subDetId = detId.subdetId();

    if (treeMem.subDetId == PixelSubdetector::PixelBarrel) {
      unsigned int whichHalfBarrel(1), rawId(detId.rawId());  //DetId does not know about halfBarrels is PXB ...
      if ((rawId >= 302056964 && rawId < 302059300) || (rawId >= 302123268 && rawId < 302127140) ||
          (rawId >= 302189572 && rawId < 302194980))
        whichHalfBarrel = 2;
      treeMem.layer = tTopo->pxbLayer(detId);
      treeMem.half = whichHalfBarrel;
      treeMem.rod = tTopo->pxbLadder(
          detId);  // ... so, ladder is not per halfBarrel-Layer, but per barrel-layer! Needs complicated calculation in associateModuleHistsWithTree()
      treeMem.module = tTopo->pxbModule(detId);
    } else if (treeMem.subDetId == PixelSubdetector::PixelEndcap) {
      unsigned int whichHalfCylinder(1), rawId(detId.rawId());  //DetId does not kmow about halfCylinders in PXF
      if ((rawId >= 352394500 && rawId < 352406032) || (rawId >= 352460036 && rawId < 352471568) ||
          (rawId >= 344005892 && rawId < 344017424) || (rawId >= 344071428 && rawId < 344082960))
        whichHalfCylinder = 2;
      treeMem.layer = tTopo->pxfDisk(detId);
      treeMem.side = tTopo->pxfSide(detId);
      treeMem.half = whichHalfCylinder;
      treeMem.blade = tTopo->pxfBlade(detId);
      treeMem.panel = tTopo->pxfPanel(detId);
      treeMem.module = tTopo->pxfModule(detId);
    } else if (treeMem.subDetId == StripSubdetector::TIB) {
      unsigned int whichHalfShell(1), rawId(detId.rawId());  //DetId does not kmow about halfShells in TIB
      if ((rawId >= 369120484 && rawId < 369120688) || (rawId >= 369121540 && rawId < 369121776) ||
          (rawId >= 369136932 && rawId < 369137200) || (rawId >= 369137988 && rawId < 369138288) ||
          (rawId >= 369153396 && rawId < 369153744) || (rawId >= 369154436 && rawId < 369154800) ||
          (rawId >= 369169844 && rawId < 369170256) || (rawId >= 369170900 && rawId < 369171344) ||
          (rawId >= 369124580 && rawId < 369124784) || (rawId >= 369125636 && rawId < 369125872) ||
          (rawId >= 369141028 && rawId < 369141296) || (rawId >= 369142084 && rawId < 369142384) ||
          (rawId >= 369157492 && rawId < 369157840) || (rawId >= 369158532 && rawId < 369158896) ||
          (rawId >= 369173940 && rawId < 369174352) || (rawId >= 369174996 && rawId < 369175440))
        whichHalfShell = 2;
      treeMem.layer = tTopo->tibLayer(detId);
      treeMem.side = tTopo->tibStringInfo(detId)[0];
      treeMem.half = whichHalfShell;
      treeMem.rod = tTopo->tibStringInfo(detId)[2];
      treeMem.outerInner = tTopo->tibStringInfo(detId)[1];
      treeMem.module = tTopo->tibModule(detId);
      treeMem.isStereo = tTopo->tibStereo(detId);
      treeMem.isDoubleSide = tTopo->tibIsDoubleSide(detId);
    } else if (treeMem.subDetId == StripSubdetector::TID) {
      treeMem.layer = tTopo->tidWheel(detId);
      treeMem.side = tTopo->tidSide(detId);
      treeMem.ring = tTopo->tidRing(detId);
      treeMem.outerInner = tTopo->tidModuleInfo(detId)[0];
      treeMem.module = tTopo->tidModuleInfo(detId)[1];
      treeMem.isStereo = tTopo->tidStereo(detId);
      treeMem.isDoubleSide = tTopo->tidIsDoubleSide(detId);
    } else if (treeMem.subDetId == StripSubdetector::TOB) {
      treeMem.layer = tTopo->tobLayer(detId);
      treeMem.side = tTopo->tobRodInfo(detId)[0];
      treeMem.rod = tTopo->tobRodInfo(detId)[1];
      treeMem.module = tTopo->tobModule(detId);
      treeMem.isStereo = tTopo->tobStereo(detId);
      treeMem.isDoubleSide = tTopo->tobIsDoubleSide(detId);
    } else if (treeMem.subDetId == StripSubdetector::TEC) {
      treeMem.layer = tTopo->tecWheel(detId);
      treeMem.side = tTopo->tecSide(detId);
      treeMem.ring = tTopo->tecRing(detId);
      treeMem.petal = tTopo->tecPetalInfo(detId)[1];
      treeMem.outerInner = tTopo->tecPetalInfo(detId)[0];
      treeMem.module = tTopo->tecModule(detId);
      treeMem.isStereo = tTopo->tecStereo(detId);
      treeMem.isDoubleSide = tTopo->tecIsDoubleSide(detId);
    }

    //variables concerning the tracker geometry

    const Surface::PositionType& gPModule = tkgeom.idToDet(detId)->position();
    treeMem.posPhi = gPModule.phi();
    treeMem.posEta = gPModule.eta();
    treeMem.posR = gPModule.perp();
    treeMem.posX = gPModule.x();
    treeMem.posY = gPModule.y();
    treeMem.posZ = gPModule.z();

    const Surface& surface = tkgeom.idToDet(detId)->surface();

    //global Orientation of local coordinate system of dets/detUnits
    LocalPoint lUDirection(1., 0., 0.), lVDirection(0., 1., 0.), lWDirection(0., 0., 1.);
    GlobalPoint gUDirection = surface.toGlobal(lUDirection), gVDirection = surface.toGlobal(lVDirection),
                gWDirection = surface.toGlobal(lWDirection);
    double dR(999.), dPhi(999.), dZ(999.);
    if (treeMem.subDetId == PixelSubdetector::PixelBarrel || treeMem.subDetId == StripSubdetector::TIB ||
        treeMem.subDetId == StripSubdetector::TOB) {
      dR = gWDirection.perp() - gPModule.perp();
      dPhi = deltaPhi(gUDirection.barePhi(), gPModule.barePhi());
      dZ = gVDirection.z() - gPModule.z();
      if (dZ >= 0.)
        treeMem.rOrZDirection = 1;
      else
        treeMem.rOrZDirection = -1;
    } else if (treeMem.subDetId == PixelSubdetector::PixelEndcap) {
      dR = gUDirection.perp() - gPModule.perp();
      dPhi = deltaPhi(gVDirection.barePhi(), gPModule.barePhi());
      dZ = gWDirection.z() - gPModule.z();
      if (dR >= 0.)
        treeMem.rOrZDirection = 1;
      else
        treeMem.rOrZDirection = -1;
    } else if (treeMem.subDetId == StripSubdetector::TID || treeMem.subDetId == StripSubdetector::TEC) {
      dR = gVDirection.perp() - gPModule.perp();
      dPhi = deltaPhi(gUDirection.barePhi(), gPModule.barePhi());
      dZ = gWDirection.z() - gPModule.z();
      if (dR >= 0.)
        treeMem.rOrZDirection = 1;
      else
        treeMem.rOrZDirection = -1;
    }
    if (dR >= 0.)
      treeMem.rDirection = 1;
    else
      treeMem.rDirection = -1;
    if (dPhi >= 0.)
      treeMem.phiDirection = 1;
    else
      treeMem.phiDirection = -1;
    if (dZ >= 0.)
      treeMem.zDirection = 1;
    else
      treeMem.zDirection = -1;

    // Assign histos from first step (TrackerOfflineValidation) to the module's entry in the TTree for retrieving mean, rms, median ...
    const std::string histDir = associateModuleHistsWithTree(treeMem, it->second, substructureName);

    //mean and RMS values (extracted from histograms Xprime on module level)
    treeMem.entries = static_cast<UInt_t>(it->second.ResXprimeHisto->GetEntries());
    treeMem.meanX = it->second.ResXprimeHisto->GetMean();
    treeMem.rmsX = it->second.ResXprimeHisto->GetRMS();
    //treeMem.sigmaX = Fwhm(it->second.ResXprimeHisto)/2.355;
    if (useFit_) {
      //call fit function which returns mean and sigma from the fit
      //for absolute residuals
      std::pair<float, float> fitResult1 = this->fitResiduals(it->second.ResXprimeHisto);
      treeMem.fitMeanX = fitResult1.first;
      treeMem.fitSigmaX = fitResult1.second;
      //for normalized residuals
      std::pair<float, float> fitResult2 = this->fitResiduals(it->second.NormResXprimeHisto);
      treeMem.fitMeanNormX = fitResult2.first;
      treeMem.fitSigmaNormX = fitResult2.second;
    }

    //get median for absolute residuals
    treeMem.medianX = this->getMedian(it->second.ResXprimeHisto);

    int numberOfBins = it->second.ResXprimeHisto->GetNbinsX();
    treeMem.numberOfUnderflows = it->second.ResXprimeHisto->GetBinContent(0);
    treeMem.numberOfOverflows = it->second.ResXprimeHisto->GetBinContent(numberOfBins + 1);
    treeMem.numberOfOutliers =
        it->second.ResXprimeHisto->GetBinContent(0) + it->second.ResXprimeHisto->GetBinContent(numberOfBins + 1);
    //mean and RMS values (extracted from histograms(normalized Xprime on module level)
    treeMem.meanNormX = it->second.NormResXprimeHisto->GetMean();
    treeMem.rmsNormX = it->second.NormResXprimeHisto->GetRMS();

    double stats[20];
    it->second.NormResXprimeHisto->GetStats(stats);
    // GF  treeMem.chi2PerDofX = stats[3]/(stats[0]-1);
    if (stats[0])
      treeMem.chi2PerDofX = stats[3] / stats[0];

    //treeMem.sigmaNormX = Fwhm(it->second.NormResXprimeHisto)/2.355;
    treeMem.histNameX = it->second.ResXprimeHisto->GetName();
    treeMem.histNameNormX = it->second.NormResXprimeHisto->GetName();

    // fill tree variables in local coordinates if set in cfg of TrackerOfllineValidation
    if (it->second.ResHisto && it->second.NormResHisto) {  // if(lCoorHistOn_) {
      treeMem.meanLocalX = it->second.ResHisto->GetMean();
      treeMem.rmsLocalX = it->second.ResHisto->GetRMS();
      treeMem.meanNormLocalX = it->second.NormResHisto->GetMean();
      treeMem.rmsNormLocalX = it->second.NormResHisto->GetRMS();
      treeMem.histNameLocalX = it->second.ResHisto->GetName();
      treeMem.histNameNormLocalX = it->second.NormResHisto->GetName();
    }

    // mean and RMS values in local y (extracted from histograms Yprime on module level)
    // might exist in pixel only
    if (it->second.ResYprimeHisto) {  //(stripYResiduals_){
      TH1* h = it->second.ResYprimeHisto;
      treeMem.meanY = h->GetMean();
      treeMem.rmsY = h->GetRMS();

      if (useFit_) {  // fit function which returns mean and sigma from the fit
        std::pair<float, float> fitMeanSigma = this->fitResiduals(h);
        treeMem.fitMeanY = fitMeanSigma.first;
        treeMem.fitSigmaY = fitMeanSigma.second;
      }

      //get median for absolute residuals
      treeMem.medianY = this->getMedian(h);

      treeMem.histNameY = h->GetName();
    }

    if (it->second.NormResYprimeHisto) {
      TH1* h = it->second.NormResYprimeHisto;
      treeMem.meanNormY = h->GetMean();
      treeMem.rmsNormY = h->GetRMS();
      h->GetStats(stats);  // stats buffer defined above
      if (stats[0])
        treeMem.chi2PerDofY = stats[3] / stats[0];

      if (useFit_) {  // fit function which returns mean and sigma from the fit
        std::pair<float, float> fitMeanSigma = this->fitResiduals(h);
        treeMem.fitMeanNormY = fitMeanSigma.first;
        treeMem.fitSigmaNormY = fitMeanSigma.second;
      }
      treeMem.histNameNormY = h->GetName();
    }

    tree.Fill();
  }
}

fitResiduals()

std::pair< float, float > TrackerOfflineValidationSummary::fitResiduals ( TH1 *  hist ) const

private

Definition at line 687 of file TrackerOfflineValidationSummary.cc.

References MillePedeFileConverter_cfg::e, EcalMonitorTask_cff::func, compareTotals::hist, and SiStripPI::mean.

Referenced by fillTree().

                                                                                   {
  std::pair<float, float> fitResult(9999., 9999.);
  if (!hist || hist->GetEntries() < 20)
    return fitResult;

  float mean = hist->GetMean();
  float sigma = hist->GetRMS();

  try {  // for < CMSSW_2_2_0 since ROOT warnings from fit are converted to exceptions
    // Remove the try/catch for more recent CMSSW!
    // first fit: two RMS around mean
    TF1 func("tmp", "gaus", mean - 2. * sigma, mean + 2. * sigma);
    if (0 == hist->Fit(&func, "QNR")) {  // N: do not blow up file by storing fit!
      mean = func.GetParameter(1);
      sigma = func.GetParameter(2);
      // second fit: three sigma of first fit around mean of first fit
      func.SetRange(mean - 3. * sigma, mean + 3. * sigma);
      // I: integral gives more correct results if binning is too wide
      // L: Likelihood can treat empty bins correctly (if hist not weighted...)
      if (0 == hist->Fit(&func, "Q0LR")) {
        if (hist->GetFunction(func.GetName())) {  // Take care that it is later on drawn:
          hist->GetFunction(func.GetName())->ResetBit(TF1::kNotDraw);
        }
        fitResult.first = func.GetParameter(1);
        fitResult.second = func.GetParameter(2);
      }
    }
  } catch (cms::Exception const& e) {
    edm::LogWarning("Alignment") << "@SUB=TrackerOfflineValidation::fitResiduals"
                                 << "Caught this exception during ROOT fit: " << e.what();
  }
  return fitResult;
}

getBinning()

void TrackerOfflineValidationSummary::getBinning ( const std::string &  binningPSetName, int &  nBinsX, double &  lowerBoundX, double &  upperBoundX  ) const

private

Definition at line 834 of file TrackerOfflineValidationSummary.cc.

References edm::ParameterSet::getParameter(), and parSet_.

Referenced by bookHarvestingHists().

                                                                            {
  const edm::ParameterSet& binningPSet = parSet_.getParameter<edm::ParameterSet>(binningPSetName);
  nBinsX = binningPSet.getParameter<int>("Nbinx");
  lowerBoundX = binningPSet.getParameter<double>("xmin");
  upperBoundX = binningPSet.getParameter<double>("xmax");
}

getMedian()

float TrackerOfflineValidationSummary::getMedian ( const TH1 *  hist ) const

private

Definition at line 721 of file TrackerOfflineValidationSummary.cc.

References timingPdfMaker::histo, dqmiolumiharvest::j, median(), LaserClient_cfi::nbins, x, and y.

Referenced by fillTree().

                                                                       {
  float median = 999;
  const int nbins = histo->GetNbinsX();

  //extract median from histogram
  double* x = new double[nbins];
  double* y = new double[nbins];
  for (int j = 0; j < nbins; j++) {
    x[j] = histo->GetBinCenter(j + 1);
    y[j] = histo->GetBinContent(j + 1);
  }
  median = TMath::Median(nbins, x, y);

  delete[] x;
  x = nullptr;
  delete[] y;
  y = nullptr;

  return median;
}

Member Data Documentation

dbe_

DQMStore* TrackerOfflineValidationSummary::dbe_

private

Definition at line 136 of file TrackerOfflineValidationSummary.cc.

Referenced by associateModuleHistsWithTree(), bookHarvestingHists(), and TrackerOfflineValidationSummary().

geomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TrackerOfflineValidationSummary::geomToken_

private

Definition at line 67 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze().

moduleDirectory_

const std::string TrackerOfflineValidationSummary::moduleDirectory_

private

Definition at line 133 of file TrackerOfflineValidationSummary.cc.

Referenced by associateModuleHistsWithTree(), and bookHarvestingHists().

moduleMapsInitialized

bool TrackerOfflineValidationSummary::moduleMapsInitialized

private

Definition at line 138 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze().

mPxbResiduals_

std::map<int, TrackerOfflineValidationSummary::ModuleHistos> TrackerOfflineValidationSummary::mPxbResiduals_

private

Definition at line 140 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze(), and endJob().

mPxeResiduals_

std::map<int, TrackerOfflineValidationSummary::ModuleHistos> TrackerOfflineValidationSummary::mPxeResiduals_

private

Definition at line 141 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze(), and endJob().

mTecResiduals_

std::map<int, TrackerOfflineValidationSummary::ModuleHistos> TrackerOfflineValidationSummary::mTecResiduals_

private

Definition at line 145 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze(), and endJob().

mTibResiduals_

std::map<int, TrackerOfflineValidationSummary::ModuleHistos> TrackerOfflineValidationSummary::mTibResiduals_

private

Definition at line 142 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze(), and endJob().

mTidResiduals_

std::map<int, TrackerOfflineValidationSummary::ModuleHistos> TrackerOfflineValidationSummary::mTidResiduals_

private

Definition at line 143 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze(), and endJob().

mTobResiduals_

std::map<int, TrackerOfflineValidationSummary::ModuleHistos> TrackerOfflineValidationSummary::mTobResiduals_

private

Definition at line 144 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze(), and endJob().

parSet_

const edm::ParameterSet TrackerOfflineValidationSummary::parSet_

private

Definition at line 128 of file TrackerOfflineValidationSummary.cc.

Referenced by bookHarvestingHists(), fillHarvestingHists(), and getBinning().

tkGeom_

edm::ESHandle<TrackerGeometry> TrackerOfflineValidationSummary::tkGeom_

private

Definition at line 129 of file TrackerOfflineValidationSummary.cc.

Referenced by analyze(), and endJob().

tTopo_

std::unique_ptr<TrackerTopology> TrackerOfflineValidationSummary::tTopo_

private

Definition at line 130 of file TrackerOfflineValidationSummary.cc.

Referenced by endJob(), and endRun().

tTopoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> TrackerOfflineValidationSummary::tTopoToken_

private

Definition at line 68 of file TrackerOfflineValidationSummary.cc.

Referenced by endRun().

useFit_

const bool TrackerOfflineValidationSummary::useFit_

private

Definition at line 134 of file TrackerOfflineValidationSummary.cc.

Referenced by fillTree().

vHarvestingHierarchy_

std::vector<HarvestingHierarchy> TrackerOfflineValidationSummary::vHarvestingHierarchy_

private

Definition at line 147 of file TrackerOfflineValidationSummary.cc.

Referenced by applyHarvestingHierarchy(), bookHarvestingHists(), and fillHarvestingHists().