Clusterer.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_CPV_CLUSTERER_H
#define ALICEO2_CPV_CLUSTERER_H
#include "DataFormatsCPV/Digit.h"
#include "DataFormatsCPV/Cluster.h"
#include "CPVReconstruction/FullCluster.h"
#include "SimulationDataFormat/MCTruthContainer.h"
#include "DataFormatsCPV/TriggerRecord.h"
 
namespace o2
{
namespace cpv
{
 
class Clusterer
{
 public:
  Clusterer() = default;
  ~Clusterer() = default;
 
  void initialize();
  void process(gsl::span<const Digit> digits, gsl::span<const TriggerRecord> dtr,
               const o2::dataformats::MCTruthContainer<o2::MCCompLabel>* dmc,
               std::vector<Cluster>* clusters, std::vector<TriggerRecord>* trigRec,
               o2::dataformats::MCTruthContainer<o2::MCCompLabel>* cluMC,
               std::vector<Digit>* calibDigits);
 
  void makeClusters(gsl::span<const Digit> digits);
  void evalCluProperties(gsl::span<const Digit> digits, std::vector<Cluster>* clusters,
                         const o2::dataformats::MCTruthContainer<o2::MCCompLabel>* dmc,
                         o2::dataformats::MCTruthContainer<o2::MCCompLabel>* cluMC);
 
  float responseShape(float dx, float dz); // Parameterization of EM shower
  void propagateMC(bool toRun = true) { mRunMC = toRun; }
 
  void makeUnfoldingsAndCalibDigits(gsl::span<const Digit> digits, std::vector<Digit>* calibDigits); // Find and unfold clusters with few local maxima
  void makeCalibDigits(std::vector<Digit>* calibDigits);                                             // Find clusters with 1 local maximum and make calibDigits using them
  void unfoldOneCluster(FullCluster& iniClu, char nMax, gsl::span<int> digitId, gsl::span<const Digit> digits);
 
 protected:
  static constexpr short NLMMax = 10; 
 
  bool mRunMC = false;                
  int mFirstDigitInEvent;             
  int mLastDigitInEvent;              
  std::vector<FullCluster> mClusters; 
  std::vector<Digit> mDigits;         
 
  std::vector<std::vector<float>> meInClusters = std::vector<std::vector<float>>(10, std::vector<float>(NLMMax));
  std::vector<std::vector<float>> mfij = std::vector<std::vector<float>>(10, std::vector<float>(NLMMax));
};
} // namespace cpv
} // namespace o2
 
#endif