FSimTrack.h File Reference

#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"
#include "SimDataFormats/Track/interface/SimTrack.h"
#include "SimDataFormats/Vertex/interface/SimVertex.h"
#include "FastSimulation/Particle/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

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

Definition at line 118 of file FSimTrack.cc.

References FSimTrack::daughter(), FSimTrack::ecalEntrance(), FSimTrack::endVertex(), alignBH_cfg::fixed, SimTrack::genpartIndex(), FSimTrack::hcalEntrance(), mps_fire::i, FSimVertex::id(), FSimTrack::id(), FSimTrack::layer1Entrance(), FSimTrack::layer2Entrance(), FSimTrack::momentum(), FSimTrack::mother(), dataset::name, FSimTrack::nDaughters(), FSimTrack::noEndVertex(), FSimTrack::noMother(), connectstrParser::o, FSimTrack::onEcal(), FSimTrack::onHcal(), FSimTrack::onLayer1(), FSimTrack::onLayer2(), FSimTrack::onVFcal(), FSimTrack::particleInfo(), FSimVertex::position(), AlCaHLTBitMon_QueryRunRegistry::string, FSimTrack::vertex(), RawParticle::vertex(), and FSimTrack::vfcalEntrance().

                                                            {
  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;
}