BasicHepMCValidation::ParticleMonitor Class Reference

Public Member Functions
bool	Fill (const HepMC::GenParticle *p, double weight)
void	FillCount (double weight)
	ParticleMonitor (std::string name_, int pdgid_, DQMStore::IBooker &i, bool nlog_=false)
int	PDGID ()
	~ParticleMonitor ()

Private Member Functions
const HepMC::GenParticle *	GetFinal (const HepMC::GenParticle *p)
bool	isFirst (const HepMC::GenParticle *p)

Private Attributes
unsigned int	count
MonitorElement *	eta_init
MonitorElement *	lifetime_final
MonitorElement *	lifetime_init
std::string	name
bool	nlog
MonitorElement *	numberPerEvent
MonitorElement *	p_final
MonitorElement *	p_init
int	pdgid

Detailed Description

Definition at line 50 of file BasicHepMCValidation.h.

Constructor & Destructor Documentation

ParticleMonitor()

BasicHepMCValidation::ParticleMonitor::ParticleMonitor ( std::string  name_, int  pdgid_, DQMStore::IBooker &  i, bool  nlog_ = false  )

inline

Definition at line 52 of file BasicHepMCValidation.h.

        : name(name_), pdgid(pdgid_), count(0), nlog(nlog_) {
      DQMHelper dqm(&i);
      // Number of analyzed events
      if (!nlog) {
        numberPerEvent = dqm.book1dHisto(
            name + "Number", "Number of  " + name + "'s per event", 20, 0, 20, "No. of " + name, "Number of Events");
      } else {
        numberPerEvent = dqm.book1dHisto(name + "Number",
                                         "Number of  " + name + "'s per event",
                                         20,
                                         0,
                                         20,
                                         "log_{10}(No. of " + name + ")",
                                         "Number of Events");
      }
      p_init = dqm.book1dHisto(name + "Momentum",
                               "log_{10}(P) of the " + name + "s",
                               60,
                               -2,
                               4,
                               "log_{10}(P) (log_{10}(GeV))",
                               "Number of " + name);
 
      eta_init = dqm.book1dHisto(name + "Eta", "#eta of the " + name + "s", 100, -5., 5., "#eta", "Number of " + name);
 
      lifetime_init = dqm.book1dHisto(name + "LifeTime",
                                      "#phi of the " + name + "s",
                                      100,
                                      -15,
                                      -5,
                                      "Log_{10}(life-time^{final}) (log_{10}(s))",
                                      "Number of " + name);
 
      p_final = dqm.book1dHisto(name + "MomentumFinal",
                                "log_{10}(P^{final}) of " + name + "s at end of decay chain",
                                60,
                                -2,
                                4,
                                "log_{10}(P^{final}) (log_{10}(GeV))",
                                "Number of " + name);
 
      lifetime_final = dqm.book1dHisto(name + "LifeTimeFinal",
                                       "Log_{10}(life-time^{final}) of " + name + "s at end of decay chain",
                                       100,
                                       -15,
                                       -5,
                                       "Log_{10}(life-time^{final}) (log_{10}(s))",
                                       "Number of " + name);
    }

References eta_init, mps_fire::i, lifetime_final, lifetime_init, name, nlog, numberPerEvent, p_final, and p_init.

~ParticleMonitor()

BasicHepMCValidation::ParticleMonitor::~ParticleMonitor ( )

inline

Definition at line 103 of file BasicHepMCValidation.h.

{};

Member Function Documentation

Fill()

bool BasicHepMCValidation::ParticleMonitor::Fill ( const HepMC::GenParticle *  p, double  weight  )

inline

Definition at line 105 of file BasicHepMCValidation.h.

                                                        {
      if (p->pdg_id() == pdgid) {
        if (isFirst(p)) {
          p_init->Fill(log10(p->momentum().rho()), weight);
          eta_init->Fill(p->momentum().eta(), weight);
          const HepMC::GenParticle *pf = GetFinal(p);  // inlcude mixing
          p_final->Fill(log10(pf->momentum().rho()), weight);
          // compute lifetime...
          if (p->production_vertex() && p->end_vertex()) {
            TVector3 PV(p->production_vertex()->point3d().x(),
                        p->production_vertex()->point3d().y(),
                        p->production_vertex()->point3d().z());
            TVector3 SV(p->end_vertex()->point3d().x(), p->end_vertex()->point3d().y(), p->end_vertex()->point3d().z());
            TVector3 DL = SV - PV;
            double c(2.99792458E8), Ltau(DL.Mag() / 100) /*cm->m*/, beta(p->momentum().rho() / p->momentum().m());
            double lt = Ltau / (c * beta);
            if (lt > 1E-16)
              lifetime_init->Fill(log10(lt), weight);
            if (pf->end_vertex()) {
              TVector3 SVf(
                  pf->end_vertex()->point3d().x(), pf->end_vertex()->point3d().y(), pf->end_vertex()->point3d().z());
              DL = SVf - PV;
              Ltau = DL.Mag() / 100;
              lt = Ltau / (c * beta);
              if (lt > 1E-16)
                lifetime_final->Fill(log10(lt), weight);
            }
          }
          count++;
        }
        return true;
      }
      return false;
    }

References HLT_FULL_cff::beta, c, count, eta_init, dqm::impl::MonitorElement::Fill(), GenParticle::GenParticle, GetFinal(), isFirst(), lifetime_final, lifetime_init, AlCaHLTBitMon_ParallelJobs::p, p_final, p_init, pdgid, packedPFCandidateRefMixer_cfi::pf, nanoDQM_cfi::PV, nanoDQM_cfi::SV, and mps_merge::weight.

FillCount()

void BasicHepMCValidation::ParticleMonitor::FillCount ( double  weight )

inline

Definition at line 140 of file BasicHepMCValidation.h.

                                  {
      if (nlog)
        numberPerEvent->Fill(log10(count), weight);
      else
        numberPerEvent->Fill(count, weight);
      count = 0;
    }

References count, dqm::impl::MonitorElement::Fill(), nlog, and numberPerEvent.

GetFinal()

const HepMC::GenParticle* BasicHepMCValidation::ParticleMonitor::GetFinal ( const HepMC::GenParticle *  p )

inlineprivate

Definition at line 164 of file BasicHepMCValidation.h.

                                                                {
      const HepMC::GenParticle *aPart = p;
      for (unsigned int iMix = 0; iMix < 10; iMix++) {
        bool foundSimilar = false;
        if (aPart->end_vertex()) {
          if (aPart->end_vertex()->particles_out_size() != 0) {
            for (HepMC::GenVertex::particles_out_const_iterator d = aPart->end_vertex()->particles_out_const_begin();
                 d != aPart->end_vertex()->particles_out_const_end();
                 d++) {
              if (abs((*d)->pdg_id()) == abs(aPart->pdg_id())) {
                aPart = *d;
                foundSimilar = true;
                break;
              }
            }
          }
          if (!foundSimilar)
            break;
        }
      }
      return aPart;
    }

References funct::abs(), ztail::d, GenParticle::GenParticle, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by Fill().

isFirst()

bool BasicHepMCValidation::ParticleMonitor::isFirst ( const HepMC::GenParticle *  p )

inlineprivate

Definition at line 151 of file BasicHepMCValidation.h.

                                            {
      if (p->production_vertex()) {
        for (HepMC::GenVertex::particles_in_const_iterator m = p->production_vertex()->particles_in_const_begin();
             m != p->production_vertex()->particles_in_const_end();
             m++) {
          if (abs((*m)->pdg_id()) == abs(p->pdg_id()))
            return false;
        }
      }
      return true;
    }

References funct::abs(), visualization-live-secondInstance_cfg::m, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by Fill().

PDGID()

int BasicHepMCValidation::ParticleMonitor::PDGID ( )

inline

Definition at line 148 of file BasicHepMCValidation.h.

{ return pdgid; }

References pdgid.

Member Data Documentation

count

unsigned int BasicHepMCValidation::ParticleMonitor::count

private

Definition at line 189 of file BasicHepMCValidation.h.

Referenced by Fill(), and FillCount().

eta_init

MonitorElement * BasicHepMCValidation::ParticleMonitor::eta_init

private

Definition at line 191 of file BasicHepMCValidation.h.

Referenced by Fill(), and ParticleMonitor().

lifetime_final

MonitorElement * BasicHepMCValidation::ParticleMonitor::lifetime_final

private

Definition at line 191 of file BasicHepMCValidation.h.

Referenced by Fill(), and ParticleMonitor().

lifetime_init

MonitorElement * BasicHepMCValidation::ParticleMonitor::lifetime_init

private

Definition at line 191 of file BasicHepMCValidation.h.

Referenced by Fill(), and ParticleMonitor().

name

std::string BasicHepMCValidation::ParticleMonitor::name

private

Definition at line 187 of file BasicHepMCValidation.h.

Referenced by config.CFG::__str__(), validation.Sample::digest(), VIDSelectorBase.VIDSelectorBase::initialize(), and ParticleMonitor().

nlog

bool BasicHepMCValidation::ParticleMonitor::nlog

private

Definition at line 190 of file BasicHepMCValidation.h.

Referenced by FillCount(), and ParticleMonitor().

numberPerEvent

MonitorElement * BasicHepMCValidation::ParticleMonitor::numberPerEvent

private

Definition at line 191 of file BasicHepMCValidation.h.

Referenced by FillCount(), and ParticleMonitor().

p_final

MonitorElement * BasicHepMCValidation::ParticleMonitor::p_final

private

Definition at line 191 of file BasicHepMCValidation.h.

Referenced by Fill(), and ParticleMonitor().

p_init

MonitorElement* BasicHepMCValidation::ParticleMonitor::p_init

private

Definition at line 191 of file BasicHepMCValidation.h.

Referenced by Fill(), and ParticleMonitor().

pdgid

int BasicHepMCValidation::ParticleMonitor::pdgid

private

Definition at line 188 of file BasicHepMCValidation.h.

Referenced by Fill(), and PDGID().