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 51 of file BasicHepMCValidation.h.

Constructor & Destructor Documentation

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

inline

Definition at line 53 of file BasicHepMCValidation.h.

References DQMHelper::book1dHisto(), eta_init, lifetime_final, lifetime_init, name, nlog, numberPerEvent, p_final, and p_init.

         : 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);
     }

BasicHepMCValidation::ParticleMonitor::~ParticleMonitor ( )

inline

Definition at line 104 of file BasicHepMCValidation.h.

104 {};

Member Function Documentation

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

inline

Definition at line 106 of file BasicHepMCValidation.h.

References HLT_FULL_cff::beta, c, count, eta_init, dqm::impl::MonitorElement::Fill(), GenParticle::GenParticle, GetFinal(), isFirst(), lifetime_final, lifetime_init, lt, p_final, p_init, pdgid, and histoStyle::weight.

                                                         {
       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;
     }

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

inline

Definition at line 141 of file BasicHepMCValidation.h.

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

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

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

inlineprivate

Definition at line 165 of file BasicHepMCValidation.h.

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

Referenced by Fill().

                                                                 {
       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;
     }

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

inlineprivate

Definition at line 152 of file BasicHepMCValidation.h.

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

Referenced by Fill().

                                             {
       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;
     }