27 PHG4InEvent *ineve = findNode::getClass<PHG4InEvent>(topNode,
"PHG4INEVENT");
30 cout <<
PHWHERE <<
"no PHG4INEVENT node" << endl;
33 pair<multimap<int, PHG4Particle *>::iterator, multimap<int, PHG4Particle *>::iterator > beginend = ineve->
GetParticles_Modify();
34 multimap<int, PHG4Particle *>::iterator particleiter;
37 cout <<
"PHG4InputFilter before filter" << endl;
40 particleiter = beginend.first;
41 while (particleiter != beginend.second)
45 double eta =
get_eta((particleiter->second)->get_px(), (particleiter->second)->get_py(), (particleiter->second)->get_pz());
48 multimap<int, PHG4Particle *>::iterator particleiter_cache = particleiter ;
51 particleiter = particleiter_cache;
57 double eta =
get_eta((particleiter->second)->get_px(), (particleiter->second)->get_py(), (particleiter->second)->get_pz());
60 multimap<int, PHG4Particle *>::iterator particleiter_cache = particleiter ;
63 particleiter = particleiter_cache;
69 double pt = sqrt((particleiter->second)->get_px()*(particleiter->second)->get_px()+ (particleiter->second)->get_py()*(particleiter->second)->get_py());
72 multimap<int, PHG4Particle *>::iterator particleiter_cache = particleiter ;
75 particleiter = particleiter_cache;
81 double pt = sqrt((particleiter->second)->get_px()*(particleiter->second)->get_px()+ (particleiter->second)->get_py()*(particleiter->second)->get_py());
84 multimap<int, PHG4Particle *>::iterator particleiter_cache = particleiter ;
87 particleiter = particleiter_cache;
95 cout <<
"PHG4InputFilter: after filter" << endl;
107 radius = sqrt(x * x + y * y);
108 theta = atan2(radius, z);
109 eta = -log(tan(theta / 2.));