HGCalTestGuardRing Class Reference

Inheritance diagram for HGCalTestGuardRing:

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

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)

Protected Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	beginJob () override
void	endJob () override
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)

Private Attributes
const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	geomToken_
const double	guardRingOffset_
const std::string	nameSense_
const std::string	waferFile_
std::map< HGCSiliconDetId, int >	waferID_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 37 of file HGCalTestGuardRing.cc.

Constructor & Destructor Documentation

HGCalTestGuardRing()

HGCalTestGuardRing::HGCalTestGuardRing

(

const edm::ParameterSet &

ps

)

Definition at line 55 of file HGCalTestGuardRing.cc.

References nameSense_.

    : nameSense_(ps.getParameter<std::string>("nameSense")),
      waferFile_(ps.getParameter<std::string>("waferFile")),
      guardRingOffset_(ps.getParameter<double>("guardRingOffset")),
      geomToken_(esConsumes<HGCalGeometry, IdealGeometryRecord>(edm::ESInputTag{"", nameSense_})) {
  DetId::Detector det = (nameSense_ != "HGCalHESiliconSensitive") ? DetId::HGCalEE : DetId::HGCalHSi;
  edm::LogVerbatim("HGCalSim") << "Test Guard Ring Offset " << guardRingOffset_ << " for " << nameSense_ << ":" << det
                               << " for wafers read from file " << waferFile_;
  if (!waferFile_.empty()) {
    std::string thick[4] = {"h120", "l200", "l300", "h200"};
    int addType[4] = {HGCalTypes::WaferHD120, HGCalTypes::WaferLD200, HGCalTypes::WaferLD300, HGCalTypes::WaferHD200};
    const int partTypeH[6] = {HGCalTypes::WaferFull,
                              HGCalTypes::WaferHalf2,
                              HGCalTypes::WaferChopTwoM,
                              HGCalTypes::WaferSemi2,
                              HGCalTypes::WaferSemi2,
                              HGCalTypes::WaferFive2};
    const int partTypeL[7] = {HGCalTypes::WaferFull,
                              HGCalTypes::WaferHalf,
                              HGCalTypes::WaferHalf,
                              HGCalTypes::WaferSemi,
                              HGCalTypes::WaferSemi,
                              HGCalTypes::WaferFive,
                              HGCalTypes::WaferThree};
    edm::FileInPath filetmp("SimG4CMS/Calo/data/" + waferFile_);
    std::string fileName = filetmp.fullPath();
    std::ifstream fInput(fileName.c_str());
    if (!fInput.good()) {
      edm::LogVerbatim("HGCalSim") << "Cannot open file " << fileName;
    } else {
      char buffer[80];
      while (fInput.getline(buffer, 80)) {
        std::vector<std::string> items = CaloSimUtils::splitString(std::string(buffer));
        if (items.size() > 6) {
          int layer = std::atoi(items[0].c_str());
          int waferU = std::atoi(items[4].c_str());
          int waferV = std::atoi(items[5].c_str());
          int thck = static_cast<int>(std::find(thick, thick + 4, items[2]) - thick);
          int type = (thck < 4) ? addType[thck] : 0;
          HGCSiliconDetId id(det, -1, type, layer, waferU, waferV, 0, 0);
          int orient = std::atoi(items[5].c_str());
          int part = std::atoi(items[1].c_str());
          if (part >= 0) {
            if ((type == HGCalTypes::WaferHD120) || (type == HGCalTypes::WaferHD200))
              part = partTypeH[part];
            else
              part = partTypeL[part];
          }
          waferID_[id] = orient * 100 + part;
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("HGCalSim") << "HGCalTestGuardRing::Reads " << id << " Orientation:Partial " << orient << ":"
                                       << part;
#endif
        }
      }
      edm::LogVerbatim("HGCalSim") << "HGCalTestGuardRing::Reads in " << waferID_.size() << " wafer information from "
                                   << fileName;
      fInput.close();
    }
  }
}

~HGCalTestGuardRing()

HGCalTestGuardRing::~HGCalTestGuardRing ( )

overridedefault

Member Function Documentation

analyze()

void HGCalTestGuardRing::analyze ( edm::Event const &  e, edm::EventSetup const &  iS  )

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 125 of file HGCalTestGuardRing.cc.

References python.cmstools::all(), dqmdumpme::first, relativeConstraints::geom, geomToken_, edm::EventSetup::getData(), caHitNtupletGeneratorKernels::good, guardRingOffset_, DetId::HGCalEE, DetId::HGCalHSi, HGCGuardRing::insidePolygon(), nano_mu_digi_cff::layer, nameSense_, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::nCells, hgcalPerformanceValidation::orient, edm::second(), findQualityFiles::v, validateGeometry_cfg::valid, HGCalTypes::waferFrontBack(), waferID_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), HGCalWaferMask::waferXY(), geometryCSVtoXML::xy, and ecaldqm::zside().

                                                                           {
  // get HGCalGeometry
  const HGCalGeometry* geom = &iS.getData(geomToken_);
  const HGCalDDDConstants& hgc = geom->topology().dddConstants();
  double waferSize = hgc.waferSize(false);
  HGCalCell wafer(waferSize, hgc.getUVMax(0), hgc.getUVMax(1));
  const bool v17OrLess = hgc.v17OrLess();

  // get the hit collection
  edm::LogVerbatim("HGCalSim") << "HGCalTestGuardRing: Wafer Szie " << waferSize << " v17OrLess " << v17OrLess;

  // Loop over all IDs
  int all(0), allSi(0), good(0);
  for (std::map<HGCSiliconDetId, int>::const_iterator itr = waferID_.begin(); itr != waferID_.end(); ++itr) {
    HGCSiliconDetId id = itr->first;
    ++all;
    if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
      ++allSi;
      if (((id.det() == DetId::HGCalEE) && (nameSense_ == "HGCalEESensitive")) ||
          ((id.det() == DetId::HGCalHSi) && (nameSense_ == "HGCalHESiliconSensitive"))) {
        int partial = ((itr->second) % 100);
        int orient = (((itr->second) / 100) % 100);
        int type = id.type();
        int nCells = hgc.getUVMax(type);
        for (int u = 0; u < 2 * nCells; ++u) {
          for (int v = 0; v < 2 * nCells; ++v) {
            if (((v - u) < nCells) && ((u - v) <= nCells)) {
              HGCSiliconDetId hid(id.det(), id.zside(), id.type(), id.layer(), id.waferU(), id.waferV(), u, v);
              bool valid = (geom->topology()).valid(static_cast<DetId>(hid));
              if (valid) {
                ++good;
                int placeIndex = wafer.cellPlacementIndex(1, HGCalTypes::waferFrontBack(0), orient);
                std::pair<double, double> xy = wafer.cellUV2XY1(u, v, placeIndex, type);
                std::vector<std::pair<double, double> > wxy1 =
                    HGCalWaferMask::waferXY(partial, orient, -1, waferSize, 0.0, 0.0, 0.0, v17OrLess);
                bool check1 = HGCGuardRing::insidePolygon(xy.first, xy.second, wxy1);
                std::ostringstream st1;
                for (unsigned int k1 = 0; k1 < wxy1.size(); ++k1)
                  st1 << " (" << wxy1[k1].first << ", " << wxy1[k1].second << ")";
                edm::LogVerbatim("HGCSim")
                    << "First " << hid << " Type:Partial:Orient:Place " << type << ":" << partial << ":" << orient
                    << ":" << placeIndex << " Boundary with " << wxy1.size() << " points: " << st1.str() << " check "
                    << check1 << " for (" << xy.first << ", " << xy.second << ")";

                std::vector<std::pair<double, double> > wxy2 =
                    HGCalWaferMask::waferXY(partial, orient, -1, waferSize, guardRingOffset_, 0.0, 0.0, v17OrLess);
                bool check2 = HGCGuardRing::insidePolygon(xy.first, xy.second, wxy2);
                std::ostringstream st2;
                for (unsigned int k1 = 0; k1 < wxy2.size(); ++k1)
                  st2 << " (" << wxy2[k1].first << ", " << wxy2[k1].second << ")";
                edm::LogVerbatim("HGCSim") << "Second Offset " << guardRingOffset_ << " Boundary with " << wxy2.size()
                                           << " points: " << st2.str() << " check " << check2 << " for (" << xy.first
                                           << ", " << xy.second << ")";
              }
            }
          }
        }
      }
    }
  }
  edm::LogVerbatim("HGCalSim") << "Total hits = " << all << " Good Silicon DetIds = " << allSi << ":" << good;
}

beginJob()

void HGCalTestGuardRing::beginJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 45 of file HGCalTestGuardRing.cc.

{}

endJob()

void HGCalTestGuardRing::endJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 46 of file HGCalTestGuardRing.cc.

{}

fillDescriptions()

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

static

Definition at line 117 of file HGCalTestGuardRing.cc.

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

                                                                                    {
  edm::ParameterSetDescription desc;
  desc.add<std::string>("nameSense", "HGCalEESensitive");
  desc.add<std::string>("waferFile", "testWafersEE.txt");
  desc.add<double>("guardRingOffset", 1.0);
  descriptions.add("hgcalTestGuardRingEE", desc);
}

Member Data Documentation

geomToken_

const edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGCalTestGuardRing::geomToken_

private

Definition at line 51 of file HGCalTestGuardRing.cc.

Referenced by analyze().

guardRingOffset_

const double HGCalTestGuardRing::guardRingOffset_

private

Definition at line 50 of file HGCalTestGuardRing.cc.

Referenced by analyze().

nameSense_

const std::string HGCalTestGuardRing::nameSense_

private

Definition at line 49 of file HGCalTestGuardRing.cc.

Referenced by analyze(), and HGCalTestGuardRing().

waferFile_

const std::string HGCalTestGuardRing::waferFile_

private

Definition at line 49 of file HGCalTestGuardRing.cc.

waferID_

std::map<HGCSiliconDetId, int> HGCalTestGuardRing::waferID_

private

Definition at line 52 of file HGCalTestGuardRing.cc.

Referenced by analyze().