| 471 |
|
nVertex = event->vertices.size(); |
| 472 |
|
|
| 473 |
|
tree::Particle thisGenParticle; |
| 474 |
+ |
cout << " new event " << endl; |
| 475 |
|
for( std::vector<susy::Particle>::iterator it = event->genParticles.begin(); it != event->genParticles.end(); ++it ) { |
| 476 |
|
if( it->status == 3 ) { // hard interaction |
| 477 |
+ |
cout << "hard interaction mit status " << it->pdgId << endl; |
| 478 |
|
switch( std::abs(it->pdgId) ) { |
| 479 |
|
thisGenParticle.pt = it->momentum.Pt(); |
| 480 |
|
thisGenParticle.eta = it->momentum.Eta(); |
| 481 |
|
thisGenParticle.phi = it->momentum.Phi(); |
| 482 |
|
case 22: // photon |
| 483 |
+ |
cout << "found photon with pt " << thisGenParticle.pt << endl; |
| 484 |
|
if( thisGenParticle.pt > 75 ) |
| 485 |
|
genPhoton.push_back( thisGenParticle ); |
| 486 |
|
case 11: // electron |
| 487 |
+ |
cout << "found electron with pt " << thisGenParticle.pt << endl; |
| 488 |
|
if( thisGenParticle.pt > 20 ) // pt cut is lower to estimate fake rate in all pt bins |
| 489 |
|
genElectron.push_back( thisGenParticle ); |
| 490 |
|
} |
| 491 |
|
} |
| 492 |
+ |
else if (it->status == 2) |
| 493 |
+ |
cout <<"zwischending mit status " << it->pdgId << endl; |
| 494 |
+ |
else if (it->status == 1) |
| 495 |
+ |
cout <<"endzustand mit status " << it->pdgId << endl; |
| 496 |
+ |
|
| 497 |
|
} |
| 498 |
|
|
| 499 |
|
tree->Fill(); |
