DemuxFast.h

Go to the documentation of this file.

#ifndef DEMUX_FAST_H
#define DEMUX_FAST_H

class DemuxFast;             // Forward declaration of THIS class 
class MomNavigator;        // Forward declaration of this class
class CandDigitListHandle; // Forward declaration of this class

//class PlanePossibility;
#include "Plex/PlexSEIdAltL.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "UgliGeometry/UgliStripHandle.h"
#include "JobControl/JobCModule.h"
#include "MyDemuxPatternMaster.h"
#include "TTree.h"
#include "TNamed.h"
#include "TCanvas.h" // ROOT include file for drawing canvases
#include "TH2.h"     // ROOT include file for 2d histograms

class TCanvas;
class TPad;
class TPolyMarker;
class TH2F;

const int DISPLAY_TYPE = 4;

enum DemuxSearchType{SEARCH_BACKWARDS,SEARCH_GAP,SEARCH_GAP_F,SEARCH_GAP_B,SEARCH_FORWARDS,SEARCH_SINGLES}; 



typedef struct {
  PlexSEIdAltL* altList1;
  Int_t         alt1;
  PlexSEIdAltL* altList2;
  Int_t         alt2;
  Int_t         status;
}DemuxPossibility;

typedef struct {
  Int_t         iplane;
  Int_t         lowplane;
  Int_t         highplane;
  DemuxSearchType searchType;
  bool            goodTactic;
}DemuxSearchTactic;



class DemuxFast : public JobCModule, public TObject
{

 public:

  DemuxFast();
  void BeginJob();                             
  JobCResult Ana(const MomNavigator *mom);
  JobCResult Reco(MomNavigator *mom);  
  void EndJob();                               
  void ResetMap(int ne, int nw);
  void PrintMap(int ne, int nw);
  void PrintWhatRemains();
  void MakeMap(int iplane,bool useTargets=false);
  void MakeMap(DemuxSearchTactic);
  void GroupHits(int ecount, int wcount);
  void MakeSearchTactics();
  void MakeSearchTacticsX();
  void MakeSearchTacticsY();
  void DemuxSingles(DemuxSearchTactic);
  bool SelectHits(bool makeGold, bool useGold=false);
  void DemuxHits(int ie, int iw , int is);
  void GoldHits(int ie, int iw , int is);
  void MakeGoldGuide();
  void DemuxHitE(PlexSEIdAltL* pAlt, int is);
  void DemuxHitW(PlexSEIdAltL* pAlt, int is);
  void DemuxHitE(int ie, int is);
  void DemuxHitW(int iw, int is);
  bool DemuxBigGroup(int ig,bool gold);
  void NewTactic(DemuxSearchTactic oldTactic);
  float IsItXTalk(int iew, int iplane, int iva, int ipixel);

 private:

  MyDemuxPatternMaster* pMaster;
  std::vector<PlexSEIdAltL*>PlanesAltListsE[500];
  std::vector<PlexSEIdAltL*>PlanesAltListsW[500];
  std::vector<DemuxSearchTactic>searchTactics;
  PlaneView::PlaneView_t  kView; 
  UgliGeomHandle* _pUgh;

  /* -------------Member Functions---------------- */

                               
  /* --------------Cuts--------------------------- */

  float _cutRatioLow;
  float _cutRatioHigh;
  float _cutPEmin;
  float _cutFractionMin;
  float _cutSigmaQ;
  float _cutTrackTargetQ;

  /* --------------Data Structures---------------- */

  bool _amap[100][100];
  int  _smap[100][100];
  bool _bmap[100][100];
  int  _UVmap[500];
  bool amWriting;
  bool amPSing;
  bool amHisting;
  bool _searched[500];

  float pixelE[500][3][16];
  float pixelW[500][3][16];
  float _qTotE[500];
  float _qTotW[500];
  float _qMaxE[500];
  float _qMaxW[500];

  float _goldGuide[500];
  float _sgoldGuide[500];
  std::vector <int> targets;

  PlexSEIdAltL* _pAmapE[100];
  PlexSEIdAltL* _pAmapW[100];

  int _vaChip0[18];
    
  int  _rangeWindow;

  bool  efound[100];
  bool  wfound[100];
  int   ngroup;
  int   ningroupe[100];
  int   ningroupw[100];
  int   egroup[100][100];
  int   wgroup[100][100];

  int   ecount;
  int   wcount;
  int   currentPlane;
  int   planehit[500];
  int   hitmap[500][100];
  float qhitmapE[500][100];
  float qhitmapW[500][100];
  
  int   NdemuxedHitsU;
  int   NdemuxedHitsV;
  int   NdemuxedHits;
  std::vector<int>   demuxedHitPlane;
  std::vector<int>   demuxedHitStrip;
  std::vector<float> demuxedHitQ;
 
  int   demuxedmap[500][100];
  int   goldplanehit[500];
  int   goldhitmap[500][100];

  ClassDef(DemuxFast,0);
}; 
#endif // DEMUX_FAST_H

---