#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 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< 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 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 ()

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

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

private

Definition at line 57 of file SiStripLAProfileBooker.h.

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

Definition at line 36 of file SiStripLAProfileBooker.h.

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

private

Definition at line 58 of file SiStripLAProfileBooker.h.

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

Definition at line 35 of file SiStripLAProfileBooker.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 55 of file SiStripLAProfileBooker.cc.

References detCablingToken_.

     : 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>()}

SiStripLAProfileBooker::~SiStripLAProfileBooker ( )

override

Definition at line 242 of file SiStripLAProfileBooker.cc.

References detmap, 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

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, SplitLinear::begin, better_double_hit, Surface::bounds(), bw_fw, utilities::cache(), chi2norm, ClSize, SiStripRecHit2D::cluster(), conf_, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, detmap, EtaTrack, edm::EventID::event(), eventcounter, EventNumber, EventTree, Ext_Int, TrackingRecHit::geographicalId(), edm::Event::getByLabel(), edm::EventSetup::getData(), getlayer(), edm::ParameterSet::getParameter(), histos, hit_std_dev, HitCharge, hitcounter, hitcounter_2ndloop, HitNr, HitPerTrack, HitsTree, mps_fire::i, edm::EventBase::id(), id_detector, TrackerGeometry::idToDet(), Layer, TrajectoryStateOnSurface::localDirection(), StripTopology::localPitch(), BaseTrackerRecHit::localPosition(), LogDebug, PV3DBase< T, PVType, FrameType >::mag(), MagField, Momentum, GluedGeomDet::monoDet(), SiStripMatchedRecHit2D::monoHit(), MonoStereo, mergeVDriftHistosByStation::name, nstrip, ProjectedSiStripRecHit2D::originalHit(), ParticleCharge, PhiTrack, pitch_detector, position, funct::pow(), pt, DetId::rawId(), edm::EventID::run(), 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_, StripSubdetector::TOB, TrackerTopology::tobLayer(), TrackerTopology::tobRodInfo(), GeomDet::toGlobal(), GeomDet::toLocal(), GeomDet::topology(), HLT_FULL_cff::trackCollection, trackcollsize, tracks, TrackTree, trajsize, tTopoToken_, Type, 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);
   }
 }

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

override

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();
 }

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

override

Definition at line 240 of file SiStripLAProfileBooker.cc.

240 {}

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

Definition at line 705 of file SiStripLAProfileBooker.cc.

References phase1PixelTopology::layer, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecRing(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

Referenced by analyze(), and for().

                                                              {
   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

uint8_t SiStripLAProfileBooker::Amplitudes[100]

private

Definition at line 68 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::barycenter

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::better_double_hit

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::bw_fw

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::chi2norm

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::ClSize

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

edm::ParameterSet SiStripLAProfileBooker::conf_

private

Definition at line 82 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and endJob().

DQMStore* SiStripLAProfileBooker::dbe_

private

Definition at line 78 of file SiStripLAProfileBooker.h.

Referenced by endJob(), and for().

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

private

Definition at line 89 of file SiStripLAProfileBooker.h.

Referenced by if(), and SiStripLAProfileBooker().

detparmap SiStripLAProfileBooker::detmap

private

Definition at line 80 of file SiStripLAProfileBooker.h.

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

float SiStripLAProfileBooker::EtaTrack

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

TDirectory * SiStripLAProfileBooker::Event_Tree

private

Definition at line 73 of file SiStripLAProfileBooker.h.

int SiStripLAProfileBooker::eventcounter

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::EventCounter

private

Definition at line 66 of file SiStripLAProfileBooker.h.

int SiStripLAProfileBooker::EventNumber

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

TTree * SiStripLAProfileBooker::EventTree

private

Definition at line 70 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::Ext_Int

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

TFile* SiStripLAProfileBooker::hFile

private

Definition at line 71 of file SiStripLAProfileBooker.h.

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

histomap SiStripLAProfileBooker::histos

private

Definition at line 75 of file SiStripLAProfileBooker.h.

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

float SiStripLAProfileBooker::hit_std_dev

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

TDirectory* SiStripLAProfileBooker::Hit_Tree

private

Definition at line 73 of file SiStripLAProfileBooker.h.

int SiStripLAProfileBooker::HitCharge

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::hitcounter

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::hitcounter_2ndloop

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::HitNr

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::HitPerTrack

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

TTree* SiStripLAProfileBooker::HitsTree

private

Definition at line 70 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::id_detector

private

Definition at line 66 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::Layer

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::MagField

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

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

private

Definition at line 88 of file SiStripLAProfileBooker.h.

float SiStripLAProfileBooker::Momentum

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::MonoStereo

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::nstrip

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::ParticleCharge

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::PhiTrack

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::pitch_detector

private

Definition at line 67 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::pt

private

Definition at line 62 of file SiStripLAProfileBooker.h.

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

int SiStripLAProfileBooker::RunNumber

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::SignCorrection

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

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().

detparmap SiStripLAProfileBooker::summarydetmap

private

Definition at line 81 of file SiStripLAProfileBooker.h.

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

histomap SiStripLAProfileBooker::summaryhisto

private

Definition at line 76 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and for().

float SiStripLAProfileBooker::sumx

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::TanTrackAngle

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::TanTrackAngleParallel

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::thick_detector

private

Definition at line 67 of file SiStripLAProfileBooker.h.

Referenced by analyze().

const TrackerGeometry* SiStripLAProfileBooker::tkGeom_ = nullptr

private

Definition at line 85 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and for().

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

private

Definition at line 87 of file SiStripLAProfileBooker.h.

TDirectory * SiStripLAProfileBooker::Track_Tree

private

Definition at line 73 of file SiStripLAProfileBooker.h.

int SiStripLAProfileBooker::trackcollsize

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::TrackCounter

private

Definition at line 66 of file SiStripLAProfileBooker.h.

TTree * SiStripLAProfileBooker::TrackTree

private

Definition at line 70 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::trajsize

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

std::string SiStripLAProfileBooker::treename_

private

Definition at line 83 of file SiStripLAProfileBooker.h.

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

private

Definition at line 86 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::Type

private

Definition at line 60 of file SiStripLAProfileBooker.h.

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

int SiStripLAProfileBooker::Wheel

private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::worse_double_hit

private

Definition at line 64 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::XGlobal

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::YGlobal

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::ZGlobal

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by analyze().

Classes

Public Types

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation