#include <EJTerm/SiPixelPhase1Analyzer/plugins/SiPixelPhase1Analyzer.cc>

Inheritance diagram for SiPixelPhase1Analyzer:

Classes
struct	complementaryElements

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

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

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

void	BookBarrelBins (const TrackerGeometry &theTrackerGeometry, const TrackerTopology *tt)

void	BookBarrelHistograms (TDirectory *currentDir, const std::string &currentHistoName)

void	BookBins (const TrackerGeometry &theTrackerGeometry, const TrackerTopology *tt)

void	BookForwardBins (const TrackerGeometry &theTrackerGeometry, const TrackerTopology *tt)

void	BookForwardHistograms (TDirectory *currentDir, const std::string &currentHistoName)

void	BookHistograms ()

void	FillBarrelBinsAnalyze (const TrackerGeometry &theTrackerGeometry, const TrackerTopology *tt, unsigned rawId, const GlobalPoint &globalPoint)

void	FillBarrelBinsRemap (const TrackerGeometry &theTrackerGeometry, const TrackerTopology *tt)

void	FillBins (edm::Handle< reco::TrackCollection > tracks, const TrackerGeometry &theTrackerGeometry, const TrackerTopology tt)

void	FillForwardBinsAnalyze (const TrackerGeometry &theTrackerGeometry, const TrackerTopology *tt, unsigned rawId, const GlobalPoint &globalPoint)

void	FillForwardBinsRemap (const TrackerGeometry &theTrackerGeometry, const TrackerTopology *tt)

void	SaveDetectorData (bool isBarrel, unsigned rawId, int shell_hc, int layer_disk, int ladder_blade)

void	SaveDetectorVertices (const TrackerTopology *tt)

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

std::vector< std::string >	baseHistogramName

std::map< uint32_t, TGraph * >	bins

std::map< uint32_t, TGraph * >	binsSummary

std::ofstream	debugFile

std::string	debugFileName

bool	firstEvent

edm::Service< TFileService >	fs

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

std::vector< unsigned >	isBarrelSource

std::map< unsigned short, complementaryElements >	mapOfComplementaryElements

OperationMode	opMode

mat4	orthoProjectionMatrix

std::vector< std::string >	pathToHistograms

std::map< std::string, std::vector< TH2Poly * > >	pxfTh2PolyForward

std::map< std::string, TH2Poly * >	pxfTh2PolyForwardSummary

TFile *	rootFileHandle

std::map< std::string, std::vector< TH2Poly * > >	th2PolyBarrel

std::map< std::string, TH2Poly * >	th2PolyBarrelSummary

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

edm::EDGetTokenT< reco::TrackCollection >	tracksToken

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 70 of file SiPixelPhase1Analyzer.h.

Constructor & Destructor Documentation

◆ SiPixelPhase1Analyzer()

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

explicit

Definition at line 12 of file SiPixelPhase1Analyzer.cc.

References analazedRootFileName, baseHistogramName, mat4::BuildOrthographicMatrix(), debugFile, debugFileName, f, edm::ParameterSet::getParameter(), MODE_ANALYZE, MODE_REMAP, opMode, orthoProjectionMatrix, MillePedeFileConverter_cfg::out, pathToHistograms, and tracksToken.

     : geomToken_(esConsumes()),
       topoToken_(esConsumes()),
       opMode(static_cast<OperationMode>(iConfig.getUntrackedParameter<unsigned int>("opMode"))),
       debugFileName(iConfig.getUntrackedParameter<string>("debugFileName")),
       firstEvent(true),
       rootFileHandle(nullptr),
       isBarrelSource(iConfig.getUntrackedParameter<vector<unsigned>>("isBarrelSource")),
       analazedRootFileName(iConfig.getUntrackedParameter<vector<string>>("remapRootFileName")),
       pathToHistograms(iConfig.getUntrackedParameter<vector<string>>("pathToHistograms")),
       baseHistogramName(iConfig.getUntrackedParameter<vector<string>>("baseHistogramName")) {
 #ifdef DEBUG_MODE
   debugFile = std::ofstream(debugFileName.c_str(), std::ofstream::out);
 #endif
   usesResource("TFileService");
 
   orthoProjectionMatrix.BuildOrthographicMatrix(1.0f, -1.0f, 1.0f, -1.0f, -10.0f, 10.0f);
 
   switch (opMode) {
     case MODE_ANALYZE:
 
       tracksToken = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("src"));
 
       analazedRootFileName.clear();
 
       pathToHistograms.clear();
       pathToHistograms.push_back("RecHits/");
 
       baseHistogramName.clear();
       baseHistogramName.push_back("RecHits");
 
       break;
     case MODE_REMAP:
       break;
     default:
       break;
   }
 }

◆ ~SiPixelPhase1Analyzer()

SiPixelPhase1Analyzer::~SiPixelPhase1Analyzer ( )

override

Definition at line 51 of file SiPixelPhase1Analyzer.cc.

References bins, binsSummary, debugFile, and mps_fire::i.

                                               {
   for (auto& i : bins) {
     delete i.second;
   }
 
   for (auto& i : binsSummary) {
     delete i.second;
   }
 
 #ifdef DEBUG_MODE
   debugFile.close();
 #endif
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 66 of file SiPixelPhase1Analyzer.cc.

References BookBins(), BookHistograms(), FillBins(), firstEvent, geomToken_, edm::EventSetup::getData(), iEvent, MODE_ANALYZE, MODE_REMAP, opMode, topoToken_, DiMuonV_cfg::tracks, and tracksToken.

                                                                                        {
   const auto& theTrackerGeometry = iSetup.getData(geomToken_);
   const auto& tt = &iSetup.getData(topoToken_);
 
   if (firstEvent) {
     BookHistograms();
     BookBins(theTrackerGeometry, tt);
     if (opMode == MODE_REMAP) {
       FillBins(nullptr, theTrackerGeometry, tt);
     }
     firstEvent = false;
   }
   if (opMode == MODE_ANALYZE) {
     edm::Handle<reco::TrackCollection> tracks;
     iEvent.getByToken(tracksToken, tracks);
     if (!tracks.isValid()) {
       LogInfo("Analyzer") << "reco::TrackCollection not found... Aborting...\n";
       return;
     }
     FillBins(&tracks, theTrackerGeometry, tt);
   }
   // debugFile << "SiPixelPhase1Analyzer::analyze() - Event " << iEvent.run() << "/" << iEvent.id().event() << endl;
 }

◆ BookBarrelBins()

void SiPixelPhase1Analyzer::BookBarrelBins	(	const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt
	)

private

Definition at line 268 of file SiPixelPhase1Analyzer.cc.

References b, baseHistogramName, bins, binsSummary, Surface::bounds(), debugFile, TrackerGeometry::detsPXB(), f, ALCARECOEcalPhiSym_cff::float, GeomDet::geographicalId(), mps_fire::i, isBarrelSource, isotrackApplyRegressor::k, PVValHelper::ladder, PixelBarrelName::ladderName(), nano_mu_digi_cff::layer, PixelBarrelName::layerName(), MODE_ANALYZE, opMode, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, nano_mu_digi_cff::rawId, SaveDetectorData(), PixelBarrelName::shell(), mathSSE::sqrt(), GeomDet::surface(), th2PolyBarrel, th2PolyBarrelSummary, and Surface::toGlobal().

Referenced by BookBins().

                                                                                                                {
   TrackingGeometry::DetContainer pxb = theTrackerGeometry.detsPXB();
 #ifdef DEBUG_MODE
   debugFile << "There are " << pxb.size() << " detector elements in the PXB." << endl;
 #endif
   for (auto& i : pxb) {
     const GeomDet* det = i;
 
     PXBDetId id = det->geographicalId();
 
     Local2DPoint origin;
     GlobalPoint p = det->surface().toGlobal(origin);
 
     int layer = tt->pxbLayer(id);
     int ladder = tt->pxbLadder(id);
 
 #ifdef DEBUG_MODE
     int module = tt->pxbModule(id);
     PixelBarrelName pixelBarrelName(id, tt, true);
     SaveDetectorData(
         true, id.rawId(), pixelBarrelName.shell(), pixelBarrelName.layerName(), pixelBarrelName.ladderName());
 #endif
 
 #ifdef DEBUG_MODE
     float r = sqrt(p.x() * p.x() + p.y() * p.y());
 
     debugFile << "Layer: " << layer << "\tLadder: " << ladder << "\tModule: " << module << "\t(x, y, z, r2): (" << p.x()
               << ", " << p.y() << ", " << p.z() << ", " << r << ")" << endl;
 #endif
 
     const Bounds& b = (det->surface().bounds());
     float bl = b.length();
 
 #ifdef DEBUG_MODE
     float bw = b.width();
     float bt = b.thickness();
 
     debugFile << "Length: " << bl << "\tWidth: " << bw << "\tThickness: " << bt << endl;
 #endif
 
     float vertX[] = {p.z() - bl * 0.5f, p.z() + bl * 0.5f, p.z() + bl * 0.5f, p.z() - bl * 0.5f, p.z() - bl * 0.5f};
     float vertY[] = {(ladder - 1.0f), (ladder - 1.0f), (float)ladder, (float)ladder, (ladder - 1.0f)};
 
     bins[id.rawId()] = new TGraph(5, vertX, vertY);
     bins[id.rawId()]->SetName(TString::Format("%u", id.rawId()));
 
     // Summary plot
     for (unsigned k = 0; k < 5; ++k) {
       vertX[k] += ((layer == 2 || layer == 3) ? 0.0f : -60.0f);
       vertY[k] += ((layer > 2) ? 30.0f : 0.0f);
     }
 
     binsSummary[id.rawId()] = new TGraph(5, vertX, vertY);
     binsSummary[id.rawId()]->SetName(TString::Format("%u", id.rawId()));
 
     for (unsigned nameNum = 0; nameNum < baseHistogramName.size(); ++nameNum) {
       if (isBarrelSource[nameNum] || opMode == MODE_ANALYZE) {
         const string& strName = baseHistogramName[nameNum];
         th2PolyBarrel[strName][layer - 1]->AddBin(bins[id.rawId()]->Clone());
         th2PolyBarrelSummary[strName]->AddBin(binsSummary[id.rawId()]->Clone());
       }
     }
   }
 }

◆ BookBarrelHistograms()

void SiPixelPhase1Analyzer::BookBarrelHistograms	(	TDirectory *	currentDir,
		const std::string &	currentHistoName
	)

private

Definition at line 128 of file SiPixelPhase1Analyzer.cc.

References runEdmFileComparison::currentDir, f, mps_fire::i, MODE_ANALYZE, opMode, th2PolyBarrel, th2PolyBarrelSummary, and to_string().

Referenced by BookHistograms().

                                                                                                        {
   string histName;
   TH2Poly* th2p;
 
 #ifdef DEBUG_MODE
   TH2* th2;
 #endif
 
   for (unsigned i = 0; i < 4; ++i) {
     histName = "barrel_layer_";
 
     th2p = new TH2Poly((histName + std::to_string(i + 1)).c_str(), "PXBMap", -15.0, 15.0, 0.0, 5.0);
 
     th2p->SetFloat();
 
     th2p->GetXaxis()->SetTitle("z [cm]");
     th2p->GetYaxis()->SetTitle("ladder");
 
 #ifdef DEBUG_MODE
     th2p->SetOption("COLZ 0 TEXT");
 #else
     th2p->SetOption("COLZ L");
 #endif
 
     currentDir->Add(th2p);
     th2PolyBarrel[currentHistoName].push_back(th2p);
 
 #ifdef DEBUG_MODE
     if (opMode == MODE_ANALYZE) {
       th2 = new TH2I((histName + std::to_string(i + 1) + "_DEBUG").c_str(),
                      "position",
                      3000,
                      -30.0f,
                      30.0f,
                      1000,
                      -4.5f * (i + 1),
                      4.5f * (i + 1));
 
       th2->GetXaxis()->SetTitle("z [cm]");
       th2->GetYaxis()->SetTitle("-x [?]");
 
       th2->SetOption("COLZ 0 TEXT");
 
       currentDir->Add(th2);
       th2PolyBarrelDebug[currentHistoName].push_back(th2);
     }
 #endif
   }
 
   th2p = new TH2Poly("barrel_summary", "PXBMap", -5.0, 5.0, 0.0, 5.0);
   th2p->SetFloat();
 
   th2p->GetXaxis()->SetTitle("");
   th2p->GetYaxis()->SetTitle("~ladder");
 
   th2p->SetOption("COLZ L");
 
   currentDir->Add(th2p);
   th2PolyBarrelSummary[currentHistoName] = th2p;
 }

◆ BookBins()

void SiPixelPhase1Analyzer::BookBins	(	const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt
	)

private

Definition at line 259 of file SiPixelPhase1Analyzer.cc.

References BookBarrelBins(), BookForwardBins(), and SaveDetectorVertices().

Referenced by analyze().

                                                                                                          {
   BookBarrelBins(theTrackerGeometry, tt);
   BookForwardBins(theTrackerGeometry, tt);
 
 #ifdef DEBUG_MODE
   SaveDetectorVertices(tt);
 #endif
 }

◆ BookForwardBins()

void SiPixelPhase1Analyzer::BookForwardBins	(	const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt
	)

private

Definition at line 333 of file SiPixelPhase1Analyzer.cc.

References b, baseHistogramName, bins, binsSummary, PixelEndcapName::bladeName(), Surface::bounds(), CODE_FORWARD, debugFile, TrackerGeometry::detsPXF(), PixelEndcapName::diskName(), f, ALCARECOEcalPhiSym_cff::float, GeomDet::geographicalId(), PixelEndcapName::halfCylinder(), mps_fire::i, isBarrelSource, dqmiolumiharvest::j, isotrackApplyRegressor::k, mapOfComplementaryElements, MODE_ANALYZE, mat4::MulVec(), opMode, orthoProjectionMatrix, AlCaHLTBitMon_ParallelJobs::p, pxfTh2PolyForward, pxfTh2PolyForwardSummary, alignCSCRings::r, nano_mu_digi_cff::rawId, makeMuonMisalignmentScenario::rot, GeomDet::rotation(), SaveDetectorData(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), mathSSE::sqrt(), GeomDet::surface(), edm::swap(), and Surface::toGlobal().

Referenced by BookBins().

                                                                                                                 {
   TrackingGeometry::DetContainer pxf = theTrackerGeometry.detsPXF();
 #ifdef DEBUG_MODE
   debugFile << "There are " << pxf.size() << " detector elements in the PXF." << endl;
 #endif
   bool firstForwardElem = true;
 
   float elemWidth = 1.0f;
   float elemLength = 1.0f;
 
   // FIRST PASS -> MAKE MAP OF CORRESPONDING ELEMENTS (BLADES ON BOTH PANELS)
   for (auto& i : pxf) {
     const GeomDet* det = i;
 
     PXFDetId id = det->geographicalId();
 
     Local2DPoint origin;
     GlobalPoint p = det->surface().toGlobal(origin);
 
     int panel = tt->pxfPanel(id);
     int side = tt->side(id);   //tt->pxfSide(id);
     int disk = tt->layer(id);  //tt->pxfDisk(id);
     int blade = tt->pxfBlade(id);
 
 #ifdef DEBUG_MODE
     int module = tt->module(id);  //tt->pxfModule(id);
     PixelEndcapName pixelEndcapName(id, tt, true);
     SaveDetectorData(
         false, id.rawId(), pixelEndcapName.halfCylinder(), pixelEndcapName.diskName(), pixelEndcapName.bladeName());
 #endif
 
 #ifdef DEBUG_MODE
     float r = sqrt(p.x() * p.x() + p.y() * p.y());
 
     debugFile << "Panel: " << panel << "\tSide: " << side << "\tDisk: " << disk << "\tBlade: " << blade
               << "\tModule: " << module << "\t(x, y, z, r): (" << p.x() << ", " << p.y() << ", " << p.z() << ", " << r
               << ")" << endl;
 #endif
     if (firstForwardElem) {
       const Bounds& b = det->surface().bounds();  //const RectangularPlaneBounds& b
 
       elemLength = b.length();
       elemWidth = b.width();
 
       firstForwardElem = false;
     }
 
     const auto& rot = det->rotation();
 
     mat4 transMat(
         rot.xx(), rot.xy(), rot.xz(), rot.yx(), rot.yy(), rot.yz(), rot.zx(), rot.zy(), rot.zz(), p.x(), p.y(), p.z());
 
     mapOfComplementaryElements[CODE_FORWARD(side, disk, blade)].mat[panel - 1] = transMat;
     mapOfComplementaryElements[CODE_FORWARD(side, disk, blade)].rawId[panel - 1] = id.rawId();
   }
 
   // SECOND PASS -> USE INFORMATION FROM MAP TO MAKE GEOMETRIC OBJECTS (BINS)
   for (auto& i : mapOfComplementaryElements) {
     // decode side&disk from the map key
     unsigned side = i.first & 0xF;
     unsigned disk = (i.first & 0xF0) >> 4;
     // unsigned blade = (i.first & 0xFF00) >> 8;
 
     unsigned mapIdx = disk + (side - 1) * 3 - 1;
 
     // normal vectors of elements point to the (almost) opposite direction, so correction is needed before interploation (probably not 100% correct but fast)
     i.second.mat[1].data[0] = -i.second.mat[1].data[0];
     i.second.mat[1].data[1] = -i.second.mat[1].data[1];
     i.second.mat[1].data[2] = -i.second.mat[1].data[2];
 
     i.second.mat[1].data[6] = -i.second.mat[1].data[6];
     i.second.mat[1].data[7] = -i.second.mat[1].data[7];
     i.second.mat[1].data[8] = -i.second.mat[1].data[8];
 
     mat4 meanTransform = (i.second.mat[0] + i.second.mat[1]) * 0.5f;
     // mat4 meanTransform = i.second.mat[0];
 
     static const float baseVertX[4] = {-elemWidth * 0.8f, -elemWidth * 0.5f, elemWidth * 0.8f, elemWidth * 0.5f};
     static const float baseVertY[4] = {
         elemLength * 0.38f, -elemLength * 0.38f, elemLength * 0.38f, -elemLength * 0.38f};
 
     float vertXPanel[2][4], vertYPanel[2][4];
     float vertIn[3], vertOut[3];
 
     /*              
         (1)  __________________ (3)
              \             /
               \           /
                \             /
                 \       /
               (2)\_________/(4)
              
              - division line: (2) - (3)
              
              - panel 1: triangle of lower area (2, 3, 4)
              - panel 2: triangle of bigger area (1, 2, 3)
     */
 
     // obtain transformed vertices
     for (unsigned j = 0; j < 4; ++j) {
       vertIn[0] = baseVertX[j];
       vertIn[1] = baseVertY[j];
       vertIn[2] = 0.0f;
 
       meanTransform.MulVec(vertIn, vertOut);
       std::swap(vertIn, vertOut);
       orthoProjectionMatrix.MulVec(vertIn, vertOut);
 
       // vertical flip
       vertOut[0] = -vertOut[0];  // so that inner elements have positive x-coordinate
 
       if (j > 0) {
         vertXPanel[0][j - 1] = vertOut[0];
         vertYPanel[0][j - 1] = vertOut[1];  // for panel 2
       }
       if (j < 3) {
         vertXPanel[1][j] = vertOut[0];
         vertYPanel[1][j] = vertOut[1];  // for panel 1
       }
     }
 
     for (unsigned j = 0; j < 2; ++j) {
       vertXPanel[j][3] = vertXPanel[j][0];
       vertYPanel[j][3] = vertYPanel[j][0];
 
       bins[i.second.rawId[j]] = new TGraph(4, vertXPanel[j], vertYPanel[j]);
       bins[i.second.rawId[j]]->SetName(TString::Format("%u", i.second.rawId[j]));
 
       // for (auto strName: baseHistogramName)
       // {
       // pxfTh2PolyForward[strName][mapIdx]->AddBin(bins[i.second.rawId[j]]->Clone());
       // }
 
       // Summary plot
       for (unsigned k = 0; k < 4; ++k) {
         vertXPanel[j][k] += (float(side) - 1.5f) * 40.0f;
         vertYPanel[j][k] += (disk - 1) * 35.0f;
       }
 
       binsSummary[i.second.rawId[j]] = new TGraph(4, vertXPanel[j], vertYPanel[j]);
       binsSummary[i.second.rawId[j]]->SetName(TString::Format("%u", i.second.rawId[j]));
 
       for (unsigned nameNum = 0; nameNum < baseHistogramName.size(); ++nameNum) {
         if (!isBarrelSource[nameNum] || opMode == MODE_ANALYZE) {
           const string& strName = baseHistogramName[nameNum];
           pxfTh2PolyForward[strName][mapIdx]->AddBin(bins[i.second.rawId[j]]->Clone());
           pxfTh2PolyForwardSummary[strName]->AddBin(binsSummary[i.second.rawId[j]]->Clone());
         }
       }
     }
   }
 }

◆ BookForwardHistograms()

void SiPixelPhase1Analyzer::BookForwardHistograms	(	TDirectory *	currentDir,
		const std::string &	currentHistoName
	)

private

Definition at line 189 of file SiPixelPhase1Analyzer.cc.

References runEdmFileComparison::currentDir, f, MODE_ANALYZE, opMode, pxfTh2PolyForward, pxfTh2PolyForwardSummary, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), and to_string().

Referenced by BookHistograms().

                                                                                                         {
   string histName;
   TH2Poly* th2p;
 #ifdef DEBUG_MODE
   TH2* th2;
 #endif
 
   for (unsigned side = 1; side <= 2; ++side) {
     for (unsigned disk = 1; disk <= 3; ++disk) {
       histName = "forward_disk_";
 
       th2p = new TH2Poly((histName + std::to_string((side == 1 ? -(int(disk)) : (int)disk))).c_str(),
                          "PXFMap",
                          -15.0,
                          15.0,
                          -15.0,
                          15.0);
 
       th2p->SetFloat();
 
       th2p->GetXaxis()->SetTitle("x [cm]");
       th2p->GetYaxis()->SetTitle("y [cm]");
 
 #ifdef DEBUG_MODE
       th2p->SetOption("COLZ 0 TEXT");
 #else
       th2p->SetOption("COLZ L");
 #endif
       currentDir->Add(th2p);
       pxfTh2PolyForward[currentHistoName].push_back(th2p);
 
 #ifdef DEBUG_MODE
       if (opMode == MODE_ANALYZE) {
         th2 = new TH2I((histName + std::to_string((side == 1 ? -(int(disk)) : (int)disk)) + "_DEBUG").c_str(),
                        "position",
                        1000,
                        -15.0f,
                        15.0f,
                        1000,
                        -15.0f,
                        15.0f);
 
         th2->GetXaxis()->SetTitle("x [cm]");
         th2->GetYaxis()->SetTitle("y [cm]");
 
         th2->SetOption("COLZ 0 TEXT");
 
         currentDir->Add(th2);
         pxfTh2PolyForwardDebug[currentHistoName].push_back(th2);
       }
 #endif
     }
   }
 
   th2p = new TH2Poly("forward_summary", "PXFMap", -40.0, 50.0, -20.0, 90.0);
   th2p->SetFloat();
 
   th2p->GetXaxis()->SetTitle("");
   th2p->GetYaxis()->SetTitle("");
 
 #ifdef DEBUG_MODE
   th2p->SetOption("COLZ 0 TEXT");
 #else
   th2p->SetOption("COLZ L");
 #endif
 
   currentDir->Add(th2p);
   pxfTh2PolyForwardSummary[currentHistoName] = th2p;
 }

◆ BookHistograms()

void SiPixelPhase1Analyzer::BookHistograms ( )

private

Definition at line 94 of file SiPixelPhase1Analyzer.cc.

References baseHistogramName, BookBarrelHistograms(), BookForwardHistograms(), runEdmFileComparison::currentDir, debugFile, fs, TFileService::fullPath(), TFileService::getBareDirectory(), isBarrelSource, dqmiolumiharvest::j, TFileService::make(), MODE_REMAP, and opMode.

Referenced by analyze().

                                            {
   // ghost object <initializes> FileService (without it, it crashes and creation of directories would not be possible)
   TObject* ghostObj = fs->make<TH2Poly>("ghost", "ghost", -1, 1, -1, 1);
 
   TDirectory* topDir = fs->getBareDirectory();
   topDir->cd();
 
 #ifdef DEBUG_MODE
   debugFile << "Full path: " << fs->fullPath() << endl << endl;
 #endif
   string histName;
   for (unsigned j = 0; j < baseHistogramName.size(); ++j) {
     string currentHistoName = baseHistogramName[j];
 
     TDirectory* currentDir = topDir->mkdir(currentHistoName.c_str());
     currentDir->cd();
 
     if (opMode == MODE_REMAP) {
       if (isBarrelSource[j]) {
         BookBarrelHistograms(currentDir, currentHistoName);
       } else {
         BookForwardHistograms(currentDir, currentHistoName);
       }
     } else {
       BookBarrelHistograms(currentDir, currentHistoName);
       BookForwardHistograms(currentDir, currentHistoName);
     }
 
     topDir->cd();
   }
 
   ghostObj->Delete();  //not needed any more
 }

◆ FillBarrelBinsAnalyze()

void SiPixelPhase1Analyzer::FillBarrelBinsAnalyze	(	const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt,
		unsigned	rawId,
		const GlobalPoint &	globalPoint
	)

private

Definition at line 595 of file SiPixelPhase1Analyzer.cc.

References baseHistogramName, EcalPhiSymFlatTableProducers_cfi::id, nano_mu_digi_cff::layer, nano_mu_digi_cff::rawId, th2PolyBarrel, th2PolyBarrelSummary, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by FillBins().

                                                                                   {
   for (unsigned nameNum = 0; nameNum < baseHistogramName.size(); ++nameNum) {
     string strName = baseHistogramName[nameNum];
 
     PXBDetId id(rawId);
 
     int layer = tt->pxbLayer(id);
 
     th2PolyBarrel[strName][layer - 1]->Fill(TString::Format("%u", rawId), 1);
     th2PolyBarrelSummary[strName]->Fill(TString::Format("%u", rawId), 1);
 #ifdef DEBUG_MODE
     th2PolyBarrelDebug[strName][layer - 1]->Fill((globalPoint.y() < 0 ? globalPoint.z() + 0.5f : globalPoint.z()),
                                                  -globalPoint.x(),
                                                  (globalPoint.y() < 0 ? -1 : 1));
 #endif
   }
 }

◆ FillBarrelBinsRemap()

void SiPixelPhase1Analyzer::FillBarrelBinsRemap	(	const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt
	)

private

Definition at line 639 of file SiPixelPhase1Analyzer.cc.

References analazedRootFileName, baseHistogramName, debugFile, TrackerGeometry::detsPXB(), fs, GeomDet::geographicalId(), TFileService::getBareDirectory(), h, mps_fire::i, isBarrelSource, PixelBarrelName::ladderName(), nano_mu_digi_cff::layer, pathToHistograms, nano_mu_digi_cff::rawId, rootFileHandle, PixelBarrelName::shell(), th2PolyBarrel, th2PolyBarrelSummary, to_string(), and heppy_batch::val.

Referenced by FillBins().

                                                                                                                     {
   rootFileHandle = new TFile(analazedRootFileName[0].c_str());
 
   if (!rootFileHandle) {
     LogInfo("Analyzer") << "Could not open file: " << analazedRootFileName[0] << "..." << endl;
     return;
   }
 
 #ifdef DEBUG_MODE
   rootFileHandle->ls();
   LogInfo("Analyzer") << "\n\n";
   rootFileHandle->pwd();
   LogInfo("Analyzer") << "\n\n";
 #endif
 
   for (unsigned nameNum = 0; nameNum < baseHistogramName.size(); ++nameNum) {
     if (!isBarrelSource[nameNum])
       continue;
 
     // if (pathToHistograms[nameNum][pathToHistograms.size() - 1] != '/') pathToHistograms[nameNum] += "/";
     string baseHistogramNameWithPath = pathToHistograms[nameNum] + baseHistogramName[nameNum] + "_";
 
     const TProfile2D* handles[4];
 #ifndef DEBUG_MODE
     const TProfile2D* h;
 #endif
     bool problemWithHandles = false;
 
     for (unsigned i = 0; i < 4; ++i) {
       string fullFileName = (baseHistogramNameWithPath + std::to_string(i + 1) + ";1 ");
       handles[i] = (TProfile2D*)rootFileHandle->Get(fullFileName.c_str());
       if (!handles[i]) {
         problemWithHandles = true;
         LogInfo("Analyzer") << "Histogram: " << fullFileName << " does not exist!\n";
 
         break;
       }
     }
 
     if (!problemWithHandles) {
       LogInfo("Analyzer") << "\nInput histograms: " << baseHistogramNameWithPath << " opened successfully\n";
 
       //Add original histograms to this file
 
       TDirectory* currDir = fs->getBareDirectory()->GetDirectory(baseHistogramName[nameNum].c_str());
       currDir->cd();
 
       for (unsigned i = 0; i < 4; ++i) {
         currDir->Add(handles[i]->Clone());
       }
 
       TrackingGeometry::DetContainer pxb = theTrackerGeometry.detsPXB();
 #ifdef DEBUG_MODE
       debugFile << "There are " << pxb.size() << " detector elements in the PXB." << endl;
 #endif
       for (auto& i : pxb) {
         const GeomDet* det = i;
 
         PXBDetId id = det->geographicalId();
         unsigned rawId = id.rawId();
 
         int module = tt->pxbModule(id);
         //int ladder = tt->pxbLadder(id);
         int layer = tt->pxbLayer(id);
 
         int signedOnlineModule = module - 4;
         if (signedOnlineModule <= 0)
           --signedOnlineModule;
 
         PixelBarrelName pixelBarrelName = PixelBarrelName(id, tt, true);
         int onlineShell = pixelBarrelName.shell();
 
         int signedOnlineLadder = ((onlineShell & 1) ? -pixelBarrelName.ladderName() : pixelBarrelName.ladderName());
         string strName = baseHistogramName[nameNum];
 
 #ifdef DEBUG_MODE
         th2PolyBarrel[strName][layer - 1]->Fill(TString::Format("%u", rawId), signedOnlineLadder);
         th2PolyBarrelSummary[strName]->Fill(TString::Format("%u", rawId), signedOnlineLadder);
 #else
         h = handles[layer - 1];
         unsigned nx = h->GetNbinsX();
         unsigned ny = h->GetNbinsY();
         unsigned binX = signedOnlineModule + ((nx + 1) >> 1);
         unsigned binY = (signedOnlineLadder) + ((ny + 1) >> 1);
         double val = h->GetBinContent(binX, binY);
         th2PolyBarrel[strName][layer - 1]->Fill(TString::Format("%u", rawId), val);
         th2PolyBarrelSummary[strName]->Fill(TString::Format("%u", rawId), val);
 #endif
       }
     }
   }
 
   rootFileHandle->Close();
   delete rootFileHandle;
 }

◆ FillBins()

void SiPixelPhase1Analyzer::FillBins	(	edm::Handle< reco::TrackCollection > *	tracks,
		const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt
	)

private

Definition at line 554 of file SiPixelPhase1Analyzer.cc.

References FillBarrelBinsAnalyze(), FillBarrelBinsRemap(), FillForwardBinsAnalyze(), FillForwardBinsRemap(), mps_fire::i, TrackerGeometry::idToDet(), MODE_ANALYZE, MODE_REMAP, opMode, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, nano_mu_digi_cff::rawId, rpcPointValidation_cfi::recHit, GeomDet::surface(), Surface::toGlobal(), HLT_2024v14_cff::track, and DiMuonV_cfg::tracks.

Referenced by analyze().

                                                                 {
   switch (opMode) {
     case MODE_ANALYZE:
       for (auto const& track : *(*tracks)) {
         auto recHitsBegin = track.recHitsBegin();
         for (unsigned i = 0; i < track.recHitsSize(); ++i) {
           auto recHit = *(recHitsBegin + i);
           if (!recHit->isValid())
             continue;
 
           DetId id = recHit->geographicalId();
           unsigned subdetId = id.subdetId();
 
           if (subdetId != PixelSubdetector::PixelBarrel && subdetId != PixelSubdetector::PixelEndcap)
             continue;
 
           const PixelGeomDetUnit* geomdetunit =
               dynamic_cast<const PixelGeomDetUnit*>(theTrackerGeometry.idToDet(id));  // theTrackerGeometry ?????
           //const PixelTopology& topol = geomdetunit->specificTopology();
 
           LocalPoint localPoint = recHit->localPosition();
           GlobalPoint globalPoint = geomdetunit->surface().toGlobal(localPoint);
 
           if (subdetId == PixelSubdetector::PixelBarrel)
             FillBarrelBinsAnalyze(theTrackerGeometry, tt, id.rawId(), globalPoint);
           else
             FillForwardBinsAnalyze(theTrackerGeometry, tt, id.rawId(), globalPoint);
         }
       }
       break;
     case MODE_REMAP:
       FillBarrelBinsRemap(theTrackerGeometry, tt);
       FillForwardBinsRemap(theTrackerGeometry, tt);
       break;
     default:
       break;
   }
 }

◆ fillDescriptions()

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

static

Definition at line 823 of file SiPixelPhase1Analyzer.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                        {
   //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.setComment(
       "Creates TH2Poly Pixel Tracker maps by either analyzing the event or remapping exising DQM historams");
   desc.add<edm::InputTag>("src", edm::InputTag("generalTracks"));
   desc.addUntracked<unsigned int>("opMode", 1);
   desc.addUntracked<std::string>("debugFileName", "debug.txt");
   desc.addUntracked<std::vector<unsigned int>>("isBarrelSource", {0, 0, 1});
   desc.addUntracked<std::vector<std::string>>("remapRootFileName", {"dqmFile.root"});
   desc.addUntracked<std::vector<std::string>>(
       "pathToHistograms",
       {"DQMData/Run 1/PixelPhase1/Run summary/Phase1_MechanicalView/PXForward/",
        "DQMData/Run 1/PixelPhase1/Run summary/Phase1_MechanicalView/PXForward/",
        "DQMData/Run 1/PixelPhase1/Run summary/Phase1_MechanicalView/PXBarrel/"});
   desc.addUntracked<std::vector<std::string>>("baseHistogramName",
                                               {"num_clusters_per_PXDisk_per_SignedBladePanel_PXRing",
                                                "num_digis_per_PXDisk_per_SignedBladePanel_PXRing",
                                                "num_digis_per_SignedModule_per_SignedLadder_PXLayer"});
   descriptions.addWithDefaultLabel(desc);
 }

◆ FillForwardBinsAnalyze()

void SiPixelPhase1Analyzer::FillForwardBinsAnalyze	(	const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt,
		unsigned	rawId,
		const GlobalPoint &	globalPoint
	)

private

Definition at line 616 of file SiPixelPhase1Analyzer.cc.

References baseHistogramName, EcalPhiSymFlatTableProducers_cfi::id, pxfTh2PolyForward, pxfTh2PolyForwardSummary, nano_mu_digi_cff::rawId, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by FillBins().

                                                                                    {
   for (unsigned nameNum = 0; nameNum < baseHistogramName.size(); ++nameNum) {
     string strName = baseHistogramName[nameNum];
 
     PXFDetId id(rawId);
 
     int side = tt->side(id);   //tt->pxfSide(id);
     int disk = tt->layer(id);  //tt->pxfDisk(id);
     // int blade = tt->pxfBlade(id);
     unsigned mapIdx = disk + (side - 1) * 3 - 1;
 
     pxfTh2PolyForward[strName][mapIdx]->Fill(TString::Format("%u", rawId), 1);
     pxfTh2PolyForwardSummary[strName]->Fill(TString::Format("%u", rawId), 1);
 
 #ifdef DEBUG_MODE
     pxfTh2PolyForwardDebug[strName][mapIdx]->Fill(globalPoint.x(), globalPoint.y(), 1);
 #endif
   }
 }

◆ FillForwardBinsRemap()

void SiPixelPhase1Analyzer::FillForwardBinsRemap	(	const TrackerGeometry &	theTrackerGeometry,
		const TrackerTopology *	tt
	)

private

Definition at line 735 of file SiPixelPhase1Analyzer.cc.

References analazedRootFileName, baseHistogramName, PixelEndcapName::bladeName(), TrackerGeometry::detsPXF(), fs, GeomDet::geographicalId(), TFileService::getBareDirectory(), h, PixelEndcapName::halfCylinder(), mps_fire::i, isBarrelSource, PixelEndcapName::pannelName(), pathToHistograms, pxfTh2PolyForward, pxfTh2PolyForwardSummary, nano_mu_digi_cff::rawId, PixelEndcapName::ringName(), rootFileHandle, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), and heppy_batch::val.

Referenced by FillBins().

                                                                                                                      {
   rootFileHandle = new TFile(analazedRootFileName[0].c_str());
 
   if (!rootFileHandle) {
     return;
   }
 
   TrackingGeometry::DetContainer pxf = theTrackerGeometry.detsPXF();
 
 #ifdef DEBUG_MODE
   rootFileHandle->ls();
   LogInfo("Analyzer") << "\n\n";
   rootFileHandle->pwd();
   LogInfo("Analyzer") << "\n\n";
 #endif
 
   for (unsigned nameNum = 0; nameNum < baseHistogramName.size(); ++nameNum) {
     if (isBarrelSource[nameNum])
       continue;
 
     string baseHistogramNameWithPath = pathToHistograms[nameNum] + baseHistogramName[nameNum] + "_";
 
     const TProfile2D* h_1 = (TProfile2D*)rootFileHandle->Get((baseHistogramNameWithPath + "1;1 ").c_str());
     const TProfile2D* h_2 = (TProfile2D*)rootFileHandle->Get((baseHistogramNameWithPath + "2;1 ").c_str());
 #ifndef DEBUG_MODE
     const TProfile2D* h;
 #endif
     if (h_2 && h_1) {
       LogInfo("Analyzer") << "\nInput histograms: " << baseHistogramNameWithPath << " opened successfully\n";
 
       //Add original histograms to this file
       TDirectory* currDir = fs->getBareDirectory()->GetDirectory(baseHistogramName[nameNum].c_str());
       currDir->cd();
       currDir->Add(h_1->Clone());
       currDir->Add(h_2->Clone());
 
       for (auto& i : pxf) {
         const GeomDet* det = i;
 
         PXFDetId id = det->geographicalId();
 
         int side = tt->side(id);   //tt->pxfSide(id);
         int disk = tt->layer(id);  //tt->pxfDisk(id);
 
         unsigned rawId = id.rawId();
         PixelEndcapName pixelEndcapName = PixelEndcapName(PXFDetId(rawId), tt, true);
 
 #ifdef DEBUG_MODE
         int blade = tt->pxfBlade(id);
 #else
         int onlineBlade = pixelEndcapName.bladeName();
         bool isInnerOnlineBlade = !(pixelEndcapName.halfCylinder() & 1);  // inner -> blade > 0 (?)
 
         int signedOnlineBlade = (isInnerOnlineBlade) ? onlineBlade : -onlineBlade;
         int signedDisk = (side == 2) ? disk : -disk;
 
         int pannel = pixelEndcapName.pannelName() - 1;
 
 #endif
         unsigned mapIdx = disk + (side - 1) * 3 - 1;
         string strName = baseHistogramName[nameNum];
 
 #ifdef DEBUG_MODE
         pxfTh2PolyForward[strName][mapIdx]->Fill(TString::Format("%u", rawId), blade);
         pxfTh2PolyForwardSummary[strName]->Fill(TString::Format("%u", rawId), blade);
 #else
         if (pixelEndcapName.ringName() == 1)
           h = h_1;
         else
           h = h_2;
         // ---- REMAP (Online -> Offline)
         unsigned nx = h->GetNbinsX();
         unsigned ny = h->GetNbinsY();
         unsigned binX = signedDisk + ((nx + 1) >> 1);
         unsigned binY = (signedOnlineBlade * 2) + (ny >> 1);
         double val = h->GetBinContent(binX, binY + pannel);
         pxfTh2PolyForward[strName][mapIdx]->Fill(TString::Format("%u", rawId), val);
         pxfTh2PolyForwardSummary[strName]->Fill(TString::Format("%u", rawId), val);
 #endif
       }
     }
   }
 
   rootFileHandle->Close();
   delete rootFileHandle;
 }

◆ SaveDetectorData()

void SiPixelPhase1Analyzer::SaveDetectorData	(	bool	isBarrel,
		unsigned	rawId,
		int	shell_hc,
		int	layer_disk,
		int	ladder_blade
	)

inlineprivate

Definition at line 157 of file SiPixelPhase1Analyzer.h.

References geometryDiff::file, PixelPluginsPhase0_cfi::isBarrel, MillePedeFileConverter_cfg::out, and nano_mu_digi_cff::rawId.

Referenced by BookBarrelBins(), and BookForwardBins().

                                                                                                        {
     std::ofstream file("det.data", std::ofstream::out);
 
     file << isBarrel << "\t" << rawId << "\t" << shell_hc << "\t" << layer_disk << "\t" << ladder_blade << std::endl;
   }

◆ SaveDetectorVertices()

void SiPixelPhase1Analyzer::SaveDetectorVertices ( const TrackerTopology * tt )

private

Definition at line 486 of file SiPixelPhase1Analyzer.cc.

References newFWLiteAna::bin, bins, f, mps_fire::i, EcalPhiSymFlatTableProducers_cfi::id, dqmiolumiharvest::j, nano_mu_digi_cff::layer, PixelBarrelName::name(), PixelEndcapName::name(), MillePedeFileConverter_cfg::out, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, nano_mu_digi_cff::rawId, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), and to_string().

Referenced by BookBins().

                                                                           {
   vector<std::ofstream*> verticesFiles[2];
   for (unsigned i = 0; i < 4; ++i) {
     std::ofstream* f = new std::ofstream(("vertices_barrel_" + std::to_string(i + 1)).c_str(), std::ofstream::out);
 
     verticesFiles[0].push_back(f);
   }
 
   for (unsigned side = 1; side <= 2; ++side) {
     for (unsigned disk = 1; disk <= 3; ++disk) {
       std::ofstream* f = new std::ofstream(
           ("vertices_forward_" + std::to_string((side == 1 ? -(int(disk)) : (int)disk))).c_str(), std::ofstream::out);
 
       verticesFiles[1].push_back(f);
     }
   }
 
   for (auto& bin : bins) {
     unsigned rawId = bin.first;
     DetId id(rawId);
     unsigned subdetId = id.subdetId();
 
     if (subdetId != PixelSubdetector::PixelBarrel && subdetId != PixelSubdetector::PixelEndcap)
       continue;
 
     double* vertX = bin.second->GetX();
     double* vertY = bin.second->GetY();
 
     if (subdetId == PixelSubdetector::PixelBarrel) {
       PXBDetId pxbId(rawId);
       unsigned layer = tt->pxbLayer(pxbId);
       string onlineName = PixelBarrelName(pxbId, tt, true).name();
 
       *(verticesFiles[0][layer - 1]) << rawId << " " << onlineName << " \"";
       for (unsigned i = 0; i < 4; ++i) {
         *(verticesFiles[0][layer - 1]) << vertX[i] << "," << vertY[i];
         if (i == 3)
           *(verticesFiles[0][layer - 1]) << "\"\n";
         else
           *(verticesFiles[0][layer - 1]) << " ";
       }
     } else {
       PXFDetId pxfId(rawId);
       unsigned side = tt->pxfSide(pxfId);
       unsigned disk = tt->pxfDisk(pxfId);
       string onlineName = PixelEndcapName(pxfId, tt, true).name();
       unsigned mapIdx = disk + (side - 1) * 3 - 1;
 
       *(verticesFiles[1][mapIdx]) << rawId << " " << onlineName << " \"";
       for (unsigned i = 0; i < 3; ++i) {
         *(verticesFiles[1][mapIdx]) << vertX[i] << "," << vertY[i];
 
         if (i == 2)
           *(verticesFiles[1][mapIdx]) << "\"\n";
         else
           *(verticesFiles[1][mapIdx]) << " ";
       }
     }
   }
 
   for (unsigned i = 0; i < 2; ++i) {
     for (auto& j : verticesFiles[i]) {
       j->close();
       delete j;
     }
   }
 }