Inheritance diagram for AlignmentStats:

Public Member Functions
	AlignmentStats (const edm::ParameterSet &iConfig)

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

void	beginJob () override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	endJob () override

void	endRun (edm::Run const &iRun, edm::EventSetup const &iSetup) override

	~AlignmentStats () override=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	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 Types
typedef std::map< uint32_t, uint32_t >	DetHitMap

Private Attributes
float	Chi2n [MAXTRKS_]

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	esTokenTkGeoER_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	esTokenTTopoER_

float	Eta [MAXTRKS_]

int	event_

std::vector< aliStats::GeoInfo >	geomInfoList_

DetHitMap	hitmap_

TFile *	hitsfile_

TTree *	hitstree_

const std::string	hitsTreeName_

unsigned int	id_

bool	is2D_

bool	isStereo_

const bool	keepHitPopulation_

const bool	keepTrackStats_

unsigned int	layer_

const edm::EDGetTokenT< AliClusterValueMap >	mapToken_

int	Nhits [MAXTRKS_][7]

unsigned int	nhits_

unsigned int	noverlaps_

unsigned int	ntracks

TTree *	outtree_

const edm::InputTag	overlapAM_

DetHitMap	overlapmap_

float	P [MAXTRKS_]

float	Phi [MAXTRKS_]

float	posEta_

float	posPhi_

float	posR_

float	posX_

float	posY_

float	posZ_

const uint32_t	prescale_

float	Pt [MAXTRKS_]

int	run_

const edm::InputTag	src_

const std::string	statsTreeName_

int	subdet_

uint32_t	tmpPresc_

std::unique_ptr< TrackerTopology >	trackerTopology_

const edm::EDGetTokenT< reco::TrackCollection >	trackToken_

TFile *	treefile_

Static Private Attributes
static const int	MAXTRKS_ = 200

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

Definition at line 69 of file AlignmentStats.cc.

Member Typedef Documentation

◆ DetHitMap

typedef std::map<uint32_t, uint32_t> AlignmentStats::DetHitMap

private

Definition at line 122 of file AlignmentStats.cc.

Constructor & Destructor Documentation

◆ AlignmentStats()

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

Definition at line 131 of file AlignmentStats.cc.

References outtree_.

     : esTokenTTopoER_(esConsumes<edm::Transition::EndRun>()),
       esTokenTkGeoER_(esConsumes<edm::Transition::EndRun>()),
       src_(iConfig.getParameter<edm::InputTag>("src")),
       overlapAM_(iConfig.getParameter<edm::InputTag>("OverlapAssoMap")),
       keepTrackStats_(iConfig.getParameter<bool>("keepTrackStats")),
       keepHitPopulation_(iConfig.getParameter<bool>("keepHitStats")),
       statsTreeName_(iConfig.getParameter<string>("TrkStatsFileName")),
       hitsTreeName_(iConfig.getParameter<string>("HitStatsFileName")),
       prescale_(iConfig.getParameter<uint32_t>("TrkStatsPrescale")),
       trackToken_(consumes<reco::TrackCollection>(src_)),
       mapToken_(consumes<AliClusterValueMap>(overlapAM_)) {
   //sanity checks
 
   //init
   outtree_ = nullptr;
 
 }  //end constructor

◆ ~AlignmentStats()

AlignmentStats::~AlignmentStats ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 278 of file AlignmentStats.cc.

References Chi2n, className(), SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), Eta, event_, Exception, hitmap_, iEvent, AlignmentClusterFlag::isOverlap(), dqmiolumiharvest::j, dqmdumpme::k, mapToken_, MAXTRKS_, Nhits, ntracks, outtree_, overlapmap_, Phi, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixelClusterTagInfos_cfi::pixelhit, prescale_, Pt, nano_mu_digi_cff::rawId, run_, DetId::subdetId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, tmpPresc_, SiStripDetId::TOB, gather_cfg::Tracks, and trackToken_.

                                                                                 {
   //take trajectories and tracks to loop on
   // edm::Handle<TrajTrackAssociationCollection> TrackAssoMap;
   const edm::Handle<reco::TrackCollection> &Tracks = iEvent.getHandle(trackToken_);
 
   //take overlap HitAssomap
   const edm::Handle<AliClusterValueMap> &hMap = iEvent.getHandle(mapToken_);
   const AliClusterValueMap &OverlapMap = *hMap;
 
   // Initialise
   run_ = 1;
   event_ = 1;
   ntracks = 0;
   run_ = iEvent.id().run();
   event_ = iEvent.id().event();
   ntracks = Tracks->size();
   if (ntracks > 1)
     edm::LogVerbatim("AlignmenStats") << "~~~~~~~~~~~~\n For this event processing " << ntracks << " tracks";
 
   unsigned int trk_cnt = 0;
 
   for (int j = 0; j < MAXTRKS_; j++) {
     Eta[j] = -9999.0;
     Phi[j] = -8888.0;
     P[j] = -7777.0;
     Pt[j] = -6666.0;
     Chi2n[j] = -2222.0;
     for (int k = 0; k < 7; k++) {
       Nhits[j][k] = 0;
     }
   }
 
   // int npxbhits=0;
 
   //loop on tracks
   for (const auto &ittrk : *Tracks) {
     Eta[trk_cnt] = ittrk.eta();
     Phi[trk_cnt] = ittrk.phi();
     Chi2n[trk_cnt] = ittrk.normalizedChi2();
     P[trk_cnt] = ittrk.p();
     Pt[trk_cnt] = ittrk.pt();
     Nhits[trk_cnt][0] = ittrk.numberOfValidHits();
 
     if (ntracks > 1)
       edm::LogVerbatim("AlignmenStats") << "Track #" << trk_cnt + 1 << " params:    Eta=" << Eta[trk_cnt]
                                         << "  Phi=" << Phi[trk_cnt] << "  P=" << P[trk_cnt]
                                         << "   Nhits=" << Nhits[trk_cnt][0];
 
     //loop on tracking rechits
     //edm::LogVerbatim("AlignmenStats") << "   loop on hits of track #" << (itt - tracks->begin());
     for (auto const &hit : ittrk.recHits()) {
       if (!hit->isValid())
         continue;
       DetId detid = hit->geographicalId();
       int subDet = detid.subdetId();
       uint32_t rawId = hit->geographicalId().rawId();
 
       //  if(subDet==1)npxbhits++;
 
       //look if you find this detid in the map
       DetHitMap::iterator mapiter;
       mapiter = hitmap_.find(rawId);
       if (mapiter != hitmap_.end()) {  //present, increase its value by one
         //  hitmap_[rawId]=hitmap_[rawId]+1;
         ++(hitmap_[rawId]);
       } else {  //not present, let's add this key to the map with value=1
         hitmap_.insert(pair<uint32_t, uint32_t>(rawId, 1));
       }
 
       AlignmentClusterFlag inval;
 
       bool hitInPixel = (subDet == PixelSubdetector::PixelBarrel) || (subDet == PixelSubdetector::PixelEndcap);
       bool hitInStrip = (subDet == SiStripDetId::TIB) || (subDet == SiStripDetId::TID) ||
                         (subDet == SiStripDetId::TOB) || (subDet == SiStripDetId::TEC);
 
       if (!(hitInPixel || hitInStrip)) {
         //skip only this hit, don't stop everything throwing an exception
         edm::LogError("AlignmentStats") << "Hit not belonging neither to pixels nor to strips! Skipping it. SubDet= "
                                         << subDet;
         continue;
       }
 
       //check also if the hit is an overlap. If yes fill a dedicated hitmap
       if (hitInStrip) {
         const std::type_info &type = typeid(*hit);
 
         if (type == typeid(SiStripRecHit1D)) {
           //Notice the difference respect to when one loops on Trajectories: the recHit is a TrackingRecHit and not a TransientTrackingRecHit
           const SiStripRecHit1D *striphit = dynamic_cast<const SiStripRecHit1D *>(hit);
           if (striphit != nullptr) {
             SiStripRecHit1D::ClusterRef stripclust(striphit->cluster());
             inval = OverlapMap[stripclust];
           } else {
             //  edm::LogError("AlignmentStats")<<"ERROR in <AlignmentStats::analyze>: Dynamic cast of Strip RecHit1D failed!   TypeId of the RecHit: "<<className(*hit);
             throw cms::Exception("NullPointerError")
                 << "ERROR in <AlignmentStats::analyze>: Dynamic cast of Strip RecHit1D failed!   TypeId of the RecHit: "
                 << className(*hit);
           }
         }  //end if sistriprechit1D
         if (type == typeid(SiStripRecHit2D)) {
           const SiStripRecHit2D *striphit = dynamic_cast<const SiStripRecHit2D *>(hit);
           if (striphit != nullptr) {
             SiStripRecHit2D::ClusterRef stripclust(striphit->cluster());
             inval = OverlapMap[stripclust];
             //edm::LogVerbatim("AlignmenStats")<<"Taken the Strip Cluster with ProdId "<<stripclust.id() <<"; the Value in the map is "<<inval<<"  (DetId is "<<hit->geographicalId().rawId()<<")";
           } else {
             throw cms::Exception("NullPointerError")
                 << "ERROR in <AlignmentStats::analyze>: Dynamic cast of Strip RecHit2D failed!   TypeId of the RecHit: "
                 << className(*hit);
             //    edm::LogError("AlignmentStats")<<"ERROR in <AlignmentStats::analyze>: Dynamic cast of Strip RecHit2D failed!   TypeId of the RecHit: "<<className(*hit);
           }
         }  //end if sistriprechit2D
 
       }  //end if hit in Strips
       else {
         const SiPixelRecHit *pixelhit = dynamic_cast<const SiPixelRecHit *>(hit);
         if (pixelhit != nullptr) {
           SiPixelClusterRefNew pixclust(pixelhit->cluster());
           inval = OverlapMap[pixclust];
           //edm::LogVerbatim("AlignmenStats")<<"Taken the Pixel Cluster with ProdId "<<pixclust.id() <<"; the Value in the map is "<<inval<<"  (DetId is "<<hit->geographicalId().rawId()<<")";
         } else {
           edm::LogError("AlignmentStats")
               << "ERROR in <AlignmentStats::analyze>: Dynamic cast of Pixel RecHit failed!   TypeId of the RecHit: "
               << className(*hit);
         }
       }  //end else hit is in Pixel
 
       bool isOverlapHit(inval.isOverlap());
 
       if (isOverlapHit) {
         edm::LogVerbatim("AlignmenStats") << "This hit is an overlap !";
         DetHitMap::iterator overlapiter;
         overlapiter = overlapmap_.find(rawId);
 
         if (overlapiter !=
             overlapmap_.end()) {  //the det already collected at least an overlap, increase its value by one
           overlapmap_[rawId] = overlapmap_[rawId] + 1;
         } else {  //first overlap on det unit, let's add it to the map
           overlapmap_.insert(pair<uint32_t, uint32_t>(rawId, 1));
         }
       }  //end if the hit is an overlap
 
       int subdethit = static_cast<int>(hit->geographicalId().subdetId());
       if (ntracks > 1)
         edm::LogVerbatim("AlignmenStats") << "Hit in SubDet=" << subdethit;
       Nhits[trk_cnt][subdethit] = Nhits[trk_cnt][subdethit] + 1;
     }  //end loop on trackingrechits
     trk_cnt++;
 
   }  //end loop on tracks
 
   //edm::LogVerbatim("AlignmenStats")<<"Total number of pixel hits is " << npxbhits;
 
   tmpPresc_--;
   if (tmpPresc_ < 1) {
     outtree_->Fill();
     tmpPresc_ = prescale_;
   }
   if (trk_cnt != ntracks)
     edm::LogError("AlignmentStats") << "\nERROR! trk_cnt=" << trk_cnt << "   ntracks=" << ntracks;
 
   return;
 }

◆ beginJob()

void AlignmentStats::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 162 of file AlignmentStats.cc.

References Chi2n, Eta, event_, hitsfile_, hitstree_, hitsTreeName_, id_, is2D_, isStereo_, layer_, Nhits, nhits_, noverlaps_, ntracks, outtree_, Phi, posEta_, posPhi_, posR_, posX_, posY_, posZ_, prescale_, Pt, run_, statsTreeName_, subdet_, tmpPresc_, and treefile_.

                               {  // const edm::EventSetup &iSetup
 
   //book track stats tree
   treefile_ = new TFile(statsTreeName_.c_str(), "RECREATE");
   treefile_->cd();
   outtree_ = new TTree("AlignmentTrackStats", "Statistics of Tracks used for Alignment");
   // int nHitsinPXB[MAXTRKS_], nHitsinPXE[MAXTRKS_], nHitsinTEC[MAXTRKS_], nHitsinTIB[MAXTRKS_],nHitsinTOB[MAXTRKS_],nHitsinTID[MAXTRKS_];
 
   outtree_->Branch("Ntracks", &ntracks, "Ntracks/i");
   outtree_->Branch("Run_", &run_, "Run_Nr/I");
   outtree_->Branch("Event", &event_, "EventNr/I");
   outtree_->Branch("Eta", &Eta, "Eta[Ntracks]/F");
   outtree_->Branch("Phi", &Phi, "Phi[Ntracks]/F");
   outtree_->Branch("P", &P, "P[Ntracks]/F");
   outtree_->Branch("Pt", &Pt, "Pt[Ntracks]/F");
   outtree_->Branch("Chi2n", &Chi2n, "Chi2n[Ntracks]/F");
   outtree_->Branch("Nhits", &Nhits, "Nhits[Ntracks][7]/I");
   /*
     outtree_->Branch("NhitsPXB"       , ,);
     outtree_->Branch("NhitsPXE"       , ,);
     outtree_->Branch("NhitsTIB"       , ,);
     outtree_->Branch("NhitsTID"       , ,);
     outtree_->Branch("NhitsTOB"       , ,);
     outtree_->Branch("NhitsTOB"       , ,);
   */
 
   tmpPresc_ = prescale_;
 
   // create tree with hit maps (hitstree)
   // book hits stats tree
   hitsfile_ = new TFile(hitsTreeName_.c_str(), "RECREATE");
   hitsfile_->cd();
   hitstree_ = new TTree("AlignmentHitMap", "Maps of Hits used for Alignment");
   hitstree_->Branch("DetId", &id_, "DetId/i");
   hitstree_->Branch("Nhits", &nhits_, "Nhits/i");
   hitstree_->Branch("Noverlaps", &noverlaps_, "Noverlaps/i");
   hitstree_->Branch("SubDet", &subdet_, "SubDet/I");
   hitstree_->Branch("Layer", &layer_, "Layer/i");
   hitstree_->Branch("is2D", &is2D_, "is2D/B");
   hitstree_->Branch("isStereo", &isStereo_, "isStereo/B");
   hitstree_->Branch("posX", &posX_, "posX/F");
   hitstree_->Branch("posY", &posY_, "posY/F");
   hitstree_->Branch("posZ", &posZ_, "posZ/F");
   hitstree_->Branch("posR", &posR_, "posR/F");
   hitstree_->Branch("posEta", &posEta_, "posEta/F");
   hitstree_->Branch("posPhi", &posPhi_, "posPhi/F");
 }  //end beginJob

◆ beginRun()

void AlignmentStats::beginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

inlineoverride

Definition at line 75 of file AlignmentStats.cc.

75 {};

◆ endJob()

void AlignmentStats::endJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 442 of file AlignmentStats.cc.

References geomInfoList_, hitmap_, hitsfile_, hitstree_, id_, is2D_, isStereo_, layer_, nhits_, noverlaps_, outtree_, overlapmap_, posEta_, posPhi_, posR_, posX_, posY_, posZ_, subdet_, and treefile_.

Referenced by o2olib.O2ORunMgr::executeJob().

                             {
   treefile_->cd();
   edm::LogInfo("AlignmentStats") << "Writing out the TrackStatistics in " << gDirectory->GetPath();
   outtree_->Write();
   delete outtree_;
 
   int ndetunits = geomInfoList_.size();
   edm::LogInfo("AlignmentStats") << __PRETTY_FUNCTION__ << "Number of DetUnits in the AlignableTracker: " << ndetunits
                                  << std::endl;
 
   hitsfile_->cd();
   // save the information on hits and overlaps
   for (const auto &detUnitInfo : geomInfoList_) {
     detUnitInfo.printAll();
 
     // copy the values from the struct to the TTree
     id_ = detUnitInfo.id_;
     posX_ = detUnitInfo.posX_;
     posY_ = detUnitInfo.posY_;
     posZ_ = detUnitInfo.posZ_;
     posEta_ = detUnitInfo.posEta_;
     posPhi_ = detUnitInfo.posPhi_;
     posR_ = detUnitInfo.posR_;
     subdet_ = detUnitInfo.subdet_;
     layer_ = detUnitInfo.layer_;
     is2D_ = detUnitInfo.is2D_;
     isStereo_ = detUnitInfo.isStereo_;
 
     if (hitmap_.find(id_) != hitmap_.end()) {
       nhits_ = hitmap_[id_];
     }
     //if not, save nhits=0
     else {
       nhits_ = 0;
       hitmap_.insert(pair<uint32_t, uint32_t>(id_, 0));
     }
 
     if (overlapmap_.find(id_) != overlapmap_.end()) {
       noverlaps_ = overlapmap_[id_];
     }
     //if not, save nhits=0
     else {
       noverlaps_ = 0;
       overlapmap_.insert(pair<uint32_t, uint32_t>(id_, 0));
     }
     //write in the hitstree
     hitstree_->Fill();
   }  //end loop over detunits
 
   //save hitstree
   hitstree_->Write();
   delete hitstree_;
   hitmap_.clear();
   overlapmap_.clear();
   delete hitsfile_;
 }

◆ endRun()

void AlignmentStats::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

override

Definition at line 211 of file AlignmentStats.cc.

References esTokenTkGeoER_, esTokenTTopoER_, PV3DBase< T, PVType, FrameType >::eta(), geomInfoList_, edm::EventSetup::getData(), hitsfile_, aliStats::GeoInfo::id_, aliStats::GeoInfo::is2D_, aliStats::GeoInfo::isStereo_, aliStats::GeoInfo::layer_, LogDebug, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, aliStats::GeoInfo::posEta_, aliStats::GeoInfo::posPhi_, aliStats::GeoInfo::posR_, aliStats::GeoInfo::posX_, aliStats::GeoInfo::posY_, aliStats::GeoInfo::posZ_, DetId::rawId(), aliStats::GeoInfo::subdet_, DetId::subdetId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, and trackerTopology_.

                                                                            {
   if (!trackerTopology_) {
     trackerTopology_ = std::make_unique<TrackerTopology>(iSetup.getData(esTokenTTopoER_));
     const TrackerGeometry *trackerGeometry_ = &iSetup.getData(esTokenTkGeoER_);
     auto theAliTracker = std::make_unique<AlignableTracker>(trackerGeometry_, trackerTopology_.get());
 
     hitsfile_->cd();
     for (const auto &detUnit : theAliTracker->deepComponents()) {
       aliStats::GeoInfo detUnitInfo;
       detUnitInfo.id_ = static_cast<uint32_t>(detUnit->id());
       DetId detid(detUnitInfo.id_);
       detUnitInfo.subdet_ = detid.subdetId();
 
       //take other geometrical infos from the det
       detUnitInfo.posX_ = detUnit->globalPosition().x();
       detUnitInfo.posY_ = detUnit->globalPosition().y();
       detUnitInfo.posZ_ = detUnit->globalPosition().z();
 
       align::GlobalVector vec(detUnitInfo.posX_, detUnitInfo.posY_, detUnitInfo.posZ_);
       detUnitInfo.posR_ = vec.perp();
       detUnitInfo.posPhi_ = vec.phi();
       detUnitInfo.posEta_ = vec.eta();
       // detUnitInfo.posPhi_ = atan2(posY_,posX_);
 
       //get layers, petals, etc...
       if (detUnitInfo.subdet_ == PixelSubdetector::PixelBarrel) {  //PXB
         detUnitInfo.layer_ = trackerTopology_->pxbLayer(detUnitInfo.id_);
         detUnitInfo.is2D_ = true;
         detUnitInfo.isStereo_ = false;
       } else if (detUnitInfo.subdet_ == PixelSubdetector::PixelEndcap) {
         detUnitInfo.layer_ = trackerTopology_->pxfDisk(detUnitInfo.id_);
         detUnitInfo.is2D_ = true;
         detUnitInfo.isStereo_ = false;
       } else if (detUnitInfo.subdet_ == SiStripDetId::TIB) {
         detUnitInfo.layer_ = trackerTopology_->tibLayer(detUnitInfo.id_);
         detUnitInfo.is2D_ = trackerTopology_->tibIsDoubleSide(detUnitInfo.id_);
         detUnitInfo.isStereo_ = trackerTopology_->tibIsStereo(detUnitInfo.id_);
       } else if (detUnitInfo.subdet_ == SiStripDetId::TID) {
         detUnitInfo.layer_ = trackerTopology_->tidWheel(detUnitInfo.id_);
         detUnitInfo.is2D_ = trackerTopology_->tidIsDoubleSide(detUnitInfo.id_);
         detUnitInfo.isStereo_ = trackerTopology_->tidIsStereo(detUnitInfo.id_);
       } else if (detUnitInfo.subdet_ == SiStripDetId::TOB) {
         detUnitInfo.layer_ = trackerTopology_->tobLayer(detUnitInfo.id_);
         detUnitInfo.is2D_ = trackerTopology_->tobIsDoubleSide(detUnitInfo.id_);
         detUnitInfo.isStereo_ = trackerTopology_->tobIsStereo(detUnitInfo.id_);
       } else if (detUnitInfo.subdet_ == SiStripDetId::TEC) {
         detUnitInfo.layer_ = trackerTopology_->tecWheel(detUnitInfo.id_);
         detUnitInfo.is2D_ = trackerTopology_->tecIsDoubleSide(detUnitInfo.id_);
         detUnitInfo.isStereo_ = trackerTopology_->tecIsStereo(detUnitInfo.id_);
       } else {
         edm::LogError("AlignmentStats")
             << "Detector not belonging neither to pixels nor to strips! Skipping it. SubDet= " << detUnitInfo.subdet_;
       }
 
       LogDebug("AlignmentStats") << "id " << detUnitInfo.id_ << " detid.rawId()" << detid.rawId() << " subdet "
                                  << detUnitInfo.subdet_;
 
       // push back in the list
       geomInfoList_.push_back(detUnitInfo);
     }  // end loop over detunits
 
     int ndetunits = geomInfoList_.size();
     edm::LogInfo("AlignmentStats") << __PRETTY_FUNCTION__
                                    << " Number of DetUnits in the AlignableTracker: " << ndetunits;
   }
 }

◆ fillDescriptions()

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

static

Definition at line 150 of file AlignmentStats.cc.

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

                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("src", edm::InputTag("generalTracks"));
   desc.add<edm::InputTag>("OverlapAssoMap", edm::InputTag("OverlapAssoMap"));
   desc.add<bool>("keepTrackStats", false);
   desc.add<bool>("keepHitStats", false);
   desc.add<std::string>("TrkStatsFileName", "tracks_statistics.root");
   desc.add<std::string>("HitStatsFileName", "HitMaps.root");
   desc.add<unsigned int>("TrkStatsPrescale", 1);
   descriptions.add("AlignmentStats", desc);
 }

Member Data Documentation

◆ Chi2n

float AlignmentStats::Chi2n[MAXTRKS_]

private

Definition at line 109 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ esTokenTkGeoER_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> AlignmentStats::esTokenTkGeoER_

private

Definition at line 86 of file AlignmentStats.cc.

Referenced by endRun().

◆ esTokenTTopoER_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> AlignmentStats::esTokenTTopoER_

private

Definition at line 85 of file AlignmentStats.cc.

Referenced by endRun().

◆ Eta

float AlignmentStats::Eta[MAXTRKS_]

private

Definition at line 109 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ event_

int AlignmentStats::event_

private

Definition at line 107 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ geomInfoList_

std::vector<aliStats::GeoInfo> AlignmentStats::geomInfoList_

private

Definition at line 82 of file AlignmentStats.cc.

Referenced by endJob(), and endRun().

◆ hitmap_

DetHitMap AlignmentStats::hitmap_

private

Definition at line 123 of file AlignmentStats.cc.

Referenced by analyze(), and endJob().

◆ hitsfile_

TFile* AlignmentStats::hitsfile_

private

Definition at line 113 of file AlignmentStats.cc.

Referenced by beginJob(), endJob(), and endRun().

◆ hitstree_

TTree* AlignmentStats::hitstree_

private

Definition at line 114 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ hitsTreeName_

const std::string AlignmentStats::hitsTreeName_

private

Definition at line 94 of file AlignmentStats.cc.

Referenced by beginJob().

◆ id_

unsigned int AlignmentStats::id_

private

Definition at line 115 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ is2D_

bool AlignmentStats::is2D_

private

Definition at line 120 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ isStereo_

bool AlignmentStats::isStereo_

private

Definition at line 120 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ keepHitPopulation_

const bool AlignmentStats::keepHitPopulation_

private

Definition at line 92 of file AlignmentStats.cc.

◆ keepTrackStats_

const bool AlignmentStats::keepTrackStats_

private

Definition at line 91 of file AlignmentStats.cc.

◆ layer_

unsigned int AlignmentStats::layer_

private

Definition at line 119 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ mapToken_

const edm::EDGetTokenT<AliClusterValueMap> AlignmentStats::mapToken_

private

Definition at line 98 of file AlignmentStats.cc.

Referenced by analyze().

◆ MAXTRKS_

const int AlignmentStats::MAXTRKS_ = 200

staticprivate

Definition at line 106 of file AlignmentStats.cc.

Referenced by analyze().

◆ Nhits

int AlignmentStats::Nhits[MAXTRKS_][7]

private

Definition at line 110 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ nhits_

unsigned int AlignmentStats::nhits_

private

Definition at line 115 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ noverlaps_

unsigned int AlignmentStats::noverlaps_

private

Definition at line 115 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ ntracks

unsigned int AlignmentStats::ntracks

private

Definition at line 108 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ outtree_

TTree* AlignmentStats::outtree_

private

Definition at line 105 of file AlignmentStats.cc.

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

◆ overlapAM_

const edm::InputTag AlignmentStats::overlapAM_

private

Definition at line 90 of file AlignmentStats.cc.

◆ overlapmap_

DetHitMap AlignmentStats::overlapmap_

private

Definition at line 124 of file AlignmentStats.cc.

Referenced by analyze(), and endJob().

◆ P

float AlignmentStats::P[MAXTRKS_]

private

Definition at line 109 of file AlignmentStats.cc.

◆ Phi

float AlignmentStats::Phi[MAXTRKS_]

private

Definition at line 109 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ posEta_

float AlignmentStats::posEta_

private

Definition at line 117 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ posPhi_

float AlignmentStats::posPhi_

private

Definition at line 117 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ posR_

float AlignmentStats::posR_

private

Definition at line 117 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ posX_

float AlignmentStats::posX_

private

Definition at line 116 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ posY_

float AlignmentStats::posY_

private

Definition at line 116 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ posZ_

float AlignmentStats::posZ_

private

Definition at line 116 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ prescale_

const uint32_t AlignmentStats::prescale_

private

Definition at line 95 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ Pt

float AlignmentStats::Pt[MAXTRKS_]

private

Definition at line 109 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ run_

int AlignmentStats::run_

private

Definition at line 107 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ src_

const edm::InputTag AlignmentStats::src_

private

Definition at line 89 of file AlignmentStats.cc.

◆ statsTreeName_

const std::string AlignmentStats::statsTreeName_

private

Definition at line 93 of file AlignmentStats.cc.

Referenced by beginJob().

◆ subdet_

int AlignmentStats::subdet_

private

Definition at line 118 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

◆ tmpPresc_

uint32_t AlignmentStats::tmpPresc_

private

Definition at line 101 of file AlignmentStats.cc.

Referenced by analyze(), and beginJob().

◆ trackerTopology_

std::unique_ptr<TrackerTopology> AlignmentStats::trackerTopology_

private

Definition at line 126 of file AlignmentStats.cc.

Referenced by endRun().

◆ trackToken_

const edm::EDGetTokenT<reco::TrackCollection> AlignmentStats::trackToken_

private

Definition at line 97 of file AlignmentStats.cc.

Referenced by analyze().

◆ treefile_

TFile* AlignmentStats::treefile_

private

Definition at line 104 of file AlignmentStats.cc.

Referenced by beginJob(), and endJob().

Public Member Functions

Static Public Member Functions

Private Types

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ DetHitMap

Constructor & Destructor Documentation

◆ AlignmentStats()

◆ ~AlignmentStats()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ fillDescriptions()

Member Data Documentation

◆ Chi2n

◆ esTokenTkGeoER_

◆ esTokenTTopoER_

◆ Eta

◆ event_

◆ geomInfoList_

◆ hitmap_

◆ hitsfile_

◆ hitstree_

◆ hitsTreeName_

◆ id_

◆ is2D_

◆ isStereo_

◆ keepHitPopulation_

◆ keepTrackStats_

◆ layer_

◆ mapToken_

◆ MAXTRKS_

◆ Nhits

◆ nhits_

◆ noverlaps_

◆ ntracks

◆ outtree_

◆ overlapAM_

◆ overlapmap_

◆ P

◆ Phi

◆ posEta_

◆ posPhi_

◆ posR_

◆ posX_

◆ posY_

◆ posZ_

◆ prescale_

◆ Pt

◆ run_

◆ src_

◆ statsTreeName_

◆ subdet_

◆ tmpPresc_

◆ trackerTopology_

◆ trackToken_

◆ treefile_