Pythia8::mergingDJRs Class Reference

Inheritance diagram for Pythia8::mergingDJRs:

Collaboration diagram for Pythia8::mergingDJRs:

Public Member Functions
	mergingDJRs ()
	~mergingDJRs ()
virtual bool	canVetoPartonLevel () override
virtual bool	doVetoPartonLevel (const Event &event) override
virtual void	getDJR (const Event &event)
virtual bool	initAfterBeams () override

Private Attributes
int	m_power = 1
double	m_etaMax = 10.
double	m_radius = 1.0
double	m_pTjetMin = 10.
Pythia8::SlowJetHook *	m_sjHookPtrIn = 0
bool	m_useFJcoreIn = false
bool	m_useStandardRin = true
int	m_nSel = 2
int	m_massSetIn = 1
std::unique_ptr< Pythia8::SlowJet >	m_slowJet
std::unique_ptr< TH1F >	m_HistDjr
std::unique_ptr< TH1F >	m_HistDjr2
std::unique_ptr< TFile >	m_outFile
vector< double >	m_result
vector< double >	m_DJR
Event	m_workEventJet

Detailed Description

Definition at line 20 of file mergingDJRs.cxx.

Constructor & Destructor Documentation

mergingDJRs()

Pythia8::mergingDJRs::mergingDJRs ( )

inline

Definition at line 44 of file mergingDJRs.cxx.

                {
    // Slowjet pointer
    m_slowJet = std::make_unique<Pythia8::SlowJet>(
        m_power, m_radius, m_pTjetMin, m_etaMax, m_nSel, m_massSetIn,
        m_sjHookPtrIn, m_useFJcoreIn, m_useStandardRin);
 
    // ROOT histograms and the output file where to save them
    m_HistDjr = std::make_unique<TH1F>("HistDjr", "The first DJR", 100, 0.0, 3.0);
    m_HistDjr2 =
        std::make_unique<TH1F>("HistDjr2", "The second DJR", 100, 0.0, 3.0);
    m_outFile = std::make_unique<TFile>("hist-DJR.root", "RECREATE");
 
    std::cout << "**********************************************************"
              << std::endl;
    std::cout << "*                                                        *"
              << std::endl;
    std::cout << "*   the jet merging userhook CKKWL DJRS is working   *"
              << std::endl;
    std::cout << "*                                                        *"
              << std::endl;
    std::cout << "**********************************************************"
              << std::endl;
  }

~mergingDJRs()

Pythia8::mergingDJRs::~mergingDJRs ( )

inline

Definition at line 69 of file mergingDJRs.cxx.

                 {
    m_HistDjr->Write();
    m_HistDjr2->Write();
  }

Member Function Documentation

canVetoPartonLevel()

virtual bool Pythia8::mergingDJRs::canVetoPartonLevel ( )

inlineoverridevirtual

Definition at line 75 of file mergingDJRs.cxx.

{ return true; }

doVetoPartonLevel()

bool Pythia8::mergingDJRs::doVetoPartonLevel ( const Event &  event )

inlineoverridevirtual

Definition at line 91 of file mergingDJRs.cxx.

                                                             {
 
  // subEvent method extract a list of the current partons list and save the
  // output in workEvent.
  subEvent(event);
  m_workEventJet = workEvent;
 
  // The selected particles to pass to slowjet
  for (int j = 0; j < m_workEventJet.size(); ++j) {
    if (!(m_workEventJet[j].isFinal()) || m_workEventJet[j].isLepton() ||
        m_workEventJet[j].id() == 23 || std::abs(m_workEventJet[j].id()) == 24 ||
        m_workEventJet[j].id() == 22) {
      m_workEventJet[j].statusNeg();
      continue;
    }
  }
 
  // slowjet analyze the events
  if (not m_slowJet->setup(m_workEventJet)){
    std::cout<< "setup failed in  mergingDJRs::doVetoPartonLevel\n";
    return false;
  }
 
  // Call getDJR and store the DJRs vector
  getDJR(m_workEventJet);
 
  // if we reached here then no veto
  return false;
}

getDJR()

void Pythia8::mergingDJRs::getDJR ( const Event &  event )

inlinevirtual

Definition at line 122 of file mergingDJRs.cxx.

                                                  {
 
  // setup slowjet pointer
  if (not m_slowJet->setup(event)){
    std::cout<< "setup failed in  mergingDJRs::getDJR\n";
    return;
  };
 
  // Clear members.
  m_DJR.clear();
  m_result.clear();
 
  while (m_slowJet->sizeAll() - m_slowJet->sizeJet() > 0) {
    // Save the next clustering scale.
    m_result.push_back(sqrt(m_slowJet->dNext()));
    // Perform step.
    m_slowJet->doStep();
  }
 
  // Save clustering scales in reverse order
  for (int i = m_result.size() - 1; i >= 0; --i) {
    m_DJR.push_back(m_result[i]);
  }
 
  // Fill the histogram and Normalize them
  double eventWeight = infoPtr->mergingWeight() * infoPtr->weight();
 
  if (m_DJR.size() > 0) {
    m_HistDjr->Fill(log10(m_DJR[0]), eventWeight);
    m_HistDjr2->Fill(log10(m_DJR[1]), eventWeight);
  }
}

initAfterBeams()

virtual bool Pythia8::mergingDJRs::initAfterBeams ( )

inlineoverridevirtual

Definition at line 82 of file mergingDJRs.cxx.

                                         {
 
    m_workEventJet.init("(workEventJet)", particleDataPtr);
 
    return true;
  }

Member Data Documentation

m_DJR

vector<double> Pythia8::mergingDJRs::m_DJR

private

Definition at line 40 of file mergingDJRs.cxx.

m_etaMax

double Pythia8::mergingDJRs::m_etaMax = 10.

private

Definition at line 25 of file mergingDJRs.cxx.

m_HistDjr

std::unique_ptr<TH1F> Pythia8::mergingDJRs::m_HistDjr

private

Definition at line 37 of file mergingDJRs.cxx.

m_HistDjr2

std::unique_ptr<TH1F> Pythia8::mergingDJRs::m_HistDjr2

private

Definition at line 37 of file mergingDJRs.cxx.

m_massSetIn

int Pythia8::mergingDJRs::m_massSetIn = 1

private

Definition at line 33 of file mergingDJRs.cxx.

m_nSel

int Pythia8::mergingDJRs::m_nSel = 2

private

Definition at line 32 of file mergingDJRs.cxx.

m_outFile

std::unique_ptr<TFile> Pythia8::mergingDJRs::m_outFile

private

Definition at line 38 of file mergingDJRs.cxx.

m_power

int Pythia8::mergingDJRs::m_power = 1

private

Definition at line 24 of file mergingDJRs.cxx.

m_pTjetMin

double Pythia8::mergingDJRs::m_pTjetMin = 10.

private

Definition at line 27 of file mergingDJRs.cxx.

m_radius

double Pythia8::mergingDJRs::m_radius = 1.0

private

Definition at line 26 of file mergingDJRs.cxx.

m_result

vector<double> Pythia8::mergingDJRs::m_result

private

Definition at line 40 of file mergingDJRs.cxx.

m_sjHookPtrIn

Pythia8::SlowJetHook* Pythia8::mergingDJRs::m_sjHookPtrIn = 0

private

Definition at line 28 of file mergingDJRs.cxx.

m_slowJet

std::unique_ptr<Pythia8::SlowJet> Pythia8::mergingDJRs::m_slowJet

private

Definition at line 34 of file mergingDJRs.cxx.

m_useFJcoreIn

bool Pythia8::mergingDJRs::m_useFJcoreIn = false

private

Definition at line 29 of file mergingDJRs.cxx.

m_useStandardRin

bool Pythia8::mergingDJRs::m_useStandardRin = true

private

Definition at line 31 of file mergingDJRs.cxx.

m_workEventJet

Event Pythia8::mergingDJRs::m_workEventJet

private

Definition at line 41 of file mergingDJRs.cxx.

---

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

• mergingDJRs.cxx