24 = consumes<GenParticleCollection>(inputTagGenParticles_);
32 produces<GenJetCollection>();
47 std::ostringstream err;
48 err<<
" cannot get collection: "
49 <<inputTagGenParticles_<<std::endl;
54 std::auto_ptr<GenJetCollection>
57 using namespace GenParticlesHelper;
61 allStatus2Taus, 15, 2);
63 for (
IGR iTau=allStatus2Taus.begin(); iTau!=allStatus2Taus.end(); ++iTau ) {
73 if ( status2TauDaughters.
size() > 0 )
continue;
81 cout<<
"tau "<<(*iTau)<<endl;
84 igr!= descendents.
end(); ++igr ) {
86 int absPdgId =
abs((*igr)->pdgId());
89 if(!includeNeutrinos_ ) {
97 cout<<
"\t"<<(*igr)<<endl;
99 charge += (*igr)->charge();
100 sumVisMom += (*igr)->p4();
104 constituents.push_back(
refToPtr( *igr) );
110 GenJet jet( sumVisMom, vertex, specific, constituents);
112 if (charge != (*iTau)->charge() )
113 std::cout<<
" charge of Tau: " << (*iTau) <<
" not equal to charge of sum of charge of all descendents. " <<
std::cout;
116 pOutVisTaus->push_back( jet );
119 iEvent.
put( pOutVisTaus );
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
virtual void setCharge(Charge q) GCC11_FINAL
set electric charge
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
bool getByToken(EDGetToken token, Handle< PROD > &result) const
#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.
Abs< T >::type abs(const T &t)
Jets made from MC generator particles.
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)