LIPlexMaps.h

Go to the documentation of this file.

// Program name: LIPlexMaps.h
//                                                                     
// Package: LISummary                                         
//                                                                    
// Coded by Jeff Hartnell Oct/2002-Sep/2003          
//                                                                    
// Purpose: To create wiring maps of LI cabling using the plex 
//                                                                    
// Contact: jeffrey.hartnell@physics.ox.ac.uk                         

#ifndef LIPLEXMAPS_H
#define LIPLEXMAPS_H

#include <vector>
#include <string>
#include <map>

#include "TStyle.h"
#include "TRandom.h"
#include "TGraph.h"

#include "DatabaseInterface/DbiResultPtr.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexPinDiodeToLed.h"
#include "Plex/PlexStripEndToLed.h"
#include "PulserCalibration/PulserRawDrift.h"
#include "PulserCalibration/PulserRawDriftF.h"
#include "Validity/VldContext.h"
#include "Validity/VldRange.h"

#include "LISummary/LILookup.h"
#include "LISummary/LIPlane.h"

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

class LIPulserPin
{
  
  //a little class to hold the pulser raw drift pin table info

public:
  
  LIPulserPin() { 
    mean=-1;
    numEnt=-1;
    error=-1;
  };

  Float_t mean;
  Float_t numEnt;
  Float_t error;

};

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

class myCalStripToStrip
{
  //a little class to hold the CalStripToStrip objects data since they
  //can't be copy constructed

public:
  
  myCalStripToStrip() { 
    seid_enc=-1;
    response=-1;
    responseErr=-1;
  };

  myCalStripToStrip(Int_t s,Float_t r,Float_t rErr) { 
    seid_enc=s;
    response=r;
    responseErr=rErr;
  };

  Int_t seid_enc;
  Float_t response;
  Float_t responseErr;
};

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

class LIPlexMaps : public LILookup
{

 public:  
  
  LIPlexMaps();
  ~LIPlexMaps();

  void CompareStripToStrip(Int_t det,VldTimeStamp ts,Int_t taskA,
                           Int_t taskB);
  void DriftWithTime(VldRange vldRng,Float_t pointSpacingMinutes=20);
  std::map<Int_t,Float_t> GetRawPinMap(VldContext vldCtx);
  std::map<PlexLedId,LIPulserPin> GetRawPinLedMap(VldContext vldCtx); 
  std::map<Int_t,LIPulserPin> GetRefAdcPinRatio(VldContext vldCtx);
  //map<Int_t,LIPulserPin> GetPulserRawDrift
  //(DbiResultPtr<PulserRawDrift> dbiPulserRawDrift,VldContext vldCtx);
  void DumpCalStripToStrip(std::string sDataFileName,
                           Detector::Detector_t det,
                           SimFlag::ESimFlag simFlag,
                           VldTimeStamp ts,Int_t task) const;
  void DumpCalStripAtten(std::string sDataFileName,
                         Detector::Detector_t det,
                         SimFlag::ESimFlag simFlag,
                         VldTimeStamp ts,Int_t task);

//  void DumpGains(Detector::Detector_t det, VldTimeStamp ts,
//               std::string gaindump_filename);

  Bool_t IsPermitted() const;

  void MakeCalAdcToPe(std::string gainsFileName,
                      Detector::Detector_t det,
                      SimFlag::ESimFlag simFlag,
                      VldTimeStamp vldStart,VldTimeStamp vldEnd,
                      Int_t aggNo,Int_t task,
                      Int_t aggStat[484],
                      VldTimeStamp creationDate,std::string dbName,
                      std::string sLogComment,
                      Int_t gainTooHigh,Int_t gainCeiling,
                      Bool_t writeToDb=false);
  void MakeCalStripToStrip(std::string sDataFileName,
                           Detector::Detector_t det,
                           SimFlag::ESimFlag simFlag,
                           VldTimeStamp vldStart,VldTimeStamp vldEnd,
                           Int_t aggNo,Int_t task,
                           VldTimeStamp creationDate,std::string dbName,
                           std::string sLogComment,
                           Bool_t writeToDb=false);
  void MakeCalStripAtten(std::string sDataFileName,
                         Detector::Detector_t det,
                         SimFlag::ESimFlag simFlag,
                         VldTimeStamp vldStart,VldTimeStamp vldEnd,
                         Int_t aggNo,Int_t task,
                         VldTimeStamp creationDate,std::string dbName,
                         std::string sLongComment,
                         Int_t iFileFormat,
                         Bool_t writeToDb=false);
  void MakeCalMIPCalibration(Float_t scale,
                             Detector::Detector_t det,
                             SimFlag::ESimFlag simFlag,
                             VldTimeStamp vldStart,
                             VldTimeStamp vldEnd,
                             Int_t aggNo,Int_t task,
                             VldTimeStamp creationDate,
                             std::string dbName,std::string sLogComment,
                             Bool_t writeToDb);
  void MakePinPlexTables();
  void MakePlexStripEndToLed(Bool_t writeToDb=false); 
  void NormaliseVector(std::vector<Double_t>& v,Int_t mode=1);
  void PlotPlexPinMap();
  void PlotPlexStripMap(Int_t detectorType);
  void PrintPlexPinDiodeInfo(Int_t detectorType);
  void PrintDbPinDiodeInfo();
  void RatioRawToRefDrift();
  void ReadDbCalStripToStrip(Int_t det,VldTimeStamp ts,Int_t task);
  void ReadDbCalStripAtten(Int_t det,VldTimeStamp ts, Int_t task);
  void ReadDbGains(Int_t det,SimFlag::ESimFlag simFlag,VldTimeStamp ts,std::string runNumber,std::string sAgg,Int_t task=0);
  // Remake this as a less bloated function, to grab the detector gains avearge
  Double_t ReadDbGainAverage(Detector::Detector_t det, VldTimeStamp ts, VldTimeStamp *end = 0, VldTimeStamp *start = 0);
  TH1D *ReadDbGains(Detector::Detector_t det, VldTimeStamp ts);
  std::map<Int_t,Double_t> GetGainMap(Detector::Detector_t det, VldTimeStamp ts);
  TH1D GainDifference(Detector::Detector_t det, VldTimeStamp t1, VldTimeStamp t2);
  TH2D GainDifferenceSvP(Detector::Detector_t det, VldTimeStamp t1, VldTimeStamp t2);
  void DumpGains(Detector::Detector_t det, VldTimeStamp ts, std::string filename);
  
  void ReadDbPulserDrift(Int_t det,VldTimeStamp ts);
  void ReadDbPulserDriftPin();
  void ReadDbPulserRawGain(); 
  void SetGraphAxisEtc(TGraph* g,Int_t startTimeSecs,Int_t endTimeSecs);

  //MC methods
  void ADCRailEffect();
  void Attenuation();
  void BinomialGen(std::vector<Double_t>& v,Int_t N,Double_t p);
  void FactorialGen(std::vector<Double_t>& v);
  void LIErrors();
  void NonLinearity();
  void PeGen();
  Int_t SampleDist(std::vector<Double_t> v,TRandom* tRand);
  void ScaleVector(std::vector<Double_t>& v,Double_t scaleFactor);

 private:

  std::string GetElecString(RawChannelId rChId);
  std::string fS;//a general purpose string
  Bool_t IsScintStrip(PlexStripEndId& stripEndId,
                      PlexHandle* plexHandle);

  //data members:
  LILookup fLookup;

  ClassDef(LIPlexMaps,0)
};

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

#endif // LIPLEXMAPS_H

---