CbmAnaConversionRich.h

Go to the documentation of this file.

 
#ifndef CBM_ANA_CONVERSION_RICH
#define CBM_ANA_CONVERSION_RICH
 
// included from ROOT
#include <TClonesArray.h>
#include <TH1.h>
#include <TH2.h>
#include <TH3.h>
#include <TStopwatch.h>
 
// included from CbmRoot
#include "CbmMCTrack.h"
#include "CbmVertex.h"
 
#include <vector>
 
class CbmAnaConversionRich {
 
public:
  CbmAnaConversionRich();
  virtual ~CbmAnaConversionRich();
 
  void Init();
  void InitHistos();
  void Finish();
 
  void AnalyseRICHdata();
  void FillAdditionalPDGhisto(Int_t pdg, Int_t motherpdg);
  void Protons(CbmMCTrack* mcTrack);
 
  void CheckMC();
 
 
private:
  TClonesArray* fRichPoints;
  TClonesArray* fRichRings;
  TClonesArray* fRichRingMatches;
  TClonesArray* fMcTracks;
  TClonesArray* fStsTracks;
  TClonesArray* fStsTrackMatches;
  TClonesArray* fGlobalTracks;
  CbmVertex* fPrimVertex;
 
  std::vector<TH1*>
    fHistoList_richrings;  // list of all histograms related to rich rings
 
 
  TH2D* fTest;
 
  // histograms containing data about rich rings
  TH1D* fRichRings_nofRings;
  TH1D* fRichRings_motherpdg;
  TH1D* fRichRings_richpdg;
  TH1D* fRichRings_electronspE;
  TH1D* fRichRings_sourcePI0;
  TH1D* fRichRings_test;
 
  TH1D* fRichRings_detectedParticles;
  TH1D* fRichRings_detParticlesMother;
 
  TH1D* fRichRings_Aaxis;
  TH1D* fRichRings_Aaxis_part1;
  TH1D* fRichRings_Aaxis_part2;
  TH1D* fRichRings_Aaxis_part3;
  TH1D* fRichRings_Aaxis_electrons;
  TH1D* fRichRings_Baxis;
  TH1D* fRichRings_Baxis_part1;
  TH1D* fRichRings_Baxis_part2;
  TH1D* fRichRings_Baxis_part3;
  TH1D* fRichRings_Baxis_electrons;
  TH1D* fRichRings_radius;
  TH1D* fRichRings_radius_electrons;
  TH2D* fRichRings_radius_vs_momentum;
  TH2D* fRichRings_radius_vs_pt;
  TH1D* fRichRings_distance;
  TH1D* fRichRings_distance_electrons;
  TH2D* fhRingtest;
  TH2D* fhRichRings_AaxisVSmom;
  TH2D* fhRichRings_BaxisVSmom;
 
  TH1D* fhRichRings_test1;
  TH1D* fhRichRings_test2;
  TH1D* fhRichRings_test3;
  TH1D* fhRichRings_test4;
  TH1D* fhRichRings_test5;
  TH1D* fhRichRings_test6;
  TH2D* fhRichRings_pos;
 
 
  TH2D* fhRichRings_protons;
  TH1D* fhRichRings_protons2;
 
  TH1D* fhRichRings_start;
 
  // timer
  TStopwatch timer;
  Double_t fTime;
 
  CbmAnaConversionRich(const CbmAnaConversionRich&);
  CbmAnaConversionRich operator=(const CbmAnaConversionRich&);
 
  ClassDef(CbmAnaConversionRich, 1)
};
 
#endif