HarpoAnalyseBasic.h

Go to the documentation of this file.

#ifndef HARPOANALYSEBASIC_H
#define HARPOANALYSEBASIC_H

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

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

class HarpoAnalyseBasic : public HarpoAnalyse {

 public:
  
  HarpoAnalyseBasic() {}

  void   Init();
  void   Save(char *mode=NULL);
  void   process(); 
  void   print(); 
  void   FillHistograms();
  void   FillHistogramsPMM2();
  void   FillHistogramsTrigger();
 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

  Double_t fNormTTz;
  Double_t fNormTTzX;
  Double_t fNormTTzY;

  //position of signal
  TH1F *h1_tmin;
  TH1F *h1_tmax;
  TH1F *h1_xmin;
  TH1F *h1_xmax;
  TH1F *h1_ymin;
  TH1F *h1_ymax;
  TH1F *h1_dtime;


  //
  TH1F *h1_Qpix_x;
  TH1F *h1_Q_time_x;
  TH1F *h1_Q_channel_x;

  TH1F *h1_Qpix_y;
  TH1F *h1_Q_time_y;
  TH1F *h1_Q_channel_y;

  //clustering
  TH1F *h1_QTcluster_x;
  TH1F *h1_QCcluster_x;
  TH1F *h1_SigTcluster_x;
  TH1F *h1_SigCcluster_x;
  TH1F *h1_WidthTcluster_x;
  TH1F *h1_WidthCcluster_x;
  TH1F *h1_IndexTcluster_x;
  TH1F *h1_IndexCcluster_x;
  TH1F *h1_MeanTcluster_x;
  TH1F *h1_MeanCcluster_x;

  TH1F *h1_QTcluster_y;
  TH1F *h1_QCcluster_y;
  TH1F *h1_SigTcluster_y;
  TH1F *h1_SigCcluster_y;
  TH1F *h1_WidthTcluster_y;
  TH1F *h1_WidthCcluster_y;
  TH1F *h1_IndexTcluster_y;
  TH1F *h1_IndexCcluster_y;
  TH1F *h1_MeanTcluster_y;
  TH1F *h1_MeanCcluster_y;

  //Tracking
  TH1F *h1_theta_x;
  TH1F *h1_rho_x;
  TH1F *h1_sigtheta_x;
  TH1F *h1_sigrho_x;
  TH1F *h1_QcAssoc_x;
  TH1F *h1_QcNAssoc_x;
  TH1F *h1_QtAssoc_x;
  TH1F *h1_QtNAssoc_x;
  TH1F *h1_TypeTrack_x;
  TH1F *h1_angleTrack_x;

  TH1F *h1_theta_y;
  TH1F *h1_rho_y;
  TH1F *h1_sigtheta_y;
  TH1F *h1_sigrho_y;
  TH1F *h1_QcAssoc_y;
  TH1F *h1_QcNAssoc_y;
  TH1F *h1_QtAssoc_y;
  TH1F *h1_QtNAssoc_y;
  TH1F *h1_TypeTrack_y;
  TH1F *h1_angleTrack_y;

  //distortion
  Double_t fXT, fYT, fXC, fYC, fXdT, fXdC,fYdT,fYdC;

  TH2F *h2_distortion_dc_t_x;
  TH2F *h2_distortion_dc_c_x;
  TH2F *h2_distortion_dt_t_x;
  TH2F *h2_distortion_dt_c_x;
  TProfile *h2_prof_dc_t_x;
  TProfile *h2_prof_dc_c_x;
  TProfile *h2_prof_dt_t_x;
  TProfile *h2_prof_dt_c_x;
  TH1F *h1_distortion_dt_x;
  TH1F *h1_distortion_dc_x;

  TH2F *h2_distortion_dc_t_y;
  TH2F *h2_distortion_dc_c_y;
  TH2F *h2_distortion_dt_t_y;
  TH2F *h2_distortion_dt_c_y;
  TProfile *h2_prof_dc_t_y;
  TProfile *h2_prof_dc_c_y;
  TProfile *h2_prof_dt_t_y;
  TProfile *h2_prof_dt_c_y;
  TH1F *h1_distortion_dt_y;
  TH1F *h1_distortion_dc_y;

  //
  TH1F *h1_Qassoc_time_x;
  TH1F *h1_QtAssocN_x;
  TH1F *h1_QtAssocN_time_x;
  TH1F *h1_QtAssocN_channel_x;
  
  TH1F *h1_Qassoc_time_y;
  TH1F *h1_QtAssocN_y;
  TH1F *h1_QtAssocN_time_y;
  TH1F *h1_QtAssocN_channel_y;
  //
  TH2F *h2_QtAssocN_time_x;
  TH2F *h2_QtAssocN_time_y;
  TProfile *h2_Qprof_x;
  TProfile *h2_Qprof_y;
 

  //pmm2
  TH1F * helen;   
  TH1F * heprof;  
  TH1F * hadc[2]; 
  TH2F * htvse;  
  TH2F * hpps;    
  TH2F * hramp;   
  TH1F * htdiff;  
  TH1F * htdiff2;  

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

  TH1F* hTriggerCorrection;

  TH1F* hTriggerCreated;
  TH1F* hTriggerEmpty;
  TH1F* hTriggerTraversing;
  TH1F* hTriggerShower;
  TH1F* hTriggerOther;
  TH1F* hTriggerTTZ;

  TH1F* hTriggerContained;
  TH1F* hTriggerDepot;
  TH1F* hTriggerTwoTrGamma;
  TH1F* hTriggerTwoTrDelta;
  TH1F* hTriggerTwoTrOther;

  static const Int_t kNtriggerTypes = 16;

  TGraphErrors* gRateCreated[kNtriggerTypes+1];
  TGraphErrors* gRateEmpty[kNtriggerTypes+1];
  TGraphErrors* gRateTraversing[kNtriggerTypes+1];
  TGraphErrors* gRateSaturated[kNtriggerTypes+1];
  TGraphErrors* gRateOther[kNtriggerTypes+1];
  TGraphErrors* gRateTotal[kNtriggerTypes+1];

  TGraphErrors* gRateContained[kNtriggerTypes+1];
  TGraphErrors* gRateDepot[kNtriggerTypes+1];
  TGraphErrors* gRateTwoTrGamma[kNtriggerTypes+1];
  TGraphErrors* gRateTwoTrDelta[kNtriggerTypes+1];
  TGraphErrors* gRateTwoTrOther[kNtriggerTypes+1];
  
  Int_t fNcreated[kNtriggerTypes+1];
  Int_t fNempty[kNtriggerTypes+1];
  Int_t fNtraversing[kNtriggerTypes+1];
  Int_t fNshower[kNtriggerTypes+1];
  Int_t fNother[kNtriggerTypes+1];

  Int_t fNcontained[kNtriggerTypes+1];
  Int_t fNdepot[kNtriggerTypes+1];
  Int_t fNtwotrgamma[kNtriggerTypes+1];
  Int_t fNtwotrdelta[kNtriggerTypes+1];
  Int_t fNtwotrother[kNtriggerTypes+1];
  Int_t f1TZ[kNtriggerTypes+1];

  Int_t fTriggerTime[kNtriggerTypes];
  Int_t fTriggerDownscaling[kNtriggerTypes];
  TGraph* alignmentx;
  TGraph* alignmenty;
/* TH1F* alignmentx; */
/* TH1F* alignmenty; */
  TH1F* hTriggerType;
  int itsold;   // old time stamp value
  Double_t TruncMean(TArrayD* vect, Double_t tl, Double_t th);

  //vertexing
  TH1F* h1Doca;
  TH1F* h1Poca_x;
  TH1F* h1Poca_y;
  TH1F* h1Poca_z;
  TH2F* h2Poca_xy;
  TH2F* h2Poca_xz;
  TH2F* h2Poca_yz;

  TH1F* h1Px;
  TH1F* h1Py;
  TH1F* h1Pz;
  TH2F* h2Pxy;
  TH2F* h2Pxz;
  TH2F* h2Pyz;

  TH1F* h1Azimutal;
  TH1F* h1AOpen;
  //alignement
  TH1F* h1xd;
  TH1F* h1yd;
  TH1F* h1zxd;
  TH1F* h1zyd;
  TH2F* h2xz;
  TH2F* h2yz;
  TH2F* h2xzn;
  TH2F* h2yzn;
  TH2F* h2xyn;

  //alignement first 3 clusters
  TH1F* h1x3d;
  TH1F* h1y3d;
  TH1F* h1zx3d;
  TH1F* h1zy3d;
  TH2F* h2xz3n;
  TH2F* h2yz3n;
  TH2F* h2xy3n;

  //cluster <---> track

  TH1F* h1NTrackX;
 TH1F* h1NTrackXTest;
  TH1F* h1NTrackY; 
  TH1F* h1NClusterX;
  TH1F* h1NClusterY;
  TH1F* h1NClusterNX;
  TH1F* h1NClusterNY;
  TH1F *h1_pro_timeX;
  TH1F *h1_pro_timeY;

  TH1F *h1_NEvent;
  TH1F *h1_TimeSt;

  TH1F *h1_vertex_tr1;
  TH1F *h1_vertex_tr2;
  TH2F *h2_vertex;

  Int_t NEventEmpty;
  Int_t NEventDEPOSIT;
  Int_t NEventTTZ;
  Int_t NEventTT;
  Int_t NEventTCR;
  Int_t NEventTCN;
  Int_t NEventGAMMA;
  Int_t NEventDELTA;
  Int_t NEvent2TOTHER;
  Int_t NEventSHOWER;
  Int_t NEventOTHER;
  Int_t NEvent;
  Double_t TimeSt;

  //ULong_t TimeSt;



  TH1F* h1_sigrho_gamma_x;
  TH1F* h1_sigrho_gamma_y;
  TH1F* h1_sigrho_TTZ_x;
  TH1F* h1_sigrho_TTZ_y;

  Double_t fAlignmentX, fAlignmentY, fAlignmentZx,fAlignmentZy, fAlignmentZ;
  Int_t fNe;
  Double_t fDoca,fPoca_x,fPoca_y,fPoca_z,fAngleOpen, fAngleAzimuth,fsigrho;
  Int_t fEventType;
  TTree *fTree;
  Int_t ftriggerType;
  //TMapFile* mfile;

  ClassDef(HarpoAnalyseBasic,1)  //Run AnalysePrint
    };

#endif