#include <DebugTools/OverlapProblem/plugins/OverlapProblemTSOSAnalyzer.cc>

Inheritance diagram for OverlapProblemTSOSAnalyzer:

Public Member Functions
	OverlapProblemTSOSAnalyzer (const edm::ParameterSet &)

	~OverlapProblemTSOSAnalyzer ()

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

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override

virtual void	beginRun (const edm::Run &, const edm::EventSetup &) override

virtual void	endRun (const edm::Run &, const edm::EventSetup &) override

Private Attributes
const bool	m_debug

TH1F *	m_etatrk

TH1F *	m_ptrk

TSOSHistogramMaker	m_tsoshm

edm::EDGetTokenT < TrajTrackAssociationCollection >	m_ttacollToken

bool	m_validOnly

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

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 67 of file OverlapProblemTSOSAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 101 of file OverlapProblemTSOSAnalyzer.cc.

References m_etatrk, m_ptrk, and TFileService::make().

                                                                                     :
   m_validOnly(iConfig.getParameter<bool>("onlyValidRecHit")),
   m_ttacollToken(consumes<TrajTrackAssociationCollection>(iConfig.getParameter<edm::InputTag>("trajTrackAssoCollection"))),
   m_debug(iConfig.getUntrackedParameter<bool>("debugMode",false)),
   m_tsoshm(iConfig.getParameter<edm::ParameterSet>("tsosHMConf"))
 {
    //now do what ever initialization is needed
 
 
 
   edm::Service<TFileService> tfserv;
 
   m_ptrk = tfserv->make<TH1F>("trkmomentum","Refitted Track  momentum",100,0.,200.);
   m_etatrk = tfserv->make<TH1F>("trketa","Refitted Track pseudorapidity",100,-4.,4.);
 
 
 }

OverlapProblemTSOSAnalyzer::~OverlapProblemTSOSAnalyzer ( )

Definition at line 120 of file OverlapProblemTSOSAnalyzer.cc.

 {
 
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 135 of file OverlapProblemTSOSAnalyzer.cc.

References edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, TSOSHistogramMaker::fill(), edm::Event::getByToken(), DetIdSelector::isSelected(), edm::Ref< C, T, F >::key(), m_debug, m_etatrk, m_ptrk, m_tsoshm, m_ttacollToken, and m_validOnly.

 {
   using namespace edm;
 
   // loop on trajectories and plot TSOS local coordinate
   
   TrajectoryStateCombiner tsoscomb;
   
   // Trajectory Handle
 
   DetIdSelector selector("0x1fbff004-0x14ac1004");
   
   Handle<TrajTrackAssociationCollection> ttac;
   iEvent.getByToken(m_ttacollToken,ttac);
   
   for(TrajTrackAssociationCollection::const_iterator pair=ttac->begin();pair!=ttac->end();++pair) {
     
     const edm::Ref<std::vector<Trajectory> > & traj = pair->key;
     const reco::TrackRef & trk = pair->val;
     const std::vector<TrajectoryMeasurement> & tmcoll = traj->measurements();
     
     m_ptrk->Fill(trk->p());
     m_etatrk->Fill(trk->eta());
 
 
     for(std::vector<TrajectoryMeasurement>::const_iterator meas = tmcoll.begin() ; meas!= tmcoll.end() ; ++meas) {
       
       if(!meas->updatedState().isValid()) continue;
       
       TrajectoryStateOnSurface tsos = tsoscomb(meas->forwardPredictedState(), meas->backwardPredictedState());
       TransientTrackingRecHit::ConstRecHitPointer hit = meas->recHit();
       
       m_tsoshm.fill(tsos,hit);
 
       if(!hit->isValid() && m_validOnly) continue;
 
       if(m_debug) {
     if(selector.isSelected(hit->geographicalId())) {
       const SiPixelRecHit* pixelrh = dynamic_cast<const SiPixelRecHit*>(hit->hit());
       if(pixelrh) {
         edm::LogInfo("ClusterFound") << "Cluster reference" << pixelrh->cluster().key(); 
       }
       else {
         edm::LogInfo("NoCluster") << "No cluster found!";
       }
     }
       }
 
     }
     
   }
   
   
 }