QCoProjector.hh

Go to the documentation of this file.

#ifndef _Q_COPROJECTER_HH_
#define _Q_COPROJECTER_HH_

#include <vector>
#include <list>
#include <string>
#include <QCoSelector.hh>

class QCoincidence;

class TEventList;
class QCoincidence;


class QCoProjector
{


public:
  QCoProjector ();
  virtual ~QCoProjector ();

  //* @brief Set coincidences
  void SetCoincidences (const std::list<QCoincidence>& coincidences); 

  const std::list<QCoincidence>& GetCoincidences() const
   { return fCoincidences; }

  //* @brief Return Events in fCoincidences in form of root TEventlist,
  //   that can be used to filter a diana tree 
  TEventList * GetEventList ();

  //* @brief Projects arg variables in internal arrays in Event-based way
  //   arg: Root syntaxis, by now, 1 or two arguments, separated by : 
  //   Consider only events in evList
  //   To recover arrays, do this->GetV1(), this->GetV2() 
  //   Beware these arrays will be overwritten by the next call to Project...
  //   If an eventList is given, only this events will be considered
  //   Return the number of projected events
  int ProjectEventBased (std::string arg, TEventList * evList=0);

  //* @brief Projects arg variables in internal arrays in Co-based way
  //   arg: Root syntaxis, by now, 1 or two arguments, separated by : 
  //   To recover arrays, do this->GetV1(), this->GetV2() 
  //   Beware these arrays will be overwritten by the next call to Project...
  void ProjectCoBased (std::string arg);

  //* @brief Return pointers to internals arrays with projected variables
  //   Beware these arrays will be overwritten by the next call to Project...
  double * GetV1() {return fV1;}
  double * GetV2() {return fV2;}
  int GetNEvents() {if (fRegenerateEventList) RegenerateEventList(); return fNEvents;}
  int GetNCoincidences() {if (fRegenerateEventList) RegenerateEventList(); return fNCoincidences;}

protected:
  std::list<QCoincidence> fCoincidences;
  QCoSelector fSelector;
  TEventList * fEventList; 
  bool fRegenerateEventList;
  int fNEvents;
  int fNCoincidences; 
  double * fV1;
  double * fV2;


  // Auxiliar
protected:
  void RegenerateEventList();
  void CreateArrays();
  std::vector<std::string> ScanArg(std::string arg, int &narg);
};

#endif

---