ReconstructionUtils.h

Go to the documentation of this file.

/*
 * ReconstructionUtils.h
 *
 *  Created on: 4 Mar 2022
 *      Author: Valerio Formato
 */

#ifndef HERD_RECONSTRUCTIONUTILS_H
#define HERD_RECONSTRUCTIONUTILS_H

// HerdSoftware headers
#include "common/BinnedProfile.h"
#include "dataobjects/CaloGeoParams.h"
#include "dataobjects/CaloHitsRefFrame.h"
#include "dataobjects/CaloShower.h"
#include "dataobjects/Line.h"
#include "dataobjects/Point.h"

namespace Herd {
namespace Reconstruction {

Point CaloTwoThreshShowerStartingPoint(const CaloHits &caloHits, const Line &showerAxis,
                                       const CaloGeoParams &caloGeoParams, float thr1, float thr2, float radius,
                                       float distSecond, float threshold = 0.0f);

float CaloLocalEnergy(const CaloHits &hits, const Line &line, const CaloGeoParams &geoParams, float radius,
                      float innerRadius = 0, Point refPoint = defaultValue<Point>,
                      float startPos = std::numeric_limits<float>::lowest(),
                      float endPos = std::numeric_limits<float>::max(), float threshold = 0.0f);

float CaloLongitudinalMoment(unsigned int order, const CaloHits &hits, const Vec3D &dir, const Point &refPoint,
                             const CaloGeoParams &geoParams, float threshold = 0.0f);

float CaloLateralMoment(unsigned int order, const CaloHits &hits, const Line &dir, const CaloGeoParams &geoParams,
                        float threshold = 0.0f);

BinnedProfile<float> CaloSignalLongitudinalProfile(const CaloHits &hits, const Line &line,
                                                   const CaloGeoParams &geoParams, float maxRadius,
                                                   float longitudinalDistStep, float threshold = 0.0f);

BinnedProfile<float> CaloLateralMomentLongitudinalProfile(const CaloHits &hits, const Line &line,
                                                          const CaloGeoParams &geoParams, float maxRadius,
                                                          float longitudinalDistStep, float momentOrder,
                                                          float threshold = 0.0f);

} // namespace Reconstruction
} // namespace Herd

#endif