TopSanityCheckerModules

From T2B Wiki
Revision as of 12:29, 26 August 2015 by Maintenance script (talk | contribs) (Created page with " === TtGenEventChecker === === jetmet === === muonchecker === === vertex === === kinematics === === TruthReco === This module is executed after reconstructing ...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

TtGenEventChecker

jetmet

muonchecker

vertex

kinematics

TruthReco

This module is executed after reconstructing the ttbar generated event and after selecting the semi-leptonic ttbar decay into a muon (on generator level).

    • When matching is done, no jet is used twice and no parton is used twice!!!
    • Input is the selectedLayer1Objects and the generated ttbar event, genEvt.
    • For the calculation of the expected light and b JEC in bins of pT and eta, the default values of the ptbins and etabins can be found in python/TruthRecoChecker_cfi.py.
    • at least 1 muon, pT > 20, |eta| < 2.5 and isolated: vetoEM<4, vetoHad<6
    • at least 4 jets, pT > 30, |eta| < 2.5
subdir name of plot description/method remarks
hExpLightJECIncl (quarkEnergy-jetEnergy)/jetEnergy when the 4 quarks are matched!
hExpBJECIncl (quarkEnergy-jetEnergy)/jetEnergy when the 4 quarks are matched!
LowestPtMatchedJet_Rank jet with the lowest pT, matched to parton -> it's rank when the 4 quarks are matched!
LowestPtMatchedJet_Pt jet with the lowest pT, matched to parton -> it's pT when the 4 quarks are matched!
mWHad_Gen sum p4 of 2 light partons -> mass() when the 4 quarks are matched!
mWLep_Gen sum p4 of muon and neutrino -> mass() when the 4 quarks are matched!
mtHad_Gen sum p4 of the 3 quarks from top -> mass() when the 4 quarks are matched!
mtLep_Gen sum p4 of muon, neutrino and parton -> mass() when the 4 quarks are matched!
mWHad_Rec sum p4 of 2 light jets -> mass() when the 4 quarks are matched!
mWLep_Rec sum p4 of highest pT muon and MET -> mass() when the 4 quarks are matched!
mtHad_Rec sum p4 of the 3 jets from top -> mass() when the 4 quarks are matched!
mtLep_Rec sum p4 of highest pT muon, MET and b-jet -> mass() when the 4 quarks are matched!
Dphi_HadTops ROOT::Math::VectorUtil::DeltaPhi(tHadGen,tHadRec) when the 4 quarks are matched!
Dphi_LepTops ROOT::Math::VectorUtil::DeltaPhi(tLepGen,tLepRec) when the 4 quarks are matched!
Deta_HadTops fabs(tHadGen.eta() - tHadRec.eta()) when the 4 quarks are matched!
Deta_LepTops fabs(tLepGen.eta() - tLepRec.eta()) when the 4 quarks are matched!
DR_HadTops ROOT::Math::VectorUtil::DeltaR(tHadGen,tHadRec) when the 4 quarks are matched!
DR_LepTops ROOT::Math::VectorUtil::DeltaR(tHadGen,tHadRec) when the 4 quarks are matched!
Dphi_Tops_Rec ROOT::Math::VectorUtil::DeltaPhi(tHadRec,tLepRec) when the 4 quarks are matched!
Dphi_Tops_Gen ROOT::Math::VectorUtil::DeltaPhi(tHadGen,tLepGen) when the 4 quarks are matched!
DR_Tops_Rec ROOT::Math::VectorUtil::DeltaR(tHadRec,tLepRec) when the 4 quarks are matched!
DR_Tops_Gen ROOT::Math::VectorUtil::DeltaR(tHadGen,tLepGen) when the 4 quarks are matched!
minJetMuVSminQuarkMu smallest angle between jet and muon versus the smallest angle between the quarks and generated muon when the 4 quarks are matched!
minJetsVSminQuarks smallest angle between the 4 jets versus the smallest angle between the 4 quarks when the 4 quarks are matched!
QuarksMatchedJetRadiationMultiplicity # jets matched to radiated partons when the 4 quarks are matched!
QuarksMatchedHighestPtRadMatchedJet_Rank among the 4 highest pT jets, if at least 1 is matched with radiation, take the one with the highest pT -> rank when the 4 quarks are matched!
QuarksMatchedHighestPtRadMatchedJet_Pt among the 4 highest pT jets, if at least 1 is matched with radiation, take the one with the highest pT -> pT when the 4 quarks are matched!
LowestPtQuarkMatchedJetPt_vs_HighestPtRadMatchedJetPt 4th plot in list versus plot above when the 4 quarks are matched!
JetRadiationMultiplicity # jets matched to radiated partons when at least one of the 4 quarks is not matched!
HighestPtRadMatchedJet_Rank among the 4 highest pT jets, if at least 1 is matched with radiation, take the one with the highest pT -> rank when at least one of the 4 quarks is not matched!
HighestPtRadMatchedJet_Pt among the 4 highest pT jets, if at least 1 is matched with radiation, take the one with the highest pT -> pT when at least one of the 4 quarks is not matched!
UnmatchedQuarkMultiplicity # unmatched quarks among the 4 primary when at least one of the 4 quarks is not matched!
ExpectedJECinbins ExpectedLightJEC_etabinX_ptbinY same as inclusive but in bins of pT and eta fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
ExpectedJECinbins ExpectedBJEC_etabinX_ptbinY same as inclusive but in bins of pT and eta fitted with gaussian over 1/3 of the X-axis around the bin with the highest content


Resolutions_object

This module is executed after reconstructing the ttbar generated event.

    • object can be "lJets", "bJets", "muon", "electron", "met" (input is the selectedLayer1Objects and the generated ttbar event, genEvt)
    • resolutions are calculated as a function of pT and eta (except for MET, only as a function of pT, so by default only 1 etabin)
    • for the default values of the ptbins and etabins, look at python/ResolutionChecker_cfi.py
    • electrons and muons are required to match with the generated particle within DR 0.1
    • jets are required to match with the generated particle within DR 0.3
    • jets are required to match with the generated particle within DR 1000. (basically no requirement)
    • ares, bres, cres and dres: resolution on "a", "b", "c" & "d" with "a", "b", "c" & "d" according to the definition (CMS-NOTE-2006-023)
    • thres, phres, etres and etares: resolution on theta, phi, ET and eta
    • NOTE: more then one object can be matched to the same particle (e.g. 2 jets can be matched to the same parton)
name of plot description/method remarks
object_ares_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_bres_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_cres_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_dres_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_thres_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_phres_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_etres_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_etares_etabinX_ptbinY reconstructed value - generated value fitted with gaussian over 1/3 of the X-axis around the bin with the highest content
object_ares_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "[0]+[1]*exp(-[2]*x)", "x" is pT-dependence
object_bres_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "[0]+[1]*exp(-[2]*x)", "x" is pT-dependence
object_cres_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "[0]+[1]*exp(-[2]*x)", "x" is pT -dependence
object_dres_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "[0]+[1]*exp(-[2]*x)", "x" is pT-dependence
object_thres_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "[0]+[1]*exp(-[2]*x)", "x" is pT-dependence
object_phres_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "[0]+[1]*exp(-[2]*x)", "x" is pT-dependence
object_etres_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "pol1", "x" is pT-dependence
object_etares_etabinX for each pt-bin, mean and width of gaussian fit is used fitted with "[0]+[1]*exp(-[2]*x)", "x" is pT-dependence



Template:TracNotice