MuonTrackAnalyzer Class Reference

Analyzer of the Muon tracks. More...

#include <Validation/RecoMuon/src/MuonTrackAnalyzer.h>

Inheritance diagram for MuonTrackAnalyzer:

Detailed Description

Analyzer of the Muon tracks.

Analyzer of the StandAlone muon tracks.

Date: 2008/08/11 12:09:39

Revision: 1.7

Author:: R. Bellan - INFN Torino <riccardo.bellan@cern.ch>

Date: 2008/06/24 16:40:53

Revision: 1.4

Author:: R. Bellan - INFN Torino <riccardo.bellan@cern.ch>

Definition at line 40 of file MuonTrackAnalyzer.h.

Member Enumeration Documentation

enum MuonTrackAnalyzer::EtaRange

Enumerator:

all
barrel
endcap

Definition at line 43 of file MuonTrackAnalyzer.h.

00043 {all,barrel,endcap};

Constructor & Destructor Documentation

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

Constructor.

Definition at line 49 of file MuonTrackAnalyzer.cc.

References dbe_, dirName_, doSeedsAnalysis, doTracksAnalysis, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), MuonServiceProxy_cff::MuonServiceProxy, MuonUpdatorAtVertex_cff::MuonUpdatorAtVertex, numberOfRecTracks, numberOfSimTracks, out, theCSCSimHitLabel, theDTSimHitLabel, theEtaRange, theRPCSimHitLabel, theSeedPropagatorName, theSeedsLabel, theService, theTracksLabel, and theUpdator.

00049                                                             {
00050   
00051   // service parameters
00052   ParameterSet serviceParameters = pset.getParameter<ParameterSet>("ServiceParameters");
00053   // the services
00054   theService = new MuonServiceProxy(serviceParameters);
00055   
00056   theTracksLabel = pset.getParameter<InputTag>("Tracks");
00057   doTracksAnalysis = pset.getUntrackedParameter<bool>("DoTracksAnalysis",true);
00058 
00059   doSeedsAnalysis = pset.getUntrackedParameter<bool>("DoSeedsAnalysis",false);
00060   if(doSeedsAnalysis){
00061     theSeedsLabel = pset.getParameter<InputTag>("MuonSeed");
00062     ParameterSet updatorPar = pset.getParameter<ParameterSet>("MuonUpdatorAtVertexParameters");
00063     theSeedPropagatorName = updatorPar.getParameter<string>("Propagator");
00064 
00065     theUpdator = new MuonUpdatorAtVertex(updatorPar,theService);
00066   }
00067   
00068   theCSCSimHitLabel = pset.getParameter<InputTag>("CSCSimHit");
00069   theDTSimHitLabel = pset.getParameter<InputTag>("DTSimHit");
00070   theRPCSimHitLabel = pset.getParameter<InputTag>("RPCSimHit");
00071 
00072   theEtaRange = (EtaRange) pset.getParameter<int>("EtaRange");
00073   
00074   // number of sim tracks
00075   numberOfSimTracks=0;
00076   // number of reco tracks
00077   numberOfRecTracks=0;
00078 
00079   dbe_ = edm::Service<DQMStore>().operator->();
00080   out = pset.getUntrackedParameter<string>("rootFileName");
00081   dirName_ = pset.getUntrackedParameter<std::string>("dirName");
00082 
00083 }

MuonTrackAnalyzer::~MuonTrackAnalyzer ( ) [virtual]

Destructor.

Definition at line 86 of file MuonTrackAnalyzer.cc.

References theService.

00086                                      {
00087   if (theService) delete theService;
00088 }

Member Function Documentation

void MuonTrackAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)			`[virtual]`

Implements edm::EDAnalyzer.

Definition at line 190 of file MuonTrackAnalyzer.cc.

References doSeedsAnalysis, doTracksAnalysis, fillPlots(), LogDebug, seedsAnalysis(), theService, and tracksAnalysis().

00190                                                                                 {
00191   
00192   LogDebug("MuonTrackAnalyzer") << "Run: " << event.id().run() << " Event: " << event.id().event();
00193 
00194   // Update the services
00195   theService->update(eventSetup);
00196 
00197   Handle<SimTrackContainer> simTracks;
00198   event.getByLabel("g4SimHits",simTracks);  
00199   fillPlots(event,simTracks);
00200 
00201   
00202   if(doTracksAnalysis)
00203     tracksAnalysis(event,eventSetup,simTracks);
00204   
00205   if(doSeedsAnalysis)
00206     seedsAnalysis(event,eventSetup,simTracks);
00207   
00208 
00209 }

void MuonTrackAnalyzer::beginJob ( const edm::EventSetup & eventSetup ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 90 of file MuonTrackAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::cd(), dbe_, test_cfg::dirName, dirName_, doSeedsAnalysis, doTracksAnalysis, hChargeVsEta, hChargeVsPt, hChi2, hChi2Norm, hChi2NormVsEta, hChi2Prob, hChi2ProbVsEta, hChi2VsEta, hDeltaPt_In_Out_VsEta, hDeltaPtVsEta, hDof, hDofVsEta, hHitsPerTrack, hHitsPerTrackVsEta, hNumberOfTracks, hNumberOfTracksVsEta, hPtRecVsPtGen, hRecoSeedInner, hRecoSeedPCA, hRecoTracksInner, hRecoTracksOuter, hRecoTracksPCA, hSimTracks, edm::InputTag::instance(), edm::InputTag::label(), edm::InputTag::process(), replace(), DQMStore::setCurrentFolder(), DQMStore::showDirStructure(), and theTracksLabel.

00090                                                             {
00091   dbe_->showDirStructure();
00092   
00093   dbe_->cd();
00094   InputTag algo = theTracksLabel;
00095   string dirName=dirName_;
00096   if (algo.process()!="")
00097     dirName+=algo.process()+"_";
00098   if(algo.label()!="")
00099     dirName+=algo.label()+"_";
00100   if(algo.instance()!="")
00101     dirName+=algo.instance()+"";      
00102   if (dirName.find("Tracks")<dirName.length()){
00103     dirName.replace(dirName.find("Tracks"),6,"");
00104   }
00105   std::replace(dirName.begin(), dirName.end(), ':', '_');
00106   dbe_->setCurrentFolder(dirName.c_str());
00107   
00108   //dbe_->goUp();
00109   std::string simName = dirName;
00110   simName+="/SimTracks";
00111   hSimTracks = new HTrackVariables(simName.c_str(),"SimTracks"); 
00112 
00113   dbe_->cd();
00114   dbe_->setCurrentFolder(dirName.c_str());
00115 
00116   // Create the root file
00117   //theFile = new TFile(theRootFileName.c_str(), "RECREATE");
00118 
00119   if(doSeedsAnalysis){
00120     dbe_->cd();
00121     dbe_->setCurrentFolder(dirName.c_str());
00122     hRecoSeedInner = new HTrack(dirName.c_str(),"RecoSeed","Inner");
00123     hRecoSeedPCA = new HTrack(dirName.c_str(),"RecoSeed","PCA");
00124   }
00125   
00126   if(doTracksAnalysis){
00127     dbe_->cd();
00128     dbe_->setCurrentFolder(dirName.c_str());    
00129     hRecoTracksPCA = new HTrack(dirName.c_str(),"RecoTracks","PCA"); 
00130     hRecoTracksInner = new HTrack(dirName.c_str(),"RecoTracks","Inner"); 
00131     hRecoTracksOuter = new HTrack(dirName.c_str(),"RecoTracks","Outer"); 
00132 
00133     dbe_->cd();
00134     dbe_->setCurrentFolder(dirName.c_str());
00135     
00136     // General Histos
00137 
00138 
00139     hChi2 = dbe_->book1D("chi2","#chi^2",200,0,200);
00140     hChi2VsEta = dbe_->book2D("chi2VsEta","#chi^2 VS #eta",120,-3.,3.,200,0,200);
00141     
00142     hChi2Norm = dbe_->book1D("chi2Norm","Normalized #chi^2",400,0,100);
00143     hChi2NormVsEta = dbe_->book2D("chi2NormVsEta","Normalized #chi^2 VS #eta",120,-3.,3.,400,0,100);
00144     
00145     hHitsPerTrack  = dbe_->book1D("HitsPerTrack","Number of hits per track",55,0,55);
00146     hHitsPerTrackVsEta  = dbe_->book2D("HitsPerTrackVsEta","Number of hits per track VS #eta",
00147                                    120,-3.,3.,55,0,55);
00148     
00149     hDof  = dbe_->book1D("dof","Number of Degree of Freedom",55,0,55);
00150     hDofVsEta  = dbe_->book2D("dofVsEta","Number of Degree of Freedom VS #eta",120,-3.,3.,55,0,55);
00151     
00152     hChi2Prob = dbe_->book1D("chi2Prob","#chi^2 probability",200,0,1);
00153     hChi2ProbVsEta = dbe_->book2D("chi2ProbVsEta","#chi^2 probability VS #eta",120,-3.,3.,200,0,1);
00154     
00155     hNumberOfTracks = dbe_->book1D("NumberOfTracks","Number of reconstructed tracks per event",200,0,200);
00156     hNumberOfTracksVsEta = dbe_->book2D("NumberOfTracksVsEta",
00157                                     "Number of reconstructed tracks per event VS #eta",
00158                                     120,-3.,3.,10,0,10);
00159     
00160     hChargeVsEta = dbe_->book2D("ChargeVsEta","Charge vs #eta gen",120,-3.,3.,4,-2.,2.);
00161     hChargeVsPt = dbe_->book2D("ChargeVsPt","Charge vs P_{T} gen",250,0,200,4,-2.,2.);
00162     hPtRecVsPtGen = dbe_->book2D("PtRecVsPtGen","P_{T} rec vs P_{T} gen",250,0,200,250,0,200);
00163 
00164     hDeltaPtVsEta = dbe_->book2D("DeltaPtVsEta","#Delta P_{t} vs #eta gen",120,-3.,3.,500,-250.,250.);
00165     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.);
00166   }    
00167 
00168   //theFile->cd();
00169 }

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

Definition at line 432 of file MuonTrackAnalyzer.cc.

References funct::abs(), edm::InputTag::instance(), edm::InputTag::label(), theCSCSimHitLabel, theDTSimHitLabel, and theRPCSimHitLabel.

Referenced by fillPlots().

00433                                                                                           {
00434 
00435   // Get the SimHit collection from the event
00436   Handle<PSimHitContainer> dtSimHits;
00437   event.getByLabel(theDTSimHitLabel.instance(),theDTSimHitLabel.label(), dtSimHits);
00438   
00439   Handle<PSimHitContainer> cscSimHits;
00440   event.getByLabel(theCSCSimHitLabel.instance(),theCSCSimHitLabel.label(), cscSimHits);
00441   
00442   Handle<PSimHitContainer> rpcSimHits;
00443   event.getByLabel(theRPCSimHitLabel.instance(),theRPCSimHitLabel.label(), rpcSimHits);  
00444   
00445   map<unsigned int, vector<const PSimHit*> > mapOfMuonSimHits;
00446   
00447   for(PSimHitContainer::const_iterator simhit = dtSimHits->begin();
00448       simhit != dtSimHits->end(); ++simhit) {
00449     if (abs(simhit->particleType()) != 13) continue;
00450     mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
00451   }
00452   
00453   for(PSimHitContainer::const_iterator simhit = cscSimHits->begin();
00454       simhit != cscSimHits->end(); ++simhit) {
00455     if (abs(simhit->particleType()) != 13) continue;
00456     mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
00457   }
00458   
00459   for(PSimHitContainer::const_iterator simhit = rpcSimHits->begin();
00460       simhit != rpcSimHits->end(); ++simhit) {
00461     if (abs(simhit->particleType()) != 13) continue;
00462     mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
00463   }
00464 
00465   bool presence = false;
00466 
00467   for (SimTrackContainer::const_iterator simTrack = simTracks->begin(); 
00468        simTrack != simTracks->end(); ++simTrack){
00469     
00470     if (abs(simTrack->type()) != 13) continue;
00471     
00472     map<unsigned int, vector<const PSimHit*> >::const_iterator mapIterator = 
00473       mapOfMuonSimHits.find(simTrack->trackId());
00474     
00475     if (mapIterator != mapOfMuonSimHits.end())
00476       presence = true;
00477   }
00478   
00479   return presence;
00480 }

void MuonTrackAnalyzer::endJob ( void ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 171 of file MuonTrackAnalyzer.cc.

References HTrack::computeEfficiency(), dbe_, doSeedsAnalysis, doTracksAnalysis, hRecoSeedInner, hRecoTracksPCA, hSimTracks, numberOfRecTracks, numberOfSimTracks, out, and DQMStore::save().

00171                               {
00172   LogInfo("MuonTrackAnalyzer")<< "Number of Sim tracks: " << numberOfSimTracks;
00173 
00174   LogInfo("MuonTrackAnalyzer") << "Number of Reco tracks: " << numberOfRecTracks;
00175 
00176   
00177   if(doTracksAnalysis){
00178     double eff = hRecoTracksPCA->computeEfficiency(hSimTracks);
00179     LogInfo("MuonTrackAnalyzer") <<" *Track Efficiency* = "<< eff <<"%";
00180   }
00181 
00182   if(doSeedsAnalysis){
00183     double eff = hRecoSeedInner->computeEfficiency(hSimTracks);
00184     LogInfo("MuonTrackAnalyzer")<<" *Seed Efficiency* = "<< eff <<"%";
00185   }
00186 
00187   if ( out.size() != 0 && dbe_ ) dbe_->save(out);
00188 }

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

Definition at line 361 of file MuonTrackAnalyzer.cc.

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

00362                                                                                {
00363   
00364   LogTrace("MuonTrackAnalyzer") << debug.dumpFTS(recoFTS)<<endl;
00365   histo->Fill(recoFTS);
00366   
00367   GlobalVector ftsVect = recoFTS.momentum();
00368   math::XYZVectorD reco(ftsVect.x(), ftsVect.y(), ftsVect.z());
00369   double deltaRVal = deltaR<double>(reco.eta(),reco.phi(),
00370                                     simTrack.momentum().eta(),simTrack.momentum().phi());
00371   histo->FillDeltaR(deltaRVal);
00372 
00373   histo->computeResolutionAndPull(recoFTS,simTrack);
00374 }

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

Definition at line 345 of file MuonTrackAnalyzer.cc.

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

00346                                                                                {
00347   
00348   LogTrace("MuonTrackAnalyzer") << debug.dumpTSOS(recoTSOS)<<endl;
00349   histo->Fill(recoTSOS);
00350   
00351   GlobalVector tsosVect = recoTSOS.globalMomentum();
00352   math::XYZVectorD reco(tsosVect.x(), tsosVect.y(), tsosVect.z());
00353   double deltaRVal = deltaR<double>(reco.eta(),reco.phi(),
00354                                     simTrack.momentum().eta(),simTrack.momentum().phi());
00355   histo->FillDeltaR(deltaRVal);
00356 
00357   histo->computeResolutionAndPull(recoTSOS,simTrack);
00358 }

void MuonTrackAnalyzer::fillPlots	(	reco::TransientTrack &	track,
		SimTrack &	simTrack
	)			`[private]`

Definition at line 327 of file MuonTrackAnalyzer.cc.

References reco::TransientTrack::chi2(), ChiSquaredProbability(), lat::endl(), MonitorElement::Fill(), hChi2, hChi2Norm, hChi2NormVsEta, hChi2Prob, hChi2ProbVsEta, hChi2VsEta, hDof, hDofVsEta, hHitsPerTrack, hHitsPerTrackVsEta, LogTrace, CoreSimTrack::momentum(), reco::TransientTrack::ndof(), reco::TransientTrack::normalizedChi2(), and reco::TransientTrack::recHitsSize().

00327                                                                                {
00328 
00329   LogTrace("MuonTrackAnalyzer")<<"Analizer: New track, chi2: "<<track.chi2()<<" dof: "<<track.ndof()<<endl;
00330   hChi2->Fill(track.chi2());
00331   hDof->Fill(track.ndof());
00332   hChi2Norm->Fill(track.normalizedChi2());
00333   hHitsPerTrack->Fill(track.recHitsSize());
00334   
00335   hChi2Prob->Fill( ChiSquaredProbability(track.chi2(),track.ndof()) );
00336 
00337   hChi2VsEta->Fill(simTrack.momentum().eta(),track.chi2());
00338   hChi2NormVsEta->Fill(simTrack.momentum().eta(),track.normalizedChi2());
00339   hChi2ProbVsEta->Fill(simTrack.momentum().eta(),ChiSquaredProbability(track.chi2(),track.ndof()));     
00340   hHitsPerTrackVsEta->Fill(simTrack.momentum().eta(),track.recHitsSize());
00341   hDofVsEta->Fill(simTrack.momentum().eta(),track.ndof()); 
00342 }

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

Definition at line 296 of file MuonTrackAnalyzer.cc.

References funct::abs(), checkMuonSimHitPresence(), lat::endl(), HTrackVariables::Fill(), hSimTracks, isInTheAcceptance(), LogTrace, numberOfSimTracks, and funct::sqrt().

Referenced by analyze(), seedsAnalysis(), and tracksAnalysis().

00296                                                                                                 {
00297   
00298   if(!checkMuonSimHitPresence(event,simTracks)) return;
00299   
00300   // Loop over the Sim tracks
00301   SimTrackContainer::const_iterator simTrack;
00302   LogTrace("MuonTrackAnalyzer")<<"Simulated tracks: "<<simTracks->size()<<endl;
00303   
00304   for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
00305     if (abs((*simTrack).type()) == 13) {
00306       
00307       if( !isInTheAcceptance( (*simTrack).momentum().eta()) ) continue; // FIXME!!
00308       
00309         numberOfSimTracks++;
00310         
00311         LogTrace("MuonTrackAnalyzer")<<"Simualted muon:"<<endl;
00312         LogTrace("MuonTrackAnalyzer")<<"Sim pT: "<<sqrt((*simTrack).momentum().perp2())<<endl;
00313         LogTrace("MuonTrackAnalyzer")<<"Sim Eta: "<<(*simTrack).momentum().eta()<<endl; // FIXME
00314         
00315         hSimTracks->Fill((*simTrack).momentum().mag(), 
00316                          sqrt((*simTrack).momentum().perp2()), 
00317                          (*simTrack).momentum().eta(), 
00318                          (*simTrack).momentum().phi(), 
00319                          -(*simTrack).type()/ abs((*simTrack).type()) ); // Double FIXME  
00320         LogTrace("MuonTrackAnalyzer") << "hSimTracks filled" << endl;
00321     }
00322   
00323   LogTrace("MuonTrackAnalyzer") << endl; 
00324 }

TrajectoryStateOnSurface MuonTrackAnalyzer::getSeedTSOS ( const TrajectorySeed & seed ) [private]

Definition at line 482 of file MuonTrackAnalyzer.cc.

References DetLayer::compatibleLayers(), PTrajectoryStateOnDet::detId(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::isValid(), oppositeToMomentum, HLT_VtxMuL3::result, TrajectorySeed::startingState(), GeometricSearchDet::surface(), theSeedPropagatorName, theService, and TrajectoryStateTransform::transientState().

Referenced by seedsAnalysis().

00482                                                                                  {
00483 
00484   // Get the Trajectory State on Det (persistent version of a TSOS) from the seed
00485   PTrajectoryStateOnDet pTSOD = seed.startingState();
00486 
00487   // Transform it in a TrajectoryStateOnSurface
00488   TrajectoryStateTransform tsTransform;
00489 
00490   DetId seedDetId(pTSOD.detId());
00491 
00492   const GeomDet* gdet = theService->trackingGeometry()->idToDet( seedDetId );
00493 
00494   TrajectoryStateOnSurface initialState = tsTransform.transientState(pTSOD, &(gdet->surface()), &*theService->magneticField());
00495 
00496   // Get the layer on which the seed relies
00497   const DetLayer *initialLayer = theService->detLayerGeometry()->idToLayer( seedDetId );
00498 
00499   PropagationDirection detLayerOrder = oppositeToMomentum;
00500 
00501   // ask for compatible layers
00502   vector<const DetLayer*> detLayers;
00503   detLayers = initialLayer->compatibleLayers( *initialState.freeState(),detLayerOrder);
00504   
00505   TrajectoryStateOnSurface result = initialState;
00506   if(detLayers.size()){
00507     const DetLayer* finalLayer = detLayers.back();
00508     const TrajectoryStateOnSurface propagatedState = theService->propagator(theSeedPropagatorName)->propagate(initialState, finalLayer->surface());
00509     if(propagatedState.isValid())
00510       result = propagatedState;
00511   }
00512   
00513   return result;
00514 }

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

Definition at line 376 of file MuonTrackAnalyzer.cc.

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

Referenced by seedsAnalysis(), and tracksAnalysis().

00377                                                                                          {
00378   
00379 //   // Loop over the Sim tracks
00380 //   SimTrackContainer::const_iterator simTrack;
00381   
00382 //   SimTrack result;
00383 //   int mu=0;
00384 //   for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
00385 //     if (abs((*simTrack).type()) == 13) { 
00386 //       result = *simTrack;
00387 //       ++mu;
00388 //     }
00389   
00390 //   if(mu != 1) LogTrace("MuonTrackAnalyzer") << "WARNING!! more than 1 simulated muon!!" <<endl;
00391 //   return result;
00392 
00393 
00394   // Loop over the Sim tracks
00395   SimTrackContainer::const_iterator simTrack;
00396   
00397   SimTrack result;
00398 
00399   double bestDeltaR = 10e5;
00400   for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack){
00401     if (abs((*simTrack).type()) != 13) continue;
00402     
00403     //    double newDeltaR = tsos.globalMomentum().basicVector().deltaR(simTrack->momentum().vect());
00404     GlobalVector tsosVect = tsos.globalMomentum();
00405     math::XYZVectorD vect(tsosVect.x(), tsosVect.y(), tsosVect.z());
00406     double newDeltaR = deltaR<double>(vect.eta(),vect.phi(),
00407                                       simTrack->momentum().eta(),simTrack->momentum().phi());
00408 
00409     if (  newDeltaR < bestDeltaR ) {
00410       LogTrace("MuonTrackAnalyzer") << "Matching Track with DeltaR = " << newDeltaR<<endl;
00411       bestDeltaR = newDeltaR;
00412       result  = *simTrack;
00413     }
00414   } 
00415   return pair<SimTrack,double>(result,bestDeltaR);
00416 }

bool MuonTrackAnalyzer::isInTheAcceptance ( double eta ) [private]

Definition at line 419 of file MuonTrackAnalyzer.cc.

References funct::abs(), all, barrel, endcap, lat::endl(), LogTrace, and theEtaRange.

Referenced by fillPlots().

00419                                                    {
00420   switch(theEtaRange){
00421   case all:
00422     return ( abs(eta) <= 2.4 ) ? true : false;
00423   case barrel:
00424     return ( abs(eta) < 1.1 ) ? true : false;
00425   case endcap:
00426     return ( abs(eta) >= 1.1 && abs(eta) <= 2.4 ) ? true : false;  
00427   default:
00428     {LogTrace("MuonTrackAnalyzer")<<"No correct Eta range selected!! "<<endl; return false;}
00429   }
00430 }

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

Definition at line 211 of file MuonTrackAnalyzer.cc.

References debug, lat::endl(), fillPlots(), getSeedTSOS(), getSimTrack(), hRecoSeedInner, hRecoSeedPCA, LogTrace, MuonUpdatorAtVertex::propagateToNominalLine(), theSeedsLabel, and theUpdator.

Referenced by analyze().

00212                                                                           {
00213 
00214   MuonPatternRecoDumper debug;
00215 
00216   // Get the RecTrack collection from the event
00217   Handle<TrajectorySeedCollection> seeds;
00218   event.getByLabel(theSeedsLabel, seeds);
00219   
00220   LogTrace("MuonTrackAnalyzer")<<"Number of reconstructed seeds: " << seeds->size()<<endl;
00221 
00222   for(TrajectorySeedCollection::const_iterator seed = seeds->begin(); 
00223       seed != seeds->end(); ++seed){
00224     TrajectoryStateOnSurface seedTSOS = getSeedTSOS(*seed);
00225     pair<SimTrack,double> sim = getSimTrack(seedTSOS,simTracks);
00226     fillPlots(seedTSOS, sim.first,
00227               hRecoSeedInner, debug);
00228     
00229     std::pair<bool,FreeTrajectoryState> propSeed =
00230       theUpdator->propagateToNominalLine(seedTSOS);
00231     if(propSeed.first)
00232       fillPlots(propSeed.second, sim.first,
00233                 hRecoSeedPCA, debug);
00234     else
00235       LogTrace("MuonTrackAnalyzer")<<"Error in seed propagation"<<endl;
00236    
00237   }
00238 }

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

Definition at line 241 of file MuonTrackAnalyzer.cc.

References TrajectoryStateOnSurface::charge(), debug, lat::endl(), MonitorElement::Fill(), fillPlots(), getSimTrack(), TrajectoryStateOnSurface::globalMomentum(), hChargeVsEta, hChargeVsPt, hDeltaPt_In_Out_VsEta, hDeltaPtVsEta, hNumberOfTracks, hNumberOfTracksVsEta, hPtRecVsPtGen, hRecoTracksInner, hRecoTracksOuter, hRecoTracksPCA, LogTrace, CoreSimTrack::momentum(), numberOfRecTracks, PV3DBase< T, PVType, FrameType >::perp(), funct::sqrt(), t, theService, theTracksLabel, track, and tracks.

Referenced by analyze().

00242                                                                           {
00243   MuonPatternRecoDumper debug;
00244   
00245   
00246   // Get the RecTrack collection from the event
00247   Handle<reco::TrackCollection> tracks;
00248   event.getByLabel(theTracksLabel, tracks);
00249 
00250   LogTrace("MuonTrackAnalyzer")<<"Reconstructed tracks: " << tracks->size() << endl;
00251   hNumberOfTracks->Fill(tracks->size());
00252   
00253   if(tracks->size()) numberOfRecTracks++;
00254   
00255   // Loop over the Rec tracks
00256   for(reco::TrackCollection::const_iterator t = tracks->begin(); t != tracks->end(); ++t) {
00257     
00258     reco::TransientTrack track(*t,&*theService->magneticField(),theService->trackingGeometry()); 
00259 
00260     TrajectoryStateOnSurface outerTSOS = track.outermostMeasurementState();
00261     TrajectoryStateOnSurface innerTSOS = track.innermostMeasurementState();
00262     TrajectoryStateOnSurface pcaTSOS   = track.impactPointState();
00263 
00264     pair<SimTrack,double> sim = getSimTrack(pcaTSOS,simTracks);
00265     SimTrack simTrack = sim.first;
00266     hNumberOfTracksVsEta->Fill(simTrack.momentum().eta(), tracks->size());
00267     fillPlots(track,simTrack);
00268     
00269     LogTrace("MuonTrackAnalyzer") << "State at the outer surface: " << endl; 
00270     fillPlots(outerTSOS,simTrack,hRecoTracksOuter,debug);
00271 
00272     LogTrace("MuonTrackAnalyzer") << "State at the inner surface: " << endl; 
00273     fillPlots(innerTSOS,simTrack,hRecoTracksInner,debug);
00274 
00275     LogTrace("MuonTrackAnalyzer") << "State at PCA: " << endl; 
00276     fillPlots(pcaTSOS,simTrack,hRecoTracksPCA,debug);
00277     
00278     double deltaPt_in_out = innerTSOS.globalMomentum().perp()-outerTSOS.globalMomentum().perp();
00279     hDeltaPt_In_Out_VsEta->Fill(simTrack.momentum().eta(),deltaPt_in_out);
00280 
00281     double deltaPt_pca_sim = pcaTSOS.globalMomentum().perp()-sqrt(simTrack.momentum().Perp2());
00282     hDeltaPtVsEta->Fill(simTrack.momentum().eta(),deltaPt_pca_sim);
00283     
00284     hChargeVsEta->Fill(simTrack.momentum().eta(),pcaTSOS.charge());
00285     
00286     hChargeVsPt->Fill(sqrt(simTrack.momentum().perp2()),pcaTSOS.charge());
00287     
00288     hPtRecVsPtGen->Fill(sqrt(simTrack.momentum().perp2()),pcaTSOS.globalMomentum().perp());    
00289   }
00290   LogTrace("MuonTrackAnalyzer")<<"--------------------------------------------"<<endl;  
00291 }

Member Data Documentation

DQMStore* MuonTrackAnalyzer::dbe_ [private]

Definition at line 83 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), endJob(), and MuonTrackAnalyzer().

std::string MuonTrackAnalyzer::dirName_ [private]

Definition at line 84 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and MuonTrackAnalyzer().

bool MuonTrackAnalyzer::doSeedsAnalysis [private]

Definition at line 98 of file MuonTrackAnalyzer.h.

Referenced by analyze(), beginJob(), endJob(), and MuonTrackAnalyzer().

bool MuonTrackAnalyzer::doTracksAnalysis [private]

Definition at line 97 of file MuonTrackAnalyzer.h.

Referenced by analyze(), beginJob(), endJob(), and MuonTrackAnalyzer().

MonitorElement* MuonTrackAnalyzer::hChargeVsEta [private]

Definition at line 113 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

MonitorElement* MuonTrackAnalyzer::hChargeVsPt [private]

Definition at line 114 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

MonitorElement* MuonTrackAnalyzer::hChi2 [private]

Definition at line 105 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hChi2Norm [private]

Definition at line 106 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hChi2NormVsEta [private]

Definition at line 118 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hChi2Prob [private]

Definition at line 109 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hChi2ProbVsEta [private]

Definition at line 121 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hChi2VsEta [private]

Definition at line 117 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hDeltaPt_In_Out_VsEta [private]

Definition at line 123 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

MonitorElement* MuonTrackAnalyzer::hDeltaPtVsEta [private]

Definition at line 122 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

MonitorElement* MuonTrackAnalyzer::hDof [private]

Definition at line 108 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hDofVsEta [private]

Definition at line 120 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hHitsPerTrack [private]

Definition at line 107 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hHitsPerTrackVsEta [private]

Definition at line 119 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and fillPlots().

MonitorElement* MuonTrackAnalyzer::hNumberOfTracks [private]

Definition at line 111 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

MonitorElement* MuonTrackAnalyzer::hNumberOfTracksVsEta [private]

Definition at line 112 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

MonitorElement* MuonTrackAnalyzer::hPtRecVsPtGen [private]

Definition at line 115 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

HTrack* MuonTrackAnalyzer::hRecoSeedInner [private]

Definition at line 127 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), endJob(), and seedsAnalysis().

HTrack* MuonTrackAnalyzer::hRecoSeedPCA [private]

Definition at line 128 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and seedsAnalysis().

HTrack* MuonTrackAnalyzer::hRecoTracksInner [private]

Definition at line 130 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

HTrack* MuonTrackAnalyzer::hRecoTracksOuter [private]

Definition at line 131 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), and tracksAnalysis().

HTrack* MuonTrackAnalyzer::hRecoTracksPCA [private]

Definition at line 129 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), endJob(), and tracksAnalysis().

HTrackVariables* MuonTrackAnalyzer::hSimTracks [private]

Definition at line 125 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), endJob(), and fillPlots().

int MuonTrackAnalyzer::numberOfRecTracks [private]

Definition at line 135 of file MuonTrackAnalyzer.h.

Referenced by endJob(), MuonTrackAnalyzer(), and tracksAnalysis().

int MuonTrackAnalyzer::numberOfSimTracks [private]

Definition at line 134 of file MuonTrackAnalyzer.h.

Referenced by endJob(), fillPlots(), and MuonTrackAnalyzer().

std::string MuonTrackAnalyzer::out [private]

Definition at line 86 of file MuonTrackAnalyzer.h.

Referenced by endJob(), and MuonTrackAnalyzer().

edm::InputTag MuonTrackAnalyzer::theCSCSimHitLabel [private]

Definition at line 93 of file MuonTrackAnalyzer.h.

Referenced by checkMuonSimHitPresence(), and MuonTrackAnalyzer().

edm::InputTag MuonTrackAnalyzer::theDTSimHitLabel [private]

Definition at line 94 of file MuonTrackAnalyzer.h.

Referenced by checkMuonSimHitPresence(), and MuonTrackAnalyzer().

EtaRange MuonTrackAnalyzer::theEtaRange [private]

Definition at line 89 of file MuonTrackAnalyzer.h.

Referenced by isInTheAcceptance(), and MuonTrackAnalyzer().

edm::InputTag MuonTrackAnalyzer::theRPCSimHitLabel [private]

Definition at line 95 of file MuonTrackAnalyzer.h.

Referenced by checkMuonSimHitPresence(), and MuonTrackAnalyzer().

std::string MuonTrackAnalyzer::theSeedPropagatorName [private]

Definition at line 99 of file MuonTrackAnalyzer.h.

Referenced by getSeedTSOS(), and MuonTrackAnalyzer().

edm::InputTag MuonTrackAnalyzer::theSeedsLabel [private]

Definition at line 92 of file MuonTrackAnalyzer.h.

Referenced by MuonTrackAnalyzer(), and seedsAnalysis().

MuonServiceProxy* MuonTrackAnalyzer::theService [private]

Definition at line 101 of file MuonTrackAnalyzer.h.

Referenced by analyze(), getSeedTSOS(), MuonTrackAnalyzer(), tracksAnalysis(), and ~MuonTrackAnalyzer().

edm::InputTag MuonTrackAnalyzer::theTracksLabel [private]

Definition at line 91 of file MuonTrackAnalyzer.h.

Referenced by beginJob(), MuonTrackAnalyzer(), and tracksAnalysis().

MuonUpdatorAtVertex* MuonTrackAnalyzer::theUpdator [private]

Definition at line 102 of file MuonTrackAnalyzer.h.

Referenced by MuonTrackAnalyzer(), and seedsAnalysis().

The documentation for this class was generated from the following files:

Validation/RecoMuon/src/MuonTrackAnalyzer.h
Validation/RecoMuon/src/MuonTrackAnalyzer.cc

Generated on Tue Jun 9 18:28:52 2009 for CMSSW by

1.5.4


Public Types
enum	EtaRange { all, barrel, endcap }
Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup)
virtual void	beginJob (const edm::EventSetup &eventSetup)
virtual void	endJob ()
	MuonTrackAnalyzer (const edm::ParameterSet &pset)
	Constructor.
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.
Private Member Functions
bool	checkMuonSimHitPresence (const edm::Event &event, edm::Handle< edm::SimTrackContainer > simTracks)
void	fillPlots (FreeTrajectoryState &recoFTS, SimTrack &simTrack, HTrack *histo, MuonPatternRecoDumper &debug)
void	fillPlots (TrajectoryStateOnSurface &recoTSOS, SimTrack &simState, HTrack *, MuonPatternRecoDumper &)
void	fillPlots (reco::TransientTrack &track, SimTrack &simTrack)
void	fillPlots (const edm::Event &event, edm::Handle< edm::SimTrackContainer > &simTracks)
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::InputTag	theDTSimHitLabel
EtaRange	theEtaRange
edm::InputTag	theRPCSimHitLabel
std::string	theSeedPropagatorName
edm::InputTag	theSeedsLabel
MuonServiceProxy *	theService
edm::InputTag	theTracksLabel
MuonUpdatorAtVertex *	theUpdator

MuonTrackAnalyzer Class Reference

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation