SiStripLAProfileBooker Class Reference

#include <SiStripLAProfileBooker.h>

Inheritance diagram for SiStripLAProfileBooker:

Classes
struct	detparameters

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
void	beginRun (edm::Run const &, const edm::EventSetup &c)
void	endJob ()
void	getlayer (const DetId &detid, const TrackerTopology *tTopo, std::string &name, unsigned int &layerid)
	SiStripLAProfileBooker (const edm::ParameterSet &conf)
	~SiStripLAProfileBooker ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
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
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 &)
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

(

)

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
	)

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
	)

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(), NewPostProcessor_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=0;
  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().c_str(),"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().c_str();
      if(summaryhisto.find(layerid)==(summaryhisto.end())){
    folder_organizer.setSiStripFolder();
    MonitorElement * summaryprofile=0;
    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  )

virtual

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