WTAConeMaker Class Reference

#include <WTAConeMaker.h>

Inheritance diagram for WTAConeMaker:

Collaboration diagram for WTAConeMaker:

Public Member Functions
	WTAConeMaker (bool debug=false, bool verbose=false)
virtual	~WTAConeMaker ()=default
void	ResizeConstituentList ()
void	ResizeSeedSortingList ()
void	ResizeSeedList ()
void	ResizeSeedNConstLists ()
void	ResizeThisJetConstituents (WTAJet &jet)
void	ClearLists ()
void	PrintSeedList ()
WTAJet	WTATrigObjToWTAJet (const WTATrigObj &obj)
WTATrigObj	WTAJetToWTATrigObj (const WTAJet &jet)
const std::vector< WTATrigObj > &	GetConstituentList () const
const std::vector< WTATrigObj > &	GetSeedSortingList () const
const std::vector< WTAJet > &	GetSeedList () const
std::vector< int >	GetAssociateBits (WTATrigObj incoming_seed, int &max_pt_index)
std::vector< WTATrigObj >	LoadInputs (const std::vector< IntOrFloat > &ptVec, const std::vector< IntOrFloat > &etaVec, const std::vector< IntOrFloat > &phiVec, const std::vector< IntOrFloat > &mVec)
virtual void	FillLists (const std::vector< WTATrigObj > &InputTowers)
void	InitiateInputs (const std::vector< IntOrFloat > &ptVec, const std::vector< IntOrFloat > &etaVec, const std::vector< IntOrFloat > &phiVec, const std::vector< IntOrFloat > &mVec)
void	InitiateInputs (const std::vector< WTATrigObj > &InputTowers)
void	InsertToConstList (const WTATrigObj &obj)
virtual void	SeedCleaning ()
virtual void	MergeConstsToSeeds ()
void	SetDEBUG ()
void	SetVERBOSE ()
bool	GetDEBUG ()
bool	GetVERBOSE ()

Public Attributes
WTAParameters	m_WTAConeMakerParameter

Protected Attributes
std::vector< WTATrigObj >	m_ConstituentList
std::vector< WTATrigObj >	m_SeedSortingList
std::vector< WTAJet >	m_SeedList
bool	m_DEBUG
bool	m_VERBOSE

Detailed Description

Definition at line 54 of file WTAConeMaker.h.

Constructor & Destructor Documentation

WTAConeMaker()

WTAConeMaker::WTAConeMaker ( bool debug = false, bool verbose = false )

inline

Definition at line 56 of file WTAConeMaker.h.

                                                              :
            m_DEBUG(debug), m_VERBOSE(verbose)
            {}; // Constructor

~WTAConeMaker()

virtual WTAConeMaker::~WTAConeMaker ( )

virtualdefault

Member Function Documentation

ClearLists()

void WTAConeMaker::ClearLists ( )

inline

Definition at line 66 of file WTAConeMaker.h.

{m_ConstituentList.clear(); m_SeedSortingList.clear(); m_SeedList.clear();};

FillLists()

void WTAConeMaker::FillLists ( const std::vector< WTATrigObj > & InputTowers )

inlinevirtual

Reimplemented in WTACone2PassMaker.

Definition at line 114 of file WTAConeMaker.h.

{
    m_ConstituentList.clear(); m_SeedSortingList.clear();
    const unsigned int MaxSeedSortingN = m_WTAConeMakerParameter.GetMaxSeedSortingN();
    const unsigned int MaxConstN = m_WTAConeMakerParameter.GetMaxConstN();
    const int MaxTowersToReadPerEvent = m_WTAConeMakerParameter.GetMaxInputTowers();
    if(m_DEBUG) std::cout << "MaxTowersToReadPerEvent = " << MaxTowersToReadPerEvent << ", Size of event = " << InputTowers.size() << std::endl;
    int nTowers = 0;
    for(const auto& tower: InputTowers)
    {
        if(tower.pt() < m_WTAConeMakerParameter.GetConstEtCut())continue; // Skip Et < 2GeV
        if(tower.pt() >= m_WTAConeMakerParameter.GetSeedEtCut())m_SeedSortingList.push_back(tower); // Harmonize >=
        else m_ConstituentList.push_back(tower); // Initially, always put soft tower at the back of the m_ConstituentList
 
        // We can read only first N towers out of entier event
        nTowers++;
        if(nTowers>MaxTowersToReadPerEvent) break;
    }
    if(m_SeedSortingList.size() > MaxSeedSortingN)m_SeedSortingList.resize(MaxSeedSortingN);
    if(m_ConstituentList.size() > MaxConstN)m_ConstituentList.resize(MaxConstN); // Resize lists accordingly
}

GetAssociateBits()

std::vector< int > WTAConeMaker::GetAssociateBits ( WTATrigObj incoming_seed, int & max_pt_index )

inline

Definition at line 210 of file WTAConeMaker.h.

{
    int jet_N = m_SeedList.size();
    std::vector<int> associate_bit(jet_N, 0);
    IntOrFloat MaxPt = tower.pt(); MaxPtIndex = -1;
    for(int j = 0; j < jet_N; j++)
    {
        if(m_VERBOSE)std::cout << "deta, dphi, dR2 = " << tower.d_eta(m_SeedList.at(j)) << " , " << tower.d_phi_MPI_PI(m_SeedList.at(j)) << " , " << tower.dR2(m_SeedList.at(j)) << std::endl;
        if(tower.IsAssocdR(m_SeedList.at(j), m_WTAConeMakerParameter.GetIso_dR2())) // e.g) dR2 < 0.16, association
        {
            associate_bit.at(j) = 1;
            if(m_SeedList.at(j).pt() > MaxPt)
            {
                MaxPtIndex = j;
                MaxPt = m_SeedList.at(j).pt(); // Update the counter
            }
        }
    } // associate_bit filling done
    return associate_bit;
}

GetConstituentList()

const std::vector< WTATrigObj > & WTAConeMaker::GetConstituentList ( ) const

inline

Definition at line 70 of file WTAConeMaker.h.

{return m_ConstituentList;}; // Access these lists whenever we need

GetDEBUG()

bool WTAConeMaker::GetDEBUG ( )

inline

Definition at line 86 of file WTAConeMaker.h.

{return m_DEBUG;}; // Printout for debug

GetSeedList()

const std::vector< WTAJet > & WTAConeMaker::GetSeedList ( ) const

inline

Definition at line 72 of file WTAConeMaker.h.

{return m_SeedList;};

GetSeedSortingList()

const std::vector< WTATrigObj > & WTAConeMaker::GetSeedSortingList ( ) const

inline

Definition at line 71 of file WTAConeMaker.h.

{return m_SeedSortingList;};

GetVERBOSE()

bool WTAConeMaker::GetVERBOSE ( )

inline

Definition at line 87 of file WTAConeMaker.h.

{return m_VERBOSE;};

InitiateInputs() [1/2]

void WTAConeMaker::InitiateInputs ( const std::vector< IntOrFloat > & ptVec, const std::vector< IntOrFloat > & etaVec, const std::vector< IntOrFloat > & phiVec, const std::vector< IntOrFloat > & mVec )

inline

Definition at line 136 of file WTAConeMaker.h.

{
    std::vector<WTATrigObj> InputTowers = LoadInputs(ptVec, etaVec, phiVec, mVec);
    FillLists(InputTowers);
    if(m_VERBOSE)std::cout << "InputN, ConstN, SeedN = " << InputTowers.size() << " , " << m_ConstituentList.size() << " , " << m_SeedSortingList.size() << std::endl;
}

InitiateInputs() [2/2]

void WTAConeMaker::InitiateInputs ( const std::vector< WTATrigObj > & InputTowers )

inline

Definition at line 143 of file WTAConeMaker.h.

{
    FillLists(InputTowers);
    if(m_VERBOSE)std::cout << "InputN, ConstN, SeedN = " << InputTowers.size() << " , " << m_ConstituentList.size() << " , " << m_SeedSortingList.size() << std::endl;
}

InsertToConstList()

void WTAConeMaker::InsertToConstList ( const WTATrigObj & obj )

inline

Definition at line 161 of file WTAConeMaker.h.

{
    if(m_WTAConeMakerParameter.GetAddConstFirst())m_ConstituentList.insert(m_ConstituentList.begin(), obj); // Insert obj at the beginnning
    else m_ConstituentList.push_back(obj); // Insert obj at the end
}

LoadInputs()

std::vector< WTATrigObj > WTAConeMaker::LoadInputs ( const std::vector< IntOrFloat > & ptVec, const std::vector< IntOrFloat > & etaVec, const std::vector< IntOrFloat > & phiVec, const std::vector< IntOrFloat > & mVec )

inline

Definition at line 101 of file WTAConeMaker.h.

{
    std::vector<WTATrigObj> input_towers;
    unsigned int tower_n = ptVec.size();
    for(unsigned int t = 0; t < tower_n; t++)
    {
        WTATrigObj this_tower(ptVec.at(t), etaVec.at(t), phiVec.at(t), mVec.at(t));
        input_towers.push_back(this_tower);
    }
    return input_towers;
 
}

MergeConstsToSeeds()

void WTAConeMaker::MergeConstsToSeeds ( )

inlinevirtual

Reimplemented in WTACone2PassMaker.

Definition at line 282 of file WTAConeMaker.h.

{
    for(auto constituent: m_ConstituentList)
    {
        for(unsigned int j = 0; j < m_SeedList.size(); j++) // Assume Seeds are pT sorted, WTA means more energetic jet eats constituent first
        {
            IntOrFloat dR2 = constituent.dR2(m_SeedList.at(j));
            if(dR2 != 0 && constituent.IsAssocdR(m_SeedList.at(j), m_WTAConeMakerParameter.GetJet_dR2())) // Thistime, the condition is m_Jet_dR2, the usual R2Par, **WARNING: dR2!=0 IS TEMPORARY, NEED TO KNOW TOPOTOWER CREATION
            {
                m_SeedList.at(j).MergeConstituent(constituent);
                ResizeThisJetConstituents(m_SeedList.at(j)); // Check JetConstituent N
                break; // Break the jet loop, move to the next constituent
            }
        }
    }
}

PrintSeedList()

void WTAConeMaker::PrintSeedList ( )

inline

Definition at line 149 of file WTAConeMaker.h.

{
    int AllJetConstN = 0;
    for(const auto& jet: m_SeedList)
    {
        AllJetConstN += jet.GetConstituentCount();
        std::cout << "Jet pT, eta, phi, constN = " << jet.pt() << " , " << jet.eta() << " , " << jet.phi() << " , " << jet.GetConstituentCount() << std::endl;
    }
    std::cout << "PrintSeedList..." << " , SeedN, AllJetConstN, ConstN = " << m_SeedList.size() << " , " << AllJetConstN << " , " << m_ConstituentList.size() <<  std::endl;
    std::cout << "+++++++++++++++++++++++++" << std::endl;
}

ResizeConstituentList()

void WTAConeMaker::ResizeConstituentList ( )

inline

Definition at line 167 of file WTAConeMaker.h.

{
    const unsigned int MaxConstN = m_WTAConeMakerParameter.GetMaxConstN();
    while(m_ConstituentList.size() > MaxConstN)m_ConstituentList.pop_back();
}

ResizeSeedList()

void WTAConeMaker::ResizeSeedList ( )

inline

Definition at line 188 of file WTAConeMaker.h.

{
    const unsigned int MaxSeedN = m_WTAConeMakerParameter.GetMaxSeedN();
    while(m_SeedList.size() > MaxSeedN)
    {
        InsertToConstList(WTAJetToWTATrigObj(m_SeedList.back()));
        m_SeedList.pop_back();
    }
}

ResizeSeedNConstLists()

void WTAConeMaker::ResizeSeedNConstLists ( )

inline

Definition at line 64 of file WTAConeMaker.h.

{ResizeSeedList(); ResizeConstituentList();}

ResizeSeedSortingList()

void WTAConeMaker::ResizeSeedSortingList ( )

inline

Definition at line 173 of file WTAConeMaker.h.

{
    const unsigned int MaxSeedSortingN = m_WTAConeMakerParameter.GetMaxSeedSortingN();
    while(m_SeedSortingList.size() > MaxSeedSortingN)m_SeedSortingList.pop_back();
}

ResizeThisJetConstituents()

void WTAConeMaker::ResizeThisJetConstituents ( WTAJet & jet )

inline

Definition at line 179 of file WTAConeMaker.h.

{
    const unsigned int MaxConstPerJetN = m_WTAConeMakerParameter.GetMaxConstPerJetN();
    while(jet.GetConstituentCount() > MaxConstPerJetN)
    {
        jet.PopOutLastConstituent(); // Erase the last constituent, it is NOT appended to the Const list
    }
}

SeedCleaning()

void WTAConeMaker::SeedCleaning ( )

inlinevirtual

Reimplemented in WTACone2PassMaker.

Definition at line 231 of file WTAConeMaker.h.

{
    m_SeedList.clear();
    if(m_DEBUG)std::cout << "Baseline Seed Cleaning......" << std::endl;
    for(const auto& seed: m_SeedSortingList)
    {
        int jet_N = m_SeedList.size();
        if(jet_N == 0)
        {
            m_SeedList.push_back(WTATrigObjToWTAJet(seed)); // Take first seed as jet
        }
        else
        {   
            int MaxPtIndex = -1;
            std::vector<int> associate_bit = GetAssociateBits(seed, MaxPtIndex); // Associate_bit filling done
            if(std::find(associate_bit.begin(), associate_bit.end(), 1) != associate_bit.end())
            { // When there is at least one association
                if(MaxPtIndex == -1)
                { // Incoming Seed is the highest et, do the seed cleaning
                    WTAJet IncomingSeedAsJet = WTATrigObjToWTAJet(seed);
                    for(int j = jet_N - 1; j >= 0; j--)
                    { // It's important to read bits from the right, lower et jets first
                        if(associate_bit.at(j) == 1)
                        { // If Incoming seed is the highest, and there is an associated old jth jet, pop out jth jet
                            InsertToConstList(WTAJetToWTATrigObj(m_SeedList.at(j))); // Move jth jet to constituent list
                            m_SeedList.erase(m_SeedList.begin() + j); // Then, erase jth jet
                        }
                    }
                    m_SeedList.push_back(std::move(IncomingSeedAsJet)); // Add this new incoming seed as Jet, well localized protojet
                }
                else
                { // Incoming Seed is not the highest et, add seed to the ConstituentList
                    InsertToConstList(seed);
                }
            } // At least one association loop
            else
            { // When there is no association
                WTAJet IncomingSeedAsJet = WTATrigObjToWTAJet(seed);
                m_SeedList.push_back(std::move(IncomingSeedAsJet)); // Add seed to the SeedList
            }
        } // Main merging loop
        SortByPt(m_SeedList); // PtSort the SeedList. This is important for Baseline seed cleaning!
        ResizeSeedNConstLists(); // Resize SeedList, and ConstList
        if(m_VERBOSE)PrintSeedList(); // Print SeedList, for debug
    } // seed loop
    if(m_DEBUG){
        PrintSeedList();
        std::cout << "Baseline Seed Cleaning Done......" << std::endl;
    }
}

SetDEBUG()

void WTAConeMaker::SetDEBUG ( )

inline

Definition at line 84 of file WTAConeMaker.h.

{m_DEBUG = true;}; // Printout for debug

SetVERBOSE()

void WTAConeMaker::SetVERBOSE ( )

inline

Definition at line 85 of file WTAConeMaker.h.

{m_DEBUG = true; m_VERBOSE = true;};

WTAJetToWTATrigObj()

WTATrigObj WTAConeMaker::WTAJetToWTATrigObj ( const WTAJet & jet )

inline

Definition at line 204 of file WTAConeMaker.h.

{
    WTATrigObj thisobj(jet.pt(), jet.eta(), jet.phi(), jet.m(), jet.idx());
    return thisobj;
}

WTATrigObjToWTAJet()

WTAJet WTAConeMaker::WTATrigObjToWTAJet ( const WTATrigObj & obj )

inline

Definition at line 198 of file WTAConeMaker.h.

{
    WTAJet thisjet(obj.pt(), obj.eta(), obj.phi(), obj.m(), obj.idx());
    return thisjet;
}

Member Data Documentation

m_ConstituentList

std::vector<WTATrigObj> WTAConeMaker::m_ConstituentList

protected

Definition at line 92 of file WTAConeMaker.h.

m_DEBUG

bool WTAConeMaker::m_DEBUG

protected

Definition at line 96 of file WTAConeMaker.h.

m_SeedList

std::vector<WTAJet> WTAConeMaker::m_SeedList

protected

Definition at line 94 of file WTAConeMaker.h.

m_SeedSortingList

std::vector<WTATrigObj> WTAConeMaker::m_SeedSortingList

protected

Definition at line 93 of file WTAConeMaker.h.

m_VERBOSE

bool WTAConeMaker::m_VERBOSE

protected

Definition at line 97 of file WTAConeMaker.h.

m_WTAConeMakerParameter

WTAParameters WTAConeMaker::m_WTAConeMakerParameter

Definition at line 89 of file WTAConeMaker.h.

---

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

• WTAConeMaker.h