MCalCoefficients.hh

Go to the documentation of this file.

#ifndef _Q_CAL_COEF_H_
#define _Q_CAL_COEF_H_

#include "QCuore.hh"
#include "QModule.hh"
#include "MCalCoefficients.hh"
#include "QCalPeak.hh"
#include "QVector.hh"
#include <vector>
#include <map>
#include <set>

class MCalCoefficients : public QModule {
  
public:
  MCalCoefficients(QSequence* s);
  
  virtual ~MCalCoefficients();
  
  // Init method is called before event loop
  virtual void Init();
  
  // Doit method is called for each event, getting the event as argument
  virtual QEvent* Do( QEvent* ev);
  
  // Done method is called after event loop
  virtual void Done();

//  virtual void StoreData(QEvent* ev);

private:
        std::map<int, std::vector<double> > fAmplitudes;

        std::map<int, double> fMaxAmplitude;

        double fMaxEnergy;
        double fMinEnergy;

        std::set<int> fRuns;
        
        std::vector<QCalPeak*> fCalPeaks1;

        std::vector<QCalPeak*> fCalPeaks2;
        
        TF1* fCalFunction;

        double *fCalFunctionInit;

        std::string fTempFilename;

        std::map<double,double> fChisquares;
        std::map<double,double> fReducedChisquares;
        
        void InitCalPeaks(std::string source);

        bool AutoSeed(int channel, double resolution=-1, int iter=0, double binscale=1);

        double NearestNeighbor(double energy);
};

std::vector<std::vector<int> > GenerateCombinations(std::vector<int> object, unsigned int r);
std::vector<int> RemoveElements(std::vector<int> objects, unsigned int i) ;

#endif

---