#include <AlignmentStats.h>

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	endJob () override

	~AlignmentStats () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef std::map< uint32_t, uint32_t >	DetHitMap

Private Attributes
float	Chi2n [MAXTRKS_]

const edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	esTokenTkGeo_

const edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	esTokenTTopo_

float	Eta [MAXTRKS_]

int	event_

DetHitMap	hitmap_

std::string	hitsTreeName_

bool	keepHitPopulation_

bool	keepTrackStats_

int	Nhits [MAXTRKS_][7]

unsigned int	ntracks

TTree *	outtree_

edm::InputTag	overlapAM_

DetHitMap	overlapmap_

float	P [MAXTRKS_]

float	Phi [MAXTRKS_]

uint32_t	prescale_

float	Pt [MAXTRKS_]

int	run_

edm::InputTag	src_

std::string	statsTreeName_

uint32_t	tmpPresc_

std::unique_ptr< TrackerGeometry >	trackerGeometry_

std::unique_ptr< TrackerTopology >	trackerTopology_

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

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

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 28 of file AlignmentStats.h.

Member Typedef Documentation

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

private

Definition at line 62 of file AlignmentStats.h.

Constructor & Destructor Documentation

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

Definition at line 28 of file AlignmentStats.cc.

References outtree_.

     : esTokenTTopo_(esConsumes()),
       esTokenTkGeo_(esConsumes()),
       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")) {
   //sanity checks
 
   //init
   outtree_ = nullptr;
 
 }  //end constructor

AlignmentStats::~AlignmentStats ( )

override

Definition at line 45 of file AlignmentStats.cc.

                                 {
   //
 }

Member Function Documentation

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 78 of file AlignmentStats.cc.

References Chi2n, className(), SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), SiPixelRecHit::cluster(), esTokenTkGeo_, esTokenTTopo_, Eta, edm::EventID::event(), event_, Exception, edm::Event::getByLabel(), edm::EventSetup::getData(), hitmap_, edm::EventBase::id(), AlignmentClusterFlag::isOverlap(), dqmiolumiharvest::j, isotrackApplyRegressor::k, MAXTRKS_, Nhits, ntracks, outtree_, overlapAM_, overlapmap_, Phi, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, prescale_, Pt, edm::EventID::run(), run_, src_, DetId::subdetId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, tmpPresc_, SiStripDetId::TOB, trackerGeometry_, trackerTopology_, and HLT_FULL_cff::Tracks.

                                                                                 {
   //load list of detunits needed then in endJob
   if (!trackerGeometry_) {
     trackerGeometry_ = std::make_unique<TrackerGeometry>(iSetup.getData(esTokenTkGeo_));
   }
 
   if (!trackerTopology_) {
     trackerTopology_ = std::make_unique<TrackerTopology>(iSetup.getData(esTokenTTopo_));
   }
 
   //take trajectories and tracks to loop on
   // edm::Handle<TrajTrackAssociationCollection> TrackAssoMap;
   edm::Handle<reco::TrackCollection> Tracks;
   iEvent.getByLabel(src_, Tracks);
 
   //take overlap HitAssomap
   edm::Handle<AliClusterValueMap> hMap;
   iEvent.getByLabel(overlapAM_, hMap);
   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) std::cout<<"~~~~~~~~~~~~\n For this event processing "<<ntracks<<" tracks"<<std::endl;
 
   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 (std::vector<reco::Track>::const_iterator ittrk = Tracks->begin(), edtrk = Tracks->end(); ittrk != edtrk;
        ++ittrk) {
     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)std::cout<<"Track #"<<trk_cnt+1<<" params:    Eta="<< Eta[trk_cnt]<<"  Phi="<< Phi[trk_cnt]<<"  P="<<P[trk_cnt]<<"   Nhits="<<Nhits[trk_cnt][0]<<std::endl;
 
     int nhit = 0;
     //loop on tracking rechits
     //std::cout << "   loop on hits of track #" << (itt - tracks->begin()) << std::endl;
     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];
             //cout<<"Taken the Strip Cluster with ProdId "<<stripclust.id() <<"; the Value in the map is "<<inval<<"  (DetId is "<<hit->geographicalId().rawId()<<")"<<endl;
           } 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];
           //cout<<"Taken the Pixel Cluster with ProdId "<<pixclust.id() <<"; the Value in the map is "<<inval<<"  (DetId is "<<hit->geographicalId().rawId()<<")"<<endl;
         } 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) {
         //cout<<"This hit is an overlap !"<<endl;
         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)std::cout<<"Hit in SubDet="<<subdethit<<std::endl;
       Nhits[trk_cnt][subdethit] = Nhits[trk_cnt][subdethit] + 1;
       nhit++;
     }  //end loop on trackingrechits
     trk_cnt++;
 
   }  //end loop on tracks
 
   //  cout<<"Total number of pixel hits is "<<npxbhits<<endl;
 
   tmpPresc_--;
   if (tmpPresc_ < 1) {
     outtree_->Fill();
     tmpPresc_ = prescale_;
   }
   if (trk_cnt != ntracks)
     edm::LogError("AlignmentStats") << "\nERROR! trk_cnt=" << trk_cnt << "   ntracks=" << ntracks;
 
   return;
 }

void AlignmentStats::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 49 of file AlignmentStats.cc.

References Chi2n, Eta, event_, Nhits, ntracks, outtree_, Phi, prescale_, Pt, run_, statsTreeName_, 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_;
 }  //end beginJob

void AlignmentStats::endJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 251 of file AlignmentStats.cc.

References PV3DBase< T, PVType, FrameType >::eta(), hitmap_, hitsTreeName_, gpuClustering::id, fastTrackerRecHitType::is2D(), phase1PixelTopology::layer, nhits, outtree_, overlapmap_, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::subdetId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, trackerGeometry_, trackerTopology_, and treefile_.

Referenced by o2olib.O2ORunMgr::executeJob().

                             {
   treefile_->cd();
   edm::LogInfo("AlignmentStats") << "Writing out the TrackStatistics in " << gDirectory->GetPath() << std::endl;
   outtree_->Write();
   delete outtree_;
 
   //create tree with hit maps (hitstree)
   //book track stats tree
   TFile *hitsfile = new TFile(hitsTreeName_.c_str(), "RECREATE");
   hitsfile->cd();
   TTree *hitstree = new TTree("AlignmentHitMap", "Maps of Hits used for Alignment");
 
   unsigned int id = 0, nhits = 0, noverlaps = 0;
   float posX(-99999.0), posY(-77777.0), posZ(-88888.0);
   float posEta(-6666.0), posPhi(-5555.0), posR(-4444.0);
   int subdet = 0;
   unsigned int layer = 0;
   bool is2D = false, isStereo = false;
   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");
 
   /*
   TTree *overlapstree=new TTree("OverlapHitMap","Maps of Overlaps used for Alignment");
   hitstree->Branch("DetId",   &id ,     "DetId/i");
   hitstree->Branch("NOverlaps",   &nhits ,  "Nhits/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");
   */
 
   std::unique_ptr<AlignableTracker> theAliTracker =
       std::make_unique<AlignableTracker>(trackerGeometry_.get(), trackerTopology_.get());
   const auto &Detunitslist = theAliTracker->deepComponents();
   int ndetunits = Detunitslist.size();
   edm::LogInfo("AlignmentStats") << "Number of DetUnits in the AlignableTracker: " << ndetunits << std::endl;
 
   for (int det_cnt = 0; det_cnt < ndetunits; ++det_cnt) {
     //re-initialize for safety
     id = 0;
     nhits = 0;
     noverlaps = 0;
     posX = -99999.0;
     posY = -77777.0;
     posZ = -88888.0;
     posEta = -6666.0;
     posPhi = -5555.0;
     posR = -4444.0;
     subdet = 0;
     layer = 0;
     is2D = false;
     isStereo = false;
 
     //if detunit in vector is found also in the map, look for how many hits were collected
     //and save in the tree this number
     id = static_cast<uint32_t>(Detunitslist[det_cnt]->id());
     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));
     }
 
     //take other geometrical infos from the det
     posX = Detunitslist[det_cnt]->globalPosition().x();
     posY = Detunitslist[det_cnt]->globalPosition().y();
     posZ = Detunitslist[det_cnt]->globalPosition().z();
 
     align::GlobalVector vec(posX, posY, posZ);
     posR = vec.perp();
     posPhi = vec.phi();
     posEta = vec.eta();
     //   posPhi = atan2(posY,posX);
 
     DetId detid(id);
     subdet = detid.subdetId();
 
     //get layers, petals, etc...
     if (subdet == PixelSubdetector::PixelBarrel) {  //PXB
 
       layer = trackerTopology_->pxbLayer(id);
       is2D = true;
       isStereo = false;
     } else if (subdet == PixelSubdetector::PixelEndcap) {
       layer = trackerTopology_->pxfDisk(id);
       is2D = true;
       isStereo = false;
     } else if (subdet == SiStripDetId::TIB) {
       layer = trackerTopology_->tibLayer(id);
       is2D = trackerTopology_->tibIsDoubleSide(id);
       isStereo = trackerTopology_->tibIsStereo(id);
     } else if (subdet == SiStripDetId::TID) {
       layer = trackerTopology_->tidWheel(id);
       is2D = trackerTopology_->tidIsDoubleSide(id);
       isStereo = trackerTopology_->tidIsStereo(id);
     } else if (subdet == SiStripDetId::TOB) {
       layer = trackerTopology_->tobLayer(id);
       is2D = trackerTopology_->tobIsDoubleSide(id);
       isStereo = trackerTopology_->tobIsStereo(id);
     } else if (subdet == SiStripDetId::TEC) {
       layer = trackerTopology_->tecWheel(id);
       is2D = trackerTopology_->tecIsDoubleSide(id);
       isStereo = trackerTopology_->tecIsStereo(id);
     } else {
       edm::LogError("AlignmentStats") << "Detector not belonging neither to pixels nor to strips! Skipping it. SubDet= "
                                       << subdet;
     }
 
     //write in the hitstree
     hitstree->Fill();
   }  //end loop over detunits
 
   //save hitstree
   hitstree->Write();
   delete hitstree;
   //delete Detunitslist;
   hitmap_.clear();
   overlapmap_.clear();
   delete hitsfile;
 }

Member Data Documentation

float AlignmentStats::Chi2n[MAXTRKS_]

private

Definition at line 58 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

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

private

Definition at line 39 of file AlignmentStats.h.

Referenced by analyze().

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

private

Definition at line 38 of file AlignmentStats.h.

Referenced by analyze().

float AlignmentStats::Eta[MAXTRKS_]

private

Definition at line 58 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

int AlignmentStats::event_

private

Definition at line 56 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

DetHitMap AlignmentStats::hitmap_

private

Definition at line 63 of file AlignmentStats.h.

Referenced by analyze(), and endJob().

std::string AlignmentStats::hitsTreeName_

private

Definition at line 47 of file AlignmentStats.h.

Referenced by endJob().

bool AlignmentStats::keepHitPopulation_

private

Definition at line 45 of file AlignmentStats.h.

bool AlignmentStats::keepTrackStats_

private

Definition at line 44 of file AlignmentStats.h.

const int AlignmentStats::MAXTRKS_ = 200

staticprivate

Definition at line 55 of file AlignmentStats.h.

Referenced by analyze().

int AlignmentStats::Nhits[MAXTRKS_][7]

private

Definition at line 59 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

unsigned int AlignmentStats::ntracks

private

Definition at line 57 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

TTree* AlignmentStats::outtree_

private

Definition at line 54 of file AlignmentStats.h.

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

edm::InputTag AlignmentStats::overlapAM_

private

Definition at line 43 of file AlignmentStats.h.

Referenced by analyze().

DetHitMap AlignmentStats::overlapmap_

private

Definition at line 64 of file AlignmentStats.h.

Referenced by analyze(), and endJob().

float AlignmentStats::P[MAXTRKS_]

private

Definition at line 58 of file AlignmentStats.h.

float AlignmentStats::Phi[MAXTRKS_]

private

Definition at line 58 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

uint32_t AlignmentStats::prescale_

private

Definition at line 48 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

float AlignmentStats::Pt[MAXTRKS_]

private

Definition at line 58 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

int AlignmentStats::run_

private

Definition at line 56 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

edm::InputTag AlignmentStats::src_

private

Definition at line 42 of file AlignmentStats.h.

Referenced by analyze().

std::string AlignmentStats::statsTreeName_

private

Definition at line 46 of file AlignmentStats.h.

Referenced by beginJob().

uint32_t AlignmentStats::tmpPresc_

private

Definition at line 50 of file AlignmentStats.h.

Referenced by analyze(), and beginJob().

std::unique_ptr<TrackerGeometry> AlignmentStats::trackerGeometry_

private

Definition at line 67 of file AlignmentStats.h.

Referenced by analyze(), and endJob().

std::unique_ptr<TrackerTopology> AlignmentStats::trackerTopology_

private

Definition at line 66 of file AlignmentStats.h.

Referenced by analyze(), and endJob().

TFile* AlignmentStats::treefile_

private

Definition at line 53 of file AlignmentStats.h.

Referenced by beginJob(), and endJob().

Public Member Functions

Private Types

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation