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 52 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 54 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 81 of file BasicHepMCValidation.h.

81 {};

Member Function Documentation

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

inline

Definition at line 83 of file BasicHepMCValidation.h.

References beta, EnergyCorrector::c, count, eta_init, MonitorElement::Fill(), GenParticle::GenParticle, GetFinal(), isFirst(), lifetime_final, lifetime_init, p_final, p_init, pdgid, and mps_merge::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 113 of file BasicHepMCValidation.h.

References count, 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 131 of file BasicHepMCValidation.h.

References funct::abs(), edmIntegrityCheck::d, and GenParticle::GenParticle.

Referenced by Fill().

                                                                  { // includes mixing (assuming mixing is not occurring more than 5 times back and forth)
         HepMC::GenParticle* aPart = new HepMC::GenParticle(*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 122 of file BasicHepMCValidation.h.

References funct::abs(), and funct::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;
       }