SimCluster.cc

Go to the documentation of this file.

#include "SimDataFormats/CaloAnalysis/interface/SimCluster.h"

#include "DataFormats/HepMCCandidate/interface/GenParticle.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <numeric>

const unsigned int SimCluster::longLivedTag = 65536;

SimCluster::SimCluster()
{
  // No operation
}

SimCluster::SimCluster( const SimTrack& simtrk )
{
  addG4Track( simtrk );
  event_ = simtrk.eventId();
  particleId_ = simtrk.trackId();
  
  theMomentum_.SetPxPyPzE(simtrk.momentum().px(),
              simtrk.momentum().py(),
              simtrk.momentum().pz(),
              simtrk.momentum().E());
}

SimCluster::SimCluster( EncodedEventId eventID, uint32_t particleID )
{
  event_ = eventID;
  particleId_ = particleID;
}

SimCluster::~SimCluster()
{
}

std::ostream& operator<< (std::ostream& s, SimCluster const & tp)
{
    s << "CP momentum, q, ID, & Event #: "
    << tp.p4()                      << " " << tp.charge() << " "   << tp.pdgId() << " "
    << tp.eventId().bunchCrossing() << "." << tp.eventId().event() << std::endl;

    for (SimCluster::genp_iterator hepT = tp.genParticle_begin(); hepT !=  tp.genParticle_end(); ++hepT)
    {
        s << " HepMC Track Momentum " << (*hepT)->momentum().rho() << std::endl;
    }

    for (SimCluster::g4t_iterator g4T = tp.g4Track_begin(); g4T !=  tp.g4Track_end(); ++g4T)
    {
        s << " Geant Track Momentum  " << g4T->momentum() << std::endl;
        s << " Geant Track ID & type " << g4T->trackId() << " " << g4T->type() << std::endl;
        if (g4T->type() !=  tp.pdgId())
        {
            s << " Mismatch b/t SimCluster and Geant types" << std::endl;
        }
    }
    s << " # of cells = " << tp.hits_.size() 
      << ", effective cells = " << std::accumulate(tp.fractions_.begin(),tp.fractions_.end(),0.f) << std::endl;
    return s;
}