FSimTrack.h File Reference

#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"
#include "SimDataFormats/Track/interface/SimTrack.h"
#include "SimDataFormats/Vertex/interface/SimVertex.h"
#include "CommonTools/BaseParticlePropagator/interface/RawParticle.h"
#include <vector>
#include <iosfwd>
#include "FastSimulation/Event/interface/FSimTrack.icc"

Go to the source code of this file.

Classes
class	FSimTrack

Namespaces
	HepMC

Functions
std::ostream &	operator<< (std::ostream &o, const FSimTrack &t)

Function Documentation

operator<<()

std::ostream& operator<<	(	std::ostream &	o,
		const FSimTrack &	t
	)

Definition at line 146 of file FSimTrack.cc.

References alignBH_cfg::fixed, mps_fire::i, mergeVDriftHistosByStation::name, EcalTangentSkim_cfg::o, AlCaHLTBitMon_QueryRunRegistry::string, and submitPVValidationJobs::t.

                                                          {
  std::string name = t.particleInfo() ? t.particleInfo()->name() : "Unknown";
  const XYZTLorentzVector& momentum1 = t.momentum();
  XYZVector vertex1 = t.vertex().position().Vect();
  int vertexId1 = t.vertex().id();

  o.setf(std::ios::fixed, std::ios::floatfield);
  o.setf(std::ios::right, std::ios::adjustfield);

  o << std::setw(4) << t.id() << "; " << std::setw(4) << t.genpartIndex() << "; " << name << "\n\t\t";

  //for(unsigned int k=0;k<11-name.length() && k<12; k++) o << "; ";

  o << "Position: " << std::setw(6) << std::setprecision(2) << momentum1.eta() << "; " << std::setw(6)
    << std::setprecision(2) << momentum1.phi() << "; " << std::setw(6) << std::setprecision(2) << momentum1.pt() << "; "
    << std::setw(6) << std::setprecision(2) << momentum1.e() << "; " << std::setw(4) << vertexId1 << "; "
    << std::setw(6) << std::setprecision(1) << vertex1.x() << "; " << std::setw(6) << std::setprecision(1)
    << vertex1.y() << "; " << std::setw(6) << std::setprecision(1) << vertex1.z() << "; " << std::setw(4)
    << (t.particleInfo() ? (t.noMother() ? -1 : t.mother().id()) : -1) << "; "
    << "\n\t\t";

  if (!t.noEndVertex()) {
    XYZTLorentzVector vertex2 = t.endVertex().position();
    int vertexId2 = t.endVertex().id();

    o << "Decayed particle: " << std::setw(4) << vertexId2 << "; " << std::setw(6) << std::setprecision(2)
      << vertex2.eta() << "; " << std::setw(6) << std::setprecision(2) << vertex2.phi() << "; " << std::setw(5)
      << std::setprecision(1) << vertex2.pt() << "; " << std::setw(6) << std::setprecision(1) << vertex2.z() << "; ";
    for (int i = 0; i < t.nDaughters(); ++i)
      o << std::setw(4) << t.daughter(i).id() << "; "
        << "\n\t\t";

  } else {
    if (t.onLayer1()) {
      XYZTLorentzVector vertex2 = t.layer1Entrance().vertex();

      o << "Layer 1: " << std::setw(4) << -t.onLayer1() << "; " << std::setw(6) << std::setprecision(2) << vertex2.eta()
        << "; " << std::setw(6) << std::setprecision(2) << vertex2.phi() << "; " << std::setw(5) << std::setprecision(1)
        << vertex2.pt() << "; " << std::setw(6) << std::setprecision(1) << vertex2.z() << "; " << std::setw(6)
        << std::setprecision(2) << t.layer1Entrance().pt() << "; " << std::setw(6) << std::setprecision(2)
        << t.layer1Entrance().e() << "; "
        << "\n\t\t";
    }
    if (t.onLayer2()) {
      XYZTLorentzVector vertex2 = t.layer2Entrance().vertex();

      o << "Layer 2: " << std::setw(4) << -t.onLayer2() << "; " << std::setw(6) << std::setprecision(2) << vertex2.eta()
        << "; " << std::setw(6) << std::setprecision(2) << vertex2.phi() << "; " << std::setw(5) << std::setprecision(1)
        << vertex2.pt() << "; " << std::setw(6) << std::setprecision(1) << vertex2.z() << "; " << std::setw(6)
        << std::setprecision(2) << t.layer2Entrance().pt() << "; " << std::setw(6) << std::setprecision(2)
        << t.layer2Entrance().e() << "; "
        << "\n\t\t";
    }
    //if ( t.onEcal() ) {

    XYZTLorentzVector vertex2 = t.ecalEntrance().vertex();

    o << "ECAL: " << std::setw(4) << -t.onEcal() << "; " << std::setw(6) << std::setprecision(2) << vertex2.eta()
      << "; " << std::setw(6) << std::setprecision(2) << vertex2.phi() << "; " << std::setw(5) << std::setprecision(1)
      << vertex2.pt() << "; " << std::setw(6) << std::setprecision(1) << vertex2.z() << "; " << std::setw(6)
      << std::setprecision(2) << t.ecalEntrance().pt() << "; " << std::setw(6) << std::setprecision(2)
      << t.ecalEntrance().e() << "; "
      << "\n\t\t";
    //}
    //if ( t.onHcal() ) {

    vertex2 = t.hcalEntrance().vertex();

    o << "HCAL: " << std::setw(4) << -t.onHcal() << "; " << std::setw(6) << std::setprecision(2) << vertex2.eta()
      << "; " << std::setw(6) << std::setprecision(2) << vertex2.phi() << "; " << std::setw(5) << std::setprecision(1)
      << vertex2.pt() << "; " << std::setw(6) << std::setprecision(1) << vertex2.z() << "; " << std::setw(6)
      << std::setprecision(2) << t.hcalEntrance().pt() << "; " << std::setw(6) << std::setprecision(2)
      << t.hcalEntrance().e() << "; "
      << "\n\t\t";
    //}
    //if ( t.onVFcal() ) {

    vertex2 = t.vfcalEntrance().vertex();

    o << "VFCAL: " << std::setw(4) << -t.onVFcal() << "; " << std::setw(6) << std::setprecision(2) << vertex2.eta()
      << "; " << std::setw(6) << std::setprecision(2) << vertex2.phi() << "; " << std::setw(5) << std::setprecision(1)
      << vertex2.pt() << "; " << std::setw(6) << std::setprecision(1) << vertex2.z() << "; " << std::setw(6)
      << std::setprecision(2) << t.vfcalEntrance().pt() << "; " << std::setw(6) << std::setprecision(2)
      << t.vfcalEntrance().e() << "; "
      << "\n\t\t";
    //}
  }
  return o;
}