CbmAnaConversionTomography.h

Go to the documentation of this file.

/* Copyright (C) 2014-2016 Fakultaet fuer Mathematik und Naturwissenschaften, Bergische Universitaet Wuppertal, Wuppertal
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sascha Reinecke [committer], Florian Uhlig */
 
#ifndef CBM_ANA_CONVERSION_TOMOGRAPHY
#define CBM_ANA_CONVERSION_TOMOGRAPHY
 
// 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 <vector>
 
class CbmAnaConversionTomography {
 
public:
  CbmAnaConversionTomography();
  virtual ~CbmAnaConversionTomography();
 
  void Init();
  void InitHistos();
  void Finish();
 
  void Exec();
 
  void TomographyMC(int electronID);
  void TomographyReco(CbmMCTrack* mctrack);
  Bool_t GetNPoints(CbmMCTrack* mctrack);
 
 
private:
  TClonesArray* fMcTracks;
 
  std::vector<TH1*> fHistoList_tomography;  // list of all histograms
 
 
  TH1D* fhGammaZ;
  TH3D* fhTomography;
  TH2D* fhTomography_XZ;
  TH2D* fhTomography_YZ;
  //TH2D * fhTomography_XZ_cut;
  //TH2D * fhTomography_YZ_cut;
  TH2D* fhTomography_uptoRICH;
  TH2D* fhTomography_RICH_complete;
  TH2D* fhTomography_RICH_beampipe;
  TH2D* fhTomography_STS_end;
  TH2D* fhTomography_STS_lastStation;
  TH2D* fhTomography_RICH_frontplate;
  TH2D* fhTomography_RICH_backplate;
 
  TH1D* fhConversion;
  TH1D* fhConversion_cut;
  TH1D* fhConversion_inSTS;
  TH1D* fhConversion_prob;
  TH1D* fhConversion_energy;
  TH1D* fhConversion_p;
  TH2D* fhConversion_vs_momentum;
 
 
  TH3D* fhTomography_reco;
  TH2D* fhTomography_reco_XZ;
  TH2D* fhTomography_reco_YZ;
  TH1D* fhConversion_reco;
 
  std::vector<int> electronIDs;
  std::vector<int> electronMotherIDs;
  Int_t conversionsInDetector[5];  // 0 = magnet, 1 = sts, 2 = rich, 3 = trd, 4 = tof
  //Int_t conversionsInDetector_cut[5];         // 0 = magnet, 1 = sts, 2 = rich, 3 = trd, 4 = tof
  TH1I* fhConversionsPerDetector;
  TH2I* fhConversionsPerDetectorPE;
  //TH2I * fhConversionsPerDetectorPE_cut;
 
  // timer
  TStopwatch timer;
  Double_t fTime;
 
  CbmAnaConversionTomography(const CbmAnaConversionTomography&);
  CbmAnaConversionTomography operator=(const CbmAnaConversionTomography&);
 
  ClassDef(CbmAnaConversionTomography, 1)
};
 
#endif