RecoLocalTracker/SiPhase2VectorHitBuilder/test/ClustersValidationTGraph.cc - master - githubmirror/cmssw

RecoLocalTracker/SiPhase2VectorHitBuilder/test/ClustersValidationTGraph.cc
549 строк22 KB
Shahzad Malik Muzaffar
[RECONSTRUCTION] Drop Geometry/CommonDetUnit package
28 май 2026, 12:57
28 май 2026, 12:57685f036
// system includes
#include <map>
#include <vector>
#include <algorithm>

// user includes
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "CommonTools/Utils/interface/TFileDirectory.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/DetSetVectorNew.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/Phase2TrackerCluster/interface/Phase2TrackerCluster1D.h"
#include "DataFormats/Phase2TrackerDigi/interface/Phase2TrackerDigi.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "DataFormats/SiPixelDigi/interface/PixelDigi.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "Geometry/CommonTopologies/interface/GeomDet.h"
#include "Geometry/CommonTopologies/interface/PixelGeomDetUnit.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"
#include "SimDataFormats/TrackerDigiSimLink/interface/PixelDigiSimLink.h"
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "SimDataFormats/Vertex/interface/SimVertexContainer.h"

// ROOT includes
#include <TGraph.h>
#include <TH1F.h>
#include <THStack.h>
#include <TTree.h>

struct ClusterHistos {
  THStack* numberClustersMixed;
  TH1F* numberClusterPixel;
  TH1F* numberClusterStrip;

  THStack* clustersSizeMixed;
  TH1F* clusterSizePixel;
  TH1F* clusterSizeStrip;

  TGraph* globalPosXY[3];
  TGraph* localPosXY[3];

  TH1F* deltaXClusterSimHits[3];
  TH1F* deltaYClusterSimHits[3];

  TH1F* deltaXClusterSimHits_P[3];
  TH1F* deltaYClusterSimHits_P[3];

  TH1F* digiEfficiency[3];

  TH1F* primarySimHits;
  TH1F* otherSimHits;
};

class Phase2TrackerClusterizerValidationTGraph : public edm::one::EDAnalyzer<edm::one::SharedResources> {
public:
  typedef std::map<unsigned int, std::vector<PSimHit> > SimHitsMap;
  typedef std::map<unsigned int, SimTrack> SimTracksMap;

  explicit Phase2TrackerClusterizerValidationTGraph(const edm::ParameterSet&);
  ~Phase2TrackerClusterizerValidationTGraph();
  void beginJob();
  void analyze(const edm::Event&, const edm::EventSetup&);
  static void fillDescriptions(edm::ConfigurationDescriptions&);

private:
  std::map<unsigned int, ClusterHistos>::iterator createLayerHistograms(unsigned int);
  unsigned int getLayerNumber(const DetId&, const TrackerTopology*);
  unsigned int getModuleNumber(const DetId&, const TrackerTopology*);
  unsigned int getSimTrackId(const edm::Handle<edm::DetSetVector<PixelDigiSimLink> >&, const DetId&, unsigned int);

  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> esTokenGeom_;
  const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> esTokenTopo_;
  const edm::EDGetTokenT<edmNew::DetSetVector<Phase2TrackerCluster1D> > srcClu_;
  const edm::EDGetTokenT<edm::DetSetVector<PixelDigiSimLink> > siphase2OTSimLinksToken_;
  const edm::EDGetTokenT<edm::PSimHitContainer> simHitsToken_;
  const edm::EDGetTokenT<edm::SimTrackContainer> simTracksToken_;
  const edm::EDGetTokenT<edm::SimVertexContainer> simVerticesToken_;
  const TrackerGeometry* tkGeom_;
  const TrackerTopology* tkTopo_;

  TTree* tree_;
  TGraph* trackerLayout_[3];
  TGraph* trackerLayoutXY_[3];
  TGraph* trackerLayoutXYBar_;
  TGraph* trackerLayoutXYEC_;

  std::map<unsigned int, ClusterHistos> histograms_;
};

Phase2TrackerClusterizerValidationTGraph::Phase2TrackerClusterizerValidationTGraph(const edm::ParameterSet& conf)
    : esTokenGeom_(esConsumes()),
      esTokenTopo_(esConsumes()),
      srcClu_(consumes<edmNew::DetSetVector<Phase2TrackerCluster1D> >(
          edm::InputTag(conf.getParameter<std::string>("src")))),
      siphase2OTSimLinksToken_(
          consumes<edm::DetSetVector<PixelDigiSimLink> >(conf.getParameter<edm::InputTag>("links"))),
      simHitsToken_(consumes<edm::PSimHitContainer>(edm::InputTag("g4SimHits", "TrackerHitsPixelBarrelLowTof"))),
      simTracksToken_(consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"))),
      simVerticesToken_(consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"))) {
  usesResource(TFileService::kSharedResource);
}

Phase2TrackerClusterizerValidationTGraph::~Phase2TrackerClusterizerValidationTGraph() = default;

void Phase2TrackerClusterizerValidationTGraph::beginJob() {
  edm::Service<TFileService> fs;
  fs->file().cd("/");
  TFileDirectory td = fs->mkdir("Common");
  //Create common ntuple
  tree_ = td.make<TTree>("Phase2TrackerClusters", "Phase2TrackerClusters");
  // Create common histograms
  trackerLayout_[0] = td.make<TGraph>();
  trackerLayout_[0]->SetName("RVsZ_Mixed");
  trackerLayout_[1] = td.make<TGraph>();
  trackerLayout_[1]->SetName("RVsZ_Pixel");
  trackerLayout_[2] = td.make<TGraph>();
  trackerLayout_[2]->SetName("RVsZ_Strip");
  trackerLayoutXY_[0] = td.make<TGraph>();
  trackerLayoutXY_[0]->SetName("YVsX_Mixed");
  trackerLayoutXY_[1] = td.make<TGraph>();
  trackerLayoutXY_[1]->SetName("YVsX_Pixel");
  trackerLayoutXY_[2] = td.make<TGraph>();
  trackerLayoutXY_[2]->SetName("YVsX_Strip");
  trackerLayoutXYBar_ = td.make<TGraph>();
  trackerLayoutXYBar_->SetName("YVsXBar");
  trackerLayoutXYEC_ = td.make<TGraph>();
  trackerLayoutXYEC_->SetName("YVsXEC");
}

void Phase2TrackerClusterizerValidationTGraph::analyze(const edm::Event& event, const edm::EventSetup& eventSetup) {
  // Get the needed objects

  // Get the clusters
  edm::Handle<Phase2TrackerCluster1DCollectionNew> clusters;
  event.getByToken(srcClu_, clusters);

  // Get the Phase2 DigiSimLink
  edm::Handle<edm::DetSetVector<PixelDigiSimLink> > siphase2SimLinks;
  event.getByToken(siphase2OTSimLinksToken_, siphase2SimLinks);

  // Get the SimHits
  edm::Handle<edm::PSimHitContainer> simHitsRaw;
  event.getByToken(simHitsToken_, simHitsRaw);

  // Get the SimTracks
  edm::Handle<edm::SimTrackContainer> simTracksRaw;
  event.getByToken(simTracksToken_, simTracksRaw);

  // Get the SimVertex
  edm::Handle<edm::SimVertexContainer> simVertices;
  event.getByToken(simVerticesToken_, simVertices);

  // Get the geometry
  tkGeom_ = &eventSetup.getData(esTokenGeom_);
  tkTopo_ = &eventSetup.getData(esTokenTopo_);

  //set up for tree
  int layer_number;
  int module_id;
  int module_number;
  int module_type;  //1: pixel, 2: strip
  float x_global, y_global, z_global;
  float x_local, y_local, z_local;

  tree_->Branch("layer_number", &layer_number, "layer_number/I");
  tree_->Branch("module_id", &module_id, "module_id/I");
  tree_->Branch("module_type", &module_type, "module_type/I");
  tree_->Branch("module_number", &module_number, "module_number/I");
  tree_->Branch("x_global", &x_global, "x_global/F");
  tree_->Branch("y_global", &y_global, "y_global/F");
  tree_->Branch("z_global", &z_global, "z_global/F");
  tree_->Branch("x_local", &x_local, "x_local/F");
  tree_->Branch("y_local", &y_local, "y_local/F");
  tree_->Branch("z_local", &z_local, "z_local/F");

  // Rearrange the simTracks for ease of use <simTrackID, simTrack>
  SimTracksMap simTracks;
  for (const auto& simTrackIt : *simTracksRaw)
    simTracks.emplace(std::pair<unsigned int, SimTrack>(simTrackIt.trackId(), simTrackIt));

  // Rearrange the simHits by detUnit

  // Rearrange the simHits for ease of use
  SimHitsMap simHitsDetUnit;
  SimHitsMap simHitsTrackId;
  for (const auto& simHitIt : *simHitsRaw) {
    SimHitsMap::iterator simHitsDetUnitIt(simHitsDetUnit.find(simHitIt.detUnitId()));
    if (simHitsDetUnitIt == simHitsDetUnit.end()) {
      std::pair<SimHitsMap::iterator, bool> newIt(simHitsDetUnit.insert(
          std::pair<unsigned int, std::vector<PSimHit> >(simHitIt.detUnitId(), std::vector<PSimHit>())));
      simHitsDetUnitIt = newIt.first;
    }
    simHitsDetUnitIt->second.push_back(simHitIt);
    SimHitsMap::iterator simHitsTrackIdIt(simHitsTrackId.find(simHitIt.trackId()));
    if (simHitsTrackIdIt == simHitsTrackId.end()) {
      std::pair<SimHitsMap::iterator, bool> newIt(simHitsTrackId.insert(
          std::pair<unsigned int, std::vector<PSimHit> >(simHitIt.trackId(), std::vector<PSimHit>())));
      simHitsTrackIdIt = newIt.first;
    }
    simHitsTrackIdIt->second.push_back(simHitIt);
  }

  // ValidationTGraph
  unsigned int nClustersTot(0), nClustersPixelTot(0), nClustersStripTot(0);

  // Loop over modules
  for (const auto& DSViter : *clusters) {
    // Get the detector unit's id
    unsigned int rawid(DSViter.detId());
    module_id = rawid;
    DetId detId(rawid);

    layer_number = getLayerNumber(detId, tkTopo_);
    module_number = getModuleNumber(detId, tkTopo_);
    unsigned int layer(getLayerNumber(detId, tkTopo_));

    // Get the geometry of the tracker
    const GeomDetUnit* geomDetUnit(tkGeom_->idToDetUnit(detId));
    const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(geomDetUnit);
    const PixelTopology& topol = theGeomDet->specificTopology();

    if (!geomDetUnit)
      break;

    // Create histograms for the layer if they do not yet exist
    std::map<unsigned int, ClusterHistos>::iterator histogramLayer(histograms_.find(layer));
    if (histogramLayer == histograms_.end())
      histogramLayer = createLayerHistograms(layer);

    // Number of clusters
    unsigned int nClustersPixel(0), nClustersStrip(0);

    // Loop over the clusters in the detector unit
    for (const auto& clustIt : DSViter) {
      // Cluster related variables
      MeasurementPoint mpClu(clustIt.center(), clustIt.column() + 0.5);
      Local3DPoint localPosClu = geomDetUnit->topology().localPosition(mpClu);
      x_local = localPosClu.x();
      y_local = localPosClu.y();
      z_local = localPosClu.z();

      Global3DPoint globalPosClu = geomDetUnit->surface().toGlobal(localPosClu);
      x_global = globalPosClu.x();
      y_global = globalPosClu.y();
      z_global = globalPosClu.z();

      // Fill the position histograms
      trackerLayout_[0]->SetPoint(nClustersTot, globalPosClu.z(), globalPosClu.perp());
      trackerLayoutXY_[0]->SetPoint(nClustersTot, globalPosClu.x(), globalPosClu.y());

      if (layer < 100)
        trackerLayoutXYBar_->SetPoint(nClustersTot, globalPosClu.x(), globalPosClu.y());
      else
        trackerLayoutXYEC_->SetPoint(nClustersTot, globalPosClu.x(), globalPosClu.y());

      histogramLayer->second.localPosXY[0]->SetPoint(nClustersTot, localPosClu.x(), localPosClu.y());
      histogramLayer->second.globalPosXY[0]->SetPoint(nClustersTot, globalPosClu.x(), globalPosClu.y());

      // Pixel module
      if (topol.ncolumns() == 32) {
        module_type = 1;
        trackerLayout_[1]->SetPoint(nClustersPixelTot, globalPosClu.z(), globalPosClu.perp());
        trackerLayoutXY_[1]->SetPoint(nClustersPixelTot, globalPosClu.x(), globalPosClu.y());

        histogramLayer->second.localPosXY[1]->SetPoint(nClustersPixelTot, localPosClu.x(), localPosClu.y());
        histogramLayer->second.globalPosXY[1]->SetPoint(nClustersPixelTot, globalPosClu.x(), globalPosClu.y());
        histogramLayer->second.clusterSizePixel->Fill(clustIt.size());
        ++nClustersPixel;
        ++nClustersPixelTot;
      }
      // Strip module
      else if (topol.ncolumns() == 2) {
        module_type = 2;
        trackerLayout_[2]->SetPoint(nClustersStripTot, globalPosClu.z(), globalPosClu.perp());
        trackerLayoutXY_[2]->SetPoint(nClustersStripTot, globalPosClu.x(), globalPosClu.y());

        histogramLayer->second.localPosXY[2]->SetPoint(nClustersStripTot, localPosClu.x(), localPosClu.y());
        histogramLayer->second.globalPosXY[2]->SetPoint(nClustersStripTot, globalPosClu.x(), globalPosClu.y());
        histogramLayer->second.clusterSizeStrip->Fill(clustIt.size());
        ++nClustersStrip;
        ++nClustersStripTot;
      }

      // * Digis related variables

      std::vector<unsigned int> clusterSimTrackIds;

      // Get all the simTracks that form the cluster
      for (unsigned int i(0); i < clustIt.size(); ++i) {
        unsigned int channel(PixelDigi::pixelToChannel(
            clustIt.firstRow() + i,
            clustIt
                .column()));  // Here we have to use the old pixelToChannel function (not Phase2TrackerDigi but PixelDigi), change this when using new Digis
        unsigned int simTrackId(getSimTrackId(siphase2SimLinks, detId, channel));
        clusterSimTrackIds.push_back(simTrackId);
      }
    }

    if (nClustersPixel)
      histogramLayer->second.numberClusterPixel->Fill(nClustersPixel);
    if (nClustersStrip)
      histogramLayer->second.numberClusterStrip->Fill(nClustersStrip);
    nClustersTot++;
    tree_->Fill();
  }
}

// Create the histograms
std::map<unsigned int, ClusterHistos>::iterator Phase2TrackerClusterizerValidationTGraph::createLayerHistograms(
    unsigned int ival) {
  std::ostringstream fname1, fname2;

  edm::Service<TFileService> fs;
  fs->file().cd("/");

  std::string tag;
  unsigned int id;
  if (ival < 100) {
    id = ival;
    fname1 << "Barrel";
    fname2 << "Layer_" << id;
    tag = "_layer_";
  } else {
    int side = ival / 100;
    id = ival - side * 100;
    fname1 << "EndCap_Side_" << side;
    fname2 << "Disc_" << id;
    tag = "_disc_";
  }

  TFileDirectory td1 = fs->mkdir(fname1.str().c_str());
  TFileDirectory td = td1.mkdir(fname2.str().c_str());

  ClusterHistos local_histos;

  std::ostringstream histoName;

  /*
     * Number of clusters
     */

  histoName.str("");
  histoName << "Number_Clusters_Pixel" << tag.c_str() << id;
  local_histos.numberClusterPixel = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 20, 0., 20.);
  local_histos.numberClusterPixel->SetFillColor(kAzure + 7);

  histoName.str("");
  histoName << "Number_Clusters_Strip" << tag.c_str() << id;
  local_histos.numberClusterStrip = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 20, 0., 20.);
  local_histos.numberClusterStrip->SetFillColor(kOrange - 3);

  histoName.str("");
  histoName << "Number_Clusters_Mixed" << tag.c_str() << id;
  local_histos.numberClustersMixed = td.make<THStack>(histoName.str().c_str(), histoName.str().c_str());
  local_histos.numberClustersMixed->Add(local_histos.numberClusterPixel);
  local_histos.numberClustersMixed->Add(local_histos.numberClusterStrip);

  /*
     * Cluster size
     */

  histoName.str("");
  histoName << "Cluster_Size_Pixel" << tag.c_str() << id;
  local_histos.clusterSizePixel = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 10, 0., 10.);
  local_histos.clusterSizePixel->SetFillColor(kAzure + 7);

  histoName.str("");
  histoName << "Cluster_Size_Strip" << tag.c_str() << id;
  local_histos.clusterSizeStrip = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 10, 0., 10.);
  local_histos.clusterSizeStrip->SetFillColor(kOrange - 3);

  histoName.str("");
  histoName << "Cluster_Size_Mixed" << tag.c_str() << id;
  local_histos.clustersSizeMixed = td.make<THStack>(histoName.str().c_str(), histoName.str().c_str());
  local_histos.clustersSizeMixed->Add(local_histos.clusterSizePixel);
  local_histos.clustersSizeMixed->Add(local_histos.clusterSizeStrip);

  /*
     * Local and Global positions
     */

  histoName.str("");
  histoName << "Local_Position_XY_Mixed" << tag.c_str() << id;
  local_histos.localPosXY[0] = td.make<TGraph>();
  local_histos.localPosXY[0]->SetName(histoName.str().c_str());

  histoName.str("");
  histoName << "Local_Position_XY_Pixel" << tag.c_str() << id;
  local_histos.localPosXY[1] = td.make<TGraph>();
  local_histos.localPosXY[1]->SetName(histoName.str().c_str());

  histoName.str("");
  histoName << "Local_Position_XY_Strip" << tag.c_str() << id;
  local_histos.localPosXY[2] = td.make<TGraph>();
  local_histos.localPosXY[2]->SetName(histoName.str().c_str());

  histoName.str("");
  histoName << "Global_Position_XY_Mixed" << tag.c_str() << id;
  local_histos.globalPosXY[0] = td.make<TGraph>();
  local_histos.globalPosXY[0]->SetName(histoName.str().c_str());

  histoName.str("");
  histoName << "Global_Position_XY_Pixel" << tag.c_str() << id;
  local_histos.globalPosXY[1] = td.make<TGraph>();
  local_histos.globalPosXY[1]->SetName(histoName.str().c_str());

  histoName.str("");
  histoName << "Global_Position_XY_Strip" << tag.c_str() << id;
  local_histos.globalPosXY[2] = td.make<TGraph>();
  local_histos.globalPosXY[2]->SetName(histoName.str().c_str());

  /*
     * Delta positions with SimHits
     */

  histoName.str("");
  histoName << "Delta_X_Cluster_SimHits_Mixed" << tag.c_str() << id;
  local_histos.deltaXClusterSimHits[0] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_X_Cluster_SimHits_Pixel" << tag.c_str() << id;
  local_histos.deltaXClusterSimHits[1] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_X_Cluster_SimHits_Strip" << tag.c_str() << id;
  local_histos.deltaXClusterSimHits[2] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_Y_Cluster_SimHits_Mixed" << tag.c_str() << id;
  local_histos.deltaYClusterSimHits[0] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_Y_Cluster_SimHits_Pixel" << tag.c_str() << id;
  local_histos.deltaYClusterSimHits[1] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_Y_Cluster_SimHits_Strip" << tag.c_str() << id;
  local_histos.deltaYClusterSimHits[2] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  /*
     * Delta position with simHits for primary tracks only
     */

  histoName.str("");
  histoName << "Delta_X_Cluster_SimHits_Mixed_P" << tag.c_str() << id;
  local_histos.deltaXClusterSimHits_P[0] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_X_Cluster_SimHits_Pixel_P" << tag.c_str() << id;
  local_histos.deltaXClusterSimHits_P[1] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_X_Cluster_SimHits_Strip_P" << tag.c_str() << id;
  local_histos.deltaXClusterSimHits_P[2] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_Y_Cluster_SimHits_Mixed_P" << tag.c_str() << id;
  local_histos.deltaYClusterSimHits_P[0] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_Y_Cluster_SimHits_Pixel_P" << tag.c_str() << id;
  local_histos.deltaYClusterSimHits_P[1] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  histoName.str("");
  histoName << "Delta_Y_Cluster_SimHits_Strip_P" << tag.c_str() << id;
  local_histos.deltaYClusterSimHits_P[2] = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 100, 0., 0.);

  /*
     * Information on the Digis per cluster
     */

  histoName.str("");
  histoName << "Primary_Digis" << tag.c_str() << id;
  local_histos.primarySimHits = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 10, 0., 10.);

  histoName.str("");
  histoName << "Other_Digis" << tag.c_str() << id;
  local_histos.otherSimHits = td.make<TH1F>(histoName.str().c_str(), histoName.str().c_str(), 10, 0., 10.);

  /*
     * End
     */

  std::pair<std::map<unsigned int, ClusterHistos>::iterator, bool> insertedIt(
      histograms_.insert(std::make_pair(ival, local_histos)));
  fs->file().cd("/");

  return insertedIt.first;
}

unsigned int Phase2TrackerClusterizerValidationTGraph::getLayerNumber(const DetId& detid, const TrackerTopology* topo) {
  if (detid.det() == DetId::Tracker) {
    if (detid.subdetId() == PixelSubdetector::PixelBarrel)
      return (topo->pxbLayer(detid));
    else if (detid.subdetId() == PixelSubdetector::PixelEndcap)
      return (100 * topo->pxfSide(detid) + topo->pxfDisk(detid));
  }
  return 999;
}

unsigned int Phase2TrackerClusterizerValidationTGraph::getModuleNumber(const DetId& detid,
                                                                       const TrackerTopology* topo) {
  if (detid.det() == DetId::Tracker) {
    if (detid.subdetId() == PixelSubdetector::PixelBarrel) {
      return (topo->pxbModule(detid));
    } else if (detid.subdetId() == PixelSubdetector::PixelEndcap) {
      return (topo->pxfModule(detid));
    } else
      return 999;
  }
  return 999;
}

unsigned int Phase2TrackerClusterizerValidationTGraph::getSimTrackId(
    const edm::Handle<edm::DetSetVector<PixelDigiSimLink> >& siphase2SimLinks,
    const DetId& detId,
    unsigned int channel) {
  edm::DetSetVector<PixelDigiSimLink>::const_iterator DSViter(siphase2SimLinks->find(detId));
  if (DSViter == siphase2SimLinks->end())
    return 0;
  for (const auto& it : DSViter->data) {
    if (channel == it.channel())
      return it.SimTrackId();
  }
  return 0;
}

void Phase2TrackerClusterizerValidationTGraph::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.add<std::string>("src", "siPhase2Clusters");
  desc.add<edm::InputTag>("links", edm::InputTag("simSiPixelDigis", "Tracker"));
  descriptions.add("phase2TrackerClusterizerValidationTGraph", desc);
}

DEFINE_FWK_MODULE(Phase2TrackerClusterizerValidationTGraph);