CMSSW_4_4_3_patch1/src/SimDataFormats/HcalTestBeam/src/PHcalTB06Info.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:     HcalTestBeam
// Class  :     PHcalTB06Info
// 
// Implementation:
//     <Notes on implementation>
//
// Original Author:  Sunanda Banerjee
//         Created:  Tue Oct 10 10:25:44 CEST 2006
//

// system include files

// user include files
#include "SimDataFormats/HcalTestBeam/interface/PHcalTB06Info.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"



//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
PHcalTB06Info::PHcalTB06Info() {
  clear();
}

// PHcalTB06Info::PHcalTB06Info(const PHcalTB06Info& rhs) {
//    // do actual copying here;
// }

PHcalTB06Info::~PHcalTB06Info() {
}

//
// assignment operators
//
// const PHcalTB06Info& PHcalTB06Info::operator=(const PHcalTB06Info& rhs) {
//   //An exception safe implementation is
//   PHcalTB06Info temp(rhs);
//   swap(rhs);
//
//   return *this;
// }

//
// member functions
//

void PHcalTB06Info::clear() {
  nPrimary =  idBeam =0;
  eBeam = etaBeam = phiBeam = 0;

  simEtot = simEEc = simEHc = 0;

  hit = 0;
  hits.clear();

  v1EvNum = v1Type = v1Nsec = 0;
  v1X = v1Y = v1Z = v1U = v1V = v1W = 0;
  v1Px = v1Py = v1Pz = 0;
  v1Sec.clear();
}

void PHcalTB06Info::setPrimary(int primary, int id, double energy, double etav,
                               double phiv) {

  nPrimary = primary;
  idBeam   = id;
  eBeam    = (float)(energy);
  etaBeam  = (float)(etav);
  phiBeam  = (float)(phiv);
  LogDebug("SimHCalData") << "PHcalTB06Info::setPrimary: nPrimary " << nPrimary
                          << " partID " << idBeam << " initE " << eBeam 
                          << " eta " << etaBeam << " phi " << phiBeam;
}

void PHcalTB06Info::setEdep(double simtot, double sime, double simh) {

  simEtot = (float)simtot;
  simEEc  = (float)sime;
  simEHc  = (float)simh;

  LogDebug("SimHCalData") << "PHcalTB06Info::setEdep: simEtot " << simEtot 
                          << " simEEc " << simEEc << " simEHc " << simEHc;
}

void PHcalTB06Info::saveHit(unsigned int id, double eta, double phi, double e, 
                            double t) {

  int nh = hit;
  hit++;
  PHcalTB06Info::Hit newHit;
  newHit.id  = id;
  newHit.eta = (float)(eta);
  newHit.phi = (float)(phi);
  newHit.e   = (float)(e);
  newHit.t   = (float)(t);
  hits.push_back(newHit);
  LogDebug("SimHCalData") << "PHcalTB06Info::saveHit " << hit << " ID 0x" 
                          << std::hex << hits[nh].id << std::dec << " Eta " 
                          << hits[nh].eta << " Phi " << hits[nh].phi 
                          << " E " << hits[nh].e << " t " << hits[nh].t;
}

void PHcalTB06Info::setVtxPrim(int evNum, int type, double x, double y, 
                               double z, double u, double v, double w, 
                               double px, double py, double pz) {

  v1EvNum = evNum;
  v1Type  = type;
  double r= sqrt(x*x+y*y+z*z);
  v1X     = (float)(x);
  v1Y     = (float)(y);
  v1Z     = (float)(z);
  v1R     = (float)(r);
  v1Px    = (float)(px);
  v1Py    = (float)(py);
  v1Pz    = (float)(pz);
  v1U     = (float)(u);
  v1V     = (float)(v);
  v1W     = (float)(w);
  LogDebug("SimHCalData") << "PHcalTB06Info::setVtxPrim: EvNum " << v1EvNum 
                          << " Type " << v1Type << " X/Y/Z/R " << v1X << "/" 
                          << v1Y << "/" << v1Z << "/" << v1R << " Px/Py/Pz " 
                          << v1Px << "/" << v1Py << "/" << v1Pz << " U/V/W "
                          << v1U << "/" << v1V << "/" << v1W;
}

void PHcalTB06Info::setVtxSec(int id, int pdg, double px, double py, double pz,
                              double ek) {

  int ns = v1Nsec;
  v1Nsec++;
  PHcalTB06Info::Vtx newVtx;
  newVtx.trackID = id;
  newVtx.partID  = pdg;
  newVtx.px      = (float)(px);
  newVtx.py      = (float)(py);
  newVtx.pz      = (float)(pz);
  newVtx.eKin    = (float)(ek);
  v1Sec.push_back(newVtx);
  LogDebug("SimHCalData") << "PHcalTB06Info::setVtxSec " << v1Nsec << " ID " 
                          << v1Sec[ns].trackID << " PDG Code " 
                          << v1Sec[ns].partID << " Px/Py/Pz/Ek " 
                          << v1Sec[ns].px << "/" << v1Sec[ns].py << "/" 
                          << v1Sec[ns].pz << "/" << v1Sec[ns].eKin;
}