PndHypFullIdealAna.h

Go to the documentation of this file.

#ifndef PNDHYPFULLIDEALANA_H
#define PNDHYPFULLIDEALANA_H 1

#include "FairTask.h"
#include <map>
#include <string>
#include "TFile.h"
#include "TVector3.h"
#include "TLorentzVector.h"

class TClonesArray;
class TObjectArray;
class TH1F;
class TH2F;

class RhoChargedParticleSelector;
class RhoNeutralParticleSelector;
class RhoPlusParticleSelector;
class RhoMinusParticleSelector;
class RhoMassParticleSelector;
class RhoSimpleKaonSelector;
class RhoSimplePionSelector;
class RhoSimpleProtonSelector;

class PndHypFullIdealAna : public FairTask {

 public:
  typedef std::map<Int_t, Float_t> mapper;

  PndHypFullIdealAna();

  ~PndHypFullIdealAna();

  virtual InitStatus Init();

  virtual void Exec(Option_t *opt);

  void Finish(TString cat);

  void SetEnergySpectra(int event, int cluster);
  void SetTotESpectra(int clus);
  Int_t GetIonCharge(Int_t Z, Int_t &mass, Int_t &str);

  // void CreateStructure();

 protected:
  int evcount;

  // **** some basic selectors
  //
  // RhoChargedParticleSelector *chargedSel;
  RhoNeutralParticleSelector *neutralSel;
  RhoPlusParticleSelector *plusSel;
  RhoMinusParticleSelector *minusSel;

  // **** mass selectors for the resonances/composites
  //
  RhoMassParticleSelector *phiMSel;
  RhoMassParticleSelector *pi0MSel;
  RhoMassParticleSelector *dsMSel;
  RhoMassParticleSelector *LambMSel;
  RhoSimpleKaonSelector *kSel;
  RhoSimplePionSelector *piSel;
  RhoSimpleProtonSelector *pSel;

  TH2F *hvtx2[10];
  TH1F *spectra[10];

  TH1F *ds0mass;
  TH1F *ximass;
  TH1F *Lamb;
  TH1F *ppi2mass;
  TH1F *ppi2;
  TH1F *e;
  TH2F *pid;
  TH2F *pidh;

  TH1F *nmult;

 private:
  TClonesArray *fChargedArray;
  TClonesArray *fMcTr;
  TClonesArray *fMicroArray;
  TClonesArray *fMcCands;
  TClonesArray *fMc;
  TClonesArray *fGe;

  virtual void SetParContainers();

  ClassDef(PndHypFullIdealAna, 1);
};

#endif