SimG4HcalHitCluster.cc

Go to the documentation of this file.

// File: SimG4HcalHitCluster.cc
// Description: Cluster class in SimG4HcalValidation
#include "Validation/HcalHits/interface/SimG4HcalHitCluster.h"

#include "CLHEP/Vector/LorentzVector.h"

SimG4HcalHitCluster::SimG4HcalHitCluster(): ec(0), etac(0), phic(0) {}

SimG4HcalHitCluster::~SimG4HcalHitCluster() {}

bool SimG4HcalHitCluster::operator<(const SimG4HcalHitCluster& cluster) const {
  return (ec/cosh(etac) < cluster.e()/cosh(cluster.eta())) ? false : true ;
}

SimG4HcalHitCluster& SimG4HcalHitCluster::operator+=(const CaloHit& hit) {

  hitsc.push_back(hit);
 
  if (ec == 0. && etac == 0. && phic == 0.) {
    ec   = hit.e();
    etac = hit.eta();
    phic = hit.phi();
  } else {   
    // cluster px,py,pz
    double et = ec / my_cosh(etac);
    double px = et * cos(phic);
    double py = et * sin(phic);
    double pz = et * my_sinh(etac); 

    CLHEP::HepLorentzVector clusHLV(px,py,pz,ec);
      
    // hit px,py,pz
    double eh   = hit.e();
    double etah = hit.eta();
    double phih = hit.phi();
    et = eh / my_cosh(etah);
    px = et * cos(phih);
    py = et * sin(phih);
    pz = et * my_sinh(etah); 
      
    CLHEP::HepLorentzVector hitHLV(px,py,pz,eh);
      
    // clus + hit
    clusHLV += hitHLV;
      
    double theta  = clusHLV.theta();
    etac = -log(tan(theta/2.));
    phic = clusHLV.phi();
    ec   = clusHLV.t();
  }

  return *this;
}

double SimG4HcalHitCluster::collectEcalEnergyR() {

  double sum = 0.;
  std::vector<CaloHit>::iterator itr;

  for (itr = hitsc.begin(); itr < hitsc.end(); itr++) {
    if (itr->det() == 10 || itr->det() == 11 || itr->det() == 12) {
      sum += itr->e(); 
    }    
  }
  return sum;
}

std::ostream& operator<<(std::ostream& os, const SimG4HcalHitCluster& cluster){
  os << " SimG4HcalHitCluster:: E " << cluster.e() << "  eta " << cluster.eta()
     << "  phi " << cluster.phi();
  return os;
}