Skip to content
Snippets Groups Projects
EvtDataBaseList.hh 1.88 KiB
Newer Older
Bertram Kopf's avatar
Bertram Kopf committed
// EvtDataBaseList class definition file. -*- C++ -*-
// Copyright 2012 Bertram Kopf


#include "qft++/topincludes/relativistic-quantum-mechanics.hh"
#include <boost/shared_ptr.hpp>
Bertram Kopf's avatar
Bertram Kopf committed
#include <boost/unordered_map.hpp>
#include <map>

#include "Utils/PawianCollectionUtils.hh"
#include "PwaUtils/DataUtils.hh"

Bertram Kopf's avatar
Bertram Kopf committed
typedef std::map<int, Vector4<double> > mapInt4Vec;
typedef std::map<int, map<Spin,map<Spin,map<Spin,complex<double> > > > > mapIntSpinComplex;
typedef std::map<std::string, Vector4<double> > mapString4Vec;
typedef std::map<std::string, map<Spin,map<Spin,map<Spin,complex<double> > > > > mapStringSpinComplex;
typedef std::map<std::string, double> mapStringDouble;
Bertram Kopf's avatar
Bertram Kopf committed

struct EvtData {
Bertram Kopf's avatar
Bertram Kopf committed
  mapInt4Vec FourVecsProd;
  mapInt4Vec FourVecsDec;
  mapIntSpinComplex WignerDsProd;
  mapIntSpinComplex WignerDsDec;
  mapString4Vec FourVecsString;
  mapStringSpinComplex WignerDsString;
  mapStringDouble DoubleString;
  double evtWeight;
Bertram Kopf's avatar
Bertram Kopf committed
  int evtNo;
class AbsEnv;

class EvtDataBaseList {

public:
  EvtDataBaseList(AbsEnv* theEnv);
  virtual ~EvtDataBaseList();

Bertram Kopf's avatar
Bertram Kopf committed
  void read(EventList& evtListData, EventList& evtListMc);

  const std::vector<EvtData*> getDataVecs() const { return _evtDataList; }
  const std::vector<EvtData*> getMcVecs() const { return _mcDataList; }

  double NoOfWeightedDataEvts() const {return _noOfWeightedDataEvts;}
  double NoOfWeightedMcEvts() const {return _noOfWeightedMcEvts;}
  void ratioMcToData (double mcToDataRatio) {_mcToDataRatio=mcToDataRatio;}

protected:
  std::vector<EvtData*> _evtDataList;
  std::vector<EvtData*> _mcDataList;
  virtual void read4Vecs(EventList& evtList, std::vector<EvtData*>& theEvtList, double& evtWeightSum, int maxEvts );
  virtual std::string getName(std::vector<Particle*>& theVec);
Bertram Kopf's avatar
Bertram Kopf committed

  double _noOfWeightedDataEvts;
  double _noOfWeightedMcEvts;
  int _mcToDataRatio;
  bool _alreadyRead;
  int _evtNoAll;
  AbsEnv* _absEnv;
Bertram Kopf's avatar
Bertram Kopf committed