BmnPairFinder.h

Go to the documentation of this file.

 
 
#ifndef BMNPAIRFINDER_H
#define BMNPAIRFINDER_H
 
#include "BmnFieldMap.h"
#include "BmnGemStripHit.h"
#include "BmnGemStripStationSet.h"
#include "BmnGemStripStationSet_RunSpring2017.h"
#include "BmnGemStripStationSet_RunWinter2016.h"
#include "BmnGlobalTrack.h"
#include "BmnKalmanFilter.h"
#include "BmnNewFieldMap.h"
#include "BmnSiliconPoint.h"
#include "BmnTrackMatch.h"
#include "BmnVertex.h"
#include "CbmMCTrack.h"
#include "CbmStsPoint.h"
#include "CbmVertex.h"
#include "FairTask.h"
 
#include <BmnEventHeader.h>
#include <BmnParticlePair.h>
#include <DstEventHeader.h>
#include <TClonesArray.h>
#include <TDatabasePDG.h>
#include <TF1.h>
#include <TFitResult.h>
#include <TGraph.h>
#include <TH1F.h>
#include <TH1I.h>
#include <TH2F.h>
#include <TLorentzVector.h>
#include <TMath.h>
#include <TString.h>
#include <TVector2.h>
#include <TVector3.h>
#include <UniDetectorParameter.h>
#include <UniRun.h>
#include <iostream>
#include <vector>
 
using namespace std;
using namespace TMath;
 
// Extenion of "standard" dst event header by info to be used in dec. particle anal ...
 
class BmnPairFinder : public FairTask
{
  private:
    BmnKalmanFilter* fKalman;
 
    DstEventHeader* fDstHeader;
    TClonesArray* fStsTracks;
    CbmVertex* fVertex;
 
    TClonesArray* fParticlePair;
 
    Double_t FindV0ByVirtualPlanes(BmnGlobalTrack* track1, BmnGlobalTrack* track2, Double_t z_0, Double_t range);
    Double_t FindInvMass(FairTrackParam pi_plus_meson_V0, FairTrackParam pi_minus_meson_V0);
    Double_t GetPath(Double_t vPx, Double_t vPy, Double_t vPz, Double_t v0x, Double_t v0y, Double_t v0z);
    Double_t FindDCA0(FairTrackParam pi_plus_meson_V0,
                      FairTrackParam pi_minus_meson_V0,
                      Double_t vPx,
                      Double_t vPy,
                      Double_t vPz,
                      Double_t v0x,
                      Double_t v0y,
                      Double_t v0z);
    Double_t FindEta(FairTrackParam pi_plus_meson_V0);
    Double_t FindAngleDecayProducts(FairTrackParam proton_V0, FairTrackParam pi_minus_meson_V0);
 
    Int_t fRunId;
    Int_t evNumber;
    unique_ptr<TParticlePDG> fParticle;
    unique_ptr<TParticlePDG> fDecayFirstParticle;
    unique_ptr<TParticlePDG> fDecaySecondParticle;
    Double_t lowV0Z;
 
  public:
    BmnPairFinder(TParticlePDG* particle = TDatabasePDG::Instance()->GetParticle("Lambda0"));
    virtual ~BmnPairFinder();
 
    virtual void Exec(Option_t* option);
    virtual InitStatus Init();
    virtual void Finish();
 
    ClassDef(BmnPairFinder, 0)
};
 
#endif