AltDeMuxCalc.h

Go to the documentation of this file.

// $Id: AltDeMuxCalc.h,v 1.3 2005/04/20 13:25:56 thomson Exp $
//
// An object for calculating useful quantities in the demuxing
// algorithm. 
//
// Author:  M. Thosmson April 2003
#ifndef ALTDEMUXCALC_H
#define ALTDEMUXCALC_H

#include "TObject.h"
#include "Plex/PlexSEIdAltL.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "UgliGeometry/UgliStripHandle.h"
#include "AltDeMuxStructures.h"

//......................................................................

class AltDeMuxCalc :  public TObject
{

 public:
  AltDeMuxCalc();
  ~AltDeMuxCalc();
  void SetClearFibreC(Float_t c){fClearFibreC=c;};
  void SetWLSFibreC(Float_t c){fWLSFibreC=c;};
  void SetScintillatorC(Float_t c){fScintillatorC=c;fDt2Dstrip=fScintillatorC/2.0/0.042;};
  void SetNumberOfStrips(Int_t n){fNumberOfStrips=n;};
  void SetFibreLengthE(Int_t iplane, PlexSEIdAltL* paltlist);
  void SetFibreLengthW(Int_t iplane, PlexSEIdAltL* paltlist);
  void SetUgli(UgliGeomHandle* p){pUgh=p;};
  void SetPlane(Int_t iplane){fPlane = iplane;}; 
  void SetView(PlaneView::PlaneView_t view){fView=view ;};
  void SetEast(PlexSEIdAltL* pListE, Int_t stripE);
  void SetWest(PlexSEIdAltL* pListW, Int_t stripW);
  void SetEastToStrip(PlexSEIdAltL* pListE, Int_t stripE);
  void SetWestToStrip(PlexSEIdAltL* pListW, Int_t stripW);
  void SetAbsTime(Double_t absTime){fAbsTime = absTime;};
  Double_t GetAbsTime(){return fAbsTime;};
  void CalcEastWest();
  void CalcEast();
  void CalcWest();
  void CalcBestEastWest();
  void CalcBestEast();
  void CalcBestWest();
  void CalcDt();
  Int_t StripAim();
  Float_t CurrentQ(PlexSEIdAltL* pAltList);
  Float_t AttCorWLS(Float_t lWlsInMetres);
  Float_t AttCorClear(Float_t lClearInMetres);
  Float_t SigmaDQ();
  Float_t QRawE(){return fQE;};
  Float_t QRawW(){return fQW;};
  Float_t QSigCorE(){return fQEc;};
  Float_t QSigCorW(){return fQWc;};
  Float_t QAttCorE(){return fQAttCorE;};
  Float_t QAttCorW(){return fQAttCorW;};
  Float_t PairQCor(){return fQWc+fQEc;}
  Float_t GetClearE(){return fClearE;}
  Float_t GetClearW(){return fClearW;}
  Float_t GetWlsE(){return fWlsE;}
  Float_t GetWlsW(){return fWlsW;}
  Double_t GetTE(){return fTE;}
  Double_t GetTW(){return fTW;}
  Float_t DT(){return fDt;}
  Float_t TimeWalk(float adc);
  void FindNearestEdge(Float_t x, Float_t y, Int_t &iedge, Float_t &d);

private:

  UgliGeomHandle* pUgh;
  Int_t fPlane;
  PlaneView::PlaneView_t fView;
  PlexSEIdAltL* pAltListE;
  PlexSEIdAltL* pAltListW;
  Double_t fAbsTime;
  bool fEastIsSet;
  bool fWestIsSet;
  bool fDtIsSet;
  bool fDQIsSet;
  
  Int_t fStripE;
  Int_t fStripW;
  Float_t fQE;
  Float_t fQW;
  Float_t fQEc;
  Float_t fQWc;
  Float_t fQAttCorE;
  Float_t fQAttCorW;
  Float_t fClearE;
  Float_t fClearW;
  Float_t fWlsE;
  Float_t fWlsW;
  Double_t fTE;
  Double_t fTW;
  Double_t fDt;
  Float_t fDQ;
  Int_t   fStripAim;
  Float_t fQEcor;
  Float_t fQWcor;
  Float_t fClearFibreC;
  Float_t fWLSFibreC;
  Float_t fScintillatorC;
  Float_t fDt2Dstrip;

  Int_t   fNumberOfStrips;


// Fibre lengths for the EVENT - used often so store in large arrays 
  Float_t  fWLSFibreLengthE[MAX_NUMBER_OF_PLANES][MAX_NUMBER_OF_STRIPS];
  Float_t  fWLSFibreLengthW[MAX_NUMBER_OF_PLANES][MAX_NUMBER_OF_STRIPS];
  Float_t  fClearFibreLengthE[MAX_NUMBER_OF_PLANES][MAX_NUMBER_OF_STRIPS];
  Float_t  fClearFibreLengthW[MAX_NUMBER_OF_PLANES][MAX_NUMBER_OF_STRIPS];
// Lookup tables for Attenuation corrections 
  Float_t fAttCorWLS[2500];
  Float_t fAttCorClear[2500];

  ClassDef(AltDeMuxCalc, 0);
};

#endif // ALTDEMUXCALC_H

---