#include <SiStripLAProfileBooker.h>

Inheritance diagram for SiStripLAProfileBooker:

Classes
struct	detparameters

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	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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

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_

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

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

TDirectory *	Track_Tree

int	trackcollsize

int	TrackCounter

const TrackerGeometry *	tracker

TTree *	TrackTree

int	trajsize

std::string	treename_

int	Type

int	Wheel

int	worse_double_hit

float	XGlobal

float	YGlobal

float	ZGlobal

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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

Detailed Description

Definition at line 30 of file SiStripLAProfileBooker.h.

Member Typedef Documentation

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

private

Definition at line 49 of file SiStripLAProfileBooker.h.

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

private

Definition at line 50 of file SiStripLAProfileBooker.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 of file SiStripLAProfileBooker.cc.

                                                                           : 
   conf_(conf)
 {
 }

SiStripLAProfileBooker::~SiStripLAProfileBooker ( )

override

Definition at line 261 of file SiStripLAProfileBooker.cc.

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

override

Definition at line 270 of file SiStripLAProfileBooker.cc.

References Amplitudes, barycenter, begin, edm::AssociationMap< Tag >::begin(), better_double_hit, Surface::bounds(), bw_fw, helperFunctions::cache(), chi2norm, ClSize, SiStripRecHit2D::cluster(), conf_, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, detmap, edm::AssociationMap< Tag >::end(), EtaTrack, edm::EventID::event(), eventcounter, EventNumber, EventTree, Ext_Int, TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), 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(), createfilelist::int, Layer, TrajectoryStateOnSurface::localDirection(), StripTopology::localPitch(), BaseTrackerRecHit::localPosition(), LogDebug, PV3DBase< T, PVType, FrameType >::mag(), MagField, Momentum, GluedGeomDet::monoDet(), SiStripMatchedRecHit2D::monoHit(), MonoStereo, dataset::name, nstrip, ProjectedSiStripRecHit2D::originalHit(), ParticleCharge, PhiTrack, pitch_detector, position, funct::pow(), edm::ESHandle< T >::product(), 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(), StripSubdetector::TOB, TrackerTopology::tobLayer(), TrackerTopology::tobRodInfo(), GeomDet::toGlobal(), GeomDet::toLocal(), GeomDet::topology(), findElectronsInSiStrips_cfi::trackCollection, trackcollsize, tracker, l1t::tracks, TrackTree, trajsize, Type, Wheel, worse_double_hit, PV3DBase< T, PVType, FrameType >::x(), XGlobal, PV3DBase< T, PVType, FrameType >::y(), YGlobal, PV3DBase< T, PVType, FrameType >::z(), and ZGlobal.

 {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
   
   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 *)tracker->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 *)tracker->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

Definition at line 71 of file SiStripLAProfileBooker.cc.

References SiStripDetCabling::addActiveDetectorsRawIds(), Amplitudes, better_double_hit, DQMStore::bookProfile(), bw_fw, chi2norm, ClSize, conf_, SiStripHistoId::createHistoId(), dbe_, TrackerGeometry::detIds(), detmap, EtaTrack, Event_Tree, eventcounter, EventCounter, EventNumber, EventTree, Ext_Int, edm::EventSetup::get(), 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, MagneticField::inTesla(), edm::InputTag::label(), Layer, StripTopology::localPitch(), SiStripLAProfileBooker::detparameters::magfield, MagField, Momentum, MonoStereo, dataset::name, Utilities::operator, AlCaHLTBitMon_ParallelJobs::p, ParticleCharge, PhiTrack, SiStripLAProfileBooker::detparameters::pitch, pitch_detector, edm::ESHandle< T >::product(), PostProcessor_cff::profile, pt, DetId::rawId(), RunNumber, MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), SiStripFolderOrganizer::setSiStripFolder(), SignCorrection, AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), summarydetmap, summaryhisto, TanTrackAngle, TanTrackAngleParallel, StripSubdetector::TEC, thick_detector, SiStripLAProfileBooker::detparameters::thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, toLocal(), Track_Tree, trackcollsize, TrackCounter, tracker, TrackTree, trajsize, Type, Wheel, worse_double_hit, XGlobal, YGlobal, and ZGlobal.

                                                                           {
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   c.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
  
   //get magnetic field and geometry from ES
   edm::ESHandle<MagneticField> esmagfield;
   c.get<IdealMagneticFieldRecord>().get(esmagfield);
   const MagneticField *  magfield=&(*esmagfield);
   
   edm::ESHandle<TrackerGeometry> estracker;
   c.get<TrackerDigiGeometryRecord>().get(estracker);
   tracker=&(*estracker); 
 
   std::vector<uint32_t> activeDets;
   edm::ESHandle<SiStripDetCabling> tkmechstruct=nullptr;
   if (conf_.getParameter<bool>("UseStripCablingDB")){ 
     c.get<SiStripDetCablingRcd>().get(tkmechstruct);
     activeDets.clear();
     tkmechstruct->addActiveDetectorsRawIds(activeDets);
   }
   else {
     const TrackerGeometry::DetIdContainer& Id = estracker->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 = tracker->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;
   
 }

void SiStripLAProfileBooker::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 764 of file SiStripLAProfileBooker.cc.

References conf_, dbe_, edm::ParameterSet::getUntrackedParameter(), hFile, 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::getlayer	(	const DetId &	detid,
		const TrackerTopology *	tTopo,
		std::string &	name,
		unsigned int &	layerid
	)

Definition at line 730 of file SiStripLAProfileBooker.cc.

References createfilelist::int, 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

uint8_t SiStripLAProfileBooker::Amplitudes[100]

private

Definition at line 57 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::barycenter

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::better_double_hit

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::bw_fw

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::chi2norm

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::ClSize

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

edm::ParameterSet SiStripLAProfileBooker::conf_

private

Definition at line 71 of file SiStripLAProfileBooker.h.

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

DQMStore* SiStripLAProfileBooker::dbe_

private

Definition at line 67 of file SiStripLAProfileBooker.h.

Referenced by beginRun(), and endJob().

detparmap SiStripLAProfileBooker::detmap

private

Definition at line 69 of file SiStripLAProfileBooker.h.

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

float SiStripLAProfileBooker::EtaTrack

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TDirectory * SiStripLAProfileBooker::Event_Tree

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::eventcounter

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::EventCounter

private

Definition at line 55 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::EventNumber

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TTree * SiStripLAProfileBooker::EventTree

private

Definition at line 59 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::Ext_Int

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TFile* SiStripLAProfileBooker::hFile

private

Definition at line 60 of file SiStripLAProfileBooker.h.

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

histomap SiStripLAProfileBooker::histos

private

Definition at line 64 of file SiStripLAProfileBooker.h.

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

float SiStripLAProfileBooker::hit_std_dev

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TDirectory* SiStripLAProfileBooker::Hit_Tree

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::HitCharge

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::hitcounter

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::hitcounter_2ndloop

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::HitNr

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::HitPerTrack

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TTree* SiStripLAProfileBooker::HitsTree

private

Definition at line 59 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::id_detector

private

Definition at line 55 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::Layer

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::MagField

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::Momentum

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::MonoStereo

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::nstrip

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::ParticleCharge

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::PhiTrack

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::pitch_detector

private

Definition at line 56 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::pt

private

Definition at line 53 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(), 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 52 of file SiStripLAProfileBooker.h.

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

float SiStripLAProfileBooker::SignCorrection

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::size

private

Definition at line 54 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 70 of file SiStripLAProfileBooker.h.

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

histomap SiStripLAProfileBooker::summaryhisto

private

Definition at line 65 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::sumx

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::TanTrackAngle

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::TanTrackAngleParallel

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::thick_detector

private

Definition at line 56 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TDirectory * SiStripLAProfileBooker::Track_Tree

private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::trackcollsize

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::TrackCounter

private

Definition at line 55 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

const TrackerGeometry* SiStripLAProfileBooker::tracker

private

Definition at line 74 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TTree * SiStripLAProfileBooker::TrackTree

private

Definition at line 59 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::trajsize

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

std::string SiStripLAProfileBooker::treename_

private

Definition at line 72 of file SiStripLAProfileBooker.h.

int SiStripLAProfileBooker::Type

private

Definition at line 52 of file SiStripLAProfileBooker.h.

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

int SiStripLAProfileBooker::Wheel

private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::worse_double_hit

private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::XGlobal

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::YGlobal

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::ZGlobal

private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

Classes

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation