SiStripLAProfileBooker Class Reference

#include <SiStripLAProfileBooker.h>

Inheritance diagram for SiStripLAProfileBooker:

Classes
struct	detparameters

Public Types
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginRun (edm::Run const &, const edm::EventSetup &c) override
void	endJob () override
void	endRun (edm::Run const &, const edm::EventSetup &c) override
void	getlayer (const DetId &detid, const TrackerTopology *tTopo, std::string &name, unsigned int &layerid)
	SiStripLAProfileBooker (const edm::ParameterSet &conf)
	~SiStripLAProfileBooker () override
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)

Private Types
typedef std::map< unsigned int, detparameters * >	detparmap
typedef std::map< unsigned int, MonitorElement * >	histomap

Private Attributes
uint8_t	Amplitudes [100]
float	barycenter
int	better_double_hit
int	bw_fw
float	chi2norm
int	ClSize
edm::ParameterSet	conf_
DQMStore *	dbe_
edm::ESGetToken< SiStripDetCabling, SiStripDetCablingRcd >	detCablingToken_
detparmap	detmap
float	EtaTrack
TDirectory *	Event_Tree
int	eventcounter
int	EventCounter
int	EventNumber
TTree *	EventTree
int	Ext_Int
TFile *	hFile
histomap	histos
float	hit_std_dev
TDirectory *	Hit_Tree
int	HitCharge
int	hitcounter
int	hitcounter_2ndloop
int	HitNr
int	HitPerTrack
TTree *	HitsTree
int	id_detector
int	Layer
float	MagField
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_
float	Momentum
int	MonoStereo
int	nstrip
int	ParticleCharge
float	PhiTrack
float	pitch_detector
float	pt
int	RunNumber
float	SignCorrection
int	size
detparmap	summarydetmap
histomap	summaryhisto
float	sumx
float	TanTrackAngle
float	TanTrackAngleParallel
float	thick_detector
const TrackerGeometry *	tkGeom_ = nullptr
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_
TDirectory *	Track_Tree
int	trackcollsize
int	TrackCounter
TTree *	TrackTree
int	trajsize
std::string	treename_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_
int	Type
int	Wheel
int	worse_double_hit
float	XGlobal
float	YGlobal
float	ZGlobal

Additional Inherited Members
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 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 33 of file SiStripLAProfileBooker.h.

Member Typedef Documentation

detparmap

typedef std::map<unsigned int, detparameters*> SiStripLAProfileBooker::detparmap

private

Definition at line 57 of file SiStripLAProfileBooker.h.

DQMStore

typedef dqm::legacy::DQMStore SiStripLAProfileBooker::DQMStore

Definition at line 36 of file SiStripLAProfileBooker.h.

histomap

typedef std::map<unsigned int, MonitorElement*> SiStripLAProfileBooker::histomap

private

Definition at line 58 of file SiStripLAProfileBooker.h.

MonitorElement

typedef dqm::legacy::MonitorElement SiStripLAProfileBooker::MonitorElement

Definition at line 35 of file SiStripLAProfileBooker.h.

Constructor & Destructor Documentation

SiStripLAProfileBooker()

SiStripLAProfileBooker::SiStripLAProfileBooker ( const edm::ParameterSet &  conf )

explicit

Definition at line 55 of file SiStripLAProfileBooker.cc.

    : conf_(conf),
      tTopoToken_(esConsumes<edm::Transition::BeginRun>()),
      tkGeomToken_(esConsumes<edm::Transition::BeginRun>()),
      magFieldToken_(esConsumes<edm::Transition::BeginRun>()),
      detCablingToken_(conf_.getParameter<bool>("UseStripCablingDB")
                           ? decltype(detCablingToken_){esConsumes<edm::Transition::BeginRun>()}
                           : decltype(detCablingToken_){}) {}

~SiStripLAProfileBooker()

SiStripLAProfileBooker::~SiStripLAProfileBooker ( )

override

Definition at line 242 of file SiStripLAProfileBooker.cc.

References detmap, cms::cuda::for(), hFile, and summarydetmap.

                                                {
  detparmap::iterator detpariter;
  for (detpariter = detmap.begin(); detpariter != detmap.end(); ++detpariter)
    delete detpariter->second;
  for (detpariter = summarydetmap.begin(); detpariter != summarydetmap.end(); ++detpariter)
    delete detpariter->second;
  delete hFile;
}

Member Function Documentation

analyze()

void SiStripLAProfileBooker::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 253 of file SiStripLAProfileBooker.cc.

References Amplitudes, barycenter, edm::AssociationMap< Tag >::begin(), better_double_hit, Surface::bounds(), bw_fw, utilities::cache(), chi2norm, ClSize, SiStripRecHit2D::cluster(), conf_, detmap, MillePedeFileConverter_cfg::e, edm::AssociationMap< Tag >::end(), EtaTrack, eventcounter, EventNumber, EventTree, Ext_Int, TrackingRecHit::geographicalId(), edm::EventSetup::getData(), getlayer(), edm::ParameterSet::getParameter(), histos, hit_std_dev, HitCharge, hitcounter, hitcounter_2ndloop, HitNr, HitPerTrack, HitsTree, mps_fire::i, id_detector, TrackerGeometry::idToDet(), createfilelist::int, edm::helpers::KeyVal< K, V >::key, Layer, TrajectoryStateOnSurface::localDirection(), StripTopology::localPitch(), BaseTrackerRecHit::localPosition(), LogDebug, PV3DBase< T, PVType, FrameType >::mag(), MagField, Momentum, GluedGeomDet::monoDet(), SiStripMatchedRecHit2D::monoHit(), MonoStereo, Skims_PA_cff::name, nstrip, ProjectedSiStripRecHit2D::originalHit(), ParticleCharge, PhiTrack, pitch_detector, position, funct::pow(), pt, DetId::rawId(), RunNumber, SignCorrection, size, GeomDet::specificSurface(), mathSSE::sqrt(), StripSubdetector::stereo(), GluedGeomDet::stereoDet(), SiStripMatchedRecHit2D::stereoHit(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), summaryhisto, sumx, GeomDet::surface(), TanTrackAngle, TanTrackAngleParallel, StripSubdetector::TEC, TrackerTopology::tecPetalInfo(), TrackerTopology::tecWheel(), thick_detector, Bounds::thickness(), StripSubdetector::TIB, TrackerTopology::tibLayer(), TrackerTopology::tibStringInfo(), StripSubdetector::TID, TrackerTopology::tidModuleInfo(), TrackerTopology::tidWheel(), tkGeom_, Zmumu_cfg::TkTag, StripSubdetector::TOB, TrackerTopology::tobLayer(), TrackerTopology::tobRodInfo(), GeomDet::toGlobal(), GeomDet::toLocal(), GeomDet::topology(), JetHT_cfg::trackCollection, trackcollsize, DiMuonV_cfg::tracks, TrackTree, trajsize, tTopoToken_, Type, edm::helpers::KeyVal< K, V >::val, Wheel, worse_double_hit, PV3DBase< T, PVType, FrameType >::x(), XGlobal, PV3DBase< T, PVType, FrameType >::y(), YGlobal, PV3DBase< T, PVType, FrameType >::z(), and ZGlobal.

                                                                               {
  const TrackerTopology* const tTopo = &es.getData(tTopoToken_);

  RunNumber = e.id().run();
  EventNumber = e.id().event();

  eventcounter++;

  EventTree->Fill();

  //Analysis of Trajectory-RecHits

  edm::InputTag TkTag = conf_.getParameter<edm::InputTag>("Tracks");

  edm::Handle<reco::TrackCollection> trackCollection;
  e.getByLabel(TkTag, trackCollection);

  edm::Handle<std::vector<Trajectory> > TrajectoryCollection;
  e.getByLabel(TkTag, TrajectoryCollection);

  edm::Handle<TrajTrackAssociationCollection> TrajTrackMap;
  e.getByLabel(TkTag, TrajTrackMap);

  const reco::TrackCollection* tracks = trackCollection.product();

  // FIXME this has to be changed to use pointers to clusters...
  std::map<const SiStripRecHit2D*, std::pair<float, float>, DetIdLess> hitangleassociation;
  std::list<SiStripRecHit2D> cache;  // ugly, inefficient, effective in making the above working

  trackcollsize = tracks->size();
  trajsize = TrajectoryCollection->size();

  edm::LogInfo("SiStripLAProfileBooker::analyze") << " Number of tracks in event = " << trackcollsize << "\n";
  edm::LogInfo("SiStripLAProfileBooker::analyze") << " Number of trajectories in event = " << trajsize << "\n";

  TrajTrackAssociationCollection::const_iterator TrajTrackIter;

  for (TrajTrackIter = TrajTrackMap->begin(); TrajTrackIter != TrajTrackMap->end();
       TrajTrackIter++) {  //loop on trajectories

    if (TrajTrackIter->key->foundHits() >= 5) {
      TrackTree->Fill();

      ParticleCharge = -99;
      Momentum = -99;
      pt = -99;
      chi2norm = -99;
      HitPerTrack = -99;
      EtaTrack = -99;
      PhiTrack = -99;

      ParticleCharge = TrajTrackIter->val->charge();
      pt = TrajTrackIter->val->pt();
      Momentum = TrajTrackIter->val->p();
      chi2norm = TrajTrackIter->val->normalizedChi2();
      EtaTrack = TrajTrackIter->val->eta();
      PhiTrack = TrajTrackIter->val->phi();
      HitPerTrack = TrajTrackIter->key->foundHits();

      std::vector<TrajectoryMeasurement> TMeas = TrajTrackIter->key->measurements();
      std::vector<TrajectoryMeasurement>::iterator itm;

      for (itm = TMeas.begin(); itm != TMeas.end(); itm++) {  //loop on hits

        int i;
        for (i = 0; i < 100; i++) {
          Amplitudes[i] = 0;
        }

        TanTrackAngle = -99;
        TanTrackAngleParallel = -99;
        ClSize = -99;
        HitCharge = 0;
        Type = -99;
        Layer = -99;
        Wheel = -99;
        bw_fw = -99;
        Ext_Int = -99;
        MonoStereo = -99;
        MagField = -99;
        SignCorrection = -99;
        XGlobal = -99;
        YGlobal = -99;
        ZGlobal = -99;
        barycenter = -99;
        hit_std_dev = -99;
        sumx = 0;
        id_detector = -1;
        thick_detector = -1;
        pitch_detector = -1;
        HitNr = 1;

        SiStripRecHit2D lhit;
        TrajectoryStateOnSurface tsos = itm->updatedState();
        const TransientTrackingRecHit::ConstRecHitPointer thit = itm->recHit();
        if ((thit->geographicalId().subdetId() == int(StripSubdetector::TIB)) ||
            thit->geographicalId().subdetId() == int(StripSubdetector::TOB)) {  //include only barrel
          const SiStripMatchedRecHit2D* matchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>((*thit).hit());
          const ProjectedSiStripRecHit2D* phit = dynamic_cast<const ProjectedSiStripRecHit2D*>((*thit).hit());
          const SiStripRecHit2D* hit = dynamic_cast<const SiStripRecHit2D*>((*thit).hit());
          if (phit) {
            lhit = phit->originalHit();
            hit = &lhit;
          }

          LocalVector trackdirection = tsos.localDirection();

          if (matchedhit) {  //if matched hit...

            GluedGeomDet* gdet = (GluedGeomDet*)tkGeom_->idToDet(matchedhit->geographicalId());

            GlobalVector gtrkdir = gdet->toGlobal(trackdirection);

            // THIS THE POINTER TO THE MONO HIT OF A MATCHED HIT

            // top be migrated to the more direct interface of matchedhit
            cache.push_back(matchedhit->monoHit());
            const SiStripRecHit2D* monohit = &cache.back();
            const SiStripRecHit2D::ClusterRef& monocluster = monohit->cluster();
            const GeomDetUnit* monodet = gdet->monoDet();
            // this does not exists anymore! either project the matched or use CPE
            const LocalPoint monoposition = monohit->localPosition();

            StripSubdetector detid = (StripSubdetector)monohit->geographicalId();
            id_detector = detid.rawId();
            thick_detector = monodet->specificSurface().bounds().thickness();
            const StripTopology& mtopol = (StripTopology&)monodet->topology();
            pitch_detector = mtopol.localPitch(monoposition);
            const GlobalPoint monogposition = (monodet->surface()).toGlobal(monoposition);
            ClSize = (monocluster->amplitudes()).size();

            const auto& amplitudes = monocluster->amplitudes();

            barycenter = monocluster->barycenter() - 0.5;
            uint16_t FirstStrip = monocluster->firstStrip();
            auto begin = amplitudes.begin();
            nstrip = 0;
            for (auto idigi = begin; idigi != amplitudes.end(); idigi++) {
              Amplitudes[nstrip] = *idigi;
              sumx += pow(((FirstStrip + idigi - begin) - barycenter), 2) * (*idigi);
              HitCharge += *idigi;
            }
            hit_std_dev = sqrt(sumx / HitCharge);

            XGlobal = monogposition.x();
            YGlobal = monogposition.y();
            ZGlobal = monogposition.z();

            Type = detid.subdetId();
            MonoStereo = detid.stereo();

            if (detid.subdetId() == int(StripSubdetector::TIB)) {
              Layer = tTopo->tibLayer(detid);
              bw_fw = tTopo->tibStringInfo(detid)[0];
              Ext_Int = tTopo->tibStringInfo(detid)[1];
            }
            if (detid.subdetId() == int(StripSubdetector::TOB)) {
              Layer = tTopo->tobLayer(detid);
              bw_fw = tTopo->tobRodInfo(detid)[0];
            }
            if (detid.subdetId() == int(StripSubdetector::TID)) {
              Wheel = tTopo->tidWheel(detid);
              bw_fw = tTopo->tidModuleInfo(detid)[0];
            }
            if (detid.subdetId() == int(StripSubdetector::TEC)) {
              Wheel = tTopo->tecWheel(detid);
              bw_fw = tTopo->tecPetalInfo(detid)[0];
            }

            LocalVector monotkdir = monodet->toLocal(gtrkdir);

            if (monotkdir.z() != 0) {
              // THE LOCAL ANGLE (MONO)
              float tanangle = monotkdir.x() / monotkdir.z();
              TanTrackAngleParallel = monotkdir.y() / monotkdir.z();
              TanTrackAngle = tanangle;
              detparmap::iterator TheDet = detmap.find(detid.rawId());
              LocalVector localmagdir;
              if (TheDet != detmap.end())
                localmagdir = TheDet->second->magfield;
              MagField = localmagdir.mag();
              if (MagField != 0.) {
                LocalVector monoylocal(0, 1, 0);
                float signcorrection = (localmagdir * monoylocal) / (MagField);
                if (signcorrection != 0)
                  SignCorrection = 1 / signcorrection;
              }

              std::map<const SiStripRecHit2D*, std::pair<float, float>, DetIdLess>::iterator alreadystored =
                  hitangleassociation.find(monohit);

              if (alreadystored != hitangleassociation.end()) {  //decide which hit take
                if (itm->estimate() > alreadystored->second.first) {
                  worse_double_hit++;
                }
                if (itm->estimate() < alreadystored->second.first) {
                  better_double_hit++;
                  hitangleassociation.insert(std::make_pair(monohit, std::make_pair(itm->estimate(), tanangle)));
                }
              } else {
                hitangleassociation.insert(make_pair(monohit, std::make_pair(itm->estimate(), tanangle)));
                HitsTree->Fill();
                hitcounter++;
              }

              // THIS THE POINTER TO THE STEREO HIT OF A MATCHED HIT

              // top be migrated to the more direct interface of matchedhit
              cache.push_back(matchedhit->stereoHit());
              const SiStripRecHit2D* stereohit = &cache.back();
              const SiStripRecHit2D::ClusterRef& stereocluster = stereohit->cluster();
              const GeomDetUnit* stereodet = gdet->stereoDet();
              // this does not exists anymore! either project the matched or use CPE
              const LocalPoint stereoposition = stereohit->localPosition();
              StripSubdetector detid = (StripSubdetector)stereohit->geographicalId();
              id_detector = detid.rawId();
              thick_detector = stereodet->specificSurface().bounds().thickness();
              const StripTopology& stopol = (StripTopology&)stereodet->topology();
              pitch_detector = stopol.localPitch(stereoposition);
              const GlobalPoint stereogposition = (stereodet->surface()).toGlobal(stereoposition);

              ClSize = (stereocluster->amplitudes()).size();

              const auto& amplitudes = stereocluster->amplitudes();

              barycenter = stereocluster->barycenter() - 0.5;
              uint16_t FirstStrip = stereocluster->firstStrip();
              auto begin = amplitudes.begin();
              nstrip = 0;
              for (auto idigi = begin; idigi != amplitudes.end(); idigi++) {
                Amplitudes[nstrip] = *idigi;
                sumx += pow(((FirstStrip + idigi - begin) - barycenter), 2) * (*idigi);
                HitCharge += *idigi;
              }
              hit_std_dev = sqrt(sumx / HitCharge);

              XGlobal = stereogposition.x();
              YGlobal = stereogposition.y();
              ZGlobal = stereogposition.z();

              Type = detid.subdetId();
              MonoStereo = detid.stereo();

              if (detid.subdetId() == int(StripSubdetector::TIB)) {
                Layer = tTopo->tibLayer(detid);
                bw_fw = tTopo->tibStringInfo(detid)[0];
                Ext_Int = tTopo->tibStringInfo(detid)[1];
              }
              if (detid.subdetId() == int(StripSubdetector::TOB)) {
                Layer = tTopo->tobLayer(detid);
                bw_fw = tTopo->tobRodInfo(detid)[0];
              }
              if (detid.subdetId() == int(StripSubdetector::TID)) {
                Wheel = tTopo->tidWheel(detid);
                bw_fw = tTopo->tidModuleInfo(detid)[0];
              }
              if (detid.subdetId() == int(StripSubdetector::TEC)) {
                Wheel = tTopo->tecWheel(detid);
                bw_fw = tTopo->tecPetalInfo(detid)[0];
              }

              LocalVector stereotkdir = stereodet->toLocal(gtrkdir);

              if (stereotkdir.z() != 0) {
                // THE LOCAL ANGLE (STEREO)
                float tanangle = stereotkdir.x() / stereotkdir.z();
                TanTrackAngleParallel = stereotkdir.y() / stereotkdir.z();
                TanTrackAngle = tanangle;
                detparmap::iterator TheDet = detmap.find(detid.rawId());
                LocalVector localmagdir;
                if (TheDet != detmap.end())
                  localmagdir = TheDet->second->magfield;
                MagField = localmagdir.mag();
                LocalVector stereoylocal(0, 1, 0);
                if (MagField != 0.) {
                  float signcorrection = (localmagdir * stereoylocal) / (MagField);
                  if (signcorrection != 0)
                    SignCorrection = 1 / signcorrection;
                }

                std::map<const SiStripRecHit2D*, std::pair<float, float>, DetIdLess>::iterator alreadystored =
                    hitangleassociation.find(stereohit);

                if (alreadystored != hitangleassociation.end()) {  //decide which hit take
                  if (itm->estimate() > alreadystored->second.first) {
                    worse_double_hit++;
                  }
                  if (itm->estimate() < alreadystored->second.first) {
                    better_double_hit++;
                    hitangleassociation.insert(std::make_pair(stereohit, std::make_pair(itm->estimate(), tanangle)));
                  }
                } else {
                  hitangleassociation.insert(std::make_pair(stereohit, std::make_pair(itm->estimate(), tanangle)));
                  HitsTree->Fill();
                  hitcounter++;
                }
              }
            }
          } else if (hit) {
            //  hit= POINTER TO THE RECHIT

            const SiStripRecHit2D::ClusterRef& cluster = hit->cluster();

            GeomDetUnit* gdet = (GeomDetUnit*)tkGeom_->idToDet(hit->geographicalId());
            const LocalPoint position = hit->localPosition();
            StripSubdetector detid = (StripSubdetector)hit->geographicalId();
            id_detector = detid.rawId();
            thick_detector = gdet->specificSurface().bounds().thickness();
            const StripTopology& topol = (StripTopology&)gdet->topology();
            pitch_detector = topol.localPitch(position);
            const GlobalPoint gposition = (gdet->surface()).toGlobal(position);

            ClSize = (cluster->amplitudes()).size();

            const auto& amplitudes = cluster->amplitudes();

            barycenter = cluster->barycenter() - 0.5;
            uint16_t FirstStrip = cluster->firstStrip();
            nstrip = 0;
            auto begin = amplitudes.begin();
            for (auto idigi = amplitudes.begin(); idigi != amplitudes.end(); idigi++) {
              Amplitudes[nstrip] = *idigi;
              sumx += pow(((FirstStrip + idigi - begin) - barycenter), 2) * (*idigi);
              HitCharge += *idigi;
            }
            hit_std_dev = sqrt(sumx / HitCharge);

            XGlobal = gposition.x();
            YGlobal = gposition.y();
            ZGlobal = gposition.z();

            Type = detid.subdetId();
            MonoStereo = detid.stereo();

            if (detid.subdetId() == int(StripSubdetector::TIB)) {
              Layer = tTopo->tibLayer(detid);
              bw_fw = tTopo->tibStringInfo(detid)[0];
              Ext_Int = tTopo->tibStringInfo(detid)[1];
            }
            if (detid.subdetId() == int(StripSubdetector::TOB)) {
              Layer = tTopo->tobLayer(detid);
              bw_fw = tTopo->tobRodInfo(detid)[0];
            }
            if (detid.subdetId() == int(StripSubdetector::TID)) {
              Wheel = tTopo->tidWheel(detid);
              bw_fw = tTopo->tidModuleInfo(detid)[0];
            }
            if (detid.subdetId() == int(StripSubdetector::TEC)) {
              Wheel = tTopo->tecWheel(detid);
              bw_fw = tTopo->tecPetalInfo(detid)[0];
            }

            if (trackdirection.z() != 0) {
              // THE LOCAL ANGLE
              float tanangle = trackdirection.x() / trackdirection.z();
              TanTrackAngleParallel = trackdirection.y() / trackdirection.z();
              TanTrackAngle = tanangle;
              detparmap::iterator TheDet = detmap.find(detid.rawId());
              LocalVector localmagdir;
              if (TheDet != detmap.end())
                localmagdir = TheDet->second->magfield;
              MagField = localmagdir.mag();
              if (MagField != 0.) {
                LocalVector ylocal(0, 1, 0);
                float signcorrection = (localmagdir * ylocal) / (MagField);
                if (signcorrection != 0)
                  SignCorrection = 1 / signcorrection;
              }

              std::map<const SiStripRecHit2D*, std::pair<float, float>, DetIdLess>::iterator alreadystored =
                  hitangleassociation.find(hit);

              if (alreadystored != hitangleassociation.end()) {  //decide which hit take
                if (itm->estimate() > alreadystored->second.first) {
                  worse_double_hit++;
                }
                if (itm->estimate() < alreadystored->second.first) {
                  better_double_hit++;
                  hitangleassociation.insert(std::make_pair(hit, std::make_pair(itm->estimate(), tanangle)));
                }
              } else {
                hitangleassociation.insert(std::make_pair(hit, std::make_pair(itm->estimate(), tanangle)));
                HitsTree->Fill();
                hitcounter++;
              }
            }
          }
        }
      }
    }
  }
  std::map<const SiStripRecHit2D*, std::pair<float, float>, DetIdLess>::iterator hitsiter;

  for (hitsiter = hitangleassociation.begin(); hitsiter != hitangleassociation.end(); hitsiter++) {
    hitcounter_2ndloop++;

    const SiStripRecHit2D* hit = hitsiter->first;
    const SiStripRecHit2D::ClusterRef& cluster = hit->cluster();

    size = (cluster->amplitudes()).size();

    StripSubdetector detid = (StripSubdetector)hit->geographicalId();

    float tangent = hitsiter->second.second;

    //Sign and XZ plane projection correction applied in TrackLocalAngle (TIB|TOB layers)

    detparmap::iterator thedet = detmap.find(detid.rawId());
    LocalVector localmagdir;
    if (thedet != detmap.end())
      localmagdir = thedet->second->magfield;
    float localmagfield = localmagdir.mag();

    if (localmagfield != 0.) {
      LocalVector ylocal(0, 1, 0);

      float normprojection = (localmagdir * ylocal) / (localmagfield);

      if (normprojection == 0.)
        LogDebug("SiStripLAProfileBooker::analyze") << "Error: YBprojection = 0";

      else {
        float signprojcorrection = 1 / normprojection;
        tangent *= signprojcorrection;
      }
    }

    //Filling histograms

    histomap::iterator thehisto = histos.find(detid.rawId());

    if (thehisto == histos.end())
      edm::LogError("SiStripLAProfileBooker::analyze")
          << "Error: the profile associated to" << detid.rawId() << "does not exist! ";
    else
      thehisto->second->Fill(tangent, size);

    //Summary histograms
    std::string name;
    unsigned int layerid;
    getlayer(detid, tTopo, name, layerid);
    histomap::iterator thesummaryhisto = summaryhisto.find(layerid);
    if (thesummaryhisto == summaryhisto.end())
      edm::LogError("SiStripLAProfileBooker::analyze")
          << "Error: the profile associated to subdet " << name << "does not exist! ";
    else
      thesummaryhisto->second->Fill(tangent, size);
  }
}

beginRun()

void SiStripLAProfileBooker::beginRun ( edm::Run const &  , const edm::EventSetup &  c  )

override

Definition at line 66 of file SiStripLAProfileBooker.cc.

References Amplitudes, better_double_hit, dqm::implementation::IBooker::bookProfile(), bw_fw, DummyCfis::c, chi2norm, ClSize, conf_, SiStripHistoId::createHistoId(), dbe_, detCablingToken_, TrackerGeometry::detIds(), detmap, EtaTrack, Event_Tree, eventcounter, EventCounter, EventNumber, EventTree, Ext_Int, getlayer(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hFile, histos, hit_std_dev, Hit_Tree, HitCharge, hitcounter, hitcounter_2ndloop, HitNr, HitPerTrack, HitsTree, id_detector, TrackerGeometry::idToDet(), createfilelist::int, Layer, StripTopology::localPitch(), SiStripLAProfileBooker::detparameters::magfield, MagField, magFieldToken_, Momentum, MonoStereo, Skims_PA_cff::name, Utilities::operator, AlCaHLTBitMon_ParallelJobs::p, ParticleCharge, PhiTrack, SiStripLAProfileBooker::detparameters::pitch, pitch_detector, PostProcessor_cff::profile, pt, DetId::rawId(), RunNumber, dqm::impl::MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), SiStripFolderOrganizer::setSiStripFolder(), SignCorrection, AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), summarydetmap, summaryhisto, TanTrackAngle, TanTrackAngleParallel, StripSubdetector::TEC, thick_detector, SiStripLAProfileBooker::detparameters::thickness, Calorimetry_cff::thickness, StripSubdetector::TIB, StripSubdetector::TID, tkGeom_, tkGeomToken_, Zmumu_cfg::TkTag, StripSubdetector::TOB, toLocal(), Track_Tree, trackcollsize, TrackCounter, TrackTree, trajsize, tTopoToken_, Type, Wheel, worse_double_hit, XGlobal, YGlobal, and ZGlobal.

                                                                           {
  const TrackerTopology* const tTopo = &c.getData(tTopoToken_);
  tkGeom_ = &c.getData(tkGeomToken_);
  const auto& magField = c.getData(magFieldToken_);

  std::vector<uint32_t> activeDets;
  if (conf_.getParameter<bool>("UseStripCablingDB")) {
    activeDets.clear();
    c.getData(detCablingToken_).addActiveDetectorsRawIds(activeDets);
  } else {
    const TrackerGeometry::DetIdContainer& Id = tkGeom_->detIds();
    TrackerGeometry::DetIdContainer::const_iterator Iditer;
    activeDets.clear();
    for (Iditer = Id.begin(); Iditer != Id.end(); Iditer++) {
      activeDets.push_back(Iditer->rawId());
    }
  }

  edm::InputTag TkTag = conf_.getParameter<edm::InputTag>("Tracks");
  //Get Ids;
  double ModuleRangeMin = conf_.getParameter<double>("ModuleXMin");
  double ModuleRangeMax = conf_.getParameter<double>("ModuleXMax");
  double TIBRangeMin = conf_.getParameter<double>("TIBXMin");
  double TIBRangeMax = conf_.getParameter<double>("TIBXMax");
  double TOBRangeMin = conf_.getParameter<double>("TOBXMin");
  double TOBRangeMax = conf_.getParameter<double>("TOBXMax");
  int TIB_bin = conf_.getParameter<int>("TIB_bin");
  int TOB_bin = conf_.getParameter<int>("TOB_bin");
  int SUM_bin = conf_.getParameter<int>("SUM_bin");

  hFile = new TFile(conf_.getUntrackedParameter<std::string>("treeName").c_str(), "RECREATE");

  Hit_Tree = hFile->mkdir("Hit_Tree");
  Track_Tree = hFile->mkdir("Track_Tree");
  Event_Tree = hFile->mkdir("Event_Tree");

  HitsTree = new TTree("HitsTree", "HitsTree");

  HitsTree->Branch("RunNumber", &RunNumber, "RunNumber/I");
  HitsTree->Branch("EventNumber", &EventNumber, "EventNumber/I");
  HitsTree->Branch("TanTrackAngle", &TanTrackAngle, "TanTrackAngle/F");
  HitsTree->Branch("TanTrackAngleParallel", &TanTrackAngleParallel, "TanTrackAngleParallel/F");
  HitsTree->Branch("ClSize", &ClSize, "ClSize/I");
  HitsTree->Branch("HitCharge", &HitCharge, "HitCharge/I");
  HitsTree->Branch("Hit_Std_Dev", &hit_std_dev, "hit_std_dev/F");
  HitsTree->Branch("Type", &Type, "Type/I");
  HitsTree->Branch("Layer", &Layer, "Layer/I");
  HitsTree->Branch("Wheel", &Wheel, "Wheel/I");
  HitsTree->Branch("bw_fw", &bw_fw, "bw_fw/I");
  HitsTree->Branch("Ext_Int", &Ext_Int, "Ext_Int/I");
  HitsTree->Branch("MonoStereo", &MonoStereo, "MonoStereo/I");
  HitsTree->Branch("MagField", &MagField, "MagField/F");
  HitsTree->Branch("SignCorrection", &SignCorrection, "SignCorrection/F");
  HitsTree->Branch("XGlobal", &XGlobal, "XGlobal/F");
  HitsTree->Branch("YGlobal", &YGlobal, "YGlobal/F");
  HitsTree->Branch("ZGlobal", &ZGlobal, "ZGlobal/F");
  HitsTree->Branch("ParticleCharge", &ParticleCharge, "ParticleCharge/I");
  HitsTree->Branch("Momentum", &Momentum, "Momentum/F");
  HitsTree->Branch("pt", &pt, "pt/F");
  HitsTree->Branch("chi2norm", &chi2norm, "chi2norm/F");
  HitsTree->Branch("EtaTrack", &EtaTrack, "EtaTrack/F");
  HitsTree->Branch("PhiTrack", &PhiTrack, "PhiTrack/F");
  HitsTree->Branch("TrajSize", &trajsize, "trajsize/I");
  HitsTree->Branch("HitNr", &HitNr, "HitNr/I");
  HitsTree->Branch("HitPerTrack", &HitPerTrack, "HitPerTrack/I");
  HitsTree->Branch("id_detector", &id_detector, "id_detector/I");
  HitsTree->Branch("thick_detector", &thick_detector, "thick_detector/F");
  HitsTree->Branch("pitch_detector", &pitch_detector, "pitch_detector/F");
  HitsTree->Branch("Amplitudes", Amplitudes, "Amplitudes[ClSize]/I");

  HitsTree->SetDirectory(Hit_Tree);

  TrackTree = new TTree("TrackTree", "TrackTree");

  TrackTree->Branch("TrackCounter", &TrackCounter, "TrackCounter/I");

  TrackTree->SetDirectory(Track_Tree);

  EventTree = new TTree("EventTree", "EventTree");

  EventTree->Branch("EventCounter", &EventCounter, "EventCounter/I");

  EventTree->SetDirectory(Event_Tree);

  // use SistripHistoId for producing histogram id (and title)
  SiStripHistoId hidmanager;

  // create SiStripFolderOrganizer
  SiStripFolderOrganizer folder_organizer;

  dbe_ = edm::Service<DQMStore>().operator->();

  //get all detids

  for (std::vector<uint32_t>::const_iterator Id = activeDets.begin(); Id != activeDets.end(); Id++) {
    //  for(Iditer=Id.begin();Iditer!=Id.end();Iditer++){ //loop on detids
    DetId Iditero = DetId(*Id);
    DetId* Iditer = &Iditero;
    if ((Iditer->subdetId() == int(StripSubdetector::TIB)) ||
        (Iditer->subdetId() == int(StripSubdetector::TOB))) {  //include only barrel

      int module_bin = 0;
      if (Iditer->subdetId() == int(StripSubdetector::TIB)) {
        module_bin = TIB_bin;
      } else {
        module_bin = TOB_bin;
      }

      // create a TProfile for each module
      StripSubdetector subid(*Iditer);
      std::string hid;
      //Mono single sided detectors
      LocalPoint p;
      auto stripdet = tkGeom_->idToDet(subid);
      if (!stripdet->isLeaf())
        continue;
      const StripTopology& topol = (const StripTopology&)stripdet->topology();
      float thickness = stripdet->specificSurface().bounds().thickness();

      folder_organizer.setDetectorFolder(Iditer->rawId(), tTopo);
      hid = hidmanager.createHistoId(TkTag.label(), "det", Iditer->rawId());
      MonitorElement* profile = dbe_->bookProfile(hid, hid, module_bin, ModuleRangeMin, ModuleRangeMax, 20, 0, 5, "");
      detparameters* param = new detparameters;
      histos[Iditer->rawId()] = profile;
      detmap[Iditer->rawId()] = param;
      param->thickness = thickness * 10000;
      param->pitch = topol.localPitch(p) * 10000;

      const GlobalPoint globalp = (stripdet->surface()).toGlobal(p);
      GlobalVector globalmagdir = magField.inTesla(globalp);
      param->magfield = (stripdet->surface()).toLocal(globalmagdir);

      profile->setAxisTitle("tan(#theta_{t})", 1);
      profile->setAxisTitle("Cluster size", 2);

      // create a summary histo if it does not exist already
      std::string name;
      unsigned int layerid;
      getlayer(subid, tTopo, name, layerid);
      name += TkTag.label();
      if (summaryhisto.find(layerid) == (summaryhisto.end())) {
        folder_organizer.setSiStripFolder();
        MonitorElement* summaryprofile = nullptr;
        if (subid.subdetId() == int(StripSubdetector::TIB) || subid.subdetId() == int(StripSubdetector::TID))
          summaryprofile = dbe_->bookProfile(name, name, SUM_bin, TIBRangeMin, TIBRangeMax, 20, 0, 5, "");
        else if (subid.subdetId() == int(StripSubdetector::TOB) || subid.subdetId() == int(StripSubdetector::TEC))
          summaryprofile = dbe_->bookProfile(name, name, SUM_bin, TOBRangeMin, TOBRangeMax, 20, 0, 5, "");
        if (summaryprofile) {
          detparameters* summaryparam = new detparameters;
          summaryhisto[layerid] = summaryprofile;
          summarydetmap[layerid] = summaryparam;
          summaryparam->thickness = thickness * 10000;
          summaryparam->pitch = topol.localPitch(p) * 10000;
          summaryprofile->setAxisTitle("tan(#theta_{t})", 1);
          summaryprofile->setAxisTitle("Cluster size", 2);
        }
      }
    }
  }

  trackcollsize = 0;
  trajsize = 0;
  RunNumber = 0;
  EventNumber = 0;
  hitcounter = 0;
  hitcounter_2ndloop = 0;
  worse_double_hit = 0;
  better_double_hit = 0;
  eventcounter = 0;

  EventCounter = 1;
  TrackCounter = 1;
}

endJob()

void SiStripLAProfileBooker::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 737 of file SiStripLAProfileBooker.cc.

References conf_, dbe_, edm::ParameterSet::getUntrackedParameter(), hFile, dqm::implementation::DQMStore::save(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by o2olib.O2ORunMgr::executeJob().

                                    {
  std::string outputFile_ = conf_.getUntrackedParameter<std::string>("fileName", "LorentzAngle.root");
  dbe_->save(outputFile_);

  hFile->Write();
  hFile->Close();
}

endRun()

void SiStripLAProfileBooker::endRun ( edm::Run const &  , const edm::EventSetup &  c  )

override

Definition at line 240 of file SiStripLAProfileBooker.cc.

{}

getlayer()

void SiStripLAProfileBooker::getlayer	(	const DetId &	detid,
		const TrackerTopology *	tTopo,
		std::string &	name,
		unsigned int &	layerid
	)

Definition at line 705 of file SiStripLAProfileBooker.cc.

References createfilelist::int, nano_mu_digi_cff::layer, Skims_PA_cff::name, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecRing(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

Referenced by analyze(), and beginRun().

                                                             {
  int layer = 0;
  std::stringstream layernum;

  if (detid.subdetId() == int(StripSubdetector::TIB)) {
    name += "TIB_Layer_";
    layer = tTopo->tibLayer(detid);
  }

  else if (detid.subdetId() == int(StripSubdetector::TID)) {
    name += "TID_Ring_";
    layer = tTopo->tidRing(detid);
  }

  else if (detid.subdetId() == int(StripSubdetector::TOB)) {
    name += "TOB_Layer_";
    layer = tTopo->tobLayer(detid);

  }

  else if (detid.subdetId() == int(StripSubdetector::TEC)) {
    name += "TEC_Ring_";
    layer = tTopo->tecRing(detid);
  }
  layernum << layer;
  name += layernum.str();
  layerid = detid.subdetId() * 10 + layer;
}

Member Data Documentation

Amplitudes

uint8_t SiStripLAProfileBooker::Amplitudes[100]

private

Definition at line 68 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

barycenter

float SiStripLAProfileBooker::barycenter

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

better_double_hit

int SiStripLAProfileBooker::better_double_hit

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

bw_fw

int SiStripLAProfileBooker::bw_fw

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

chi2norm

float SiStripLAProfileBooker::chi2norm

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

ClSize

int SiStripLAProfileBooker::ClSize

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

conf_

edm::ParameterSet SiStripLAProfileBooker::conf_

private

Definition at line 82 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and endJob().

dbe_

DQMStore* SiStripLAProfileBooker::dbe_

private

Definition at line 78 of file SiStripLAProfileBooker.h.

Referenced by beginRun(), and endJob().

detCablingToken_

edm::ESGetToken<SiStripDetCabling, SiStripDetCablingRcd> SiStripLAProfileBooker::detCablingToken_

private

Definition at line 89 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

detmap

detparmap SiStripLAProfileBooker::detmap

private

Definition at line 80 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and ~SiStripLAProfileBooker().

EtaTrack

float SiStripLAProfileBooker::EtaTrack

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

Event_Tree

TDirectory * SiStripLAProfileBooker::Event_Tree

private

Definition at line 73 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

eventcounter

int SiStripLAProfileBooker::eventcounter

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

EventCounter

int SiStripLAProfileBooker::EventCounter

private

Definition at line 66 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

EventNumber

int SiStripLAProfileBooker::EventNumber

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

EventTree

TTree * SiStripLAProfileBooker::EventTree

private

Definition at line 70 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

Ext_Int

int SiStripLAProfileBooker::Ext_Int

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

hFile

TFile* SiStripLAProfileBooker::hFile

private

Definition at line 71 of file SiStripLAProfileBooker.h.

Referenced by beginRun(), endJob(), and ~SiStripLAProfileBooker().

histos

histomap SiStripLAProfileBooker::histos

private

Definition at line 75 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and DisplayManager.DisplayManager::Draw().

hit_std_dev

float SiStripLAProfileBooker::hit_std_dev

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

Hit_Tree

TDirectory* SiStripLAProfileBooker::Hit_Tree

private

Definition at line 73 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

HitCharge

int SiStripLAProfileBooker::HitCharge

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

hitcounter

int SiStripLAProfileBooker::hitcounter

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

hitcounter_2ndloop

int SiStripLAProfileBooker::hitcounter_2ndloop

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

HitNr

int SiStripLAProfileBooker::HitNr

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

HitPerTrack

int SiStripLAProfileBooker::HitPerTrack

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

HitsTree

TTree* SiStripLAProfileBooker::HitsTree

private

Definition at line 70 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

id_detector

int SiStripLAProfileBooker::id_detector

private

Definition at line 66 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

Layer

int SiStripLAProfileBooker::Layer

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

MagField

float SiStripLAProfileBooker::MagField

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

magFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SiStripLAProfileBooker::magFieldToken_

private

Definition at line 88 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

Momentum

float SiStripLAProfileBooker::Momentum

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

MonoStereo

int SiStripLAProfileBooker::MonoStereo

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

nstrip

int SiStripLAProfileBooker::nstrip

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

ParticleCharge

int SiStripLAProfileBooker::ParticleCharge

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

PhiTrack

float SiStripLAProfileBooker::PhiTrack

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

pitch_detector

float SiStripLAProfileBooker::pitch_detector

private

Definition at line 67 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

pt

float SiStripLAProfileBooker::pt

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by Particle.Particle::__str__(), ZMuMuRochCorAnalyzer.DiMuon::__str__(), DiObject.DiMuon::__str__(), analyze(), beginRun(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), datamodel.Object::p4(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

RunNumber

int SiStripLAProfileBooker::RunNumber

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and Page1Parser.Page1Parser::ParsePage1().

SignCorrection

float SiStripLAProfileBooker::SignCorrection

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

size

int SiStripLAProfileBooker::size

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), analyze(), esMonitoring.LineHistoryEnd::pop(), esMonitoring.LineHistoryEnd::push(), cuy.FindIssue::startElement(), esMonitoring.LineHistoryEnd::write(), and esMonitoring.LineHistoryStart::write().

summarydetmap

detparmap SiStripLAProfileBooker::summarydetmap

private

Definition at line 81 of file SiStripLAProfileBooker.h.

Referenced by beginRun(), and ~SiStripLAProfileBooker().

summaryhisto

histomap SiStripLAProfileBooker::summaryhisto

private

Definition at line 76 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

sumx

float SiStripLAProfileBooker::sumx

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

TanTrackAngle

float SiStripLAProfileBooker::TanTrackAngle

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TanTrackAngleParallel

float SiStripLAProfileBooker::TanTrackAngleParallel

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

thick_detector

float SiStripLAProfileBooker::thick_detector

private

Definition at line 67 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

tkGeom_

const TrackerGeometry* SiStripLAProfileBooker::tkGeom_ = nullptr

private

Definition at line 85 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

tkGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiStripLAProfileBooker::tkGeomToken_

private

Definition at line 87 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

Track_Tree

TDirectory * SiStripLAProfileBooker::Track_Tree

private

Definition at line 73 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

trackcollsize

int SiStripLAProfileBooker::trackcollsize

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TrackCounter

int SiStripLAProfileBooker::TrackCounter

private

Definition at line 66 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

TrackTree

TTree * SiStripLAProfileBooker::TrackTree

private

Definition at line 70 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

trajsize

int SiStripLAProfileBooker::trajsize

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

treename_

std::string SiStripLAProfileBooker::treename_

private

Definition at line 83 of file SiStripLAProfileBooker.h.

tTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiStripLAProfileBooker::tTopoToken_

private

Definition at line 86 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

Type

int SiStripLAProfileBooker::Type

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and BeamSpotObj.BeamSpot::Reset().

Wheel

int SiStripLAProfileBooker::Wheel

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

worse_double_hit

int SiStripLAProfileBooker::worse_double_hit

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

XGlobal

float SiStripLAProfileBooker::XGlobal

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

YGlobal

float SiStripLAProfileBooker::YGlobal

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

ZGlobal

float SiStripLAProfileBooker::ZGlobal

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().