poissonEM.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.
 
// 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 O2_MCH_POISSONEM_H_
#define O2_MCH_POISSONEM_H_
 
namespace o2
{
namespace mch
{
 
// namespace  PEM {
// public :
static constexpr int nMacroIterations = 8;
static int nIterations[nMacroIterations] = {5, 10, 10, 10,
                                            10, 10, 10, 30};
/*
static constexpr double minPadResidues[nMacroIterations] = {2.0, 2.0, 1.5, 1.5,
                                                            1.0, 1.0, 0.5, 0.5};
*/
static constexpr double minPadResidues[nMacroIterations] = {0.8, 0.6, 0.5, 0.4,
                                                            0.3, 0.25, 0.2, 0.15};
 
std::pair<double, double> PoissonEMLoop(const Pads& pads, Pads& pixels,
                                        const double* Cij, Mask_t* maskCij,
                                        int qCutMode, double minPadResidu,
                                        int nItMax);
// static double computeChiSquare( const Pads &pads, const double
// *qPredictedPads);
//};
 
} // namespace mch
} // namespace o2
 
#endif // O2_MCH_POISSONEM_H_