PndAnaWithTrigger Class Reference

#include <PndAnaWithTrigger.h>

Inheritance diagram for PndAnaWithTrigger:

Public Member Functions
	PndAnaWithTrigger (double pbarmom, TString outname)
	~PndAnaWithTrigger ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()

Private Member Functions
	ClassDef (PndAnaWithTrigger, 1)

Private Attributes
int	fEvtCount
TLorentzVector	fIni
PndAnalysis *	fAnalysis
TDatabasePDG *	fPdg
TFile *	fFile
TString	fOutName
RhoTuple *	ntp1
RhoTuple *	ntp2
RhoTuple *	nmc
RhoMassParticleSelector *	jpsiMassSel
RhoMassParticleSelector *	jpsiPreMassSel
TClonesArray *	fOnlineFilterInfo

Detailed Description

Definition at line 24 of file PndAnaWithTrigger.h.

Constructor & Destructor Documentation

PndAnaWithTrigger::PndAnaWithTrigger	(	double	pbarmom,
		TString	outname
	)

Definition at line 63 of file PndAnaWithTrigger.cxx.

References fIni, fOutName, mp, and sqrt().

                                                                    :
  FairTask("Panda Analysis Task using Trigger info") 
{ 
        double mp=0.938272;
        fIni.SetXYZT(0,0,pbarmom, sqrt(pbarmom*pbarmom+mp*mp)+mp);

        fOutName = outname;
}

PndAnaWithTrigger::~PndAnaWithTrigger

(

)

Definition at line 75 of file PndAnaWithTrigger.cxx.

References fFile.

{ 
        delete fFile;
}

Member Function Documentation

PndAnaWithTrigger::ClassDef ( PndAnaWithTrigger  , 1    )

private

void PndAnaWithTrigger::Exec ( Option_t *  opt )

virtual

Definition at line 132 of file PndAnaWithTrigger.cxx.

References RhoKinFitter::Add4MomConstraint(), all, RhoTuple::Column(), RhoTuple::DumpData(), fAnalysis, fEvtCount, PndAnalysis::FillList(), fIni, RhoKinFitter::Fit(), fOnlineFilterInfo, RhoFitterBase::GetChi2(), PndAnalysis::GetEventInTask(), RhoCandidate::GetFit(), RhoCandidate::GetMcTruth(), PndOnlineFilterInfo::GetNTag(), PndOnlineFilterInfo::GetNTagTotal(), i, jpsiMassSel, jpsiPreMassSel, mclist, PndAnalysis::McTruthMatch(), muminus, muplus, nmc, ntp1, ntp2, RhoCandidate::P4(), piminus, piplus, PndOnlineFilterInfo::Tagged(), and TString.

{
        // *** some variables
        int i=0,j=0;
        
        // *** necessary to read the next event
        fAnalysis->GetEventInTask();
        
        // *** print event counter
        if (!(++fEvtCount%100)) cout << "evt "<<fEvtCount<<endl;
        
        // *******
        // ******* PUT ANALYSIS CODE HERE                       
        // *******
        
        // *** fetch trigger info
        PndOnlineFilterInfo *stInfo = 0;
        if (fOnlineFilterInfo) stInfo = ( PndOnlineFilterInfo* ) fOnlineFilterInfo->At ( 0 );   
        
        TString pidalg = "PidChargedProbability";
        
        // *** 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 muplus, muminus, piplus, piminus, jpsi, psi2s, all, mclist;
        
        
        // *** Select with PID info pidalg and ('All'); type and mass are set           
        fAnalysis->FillList(muplus,  "MuonAllPlus",  pidalg);
        fAnalysis->FillList(muminus, "MuonAllMinus", pidalg);
        fAnalysis->FillList(piplus,  "PionAllPlus",  pidalg);
        fAnalysis->FillList(piminus, "PionAllMinus", pidalg);

        
        // *** Setup event shape object
        fAnalysis->FillList(all,   "All", pidalg);
        PndEventShape evsh(all, fIni, 0.05, 0.1);       
        
        
        // *** store MC info in ntuple
        fAnalysis->FillList(mclist, "McTruth");
        
        nmc->Column("ev", (Int_t) fEvtCount);
        qa.qaMcList("",   mclist, nmc);
        nmc->DumpData();
        
        
        // *** combinatorics for J/psi -> mu+ mu-
        jpsi.Combine(muplus, muminus);
        jpsi.Select(jpsiPreMassSel);
        jpsi.SetType(443);
        int njmct = fAnalysis->McTruthMatch(jpsi); // match the whole list to count #matches (should be only 1)
                        
        // *** write ntuple for the jpsi reconstruction
        for (j=0;j<jpsi.GetLength();++j) 
        {
                // some general info about event (actually written for each candidate!)
                ntp1->Column("ev",              (Float_t) i);
                ntp1->Column("cand",    (Float_t) j);
                ntp1->Column("ncand",   (Float_t) jpsi.GetLength());
                ntp1->Column("nmct",    (Float_t) njmct);
                
                // *** store info from trigger
                if (stInfo)
                {
                        ntp1->Column("sttrig",   (Int_t) stInfo->Tagged());       // event triggered
                        ntp1->Column("stntot",   (Int_t) stInfo->GetNTagTotal()); // total number of triggered candidates from all active lines
                        ntp1->Column("stn200",   (Int_t) stInfo->GetNTag(200));   // number of triggered candidates from J/psi->e+ e- line
                        ntp1->Column("stn201",   (Int_t) stInfo->GetNTag(201));   // number of triggered candidates from J/psi->mu+ mu- line
                }
        
                // store info about initial 4-vector
                qa.qaP4("beam", fIni, ntp1);
                        
                // dump information about composite candidate tree recursively (see analysis/AnalysisTools/PndRhoTupleQA)
                qa.qaComp("j", jpsi[j], ntp1);
                
                // dump info about event shapes
                qa.qaEventShapeShort("es",&evsh, ntp1);
                
                // *** store the 4-vector of the truth matched candidate (or a dummy, if not matched to keep ntuple consistent)
                RhoCandidate *truth = jpsi[j]->GetMcTruth();
                TLorentzVector lv;
                if (truth) lv = truth->P4();
                qa.qaP4("trj", lv, ntp1);
                
                ntp1->DumpData();
        }
        
        
        // *** some rough mass selection on J/psi
        jpsi.Select(jpsiMassSel);
        
        // *** combinatorics for psi(2S) -> J/psi pi+ pi-
        psi2s.Combine(jpsi, piplus, piminus);
        psi2s.SetType(88880);
        int npsimct = fAnalysis->McTruthMatch(psi2s); // match the whole list to count #matches (should be only 1)

        // *** write ntuple for the psi(2S) reconstruction
        for (j=0;j<psi2s.GetLength();++j) 
        {
                // some general info about event (actually written for each candidate!)
                ntp2->Column("ev",              (Float_t) i);
                ntp2->Column("cand",    (Float_t) j);
                ntp2->Column("ncand",   (Float_t) psi2s.GetLength());
                ntp2->Column("nmct",    (Float_t) npsimct);

                // *** store info from trigger
                if (stInfo)
                {
                        ntp2->Column("sttrig",   (Int_t) stInfo->Tagged());       // event triggered
                        ntp2->Column("stntot",   (Int_t) stInfo->GetNTagTotal()); // total number of triggered candidates from all active lines
                        ntp2->Column("stn200",   (Int_t) stInfo->GetNTag(200));   // number of triggered candidates from J/psi->e+ e- line
                        ntp2->Column("stn201",   (Int_t) stInfo->GetNTag(201));   // number of triggered candidates from J/psi->mu+ mu- line
                }
                
                RhoKinFitter kinfit(psi2s[j]);
                kinfit.Add4MomConstraint(fIni);
                kinfit.Fit();
                
                RhoCandidate *psifit=psi2s[j]->GetFit();
                // store info about initial 4-vector
                qa.qaP4("beam", fIni, ntp2);
                
                // dump information about composite candidate tree recursively (see analysis/AnalysisTools/PndRhoTupleQA)
                qa.qaComp("psi", psi2s[j], ntp2);
                qa.qaComp("fpsi",psifit, ntp2);
                ntp2->Column("fchi2", (Float_t) kinfit.GetChi2());

                // dump info about event shapes
                qa.qaEventShapeShort("es",&evsh, ntp2);
                
                // *** store the 4-vector of the truth matched candidate (or a dummy, if not matched to keep ntuple consistent)
                RhoCandidate *truth = psi2s[j]->GetMcTruth();
                TLorentzVector lv;
                if (truth) lv = truth->P4();
                qa.qaP4("trpsi", lv, ntp2);
                
                ntp2->DumpData();
        }                       
        
}

void PndAnaWithTrigger::Finish ( )

virtual

Definition at line 278 of file PndAnaWithTrigger.cxx.

References fFile.

{
        
        // *******
        // ******* STORE YOUR HISTOS AND TUPLES
        // *******
        
        fFile->Write();
        fFile->Close(); 

}

InitStatus PndAnaWithTrigger::Init ( )

virtual

Definition at line 83 of file PndAnaWithTrigger.cxx.

References fAnalysis, fEvtCount, fFile, fOnlineFilterInfo, fOutName, fPdg, RhoTuple::GetInternalTree(), jpsiMassSel, jpsiPreMassSel, PndAnalysis::McMatchAllowPhotos(), nmc, ntp1, and ntp2.

{               
        // *** initialize PndAnalysis object
        fAnalysis = new PndAnalysis();
        fAnalysis->McMatchAllowPhotos();
        
        // *** reset the event counter
        fEvtCount = 0;

        // *******
        // ******* PREPARE/CREATE THE STUFF YOU NEED
        // *******
        
        fPdg = TDatabasePDG::Instance();

        // *** Mass selector for the jpsi cands
        double m0_jpsi = fPdg->GetParticle("J/psi")->Mass();   // Get nominal PDG mass of the J/psi
        jpsiMassSel    = new RhoMassParticleSelector("jpsi"   ,m0_jpsi,0.5);
        jpsiPreMassSel = new RhoMassParticleSelector("jpsipre",m0_jpsi,2.0);
        
        // ***
        // *** Prepare RhoTuple output  
        // ***
        TDirectory *dir = gDirectory;   
    fFile=new TFile(fOutName,"RECREATE");
        fFile->cd();
        
        // *** create some ntuples
        ntp1 = new RhoTuple("ntp1", "jpsi analysis");
        ntp2 = new RhoTuple("ntp2", "psi(2S) analysis");
        nmc  = new RhoTuple("nmc",  "mctruth info");

        // assign RhoTuples to outfile
        if (ntp1) ntp1->GetInternalTree()->SetDirectory(gDirectory);
        if (ntp2) ntp2->GetInternalTree()->SetDirectory(gDirectory);
        if (nmc)  nmc->GetInternalTree()->SetDirectory(gDirectory);

        // *** restore original gDirectory
        dir->cd();
        
        // *** Connect to the Online Filter Info
        fOnlineFilterInfo = ( TClonesArray* ) FairRootManager::Instance()->GetObject ( "OnlineFilterInfo" );    
        
        return kSUCCESS;
}

Member Data Documentation

PndAnalysis* PndAnaWithTrigger::fAnalysis

private

Definition at line 54 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

int PndAnaWithTrigger::fEvtCount

private

Definition at line 48 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

TFile* PndAnaWithTrigger::fFile

private

Definition at line 56 of file PndAnaWithTrigger.h.

Referenced by Finish(), Init(), and ~PndAnaWithTrigger().

TLorentzVector PndAnaWithTrigger::fIni

private

Definition at line 51 of file PndAnaWithTrigger.h.

Referenced by Exec(), and PndAnaWithTrigger().

TClonesArray* PndAnaWithTrigger::fOnlineFilterInfo

private

Definition at line 72 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

TString PndAnaWithTrigger::fOutName

private

Definition at line 57 of file PndAnaWithTrigger.h.

Referenced by Init(), and PndAnaWithTrigger().

TDatabasePDG* PndAnaWithTrigger::fPdg

private

Definition at line 55 of file PndAnaWithTrigger.h.

Referenced by Init().

RhoMassParticleSelector* PndAnaWithTrigger::jpsiMassSel

private

Definition at line 68 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

RhoMassParticleSelector* PndAnaWithTrigger::jpsiPreMassSel

private

Definition at line 69 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

RhoTuple* PndAnaWithTrigger::nmc

private

Definition at line 66 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

RhoTuple* PndAnaWithTrigger::ntp1

private

Definition at line 64 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

RhoTuple* PndAnaWithTrigger::ntp2

private

Definition at line 65 of file PndAnaWithTrigger.h.

Referenced by Exec(), and Init().

---

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

• PndAnaWithTrigger.h
• PndAnaWithTrigger.cxx