Difference: HbyHm (8 vs. 9)

Revision 92010-10-18 - CameliaMironov

Line: 1 to 1
 
META TOPICPARENT name="Works"

Hit-by-hit matching

Line: 137 to 137
 
 
  • matching decision based only on x_purity and fracHitMatch>0/0.5 and <0.75 are giving the right results ....
  • in HI, the fake rate will be higher-->
in case one wants to present it separately, should chose a working point (purity,fracHitMatch) for which this rate is small
  • (purity_X=75%, and quality>0.5) suggested working point
  • Added:
    >
    >
    mport SimMuon.MCTruth.MuonAssociatorByHits_cfi
    process.glbMuonAssociatorByHits = SimMuon.MCTruth.MuonAssociatorByHits_cfi.muonAssociatorByHits.clone()
    process.glbMuonAssociatorByHits.tracksTag                     = 'globalMuons'
    process.glbMuonAssociatorByHits.UseTracker                    = True
    process.glbMuonAssociatorByHits.UseMuon                       = True
    process.glbMuonAssociatorByHits.PurityCut_track             = 0.75
    process.glbMuonAssociatorByHits.EfficiencyCut_track      = 0.
    process.glbMuonAssociatorByHits.PurityCut_muon           = 0.75
    process.glbMuonAssociatorByHits.EfficiencyCut_muon    = 0.
    process.glbMuonAssociatorByHits.includeZeroHitMuons  = False
    process.glbMuonAssociatorByHits.acceptOneStubMatchings  = False
    

    and then, in the anlyzer, when looping over maps,

    if(sim2Reco.find(simTrack) != sim2Reco.end())     
       {
         vector<pair<edm::RefToBase<reco::Track>, double> > recTracks = sim2Reco[simTrack];
         if(recTracks.size() !=0 )
           {
             fFracShared    = recTracks.begin()->second;
             if(  fFracShared >=0.5 )  nMatches+=1;
          
           }
       }
    
    
     
    • Note: corrections wrt pT_gen: correct up to resolution effects (best seen in the sTA plots)

     
    This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
    Ideas, requests, problems regarding TWiki? Send feedback