PndProdAnaTask Class Reference

#include <PndProdAnaTask.h>

Inheritance diagram for PndProdAnaTask:

Public Member Functions
	PndProdAnaTask (int mode=0, TString pidAlgo="")
	~PndProdAnaTask ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()

Private Member Functions
int	SelectTruePid (PndAnalysis *ana, RhoCandList &l)
	ClassDef (PndProdAnaTask, 1)

Private Attributes
int	fEvtCount
int	fMode
TString	fPidAlgo
RhoTuple *	ntp
TLorentzVector	fIni
PndAnalysis *	fAnalysis

Detailed Description

Definition at line 21 of file PndProdAnaTask.h.

Constructor & Destructor Documentation

PndProdAnaTask::PndProdAnaTask	(	int	mode = 0,
		TString	pidAlgo = ""
	)

Definition at line 54 of file PndProdAnaTask.cxx.

References fMode, fPidAlgo, and mode.

                                                        :
  FairTask("Panda Tutorial Analysis Task")
{
        fMode = mode;
        fPidAlgo = pidAlgo;
        if (fPidAlgo=="") fPidAlgo = "PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc;PidAlgoStt;PidAlgoMdtHardCuts;PidAlgoSciT;PidAlgoRich";
}

PndProdAnaTask::~PndProdAnaTask

(

)

Definition at line 65 of file PndProdAnaTask.cxx.

{ }

Member Function Documentation

PndProdAnaTask::ClassDef ( PndProdAnaTask  , 1    )

private

void PndProdAnaTask::Exec ( Option_t *  opt )

virtual

Definition at line 117 of file PndProdAnaTask.cxx.

References RhoTuple::Column(), RhoTuple::DumpData(), fAnalysis, fEvtCount, PndAnalysis::FillList(), fIni, RhoFitterBase::Fit(), fMode, fPidAlgo, RhoFitterBase::GetChi2(), PndAnalysis::GetEventInTask(), RhoCandidate::GetFit(), RhoCandidate::GetMcTruth(), i, ntp, and RhoCandidate::P4().

{
        // *** some variables
        int i=0,j=0, k=0, l=0;

        // necessary to read the next event
        fAnalysis->GetEventInTask();

        if (!(++fEvtCount%100)) cout << "[PndProdAnaTask] evt "<<fEvtCount<<endl;

        // *** QA tool for simple dumping of analysis results in RhoRuple
        // *** WIKI: https://panda-wiki.gsi.de/foswiki/bin/view/Computing/PandaRootAnalysisJuly13#PndRhoTupleQA
        PndRhoTupleQA qa(fAnalysis,fIni.P());

        // *** RhoCandLists for the analysis
        RhoCandList kp, km, pip, pim, D0, D0bar;

        // *** Select with no PID info ('All'); type and mass are set
        fAnalysis->FillList(kp,  "KaonAllPlus",   fPidAlgo);
        fAnalysis->FillList(km,  "MuonAllMinus",  fPidAlgo);
        fAnalysis->FillList(pip, "PionAllPlus",   fPidAlgo);
        fAnalysis->FillList(pim, "PionAllMinus",  fPidAlgo);

        D0.Combine(km, pip, 421);
        D0bar.Combine(kp, pim, -421);

        D0.Append(D0bar);

        for (j=0; j<D0.GetLength(); ++j)
        {
                Float_t mmiss = (fIni-(D0[j]->P4())).M();
                Float_t msum  = D0[j]->M() + mmiss;

                // dump out information
                ntp->Column("ev",               (Int_t) fEvtCount);
                ntp->Column("cand",         (Int_t) j);
                ntp->Column("ncand",    (Int_t) D0.GetLength());
                ntp->Column("mode",     (Int_t) fMode);
                ntp->Column("mmiss",    (Float_t) mmiss);
                ntp->Column("msum",         (Float_t) msum);

                // store beam info
                qa.qaP4("beam", fIni, ntp);

                // store information about composite candidate tree recursively (see analysis/AnalysisTools/PndRhoTupleQA)
                qa.qaComp("x", D0[j], ntp);

                // store the 4-vector of the truth matched candidate (or a dummy, if not matched to keep ntuple consistent)
                RhoCandidate *truth = D0[j]->GetMcTruth();
                TLorentzVector lv(0,0,0,0);
                if (truth) lv = truth->P4();
                qa.qaP4("trx", lv, ntp);

                // do a vertex fit
                RhoKinVtxFitter vtxfitter(D0[j]);          // *** instantiate the vertex fitter; input is the object to be fitted
                vtxfitter.Fit();                           // *** perform fit

                RhoCandidate *cfit = D0[j]->GetFit();      // *** get the fitted candidate
                qa.qaVtx("fvxx",cfit,ntp);
                qa.qaP4("fvxx", cfit->P4(), ntp);
                double chi2_vtx = vtxfitter.GetChi2();     // *** and the chi^2 of the fit
                ntp->Column("chi2vx", (Float_t) chi2_vtx);

                // really write information to tree
                ntp->DumpData();
        }
}

void PndProdAnaTask::Finish ( )

virtual

Definition at line 186 of file PndProdAnaTask.cxx.

References RhoTuple::GetInternalTree(), and ntp.

{
        ntp->GetInternalTree()->Write();
}

InitStatus PndProdAnaTask::Init ( )

virtual

Definition at line 89 of file PndProdAnaTask.cxx.

References fAnalysis, fEvtCount, fIni, RhoTuple::GetInternalTree(), and ntp.

{
        // initialize analysis object
        fAnalysis = new PndAnalysis();

        // reset the event counter
        fEvtCount = 0;

        // *** Save current gDirectory
        TDirectory *dir = gDirectory;
        FairRootManager::Instance()->GetOutFile()->cd();

        ntp  = new RhoTuple("ntp",  "Prod Example Analysis: D0->K pi");
        ntp->GetInternalTree()->SetDirectory(gDirectory);

        // *** restore original gDirectory
        dir->cd();

        // *** the lorentz vector of the initial psi(2S)
        fIni.SetXYZT(0, 0, 6.231552, 7.240065);

        return kSUCCESS;
}

int PndProdAnaTask::SelectTruePid ( PndAnalysis *  ana, RhoCandList &  l  )

private

Definition at line 70 of file PndProdAnaTask.cxx.

References RhoCandList::GetLength(), PndAnalysis::McTruthMatch(), and RhoCandList::Remove().

{
        int removed = 0;

        for (int ii=l.GetLength()-1;ii>=0;--ii)
        {
                if ( !(ana->McTruthMatch(l[ii])) )
                {
                        l.Remove(l[ii]);
                        removed++;
                }
        }

        return removed;
}

Member Data Documentation

PndAnalysis* PndProdAnaTask::fAnalysis

private

Definition at line 64 of file PndProdAnaTask.h.

Referenced by Exec(), and Init().

int PndProdAnaTask::fEvtCount

private

Definition at line 46 of file PndProdAnaTask.h.

Referenced by Exec(), and Init().

TLorentzVector PndProdAnaTask::fIni

private

Definition at line 61 of file PndProdAnaTask.h.

Referenced by Exec(), and Init().

int PndProdAnaTask::fMode

private

Definition at line 49 of file PndProdAnaTask.h.

Referenced by Exec(), and PndProdAnaTask().

TString PndProdAnaTask::fPidAlgo

private

Definition at line 52 of file PndProdAnaTask.h.

Referenced by Exec(), and PndProdAnaTask().

RhoTuple* PndProdAnaTask::ntp

private

Definition at line 55 of file PndProdAnaTask.h.

Referenced by Exec(), Finish(), and Init().

---

The documentation for this class was generated from the following files:

• PndProdAnaTask.h
• PndProdAnaTask.cxx