31 produces<GenJetCollection>();
46 std::ostringstream err;
47 err<<
" cannot get collection: "
48 <<inputTagGenParticles_<<std::endl;
53 std::auto_ptr<GenJetCollection>
56 using namespace GenParticlesHelper;
60 allStatus2Taus, 15, 2);
62 for (
IGR iTau=allStatus2Taus.begin(); iTau!=allStatus2Taus.end(); ++iTau ) {
72 if ( status2TauDaughters.
size() > 0 )
continue;
80 cout<<
"tau "<<(*iTau)<<endl;
83 igr!= descendents.
end(); ++igr ) {
85 int absPdgId =
abs((*igr)->pdgId());
88 if(!includeNeutrinos_ ) {
96 cout<<
"\t"<<(*igr)<<endl;
98 charge += (*igr)->charge();
99 sumVisMom += (*igr)->p4();
103 constituents.push_back(
refToPtr( *igr) );
109 GenJet jet( sumVisMom, vertex, specific, constituents);
111 if (charge != (*iTau)->charge() )
112 std::cout<<
" charge of Tau: " << (*iTau) <<
" not equal to charge of sum of charge of all descendents. " <<
std::cout;
115 pOutVisTaus->push_back( jet );
118 iEvent.
put( pOutVisTaus );
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
Ptr< typename C::value_type > refToPtr(Ref< C, typename C::value_type, refhelper::FindUsingAdvance< C, typename C::value_type > > const &ref)
void findParticles(const reco::GenParticleCollection &sourceParticles, reco::GenParticleRefVector &particleRefs, int pdgId, int status)
find all particles of a given pdgId and status
virtual void setCharge(Charge q)
set electric charge
#define DEFINE_FWK_MODULE(type)
std::vector< GenJet > GenJetCollection
collection of GenJet objects
std::vector< Constituent > Constituents
int Charge
electric charge type
const_iterator end() const
Termination of iteration.
const_iterator begin() const
Initialize an iterator over the RefVector.
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Jets made from MC generator particles.
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
builds a GenJet from the visible daughters of each status 2 tau in the event.
XYZPointD XYZPoint
point in space with cartesian internal representation
virtual void produce(edm::Event &, const edm::EventSetup &)
size_type size() const
Size of the RefVector.
TauGenJetProducer(const edm::ParameterSet &)
void findDescendents(const reco::GenParticleRef &base, reco::GenParticleRefVector &descendents, int status, int pdgId=0)
find all descendents of a given status and pdgId (recursive)