Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
#ifndef HIPPLOTS_H
#define HIPPLOTS_H

#include <Riostream.h>
#include <vector>
#include <string>
//#include <stdio.h>
#include "TFile.h"
#include "TH2F.h"
#include "TChain.h"
#include "TString.h"
#include "TCanvas.h"
#include "TFileIter.h"
#include "TStyle.h"
#include "TMath.h"
#include "TCut.h"
#include "TGraph.h"

class HIPplots{
public:
  HIPplots(int IOV, char* path, char* outFile);
  void extractAlignParams(int i, int minHits = 0, int subDet = 0, int doubleSided = 0);
  void extractAlignShifts(int i, int minHits = 0, int subDet = 0);
  void plotAlignParams(string ShiftsOrParams, char* plotName = "test.png");
  void plotAlignParamsAtIter(int iter, string ShiftsOrParams, char* plotName = "test.png");
  void plotHitMap(char *outpath, int subDet, int minHits=0);
  void extractAlignableChiSquare(int minHits=0, int subDet =0, int doubleSided = 0);
  void plotAlignableChiSquare(char *plotName ="testchi2.png", float minChi2n=-1.0);
  void extractSurveyResiduals(int currentPar, int subDet =0);
  void dumpAlignedModules(int nhits=0);

  int _IOV;
  TString _path;
  TString _outFile;
  TString _inFile_params;
  TString _inFile_uservars;
  TString _inFile_truepos;
  TString _inFile_alipos;
  TString _inFile_mispos;
  TString _inFile_HIPalign;
  TString _inFile_surveys;
  enum TKdetector_id{ unknown=0, TPBid=1, TPEid=2, TIBid=3, TIDid=4, TOBid=5, TECid=6, ALLid=99 };
  TLegend * MakeLegend(double x1=0.1, double y1=0.1, double x2=0.1, double y2=0.1);
  int GetNIterations(TDirectory *f, char *tag, int startingcounter=0);
  int GetSubDet(unsigned int id);
  int GetBarrelLayer(unsigned int id);
  void SetMinMax(TH1* h);
  int FindPeaks(TH1F *h1, float *peaklist, const int maxNpeaks, int startbin=-1, int endbin=-1);
  void SetPeakThreshold(float newpeakthreshold);
  void CheckFiles(int &ierr);
  bool CheckFileExistence(TString filename);
  bool CheckHistoRising(TH1D *h);
  float peakthreshold;
  bool plotbadchi2;
};

#endif