FastMultEst.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#ifndef ALICEO2_ITS_FASTMULTEST_
#define ALICEO2_ITS_FASTMULTEST_
 
#include "ITSMFTReconstruction/ChipMappingITS.h"
#include "DataFormatsITSMFT/ROFRecord.h"
#include "DataFormatsITSMFT/CompCluster.h"
#include <DataFormatsITSMFT/PhysTrigger.h>
#include "ITSReconstruction/FastMultEstConfig.h"
#include <gsl/span>
#include <array>
 
namespace o2
{
namespace its
{
 
struct FastMultEst {
 
  static constexpr int NLayers = o2::itsmft::ChipMappingITS::NLayers;
 
  float mult = 0.;                         
  float noisePerChip = 0.;                 
  float cov[3] = {0.};                     
  float chi2 = 0.;                         
  int nLayersUsed = 0;                     
  uint32_t lastRandomSeed = 0;             
 
  std::array<int, NLayers> nClPerLayer{0}; // measured N Cl per layer selectROFs
  FastMultEst();
 
  static uint32_t getCurrentRandomSeed();
  int selectROFs(const gsl::span<const o2::itsmft::ROFRecord> rofs, const gsl::span<const o2::itsmft::CompClusterExt> clus,
                 const gsl::span<const o2::itsmft::PhysTrigger> trig, std::vector<uint8_t>& sel);
 
  void fillNClPerLayer(const gsl::span<const o2::itsmft::CompClusterExt>& clusters);
  float process(const std::array<int, NLayers> ncl)
  {
    return FastMultEstConfig::Instance().imposeNoisePerChip > 0 ? processNoiseImposed(ncl) : processNoiseFree(ncl);
  }
  float processNoiseFree(const std::array<int, NLayers> ncl);
  float processNoiseImposed(const std::array<int, NLayers> ncl);
  float process(const gsl::span<const o2::itsmft::CompClusterExt>& clusters)
  {
    fillNClPerLayer(clusters);
    return process(nClPerLayer);
  }
  static bool sSeedSet;
 
  ClassDefNV(FastMultEst, 1);
};
 
} // namespace its
} // namespace o2
 
#endif