#include <HGCalMouseBiteTester.cc>

Inheritance diagram for HGCalMouseBiteTester:

Public Member Functions
void	analyze (edm::Event const &iEvent, edm::EventSetup const &) override

void	beginJob () override

void	endJob () override

	HGCalMouseBiteTester (const edm::ParameterSet &)

	~HGCalMouseBiteTester () 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 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 Attributes
const edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord >	dddToken_

const int	layer_

const std::string	nameSense_

const int	nTrials_

std::ofstream	outputFile

const int	waferU_

const int	waferV_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

plugins/HGCalMouseBiteTester.cc

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

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

Definition at line 53 of file HGCalMouseBiteTester.cc.

Constructor & Destructor Documentation

◆ HGCalMouseBiteTester()

HGCalMouseBiteTester::HGCalMouseBiteTester ( const edm::ParameterSet & iC )

explicit

Definition at line 73 of file HGCalMouseBiteTester.cc.

References nameSense_.

     : nameSense_(iC.getParameter<std::string>("nameSense")),
       waferU_(iC.getParameter<int>("waferU")),
       waferV_(iC.getParameter<int>("waferV")),
       nTrials_(iC.getParameter<int>("numbberOfTrials")),
       layer_(iC.getParameter<int>("layer")),
       dddToken_(esConsumes<HGCalDDDConstants, IdealGeometryRecord>(edm::ESInputTag{"", nameSense_})) {
   edm::LogVerbatim("HGCalGeom") << "Test Guard_Ring for wafer in layer" << layer_ << " U " << waferU_ << " V "
                                 << waferV_ << " with " << nTrials_ << " trials";
 
   outputFile.open("full1.csv");
   if (!outputFile.is_open()) {
     edm::LogError("HGCalGeom") << "Could not open output file.";
   } else {
     outputFile << "x,y,u,v,\n";
   }
 }

◆ ~HGCalMouseBiteTester()

HGCalMouseBiteTester::~HGCalMouseBiteTester ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 102 of file HGCalMouseBiteTester.cc.

References HGCalCell::cellPlacementIndex(), HGCalCell::cellUV2XY2(), HGCalCellUV::cellUVFromXY1(), submitPVResolutionJobs::count, gather_cfg::cout, dddToken_, HGCGuardRingPartial::exclude(), HGCMouseBite::exclude(), edm::EventSetup::getData(), HGCalWaferType::getOrient(), HGCalDDDConstants::getParameter(), HGCalWaferType::getPartial(), HGCalWaferType::getType(), mps_fire::i, layer_, HGCalTypes::layerFrontBack(), HGCalDDDConstants::layerType(), caHitNtupletGeneratorKernels::nCells, submitPVValidationJobs::now, nTrials_, hgcalPerformanceValidation::orient, outputFile, point, diffTwoXMLs::r2, mathSSE::sqrt(), funct::tan(), HGCalDDDConstants::v17OrLess(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalDDDConstants::waferPosition(), HGCalDDDConstants::waferSize(), waferU_, waferV_, HGCalWaferMask::waferXY(), protons_cff::xi, geometryCSVtoXML::xy, and ecaldqm::zside().

                                                                                       {
   const HGCalDDDConstants& hgcons_ = iSetup.getData(dddToken_);
   double waferSize_(hgcons_.waferSize(false));
   int zside(1);
   int layertype = hgcons_.layerType(layer_);
   int frontBack = HGCalTypes::layerFrontBack(layertype);
   const std::vector<double> angle_{90.0, 30.0};
   int index = HGCalWaferIndex::waferIndex(layer_, waferU_, waferV_);
   int partialType_ = 0;
   int orient = HGCalWaferType::getOrient(index, hgcons_.getParameter()->waferInfoMap_);
   int placeIndex_ = HGCalCell::cellPlacementIndex(zside, frontBack, orient);
   int waferType_ = HGCalWaferType::getType(index, hgcons_.getParameter()->waferInfoMap_);
   double mouseBiteCut_ = waferSize_ * tan(30.0 * CLHEP::deg) - 5.0;
   bool v17OrLess = hgcons_.v17OrLess();
   HGCGuardRing guardRing_(hgcons_);
   HGCGuardRingPartial guardRingPartial_(hgcons_);
   HGCMouseBite mouseBite_(hgcons_, angle_, mouseBiteCut_, true);
   const int nFine(12), nCoarse(8);
   double r2 = 0.5 * waferSize_;
   double R2 = 2 * r2 / sqrt(3);
   int nCells = (waferType_ == 0) ? nFine : nCoarse;
   std::cout << "start" << std::endl;
   HGCalCellUV wafer(waferSize_, 0.0, nFine, nCoarse);
   HGCalCell wafer2(waferSize_, nFine, nCoarse);
   std::pair<double, double> xy = hgcons_.waferPosition(layer_, waferU_, waferV_, false, false);
   double x0 = (zside > 0) ? xy.first : -xy.first;
   double y0 = xy.second;
   std::ofstream guard_ring("Guard_ring.csv");
   std::ofstream guard_ring_partial("Guard_ring_partial.csv");
   std::ofstream mouse_bite("Mouse_bite.csv");
   std::ofstream selected("Selected.csv");
   edm::LogVerbatim("HGCalGeom") << "\nHGCalMouseBiteTester:: nCells " << nCells << " FrontBack " << frontBack
                                 << " Wafer Size " << waferSize_ << " and placement index " << placeIndex_
                                 << " WaferType " << waferType_ << " Partial " << partialType_ << " WaferX " << x0
                                 << " WaferY " << y0 << "\n\n";
   auto start_t = std::chrono::high_resolution_clock::now();
 
   for (int i = 0; i < nTrials_; i++) {
     double xi = (2 * r2 * static_cast<double>(rand()) / RAND_MAX) - r2;
     double yi = (2 * R2 * static_cast<double>(rand()) / RAND_MAX) - R2;
     bool goodPoint = true;
     int ug = 0;
     int vg = 0;
     if (partialType_ == 11 || partialType_ == 13 || partialType_ == 15 || partialType_ == 21 || partialType_ == 23 ||
         partialType_ == 25 || partialType_ == 0) {
       ug = 0;
       vg = 0;
     } else if (partialType_ == 12 || partialType_ == 14 || partialType_ == 16 || partialType_ == 22 ||
                partialType_ == 24) {
       ug = nCells + 1;
       vg = 2 * (nCells - 1);
     }
     std::pair<double, double> xyg = wafer2.cellUV2XY2(ug, vg, placeIndex_, waferType_);
     std::vector<std::pair<double, double> > wxy =
         HGCalWaferMask::waferXY(0, placeIndex_, waferSize_, 0.0, 0.0, 0.0, v17OrLess);
     for (unsigned int i = 0; i < (wxy.size() - 1); ++i) {
       double xp1 = wxy[i].first;
       double yp1 = wxy[i].second;
       double xp2 = wxy[i + 1].first;
       double yp2 = wxy[i + 1].second;
       if ((((xi - xp1) / (xp2 - xp1)) - ((yi - yp1) / (yp2 - yp1))) *
               (((xyg.first - xp1) / (xp2 - xp1)) - ((xyg.second - yp1) / (yp2 - yp1))) <=
           0) {
         goodPoint = false;
       }
     }
     if (goodPoint) {  //Only allowing (x, y) inside a partial wafer 11, placement index 2
       partialType_ = HGCalWaferType::getPartial(index, hgcons_.getParameter()->waferInfoMap_);
       G4ThreeVector point(xi, yi, 0.0);
       std::pair<int32_t, int32_t> uv5 = wafer.cellUVFromXY1(xi, yi, placeIndex_, waferType_, partialType_, true, false);
       if (guardRingPartial_.exclude(point, zside, frontBack, layer_, waferU_, waferV_)) {
         guard_ring_partial << xi << "," << yi << std::endl;
       } else if (mouseBite_.exclude(point, zside, layer_, waferU_, waferV_)) {
         mouse_bite << xi << "," << yi << std::endl;
       } else {
         selected << xi << "," << yi << std::endl;
         outputFile << xi << "," << yi << "," << uv5.first << "," << uv5.second << "," << std::endl;
       }
     }
   }
   guard_ring.close();
   guard_ring_partial.close();
   mouse_bite.close();
   selected.close();
   outputFile.close();
   auto end_t = std::chrono::high_resolution_clock::now();
   auto diff_t = end_t - start_t;
   edm::LogVerbatim("HGCalGeom") << "Execution time for " << nTrials_
                                 << " events = " << std::chrono::duration<double, std::milli>(diff_t).count() << " ms";
 }

◆ beginJob()

void HGCalMouseBiteTester::beginJob ( void )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 59 of file HGCalMouseBiteTester.cc.

59 {}

◆ endJob()

void HGCalMouseBiteTester::endJob ( void )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 61 of file HGCalMouseBiteTester.cc.

61 {}

◆ fillDescriptions()

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

static

Definition at line 91 of file HGCalMouseBiteTester.cc.

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

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("nameSense", "HGCalEESensitive");
   desc.add<int>("waferU", 1);
   desc.add<int>("waferV", 9);
   desc.add<int>("numbberOfTrials", 1000000);
   desc.add<int>("layer", 1);
   descriptions.add("hgcalMouseBiteTester", desc);
 }

Member Data Documentation

◆ dddToken_

const edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> HGCalMouseBiteTester::dddToken_

private

Definition at line 69 of file HGCalMouseBiteTester.cc.

Referenced by analyze().

◆ layer_

const int HGCalMouseBiteTester::layer_

private

Definition at line 68 of file HGCalMouseBiteTester.cc.

Referenced by analyze().

◆ nameSense_

const std::string HGCalMouseBiteTester::nameSense_

private

Definition at line 64 of file HGCalMouseBiteTester.cc.

Referenced by HGCalMouseBiteTester().

◆ nTrials_

const int HGCalMouseBiteTester::nTrials_

private

Definition at line 67 of file HGCalMouseBiteTester.cc.

Referenced by analyze().

◆ outputFile

std::ofstream HGCalMouseBiteTester::outputFile

private

Definition at line 70 of file HGCalMouseBiteTester.cc.

Referenced by analyze().

◆ waferU_

const int HGCalMouseBiteTester::waferU_

private

Definition at line 65 of file HGCalMouseBiteTester.cc.

Referenced by analyze().

◆ waferV_

const int HGCalMouseBiteTester::waferV_

private

Definition at line 66 of file HGCalMouseBiteTester.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalMouseBiteTester()

◆ ~HGCalMouseBiteTester()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

Member Data Documentation

◆ dddToken_

◆ layer_

◆ nameSense_

◆ nTrials_

◆ outputFile

◆ waferU_

◆ waferV_