HGCSSRecoHit.cc

Go to the documentation of this file.

#include <iomanip>
#include <cmath>
#include <stdlib.h>

#include "HGCal/TBStandaloneSimulator/interface/HGCSSRecoHit.hh"

HGCSSRecoHit::HGCSSRecoHit(const HGCSSSimHit & aSimHit,
                           const bool isScintillator,
                           const HGCSSGeometryConversion & aGeom)
{
    energy_ = aSimHit.energy();
    adcCounts_ = 0;
    zpos_ = aSimHit.get_z();


    layer_ = aSimHit.layer();
    noiseFrac_ = 0;

    std::pair<double, double> xy = aSimHit.get_xy(isScintillator, aGeom);
    xpos_ = xy.first;
    ypos_ = xy.second;

    time_ = aSimHit.time();

    //cellid encoding:
    //bool x_side = x>0 ? true : false;
    //bool y_side = y>0 ? true : false;
    //unsigned x_cell = static_cast<unsigned>(fabs(x)/(cellSize*granularity));
    //unsigned y_cell = static_cast<unsigned>(fabs(y)/(cellSize*granularity));
    //encodeCellId(x_side,y_side,x_cell,y_cell,granularity);

}


double HGCSSRecoHit::theta() const
{
    return 2 * atan(exp(-1.*eta()));
}

double HGCSSRecoHit::eta() const
{
    return position().eta();
}

double HGCSSRecoHit::phi() const
{
    return position().phi();
}

// void HGCSSRecoHit::encodeCellId(const bool x_side,const bool y_side,const unsigned x_cell,const unsigned y_cell, const unsigned granularity){
//   cellid_ =
//     x_side | ((x_cell & 0xFFF)<<1) |
//     (y_side<<13) | ((y_cell & 0xFFF)<<14) |
//     ((granularity & 0x3F) <<26) ;

//   // std::cout << " Cross-check of encoding: cellid=" << cellid_ << std::endl
//   //         << " x_side " << x_side << " " << get_x_side() << std::endl
//   //         << " y_side " << y_side << " " << get_y_side() << std::endl
//   //         << " x_cell " << x_cell << " " << get_x_cell() << std::endl
//   //         << " y_cell " << y_cell << " " << get_y_cell() << std::endl
//   //        << " granularity " << granularity << " " << getGranularity() << std::endl
//   //   ;
// }

void HGCSSRecoHit::Add(const HGCSSSimHit & aSimHit)
{
    time_ = time_ * energy_;
    energy_ += aSimHit.energy();
    time_ += aSimHit.energy() * aSimHit.time();
    if (energy_ > 0) time_ = time_ / energy_;
}

void HGCSSRecoHit::Print(std::ostream & aOs) const
{
    aOs << "====================================" << std::endl
        << " = Layer " << layer_ //<< " cellid " << cellid_
        << std::endl
        << " = Energy " << energy_ << " noiseFrac " << noiseFrac_ << std::endl
        << " = Digi E " << adcCounts_ << " adcCounts." << std::endl
        << " = time " << time_ << std::endl
        << "====================================" << std::endl;

}