#include <MuonTrackAnalyzer.h>

Inheritance diagram for MuonTrackAnalyzer:

Public Types
enum	EtaRange { all, barrel, endcap }

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup)

virtual void	beginJob ()

virtual void	beginRun ()

virtual void	endJob ()

virtual void	endRun ()

	MuonTrackAnalyzer (const edm::ParameterSet &pset)
	Constructor. More...

void	seedsAnalysis (const edm::Event &event, const edm::EventSetup &eventSetup, edm::Handle< edm::SimTrackContainer > simTracks)

void	tracksAnalysis (const edm::Event &event, const edm::EventSetup &eventSetup, edm::Handle< edm::SimTrackContainer > simTracks)

virtual	~MuonTrackAnalyzer ()
	Destructor. More...

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
bool	checkMuonSimHitPresence (const edm::Event &event, edm::Handle< edm::SimTrackContainer > simTracks)

void	fillPlots (const edm::Event &event, edm::Handle< edm::SimTrackContainer > &simTracks)

void	fillPlots (reco::TransientTrack &track, SimTrack &simTrack)

void	fillPlots (TrajectoryStateOnSurface &recoTSOS, SimTrack &simState, HTrack *, MuonPatternRecoDumper &)

void	fillPlots (FreeTrajectoryState &recoFTS, SimTrack &simTrack, HTrack *histo, MuonPatternRecoDumper &debug)

TrajectoryStateOnSurface	getSeedTSOS (const TrajectorySeed &seed)

std::pair< SimTrack, double >	getSimTrack (TrajectoryStateOnSurface &tsos, edm::Handle< edm::SimTrackContainer > simTracks)

bool	isInTheAcceptance (double eta)

Private Attributes
DQMStore *	dbe_

std::string	dirName_

bool	doSeedsAnalysis

bool	doTracksAnalysis

MonitorElement *	hChargeVsEta

MonitorElement *	hChargeVsPt

MonitorElement *	hChi2

MonitorElement *	hChi2Norm

MonitorElement *	hChi2NormVsEta

MonitorElement *	hChi2Prob

MonitorElement *	hChi2ProbVsEta

MonitorElement *	hChi2VsEta

MonitorElement *	hDeltaPt_In_Out_VsEta

MonitorElement *	hDeltaPtVsEta

MonitorElement *	hDof

MonitorElement *	hDofVsEta

MonitorElement *	hHitsPerTrack

MonitorElement *	hHitsPerTrackVsEta

MonitorElement *	hNumberOfTracks

MonitorElement *	hNumberOfTracksVsEta

MonitorElement *	hPtRecVsPtGen

HTrack *	hRecoSeedInner

HTrack *	hRecoSeedPCA

HTrack *	hRecoTracksInner

HTrack *	hRecoTracksOuter

HTrack *	hRecoTracksPCA

HTrackVariables *	hSimTracks

int	numberOfRecTracks

int	numberOfSimTracks

std::string	out

edm::InputTag	theCSCSimHitLabel

edm::EDGetTokenT< std::vector < PSimHit > >	theCSCSimHitToken

edm::InputTag	theDTSimHitLabel

edm::EDGetTokenT< std::vector < PSimHit > >	theDTSimHitToken

EtaRange	theEtaRange

edm::InputTag	theRPCSimHitLabel

edm::EDGetTokenT< std::vector < PSimHit > >	theRPCSimHitToken

std::string	theSeedPropagatorName

edm::InputTag	theSeedsLabel

edm::EDGetTokenT < TrajectorySeedCollection >	theSeedsToken

MuonServiceProxy *	theService

edm::InputTag	theSimTracksLabel

edm::EDGetTokenT < edm::SimTrackContainer >	theSimTracksToken

edm::InputTag	theTracksLabel

edm::EDGetTokenT < reco::TrackCollection >	theTracksToken

MuonUpdatorAtVertex *	theUpdator

Additional Inherited Members
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

Analyzer of the Muon tracks

Author: R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Analyzer of the StandAlone muon tracks

Author: R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 45 of file MuonTrackAnalyzer.h.

Member Enumeration Documentation

enum MuonTrackAnalyzer::EtaRange

Enumerator
all
barrel
endcap

Definition at line 48 of file MuonTrackAnalyzer.h.

48 {all,barrel,endcap};

MuonTrackAnalyzer::endcap

Definition: MuonTrackAnalyzer.h:48

MuonTrackAnalyzer::barrel

Definition: MuonTrackAnalyzer.h:48

MuonTrackAnalyzer::all

Definition: MuonTrackAnalyzer.h:48

Constructor & Destructor Documentation

MuonTrackAnalyzer::MuonTrackAnalyzer ( const edm::ParameterSet & pset )

Constructor.

Definition at line 40 of file MuonTrackAnalyzer.cc.

References dbe_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), MuonServiceProxy_cff::MuonServiceProxy, MuonUpdatorAtVertex_cff::MuonUpdatorAtVertex, cppFunctionSkipper::operator, dbtoconf::out, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                             {
   
   // service parameters
   ParameterSet serviceParameters = pset.getParameter<ParameterSet>("ServiceParameters");
   // the services
   theService = new MuonServiceProxy(serviceParameters);
   
   theSimTracksLabel = edm::InputTag("g4SimHits");
   theSimTracksToken = consumes<edm::SimTrackContainer>(theSimTracksLabel);
 
   theTracksLabel = pset.getParameter<InputTag>("Tracks");
   theTracksToken = consumes<reco::TrackCollection>(theTracksLabel);
   doTracksAnalysis = pset.getUntrackedParameter<bool>("DoTracksAnalysis",true);
 
   doSeedsAnalysis = pset.getUntrackedParameter<bool>("DoSeedsAnalysis",false);
   if(doSeedsAnalysis){
     theSeedsLabel = pset.getParameter<InputTag>("MuonSeed");
     theSeedsToken =  consumes<TrajectorySeedCollection>(theSeedsLabel);
     ParameterSet updatorPar = pset.getParameter<ParameterSet>("MuonUpdatorAtVertexParameters");
     theSeedPropagatorName = updatorPar.getParameter<string>("Propagator");
 
     theUpdator = new MuonUpdatorAtVertex(updatorPar,theService);
   }
   
   theCSCSimHitLabel = pset.getParameter<InputTag>("CSCSimHit");
   theDTSimHitLabel = pset.getParameter<InputTag>("DTSimHit");
   theRPCSimHitLabel = pset.getParameter<InputTag>("RPCSimHit");
   theCSCSimHitToken = consumes<std::vector<PSimHit> >(theCSCSimHitLabel);
   theDTSimHitToken = consumes<std::vector<PSimHit> >(theDTSimHitLabel);
   theRPCSimHitToken = consumes<std::vector<PSimHit> >(theRPCSimHitLabel);
 
   theEtaRange = (EtaRange) pset.getParameter<int>("EtaRange");
   
   // number of sim tracks
   numberOfSimTracks=0;
   // number of reco tracks
   numberOfRecTracks=0;
 
   dbe_ = edm::Service<DQMStore>().operator->();
   out = pset.getUntrackedParameter<string>("rootFileName");
   dirName_ = pset.getUntrackedParameter<std::string>("dirName");
 
 }

MuonTrackAnalyzer::~MuonTrackAnalyzer ( )

virtual

Destructor.

Definition at line 85 of file MuonTrackAnalyzer.cc.

                                      {
   if (theService) delete theService;
 }

Member Function Documentation

void MuonTrackAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

virtual

Implements edm::EDAnalyzer.

Definition at line 193 of file MuonTrackAnalyzer.cc.

References LogDebug.

                                                                                 {
   
   LogDebug("MuonTrackAnalyzer") << "Run: " << event.id().run() << " Event: " << event.id().event();
 
   // Update the services
   theService->update(eventSetup);
 
   Handle<SimTrackContainer> simTracks;
   event.getByToken(theSimTracksToken,simTracks);  
   fillPlots(event,simTracks);
 
   
   if(doTracksAnalysis)
     tracksAnalysis(event,eventSetup,simTracks);
   
   if(doSeedsAnalysis)
     seedsAnalysis(event,eventSetup,simTracks);
   
 
 }

void MuonTrackAnalyzer::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 89 of file MuonTrackAnalyzer.cc.

                                 {
 
   //theFile->cd();
 }

void MuonTrackAnalyzer::beginRun ( void )

virtual

Definition at line 93 of file MuonTrackAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::cd(), dbe_, TrackerOfflineValidation_Dqm_cff::dirName, hChi2, edm::InputTag::instance(), edm::InputTag::label(), edm::InputTag::process(), python.rootplot.root2matplotlib::replace(), DQMStore::setCurrentFolder(), DQMStore::showDirStructure(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                 {
   dbe_->showDirStructure();
   
   dbe_->cd();
   InputTag algo = theTracksLabel;
   string dirName=dirName_;
   if (algo.process()!="")
     dirName+=algo.process()+"_";
   if(algo.label()!="")
     dirName+=algo.label()+"_";
   if(algo.instance()!="")
     dirName+=algo.instance()+"";      
   if (dirName.find("Tracks")<dirName.length()){
     dirName.replace(dirName.find("Tracks"),6,"");
   }
   std::replace(dirName.begin(), dirName.end(), ':', '_');
   dbe_->setCurrentFolder(dirName.c_str());
   
   //dbe_->goUp();
   std::string simName = dirName;
   simName+="/SimTracks";
   hSimTracks = new HTrackVariables(simName.c_str(),"SimTracks"); 
 
   dbe_->cd();
   dbe_->setCurrentFolder(dirName.c_str());
 
   // Create the root file
   //theFile = new TFile(theRootFileName.c_str(), "RECREATE");
 
   if(doSeedsAnalysis){
     dbe_->cd();
     dbe_->setCurrentFolder(dirName.c_str());
     hRecoSeedInner = new HTrack(dirName.c_str(),"RecoSeed","Inner");
     hRecoSeedPCA = new HTrack(dirName.c_str(),"RecoSeed","PCA");
   }
   
   if(doTracksAnalysis){
     dbe_->cd();
     dbe_->setCurrentFolder(dirName.c_str());    
     hRecoTracksPCA = new HTrack(dirName.c_str(),"RecoTracks","PCA"); 
     hRecoTracksInner = new HTrack(dirName.c_str(),"RecoTracks","Inner"); 
     hRecoTracksOuter = new HTrack(dirName.c_str(),"RecoTracks","Outer"); 
 
     dbe_->cd();
     dbe_->setCurrentFolder(dirName.c_str());
     
     // General Histos
 
 
     hChi2 = dbe_->book1D("chi2","#chi^2",200,0,200);
     hChi2VsEta = dbe_->book2D("chi2VsEta","#chi^2 VS #eta",120,-3.,3.,200,0,200);
     
     hChi2Norm = dbe_->book1D("chi2Norm","Normalized #chi^2",400,0,100);
     hChi2NormVsEta = dbe_->book2D("chi2NormVsEta","Normalized #chi^2 VS #eta",120,-3.,3.,400,0,100);
     
     hHitsPerTrack  = dbe_->book1D("HitsPerTrack","Number of hits per track",55,0,55);
     hHitsPerTrackVsEta  = dbe_->book2D("HitsPerTrackVsEta","Number of hits per track VS #eta",
                    120,-3.,3.,55,0,55);
     
     hDof  = dbe_->book1D("dof","Number of Degree of Freedom",55,0,55);
     hDofVsEta  = dbe_->book2D("dofVsEta","Number of Degree of Freedom VS #eta",120,-3.,3.,55,0,55);
     
     hChi2Prob = dbe_->book1D("chi2Prob","#chi^2 probability",200,0,1);
     hChi2ProbVsEta = dbe_->book2D("chi2ProbVsEta","#chi^2 probability VS #eta",120,-3.,3.,200,0,1);
     
     hNumberOfTracks = dbe_->book1D("NumberOfTracks","Number of reconstructed tracks per event",200,0,200);
     hNumberOfTracksVsEta = dbe_->book2D("NumberOfTracksVsEta",
                     "Number of reconstructed tracks per event VS #eta",
                     120,-3.,3.,10,0,10);
     
     hChargeVsEta = dbe_->book2D("ChargeVsEta","Charge vs #eta gen",120,-3.,3.,4,-2.,2.);
     hChargeVsPt = dbe_->book2D("ChargeVsPt","Charge vs P_{T} gen",250,0,200,4,-2.,2.);
     hPtRecVsPtGen = dbe_->book2D("PtRecVsPtGen","P_{T} rec vs P_{T} gen",250,0,200,250,0,200);
 
     hDeltaPtVsEta = dbe_->book2D("DeltaPtVsEta","#Delta P_{t} vs #eta gen",120,-3.,3.,500,-250.,250.);
     hDeltaPt_In_Out_VsEta = dbe_->book2D("DeltaPt_In_Out_VsEta_","P^{in}_{t} - P^{out}_{t} vs #eta gen",120,-3.,3.,500,-250.,250.);
   }    
 
 }

bool MuonTrackAnalyzer::checkMuonSimHitPresence	(	const edm::Event &	event,
		edm::Handle< edm::SimTrackContainer >	simTracks
	)

private

Definition at line 435 of file MuonTrackAnalyzer.cc.

References funct::abs().

                                                                   {
 
   // Get the SimHit collection from the event
   Handle<PSimHitContainer> dtSimHits;
   event.getByToken(theDTSimHitToken, dtSimHits);
   //  event.getByToken(theDTSimHitLabel.instance(),theDTSimHitLabel.label(), dtSimHits);
   
   Handle<PSimHitContainer> cscSimHits;
   event.getByToken(theCSCSimHitToken, cscSimHits);
   
   Handle<PSimHitContainer> rpcSimHits;
   event.getByToken(theRPCSimHitToken, rpcSimHits);  
   
   map<unsigned int, vector<const PSimHit*> > mapOfMuonSimHits;
   
   for(PSimHitContainer::const_iterator simhit = dtSimHits->begin();
       simhit != dtSimHits->end(); ++simhit) {
     if (abs(simhit->particleType()) != 13) continue;
     mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
   }
   
   for(PSimHitContainer::const_iterator simhit = cscSimHits->begin();
       simhit != cscSimHits->end(); ++simhit) {
     if (abs(simhit->particleType()) != 13) continue;
     mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
   }
   
   for(PSimHitContainer::const_iterator simhit = rpcSimHits->begin();
       simhit != rpcSimHits->end(); ++simhit) {
     if (abs(simhit->particleType()) != 13) continue;
     mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
   }
 
   bool presence = false;
 
   for (SimTrackContainer::const_iterator simTrack = simTracks->begin(); 
        simTrack != simTracks->end(); ++simTrack){
     
     if (abs(simTrack->type()) != 13) continue;
     
     map<unsigned int, vector<const PSimHit*> >::const_iterator mapIterator = 
       mapOfMuonSimHits.find(simTrack->trackId());
     
     if (mapIterator != mapOfMuonSimHits.end())
       presence = true;
   }
   
   return presence;
 }

void MuonTrackAnalyzer::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 173 of file MuonTrackAnalyzer.cc.

                               {
 
 }

void MuonTrackAnalyzer::endRun ( void )

virtual

Definition at line 176 of file MuonTrackAnalyzer.cc.

References dbe_, dbtoconf::out, and DQMStore::save().

                               {
   LogInfo("MuonTrackAnalyzer")<< "Number of Sim tracks: " << numberOfSimTracks;
 
   LogInfo("MuonTrackAnalyzer") << "Number of Reco tracks: " << numberOfRecTracks;
 
   
   if(doTracksAnalysis){
     double eff = hRecoTracksPCA->computeEfficiency(hSimTracks);
     LogInfo("MuonTrackAnalyzer") <<" *Track Efficiency* = "<< eff <<"%";
   }
 
   if(doSeedsAnalysis){
     double eff = hRecoSeedInner->computeEfficiency(hSimTracks);
     LogInfo("MuonTrackAnalyzer")<<" *Seed Efficiency* = "<< eff <<"%";
   }
   if ( out.size() != 0 && dbe_ ) dbe_->save(out);
 }

void MuonTrackAnalyzer::fillPlots	(	const edm::Event &	event,
		edm::Handle< edm::SimTrackContainer > &	simTracks
	)

private

Definition at line 299 of file MuonTrackAnalyzer.cc.

References funct::abs(), LogTrace, and mathSSE::sqrt().

                                                                                                 {
   
   if(!checkMuonSimHitPresence(event,simTracks)) return;
   
   // Loop over the Sim tracks
   SimTrackContainer::const_iterator simTrack;
   LogTrace("MuonTrackAnalyzer")<<"Simulated tracks: "<<simTracks->size()<<endl;
   
   for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
     if (abs((*simTrack).type()) == 13) {
       
       if( !isInTheAcceptance( (*simTrack).momentum().eta()) ) continue; // FIXME!!
       
     numberOfSimTracks++;
     
     LogTrace("MuonTrackAnalyzer")<<"Simualted muon:"<<endl;
     LogTrace("MuonTrackAnalyzer")<<"Sim pT: "<<sqrt((*simTrack).momentum().perp2())<<endl;
     LogTrace("MuonTrackAnalyzer")<<"Sim Eta: "<<(*simTrack).momentum().eta()<<endl; // FIXME
     
     hSimTracks->Fill((*simTrack).momentum().mag(), 
              sqrt((*simTrack).momentum().perp2()), 
              (*simTrack).momentum().eta(), 
              (*simTrack).momentum().phi(), 
              -(*simTrack).type()/ abs((*simTrack).type()) ); // Double FIXME  
     LogTrace("MuonTrackAnalyzer") << "hSimTracks filled" << endl;
     }
   
   LogTrace("MuonTrackAnalyzer") << endl; 
 }

void MuonTrackAnalyzer::fillPlots	(	reco::TransientTrack &	track,
		SimTrack &	simTrack
	)

private

Definition at line 330 of file MuonTrackAnalyzer.cc.

References reco::TransientTrack::chi2(), ChiSquaredProbability(), hChi2, LogTrace, CoreSimTrack::momentum(), reco::TransientTrack::ndof(), reco::TransientTrack::normalizedChi2(), and reco::TransientTrack::recHitsSize().

                                                                                {
 
   LogTrace("MuonTrackAnalyzer")<<"Analizer: New track, chi2: "<<track.chi2()<<" dof: "<<track.ndof()<<endl;
   hChi2->Fill(track.chi2());
   hDof->Fill(track.ndof());
   hChi2Norm->Fill(track.normalizedChi2());
   hHitsPerTrack->Fill(track.recHitsSize());
   
   hChi2Prob->Fill( ChiSquaredProbability(track.chi2(),track.ndof()) );
 
   hChi2VsEta->Fill(simTrack.momentum().eta(),track.chi2());
   hChi2NormVsEta->Fill(simTrack.momentum().eta(),track.normalizedChi2());
   hChi2ProbVsEta->Fill(simTrack.momentum().eta(),ChiSquaredProbability(track.chi2(),track.ndof()));     
   hHitsPerTrackVsEta->Fill(simTrack.momentum().eta(),track.recHitsSize());
   hDofVsEta->Fill(simTrack.momentum().eta(),track.ndof()); 
 }

void MuonTrackAnalyzer::fillPlots	(	TrajectoryStateOnSurface &	recoTSOS,
		SimTrack &	simState,
		HTrack *	histo,
		MuonPatternRecoDumper &	debug
	)

private

Definition at line 348 of file MuonTrackAnalyzer.cc.

References HTrack::computeResolutionAndPull(), MuonPatternRecoDumper::dumpTSOS(), HTrack::Fill(), HTrack::FillDeltaR(), TrajectoryStateOnSurface::globalMomentum(), LogTrace, CoreSimTrack::momentum(), dt_dqm_sourceclient_common_cff::reco, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                {
   
   LogTrace("MuonTrackAnalyzer") << debug.dumpTSOS(recoTSOS)<<endl;
   histo->Fill(recoTSOS);
   
   GlobalVector tsosVect = recoTSOS.globalMomentum();
   math::XYZVectorD reco(tsosVect.x(), tsosVect.y(), tsosVect.z());
   double deltaRVal = deltaR<double>(reco.eta(),reco.phi(),
                     simTrack.momentum().eta(),simTrack.momentum().phi());
   histo->FillDeltaR(deltaRVal);
 
   histo->computeResolutionAndPull(recoTSOS,simTrack);
 }

void MuonTrackAnalyzer::fillPlots	(	FreeTrajectoryState &	recoFTS,
		SimTrack &	simTrack,
		HTrack *	histo,
		MuonPatternRecoDumper &	debug
	)

private

Definition at line 364 of file MuonTrackAnalyzer.cc.

References HTrack::computeResolutionAndPull(), MuonPatternRecoDumper::dumpFTS(), HTrack::Fill(), HTrack::FillDeltaR(), LogTrace, CoreSimTrack::momentum(), FreeTrajectoryState::momentum(), dt_dqm_sourceclient_common_cff::reco, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                {
   
   LogTrace("MuonTrackAnalyzer") << debug.dumpFTS(recoFTS)<<endl;
   histo->Fill(recoFTS);
   
   GlobalVector ftsVect = recoFTS.momentum();
   math::XYZVectorD reco(ftsVect.x(), ftsVect.y(), ftsVect.z());
   double deltaRVal = deltaR<double>(reco.eta(),reco.phi(),
                     simTrack.momentum().eta(),simTrack.momentum().phi());
   histo->FillDeltaR(deltaRVal);
 
   histo->computeResolutionAndPull(recoFTS,simTrack);
 }

TrajectoryStateOnSurface MuonTrackAnalyzer::getSeedTSOS ( const TrajectorySeed & seed )

private

Definition at line 486 of file MuonTrackAnalyzer.cc.

References PTrajectoryStateOnDet::detId(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::isValid(), oppositeToMomentum, query::result, TrajectorySeed::startingState(), GeometricSearchDet::surface(), and trajectoryStateTransform::transientState().

                                                                                  {
 
   // Get the Trajectory State on Det (persistent version of a TSOS) from the seed
   PTrajectoryStateOnDet pTSOD = seed.startingState();
 
   // Transform it in a TrajectoryStateOnSurface
   
 
   DetId seedDetId(pTSOD.detId());
 
   const GeomDet* gdet = theService->trackingGeometry()->idToDet( seedDetId );
 
   TrajectoryStateOnSurface initialState = trajectoryStateTransform::transientState(pTSOD, &(gdet->surface()), &*theService->magneticField());
 
   // Get the layer on which the seed relies
   const DetLayer *initialLayer = theService->detLayerGeometry()->idToLayer( seedDetId );
 
   PropagationDirection detLayerOrder = oppositeToMomentum;
 
   // ask for compatible layers
   vector<const DetLayer*> detLayers;
   detLayers = initialLayer->compatibleLayers( *initialState.freeState(),detLayerOrder);
   
   TrajectoryStateOnSurface result = initialState;
   if(detLayers.size()){
     const DetLayer* finalLayer = detLayers.back();
     const TrajectoryStateOnSurface propagatedState = theService->propagator(theSeedPropagatorName)->propagate(initialState, finalLayer->surface());
     if(propagatedState.isValid())
       result = propagatedState;
   }
   
   return result;
 }

pair< SimTrack, double > MuonTrackAnalyzer::getSimTrack	(	TrajectoryStateOnSurface &	tsos,
		edm::Handle< edm::SimTrackContainer >	simTracks
	)

private

Definition at line 379 of file MuonTrackAnalyzer.cc.

References funct::abs(), TrajectoryStateOnSurface::globalMomentum(), LogTrace, query::result, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                  {
   
 //   // Loop over the Sim tracks
 //   SimTrackContainer::const_iterator simTrack;
   
 //   SimTrack result;
 //   int mu=0;
 //   for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
 //     if (abs((*simTrack).type()) == 13) { 
 //       result = *simTrack;
 //       ++mu;
 //     }
   
 //   if(mu != 1) LogTrace("MuonTrackAnalyzer") << "WARNING!! more than 1 simulated muon!!" <<endl;
 //   return result;
 
 
   // Loop over the Sim tracks
   SimTrackContainer::const_iterator simTrack;
   
   SimTrack result;
 
   double bestDeltaR = 10e5;
   for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack){
     if (abs((*simTrack).type()) != 13) continue;
     
     //    double newDeltaR = tsos.globalMomentum().basicVector().deltaR(simTrack->momentum().vect());
     GlobalVector tsosVect = tsos.globalMomentum();
     math::XYZVectorD vect(tsosVect.x(), tsosVect.y(), tsosVect.z());
     double newDeltaR = deltaR<double>(vect.eta(),vect.phi(),
                       simTrack->momentum().eta(),simTrack->momentum().phi());
 
     if (  newDeltaR < bestDeltaR ) {
       LogTrace("MuonTrackAnalyzer") << "Matching Track with DeltaR = " << newDeltaR<<endl;
       bestDeltaR = newDeltaR;
       result  = *simTrack;
     }
   } 
   return pair<SimTrack,double>(result,bestDeltaR);
 }

bool MuonTrackAnalyzer::isInTheAcceptance ( double eta )

private

Definition at line 422 of file MuonTrackAnalyzer.cc.

References funct::abs(), cond::ecalcond::all, Reference_intrackfit_cff::barrel, Reference_intrackfit_cff::endcap, and LogTrace.

                                                    {
   switch(theEtaRange){
   case all:
     return ( abs(eta) <= 2.4 ) ? true : false;
   case barrel:
     return ( abs(eta) < 1.1 ) ? true : false;
   case endcap:
     return ( abs(eta) >= 1.1 && abs(eta) <= 2.4 ) ? true : false;  
   default:
     {LogTrace("MuonTrackAnalyzer")<<"No correct Eta range selected!! "<<endl; return false;}
   }
 }

void MuonTrackAnalyzer::seedsAnalysis	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup,
		edm::Handle< edm::SimTrackContainer >	simTracks
	)

Definition at line 214 of file MuonTrackAnalyzer.cc.

References debug, and LogTrace.

                                                           {
 
   MuonPatternRecoDumper debug;
 
   // Get the RecTrack collection from the event
   Handle<TrajectorySeedCollection> seeds;
   event.getByToken(theSeedsToken, seeds);
   
   LogTrace("MuonTrackAnalyzer")<<"Number of reconstructed seeds: " << seeds->size()<<endl;
 
   for(TrajectorySeedCollection::const_iterator seed = seeds->begin(); 
       seed != seeds->end(); ++seed){
     TrajectoryStateOnSurface seedTSOS = getSeedTSOS(*seed);
     pair<SimTrack,double> sim = getSimTrack(seedTSOS,simTracks);
     fillPlots(seedTSOS, sim.first,
           hRecoSeedInner, debug);
     
     std::pair<bool,FreeTrajectoryState> propSeed =
       theUpdator->propagateToNominalLine(seedTSOS);
     if(propSeed.first)
       fillPlots(propSeed.second, sim.first,
         hRecoSeedPCA, debug);
     else
       LogTrace("MuonTrackAnalyzer")<<"Error in seed propagation"<<endl;
    
   }
 }

void MuonTrackAnalyzer::tracksAnalysis	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup,
		edm::Handle< edm::SimTrackContainer >	simTracks
	)

Definition at line 244 of file MuonTrackAnalyzer.cc.

References TrajectoryStateOnSurface::charge(), debug, TrajectoryStateOnSurface::globalMomentum(), LogTrace, CoreSimTrack::momentum(), PV3DBase< T, PVType, FrameType >::perp(), mathSSE::sqrt(), lumiQTWidget::t, and testEve_cfg::tracks.

                                                           {
   MuonPatternRecoDumper debug;
   
   
   // Get the RecTrack collection from the event
   Handle<reco::TrackCollection> tracks;
   event.getByToken(theTracksToken, tracks);
 
   LogTrace("MuonTrackAnalyzer")<<"Reconstructed tracks: " << tracks->size() << endl;
   hNumberOfTracks->Fill(tracks->size());
   
   if(tracks->size()) numberOfRecTracks++;
   
   // Loop over the Rec tracks
   for(reco::TrackCollection::const_iterator t = tracks->begin(); t != tracks->end(); ++t) {
     
     reco::TransientTrack track(*t,&*theService->magneticField(),theService->trackingGeometry()); 
 
     TrajectoryStateOnSurface outerTSOS = track.outermostMeasurementState();
     TrajectoryStateOnSurface innerTSOS = track.innermostMeasurementState();
     TrajectoryStateOnSurface pcaTSOS   = track.impactPointState();
 
     pair<SimTrack,double> sim = getSimTrack(pcaTSOS,simTracks);
     SimTrack simTrack = sim.first;
     hNumberOfTracksVsEta->Fill(simTrack.momentum().eta(), tracks->size());
     fillPlots(track,simTrack);
     
     LogTrace("MuonTrackAnalyzer") << "State at the outer surface: " << endl; 
     fillPlots(outerTSOS,simTrack,hRecoTracksOuter,debug);
 
     LogTrace("MuonTrackAnalyzer") << "State at the inner surface: " << endl; 
     fillPlots(innerTSOS,simTrack,hRecoTracksInner,debug);
 
     LogTrace("MuonTrackAnalyzer") << "State at PCA: " << endl; 
     fillPlots(pcaTSOS,simTrack,hRecoTracksPCA,debug);
     
     double deltaPt_in_out = innerTSOS.globalMomentum().perp()-outerTSOS.globalMomentum().perp();
     hDeltaPt_In_Out_VsEta->Fill(simTrack.momentum().eta(),deltaPt_in_out);
 
     double deltaPt_pca_sim = pcaTSOS.globalMomentum().perp()-sqrt(simTrack.momentum().Perp2());
     hDeltaPtVsEta->Fill(simTrack.momentum().eta(),deltaPt_pca_sim);
     
     hChargeVsEta->Fill(simTrack.momentum().eta(),pcaTSOS.charge());
     
     hChargeVsPt->Fill(sqrt(simTrack.momentum().perp2()),pcaTSOS.charge());
     
     hPtRecVsPtGen->Fill(sqrt(simTrack.momentum().perp2()),pcaTSOS.globalMomentum().perp());    
   }
   LogTrace("MuonTrackAnalyzer")<<"--------------------------------------------"<<endl;  
 }

Member Data Documentation

DQMStore* MuonTrackAnalyzer::dbe_

private

Definition at line 90 of file MuonTrackAnalyzer.h.

std::string MuonTrackAnalyzer::dirName_

private

Definition at line 91 of file MuonTrackAnalyzer.h.

bool MuonTrackAnalyzer::doSeedsAnalysis

private

Definition at line 114 of file MuonTrackAnalyzer.h.

bool MuonTrackAnalyzer::doTracksAnalysis

private

Definition at line 113 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChargeVsEta

private

Definition at line 129 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChargeVsPt

private

Definition at line 130 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChi2

private

Definition at line 121 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChi2Norm

private

Definition at line 122 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChi2NormVsEta

private

Definition at line 134 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChi2Prob

private

Definition at line 125 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChi2ProbVsEta

private

Definition at line 137 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hChi2VsEta

private

Definition at line 133 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hDeltaPt_In_Out_VsEta

private

Definition at line 139 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hDeltaPtVsEta

private

Definition at line 138 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hDof

private

Definition at line 124 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hDofVsEta

private

Definition at line 136 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hHitsPerTrack

private

Definition at line 123 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hHitsPerTrackVsEta

private

Definition at line 135 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hNumberOfTracks

private

Definition at line 127 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hNumberOfTracksVsEta

private

Definition at line 128 of file MuonTrackAnalyzer.h.

MonitorElement* MuonTrackAnalyzer::hPtRecVsPtGen

private

Definition at line 131 of file MuonTrackAnalyzer.h.

HTrack* MuonTrackAnalyzer::hRecoSeedInner

private

Definition at line 143 of file MuonTrackAnalyzer.h.

HTrack* MuonTrackAnalyzer::hRecoSeedPCA

private

Definition at line 144 of file MuonTrackAnalyzer.h.

HTrack* MuonTrackAnalyzer::hRecoTracksInner

private

Definition at line 146 of file MuonTrackAnalyzer.h.

HTrack* MuonTrackAnalyzer::hRecoTracksOuter

private

Definition at line 147 of file MuonTrackAnalyzer.h.

HTrack* MuonTrackAnalyzer::hRecoTracksPCA

private

Definition at line 145 of file MuonTrackAnalyzer.h.

HTrackVariables* MuonTrackAnalyzer::hSimTracks

private

Definition at line 141 of file MuonTrackAnalyzer.h.

int MuonTrackAnalyzer::numberOfRecTracks

private

Definition at line 151 of file MuonTrackAnalyzer.h.

int MuonTrackAnalyzer::numberOfSimTracks

private

Definition at line 150 of file MuonTrackAnalyzer.h.

std::string MuonTrackAnalyzer::out

private

Definition at line 93 of file MuonTrackAnalyzer.h.

edm::InputTag MuonTrackAnalyzer::theCSCSimHitLabel

private

Definition at line 101 of file MuonTrackAnalyzer.h.

edm::EDGetTokenT<std::vector<PSimHit> > MuonTrackAnalyzer::theCSCSimHitToken

private

Definition at line 108 of file MuonTrackAnalyzer.h.

edm::InputTag MuonTrackAnalyzer::theDTSimHitLabel

private

Definition at line 102 of file MuonTrackAnalyzer.h.

edm::EDGetTokenT<std::vector<PSimHit> > MuonTrackAnalyzer::theDTSimHitToken

private

Definition at line 109 of file MuonTrackAnalyzer.h.

EtaRange MuonTrackAnalyzer::theEtaRange

private

Definition at line 96 of file MuonTrackAnalyzer.h.

edm::InputTag MuonTrackAnalyzer::theRPCSimHitLabel

private

Definition at line 103 of file MuonTrackAnalyzer.h.

edm::EDGetTokenT<std::vector<PSimHit> > MuonTrackAnalyzer::theRPCSimHitToken

private

Definition at line 110 of file MuonTrackAnalyzer.h.

std::string MuonTrackAnalyzer::theSeedPropagatorName

private

Definition at line 115 of file MuonTrackAnalyzer.h.

edm::InputTag MuonTrackAnalyzer::theSeedsLabel

private

Definition at line 99 of file MuonTrackAnalyzer.h.

edm::EDGetTokenT<TrajectorySeedCollection> MuonTrackAnalyzer::theSeedsToken

private

Definition at line 106 of file MuonTrackAnalyzer.h.

MuonServiceProxy* MuonTrackAnalyzer::theService

private

Definition at line 117 of file MuonTrackAnalyzer.h.

edm::InputTag MuonTrackAnalyzer::theSimTracksLabel

private

Definition at line 98 of file MuonTrackAnalyzer.h.

edm::EDGetTokenT<edm::SimTrackContainer> MuonTrackAnalyzer::theSimTracksToken

private

Definition at line 105 of file MuonTrackAnalyzer.h.

edm::InputTag MuonTrackAnalyzer::theTracksLabel

private

Definition at line 100 of file MuonTrackAnalyzer.h.

edm::EDGetTokenT<reco::TrackCollection> MuonTrackAnalyzer::theTracksToken

private

Definition at line 107 of file MuonTrackAnalyzer.h.

MuonUpdatorAtVertex* MuonTrackAnalyzer::theUpdator

private

Definition at line 118 of file MuonTrackAnalyzer.h.

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation