Single-top s-Channel Analysis Author: Soeren Stamm <mailto: stamm@physik.hu-berlin.de> Update: $Id: AtlSgTop_sChannelAnalysis.cxx,v 1.29 2017/07/28 17:20:35 kaphle Exp $ Copyright: 2013 (C) Soeren Stamm
void | BookHistogramsBDT() |
void | BookHistogramsMCTruthStudy() |
void | BookHistogramsPreSel(AtlHistogramTool* htool) |
void | BookHistogramsPretag(AtlHistogramTool* htool) |
void | BookHistogramsSel(AtlHistogramTool* htool) |
void | AtlSelector::BookJobInfoHistograms() |
void | AtlSelector::ChangeOutputFile() |
void | AtlSelector::DoBookkeeping(TFile* InputFile) |
virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
void | FillHistogramsBDT() |
void | FillHistogramsElectrons(AtlHistogramTool* htool, Double_t W) |
void | FillHistogramsLeptons(AtlHistogramTool* htool, Double_t W) |
void | FillHistogramsMCTruthStudy() |
void | FillHistogramsMuons(AtlHistogramTool* htool, Double_t W) |
void | FillHistogramsPreSel(AtlHistogramTool* htool, Double_t W) |
void | FillHistogramsPretag(AtlHistogramTool* htool, Double_t W) |
void | FillHistogramsSel(AtlHistogramTool* htool, Double_t W) |
virtual void | InitEvent() |
void | TObject::MakeZombie() |
void | AtlSelector::SetSumw2(TDirectory* dir) |
void | ComputeBDTVariables() |
void | DoTruthMatching() |
Bool_t | ReconstructNeutrino() |
void | ReconstructTopDecays() |
void | ReconstructWDecay() |
void | SetBDTVariables() |
Bool_t | fAdditionalJetVeto | Remove events with additional jets outside the given pt-eta range |
Bool_t | fApplyDileptonVeto | Apply dilepton veto? |
Bool_t | fApplyTtbar4JetVeto | Apply ttbar 4 jet veto? |
Bool_t | fDoReconstructTop | Reconstruct the neutrino, the W-boson and the top-quark in the old-fashioned way and compute the 16 LPSC BDT variables (default=off). !!! Note that the neutrino reconstruction might fail in a very few cases leading to a slight loss of events !!! |
Bool_t | AtlSelector::fDoTruthTree | Read and process TruthTree information. Caution: this is a very time-consuming operation (Run2 SgTop-D3PDs only; default=false) |
Bool_t | fLeadingJetIsBTagged | Require the leading jet to be b-tagged |
HepMCParticle* | fMCAntibQuark | |
HepMCParticle* | fMCLeptonMinus | |
HepMCParticle* | fMCLeptonPlus | |
HepMCParticle* | fMCLooseLepton | |
HepMCParticle* | fMCTightLepton | |
Bool_t | fMCTruthStudy | Flag for filling mc truth info for dilepton study |
HepMCParticle* | fMCbQuark | |
Float_t | fMET_max | Maximum missing Et (GeV) |
Float_t | fMET_min | Minimum missing Et (GeV) |
Float_t | fMtW_max | Transverse W mass maximum (GeV), final event selection |
Float_t | fMtW_maxPre | Transverse W mass maximum (GeV), event pre-selection |
Float_t | fMtW_min | Transverse W mass minimum (GeV), final event selection |
Float_t | fMtW_minPre | Transverse W mass minimum (GeV), event pre-selection |
Int_t | fNBTags_max | Maximum no. of b-tagged jets |
Int_t | fNBTags_min | Minimum no. of b-tagged jets |
Float_t | fPt_BJetTop_min | Minimium Pt of b-jet from top decay |
Float_t | fPt_Jet_max | Additional maximum Pt cut for all jets (GeV |
Float_t | fPt_Jet_min | Additional minimum Pt cut for all jets (GeV) |
Float_t | fPt_LeadingJet_max | Additional maximum Pt cut for the leading jet (GeV) |
Float_t | fPt_LeadingJet_min | Additional minimum Pt cut for the leading jet (GeV) |
AtlSgTop_sChannelAnalysis::ETtbarDecayMode | fTtbarDecayMode | Top Quark decay products (MC truth, ttbar only) |
Float_t | fVetoBTagWeight_min | Minimum weight for b-tag veto (ctrl region); use 0.5 in case of run-2 data if the veto should be used (for run-2 no tag weights are stored - only boolean numbers |
AtlBTag::ETagger | fVetoBTagger | B-tagger used for b-tag veto (ctrl region) |
static const Int_t | AtlSelector::fgNumJetMults | |
static const Int_t | AtlSelector::fgNumLepChannels | |
static AtlSelector::EJetMult | AtlSelector::k4To6Jets | |
static TSelector::EAbort | TSelector::kAbortFile | |
static TSelector::EAbort | TSelector::kAbortProcess | |
static AtlSelector::EJetMult | AtlSelector::kAllJets | |
static AtlSelector::EIOMode | AtlSelector::kApp | |
static TObject::(anonymous) | TObject::kBitMask | |
static TObject::EStatusBits | TObject::kCanDelete | |
static TObject::EStatusBits | TObject::kCannotPick | |
static TSelector::EAbort | TSelector::kContinue | |
static AtlSelector::EIOMode | AtlSelector::kCustom | |
static AtlSelector::EIOMode | AtlSelector::kD3PDCKM | |
static AtlSelector::EIOMode | AtlSelector::kD3PDSgTop | |
static AtlSelector::EIOMode | AtlSelector::kD3PDSgTop_v2 | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kElEl | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kElMu | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kElQQ | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kElTau | |
static AtlSelector::ELepChannel | AtlSelector::kElectron | |
static AtlSelector::EJetMult | AtlSelector::kFiveJet | |
static AtlSelector::EJetMult | AtlSelector::kFourJet | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kHadronic | |
static TObject::EStatusBits | TObject::kHasUUID | |
static TObject::EStatusBits | TObject::kInvalidObject | |
static TObject::(anonymous) | TObject::kIsOnHeap | |
static TObject::EStatusBits | TObject::kIsReferenced | |
static AtlSelector::ELepChannel | AtlSelector::kLepton | |
static AtlSelector::EIOMode | AtlSelector::kMemNtuple | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kMuMu | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kMuQQ | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kMuTau | |
static AtlSelector::ELepChannel | AtlSelector::kMuon | |
static TObject::EStatusBits | TObject::kMustCleanup | |
static TObject::EStatusBits | TObject::kNoContextMenu | |
static TObject::(anonymous) | TObject::kNotDeleted | |
static TObject::EStatusBits | TObject::kObjInCanvas | |
static AtlSelector::EJetMult | AtlSelector::kOneJet | |
static AtlSelector::EJetMult | AtlSelector::kOneOrMoreJets | |
static TObject::(anonymous) | TObject::kOverwrite | |
static TObject::(anonymous) | TObject::kSingleKey | |
static AtlSelector::EJetMult | AtlSelector::kSixJet | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kTauQQ | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kTauTau | |
static AtlSelector::EJetMult | AtlSelector::kThreeJet | |
static AtlSelector::EJetMult | AtlSelector::kThreeOrMoreJets | |
static AtlSelector::EJetMult | AtlSelector::kTwoJet | |
static AtlSelector::EJetMult | AtlSelector::kTwoOrMoreJets | |
static AtlSgTop_sChannelAnalysis::ETtbarDecayMode | kUnknown | |
static TObject::(anonymous) | TObject::kWriteDelete | |
static TObject::(anonymous) | TObject::kZombie |
TSelector::EAbort | TSelector::fAbort | Abort status |
Int_t | AtlSelector::fAcceptedEvents | No. of accepted events (w/o weights) selected by ProcessCut() |
Float_t | AtlSelector::fAcceptedEventsB | No. of accepted events (weighted) selected by ProcessCut() |
Float_t | AtlSelector::fAcceptedEventsW | No. of accepted events (weighted) selected by ProcessCut() |
TObjArray* | AtlSelector::fBookkeepingList | List of chain elements needed for bookkeeping |
Bool_t | AtlSelector::fCountUnfilteredEvents | Count the no. of all unfiltered events (in case the input A++ file does not contain a JobInfo histogram) |
TTree* | AtlSelector::fCurrentTree | Currently analysed tree |
Long64_t | AtlSelector::fEntry | Global entry number of the current event in input chain |
TEntryList* | AtlSelector::fEntryList | Global entry list for all files in a chain |
AtlEvent* | AtlSelector::fEvent | Point to event object in memory |
AtlEvtReaderBase* | AtlSelector::fEvtReader | Reader for input events |
AtlEvtReaderBase* | AtlSelector::fEvtReaderUser | User-provided reader for input events (owned!) |
AtlEvtWriterBase* | AtlSelector::fEvtWriter | Writer for output event |
TH1D* | AtlSelector::fHistCutflow_NoWeights | SgTopD3PD CutFlow Histogram |
TH1D* | AtlSelector::fHistCutflow_allEventsHFOR_NoWeights | SgTopD3PD CutFlow Histogram |
TH1D* | AtlSelector::fHistCutflow_allEventsHFOR_genWeights | SgTopD3PD CutFlow Histogram |
TH1D* | AtlSelector::fHistCutflow_genWeights | SgTopD3PD CutFlow Histogram |
TH1D* | AtlSelector::fHistCutflow_genXpileupWeights | SgTopD3PD CutFlow Histogram |
TH1D* | AtlSelector::fHistCutflow_genXpileupXzvertexWeights | SgTopD3PD CutFlow Histogram |
TH1F* | AtlSelector::fHistEvents | Histogram for no. of processed and accepted events etc. |
TH1F* | AtlSelector::fHistEvtWeights | Histogram of event weights |
TList* | TSelector::fInput | List of objects available during processing |
AtlSelector::EIOMode | AtlSelector::fInputMode | Defines input file mode |
Bool_t | AtlSelector::fIsChain | Input is chain or tree ? |
Int_t | AtlSelector::fJetBinEvents | No. of events in jet bin (w/o weights) |
Float_t | AtlSelector::fJetBinEventsW | No. of events in jet bin (weighted) |
TList* | AtlSelector::fListOfTools | List of analysis tools |
MemInfo_t* | AtlSelector::fMemInfo | Current memory usage information |
Int_t | AtlSelector::fNBookkeeping | No. of bookkeeping files |
TObject* | TSelector::fObject | ! Current object if processing object (vs. TTree) |
TString | TSelector::fOption | Option given to TTree::Process |
TSelectorList* | TSelector::fOutput | ! List of objects created during processing |
TFile* | AtlSelector::fOutputFile | Output file |
TString* | AtlSelector::fOutputFilename | Name of the output file |
AtlSelector::EIOMode | AtlSelector::fOutputMode | Defines mode for the output tree (if any) |
TTree* | AtlSelector::fOutputTree | A++ output tree (if any) |
TString* | AtlSelector::fOutputTreeName | Name of the A++ output tree (if any) |
TString* | AtlSelector::fOutputTreeTitle | Title of the A++ output tree (if any) |
TTree* | AtlSelector::fOutputTriggerConfTree | Needed for copying the trigger config dbases into the output file in case of written out A++ events |
Bool_t | AtlSelector::fPassedSelection | variable for storing the selection decisions for the event |
Int_t | AtlSelector::fPreAcceptedEvents | No. of pre-accepted events (w/o weights) selected by ProcessPreCut() |
Float_t | AtlSelector::fPreAcceptedEventsW | No. of pre-accepted events (weighted) selected by ProcessPreCut() |
Float_t | AtlSelector::fPreTaggedEvents | No. of events before b-tagging |
Float_t | AtlSelector::fPreTaggedEventsW | No. of events before b-tagging (weighted) |
Bool_t | AtlSelector::fPrintEvent | Print the content of every event (useful for testing/debugging) |
Bool_t | AtlSelector::fPrintObjectTable | Print Root's object table for every event (useful for tracking memory leaks) |
ProcInfo_t* | AtlSelector::fProcInfo | Current process information |
Int_t | AtlSelector::fProcessedEvents | No. of processed events (w/o weights) |
Float_t | AtlSelector::fProcessedEventsW | No. of processed events (weighted) |
Long64_t | TSelector::fStatus | Selector status |
TStopwatch | AtlSelector::fStopwatch | Stopwatch used for process monitoring |
TTree* | AtlSelector::fTree | Pointer to the currently analysed tree or chain |
Bool_t | AtlSelector::fWriteEntryList | Switch for writing an entry list |
TList* | fBJets | List of selected b-jets |
Float_t | fCosLepTopBJet1_TopBJet1CMS | Cos of angle between the lepton momentum in the top-quark rest-frame and the top-quark momentum in the c.m.s. The used top-quark is the candidate reconstructed using the leading b-jet. The c.m.s vector is the sum of the top momentum and the sub-leading b-jet momentum |
Float_t | fCosLepTopBJet2_TopBJet2CMS | Cos of angle between the lepton momentum in the top-quark rest-frame and the top-quark momentum in the c.m.s. The used top-quark is the candidate reconstructed using the sub-leading b-jet. The c.m.s vector is the sum of the top momentum and the leading b-jet momentum |
Float_t | fCosLepW_WTopBJet1 | Cos of angle between lepton momentum in W rest-frame and the W momentum in the rest-frame of the top-quark candidate with the leading b-jet |
Float_t | fCosMET_BJet2 | Cos of angle between MET and sub-leading b-jet |
AtlCutFlowTool* | fCutFlowTool | |
Float_t | fDeltaEtaLep_BJet1 | Delta eta between the lepton and the leading b-jet |
Float_t | fDeltaEtaLep_BJet2 | Delta eta between the lepton and the sub-leading b-jet |
Float_t | fDeltaEtaNu_BJetNonTop | Delta eta between the rec. neutrino and the b-jet from the non-top candidate |
Float_t | fDeltaEtaTopBJet2_BJet1 | Delta eta between top decay reconstructed using the sub-leading b-jet and the leading b-jet |
Float_t | fDeltaPhiBJet1_BJet2 | Delta phi between leading and sub-leading b-jet |
Float_t | fDeltaPhiLep_MET | Delta phi between lepton and MET |
Float_t | fDeltaPhiTopBJet1_BJet2 | Delta phi between top-quark candidate reconstructed using the leading b-jet and the sub-leading b-jet |
Float_t | fDeltaPhiTopBJet2_BJet1 | Delta phi between top-quark candidate reconstructed using the sub-leading b-jet and the leading b-jet |
TList* | fElectrons | List of selected electron objects |
TH1F* | fHist_BDT_CosLepTopBJet1_TopBJet1CMS | |
TH1F* | fHist_BDT_CosLepTopBJet2_TopBJet2CMS | |
TH1F* | fHist_BDT_CosLepW_WTopBJet1 | |
TH1F* | fHist_BDT_CosMET_BJet2 | |
TH1F* | fHist_BDT_DeltaEtaLep_BJet1 | |
TH1F* | fHist_BDT_DeltaEtaLep_BJet2 | |
TH1F* | fHist_BDT_DeltaEtaNu_BJetNonTop | |
TH1F* | fHist_BDT_DeltaEtaTopBJet2_BJet1 | |
TH1F* | fHist_BDT_DeltaPhiBJet1_BJet2 | |
TH1F* | fHist_BDT_DeltaPhiLep_MET | |
TH1F* | fHist_BDT_DeltaPhiTopBJet1_BJet2 | |
TH1F* | fHist_BDT_DeltaPhiTopBJet2_BJet1 | |
TH1F* | fHist_BDT_Lep_Pt | |
TH1F* | fHist_BDT_MET | |
TH1F* | fHist_BDT_MLep_BJet2 | |
TH1F* | fHist_BDT_MtW | |
TH1F* | fHist_BDT_PtBJet_Top | |
TH1F* | fHist_BDT_SumEtLep_MET | |
TH1F* | fHist_BDT_SumEtLep_MET_BJets | |
TH1F* | fHist_BDT_SumPt | |
TH1F* | fHist_BDT_SumPtBJet1_BJet2 | |
TH1F* | fHist_mc_DecayChannels | |
TH1F* | fHist_mc_DecayChannels_AllVeto | |
TH1F* | fHist_mc_DecayChannels_TauVeto | |
TH1F* | fHist_mc_DecayChannels_Veto | |
TH2F* | fHist_mc_DeltaR_2ndLeadLepton_bQuarkOS | |
TH2F* | fHist_mc_DeltaR_2ndLeadLepton_bQuarkSS | |
TH2F* | fHist_mc_DeltaR_LeadLepton_2ndLeadLepton | |
TH2F* | fHist_mc_DeltaR_LeadLepton_bQuarkOS | |
TH2F* | fHist_mc_DeltaR_LeadLepton_bQuarkSS | |
TH2F* | fHist_mc_DeltaR_loose_lepton | |
TH2F* | fHist_mc_DeltaR_tight_lepton | |
TH2F* | fHist_mc_DeltaR_unmatchedMCLepton_truthjet | |
TH2F* | fHist_mc_MET | |
TH2F* | fHist_mc_MET_nonVeto | |
TH2F* | fHist_mc_RecoEff | |
TH2F* | fHist_mc_TriggerMatch | |
TH1F* | fHist_mc_jet_truth_label | |
TH2F* | fHist_mc_reco_2ndLeadLepton_Pt | |
TH2F* | fHist_mc_reco_2ndLeadLepton_Pt_res | |
TH2F* | fHist_mc_reco_LeadLepton_Pt | |
TH2F* | fHist_mc_reco_LeadLepton_Pt_res | |
TH2F* | fHist_mc_truth_2ndLeadLepton_Pt | |
TH2F* | fHist_mc_truth_2ndLeadLepton_Pt_matched | |
TH2F* | fHist_mc_truth_2ndLeadLepton_Pt_unmatched | |
TH2F* | fHist_mc_truth_LeadLepton_Pt | |
TH2F* | fHist_mc_unmatchedMCLepton_truthjet_pt | |
AtlHistogramTool* | fHistsBTag | |
AtlHistogramTool* | fHistsBTagElectron | |
AtlHistogramTool* | fHistsBTagMuon | |
AtlHistogramTool* | fHistsPretag | |
TList* | fJets | List of selected jets |
HepParticle* | fLepton | The leading lepton |
TList* | fLeptons | List of selected leptons |
TList* | fLooseLeptons | List of selected loose leptons |
TList* | fLooseTaus | List of selected loose taus |
Float_t | fMET | MissingEt of current event |
Float_t | fMLep_BJet2 | Invariant mass of lepton and sub-leading b-jet (GeV) |
Float_t | fMtW | Transverse W mass (GeV) |
TList* | fMuons | List of selected muon objects |
AtlObjectsToolD3PDSgTop* | fObjTool | |
AtlPdfReweightingTool* | fPdfTool | |
Float_t | fPtLep | Transverse momentum (GeV) of leading lepton |
Float_t | fPt_BJetTop | Pt of b-jet from best recosntructed top decay (GeV) |
Bool_t | fQCDMMToolIsOn | Flag for QCD Matrix Method tasks |
AtlObjRecoScaleFactorTool* | fSFTool | |
Float_t | fSumEtLep_MET | Sum Et of lepton and MET |
Float_t | fSumEtLep_MET_BJets | Sum Et of lepton, MET and all b-jets |
Float_t | fSumPt | Vectoriell Pt sum of the lepton, MET and all b-jets |
Float_t | fSumPtBJet1_BJet2 | Vectoriell Pt sum of the leading and the sub-leading b-jet (GeV) |
Inheritance Chart: | |||||||||||||||||
|
Fill histograms for final event selection
Dilepton veto (against ttbar dilepton bkg) A return value of 'true' means the veto is fulfilled, ie. the event ought to be discarded. Note: for run-1 events, the veto uses loose lepton information, which is at present not available for run-2 events. Therefore here the tight leptons are used for the moment.
Fill histograms This method will be called only for events surviving the ProcessCut() routine
Rconstruct t->Wb decay candidates. ReconstructWDecay() and ReconstructNeutrino() need to be called before. The candidates for all b-jets are added to the list of top decays of the event. The list is sorted in descendent order (!) w.r.t. the proximity of the top pole-mass of 172.5 GeV.
- Assign true particles to data members - Determine decay channel of top pair based on MC truth info
Get best top decay candidate (invariant mass closest to pole-mass.