Inheritance diagram for PPSPixelDigiAnalyzer:

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

void	beginJob () override

void	endJob () override

	PPSPixelDigiAnalyzer (const edm::ParameterSet &pset)

	~PPSPixelDigiAnalyzer () 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
unsigned int	cumulative_cluster_size_ [3]

unsigned int	found_corresponding_digi_count_

TH2D *	hAllHits

TH2D *	hOneHitperEvent

TH2D *	hOneHitperEvent2

TH2D *	hOneHitperEvent2Center

TH2D *	hOneHitperEventCenter

std::string	label_

edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelDigi > >	pixel_token

edm::ESGetToken< PPSPixelTopology, PPSPixelTopologyRcd >	pixelTopologyToken_

edm::EDGetTokenT< edm::PSimHitContainer >	psim_token

int	verbosity_

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 51 of file PPSPixelDigiAnalyzer.cc.

Constructor & Destructor Documentation

◆ PPSPixelDigiAnalyzer()

PPSPixelDigiAnalyzer::PPSPixelDigiAnalyzer ( const edm::ParameterSet & pset )

explicit

Definition at line 77 of file PPSPixelDigiAnalyzer.cc.

References geometryDiff::file, hAllHits, hOneHitperEvent, hOneHitperEvent2, hOneHitperEvent2Center, hOneHitperEventCenter, HLT_2022v15_cff::InputTag, label_, pixel_token, pixelTopologyToken_, muonDTDigis_cfi::pset, psim_token, and verbosity_.

     : hAllHits(nullptr),
       hOneHitperEvent(nullptr),
       hOneHitperEvent2(nullptr),
       hOneHitperEventCenter(nullptr),
       hOneHitperEvent2Center(nullptr) {
   label_ = pset.getUntrackedParameter<string>("label");
   verbosity_ = pset.getParameter<int>("Verbosity");
   edm::Service<TFileService> file;
 #ifdef USE_MIDDLE_OF_PIXEL
   hOneHitperEvent = file->make<TH2D>("OneHitperEvent", "One Hit per Event", 30, -8.511, -8.361, 20, 1, 1.1);
   hOneHitperEvent2 = file->make<TH2D>("OneHitperEvent2", "One Hit per Event 2", 30, -8.511, -8.361, 20, 1, 1.1);
 #else
   hOneHitperEvent = file->make<TH2D>("OneHitperEvent", "One Hit per Event", 30, -8.55, -8.4, 20, 1, 1.1);
   hOneHitperEvent2 = file->make<TH2D>("OneHitperEvent2", "One Hit per Event 2", 30, -8.55, -8.4, 20, 1, 1.1);
   hOneHitperEventCenter =
       file->make<TH2D>("OneHitperEventCenter", "One Hit per Event Center", 30, -0.075, 0.075, 20, -0.05, 0.05);
   hOneHitperEvent2Center =
       file->make<TH2D>("OneHitperEvent2Center", "Cluster Size 2", 30, -0.075, 0.075, 20, -0.05, 0.05);
 #endif
   file->cd();
   hAllHits = file->make<TH2D>("AllHits", "All Hits", 10, 1, 1.1, 10, -8.55, -8.4);
 
   psim_token = consumes<PSimHitContainer>(edm::InputTag("g4SimHits", "CTPPSPixelHits"));
   pixel_token = consumes<edm::DetSetVector<CTPPSPixelDigi>>(edm::InputTag(label_, ""));  //label=RPixDetDigitizer???
   pixelTopologyToken_ = esConsumes<PPSPixelTopology, PPSPixelTopologyRcd>();
 }

◆ ~PPSPixelDigiAnalyzer()

PPSPixelDigiAnalyzer::~PPSPixelDigiAnalyzer ( )

override

Definition at line 105 of file PPSPixelDigiAnalyzer.cc.

105 {}

Member Function Documentation

◆ analyze()

void PPSPixelDigiAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 118 of file PPSPixelDigiAnalyzer.cc.

References gpuPixelDoublets::cc, CENTERX, CENTERY, cumulative_cluster_size_, mps_fire::end, options_cfi::eventSetup, found_corresponding_digi_count_, hAllHits, hOneHitperEvent, hOneHitperEvent2, hOneHitperEvent2Center, hOneHitperEventCenter, pixel_token, PPSPixelTopology::pixelRange(), pixelTopologyToken_, psim_token, findQualityFiles::rr, SELECTED_PIXEL_COLUMN, SELECTED_PIXEL_ROW, SELECTED_UNITID, FastTrackerRecHitCombiner_cfi::simHits, TG184, USE_MIDDLE_OF_PIXEL_2, verbosity_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                    {
   if (verbosity_ > 0)
     edm::LogInfo("PPSPixelDigiAnalyzer") << "--- Run: " << event.id().run() << " Event: " << event.id().event();
 
   edm::LogInfo("PPSPixelDigiAnalyzer")
       << "                                                            I do love Pixels     ";
   Handle<PSimHitContainer> simHits;
   event.getByToken(psim_token, simHits);
 
   edm::Handle<edm::DetSetVector<CTPPSPixelDigi>> CTPPSPixelDigis;
   event.getByToken(pixel_token, CTPPSPixelDigis);
 
   edm::ESHandle<PPSPixelTopology> thePixelTopology = eventSetup.getHandle(pixelTopologyToken_);
 
   if (verbosity_ > 0)
     edm::LogInfo("PPSPixelDigiAnalyzer") << "\n=================== RPDA Starting SimHit access"
                                          << "  ===================";
 
   if (verbosity_ > 1)
     edm::LogInfo("PPSPixelDigiAnalyzer") << simHits->size();
 
   double selected_pixel_lower_x;
   double selected_pixel_lower_y;
   double selected_pixel_upper_x;
   double selected_pixel_upper_y;
   double myX = 0;
   double myY = 0;
 
   thePixelTopology->pixelRange(SELECTED_PIXEL_ROW,
                                SELECTED_PIXEL_COLUMN,
                                selected_pixel_lower_x,
                                selected_pixel_upper_x,
                                selected_pixel_lower_y,
                                selected_pixel_upper_y);
 
   double hit_inside_selected_pixel[2];
   bool found_hit_inside_selected_pixel = false;
 
   for (vector<PSimHit>::const_iterator hit = simHits->begin(); hit != simHits->end(); hit++) {
     LocalPoint entryP = hit->entryPoint();
     LocalPoint exitP = hit->exitPoint();
     LocalPoint midP((entryP.x() + exitP.x()) / 2., (entryP.y() + exitP.y()) / 2.);
 
 #ifdef USE_MIDDLE_OF_PIXEL
     if (entryP.x() > selected_pixel_lower_x && entryP.x() < selected_pixel_upper_x &&
         entryP.y() > (selected_pixel_lower_y + 0.115 * TG184) && entryP.y() < (selected_pixel_upper_y + 0.115 * TG184)
 #else
 #ifdef USE_MIDDLE_OF_PIXEL_2
     if (midP.x() > selected_pixel_lower_x && midP.x() < selected_pixel_upper_x && midP.y() > selected_pixel_lower_y &&
         midP.y() < selected_pixel_upper_y
 #else
     if (entryP.x() > selected_pixel_lower_x && entryP.x() < selected_pixel_upper_x &&
         entryP.y() > selected_pixel_lower_y && entryP.y() < selected_pixel_upper_y
 #endif
 #endif
         && hit->detUnitId() == SELECTED_UNITID) {
       hit_inside_selected_pixel[0] = entryP.x();
       hit_inside_selected_pixel[1] = entryP.y();
       found_hit_inside_selected_pixel = true;
 #ifdef USE_MIDDLE_OF_PIXEL_2
       hAllHits->Fill(midP.x(), midP.y());
       myX = midP.x();
       myY = midP.y();
 #else
       hAllHits->Fill(entryP.x(), entryP.y());
       myX = entryP.x();
       myY = entryP.y();
 #endif
       if (verbosity_ > 2)
         edm::LogInfo("PPSPixelDigiAnalyzer") << hit_inside_selected_pixel[0] << " " << hit_inside_selected_pixel[1];
     }
 
     //--------------
 
     if (verbosity_ > 1)
       if (hit->timeOfFlight() > 0) {
         edm::LogInfo("PPSPixelDigiAnalyzer")
             << "DetId: " << hit->detUnitId() << "PID: " << hit->particleType() << " TOF: " << hit->timeOfFlight()
             << " Proc Type: " << hit->processType() << " p: " << hit->pabs() << " x = " << entryP.x()
             << "   y = " << entryP.y() << "  z = " << entryP.z();
       }
   }
 
   if (verbosity_ > 0)
     edm::LogInfo("PPSPixelDigiAnalyzer") << "\n=================== RPDA Starting Digi access"
                                          << "  ===================";
   int numberOfDetUnits = 0;
 
   // Iterate on detector units
   edm::DetSetVector<CTPPSPixelDigi>::const_iterator DSViter = CTPPSPixelDigis->begin();
 
   for (; DSViter != CTPPSPixelDigis->end(); DSViter++) {
     ++numberOfDetUnits;
 
     DetId detIdObject(DSViter->detId());
     if (verbosity_ > 1)
       edm::LogInfo("PPSPixelDigiAnalyzer") << "DetId: " << DSViter->detId();
 
     bool found_corresponding_digi = false;
     unsigned int corresponding_digi_cluster_size = 0;
 
     // looping over digis in a unit id
     edm::DetSet<CTPPSPixelDigi>::const_iterator begin = (*DSViter).begin();
     edm::DetSet<CTPPSPixelDigi>::const_iterator end = (*DSViter).end();
 
     if (verbosity_ > 2) {
       edm::LogInfo("PPSPixelDigiAnalyzer") << "FF  " << DSViter->detId();
       for (edm::DetSet<CTPPSPixelDigi>::const_iterator di = begin; di != end; di++) {
         edm::LogInfo("PPSPixelDigiAnalyzer") << "           Digi row  " << di->row() << ", col " << di->column();
 
         // reconvert the digi to local coordinates
         double lx;
         double ly;
         double ux;
         double uy;
         unsigned int rr = di->row();
         unsigned int cc = di->column();
         thePixelTopology->pixelRange(rr, cc, lx, ux, ly, uy);
 
         edm::LogInfo("PPSPixelDigiAnalyzer")
             << " pixel boundaries x low up, y low up " << lx << " " << ux << " " << ly << " " << uy;
       }
     }
     if (DSViter->detId() == SELECTED_UNITID && found_hit_inside_selected_pixel) {
       for (edm::DetSet<CTPPSPixelDigi>::const_iterator di = begin; di != end; di++) {
         if (verbosity_ > 1)
           edm::LogInfo("PPSPixelDigiAnalyzer") << "           Digi row  " << di->row() << ", col " << di->column();
 
         if (di->row() == SELECTED_PIXEL_ROW && di->column() == SELECTED_PIXEL_COLUMN) {
           found_corresponding_digi_count_++;
           found_corresponding_digi = true;
           corresponding_digi_cluster_size = 1;
         }
       }
       //if coresponding digi found, re-loop to look for adjacent pixels
       if (found_corresponding_digi) {
         for (edm::DetSet<CTPPSPixelDigi>::const_iterator di = begin; di != end; di++) {
           if (verbosity_ > 1)
             edm::LogInfo("PPSPixelDigiAnalyzer") << "           Digi row  " << di->row() << ", col " << di->column();
 
           if ((di->row() == SELECTED_PIXEL_ROW + 1 && di->column() == SELECTED_PIXEL_COLUMN) ||
               (di->row() == SELECTED_PIXEL_ROW - 1 && di->column() == SELECTED_PIXEL_COLUMN) ||
               (di->row() == SELECTED_PIXEL_ROW && di->column() == SELECTED_PIXEL_COLUMN + 1) ||
               (di->row() == SELECTED_PIXEL_ROW && di->column() == SELECTED_PIXEL_COLUMN - 1)) {
             corresponding_digi_cluster_size++;
             edm::LogInfo("PPSPixelDigiAnalyzer") << "           Digi row  " << di->row() << ", col " << di->column();
           }
         }
       }
     }
     if (corresponding_digi_cluster_size > 0) {
       edm::LogInfo("PPSPixelDigiAnalyzer")
           << "corresponding_digi_cluster_size in the event: " << corresponding_digi_cluster_size;
       hOneHitperEvent->Fill(myY, myX);
       hOneHitperEventCenter->Fill(myY - CENTERY, myX - CENTERX);
       if (corresponding_digi_cluster_size < 3) {
         cumulative_cluster_size_[corresponding_digi_cluster_size - 1]++;
         if (corresponding_digi_cluster_size > 1) {
           hOneHitperEvent2->Fill(myY, myX);
           hOneHitperEvent2Center->Fill(myY - CENTERY, myX - CENTERX);
         }
       } else {
         cumulative_cluster_size_[2]++;
         hOneHitperEvent2->Fill(myY, myX);
         hOneHitperEvent2Center->Fill(myY - CENTERY, myX - CENTERX);
       }
     }
   }
 
   if (verbosity_ > 1)
     edm::LogInfo("PPSPixelDigiAnalyzer") << "numberOfDetUnits in the event: " << numberOfDetUnits;
 }

◆ beginJob()

void PPSPixelDigiAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 107 of file PPSPixelDigiAnalyzer.cc.

References a, cumulative_cluster_size_, and found_corresponding_digi_count_.

                                     {
   found_corresponding_digi_count_ = 0;
   for (int a = 0; a < 3; a++)
     cumulative_cluster_size_[a] = 0;
 }

◆ endJob()

void PPSPixelDigiAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 112 of file PPSPixelDigiAnalyzer.cc.

References cumulative_cluster_size_, and found_corresponding_digi_count_.

                                   {
   edm::LogInfo("PPSPixelDigiAnalyzer") << "found_corresponding_digi_count_: " << found_corresponding_digi_count_;
   edm::LogInfo("PPSPixelDigiAnalyzer") << "Cumulative cluster size (1,2,>2) = " << cumulative_cluster_size_[0] << ", "
                                        << cumulative_cluster_size_[1] << ", " << cumulative_cluster_size_[2];
 }