Inheritance diagram for PixelBaryCentreAnalyzer:

Public Member Functions
	PixelBaryCentreAnalyzer (const edm::ParameterSet &)

	~PixelBaryCentreAnalyzer () 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)

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 &, const edm::EventSetup &) override

void	beginJob () override

void	endJob () override

void	initBC ()

void	initBS ()

Private Attributes
std::vector< std::string >	bcLabels_

std::map< std::string, TTree * >	bcTrees_

GlobalPoint	BPIX_

GlobalPoint	BPIX_DiffFlippedNonFlipped_

GlobalPoint	BPIX_Flipped_

GlobalPoint	BPIX_NonFlipped_

GlobalPoint	BPIXLayer_ [4]

GlobalPoint	BPIXLayer_DiffFlippedNonFlipped_ [4]

GlobalPoint	BPIXLayer_Flipped_ [4]

GlobalPoint	BPIXLayer_NonFlipped_ [4]

GlobalPoint	BS_

std::vector< std::string >	bsLabels_

std::map< std::string, edm::ESGetToken< BeamSpotObjects, BeamSpotObjectsRcd > >	bsTokens_

std::map< std::string, TTree * >	bsTrees_

GlobalPoint	FPIX_

GlobalPoint	FPIX_minus_

GlobalPoint	FPIX_plus_

GlobalPoint	FPIXDisks_minus_ [3]

GlobalPoint	FPIXDisks_plus_ [3]

edm::ESGetToken< Alignments, GlobalPositionRcd >	gprToken_

int	ls_

int	phase_

GlobalPoint	PIX_

int	run_

edm::ESGetToken< SiPixelQuality, SiPixelQualityFromDbRcd >	siPixelQualityToken_

edm::Service< TFileService >	tFileService

std::map< std::string, edm::ESGetToken< Alignments, TrackerAlignmentRcd > >	tkAlignTokens_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopologyToken_

bool	usePixelQuality_

edm::ESWatcher< BeamSpotObjectsRcd >	watcherBS_

edm::ESWatcher< TrackerAlignmentRcd >	watcherTkAlign_

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 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

The analyer works as the following :

Read global tracker position from global tag
Read tracker alignment constants from different ESsource with different labels
Calculate barycentres for different pixel substructures using global tracker position and alignment constants and store them in trees, one for each ESsource label.

Python script plotBaryCentre_VS_BeamSpot.py under script dir is used to plot barycentres from alignment constants used in Prompt-Reco, End-of-Year Rereco and so-called Run-2 (Ultra)Legacy Rereco. Options of the plotting script can be found from the helper in the script.

Date: 2021/01/05

Revision: 1.0

Author: Tongguang Cheng - Beihang University tongg.nosp@m.uang.nosp@m..chen.nosp@m.g@ce.nosp@m.rn.ch

Definition at line 62 of file PixelBaryCentreAnalyzer.cc.

Constructor & Destructor Documentation

◆ PixelBaryCentreAnalyzer()

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

explicit

Definition at line 123 of file PixelBaryCentreAnalyzer.cc.

References bcLabels_, bcTrees_, bsLabels_, bsTokens_, bsTrees_, label, and tkAlignTokens_.

     : usePixelQuality_(iConfig.getUntrackedParameter<bool>("usePixelQuality")),
       bcLabels_(iConfig.getUntrackedParameter<std::vector<std::string>>("tkAlignLabels")),
       bsLabels_(iConfig.getUntrackedParameter<std::vector<std::string>>("beamSpotLabels")),
       trackerGeometryToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()),
       trackerTopologyToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>()),
       siPixelQualityToken_(esConsumes<SiPixelQuality, SiPixelQualityFromDbRcd>()),
       gprToken_(esConsumes<Alignments, GlobalPositionRcd>()) {
   for (const auto& label : bcLabels_) {
     bcTrees_[label] = nullptr;
     tkAlignTokens_[label] = esConsumes<Alignments, TrackerAlignmentRcd>(edm::ESInputTag{"", label});
   }
 
   for (const auto& label : bsLabels_) {
     bsTrees_[label] = nullptr;
     bsTokens_[label] = esConsumes<BeamSpotObjects, BeamSpotObjectsRcd>(edm::ESInputTag{"", label});
   }
 
   usesResource("TFileService");
 }

◆ ~PixelBaryCentreAnalyzer()

PixelBaryCentreAnalyzer::~PixelBaryCentreAnalyzer ( )

override

Definition at line 144 of file PixelBaryCentreAnalyzer.cc.

                                                   {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 188 of file PixelBaryCentreAnalyzer.cc.

References bcLabels_, bcTrees_, BPIX_, BPIX_DiffFlippedNonFlipped_, BPIX_Flipped_, BPIX_NonFlipped_, BPIXLayer_, BPIXLayer_DiffFlippedNonFlipped_, BPIXLayer_Flipped_, BPIXLayer_NonFlipped_, BS_, bsLabels_, bsTokens_, bsTrees_, edm::ESWatcher< T >::check(), align::DetectorGlobalPosition(), FPIX_, FPIX_minus_, FPIX_plus_, FPIXDisks_minus_, FPIXDisks_plus_, edm::EventSetup::getData(), gprToken_, PixelEndcapName::halfCylinder(), triggerObjects_cff::id, iEvent, initBC(), initBS(), SiPixelQuality::IsModuleBad(), TrackerGeometry::isThere(), label, PVValHelper::ladder, phase1PixelTopology::layer, ls_, Alignments::m_align, genParticles_cff::map, GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, phase_, PIX_, PixelSubdetector::PixelBarrel, GeomDetEnumerators::PixelBarrel, PixelSubdetector::PixelEndcap, GeomDetEnumerators::PixelEndcap, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), SiPixelPI::quadrant(), DetId::rawId(), relativeConstraints::ring, PixelEndcapName::ringName(), run_, siPixelQualityToken_, DetId::subdetId(), tkAlignTokens_, DetId::Tracker, trackerGeometryToken_, trackerTopologyToken_, AlignTransform::translation(), usePixelQuality_, watcherBS_, watcherTkAlign_, PV3DBase< T, PVType, FrameType >::x(), BeamSpotObjects::x(), PV3DBase< T, PVType, FrameType >::y(), BeamSpotObjects::y(), PV3DBase< T, PVType, FrameType >::z(), and BeamSpotObjects::z().

                                                                                          {
   bool prepareTkAlign = false;
   bool prepareBS = false;
 
   // ES watcher can noly run once in the same event,
   // otherwise it will turn false whatsoever because the condition doesn't change in the second time call.
   if (watcherTkAlign_.check(iSetup))
     prepareTkAlign = true;
   if (watcherBS_.check(iSetup))
     prepareBS = true;
 
   if (!prepareTkAlign && !prepareBS)
     return;
 
   run_ = iEvent.id().run();
   ls_ = iEvent.id().luminosityBlock();
 
   if (prepareTkAlign) {  // check for new IOV for TKAlign
 
     phase_ = -1;
 
     const TrackerGeometry* tkGeo = &iSetup.getData(trackerGeometryToken_);
     const TrackerTopology* tkTopo = &iSetup.getData(trackerTopologyToken_);
 
     if (tkGeo->isThere(GeomDetEnumerators::PixelBarrel) && tkGeo->isThere(GeomDetEnumerators::PixelEndcap))
       phase_ = 0;
     else if (tkGeo->isThere(GeomDetEnumerators::P1PXB) && tkGeo->isThere(GeomDetEnumerators::P1PXEC))
       phase_ = 1;
 
     // pixel quality
     const SiPixelQuality* badPixelInfo = &iSetup.getData(siPixelQualityToken_);
 
     // Tracker global position
     const Alignments* globalAlignments = &iSetup.getData(gprToken_);
     std::unique_ptr<const Alignments> globalPositions = std::make_unique<Alignments>(*globalAlignments);
     const AlignTransform& globalCoordinates = align::DetectorGlobalPosition(*globalPositions, DetId(DetId::Tracker));
     GlobalVector globalTkPosition(
         globalCoordinates.translation().x(), globalCoordinates.translation().y(), globalCoordinates.translation().z());
 
     // loop over bclabels
     for (const auto& label : bcLabels_) {
       // init tree content
       PixelBaryCentreAnalyzer::initBC();
 
       // Get TkAlign from EventSetup:
       const Alignments* alignments = &iSetup.getData(tkAlignTokens_[label]);
       std::vector<AlignTransform> tkAlignments = alignments->m_align;
 
       // PIX
       GlobalVector barycentre_PIX(0.0, 0.0, 0.0);
       // BPIX
       GlobalVector barycentre_BPIX(0.0, 0.0, 0.0);
       float nmodules_BPIX(0.);
       // FPIX
       GlobalVector barycentre_FPIX(0.0, 0.0, 0.0);
       float nmodules_FPIX(0.);
 
       // Per-layer/ladder barycentre for BPIX
       std::map<int, std::map<int, float>> nmodules_bpix;           // layer-ladder
       std::map<int, std::map<int, GlobalVector>> barycentre_bpix;  // layer-ladder
 
       // Per-disk/ring barycentre for FPIX
       std::map<int, std::map<int, float>> nmodules_fpix;           // disk-ring
       std::map<int, std::map<int, GlobalVector>> barycentre_fpix;  // disk-ring
 
       // Loop over tracker module
       for (const auto& ali : tkAlignments) {
         //DetId
         const DetId& detId = DetId(ali.rawId());
         // remove bad module
         if (usePixelQuality_ && badPixelInfo->IsModuleBad(detId))
           continue;
 
         // alignment for a given module
         GlobalVector ali_translation(ali.translation().x(), ali.translation().y(), ali.translation().z());
 
         int subid = DetId(detId).subdetId();
         // BPIX
         if (subid == PixelSubdetector::PixelBarrel) {
           nmodules_BPIX += 1;
           barycentre_BPIX += ali_translation;
           barycentre_PIX += ali_translation;
 
           int layer = tkTopo->pxbLayer(detId);
           int ladder = tkTopo->pxbLadder(detId);
           nmodules_bpix[layer][ladder] += 1;
           barycentre_bpix[layer][ladder] += ali_translation;
 
         }  // BPIX
 
         // FPIX
         if (subid == PixelSubdetector::PixelEndcap) {
           nmodules_FPIX += 1;
           barycentre_FPIX += ali_translation;
           barycentre_PIX += ali_translation;
 
           int disk = tkTopo->pxfDisk(detId);
           int quadrant = PixelEndcapName(detId, tkTopo, phase_).halfCylinder();
           if (quadrant < 3)
             disk *= -1;
 
           int ring = -9999;
           if (phase_ == 0) {
             ring = 1 + (tkTopo->pxfPanel(detId) + tkTopo->pxfModule(detId.rawId()) > 3);
           } else if (phase_ == 1) {
             ring = PixelEndcapName(detId, tkTopo, phase_).ringName();
           }
 
           nmodules_fpix[disk][ring] += 1;
           barycentre_fpix[disk][ring] += ali_translation;
 
         }  // FPIX
 
       }  // loop over tracker module
 
       //PIX
       float nmodules_PIX = nmodules_BPIX + nmodules_FPIX;
       barycentre_PIX *= (1.0 / nmodules_PIX);
       barycentre_PIX += globalTkPosition;
       PIX_ = GlobalPoint(barycentre_PIX.x(), barycentre_PIX.y(), barycentre_PIX.z());
 
       //BPIX
       barycentre_BPIX *= (1.0 / nmodules_BPIX);
       barycentre_BPIX += globalTkPosition;
       BPIX_ = GlobalPoint(barycentre_BPIX.x(), barycentre_BPIX.y(), barycentre_BPIX.z());
       //FPIX
       barycentre_FPIX *= (1.0 / nmodules_FPIX);
       barycentre_FPIX += globalTkPosition;
       FPIX_ = GlobalPoint(barycentre_FPIX.x(), barycentre_FPIX.y(), barycentre_FPIX.z());
 
       // Pixel substructures
 
       // BPix barycentre per-layer/per-ladder
       // assuming each ladder has the same number of modules in the same layer
       // inner =  flipped; outer = non-flipped
       //
       // Phase 0: Outer ladders are odd for layer 1,3 and even for layer 2
       // Phase 1: Outer ladders are odd for layer 4 and even for layer 1,2,3
       //
 
       int nmodules_BPIX_Flipped = 0;
       int nmodules_BPIX_NonFlipped = 0;
       GlobalVector BPIX_Flipped(0.0, 0.0, 0.0);
       GlobalVector BPIX_NonFlipped(0.0, 0.0, 0.0);
 
       // loop over layers
       for (std::map<int, std::map<int, GlobalVector>>::iterator il = barycentre_bpix.begin();
            il != barycentre_bpix.end();
            ++il) {
         int layer = il->first;
 
         int nmodulesLayer = 0;
         int nmodulesLayer_Flipped = 0;
         int nmodulesLayer_NonFlipped = 0;
         GlobalVector BPIXLayer(0.0, 0.0, 0.0);
         GlobalVector BPIXLayer_Flipped(0.0, 0.0, 0.0);
         GlobalVector BPIXLayer_NonFlipped(0.0, 0.0, 0.0);
 
         // loop over ladder
         std::map<int, GlobalVector> barycentreLayer = barycentre_bpix[layer];
         for (std::map<int, GlobalVector>::iterator it = barycentreLayer.begin(); it != barycentreLayer.end(); ++it) {
           int ladder = it->first;
           //BPIXLayerLadder_[layer][ladder] = (1.0/nmodules[layer][ladder])*barycentreLayer[ladder] + globalTkPosition;
 
           nmodulesLayer += nmodules_bpix[layer][ladder];
           BPIXLayer += barycentreLayer[ladder];
 
           // Phase-1
           //
           // Phase 1: Outer ladders are odd for layer 4 and even for layer 1,2,3
           if (phase_ == 1) {
             if (layer != 4) {  // layer 1-3
 
               if (ladder % 2 != 0) {  // odd ladder = inner = flipped
                 nmodulesLayer_Flipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_Flipped += barycentreLayer[ladder];
               } else {
                 nmodulesLayer_NonFlipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_NonFlipped += barycentreLayer[ladder];
               }
             } else {  // layer-4
 
               if (ladder % 2 == 0) {  // even ladder = inner = flipped
                 nmodulesLayer_Flipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_Flipped += barycentreLayer[ladder];
               } else {  // odd ladder = outer = non-flipped
                 nmodulesLayer_NonFlipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_NonFlipped += barycentreLayer[ladder];
               }
             }
 
           }  // phase-1
 
           // Phase-0
           //
           // Phase 0: Outer ladders are odd for layer 1,3 and even for layer 2
           if (phase_ == 0) {
             if (layer == 2) {  // layer-2
 
               if (ladder % 2 != 0) {  // odd ladder = inner = flipped
                 nmodulesLayer_Flipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_Flipped += barycentreLayer[ladder];
               } else {
                 nmodulesLayer_NonFlipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_NonFlipped += barycentreLayer[ladder];
               }
             } else {  // layer-1,3
 
               if (ladder % 2 == 0) {  // even ladder = inner = flipped
                 nmodulesLayer_Flipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_Flipped += barycentreLayer[ladder];
               } else {  // odd ladder = outer = non-flipped
                 nmodulesLayer_NonFlipped += nmodules_bpix[layer][ladder];
                 BPIXLayer_NonFlipped += barycentreLayer[ladder];
               }
             }
 
           }  // phase-0
 
         }  //loop over ladders
 
         // total BPIX flipped/non-flipped
         BPIX_Flipped += BPIXLayer_Flipped;
         BPIX_NonFlipped += BPIXLayer_NonFlipped;
         nmodules_BPIX_Flipped += nmodulesLayer_Flipped;
         nmodules_BPIX_NonFlipped += nmodulesLayer_NonFlipped;
 
         //BPIX per-layer
         BPIXLayer *= (1.0 / nmodulesLayer);
         BPIXLayer += globalTkPosition;
         BPIXLayer_Flipped *= (1.0 / nmodulesLayer_Flipped);
         BPIXLayer_Flipped += globalTkPosition;
         BPIXLayer_NonFlipped *= (1.0 / nmodulesLayer_NonFlipped);
         BPIXLayer_NonFlipped += globalTkPosition;
 
         BPIXLayer_[layer - 1] = GlobalPoint(BPIXLayer.x(), BPIXLayer.y(), BPIXLayer.z());
         BPIXLayer_Flipped_[layer - 1] =
             GlobalPoint(BPIXLayer_Flipped.x(), BPIXLayer_Flipped.y(), BPIXLayer_Flipped.z());
         BPIXLayer_NonFlipped_[layer - 1] =
             GlobalPoint(BPIXLayer_NonFlipped.x(), BPIXLayer_NonFlipped.y(), BPIXLayer_NonFlipped.z());
 
         BPIXLayer_DiffFlippedNonFlipped_[layer - 1] = GlobalPoint(BPIXLayer_Flipped.x() - BPIXLayer_NonFlipped.x(),
                                                                   BPIXLayer_Flipped.y() - BPIXLayer_NonFlipped.y(),
                                                                   BPIXLayer_Flipped.z() - BPIXLayer_NonFlipped.z());
 
       }  // loop over layers
 
       BPIX_Flipped *= (1.0 / nmodules_BPIX_Flipped);
       BPIX_Flipped += globalTkPosition;
       BPIX_Flipped_ = GlobalPoint(BPIX_Flipped.x(), BPIX_Flipped.y(), BPIX_Flipped.z());
       BPIX_NonFlipped *= (1.0 / nmodules_BPIX_NonFlipped);
       BPIX_NonFlipped += globalTkPosition;
       BPIX_NonFlipped_ = GlobalPoint(BPIX_NonFlipped.x(), BPIX_NonFlipped.y(), BPIX_NonFlipped.z());
       BPIX_DiffFlippedNonFlipped_ = GlobalPoint(BPIX_Flipped.x() - BPIX_NonFlipped.x(),
                                                 BPIX_Flipped.y() - BPIX_NonFlipped.y(),
                                                 BPIX_Flipped.z() - BPIX_NonFlipped.z());
 
       // FPIX substructures per-(signed)disk/per-ring
       int nmodules_FPIX_plus = 0;
       int nmodules_FPIX_minus = 0;
       GlobalVector FPIX_plus(0.0, 0.0, 0.0);
       GlobalVector FPIX_minus(0.0, 0.0, 0.0);
       // loop over disks
       for (std::map<int, std::map<int, GlobalVector>>::iterator id = barycentre_fpix.begin();
            id != barycentre_fpix.end();
            ++id) {
         int disk = id->first;
 
         int nmodulesDisk = 0;
         GlobalVector FPIXDisk(0.0, 0.0, 0.0);
 
         std::map<int, GlobalVector> baryCentreDisk = id->second;
         for (std::map<int, GlobalVector>::iterator ir = baryCentreDisk.begin(); ir != baryCentreDisk.end();
              ++ir) {  // loop over rings
           int ring = ir->first;
           nmodulesDisk += nmodules_fpix[disk][ring];
           FPIXDisk += ir->second;
           if (disk > 0) {
             nmodules_FPIX_plus += nmodules_fpix[disk][ring];
             FPIX_plus += ir->second;
           }
           if (disk < 0) {
             nmodules_FPIX_minus += nmodules_fpix[disk][ring];
             FPIX_minus += ir->second;
           }
 
         }  // loop over rings
 
         FPIXDisk *= (1.0 / nmodulesDisk);
         FPIXDisk += globalTkPosition;
 
         if (disk > 0)
           FPIXDisks_plus_[disk - 1] = GlobalPoint(FPIXDisk.x(), FPIXDisk.y(), FPIXDisk.z());
         if (disk < 0)
           FPIXDisks_minus_[-disk - 1] = GlobalPoint(FPIXDisk.x(), FPIXDisk.y(), FPIXDisk.z());
       }  // loop over disks
 
       FPIX_plus *= (1.0 / nmodules_FPIX_plus);
       FPIX_plus += globalTkPosition;
       FPIX_plus_ = GlobalPoint(FPIX_plus.x(), FPIX_plus.y(), FPIX_plus.z());
       FPIX_minus *= (1.0 / nmodules_FPIX_minus);
       FPIX_minus += globalTkPosition;
       FPIX_minus_ = GlobalPoint(FPIX_minus.x(), FPIX_minus.y(), FPIX_minus.z());
 
       bcTrees_[label]->Fill();
 
     }  // bcLabels_
 
   }  // check for new IOV for TKAlign
 
   // beamspot
   if (prepareBS) {
     // loop over bsLabels_
     for (const auto& label : bsLabels_) {
       // init bstree content
       PixelBaryCentreAnalyzer::initBS();
 
       // Get BeamSpot from EventSetup
       const BeamSpotObjects* mybeamspot = &iSetup.getData(bsTokens_[label]);
 
       BS_ = GlobalPoint(mybeamspot->x(), mybeamspot->y(), mybeamspot->z());
 
       bsTrees_[label]->Fill();
     }  // bsLabels_
 
   }  // check for new IOV for BS
 }

◆ beginJob()

void PixelBaryCentreAnalyzer::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 517 of file PixelBaryCentreAnalyzer.cc.

References bcLabels_, bcTrees_, BPIX_, BPIX_DiffFlippedNonFlipped_, BPIX_Flipped_, BPIX_NonFlipped_, BPIXLayer_, BPIXLayer_DiffFlippedNonFlipped_, BPIXLayer_Flipped_, BPIXLayer_NonFlipped_, BS_, bsLabels_, bsTrees_, FPIX_, FPIX_minus_, FPIX_plus_, FPIXDisks_minus_, FPIXDisks_plus_, mps_fire::i, label, phase1PixelTopology::layer, ls_, TFileService::make(), PIX_, run_, AlCaHLTBitMon_QueryRunRegistry::string, and tFileService.

                                        {
   // init bc bs trees
   for (const auto& label : bsLabels_) {
     std::string treeName = "BeamSpot";
     if (!label.empty())
       treeName = "BeamSpot_";
     treeName += label;
 
     bsTrees_[label] = tFileService->make<TTree>(TString(treeName), "PixelBarycentre analyzer ntuple");
 
     bsTrees_[label]->Branch("run", &run_, "run/I");
     bsTrees_[label]->Branch("ls", &ls_, "ls/I");
 
     bsTrees_[label]->Branch("BS", &BS_);
 
   }  // bsLabels_
 
   for (const auto& label : bcLabels_) {
     std::string treeName = "PixelBarycentre";
     if (!label.empty())
       treeName = "PixelBarycentre_";
     treeName += label;
     bcTrees_[label] = tFileService->make<TTree>(TString(treeName), "PixelBarycentre analyzer ntuple");
 
     bcTrees_[label]->Branch("run", &run_, "run/I");
     bcTrees_[label]->Branch("ls", &ls_, "ls/I");
 
     bcTrees_[label]->Branch("PIX", &PIX_);
 
     bcTrees_[label]->Branch("BPIX", &BPIX_);
     bcTrees_[label]->Branch("BPIX_Flipped", &BPIX_Flipped_);
     bcTrees_[label]->Branch("BPIX_NonFlipped", &BPIX_NonFlipped_);
     bcTrees_[label]->Branch("BPIX_DiffFlippedNonFlipped", &BPIX_DiffFlippedNonFlipped_);
 
     bcTrees_[label]->Branch("FPIX", &FPIX_);
     bcTrees_[label]->Branch("FPIX_plus", &FPIX_plus_);
     bcTrees_[label]->Branch("FPIX_minus", &FPIX_minus_);
 
     //per-layer
     for (unsigned int i = 0; i < 4; i++) {
       TString structure = "BPIXLYR";
       int layer = i + 1;
       structure += layer;
 
       bcTrees_[label]->Branch(structure, &BPIXLayer_[i]);
       bcTrees_[label]->Branch(structure + "_Flipped", &BPIXLayer_Flipped_[i]);
       bcTrees_[label]->Branch(structure + "_NonFlipped", &BPIXLayer_NonFlipped_[i]);
       bcTrees_[label]->Branch(structure + "_DiffFlippedNonFlipped", &BPIXLayer_DiffFlippedNonFlipped_[i]);
     }
 
     //per-disk/ring
     for (unsigned int i = 0; i < 3; i++) {
       TString structure = "FPIXDisk+";
       int disk = i + 1;
       structure += disk;
       bcTrees_[label]->Branch(structure, &FPIXDisks_plus_[i]);
 
       structure = "FPIXDisk-";
       structure += disk;
       bcTrees_[label]->Branch(structure, &FPIXDisks_minus_[i]);
     }
 
   }  // bcLabels_
 }

◆ endJob()

void PixelBaryCentreAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 583 of file PixelBaryCentreAnalyzer.cc.

References bcLabels_, bcTrees_, bsLabels_, and bsTrees_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                      {
   bcLabels_.clear();
   bsLabels_.clear();
 
   bcTrees_.clear();
   bsTrees_.clear();
 }

◆ fillDescriptions()

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

static

Definition at line 592 of file PixelBaryCentreAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                          {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

◆ initBC()

void PixelBaryCentreAnalyzer::initBC ( )

private

Definition at line 159 of file PixelBaryCentreAnalyzer.cc.

References BPIX_, BPIX_DiffFlippedNonFlipped_, BPIX_Flipped_, BPIX_NonFlipped_, BPIXLayer_, BPIXLayer_DiffFlippedNonFlipped_, BPIXLayer_Flipped_, BPIXLayer_NonFlipped_, FPIX_, FPIX_minus_, FPIX_plus_, FPIXDisks_minus_, FPIXDisks_plus_, mps_fire::i, and PIX_.

Referenced by analyze().

                                      {
   // init to large number (unreasonable number) not zero
   double dummy_float = 999999.0;
 
   PIX_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
   BPIX_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
   FPIX_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
 
   BPIX_Flipped_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
   BPIX_NonFlipped_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
   BPIX_DiffFlippedNonFlipped_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
 
   FPIX_plus_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
   FPIX_minus_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
 
   for (unsigned int i = 0; i < 4; i++) {
     BPIXLayer_[i] = GlobalPoint(dummy_float, dummy_float, dummy_float);
     BPIXLayer_Flipped_[i] = GlobalPoint(dummy_float, dummy_float, dummy_float);
     BPIXLayer_NonFlipped_[i] = GlobalPoint(dummy_float, dummy_float, dummy_float);
     BPIXLayer_DiffFlippedNonFlipped_[i] = GlobalPoint(dummy_float, dummy_float, dummy_float);
   }
 
   for (unsigned int i = 0; i < 3; i++) {
     FPIXDisks_plus_[i] = GlobalPoint(dummy_float, dummy_float, dummy_float);
     FPIXDisks_minus_[i] = GlobalPoint(dummy_float, dummy_float, dummy_float);
   }
 }

◆ initBS()

void PixelBaryCentreAnalyzer::initBS ( )

private

Definition at line 153 of file PixelBaryCentreAnalyzer.cc.

References BS_.

Referenced by analyze().

                                      {
   double dummy_float = 999999.0;
 
   BS_ = GlobalPoint(dummy_float, dummy_float, dummy_float);
 }

Member Data Documentation

◆ bcLabels_

std::vector<std::string> PixelBaryCentreAnalyzer::bcLabels_

private

Definition at line 85 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), endJob(), and PixelBaryCentreAnalyzer().

◆ bcTrees_

std::map<std::string, TTree*> PixelBaryCentreAnalyzer::bcTrees_

private

Definition at line 116 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), endJob(), and PixelBaryCentreAnalyzer().

◆ BPIX_

GlobalPoint PixelBaryCentreAnalyzer::BPIX_

private

Definition at line 103 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BPIX_DiffFlippedNonFlipped_

GlobalPoint PixelBaryCentreAnalyzer::BPIX_DiffFlippedNonFlipped_

private

Definition at line 104 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BPIX_Flipped_

GlobalPoint PixelBaryCentreAnalyzer::BPIX_Flipped_

private

Definition at line 104 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BPIX_NonFlipped_

GlobalPoint PixelBaryCentreAnalyzer::BPIX_NonFlipped_

private

Definition at line 104 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BPIXLayer_

GlobalPoint PixelBaryCentreAnalyzer::BPIXLayer_[4]

private

Definition at line 106 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BPIXLayer_DiffFlippedNonFlipped_

GlobalPoint PixelBaryCentreAnalyzer::BPIXLayer_DiffFlippedNonFlipped_[4]

private

Definition at line 109 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BPIXLayer_Flipped_

GlobalPoint PixelBaryCentreAnalyzer::BPIXLayer_Flipped_[4]

private

Definition at line 107 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BPIXLayer_NonFlipped_

GlobalPoint PixelBaryCentreAnalyzer::BPIXLayer_NonFlipped_[4]

private

Definition at line 108 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ BS_

GlobalPoint PixelBaryCentreAnalyzer::BS_

private

Definition at line 101 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBS().

◆ bsLabels_

std::vector<std::string> PixelBaryCentreAnalyzer::bsLabels_

private

Definition at line 87 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), endJob(), and PixelBaryCentreAnalyzer().

◆ bsTokens_

std::map<std::string, edm::ESGetToken<BeamSpotObjects, BeamSpotObjectsRcd> > PixelBaryCentreAnalyzer::bsTokens_

private

Definition at line 95 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), and PixelBaryCentreAnalyzer().

◆ bsTrees_

std::map<std::string, TTree*> PixelBaryCentreAnalyzer::bsTrees_

private

Definition at line 117 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), endJob(), and PixelBaryCentreAnalyzer().

◆ FPIX_

GlobalPoint PixelBaryCentreAnalyzer::FPIX_

private

Definition at line 103 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ FPIX_minus_

GlobalPoint PixelBaryCentreAnalyzer::FPIX_minus_

private

Definition at line 111 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ FPIX_plus_

GlobalPoint PixelBaryCentreAnalyzer::FPIX_plus_

private

Definition at line 111 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ FPIXDisks_minus_

GlobalPoint PixelBaryCentreAnalyzer::FPIXDisks_minus_[3]

private

Definition at line 113 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ FPIXDisks_plus_

GlobalPoint PixelBaryCentreAnalyzer::FPIXDisks_plus_[3]

private

Definition at line 112 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ gprToken_

edm::ESGetToken<Alignments, GlobalPositionRcd> PixelBaryCentreAnalyzer::gprToken_

private

Definition at line 93 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

◆ ls_

int PixelBaryCentreAnalyzer::ls_

private

Definition at line 99 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ phase_

int PixelBaryCentreAnalyzer::phase_

private

Definition at line 78 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

◆ PIX_

GlobalPoint PixelBaryCentreAnalyzer::PIX_

private

Definition at line 103 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), beginJob(), and initBC().

◆ run_

int PixelBaryCentreAnalyzer::run_

private

Definition at line 98 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ siPixelQualityToken_

edm::ESGetToken<SiPixelQuality, SiPixelQualityFromDbRcd> PixelBaryCentreAnalyzer::siPixelQualityToken_

private

Definition at line 91 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

◆ tFileService

edm::Service<TFileService> PixelBaryCentreAnalyzer::tFileService

private

Definition at line 115 of file PixelBaryCentreAnalyzer.cc.

Referenced by beginJob().

◆ tkAlignTokens_

std::map<std::string, edm::ESGetToken<Alignments, TrackerAlignmentRcd> > PixelBaryCentreAnalyzer::tkAlignTokens_

private

Definition at line 94 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze(), and PixelBaryCentreAnalyzer().

◆ trackerGeometryToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> PixelBaryCentreAnalyzer::trackerGeometryToken_

private

Definition at line 89 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

◆ trackerTopologyToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> PixelBaryCentreAnalyzer::trackerTopologyToken_

private

Definition at line 90 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

◆ usePixelQuality_

bool PixelBaryCentreAnalyzer::usePixelQuality_

private

Definition at line 77 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

◆ watcherBS_

edm::ESWatcher<BeamSpotObjectsRcd> PixelBaryCentreAnalyzer::watcherBS_

private

Definition at line 81 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

◆ watcherTkAlign_

edm::ESWatcher<TrackerAlignmentRcd> PixelBaryCentreAnalyzer::watcherTkAlign_

private

Definition at line 82 of file PixelBaryCentreAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PixelBaryCentreAnalyzer()

◆ ~PixelBaryCentreAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

◆ initBC()

◆ initBS()

Member Data Documentation

◆ bcLabels_

◆ bcTrees_

◆ BPIX_

◆ BPIX_DiffFlippedNonFlipped_

◆ BPIX_Flipped_

◆ BPIX_NonFlipped_

◆ BPIXLayer_

◆ BPIXLayer_DiffFlippedNonFlipped_

◆ BPIXLayer_Flipped_

◆ BPIXLayer_NonFlipped_

◆ BS_

◆ bsLabels_

◆ bsTokens_

◆ bsTrees_

◆ FPIX_

◆ FPIX_minus_

◆ FPIX_plus_

◆ FPIXDisks_minus_

◆ FPIXDisks_plus_

◆ gprToken_

◆ ls_

◆ phase_

◆ PIX_

◆ run_

◆ siPixelQualityToken_

◆ tFileService

◆ tkAlignTokens_

◆ trackerGeometryToken_

◆ trackerTopologyToken_

◆ usePixelQuality_

◆ watcherBS_

◆ watcherTkAlign_