FemtoBasicAna.h

/*
 Be * FemtoBasicAna2.h
 *
 *  Created on: 26-11-2013
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
 
#ifndef HALFEMTOBASICANA_H_
#define HALFEMTOBASICANA_H_
 
#include "Array.h"
#include "FemtoCorrFunc.h"
#include "FemtoFreezoutGenerator.h"
#include "FemtoWeightGenerator.h"
#include "TwoTrackAna.h"
 
#include <Rtypes.h>
#include <RtypesCore.h>
 
 
namespace Hal {
  class FemtoBasicAna : public TwoTrackAna {
  protected:
    Int_t fPdg1, fPdg2;
    Bool_t fUseImgMomenta;
    Bool_t fIgnoreSign;
    Double_t fFsiWeight;
    FemtoPair* fFemtoPair;
    FemtoWeightGenerator* fCalc;
    FemtoFreezoutGenerator* fFreezoutGenerator;
    ObjectMatrix_2* fCFs;
    FemtoCorrFunc* fCFTemp;
    virtual void CheckCutContainerCollections();
    virtual void PreprocessFemtoPair();
    virtual void PreprocessFemtoPair_Perfect() {};
    virtual void PreprocessFemtoPair_Rotated() { PreprocessFemtoPair(); };
    virtual void PreprocessFemtoPair_Hemisphere() { PreprocessFemtoPair(); };
    virtual void PreprocessFemtoPair_Mixed() {};
    virtual void PreprocessFemtoPair_Charged() { PreprocessFemtoPair(); };
    virtual void ProcessFemtoPair();
    virtual void ProcessFemtoPair_Perfect();
    virtual void ProcessFemtoPair_Rotated();
    virtual void ProcessFemtoPair_Hemisphere();
    virtual void ProcessFemtoPair_Mixed();
    virtual void ProcessFemtoPair_Charged();
    // converting to Pair, using preprocess method to compute stuff like
    // freezout coordinates
    virtual void ProcessPair();
    virtual void ProcessPair2() { ProcessPair(); };
    virtual void ProcessPair_Perfect();
    virtual void ProcessPair2_Perfect() { ProcessPair_Perfect(); };
    virtual void ProcessPair_Mixed();
    virtual void ProcessPair2_Mixed() { ProcessPair_Mixed(); };
    virtual void ProcessPair_Rotated();
    virtual void ProcessPair2_Rotated() { ProcessPair_Rotated(); };
    virtual void ProcessPair_Hemisphere();
    virtual void ProcessPair2_Hemisphere() { ProcessPair_Hemisphere(); };
    virtual void ProcessPair_ChargedId();
    virtual void ProcessPair_ChargedNId() { ProcessPair_ChargedId(); };
    virtual Package* Report() const;
    virtual void ProcessEvent();
    virtual Bool_t InitArray();
    virtual Task::EInitFlag Init();
 
  public:
    FemtoBasicAna();
    FemtoBasicAna(const FemtoBasicAna& ana);
    FemtoBasicAna& operator=(const FemtoBasicAna& other);
    virtual void FinishTask();
#if __cplusplus >= 201402L
    [[deprecated]]
#endif
    virtual void
    SetCorrFctn(FemtoCorrFunc* h);
    virtual void SetCorrFctn(const FemtoCorrFunc& h);
    void SetWeight(const FemtoWeightGenerator& gen) { fCalc = gen.MakeCopy(); }
    void SetFreezoutGenerator(const FemtoFreezoutGenerator& gen) { fFreezoutGenerator = gen.MakeCopy(); }
    void SetPdg(Int_t PID) { fPdg1 = fPdg2 = PID; }
    void SetPdg(Int_t PID1, Int_t PID2) {
      fPdg1 = PID1;
      fPdg2 = PID2;
    };
    void IgnoreSign() { SetOption("ignore_sign"); };
    void UseImgMomenta() { SetOption("use_im_momenta"); };
    virtual void SetOption(Option_t* option);
    virtual ~FemtoBasicAna();
    ClassDef(FemtoBasicAna, 1)
  };
}  // namespace Hal
#endif /* HALFEMTOBASICANA_H_ */