CbmEvBuildSource.h

Go to the documentation of this file.

// FairSource which reads a time slice and creates events using a simple approach
// Use the same approach as CbmBuildEventSimple
// //Dr.Sys : mikhail.prokudin@cern.ch
 
#ifndef CBMEVBUILDSOURCE_H
#define CBMEVBUILDSOURCE_H
 
#include "FairSource.h"
 
#include <iostream>
 
class TChain;
class TClonesArray;
class CbmTimeSlice;
class FairEventHeader;
 
class CbmEvBuildSource : public FairSource {
public:
  CbmEvBuildSource();
  CbmEvBuildSource(const char* fname);
  virtual ~CbmEvBuildSource();
  void AddFile(const char* fname);
 
  virtual Bool_t Init();
  virtual Int_t ReadEvent(UInt_t par = 0);
  virtual void Close();
 
  virtual void Reset();
 
  //virtual Source_Type GetSourceType() {return kFILE;}
 
  virtual Bool_t ActivateObject(TObject**, const char*) { return kFALSE; }
 
  virtual Int_t CheckMaxEventNo(Int_t = 0) { return 1000000000; }
  virtual void ReadBranchEvent(const char*) { return; }
  virtual void ReadBranchEvent(const char*, Int_t) { return; }
  virtual void FillEventHeader(FairEventHeader* feh);
 
  // Setters. Time in [ns], events in counts :) Only STS information is used.
  void SetWindDur(Double_t winddur = 2) { fWindDur = winddur; }
  void SetMinusDeltaT(Double_t minusdeltat = 2) { fMinusDeltaT = minusdeltat; }
  void SetPlusDeltaT(Double_t plusdeltat = 5) { fPlusDeltaT = plusdeltat; }
  void SetMinHitStations(Int_t minhitstations = 8) {
    fMinHitStations = minhitstations;
  }
  void SetMinDigis(Int_t mindigis = 2000) { fMinDigis = mindigis; }
  void SetDeadT(Double_t deadt = 10) { fDeadT = deadt; }
 
public:
  CbmEvBuildSource(const CbmEvBuildSource& source);
  void FillEvent(Int_t st, Int_t end);  //Fill output TClonesArray
  Int_t GetNextTimeSlice();             //Load next time slice from disk
 
private:
  TChain* fCh;             
  Int_t fN;                //Number of entries in input tree
  Int_t fI;                //Current entry in tree
  TClonesArray* fSTSDigi;  
  Int_t fNSTSDigis;        //Number of STS digis in current event
  CbmTimeSlice* fSlice;    
  Int_t fISts;             //Number of STS digi under consideration
 
  UInt_t fEv;            // Number of current event **/
  Int_t fNDigis;         // Number of digis in current time window
  Int_t fSN;             // N of start of time window
  Double_t fST;          // Time of start of time window
  Int_t fNStsDigis[16];  // Number of STS digis per station
  Double_t fWindDur;     // Duration of time window [ns]
  Double_t
    fMinusDeltaT;  // [ns] before start of time window for event construction
  Double_t fPlusDeltaT;  // [ns] after end of time window for event construction
  Double_t fDeadT;       // [ns] after last accepted digi of blindless
  Int_t
    fMinHitStations;  // Minimum number of hit STS stations for event formation
  Int_t
    fMinDigis;  // Minimum number of digis in time window for event formation
  FairEventHeader* fEvHeader;  
 
  CbmEvBuildSource& operator=(const CbmEvBuildSource&);
 
  ClassDef(CbmEvBuildSource, 1)
};
 
 
#endif