Description: <one line="" class="" summary>="">
Implementation: <Notes on="" implementation>="">
Definition at line 49 of file GenHIEventProducer.cc.
Definition at line 101 of file GenHIEventProducer.cc.
References b, begin, ALCARECOTkAlJpsiMuMu_cff::charge, DEFINE_FWK_MODULE, end, PVValHelper::eta, Exception, edm::Event::getByToken(), edm::EventSetup::getData(), edm::HepMCProduct::GetEvent(), ecalTB2006H4_GenSimDigiReco_cfg::mix, eostools::move(), hpstanc_transforms::nCharged, np, npart, HiggsValidation_cfi::pdg_id, edm::Handle< T >::product(), EnergyCorrector::pt, and edm::Event::put().
116 int nChargedPtCut = 0;
117 int nChargedPtCutMR = 0;
122 double TotEnergy = 0;
130 <<
"Mixing has "<<
mix.size()<<
" sub-events, should have been at least 1" 137 HepMC::GenEvent::particle_const_iterator
begin = evt->particles_begin();
138 HepMC::GenEvent::particle_const_iterator
end = evt->particles_end();
139 for(HepMC::GenEvent::particle_const_iterator it = begin; it !=
end; ++it){
140 if((*it)->status() != 1)
continue;
141 int pdg_id = (*it)->pdg_id();
143 int charge =
static_cast<int>(part->charge());
145 if(charge == 0)
continue;
146 float pt = (*it)->momentum().perp();
147 float eta = (*it)->momentum().eta();
148 float energy = (*it)->momentum().e();
164 if(fabs(eta) > 0.5)
continue;
169 const HepMC::HeavyIon*
hi = evt->heavy_ion();
172 ncoll = ncoll + hi->Ncoll();
173 nhard = nhard + hi->Ncoll_hard();
174 int np = hi->Npart_proj() + hi->Npart_targ();
177 b = hi->impact_parameter();
178 phi = hi->event_plane_angle();
179 ecc = hi->eccentricity();
189 meanPtMR /= nChargedMR;
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::EDGetTokenT< CrossingFrame< edm::HepMCProduct > > hepmcSrc_
edm::ESHandle< ParticleDataTable > pdt
bool getData(T &iHolder) const
HepPDT::ParticleData ParticleData
const HepMC::GenEvent * GetEvent() const
T const * product() const