HarpoAnalyseMonitorNew.h

Go to the documentation of this file.

#ifndef HARPOANALYSEMONITORNEW_H
#define HARPOANALYSEMONITORNEW_H

//
// HarpoAnalyseMonitorNew
//
// Dummy analise which print all run events.
//

#include "TObject.h"
#include "HarpoAnalyse.h"
#include "TProfile.h"
#include "TGraphErrors.h"
#include "TMapFile.h"
#include "TH1F.h"
#include "TH2F.h"

const ULong_t  gkHarpoMonNhist = 18;

typedef enum HarpoHistoIdentifiers {
  H_MAP = 0,
  H_MAP_ZS,
  H_Q,
  H_Q_EV,
  H_T,
  H_CH,
  H_QTOT_T,
  H_QTOT,
  H_Q_ZS,
  H_Q_EV_ZS,
  H_T_ZS,
  H_CH_ZS,
  H_QTOTT_ZS,
  H_QTOT_ZS,
  H_T_EV,
  H_T_EV_ZS,
  H_QTOTT_EV,
  H_QTOTT_EV_ZS
} HarpoHistoIdentifier;


class HarpoAnalyseMonitorNew : public HarpoAnalyse {

 public:
  
  HarpoAnalyseMonitorNew() {}

  void   Init();
  void   Save(char *mode=NULL);
  void  process(); 
  void print(); 

  TH1* GetHist(HarpoHistoIdentifiers i, Int_t ndet);

  //  HarpoMap *m_hmap[2];
 private:

  Int_t fQmin;     // Qtot cut for empty events
  Int_t fQmax;     // Qtot cut for saturated events
  Double_t fTmin;  // cut for beginning of Tpc
  Double_t fLTpc;  // Tmax-Tmin cut for traversing tracks

  TH1F* hTdrift_TT[2];
  TH1F* hQvsT_TT[2];
  TH2F* hQvsT_TT2D[2];

  TH1F* hQ[2];
  TH1F* hQev[2];
  TH1F* hT[2];
  TH1F* hCh[2];
  TH1F* hQtotT[2];
  TH1F* hQtot[2];

  TH1F* hQZs[2];
  TH1F* hQevZs[2];
  TH1F* hTZs[2];
  TH1F* hChZs[2];
  TH1F* hQtotTZs[2];
  TH1F* hQtotZs[2];

  TH1F* hTev[2];
  TH1F* hTevZs[2];
  TH1F* hQtotTev[2];
  TH1F* hQtotTevZs[2];

  TH1F* hTmin[2];
  TH1F* hTmax[2];

  TH2F* h[2];
  TH2F* hZs[2];
  TH2F* hCumul[2];

  TH1F* hInfo;
  TH1F* hPmm2;

  TH1F* hPmm2Size;
  TH1F* hPmm2Missed;
  TH1F* hPmm2Triggered;
  TH2F* hPmm2Spectrum2D;
  TH2F* hPmm2Coincidence;

  TH1F* hTriggerCorrection;
  TH1F* hTriggerCreated;
  TH1F* hTriggerEmpty;
  TH1F* hTriggerTraversing;
  TH1F* hTriggerSaturated;
  TH1F* hTriggerOther;

  static const Int_t kNtriggerTypes = 20;

  TGraphErrors* gRateCreated[kNtriggerTypes+1];
  TGraphErrors* gRateEmpty[kNtriggerTypes+1];
  TGraphErrors* gRateTraversing[kNtriggerTypes+1];
  TGraphErrors* gRateSaturated[kNtriggerTypes+1];
  TGraphErrors* gRateOther[kNtriggerTypes+1];
  TGraphErrors* gRateTotal[kNtriggerTypes+1];
  Int_t fNcreated[kNtriggerTypes+1];
  Int_t fNempty[kNtriggerTypes+1];
  Int_t fNtraversing[kNtriggerTypes+1];
  Int_t fNsaturated[kNtriggerTypes+1];
  Int_t fNother[kNtriggerTypes+1];
  TH1F* hDeadTime[kNtriggerTypes+1];

  TGraph* gXTmin[2];

  TGraphErrors* gQtotTT_X;
  TGraphErrors* gQtotTT_Y;
  Double_t qtotX;
  Double_t qtotY;
  Double_t qtot2X;
  Double_t qtot2Y;
  Int_t nSamples;

  Double_t timestampold;
  Double_t timestampprev;
  Double_t fRateSample;
  Int_t fSaveRate;
  Double_t fTriggerTimeTot;
  Int_t fTriggerTime[kNtriggerTypes];
  Int_t fTriggerDownscaling[kNtriggerTypes];

  Int_t fIsZs;

  TMapFile* mfile;

  Double_t TruncMean(TArrayD* vect, Double_t tl, Double_t th);
 
  ClassDef(HarpoAnalyseMonitorNew,1)  //Run AnalysePrint
};

#endif