TrigAnalysisTest::Run1BStoxAODTrigger Class Reference

#include <Run1BStoxAODTrigger.h>

Inheritance diagram for TrigAnalysisTest::Run1BStoxAODTrigger:

Collaboration diagram for TrigAnalysisTest::Run1BStoxAODTrigger:

Public Member Functions
	Run1BStoxAODTrigger ()
	Run1BStoxAODTrigger (const Run1BStoxAODTrigger &)=delete
Run1BStoxAODTrigger &	operator= (const Run1BStoxAODTrigger &)=delete
void	processEvent (Trig::TrigDecisionTool &trigDecTool)
int	finalize ()
void	setEventStore (EventPtr_t store)
EventPtr_t	eventStore ()

Private Attributes
unsigned int	m_passed_l1
unsigned int	m_passed_hlt
std::map< std::string, int >	m_trigger_counts
TFile *	m_histFile
TH1 *	m_all_triggers
TH1 *	m_HLT_EFBMuMuXFex_eta
TH1 *	m_HLT_EFBMuMuXFex_mass
TH1 *	m_HLT_TrigTauJet
TH1 *	m_HLT_MuonEF_NSegments
TH1 *	m_HLT_MuonEF_pt
TH1 *	m_HLT_TrigEFMissingET
TH1 *	m_HLT_TrigEFMissingET_FEB
TH1 *	m_HLT_TrigEFMissingET_topocl
EventPtr_t	m_eventStore = nullptr

Detailed Description

Definition at line 20 of file Run1BStoxAODTrigger.h.

Constructor & Destructor Documentation

Run1BStoxAODTrigger() [1/2]

TrigAnalysisTest::Run1BStoxAODTrigger::Run1BStoxAODTrigger

(

)

Definition at line 40 of file Run1BStoxAODTrigger.cxx.

    : m_passed_l1(0),
      m_passed_hlt(0)
      //m_first_call(true)
  {
    m_histFile = new TFile("hist.root", "RECREATE");
    if (m_histFile == nullptr || !m_histFile->IsOpen()) {
      throw std::runtime_error("unable to open file for writing, exiting");
    }
 
    // bphys
    m_HLT_EFBMuMuXFex_eta = new TH1F( "HLT_EFBMuMuXFex_eta", "HLT_EFBMuMuXFex;#eta;No. of TrigEFBphys objects", 100, -4., 4. );
    m_HLT_EFBMuMuXFex_mass = new TH1F( "HLT_EFBMuMuXFex_mass", "HLT_EFBMuMuXFex;#M (GeV);No. of TrigEFBphys objects", 100, 0., 20. );
    // taus
    m_HLT_TrigTauJet = new TH1F( "HLT_TrigTauJet", "HLT_TrigTauJet;p_{T} (GeV);Jets", 100, 0., 500. );
    // muons
    m_HLT_MuonEF_NSegments = new TH1F( "HLT_MuonEF_NSegments", "HLT_MuonEF_NSegments;NSegments;Muons", 100, 0., 100.);
    m_HLT_MuonEF_pt = new TH1F( "HLT_MuonEF_pt", "HLT_MuonEF_pt;p_{T} (GeV);Muons", 100, 0., 100.);
    // MET
    m_HLT_TrigEFMissingET = new TH1F( "HLT_TrigEFMissingET", "HLT_TrigEFMissingET;E_{T} (GeV);Entries", 100, 0., 200.);
    m_HLT_TrigEFMissingET_FEB = new TH1F( "HLT_TrigEFMissingET_FEB", "HLT_TrigEFMissingET_FEB;E_{T} (GeV);Entries", 100, 0., 200. );
    m_HLT_TrigEFMissingET_topocl = new TH1F( "HLT_TrigEFMissingET_topocl", "HLT_TrigEFMissingET_topocl;E_{T} (GeV);Entries", 100, 0., 200. );
    // all trigger bits. The binning is changed in finalize() when the total number of triggers is known
    m_all_triggers = new TH1F("all_triggers", ";Trigger;No. of events", 1, 0., 1.);
  }

Run1BStoxAODTrigger() [2/2]

TrigAnalysisTest::Run1BStoxAODTrigger::Run1BStoxAODTrigger ( const Run1BStoxAODTrigger &  )

delete

Member Function Documentation

eventStore()

EventPtr_t TrigAnalysisTest::ITest::eventStore ( )

inlineinherited

Definition at line 29 of file ITest.h.

{ return m_eventStore; }

finalize()

int TrigAnalysisTest::Run1BStoxAODTrigger::finalize ( )

virtual

Implements TrigAnalysisTest::ITest.

Definition at line 166 of file Run1BStoxAODTrigger.cxx.

  {
    //athena doesn't store m_all_triggers without this
    m_histFile->cd();
 
    // Dump for debugging.
    cout << "L1 triggers passed: " << m_passed_l1 << endl;
    cout << "EF triggers passed: " << m_passed_hlt << endl;
    cout << "START SUMMARY: Run1BStoxAODTrigger" << endl;
    cout << "creating histogram with one bin per trigger, n of bins is " << m_trigger_counts.size() << endl;
    // Set the correct binning
    m_all_triggers->SetBins( m_trigger_counts.size(), 0.0, m_trigger_counts.size() );
    cout << "histogram created" << endl;
    int trigger_counter = 0;
    for (auto it = m_trigger_counts.begin(); it!=m_trigger_counts.end(); it++) {
      cout << (*it).first << "\t" << (*it).second <<endl;
      trigger_counter++;
      m_all_triggers->GetXaxis()->SetBinLabel( trigger_counter, ((*it).first).c_str() );
      m_all_triggers->SetBinContent( trigger_counter, (*it).second );
    }
    cout << "END SUMMARY  : Run1BStoxAODTrigger" << endl;
 
    m_histFile->Write();
    m_histFile->Close();
 
    // our aim here is to test the conversion of variables
    return 0;
  }

operator=()

Run1BStoxAODTrigger& TrigAnalysisTest::Run1BStoxAODTrigger::operator= ( const Run1BStoxAODTrigger &  )

delete

processEvent()

void TrigAnalysisTest::Run1BStoxAODTrigger::processEvent ( Trig::TrigDecisionTool &  trigDecTool )

virtual

Implements TrigAnalysisTest::ITest.

Definition at line 67 of file Run1BStoxAODTrigger.cxx.

                                                                       {
 
    // Did we pass a trigger level?
    if (trigDecTool.isPassed("L1_.*"))
      m_passed_l1++;
    if (trigDecTool.isPassed("EF_.*"))
      m_passed_hlt++;
 
    auto chainGroups = trigDecTool.getChainGroup(".*");
    for(auto &trig : chainGroups->getListOfTriggers()) {
      if(!m_trigger_counts.count(trig))
    m_trigger_counts[trig] = 0;
      if(trigDecTool.isPassed(trig)) {
    m_trigger_counts[trig] += 1;
    cout << "Trigger passed: " << trig << endl;
      }
    }
 
    //
    //extract trigger containers and fill histograms
    //
 
    //    bphys
    // HLT_EFBMuMuXFex
    const xAOD::TrigBphysContainer* HLT_xAOD__TrigBphysContainer_EFBMuMuXFex = 0;
    eventStore()->retrieve( HLT_xAOD__TrigBphysContainer_EFBMuMuXFex, "HLT_xAOD__TrigBphysContainer_EFBMuMuXFex" ).ignore();
    cout << "size of HLT_xAOD__TrigBphysContainer_EFBMuMuXFex is " << HLT_xAOD__TrigBphysContainer_EFBMuMuXFex->size() << endl;
    xAOD::TrigBphysContainer::const_iterator TrigBphys_iter;
    for(TrigBphys_iter = HLT_xAOD__TrigBphysContainer_EFBMuMuXFex->begin(); TrigBphys_iter != HLT_xAOD__TrigBphysContainer_EFBMuMuXFex->end(); ++TrigBphys_iter) {
      cout << "eta, mass = "<< (*TrigBphys_iter)->eta() << ", "<< (*TrigBphys_iter)->mass() << endl;
      m_HLT_EFBMuMuXFex_eta->Fill( (*TrigBphys_iter)->eta() );
      m_HLT_EFBMuMuXFex_mass->Fill( (*TrigBphys_iter)->mass() );
    }
 
    //    taus
    // HLT_TrigTauJet
    const xAOD::JetContainer* HLT_xAOD__JetContainer_TrigTauJet = 0;
    eventStore()->retrieve( HLT_xAOD__JetContainer_TrigTauJet, "HLT_xAOD__JetContainer_TrigTauJet" ).ignore();//weird name
    cout << "HLT_xAOD__JetContainer_TrigTauJet successfully retrieved - size is " << HLT_xAOD__JetContainer_TrigTauJet->size() << " jets" << endl;
    xAOD::JetContainer::const_iterator Jet_iter;
    for(Jet_iter = HLT_xAOD__JetContainer_TrigTauJet->begin(); Jet_iter != HLT_xAOD__JetContainer_TrigTauJet->end(); ++Jet_iter) {
      cout <<  "tau jet eta, pt: " << (*Jet_iter)->eta() << ", " << (*Jet_iter)->pt() << endl;
      m_HLT_TrigTauJet->Fill( (*Jet_iter)->pt() / 1000. );
    }
 
 
    //    muons
    // HLT_MuonEFInfo
    const xAOD::MuonContainer* HLT_xAOD__MuonContainer_MuonEFInfo = 0;
    eventStore()->retrieve( HLT_xAOD__MuonContainer_MuonEFInfo, "HLT_xAOD__MuonContainer_MuonEFInfo").ignore();
    cout << "size of HLT_xAOD__MuonContainer_MuonEFInfo is " << HLT_xAOD__MuonContainer_MuonEFInfo->size() << endl;
    xAOD::MuonContainer::const_iterator Muon_iter;
    for(Muon_iter = HLT_xAOD__MuonContainer_MuonEFInfo->begin(); Muon_iter!= HLT_xAOD__MuonContainer_MuonEFInfo->end(); ++Muon_iter) {
      cout << "muon n segments, pt: " << (*Muon_iter)->nMuonSegments() << ", " << (*Muon_iter)->pt() << endl;
      m_HLT_MuonEF_NSegments->Fill( float((*Muon_iter)->nMuonSegments()) );
      m_HLT_MuonEF_pt->Fill( (*Muon_iter)->pt() / 1000. );
    }
 
 
    //    MET
    // HLT_xAOD__TrigMissingETContainer_TrigEFMissingET
    const xAOD::TrigMissingETContainer* HLT_xAOD__TrigMissingETContainer_TrigEFMissingET = 0;
    eventStore()->retrieve( HLT_xAOD__TrigMissingETContainer_TrigEFMissingET, "HLT_xAOD__TrigMissingETContainer_TrigEFMissingET" ).ignore();
    xAOD::TrigMissingETContainer::const_iterator TrigMissingET_iter;
    cout << "Size of TrigMissingETContainer is " << HLT_xAOD__TrigMissingETContainer_TrigEFMissingET->size() << endl;
    for(TrigMissingET_iter = HLT_xAOD__TrigMissingETContainer_TrigEFMissingET->begin(); TrigMissingET_iter != HLT_xAOD__TrigMissingETContainer_TrigEFMissingET->end(); ++TrigMissingET_iter) {
      Double_t ex = (*TrigMissingET_iter)->ex();
      Double_t ey = (*TrigMissingET_iter)->ey();
      Double_t et = TMath::Sqrt(ex*ex + ey*ey);
      cout << " missing et = " << et << endl;
      m_HLT_TrigEFMissingET->Fill( et / 1000. );
    }
    // HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_FEB
    const xAOD::TrigMissingETContainer* HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_FEB = 0;
    eventStore()->retrieve( HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_FEB, "HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_FEB" ).ignore();
    cout << "Size of TrigMissingETContainer is " << HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_FEB->size() << endl;
    for(TrigMissingET_iter = HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_FEB->begin(); TrigMissingET_iter != HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_FEB->end(); ++TrigMissingET_iter) {
      Double_t ex = (*TrigMissingET_iter)->ex();
      Double_t ey = (*TrigMissingET_iter)->ey();
      Double_t et = TMath::Sqrt(ex*ex + ey*ey);
      cout << " missing et = " << et << endl;
      m_HLT_TrigEFMissingET_FEB->Fill( et / 1000. );
    }
    // HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl
    const xAOD::TrigMissingETContainer* HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl = 0;
    eventStore()->retrieve( HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl, "HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl" ).ignore();
    cout << "Size of TrigMissingETContainer is " << HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl->size() << endl;
    for(TrigMissingET_iter = HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl->begin(); TrigMissingET_iter != HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl->end(); ++TrigMissingET_iter) {
      Double_t ex = (*TrigMissingET_iter)->ex();
      Double_t ey = (*TrigMissingET_iter)->ey();
      Double_t et = TMath::Sqrt(ex*ex + ey*ey);
      cout << " missing et = " << et << endl;
      m_HLT_TrigEFMissingET_topocl->Fill( et / 1000. );
    }
 
    cout << "    --- Event processing done ---   \n" << endl;
  }

setEventStore()

void TrigAnalysisTest::ITest::setEventStore ( EventPtr_t  store )

inlineinherited

Definition at line 28 of file ITest.h.

{ m_eventStore = store; }

Member Data Documentation

m_all_triggers

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_all_triggers

private

Definition at line 38 of file Run1BStoxAODTrigger.h.

m_eventStore

EventPtr_t TrigAnalysisTest::ITest::m_eventStore = nullptr

privateinherited

Definition at line 31 of file ITest.h.

m_histFile

TFile* TrigAnalysisTest::Run1BStoxAODTrigger::m_histFile

private

Definition at line 35 of file Run1BStoxAODTrigger.h.

m_HLT_EFBMuMuXFex_eta

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_EFBMuMuXFex_eta

private

Definition at line 40 of file Run1BStoxAODTrigger.h.

m_HLT_EFBMuMuXFex_mass

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_EFBMuMuXFex_mass

private

Definition at line 41 of file Run1BStoxAODTrigger.h.

m_HLT_MuonEF_NSegments

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_MuonEF_NSegments

private

Definition at line 45 of file Run1BStoxAODTrigger.h.

m_HLT_MuonEF_pt

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_MuonEF_pt

private

Definition at line 46 of file Run1BStoxAODTrigger.h.

m_HLT_TrigEFMissingET

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_TrigEFMissingET

private

Definition at line 48 of file Run1BStoxAODTrigger.h.

m_HLT_TrigEFMissingET_FEB

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_TrigEFMissingET_FEB

private

Definition at line 49 of file Run1BStoxAODTrigger.h.

m_HLT_TrigEFMissingET_topocl

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_TrigEFMissingET_topocl

private

Definition at line 50 of file Run1BStoxAODTrigger.h.

m_HLT_TrigTauJet

TH1* TrigAnalysisTest::Run1BStoxAODTrigger::m_HLT_TrigTauJet

private

Definition at line 43 of file Run1BStoxAODTrigger.h.

m_passed_hlt

unsigned int TrigAnalysisTest::Run1BStoxAODTrigger::m_passed_hlt

private

Definition at line 31 of file Run1BStoxAODTrigger.h.

m_passed_l1

unsigned int TrigAnalysisTest::Run1BStoxAODTrigger::m_passed_l1

private

Definition at line 30 of file Run1BStoxAODTrigger.h.

m_trigger_counts

std::map<std::string,int> TrigAnalysisTest::Run1BStoxAODTrigger::m_trigger_counts

private

Definition at line 33 of file Run1BStoxAODTrigger.h.

---

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

• Run1BStoxAODTrigger.h
• Run1BStoxAODTrigger.cxx