PndSoftTriggerTask Class Reference

#include <PndSoftTriggerTask.h>

Inheritance diagram for PndSoftTriggerTask:

Public Types
typedef std::map< Int_t, Float_t >	mapper

Public Member Functions
	PndSoftTriggerTask (double pmom, int mode=0, int runnum=0, TString trigfilename="")
	~PndSoftTriggerTask ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()
void	SetFastSimDefaults ()
void	SetFullSimDefaults ()
void	SetConfigurationFile (TString fname)
void	ApplyFullSelection (int selmode=1)
void	SetTagNSigMode (int mode, double nsig)
void	SetTagNSigAll (double nsig)
void	SetAuxNSig (double nsig)
void	SetDstMDiffCut (double cut)
void	SetPi0QASelection (double min, double max)
void	SetKs0QASelection (double min, double max)
void	SetEtaQASelection (double min, double max)
void	SetQASelectionDefaults ()
void	SetPi0SignalParams (double mean, double sigma)
void	SetKs0SignalParams (double mean, double sigma)
void	SetEtaSignalParams (double mean, double sigma)
void	SetSignalParamsDefaults ()
void	SetGammaMinE (double min)
void	SetTrackMinP (double min)
void	SetInitialPidCut (double cut)
void	SetPidAlgoAll (TString algo)
void	SetPidAlgoElectron (TString algo)
void	SetPidAlgoMuon (TString algo)
void	SetPidAlgoPion (TString algo)
void	SetPidAlgoKaon (TString algo)
void	SetPidAlgoProton (TString algo)
void	SetTagAll (bool tag=true)
void	SetTagMode (int mode, bool tag=true)
void	SetQAPi0 (bool qa=true)
void	SetQAEta (bool qa=true)
void	SetQAKs0 (bool qa=true)
void	SetQAMc (bool qa=true)
void	SetQAMctOnly (bool qa=true)
void	SetQAEvent (bool qa=true)
void	SetQAAll (bool qa=true)
void	SetQAMode (int mode, bool qa=true)
void	McMatchAllowPhotos (int maxn=1, double thresh=0.05)
void	SetVerbose (int verb=10)
	PndSoftTriggerTask ()
	~PndSoftTriggerTask ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()

Protected Attributes
int	evcount
int	chmax
int	neutmax
int	mcmax
int	allmax
int	cntsel
int	cntjpsi
int	cntd0
int	cntdpm
int	cntds
int	cntphi
int	cntlamc

Private Member Functions
void	FillGlobalLists ()
bool	ReadConfiguration ()
bool	ReadTriggerLines ()
int	DoCombinatorics (RhoCandList &l, PndSoftTriggerLine *tl)
void	CombineList (RhoCandList &l, int mothpdg, int amothpdg, std::vector< int > &idx, std::vector< int > &aidx, bool cc=false)
int	AntiPdg (int pdg)
double	GetVarValue (RhoCandidate *c, int id)
void	FillVarArray (RhoCandidate *c, int id, Bool_t tmva=false)
void	FillEventShapeVarArray ()
bool	AcceptCandidate (int mode, RhoCandidate *c, RhoParticleSelectorBase *sel=0)
bool	AcceptDstCut (RhoCandidate *c)
double	GetPocaVtx (RhoCandidate *c, double &dist, double &ctau)
void	GetAngles (RhoCandidate *c, double &oang, double &decang)
int	SelectTruePid (RhoCandList &l)
int	SelectPidProb (RhoCandList &l, int pididx, double cut)
int	MultPidProb (RhoCandList &l, int pididx, double prob)
double	DbMass (TString name)
int	SplitString (TString s, TString delim, TString *toks, int maxtoks)
int	DetermineRecoilMode (int &mode)
TLorentzVector	BoostCms (TLorentzVector in)
int	CodeVariable (TString v)
int	CreateKs0Cands (RhoTuple *n)
int	TagMode (PndSoftTriggerLine *tl, int &npre)
virtual void	SetParContainers ()
	ClassDef (PndSoftTriggerTask, 1)
void	FillMassHisto (TH1F *h, RhoCandList &l)
int	SelectPdgCode (RhoCandList &mct, RhoCandList &l)
void	PrintList (RhoCandList &l, int max=10)
int	RemoveDoubles (RhoCandList &l, double limit=0.0001)
void	SelectPid (int type, int pdg, int chrg, RhoCandList &l, RhoCandList &lpid, double cut=0.2)
void	ConfigureHistos (TH1F *hall, TH1F *htrue, TH1F *hsel)
virtual void	SetParContainers ()
	ClassDef (PndSoftTriggerTask, 1)

Private Attributes
int	fVerbose
int	fMode
int	fRecoilMode
int	fRecoilCnt
int	fEvtCount
int	fRunNum
int	fSigCount
double	fNsigTag
double	fNsigAux
double	fDstMDiffCut
TString	fCfgFileName
TString	fTriggerFileName
Int_t	fPhotosMax
Double_t	fPhotosThresh
int	fApplyFullSelection
TLorentzVector	fIniP4
double	fEcm
double	fPbarMom
TVector3	fPrimVtx
double	fPrimVtxQa
bool	fQAPi0
bool	fQAEta
bool	fQAKs0
bool	fQAEvent
bool	fQAMc
bool	fQAMctOnly
double	fPi0QaMin
double	fPi0QaMax
double	fEtaQaMin
double	fEtaQaMax
double	fKs0QaMin
double	fKs0QaMax
double	fPi0Mean
double	fPi0Sigma
double	fEtaMean
double	fEtaSigma
double	fKs0Mean
double	fKs0Sigma
double	fGammaMinE
double	fPi0MinE
double	fEtaMinE
double	fTrackMinP
double	fIniPidCut
RhoMassParticleSelector *	fPi0Sel
RhoMassParticleSelector *	fPi0PreSel
RhoMassParticleSelector *	fEtaSel
RhoMassParticleSelector *	fEtaPreSel
RhoMassParticleSelector *	fKs0Sel
RhoMassParticleSelector *	fKs0PreSel
RhoMomentumParticleSelector *	fMomentumSel
RhoEnergyParticleSelector *	fEnergySel
TString	fAlgoElectron
TString	fAlgoMuon
TString	fAlgoPion
TString	fAlgoKaon
TString	fAlgoProton
Int_t	fPidMult_025 [5]
RhoTuple *	ntp
RhoTuple *	nks0
RhoTuple *	npi0
RhoTuple *	neta
RhoTuple *	nmc
PndAnalysis *	fAnalysis
PndEventShape *	fEventShape
PndRhoTupleQA *	fQA
RhoVtxPoca *	fPocaVertexer
TDatabasePDG *	fPdg
RhoCandList	fMcTruth
RhoCandList	fAllCands
RhoCandList	fChargedCands
RhoCandList	fNeutralCands
RhoCandList	fGammaCands
RhoCandList	fKs0Cands
RhoCandList	fPi0Cands
RhoCandList	fEtaCands
RhoCandList	fPidList [16]
FairRootManager *	fRootManager
TClonesArray *	fTcaOnlineFilterInfo
TH1F *	h_mom
TH1F *	h_momc
TH1F *	h_momn
TH1F *	h_mult
TH1F *	h_multc
TH1F *	h_multn
TH1F *	h_multrem
TH1F *	h_jpsi
TH1F *	h_jpsit
TH1F *	h_jpsisel
TH1F *	h_d0
TH1F *	h_d0t
TH1F *	h_d0sel
TH1F *	h_d0_pocx
TH1F *	h_d0_pocy
TH1F *	h_d0_pocz
TH1F *	h_d0_pocxm
TH1F *	h_d0_pocym
TH1F *	h_d0_poczm
TH1F *	h_d0_l
TH1F *	h_d0_lm
TH1F *	h_dpm
TH1F *	h_dpmt
TH1F *	h_dpmsel
TH1F *	h_ds
TH1F *	h_dst
TH1F *	h_dssel
TH1F *	h_phi
TH1F *	h_phit
TH1F *	h_phisel
TH1F *	h_lamc
TH1F *	h_lamct
TH1F *	h_lamcsel
RhoTuple *	ntp2
PndAnalysis *	theAnalysis
TStopwatch *	timer

Detailed Description

Definition at line 31 of file softrig/PndSoftTriggerTask.h.

Member Typedef Documentation

typedef std::map<Int_t, Float_t> PndSoftTriggerTask::mapper

Definition at line 25 of file tutorials/anatask/PndSoftTriggerTask.h.

Constructor & Destructor Documentation

PndSoftTriggerTask::PndSoftTriggerTask	(	double	pmom,
		int	mode = 0,
		int	runnum = 0,
		TString	trigfilename = ""
	)

Definition at line 134 of file softrig/PndSoftTriggerTask.cxx.

References fAlgoElectron, fAlgoKaon, fAlgoMuon, fAlgoPion, fAlgoProton, fApplyFullSelection, fCfgFileName, fEcm, fIniP4, fIniPidCut, fPdg, fSTMapListIndex, fSTPidIndex, fTriggerFileName, i, mp, ReadTriggerLines(), SetQASelectionDefaults(), SetSignalParamsDefaults(), sqrt(), and TString.

                                                                                              :
        FairTask("Panda Softtrigger Task"),
        fVerbose(0), fMode(mode), fEvtCount(0), fRunNum(runnum), fSigCount(0), fNsigTag(8.0), fNsigAux(3.0), fDstMDiffCut(10.),
        fTriggerFileName(trigfilename), fPhotosMax(0), fPhotosThresh(0.05), 
        fIniP4(0,0,0,0), fEcm(0.), fPbarMom(pmom),
        fQAPi0(false),fQAEta(false),fQAKs0(false),fQAEvent(false), fQAMc(false), fQAMctOnly(false),
        fGammaMinE(0.03), fPi0MinE(0.0), fEtaMinE(0.0), fTrackMinP(0.15), fIniPidCut(0.0),
        fPi0Sel(NULL), fEtaSel(NULL), fKs0Sel(NULL),
        fMomentumSel(NULL), fEnergySel(NULL),
        ntp(0), nks0(0), npi0(0), neta(0), nmc(0), 
        fEventShape(NULL),
        fQA(NULL)       
{
        fPdg = TDatabasePDG::Instance();
        
        // *** add several pbar p/n/dd Systems for MC truth match
        double ppwidth = 0.01;
        fPdg->AddParticle("pbarpSystem", "pbar p", fEcm, false, ppwidth,0,"",88888);
        fPdg->AddParticle("pbarpSystem0","pbar p", fEcm, false, ppwidth,0,"",88880);
        fPdg->AddParticle("pbarpSystem1","pbar p", fEcm, false, ppwidth,0,"",88881);
        fPdg->AddParticle("pbarpSystem2","pbar p", fEcm, false, ppwidth,0,"",88882);
        
        fPdg->AddParticle("pbarnSystem", "pbar n", fEcm, false, ppwidth,0,"",88887);
        fPdg->AddParticle("pbardSystem", "pbar d", fEcm, false, ppwidth,0,"",88889);
        
        double mp = fPdg->GetParticle("proton")->Mass(); //Proton mass for computation of p4_ini
        
        // set 4-vector of pbar-p-system
        fIniP4.SetPz(pmom);
        fIniP4.SetE(sqrt(pmom*pmom+mp*mp)+mp);
        fEcm = fIniP4.M();
        
        // set default algorithms for pid
        fAlgoElectron   = "PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc;PidAlgoStt;PidAlgoMdtHardCuts";
        fAlgoMuon               = fAlgoElectron;
        fAlgoPion               = fAlgoElectron; 
        fAlgoKaon               = fAlgoElectron;
        fAlgoProton             = fAlgoElectron; 
        
        fIniPidCut              = 0.0; 
        
        fCfgFileName    = TString(gSystem->Getenv("VMCWORKDIR"))+"/softrig/selection_10ch_tight.cfg";
        fApplyFullSelection = false;
        
        // *** set default signal parameters (mean, sigma)
        SetSignalParamsDefaults();
        SetQASelectionDefaults();       
        
        // if no trigger definition file is given, use the default one
        if (fTriggerFileName.Length()==0)
        {
                fTriggerFileName = TString(gSystem->Getenv("VMCWORKDIR"))+"/softrig/triggerlines.cfg";
                cout <<"Reading default trigger lines file "<<fTriggerFileName.Data()<<endl;
        }       
        // map pdg codes to RhoCandList index
        for (int i=0;i<14;++i) fSTMapListIndex[fSTPidIndex[i]] = i;
        
        ReadTriggerLines();
}

PndSoftTriggerTask::~PndSoftTriggerTask

(

)

Definition at line 197 of file softrig/PndSoftTriggerTask.cxx.

References fPocaVertexer.

{
        delete fPocaVertexer;
}

PndSoftTriggerTask::PndSoftTriggerTask

(

)

Default constructor

Definition at line 66 of file tutorials/anatask/PndSoftTriggerTask.cxx.

                                       :
  FairTask("Panda Analysis Task") { 
}

PndSoftTriggerTask::~PndSoftTriggerTask

(

)

Destructor

Member Function Documentation

bool PndSoftTriggerTask::AcceptCandidate ( int  mode, RhoCandidate *  c, RhoParticleSelectorBase *  sel = 0  )

private

Definition at line 1421 of file softrig/PndSoftTriggerTask.cxx.

References RhoParticleSelectorBase::Accept(), AcceptDstCut(), STCutSet::cutval, FillVarArray(), fSTencode, fSTInputValues, fSTModeIndex, fSTOps, fSTSelmap, fVerbose, i, mode, STCutSet::ncut, STCutSet::op, STCutSet::tmvacut, STCutSet::tmvameth, STCutSet::tmvaread, and STCutSet::varid.

Referenced by TagMode().

{
        int mcode = fSTencode[fSTModeIndex]*1000+mode;  
        
        if (fVerbose>1) cout <<"AcceptCandidate : mode="<<mode<<" mcode="<<mcode;

        // no selection defined for this mode
        if ( fSTSelmap.find(mcode) == fSTSelmap.end() ) {if (fVerbose>1) cout <<endl;return false;} 
        if (fVerbose>1) cout <<" found cut set";
        
        // not accepted by final (mass) selector
        if ( sel && !(sel->Accept(c)) ) {if (fVerbose>1) cout <<endl;return false;}
        // not accepted by D* mass diff cut
        if (!AcceptDstCut(c)) return false;
                
        if (fVerbose>1) cout <<" accepted by precuts"<<endl;
        
        // *** 
        // *** check for cut selection
        // *** 
        STCutSet cs = fSTSelmap[mcode];
        FillVarArray(c, mcode);
        
        bool acc = true;
        
        for (int i=0;i<cs.ncut;++i)
        {
                if (fVerbose>1) cout <<"  -> checking var["<<cs.varid[i]<<"] ("<<fSTInputValues[i]<<") "<<fSTOps[cs.op[i]]<<" "<<cs.cutval[i]<<endl;
                
                switch (cs.op[i]) // which operator is used for this cut?
                {
                case 0: // check var >= value
                        if ( !(fSTInputValues[i]>=cs.cutval[i]) ) acc=false; 
                        break;
                case 1: // check var > value
                        if ( !(fSTInputValues[i]>cs.cutval[i]) ) acc=false; 
                        break;
                case 2: // check var == value
                        if ( !(fSTInputValues[i]==cs.cutval[i]) ) acc=false; 
                        break;
                case 3: // check var <= value
                        if ( !(fSTInputValues[i]<=cs.cutval[i]) ) acc=false; 
                        break;
                case 4: // check var < value
                        if ( !(fSTInputValues[i]<cs.cutval[i]) ) acc=false; 
                        break;
                default : 
                        acc=false; 
                        break;
                }
        }
        
        // *** 
        // *** check for TMVA selector in addition
        // *** 
        if (cs.tmvaread)
        {
                FillVarArray(c, mcode, true);
                if (fVerbose>1) cout <<"  -> checking TMVA "<<cs.tmvameth<<endl;
                float mvaout = cs.tmvaread->EvaluateMVA(cs.tmvameth);
                if (mvaout<cs.tmvacut) acc=false;
        }
        
        
        if (fVerbose>1)
        {
                if (acc) cout <<" final accept"<<endl;
                else cout <<endl;
        }
        
        return acc;
}

bool PndSoftTriggerTask::AcceptDstCut ( RhoCandidate *  c )

private

Definition at line 1410 of file softrig/PndSoftTriggerTask.cxx.

References RhoCandidate::Daughter(), fabs(), fDstMDiffCut, RhoCandidate::M(), and RhoCandidate::PdgCode().

Referenced by AcceptCandidate(), and TagMode().

{
        int pdg = abs(c->PdgCode());
        
        if ( (pdg==413 || pdg==423 || pdg==433) && fabs(c->M() - c->Daughter(0)->M() - 0.143) > fDstMDiffCut ) return false;
        
        return true;
}

int PndSoftTriggerTask::AntiPdg ( int  pdg )

private

Definition at line 1497 of file softrig/PndSoftTriggerTask.cxx.

References fPdg.

Referenced by DoCombinatorics().

{
        int antipdg = pdg;
        TParticlePDG *part = fPdg->GetParticle(pdg)->AntiParticle();
        if (part) antipdg = part->PdgCode();
        
        return antipdg;
}

void PndSoftTriggerTask::ApplyFullSelection ( int  selmode = 1 )

inline

Definition at line 61 of file softrig/PndSoftTriggerTask.h.

References fApplyFullSelection.

Referenced by ana_jpsi_task_ST(), prod_fsim(), SetFastSimDefaults(), and SetFullSimDefaults().

{fApplyFullSelection = selmode;}  

TLorentzVector PndSoftTriggerTask::BoostCms ( TLorentzVector  in )

private

Definition at line 1183 of file softrig/PndSoftTriggerTask.cxx.

References fIniP4.

Referenced by GetVarValue().

{
        TLorentzVector lcm = l;
        lcm.Boost(-fIniP4.BoostVector());
        return lcm;
}

PndSoftTriggerTask::ClassDef ( PndSoftTriggerTask  , 1    )

private

PndSoftTriggerTask::ClassDef ( PndSoftTriggerTask  , 1    )

private

int PndSoftTriggerTask::CodeVariable ( TString  v )

private

Definition at line 1255 of file softrig/PndSoftTriggerTask.cxx.

References C(), fSTNQuant, fSTQuantCode, fSTVarSuff, i, r1, r2, TString, and v.

Referenced by ReadConfiguration().

{
        int i=0, code=0, A=9, B=9, C=9, D=9, EF=0;
        
        if (v=="mmiss") return 80;
        if (v=="xmdif") return 82;
        if (v=="tag")   return 99;
        
        for (i=0; i<fSTNQuant; ++i) 
                if (v.EndsWith(fSTVarSuff[i])) 
                {
                        EF=fSTQuantCode[i]; 
                        v=v(0,v.Length()-fSTVarSuff[i].Length());
                        i=1000;
                } 
        
        if (i<1000) {cout <<"[PndSoftTriggerTask] ***** Unknown variable suffix in variable '"<<v<<"'"<<endl; return -1;}
        
        // get rid of special daughter namings for pi0 gammas and KS pi+-
        v.ReplaceAll("pi1","d0");
        v.ReplaceAll("pi2","d1");
        v.ReplaceAll("g1","d0");
        v.ReplaceAll("g2","d1");
        
        TRegexp r1("xd[0-4]$");
        TRegexp r2("xd[0-4]d[0-4]$");
        TRegexp r3("xd[0-4]d[0-4]d[0-4]$");
        
        if (v!="x" && v(r1)=="" && v(r2)=="" && v(r3)=="") {cout <<"[PndSoftTriggerTask] ***** Unknown candidate prefix in variable '"<<v<<"'"<<endl; return -1;}
        
        int l = v.Length();
        
        if (l>1) B = TString(v(2,1)).Atoi();
        if (l>3) C = TString(v(4,1)).Atoi();
        if (l>5) D = TString(v(6,1)).Atoi();
        
        code=EF+100*D+1000*C+10000*B+100000*A;
        
        return code;
}

void PndSoftTriggerTask::CombineList ( RhoCandList &  l, int  mothpdg, int  amothpdg, std::vector< int > &  idx, std::vector< int > &  aidx, bool  cc = false  )

private

Definition at line 1508 of file softrig/PndSoftTriggerTask.cxx.

References RhoCandList::Cleanup(), RhoCandList::Combine(), RhoCandList::CombineAndAppend(), and fPidList.

Referenced by DoCombinatorics().

{
        l.Cleanup();
        
        int nd = idx.size();
        
        switch (nd) 
        {
        case 2: 
                l.Combine(fPidList[idx[0]], fPidList[idx[1]], mothpdg); 
                if (cc) l.CombineAndAppend(fPidList[aidx[0]], fPidList[aidx[1]], amothpdg);
                break;
                
        case 3: 
                l.Combine(fPidList[idx[0]], fPidList[idx[1]], fPidList[idx[2]], mothpdg);
                if (cc) l.CombineAndAppend(fPidList[aidx[0]], fPidList[aidx[1]], fPidList[aidx[2]], amothpdg);
                break;
                
        case 4: 
                l.Combine(fPidList[idx[0]], fPidList[idx[1]], fPidList[idx[2]], fPidList[idx[3]], mothpdg); 
                if (cc) l.CombineAndAppend(fPidList[aidx[0]], fPidList[aidx[1]], fPidList[aidx[2]], fPidList[aidx[3]], amothpdg); 
                break;

        case 5: 
                l.Combine(fPidList[idx[0]], fPidList[idx[1]], fPidList[idx[2]], fPidList[idx[3]], fPidList[idx[4]], mothpdg); 
                if (cc) l.CombineAndAppend(fPidList[aidx[0]], fPidList[aidx[1]], fPidList[aidx[2]], fPidList[aidx[3]], fPidList[aidx[4]], amothpdg); 
                break;
        default: return;
        }
        
}

void PndSoftTriggerTask::ConfigureHistos ( TH1F *  hall, TH1F *  htrue, TH1F *  hsel  )

private

Definition at line 457 of file tutorials/anatask/PndSoftTriggerTask.cxx.

{
        htrue->SetLineColor(2);
        hsel->SetFillStyle(3001);
        hsel->SetFillColor(602);
        hsel->SetLineStyle(3);
        hsel->SetLineColor(602);
        
}

int PndSoftTriggerTask::CreateKs0Cands ( RhoTuple *  n )

private

Definition at line 977 of file softrig/PndSoftTriggerTask.cxx.

References RhoTuple::Column(), RhoCandList::Combine(), RhoTuple::DumpData(), f, fabs(), fEvtCount, fIniP4, fKs0Cands, fKs0Mean, fKs0PreSel, fKs0Sel, fKs0Sigma, fMode, fNsigAux, fPidList, fQA, fRunNum, RhoCandList::GetLength(), i, Mass, nsig, PndRhoTupleQA::qaKs0(), PndRhoTupleQA::qaP4(), RhoCandList::Select(), RhoCandList::SetType(), and tag.

Referenced by FillGlobalLists().

{
        // *** basic KS0 Reco
        //fKs0Cands.Combine(fPionPlus, fPionMinus);
        fKs0Cands.Combine(fPidList[4], fPidList[5]);
        // *** pre select for qa
        fKs0Cands.Select(fKs0PreSel);
        fKs0Cands.SetType(310);
        
        if (n)
        {
                for (int i=0; i<fKs0Cands.GetLength();++i)
                {                       
                        Float_t nsig = (Float_t) fabs(fKs0Cands[i]->Mass()-fKs0Mean)/fKs0Sigma;
                        Float_t tag  = nsig<fNsigAux;  // instead of fNsigTag, fNsigAux cut is applied for KS
                        
                        n->Column("tag",        (Int_t) tag,            0);
                        n->Column("nsig",       (Float_t) nsig,         0.0f);
                        n->Column("ev",         (Int_t) fEvtCount,      0);
                        n->Column("run",        (Int_t) fRunNum,        0);
                        n->Column("num",        (Int_t) i,                      0);
                        n->Column("mode",       (Int_t) fMode,          0);
                        n->Column("xmean", (Float_t) fKs0Mean,          0.0f);
                        n->Column("xsig",       (Float_t) fKs0Sigma,    0.0f);
                        
                        fQA->qaP4("beam", fIniP4, n);
                        fQA->qaKs0("x",  fKs0Cands[i], n);
                        
                        n->DumpData();
                }
        }
        
        // *** final selection
        fKs0Cands.Select(fKs0Sel);
        
        return fKs0Cands.GetLength();
}

double PndSoftTriggerTask::DbMass ( TString  name )

inlineprivate

Definition at line 146 of file softrig/PndSoftTriggerTask.h.

References fPdg.

Referenced by SetQASelectionDefaults().

{if (fPdg->GetParticle(name)) return fPdg->GetParticle(name)->Mass(); else return 0.;}

int PndSoftTriggerTask::DetermineRecoilMode ( int &  mode )

private

Definition at line 933 of file softrig/PndSoftTriggerTask.cxx.

References fMcTruth, RhoCandList::GetLength(), and i.

Referenced by Exec().

{
        mode = 0;
        int code=-1;
        
        if ( fMcTruth.GetLength()==0 ) return -1;
        if ( fMcTruth[0]==0 ) return -1;
        
        for (int i=0; i<fMcTruth[0]->NDaughters(); ++i)
        {
                int dpdg = abs(fMcTruth[0]->Daughter(i)->PdgCode());
                
                switch (dpdg)
                {
                        case 22:  mode += 1;      break;  // gamma
                        case 111: mode += 10;     break;  // pi0 
                        case 211: mode += 100;    break;  // pi+
                        case 321: mode += 1000;   break;  // K+
                        case 130: mode += 10000;  break;  // K0
                        case 310: mode += 10000;  break;  // K0
                        case 221: mode += 100000; break;  // eta
                }
        }
        
        switch (mode) 
        {
                case       0 : code = 0;  break;  // no recoil
                case       1 : code = 1;  break;  // gamma
                case      10 : code = 2;  break;  // pi0
                case  100000 : code = 3;  break;  // eta
                case      20 : code = 4;  break;  // pi0 pi0
                case     200 : code = 5;  break;  // pi+ pi-
                case    2000 : code = 6;  break;  // K+ K-
                case   20000 : code = 7;  break;  // K0 K0b
                case  200000 : code = 8;  break;  // eta eta
                case     210 : code = 9;  break;  // pi+ pi- pi0
                default      : code =-1;  break;
        }
        
        return code;
}

int PndSoftTriggerTask::DoCombinatorics ( RhoCandList &  l, PndSoftTriggerLine *  tl  )

private

Definition at line 1542 of file softrig/PndSoftTriggerTask.cxx.

References AntiPdg(), RhoCandList::Cleanup(), CombineList(), fPidList, fSTMapListIndex, fSTPidIndex, fVerbose, PndSoftTriggerLine::GetAuxNeeded(), PndSoftTriggerLine::GetCC(), PndSoftTriggerLine::GetDaughterPdg(), RhoCandList::GetLength(), PndSoftTriggerLine::GetMotherPdg(), PndSoftTriggerLine::GetNDaughters(), i, and idx.

Referenced by TagMode().

{
        l.Cleanup();
        fPidList[14].Cleanup();
        fPidList[15].Cleanup();
        
        int mothpdg = tl->GetMotherPdg();
        if (mothpdg<0) {cout << "[PndSoftTriggerTask] **** Invalid mother for combinatorics."<<endl; return 0;}
        
        int nd = tl->GetNDaughters();    // cache for mapping: pdgcode -> fPdgList indices
        
        std::vector<int> idx, aidx, auxidx, auxaidx;
        
        
        // fetch pdg code of antiparticle for mother
        int amothpdg = AntiPdg(mothpdg);
        
        // do we need to add charged conjugate combinatorics?
        bool cc = tl->GetCC();
        
        if (fVerbose>1)
        {
                for (int i=0; i<nd; ++i) cout <<tl->GetDaughterPdg(i)<<" ";     
                if (cc) cout <<" +cc";
                cout <<endl;
        }
        
        // flag whether pdg goes to idx list or auxidx list
        bool auxmode=false; 

        // idx for the mothers of the aux list
        int auxmothpdg, auxamothpdg; 
        
        // cache the indices of the fPidList according to the pdg codes
        for (int i=0; i<nd; ++i)
        {
                int dpdg = tl->GetDaughterPdg(i);
                
                // skip the codes for '[' and ']'
                if (dpdg==-99) continue;
                if (dpdg==-98) {auxmode=false; continue;} // filling auf aux list finished

                // if no list for particle with code 'dpdg'
                if ( fSTMapListIndex.find(dpdg) == fSTMapListIndex.end() )
                {
                        // is it start of aux list definition? (e.g. D*0 -> D0 [K- pi+] pi0)
                        if (tl->GetDaughterPdg(i+1) != -99)     
                        {
                                // if not, throw error!
                                cout << "[PndSoftTriggerTask] **** Invalid daughter pdg code: "<<dpdg<<endl; 
                                return 0;
                        }
                        else // start aux mode and decide about cc of aux list
                        {
                                auxmode = true;
                                auxmothpdg  = dpdg;
                                auxamothpdg = AntiPdg(dpdg);
                                idx.push_back(14);
                                if (auxamothpdg==auxmothpdg) aidx.push_back(14);
                                else aidx.push_back(15);
                                continue;
                        }
                }
                if (auxmode)
                {
                        auxidx.push_back(fSTMapListIndex[dpdg]);
                        auxaidx.push_back(fSTMapListIndex[AntiPdg(dpdg)]); 
                }
                else
                {
                        idx.push_back(fSTMapListIndex[dpdg]);
                        aidx.push_back(fSTMapListIndex[AntiPdg(dpdg)]); 
                }
        }
        
        // does one of the lists have too many daughters?
        if (idx.size()>5)    {cout << "[PndSoftTriggerTask] **** Too many daughters ("<<idx.size()<<") for combinatorics of target resonance."<<endl;return 0;}
        if (auxidx.size()>5) {cout << "[PndSoftTriggerTask] **** Too many daughters ("<<auxidx.size()<<") for combinatorics of auxiliary resonance."<<endl;return 0;}
        
        if (fVerbose>1)
        {
                cout <<"("<<mothpdg<<"/"<<amothpdg<<") : ";
                for (size_t i=0; i<idx.size(); ++i) cout <<"("<<idx[i]<<"/"<<aidx[i]<<") "; 
                cout <<endl<<endl;
                if (auxidx.size()>0)
                {
                        cout <<"("<<auxmothpdg<<"/"<<auxamothpdg<<") : ";
                        for (size_t i=0; i<auxidx.size(); ++i) cout <<"("<<auxidx[i]<<"/"<<auxaidx[i]<<") "; cout <<endl;
                }
        }
        
        // do the combinatorics for aux list and final list
        if (tl->GetAuxNeeded())
        {
                // create different aux lists is pdg!=anti-pdg (e.g. D0 -> K- pi+ and anti-D0 -> K+ pi-)
                if (auxmothpdg!=auxamothpdg)
                {
                        CombineList(fPidList[14], auxmothpdg, 0,  auxidx, auxidx, false);
                        CombineList(fPidList[15], auxamothpdg, 0, auxaidx, auxaidx, false);
                }
                // create one aux list e.g. for eta_c -> ...
                else 
                {
                        // check whether final state is its anti-final state (e.g. pi+ pi- pi0 = pi- pi+ pi0)
                        // or not (e.g. eta_c -> KS K- pi+ != KS K+ pi-)
                        // this can be done by summing all pdg codes for particles with pdg!=anti-pdg
                        // if sum = 0, final state and anti-final state are the same (i.e. FS has for each particle the according anti-particle)
                        int pdgsum = 0;
                        for (size_t i=0;i<auxidx.size();++i) if (auxidx[i]<10) pdgsum+=fSTPidIndex[auxidx[i]];
                        
                        // if FS = anti-FS, just do combinatorics once (not adding the composites from anti-FS list)
                        if (pdgsum==0) 
                                CombineList(fPidList[14], auxmothpdg, 0,  auxidx, auxidx, false);
                        // if FS != anti-FS, append also comb from anti-FS list (e.g. etac from (KS K- pi+) + (KS K+ pi-))
                        else 
                                CombineList(fPidList[14], auxmothpdg, auxamothpdg,  auxidx, auxaidx, true);     
                }
        }
        
        // create the final list using (or not using) the aux list
        CombineList(l, mothpdg, amothpdg, idx, aidx, cc);
                
        return l.GetLength();
}

virtual void PndSoftTriggerTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

void PndSoftTriggerTask::Exec ( Option_t *  opt )

virtual

Definition at line 816 of file softrig/PndSoftTriggerTask.cxx.

References RhoTuple::Column(), DetermineRecoilMode(), RhoTuple::DumpData(), f, fAllCands, fAnalysis, fApplyFullSelection, fChargedCands, fEcm, fEventShape, fEvtCount, fGammaMinE, FillEventShapeVarArray(), FillGlobalLists(), fIniP4, fMcTruth, fMode, fPidMult_025, fPocaVertexer, fPrimVtx, fPrimVtxQa, fQA, fQAEvent, fQAMc, fRecoilCnt, fRecoilMode, fRunNum, fSTTriggers, fTcaOnlineFilterInfo, fTrackMinP, PndAnalysis::GetEventInTask(), PndSoftTriggerLine::GetModeCode(), PndSoftTriggerLine::GetNTagged(), RhoVtxPoca::GetPocaVtx(), PndSoftTriggerLine::GetTagActive(), PndSoftTriggerLine::GetThreshold(), nmc, ntp, PndRhoTupleQA::qaEventShape(), PndRhoTupleQA::qaMcList(), PndRhoTupleQA::qaP4(), PndOnlineFilterInfo::SetNTag(), PndSoftTriggerLine::SetNTagged(), and TagMode().

{
        // *** prepare TCA for writing OnlineFilterInfo
        if ( fTcaOnlineFilterInfo->GetEntriesFast() != 0 ) fTcaOnlineFilterInfo->Delete();
        
        // *** some variables
        //int i=0,j=0, k=0; //[R.K. 01/2017] unused variable
        
        if (!(++fEvtCount%100)) cout << "[PndSoftTriggerTask] evt "<<fEvtCount<<endl;
        
        // *** read the next event
        fAnalysis->GetEventInTask();
        
        // *** fill all lists necessary for combinatorics
        FillGlobalLists();
        
        // *** go through MC truth list and codify the recoil mode (pi+-, pi0, eta, gamma, K+-)
        fRecoilMode  = DetermineRecoilMode(fRecoilCnt);
        
        // *** estimate primary vertex
        fPrimVtxQa = fPocaVertexer->GetPocaVtx(fPrimVtx, fChargedCands);
        
        // *** setup eventshape object
        PndEventShape evtShape(fAllCands,fIniP4,fGammaMinE,fTrackMinP);
        fEventShape = &evtShape;
        
        if (fApplyFullSelection) FillEventShapeVarArray();
        
        int tag_glob = 0, tag_glob_pre = 0, tag_pre = 0;

        PndOnlineFilterInfo* info=new ( (*fTcaOnlineFilterInfo)[0] ) PndOnlineFilterInfo();     
        
        // *** loop through all channels and tag
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
        {
                PndSoftTriggerLine *tl = it->second;    
                
                if (tl->GetThreshold()<=fEcm && tl->GetTagActive() )
                {
                        int ntag =  TagMode( tl , tag_pre);
                        tl->SetNTagged( ntag );
                        info->SetNTag(it->first, ntag);                      
                        tag_glob += ntag;
                        tag_glob_pre += tag_pre;
                }
        }
        
        Int_t tagged  = (tag_glob>0);
        Int_t taggedm = (tag_glob_pre>0);  // tagged by mass cut only
        
        // *** write common information 
        if (fQAEvent)
        {
                ntp->Column("ev",               (Int_t)         fEvtCount,        0);
                ntp->Column("run",      (Int_t)         fRunNum,          0);
                ntp->Column("mode",             (Int_t)         fMode,            0);
                ntp->Column("recmode",  (Int_t)     fRecoilMode,  0);
                ntp->Column("reccnt",   (Int_t)     fRecoilCnt,   0);

                fQA->qaP4("beam", fIniP4, ntp);
                ntp->Column("primvx",   (Float_t)   fPrimVtx.X(), 0.0f);
                ntp->Column("primvy",   (Float_t)   fPrimVtx.Y(), 0.0f);
                ntp->Column("primvz",   (Float_t)   fPrimVtx.Z(), 0.0f);
                ntp->Column("primvqa",  (Float_t)   fPrimVtxQa  , 0.0f);
                
                // write tags of individual modes
                for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
                {
                        PndSoftTriggerLine *tl = it->second;
                        if ( tl->GetThreshold()<=fEcm && tl->GetTagActive() ) 
                        {
                                //ntp->Column("tag_"+tl->GetName(), (Int_t) tl->GetNTagged(), 0);
                                ntp->Column(TString::Format("tag%d",tl->GetModeCode()), (Int_t) tl->GetNTagged(), 0);
                        }
                }

                ntp->Column("tagall",   (Int_t) tag_glob,               0);
                ntp->Column("tagpre",   (Int_t) tag_glob_pre,   0);
                ntp->Column("tag",              (Int_t) tagged,                 0);
                ntp->Column("tagm",             (Int_t) taggedm,                0);
                
                fQA->qaEventShape("es", fEventShape, ntp);
                
                // replace PID mult values from fEventShape by actual counts with individual algorithms
                ntp->Column("eslnpide", (Int_t)  fPidMult_025[0],               0 );
                ntp->Column("eslnpidmu",(Int_t)  fPidMult_025[1],               0 );
                ntp->Column("eslnpidpi",(Int_t)  fPidMult_025[2],               0 );
                ntp->Column("eslnpidk", (Int_t)  fPidMult_025[3],               0 );
                ntp->Column("eslnpidp", (Int_t)  fPidMult_025[4],               0 );
                
                ntp->DumpData();
        }
        
        // *** write out MC info
        if (fQAMc) {fQA->qaMcList("",fMcTruth,nmc); nmc->DumpData();}
        
}

void PndSoftTriggerTask::FillEventShapeVarArray ( )

private

Definition at line 1128 of file softrig/PndSoftTriggerTask.cxx.

References PndEventShape::Aplanarity(), PndEventShape::ChrgPSumCms(), PndEventShape::ChrgPSumLab(), PndEventShape::ChrgPtSumCms(), PndEventShape::Circularity(), fEventShape, PndEventShape::FoxWolfMomR(), fPidMult_025, fSTVarEvArray, i, PndEventShape::MultChrgPminLab(), PndEventShape::MultNeutEminLab(), PndEventShape::NCharged(), PndEventShape::NeutESumCms(), PndEventShape::NeutESumLab(), PndEventShape::NeutEtSumCms(), PndEventShape::NNeutral(), PndEventShape::NParticles(), PndEventShape::Planarity(), PndEventShape::PmaxCms(), PndEventShape::PmaxLab(), PndEventShape::PminCms(), PndEventShape::PminLab(), PndEventShape::PRapmax(), PndEventShape::Ptmax(), PndEventShape::Ptmin(), PndEventShape::PtSumLab(), PndEventShape::Sphericity(), PndEventShape::SumChrgPminCms(), and PndEventShape::Thrust().

Referenced by Exec().

{
        // *** Event Shape Vars

        // (  0) eslnpide        (  1) eslnpidmu       (  2) eslnpidpi       (  3) eslnpidk        (  4) eslnpidp           
        // (  5) esfw1           (  6) esfw2           (  7) esfw3           (  8) esfw4           (  9) esfw5
        // ( 10) esapl           ( 11) escir           ( 12) espla           ( 13) esnneut         ( 14) esnpart      
        // ( 15) esnchrg         ( 16) espmax          ( 17) espmaxl         ( 18) espmin          ( 19) espminl       
        // ( 20) esprapmax       ( 21) esptmax         ( 22) essumen         ( 23) essumenl        ( 24) essumetn       
        // ( 25) essumpc         ( 26) essumpc05       ( 27) essumpcl        ( 28) essumpt         ( 29) essumptc    
        // ( 30) esthr           ( 31) essph           ( 32) esptmin         ( 33) esncp10l        ( 34) esnne10l   
        
        // don't use PID mult values from fEventShape (based on AllCands and only one algo)
        for (int i=0;i<5;++i) fSTVarEvArray[i] = fPidMult_025[i]; 
        
        fSTVarEvArray[ 5] = fEventShape->FoxWolfMomR(1);
        fSTVarEvArray[ 6] = fEventShape->FoxWolfMomR(2);
        fSTVarEvArray[ 7] = fEventShape->FoxWolfMomR(3);
        fSTVarEvArray[ 8] = fEventShape->FoxWolfMomR(4) ;
        fSTVarEvArray[ 9] = fEventShape->FoxWolfMomR(5);
        
        fSTVarEvArray[10] = fEventShape->Aplanarity();
        fSTVarEvArray[11] = fEventShape->Circularity();
        fSTVarEvArray[12] = fEventShape->Planarity();
        fSTVarEvArray[13] = fEventShape->NNeutral();
        fSTVarEvArray[14] = fEventShape->NParticles();
        
        fSTVarEvArray[15] = fEventShape->NCharged();
        fSTVarEvArray[16] = fEventShape->PmaxCms();
        fSTVarEvArray[17] = fEventShape->PmaxLab();
        fSTVarEvArray[18] = fEventShape->PminCms();
        fSTVarEvArray[19] = fEventShape->PminLab();
        
        fSTVarEvArray[20] = fEventShape->PRapmax();
        fSTVarEvArray[21] = fEventShape->Ptmax();
        fSTVarEvArray[22] = fEventShape->NeutESumCms();
        fSTVarEvArray[23] = fEventShape->NeutESumLab();
        fSTVarEvArray[24] = fEventShape->NeutEtSumCms();
        
        fSTVarEvArray[25] = fEventShape->ChrgPSumCms();
        fSTVarEvArray[26] = fEventShape->SumChrgPminCms(0.5);
        fSTVarEvArray[27] = fEventShape->ChrgPSumLab();
        fSTVarEvArray[28] = fEventShape->PtSumLab();
        fSTVarEvArray[29] = fEventShape->ChrgPtSumCms();

        fSTVarEvArray[30] = fEventShape->Thrust();
        fSTVarEvArray[31] = fEventShape->Sphericity();
        fSTVarEvArray[32] = fEventShape->Ptmin();
        fSTVarEvArray[33] = fEventShape->MultChrgPminLab(1.0);
        fSTVarEvArray[34] = fEventShape->MultNeutEminLab(1.0);
        
}

void PndSoftTriggerTask::FillGlobalLists ( )

private

Definition at line 1018 of file softrig/PndSoftTriggerTask.cxx.

References RhoTuple::Column(), RhoCandList::Combine(), CreateKs0Cands(), RhoTuple::DumpData(), f, fabs(), fAlgoElectron, fAlgoKaon, fAlgoMuon, fAlgoPion, fAlgoProton, fAllCands, fAnalysis, fChargedCands, fEnergySel, fEtaCands, fEtaMean, fEtaPreSel, fEtaSel, fEtaSigma, fEvtCount, fGammaCands, PndAnalysis::FillList(), fIniP4, fIniPidCut, fKs0Cands, fMcTruth, fMode, fMomentumSel, fNeutralCands, fNsigAux, fPi0Cands, fPi0Mean, fPi0PreSel, fPi0Sel, fPi0Sigma, fPidList, fPidMult_025, fQA, fQAEta, fQAPi0, fRunNum, fTrackMinP, RhoCandList::GetLength(), i, Mass, MultPidProb(), neta, nks0, npi0, nsig, PndRhoTupleQA::qaP4(), PndRhoTupleQA::qaPi0(), RhoCandList::Select(), SelectPidProb(), RhoCandList::SetType(), and tag.

Referenced by Exec().

{
        int i;
        
        // *** fill standard lists
        fAnalysis->FillList( fMcTruth,                  "McTruth" );
        fAnalysis->FillList( fAllCands,                 "All"  ,    fAlgoElectron);
        fAnalysis->FillList( fChargedCands,             "Charged" , fAlgoElectron);
        fAnalysis->FillList( fNeutralCands,             "Neutral" );    
        
        // *** fill PID lists

        fAnalysis->FillList( fPidList[0],       "ElectronAllPlus",      fAlgoElectron );
        fAnalysis->FillList( fPidList[1],       "ElectronAllMinus", fAlgoElectron );
        
        fAnalysis->FillList( fPidList[2],       "MuonAllPlus",          fAlgoMuon );
        fAnalysis->FillList( fPidList[3],       "MuonAllMinus",         fAlgoMuon );
        
        fAnalysis->FillList( fPidList[4],       "PionAllPlus",          fAlgoPion );
        fAnalysis->FillList( fPidList[5],       "PionAllMinus",         fAlgoPion );
        
        fAnalysis->FillList( fPidList[6],       "KaonAllPlus",          fAlgoKaon );
        fAnalysis->FillList( fPidList[7],       "KaonAllMinus",         fAlgoKaon );
        
        fAnalysis->FillList( fPidList[8],       "ProtonAllPlus",        fAlgoProton );
        fAnalysis->FillList( fPidList[9],       "ProtonAllMinus",       fAlgoProton );

        // *** apply minimum PID cut if required
        if (fIniPidCut>0) for (i=0;i<10;++i) SelectPidProb(fPidList[i],i/2, fIniPidCut);
                
        // *** select on lists
        fGammaCands.Select(fNeutralCands, fEnergySel);
        
        // *** apply minimum momentum cut if required
        if (fTrackMinP>0) for (i=0;i<10;++i) fPidList[i].Select(fMomentumSel);
                
        // *** count PID multiplicities
        for (i=0;i<5;++i) fPidMult_025[i] = MultPidProb(fPidList[2*i], i, 0.25) + MultPidProb(fPidList[2*i+1], i, 0.25);

        // *** Ks reco
        CreateKs0Cands(nks0);
        
        // *** pi0 reco
        fPi0Cands.Combine(fGammaCands, fGammaCands);
        fPi0Cands.Select(fPi0PreSel);
        fPi0Cands.SetType(111);
        if (fQAPi0)
        {
                for (i=0; i<fPi0Cands.GetLength();++i)
                {
                        Float_t nsig = (Float_t) fabs(fPi0Cands[i]->Mass()-fPi0Mean)/fPi0Sigma;
                        Float_t tag  = nsig<fNsigAux;  // instead of fNsigTag, a fNsigAux cut is applied for pi0
                        
                        npi0->Column("tag",     (Int_t) tag,            0);
                        npi0->Column("nsig",    (Float_t) nsig,         0.0f);
                        npi0->Column("ev",      (Int_t) fEvtCount,      0);
                        npi0->Column("run",     (Int_t) fRunNum,        0);
                        npi0->Column("num",     (Int_t) i,                      0);
                        npi0->Column("mode",    (Int_t) fMode,          0);
                        npi0->Column("xmean", (Float_t) fPi0Mean,               0.0f);
                        npi0->Column("xsig",    (Float_t) fPi0Sigma,    0.0f);
                        
                        fQA->qaP4("beam", fIniP4, npi0);
                        fQA->qaPi0("x", fPi0Cands[i], npi0);
                        npi0->DumpData();
                }
        }
        fPi0Cands.Select(fPi0Sel);
        
        // *** era reco
        fEtaCands.Combine(fGammaCands, fGammaCands);
        fEtaCands.Select(fEtaPreSel);
        fEtaCands.SetType(221);
        if (fQAEta)
        {
                for (i=0; i<fEtaCands.GetLength();++i)
                {
                        Float_t nsig = (Float_t) fabs(fEtaCands[i]->Mass()-fEtaMean)/fEtaSigma;
                        Float_t tag  = nsig<fNsigAux;  // instead of fNsigTag, a fNsigAux cut is applied for eta
                        
                        neta->Column("tag",     (Int_t) tag,            0);
                        neta->Column("nsig",    (Float_t) nsig,         0.0f);
                        neta->Column("ev",      (Int_t) fEvtCount,      0);
                        neta->Column("run",     (Int_t) fRunNum,        0);
                        neta->Column("num",     (Int_t) i,                      0);
                        neta->Column("mode",    (Int_t) fMode,          0);
                        neta->Column("xmean", (Float_t) fEtaMean,               0.0f);
                        neta->Column("xsig",    (Float_t) fEtaSigma,    0.0f);
                        
                        fQA->qaP4("beam", fIniP4, neta);
                        fQA->qaPi0("x", fEtaCands[i], neta);
                        
                        neta->DumpData();
                }
        }
        fEtaCands.Select(fEtaSel);
        
        // *** copy all lists to indexed ones for combinatorics

        fPidList[10] = fGammaCands;
        
        fPidList[11] = fPi0Cands;
        fPidList[12] = fKs0Cands;
        fPidList[13] = fEtaCands;
}

void PndSoftTriggerTask::FillMassHisto ( TH1F *  h, RhoCandList &  l  )

private

Definition at line 469 of file tutorials/anatask/PndSoftTriggerTask.cxx.

References RhoCandList::GetLength(), and i.

{
        Int_t i=0;
        for (i=0;i<l.GetLength();++i)
        {
                h->Fill(l[i]->M());
        }
}

void PndSoftTriggerTask::FillVarArray ( RhoCandidate *  c, int  id, Bool_t  tmva = false  )

private

Definition at line 1392 of file softrig/PndSoftTriggerTask.cxx.

References fSTInputValues, fSTSelmap, fSTVarCandArray, GetVarValue(), and i.

Referenced by AcceptCandidate().

{
        int i=0;
        
        fSTVarCandArray.clear();   // reset cand var cache (used during fill of values)
        
        // fill all requested variables
        if (tmva) 
                for (i=0; i<fSTSelmap[mcode].tmvanvar; ++i) 
                        fSTInputValues[i] = GetVarValue(c, fSTSelmap[mcode].tmvavarid[i]);  // fill variables for tmva
        else
                for (i=0; i<fSTSelmap[mcode].ncut; ++i) 
                        fSTInputValues[i] = GetVarValue(c, fSTSelmap[mcode].varid[i]);      // fill variables for cut selection
}       

virtual void PndSoftTriggerTask::Finish ( )

virtual

void PndSoftTriggerTask::Finish ( )

virtual

Definition at line 915 of file softrig/PndSoftTriggerTask.cxx.

References fSTTriggers, RhoTuple::GetInternalTree(), PndSoftTriggerLine::GetRhoTuple(), neta, nks0, nmc, npi0, and ntp.

{       
        if (ntp) ntp->GetInternalTree()->Write();                         // overall info
        if (nks0) nks0->GetInternalTree()->Write();                   // Ks0 QA
        if (npi0) npi0->GetInternalTree()->Write();                       // pi0 QA
        if (neta) neta->GetInternalTree()->Write();                       // eta QA
        if (nmc)  nmc->GetInternalTree()->Write();                        // MC QA
        
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
        {
                PndSoftTriggerLine *tl = it->second;
                if (tl->GetRhoTuple())
                        tl->GetRhoTuple()->GetInternalTree()->Write();
        }
}

void PndSoftTriggerTask::GetAngles ( RhoCandidate *  c, double &  oang, double &  decang  )

private

Definition at line 1191 of file softrig/PndSoftTriggerTask.cxx.

References d0, RhoCandidate::Daughter(), RhoCandidate::NDaughters(), RhoCandidate::P3(), and RhoCandidate::P4().

Referenced by GetVarValue().

{
        oang   = -999.;
        decang = -999.;
        
        if (c->NDaughters()!=2) return;

        RhoCandidate *d0 = c->Daughter(0);
        RhoCandidate *d1 = c->Daughter(1);

        // opening angle lab
        oang = d0->P3().Angle(d1->P3());

        // decay angle
        TLorentzVector d_cms = d0->P4();
        d_cms.Boost(-(c->P4().BoostVector()));
        decang  = d_cms.Vect().Angle(c->P3());
}

double PndSoftTriggerTask::GetPocaVtx ( RhoCandidate *  c, double &  dist, double &  ctau  )

private

Definition at line 1213 of file softrig/PndSoftTriggerTask.cxx.

References fAnalysis, PndAnalysis::FillList(), fPocaVertexer, fPrimVtx, RhoCandList::GetLength(), RhoVtxPoca::GetPocaVtx(), RhoCandidate::M(), RhoCandidate::P(), and RhoCandList::RemoveFamily().

Referenced by GetVarValue().

{
        // *** simple vtx finder
        TVector3 vtx, altvtx, primvtx;
        double qavtx = fPocaVertexer->GetPocaVtx(vtx, c);
        
        // *** determine poca of rest of tracks
        RhoCandList l;
        fAnalysis->FillList(l, "Charged");
        
        // remove c's daughters from charged tracks list
        l.RemoveFamily(c);

        // if at least 2 tracks left, compute the poca vertex of the residual tracks
        if (l.GetLength()>1) fPocaVertexer->GetPocaVtx(altvtx, l);
        else altvtx = fPrimVtx; // else use the primary vtx found before
        
        dist = 999.;
        // *** if vertex of either all charged or residual tracks was found, compute distance and ctau
        if (altvtx.Mag()>0.) dist = (vtx-altvtx).Mag();
        
        ctau = dist*c->M()/c->P();

        return qavtx;
}

double PndSoftTriggerTask::GetVarValue ( RhoCandidate *  c, int  id  )

private

Definition at line 1301 of file softrig/PndSoftTriggerTask.cxx.

References BoostCms(), C(), c, cos(), RhoCandidate::Daughter(), fIniP4, fSTVarCandArray, fSTVarEvArray, GetAngles(), RhoCandidate::GetPidInfo(), GetPocaVtx(), i, RhoCandidate::M(), RhoCandidate::P(), RhoCandidate::P4(), and RhoCandidate::Pt().

Referenced by FillVarArray().

{       
        // previously filled event shape var?
        if (id>=0 && id<35) return fSTVarEvArray[id];
        
        // is variable already in cache?
        if (fSTVarCandArray.find(id) != fSTVarCandArray.end()) return fSTVarCandArray[id];
        
        // *** now compute variable values based on code 'id'
        
        // these are special cases with simple code numbers
        if (id==80) return (fIniP4 - c->P4()).M();         // mmiss
        if (id==82) return (c->M() - c->Daughter(0)->M()); // D* mdif
        if (id==99) return 1.0;                            // artificial value for tag>0 (= empty cut)
                
        // determine quantity and candidate from id code
        int quant = id%100;
        int cand  = id/100;
        int B = (id/10000)%10, C = (id/1000)%10, D = (id/100)%10;
        
        // find particle under consideration
        RhoCandidate *cnd=0;
        
        if (B==9) cnd = c;
        else if (C==9) cnd = c->Daughter(B);
        else if (D==9) cnd = c->Daughter(B)->Daughter(C);
        else cnd = c->Daughter(B)->Daughter(C)->Daughter(D);
        
        // coding for variable suffix (to translate variable names to code)
        //   E\F    |    0    |    1    |    2    |    3    |    4    |    5 
        // ----------------------------------------------------------------------
        //  1 (kin) |    p    |   pt    |  tht    |  pcm    |  thtcm  |
        //  2 (pid) |  pide   |  pidmu  |  pidpi  |  pidk   |  pidp   |  pidmax
        //  3 (vtx) |  pocqa  | pocdist | pocctau |         |         |
        //  4 (ang) |  oang   | decang  | cdecang |         |         |
        
        //double val = -999.; //[R.K. 01/2017] unused variable
        double pocdist, pocctau, oang, decang, cdecang, pidmax; //pocqa, //[R.K. 01/2017] unused variable
        TLorentzVector bl;
        
        switch (quant)
        {
                case 10 : // p
                        fSTVarCandArray[id] = cnd->P();  break;
                
                case 11 : // pt 
                        fSTVarCandArray[id] = cnd->Pt(); break;
                
                case 12 : // tht 
                        fSTVarCandArray[id] = cnd->P4().Theta(); break;
                
                case 13 : case 14 : // pcm, thtcm
                        bl = BoostCms(cnd->P4());
                        fSTVarCandArray[cand*100+13] = bl.P();
                        fSTVarCandArray[cand*100+14] = bl.Theta();
                        break;
                        
                case 20 : case 21 : case 22 : case 23 : case 24 :  // pide,..., pidp
                        fSTVarCandArray[id] = cnd->GetPidInfo(quant%10);
                        break;
                
                case 25 : // pidmax
                        pidmax = 0;
                        for (int i=0;i<5;++i) if (cnd->GetPidInfo(i)>pidmax) pidmax = cnd->GetPidInfo(i);
                        fSTVarCandArray[id] = pidmax;
                        break;
                
                case 30 : case 31 : case 32 : // pocqa, pocdist, pocctau
                        fSTVarCandArray[cand*100+30] = GetPocaVtx(cnd, pocdist, pocctau);
                        fSTVarCandArray[cand*100+31] = pocdist; 
                        fSTVarCandArray[cand*100+32] = pocctau; 
                        break;
                
                case 40 : case 41 : case 42 : // oang, decang, cdecang
                        GetAngles(cnd, oang, decang);
                        cdecang = cos(decang);
                        fSTVarCandArray[cand*100+40] = oang; 
                        fSTVarCandArray[cand*100+41] = decang; 
                        fSTVarCandArray[cand*100+42] = cdecang; 
                        break;
                default:
                        cout <<"[PndSoftTriggerTask] ***** wrong variable code '"<<quant<<"'"<<endl;
        }
        
        return fSTVarCandArray[id];
}

virtual InitStatus PndSoftTriggerTask::Init ( )

virtual

Virtual method Init

InitStatus PndSoftTriggerTask::Init ( )

virtual

Definition at line 250 of file softrig/PndSoftTriggerTask.cxx.

References STCutSet::cutval, fabs(), fAnalysis, fApplyFullSelection, fCfgFileName, fEcm, fEnergySel, fEtaMean, fEtaPreSel, fEtaQaMax, fEtaQaMin, fEtaSel, fEtaSigma, fGammaMinE, fKs0Mean, fKs0PreSel, fKs0QaMax, fKs0QaMin, fKs0Sel, fKs0Sigma, fMomentumSel, fNsigAux, fPbarMom, fPhotosMax, fPhotosThresh, fPi0Mean, fPi0PreSel, fPi0QaMax, fPi0QaMin, fPi0Sel, fPi0Sigma, fPocaVertexer, fQA, fQAEta, fQAEvent, fQAKs0, fQAMc, fQAPi0, fRootManager, fSTencode, fSTModeIndex, fSTOps, fSTSelmap, fSTTriggers, fTcaOnlineFilterInfo, fTrackMinP, fVerbose, RhoTuple::GetInternalTree(), PndSoftTriggerLine::GetRhoTuple(), PndSoftTriggerLine::GetThreshold(), i, if(), PndSoftTriggerLine::Init(), PndAnalysis::McMatchAllowPhotos(), STCutSet::ncut, neta, nks0, nmc, npi0, ntp, STCutSet::op, PndSoftTriggerLine::Print(), ReadConfiguration(), STCutSet::tmvacut, STCutSet::tmvanvar, STCutSet::tmvaread, STCutSet::tmvavarid, and STCutSet::varid.

Referenced by qa_softtrig().

{               
        fRootManager = FairRootManager::Instance();
        
        // Register TCA for tagging info
    fTcaOnlineFilterInfo = new TClonesArray ( "PndOnlineFilterInfo" );
        if (fRootManager) 
                fRootManager->Register ( "OnlineFilterInfo","PndOnlineFolder", fTcaOnlineFilterInfo, kTRUE );
        
        // *** initialize analysis object
        fAnalysis = new PndAnalysis();
        if (fPhotosMax) fAnalysis->McMatchAllowPhotos(fPhotosMax,fPhotosThresh);
        
        // *** RhoTuple QA helper
        fQA = new PndRhoTupleQA(fAnalysis,fPbarMom);
        
        // *** Save current gDirectory
        TDirectory *dir = gDirectory;
        fRootManager->GetOutFile()->cd();
        
        // *** create ntuple
        if (fQAEvent) { ntp = new RhoTuple("ntpev","Soft Trigger Common"); ntp->GetInternalTree()->SetDirectory(gDirectory);}
        if (fQAKs0)   { nks0 = new RhoTuple("nks0","K_S -> pi+ pi-"); nks0->GetInternalTree()->SetDirectory(gDirectory);}
        if (fQAPi0)   { npi0 = new RhoTuple("npi0","pi0 -> gamma gamma"); npi0->GetInternalTree()->SetDirectory(gDirectory);}
        if (fQAEta)   { neta = new RhoTuple("neta","eta -> gamma gamma"); neta->GetInternalTree()->SetDirectory(gDirectory);}
        if (fQAMc)    { nmc  = new RhoTuple("nmc", "MC info"); nmc->GetInternalTree()->SetDirectory(gDirectory);}
        
        
        // *** create mass pre selectors for QA (formular takes into account RhoSelector definition mean +- win/2
        fPi0PreSel   = new RhoMassParticleSelector("pi0PreSel",  (fPi0QaMax + fPi0QaMin)/2.0,   fPi0QaMax - fPi0QaMin );  
        fEtaPreSel   = new RhoMassParticleSelector("etaPreSel",  (fEtaQaMax + fEtaQaMin)/2.0,   fEtaQaMax - fEtaQaMin ); 
        fKs0PreSel   = new RhoMassParticleSelector("Ks0PreSel",  (fKs0QaMax + fKs0QaMin)/2.0,   fKs0QaMax - fKs0QaMin);         
        
        // *** number of sigmas deviation for tag
        //double tagNumSig   = fNsigTag;

        // *** create final selectors for pi0, eta, KS 
        fPi0Sel     = new RhoMassParticleSelector("pi0Sel",   fPi0Mean, fPi0Sigma*2.0*fNsigAux);
        fEtaSel     = new RhoMassParticleSelector("etaSel",   fEtaMean, fEtaSigma*2.0*fNsigAux);
        fKs0Sel     = new RhoMassParticleSelector("Ks0Sel",   fKs0Mean, fKs0Sigma*2.0*fNsigAux);

        // *** basic selectors for preselection
        fMomentumSel = new RhoMomentumParticleSelector("PSel",50.+fTrackMinP,100.);
        fEnergySel   = new RhoEnergyParticleSelector("ESel",50.+fGammaMinE,100.);
        
        // *** the poca vertexer
        fPocaVertexer = new RhoVtxPoca();
                
        // *** read selection from configuration file
        fSTencode.clear();
        if (fCfgFileName!="" && fApplyFullSelection) ReadConfiguration();

        // *** set mode index for current beam momentum
        fSTModeIndex = 0;
        double diff = 1000.;
        
        for (size_t i=0; i<fSTencode.size(); ++i) 
        {
                double en = (double)fSTencode[i]/100.;
                if (fabs(fEcm-en)<diff)
                {
                        diff=fabs(fEcm-en);
                        fSTModeIndex = i;
                }
        }
        
        if (fApplyFullSelection && diff>1.0) cout <<"[PndSoftTriggerTask::Init] **** WARNING: Energy of best matching selection differs "<<diff<<" GeV from current energy!"<<endl;

        // *** initialize triggers
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
        {
                PndSoftTriggerLine *tl = it->second;
                
                if (tl->GetThreshold()<=fEcm)
                {
                        tl->Init();
                        if (tl->GetRhoTuple()) tl->GetRhoTuple()->GetInternalTree()->SetDirectory(gDirectory);
                        
                        if (fVerbose>0)
                        {
                                cout <<"*** MODE: "<<it->first<<endl;
                                tl->Print();
                                cout <<endl<<endl;
                        }
                }
        }
        
        // *** Resore original gDirectory
        dir->cd();

        if (fVerbose>0)
        {
                // print selection setup
                cout <<"[PndSoftTriggerTask] **** Selection setup:"<<endl;
                for (std::map<int,STCutSet>::iterator it=fSTSelmap.begin(); it!=fSTSelmap.end(); ++it)
                {
                        std::cout << it->first << " => ";
                        STCutSet cs = it->second;
                        for (int i=0;i<cs.ncut;++i) cout <<"v["<<cs.varid[i]<<"]"<<fSTOps[cs.op[i]]<<cs.cutval[i]<<"  ";
                        cout <<endl;
                        if (cs.tmvaread)
                        {
                                cout<<"  TMVA cut:"<<cs.tmvacut<<"  vars:";
                                for (int i=0;i<cs.tmvanvar;++i) cout <<cs.tmvavarid[i]<<" ";
                                cout<< endl;
                        }
                }
        }

        return kSUCCESS;
}

void PndSoftTriggerTask::McMatchAllowPhotos ( int  maxn = 1, double  thresh = 0.05  )

inline

Definition at line 116 of file softrig/PndSoftTriggerTask.h.

References fPhotosMax, and fPhotosThresh.

Referenced by ana_jpsi_task_ST().

{fPhotosMax=maxn;fPhotosThresh=thresh;}  

int PndSoftTriggerTask::MultPidProb ( RhoCandList &  l, int  pididx, double  prob  )

private

Definition at line 513 of file softrig/PndSoftTriggerTask.cxx.

References RhoCandList::GetLength(), and mult.

Referenced by FillGlobalLists().

{
        int mult = 0;
        
        if (pididx>=0 && pididx<5)
                for (int ii=0;ii<l.GetLength();++ii)
                        if (l[ii]->GetPidInfo(pididx)>=prob) ++mult;
                                
        return mult;
}

void PndSoftTriggerTask::PrintList ( RhoCandList &  l, int  max = 10  )

private

Definition at line 511 of file tutorials/anatask/PndSoftTriggerTask.cxx.

References RhoCandList::GetLength(), and max().

{
        int N=l.GetLength();
        if (N>max) N=max;
        for (int j=0;j<N;++j)
        {
                cout <<*(l[j])<<" PDG:"<<l[j]->PdgCode()<<" MC pointer:"<<l[j]->GetMcTruth()<<endl;
        }
        cout <<endl;
}

bool PndSoftTriggerTask::ReadConfiguration ( )

private

Definition at line 597 of file softrig/PndSoftTriggerTask.cxx.

References CodeVariable(), cuts, STCutSet::cutval, fApplyFullSelection, fCfgFileName, file, fSTenames, fSTencode, fSTInputValues, fSTNEvVars, fSTOps, fSTSelmap, fSTVarmap, fVerbose, i, if(), name, STCutSet::ncut, ok, STCutSet::op, SplitString(), STCutSet::tmvacut, STCutSet::tmvameth, STCutSet::tmvanvar, STCutSet::tmvaread, STCutSet::tmvavarid, TString, val, and STCutSet::varid.

Referenced by Init().

{
        for (int i=0;i<fSTNEvVars;++i) fSTVarmap[fSTenames[i]] = i;

        std::ifstream file(fCfgFileName.Data());
        
        if (!file.is_open()) 
        {
                cout <<"[PndSoftTriggerTask] **** Unable to open selection file "<<fCfgFileName.Data()<<endl;
                fApplyFullSelection=false;
        }
        else cout <<"[PndSoftTriggerTask] **** Reading selection from "<<fCfgFileName.Data()<<endl;
        
        char line[1500];
        
        TString toks[200];
        TString cuts[200];
        TString toks2[200];

        // loop through file line by line
        while (!file.eof())
        {
                file.getline(line,1499);
                TString sline(line);
                
                // remove tabs everywhere
                sline.ReplaceAll("\t","");

                // if lines contains a '#' somewhere, cut string from that position (=comment)
                if (sline.Contains("#")) sline = sline(0,sline.Index("#")-1);
                
                // remove whitespace at begin and end
                sline = sline.Strip(TString::kBoth);
                
                if (sline=="") continue;
                
                // split the line into tokens; token 0 ist the mode code, token 1 is the complete cut string
                // e.g. '38400 : eslnpidp>0&&abs(pcm-2.105)<0.695&&esthr>0.9&&pt>0.8'
                int N1 = SplitString(sline, ":", toks,200);
                
                // if N==2, only simple cut; N==3 includes TMVA selector
                if (N1<2 || N1>3)  {cout <<"invalid line: "<<sline.Data()<<endl; continue;}

                // extract the mode code by converting to integer
                int mcode = toks[0].Atoi();
                int en    = mcode / 1000; // energy prefix (e.g. 550 for sqs=5.5GeV)
                
                // if new energy, store the energy prefix in fSTencode
                if ( std::find(fSTencode.begin(), fSTencode.end(), en) == fSTencode.end() ) fSTencode.push_back(en);
                
                // now split the cut string into single cuts
                int N2 = SplitString(toks[1],"&&",cuts,200);
                
                // replace windows cuts 'abs(<name>-x)<y' by 2 single cuts
                for (int i=0;i<N2;++i) 
                {
                        cuts[i].ReplaceAll(" ","");
                        if (cuts[i]=="tag") cuts[i]="tag>0";
                        
                        // is cut a window cut?
                        if (cuts[i].BeginsWith("abs"))
                        {
                                int pos1 = cuts[i].First('-'); // position of '-'
                                int pos2 = cuts[i].Last(')');  // position of closing abs bracket 
                                int pos3 = cuts[i].Last('<');  // position of operator '<'
                                
                                TString name = TString(cuts[i](4,pos1-4));  // name is from first char after 'abs(' to '-' sign
                                double val = TString(cuts[i](pos1+1,(pos2-pos1)-1)).Atof(); // val to cut on is from '-' to closing bracket ')'
                                double win = TString(cuts[i](pos3+1,1000)).Atof();  // window size is from '<' to end
                                cuts[i]= TString::Format("%s>%f",name.Data(), val-win);      // construct the 1st and
                                cuts[N2++] = TString::Format("%s<%f",name.Data(), val+win);   // 2nd single cut
                        }
                }
                
                // *** now encode the string cuts into sets of (variable index, operator, cutvalue)
                // *** and put the cut sets into a map with key being the mode code
                
                STCutSet cs;
                cs.ncut = N2;
                
                cs.tmvaread = 0;
                cs.tmvanvar = 0;
                
                bool ok = true;  // if a variable is not know switch to false
                
                for (int i=0;i<N2;++i)
                {
                        int op=-1, j=0;
                        while (j<5 && op==-1) 
                        {
                                if (cuts[i].Contains(fSTOps[j])) 
                                {
                                        SplitString(cuts[i],fSTOps[j],toks2,5); 
                                        op = j;
                                }
                                j++;
                        }

                        if (op<0) // fail
                        {
                                i=N2+1;
                                ok=false;
                        }
                        
                        // event shape variable?
                        if (fSTVarmap.find(toks2[0]) != fSTVarmap.end())
                        {
                                cs.varid[i]  = fSTVarmap[toks2[0]];
                                cs.op[i]     = op;
                                cs.cutval[i] = toks2[1].Atof();
                        }
                        // candidate variable?
                        else if (CodeVariable(toks2[0])>0)
                        {
                                cs.varid[i]  = CodeVariable(toks2[0]);
                                cs.op[i]     = op;
                                cs.cutval[i] = toks2[1].Atof();
                                
                        }
                        else // fail
                        {
                                cout <<"[PndSoftTriggerTask] **** Unmapped cut variable: "<<toks2[0].Data()<<". Skipping mode "<<mcode<<"."<<endl;
                                i=N2+1;
                                ok=false;
                        }
                }
                
                // *** do we have a TMVA Selector defined?
                if (N1==3)
                {
                        // now split the TMVA string into single tokens
                        int N3 = SplitString(toks[2]," ",toks2,200);
                        
                        if (N3<4) // at least four elements expected; fail
                        {
                                cout <<"[PndSoftTriggerTask] **** Invalid TMVA configuration for mode "<<mcode<<". Skipping."<<endl;
                                ok = false;
                        }
                        
                        if (ok)
                        {
                                // first token is the weightfile name without leading path and ending '.weights.xml'
                                TString wfile   = TString(getenv("VMCWORKDIR"))+"/softrig/weights/"+toks2[0]+".weights.xml";
                                cs.tmvameth = toks2[0](toks2[0].Length()-3,3); 
                                
                                // last token is the cut to the TMVA output
                                cs.tmvacut  = toks2[N3-1].Atof();
                                
                                // tokens inbetween are the variables
                                cs.tmvanvar = N3 - 2;
                                
                                // create the reader
                                cs.tmvaread = new TMVA::Reader("Silent");
                                
                                for (int i=1; i<N3-1; ++i)
                                {
                                        if (fSTVarmap.find(toks2[i]) != fSTVarmap.end()) cs.tmvavarid[i-1] = fSTVarmap[toks2[i]];
                                        else if (CodeVariable(toks2[i])>0) cs.tmvavarid[i-1] = CodeVariable(toks2[i]);
                                        else // fail 
                                        {
                                                cout <<"[PndSoftTriggerTask] **** Unmapped TMVA variable: "<<toks2[i].Data()<<". Skipping mode "<<mcode<<"."<<endl;
                                                i=N2+1;
                                                ok=false;
                                        }
                                        if (ok) cs.tmvaread->AddVariable(toks2[i], &fSTInputValues[i-1]);
                                }
                                if (ok) cs.tmvaread->BookMVA(cs.tmvameth,wfile);
                        }
                }
                
                if (ok) fSTSelmap[mcode] = cs;
                
                if (fVerbose)
                {
                        cout <<"mcode="<<mcode<<endl;
                        STCutSet cts = fSTSelmap[mcode];
                        cout <<"  Ncut="<<cts.ncut<<" (";
                        for (int k=0;k<cts.ncut;++k) cout <<cts.varid[k]<<" ";
                        cout <<")"<<endl;
                        if (cts.tmvaread!=0)
                        {
                                cout <<"  TMVAMeth="<<cts.tmvameth<<"; Nvar="<<cts.tmvanvar<<"; Cut="<<cts.tmvacut<<"; TMVAobject="<<cts.tmvaread<<" (";                        
                                for (int k=0;k<cts.tmvanvar;++k) cout <<cts.tmvavarid[k]<<" ";
                                cout <<")"<<endl;
                        }
                }
        }
        return true;
}

bool PndSoftTriggerTask::ReadTriggerLines ( )

private

Definition at line 530 of file softrig/PndSoftTriggerTask.cxx.

References fApplyFullSelection, fEcm, file, fSTTriggers, fTriggerFileName, mean, mode, PndSoftTriggerLine::SetMeanSigma(), PndSoftTriggerLine::SetQAMassWindow(), PndSoftTriggerLine::SetTagActive(), PndSoftTriggerLine::SetTagNSig(), PndSoftTriggerLine::SetThreshold(), PndSoftTriggerLine::SetWriteQA(), SplitString(), and TString.

Referenced by PndSoftTriggerTask().

{
        std::ifstream file(fTriggerFileName.Data());
        if (!file.is_open()) 
        {
                cout <<"[PndSoftTriggerTask] **** Unable to open trigger file "<<fTriggerFileName.Data()<<endl;
                fApplyFullSelection=false;
        }
        else cout <<"[PndSoftTriggerTask] **** Reading trigger lines from "<<fTriggerFileName.Data()<<endl;
        
        char line[500];
        TString toks[13];
        
        // loop through file line by line
        while (!file.eof())
        {
                file.getline(line,499);
                TString sline(line);

                // empty lines or comments are skipped (a comment has to have a # as _first_ char)
                if (sline.BeginsWith("#")||sline=="") continue;  
        
                // split the line into tokens, has to be
                // 0-mode, 1-name, 2-decay pattern, 3-prefix, 4-ntpname, 5-qamin, 6-qamax, 7-mean, 8-sigma, 9-threshold, 10-nsig, 11-fqa, 12-ftag
                int N = SplitString(sline, ":", toks,13);
                if (N!=13)  {cout <<"[PndSoftTriggerTask] **** Invalid trigger line: "<<sline.Data()<<endl; continue;}
                //if (fVerbose>1) for (int i=0;i<N;++i) cout <<toks[i]<<":"; cout <<endl; 
                
                int mode = toks[0].Atoi();
                PndSoftTriggerLine *tl = new PndSoftTriggerLine(mode, toks[1], toks[2], toks[3], toks[4]);

                double mean=0, qamin = toks[5].Atof(), qamax = toks[6].Atof();
                
                // check whether mean = Ecm
                if (toks[7]=="Ecm") 
                        mean = fEcm;
                else 
                        mean = toks[7].Atof();
                
                // check whether qa mass window values are relative to mean 
                if (toks[5].BeginsWith("-")) qamin = mean + toks[5].Atof(); // sum, since toks[5].Atof()<0 due to '-' sign
                if (toks[6].BeginsWith("+")) qamax = mean + toks[6].Atof();
                
                tl->SetQAMassWindow(qamin, qamax);
                tl->SetMeanSigma(mean, toks[8].Atof());
                tl->SetThreshold(toks[9].Atof());
                tl->SetTagNSig(toks[10].Atof());
                tl->SetWriteQA(toks[11].Atoi());
                tl->SetTagActive(toks[12].Atoi());
                
                fSTTriggers[mode] = tl;
        }
        
        cout << "[PndSoftTriggerTask] **** Found "<<fSTTriggers.size()<<" trigger definitions: ";
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
                cout << it->second->GetName()<<" ";
        cout <<endl;
        
        return true;
}

int PndSoftTriggerTask::RemoveDoubles ( RhoCandList &  l, double  limit = 0.0001  )

private

Definition at line 479 of file tutorials/anatask/PndSoftTriggerTask.cxx.

References fabs(), RhoCandList::GetLength(), i, limit, RhoCandList::Remove(), and sim().

{
        int i,j;
        int rem=0;

        //cout <<"********* before*****"<<endl<<l<<endl;
                
        for (i=0;i<l.GetLength();++i)
        {
                bool sim=false;
                for (j=i+1;j<l.GetLength();++j)
                {
                  if(!sim && 
                     fabs(l[i]->Px()-l[j]->Px())<limit &&
                     fabs(l[i]->Py()-l[j]->Py())<limit &&
                     fabs(l[i]->Pz()-l[j]->Pz())<limit &&
                     fabs(l[i]->E()-l[j]->E())<limit)
                    {
                      sim=true;
                    }
                      //if (similarity(l[i].P4(), l[j]->P4())<limit){ sim=true; j=10000;}
                }
                if (sim) {l.Remove(l[i--]); rem++;}
        }
        //cout<<"1:" <<rem<<endl;
        //cout <<"********* after *****"<<endl<<l<<endl;
        
        return rem;
}

int PndSoftTriggerTask::SelectPdgCode ( RhoCandList &  mct, RhoCandList &  l  )

private

Definition at line 542 of file tutorials/anatask/PndSoftTriggerTask.cxx.

References RhoCandList::GetLength(), RhoCandidate::PdgCode(), and RhoCandList::Remove().

{
        int removed = 0;
        int pdgcode=0;
        int nmct = mct.GetLength();
        
        //PndMcTruthMatch mcm;
        
        if (l.GetLength()>0) pdgcode = l[0]->PdgCode();
        
        for (int ii=l.GetLength();ii>=0;--ii)
        {
                RhoCandidate* mccand = l[ii]->GetMcTruth();
        
                if (mccand)
          if(mccand->PdgCode()==pdgcode)
//      if (mcm.MctMatch(l[ii],mct))
            continue; //don't remove if good mctruth
                l.Remove(l[ii]);
                removed++;
        }
        return removed;
}

void PndSoftTriggerTask::SelectPid ( int  type, int  pdg, int  chrg, RhoCandList &  l, RhoCandList &  lpid, double  cut = 0.2  )

private

Definition at line 524 of file tutorials/anatask/PndSoftTriggerTask.cxx.

References RhoCandList::Cleanup(), cut, fabs(), RhoCandList::GetLength(), and RhoCandList::Put().

{
        lpid.Cleanup();
                        
        for (int j=0;j<l.GetLength();++j)
        {
                if (fabs(l[j]->Charge()-double(chrg))<0.001 && l[j]->GetPidInfo(type)>=cut)
                {
          //RhoCandidate c(l[j]);
                        //c.SetMass(pdgmass[type]);
                        lpid.Put(l[j]);
                }               
        }
        for (int j=0;j<lpid.GetLength();++j) lpid[j]->SetType(pdg);
}

int PndSoftTriggerTask::SelectPidProb ( RhoCandList &  l, int  pididx, double  cut  )

private

Definition at line 496 of file softrig/PndSoftTriggerTask.cxx.

References cut, RhoCandList::GetLength(), and RhoCandList::Remove().

Referenced by FillGlobalLists().

{
        int removed = 0;
        
        if (pididx>=0 && pididx<5)
                for (int ii=l.GetLength()-1;ii>=0;--ii)
                        if (l[ii]->GetPidInfo(pididx)<cut )
                        {
                                l.Remove(l[ii]);
                                removed++;
                        }
        
        return removed;
}

int PndSoftTriggerTask::SelectTruePid ( RhoCandList &  l )

private

Definition at line 478 of file softrig/PndSoftTriggerTask.cxx.

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

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

void PndSoftTriggerTask::SetAuxNSig ( double  nsig )

inline

Definition at line 66 of file softrig/PndSoftTriggerTask.h.

References fNsigAux, and nsig.

{fNsigAux = nsig;}  // n sigma cut for KS, pi0, eta

void PndSoftTriggerTask::SetConfigurationFile ( TString  fname )

inline

Definition at line 55 of file softrig/PndSoftTriggerTask.h.

References fCfgFileName.

Referenced by ana_jpsi_task_ST(), SetFastSimDefaults(), and SetFullSimDefaults().

{fCfgFileName = fname;}

void PndSoftTriggerTask::SetDstMDiffCut ( double  cut )

inline

Definition at line 67 of file softrig/PndSoftTriggerTask.h.

References cut, and fDstMDiffCut.

Referenced by SetFastSimDefaults(), and SetFullSimDefaults().

{fDstMDiffCut = cut;} // cut on delta M(D*, D); default = 10. GeV/c^2

void PndSoftTriggerTask::SetEtaQASelection ( double  min, double  max  )

inline

Definition at line 72 of file softrig/PndSoftTriggerTask.h.

References fEtaQaMax, fEtaQaMin, max(), and min().

{fEtaQaMin=min; fEtaQaMax=max;}

void PndSoftTriggerTask::SetEtaSignalParams ( double  mean, double  sigma  )

inline

Definition at line 79 of file softrig/PndSoftTriggerTask.h.

References fEtaMean, fEtaSigma, mean, and sigma.

Referenced by ana_jpsi_task_ST(), qa_softtrig(), SetFastSimDefaults(), and SetFullSimDefaults().

{fEtaMean=mean; fEtaSigma=sigma;}

void PndSoftTriggerTask::SetFastSimDefaults

(

)

Definition at line 203 of file softrig/PndSoftTriggerTask.cxx.

References ApplyFullSelection(), SetConfigurationFile(), SetDstMDiffCut(), SetEtaSignalParams(), SetGammaMinE(), SetInitialPidCut(), SetKs0SignalParams(), SetPi0SignalParams(), SetPidAlgoAll(), SetQAAll(), SetTagAll(), SetTrackMinP(), and TString.

Referenced by prod_fsim().

{
        TString selectioncfg = TString(gSystem->Getenv("VMCWORKDIR"))+"/softrig/selection_fsim_dec2014.cfg"; 
        
        SetConfigurationFile(selectioncfg);
        ApplyFullSelection(1);  // apply selection defined in 'TString selectioncfg'
        
        SetPi0SignalParams(0.134, 0.0035);  // set parameters for pi0
        SetKs0SignalParams(0.497, 0.0055);  // set parameters for KS
        SetEtaSignalParams(0.549, 0.0088);  // set parameters for eta
        
        SetGammaMinE(0.10);             // global energy pre-cut for neutrals 
        SetTrackMinP(0.10);             // global momentum pre-cut for charged  
        SetInitialPidCut(0.1);  // global PID pre-cut for charged       
        SetDstMDiffCut(0.1);    // special cut on D*-D mass difference (to reduce comb and output file size)
        
        SetPidAlgoAll("PidChargedProbability");
        
        SetTagAll(true);                // tag all modes        
        SetQAAll(false);        // don't write any QA tuple     
}

void PndSoftTriggerTask::SetFullSimDefaults

(

)

Definition at line 226 of file softrig/PndSoftTriggerTask.cxx.

References ApplyFullSelection(), SetConfigurationFile(), SetDstMDiffCut(), SetEtaSignalParams(), SetGammaMinE(), SetInitialPidCut(), SetKs0SignalParams(), SetPi0SignalParams(), SetPidAlgoAll(), SetQAAll(), SetTagAll(), SetTrackMinP(), and TString.

{
        TString selectioncfg = TString(gSystem->Getenv("VMCWORKDIR"))+"/softrig/selection_full_jan2015.cfg"; 
        
        SetConfigurationFile(selectioncfg);
        ApplyFullSelection(1);  // apply selection defined in 'TString selectioncfg'
        
    SetPi0SignalParams(0.134, 0.0049);
    SetKs0SignalParams(0.497, 0.0085);
    SetEtaSignalParams(0.549, 0.0092);

    SetGammaMinE(0.10);        // global energy pre-cut for neutrals
    SetTrackMinP(0.10);        // global momentum pre-cut for charged
    SetInitialPidCut(0.1);            // global PID pre-cut for charged
    SetDstMDiffCut(0.1);            // special cut on D*-D mass difference (to reduce comb and output file size)

        SetPidAlgoAll("PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc;PidAlgoStt;PidAlgoMdtHardCuts;PidAlgoMvd");

        SetTagAll(true);                // tag all modes        
        SetQAAll(false);        // don't write any QA tuple     
}

void PndSoftTriggerTask::SetGammaMinE ( double  min )

inline

Definition at line 84 of file softrig/PndSoftTriggerTask.h.

References fGammaMinE, and min().

Referenced by ana_jpsi_task_ST(), qa_softtrig(), SetFastSimDefaults(), and SetFullSimDefaults().

{fGammaMinE = min;}

void PndSoftTriggerTask::SetInitialPidCut ( double  cut )

inline

Definition at line 86 of file softrig/PndSoftTriggerTask.h.

References cut, and fIniPidCut.

Referenced by ana_jpsi_task_ST(), SetFastSimDefaults(), and SetFullSimDefaults().

{fIniPidCut = cut;}

void PndSoftTriggerTask::SetKs0QASelection ( double  min, double  max  )

inline

Definition at line 71 of file softrig/PndSoftTriggerTask.h.

References fKs0QaMax, fKs0QaMin, max(), and min().

{fKs0QaMin=min; fKs0QaMax=max;}

void PndSoftTriggerTask::SetKs0SignalParams ( double  mean, double  sigma  )

inline

Definition at line 78 of file softrig/PndSoftTriggerTask.h.

References fKs0Mean, fKs0Sigma, mean, and sigma.

Referenced by ana_jpsi_task_ST(), qa_softtrig(), SetFastSimDefaults(), and SetFullSimDefaults().

{fKs0Mean=mean; fKs0Sigma=sigma;}

virtual void PndSoftTriggerTask::SetParContainers ( )

privatevirtual

Get parameter containers

void PndSoftTriggerTask::SetParContainers ( )

privatevirtual

Definition at line 805 of file softrig/PndSoftTriggerTask.cxx.

References run.

{
  // Get run and runtime database
  FairRun* run = FairRun::Instance();
  if ( ! run ) Fatal("SetParContainers", "No analysis run");
}

void PndSoftTriggerTask::SetPi0QASelection ( double  min, double  max  )

inline

Definition at line 70 of file softrig/PndSoftTriggerTask.h.

References fPi0QaMax, fPi0QaMin, max(), and min().

{fPi0QaMin=min; fPi0QaMax=max;}

void PndSoftTriggerTask::SetPi0SignalParams ( double  mean, double  sigma  )

inline

Definition at line 77 of file softrig/PndSoftTriggerTask.h.

References fPi0Mean, fPi0Sigma, mean, and sigma.

Referenced by ana_jpsi_task_ST(), qa_softtrig(), SetFastSimDefaults(), and SetFullSimDefaults().

{fPi0Mean=mean; fPi0Sigma=sigma;}

void PndSoftTriggerTask::SetPidAlgoAll

(

TString

algo

)

Definition at line 400 of file softrig/PndSoftTriggerTask.cxx.

References SetPidAlgoElectron(), SetPidAlgoKaon(), SetPidAlgoMuon(), SetPidAlgoPion(), and SetPidAlgoProton().

Referenced by ana_jpsi_task_ST(), SetFastSimDefaults(), and SetFullSimDefaults().

{
        SetPidAlgoElectron(algo);
        SetPidAlgoMuon(algo);
        SetPidAlgoPion(algo); 
        SetPidAlgoKaon(algo);
        SetPidAlgoProton(algo);
}

void PndSoftTriggerTask::SetPidAlgoElectron ( TString  algo )

inline

Definition at line 91 of file softrig/PndSoftTriggerTask.h.

References fAlgoElectron.

Referenced by SetPidAlgoAll().

{fAlgoElectron=algo;}

void PndSoftTriggerTask::SetPidAlgoKaon ( TString  algo )

inline

Definition at line 94 of file softrig/PndSoftTriggerTask.h.

References fAlgoKaon.

Referenced by SetPidAlgoAll().

{fAlgoKaon=algo;}

void PndSoftTriggerTask::SetPidAlgoMuon ( TString  algo )

inline

Definition at line 92 of file softrig/PndSoftTriggerTask.h.

References fAlgoMuon.

Referenced by SetPidAlgoAll().

{fAlgoMuon=algo;}

void PndSoftTriggerTask::SetPidAlgoPion ( TString  algo )

inline

Definition at line 93 of file softrig/PndSoftTriggerTask.h.

References fAlgoPion.

Referenced by SetPidAlgoAll().

{fAlgoPion=algo;}

void PndSoftTriggerTask::SetPidAlgoProton ( TString  algo )

inline

Definition at line 95 of file softrig/PndSoftTriggerTask.h.

References fAlgoProton.

Referenced by SetPidAlgoAll().

{fAlgoProton=algo;}

void PndSoftTriggerTask::SetQAAll

(

bool

qa = true

)

Definition at line 426 of file softrig/PndSoftTriggerTask.cxx.

References fSTTriggers, SetQAEta(), SetQAEvent(), SetQAKs0(), SetQAMc(), SetQAMode(), and SetQAPi0().

Referenced by ana_jpsi_task_ST(), prod_fsim(), SetFastSimDefaults(), and SetFullSimDefaults().

{       
        SetQAPi0(qa);
        SetQAEta(qa);
        SetQAKs0(qa);
        SetQAEvent(qa);
        SetQAMc(qa);
        
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it) 
                SetQAMode(it->first, qa);
}

void PndSoftTriggerTask::SetQAEta ( bool  qa = true )

inline

Definition at line 104 of file softrig/PndSoftTriggerTask.h.

References fQAEta.

Referenced by SetQAAll().

{ fQAEta = qa;}

void PndSoftTriggerTask::SetQAEvent ( bool  qa = true )

inline

Definition at line 109 of file softrig/PndSoftTriggerTask.h.

References fQAEvent.

Referenced by ana_jpsi_task_ST(), prod_fsim(), and SetQAAll().

{ fQAEvent = qa;}

void PndSoftTriggerTask::SetQAKs0 ( bool  qa = true )

inline

Definition at line 105 of file softrig/PndSoftTriggerTask.h.

References fQAKs0.

Referenced by SetQAAll().

{ fQAKs0 = qa;}

void PndSoftTriggerTask::SetQAMc ( bool  qa = true )

inline

Definition at line 106 of file softrig/PndSoftTriggerTask.h.

References fQAMc.

Referenced by SetQAAll().

{ fQAMc  = qa;}

void PndSoftTriggerTask::SetQAMctOnly ( bool  qa = true )

inline

Definition at line 107 of file softrig/PndSoftTriggerTask.h.

References fQAMctOnly.

Referenced by prod_fsim().

{ fQAMctOnly  = qa;}

void PndSoftTriggerTask::SetQAMode	(	int	mode,
		bool	qa = true
	)

Definition at line 411 of file softrig/PndSoftTriggerTask.cxx.

References fSTTriggers.

Referenced by ana_jpsi_task_ST(), prod_fsim(), and SetQAAll().

{
        int divi = 1;
        if (mode<10) divi=100;
        else if (mode<100) divi=10;
        
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
                if (it->first/divi == mode) it->second->SetWriteQA(qa);
                
//      if (fSTTriggers.find(mode) == fSTTriggers.end()) return;        
//      fSTTriggers[mode]->SetWriteQA(qa);
}

void PndSoftTriggerTask::SetQAPi0 ( bool  qa = true )

inline

Definition at line 103 of file softrig/PndSoftTriggerTask.h.

References fQAPi0.

Referenced by SetQAAll().

{ fQAPi0 = qa;}

void PndSoftTriggerTask::SetQASelectionDefaults

(

)

Definition at line 363 of file softrig/PndSoftTriggerTask.cxx.

References DbMass(), fEtaQaMax, fEtaQaMin, fKs0QaMax, fKs0QaMin, fPi0QaMax, and fPi0QaMin.

Referenced by PndSoftTriggerTask().

{
        double Pi0Mass = DbMass("pi0");
        double EtaMass = DbMass("eta");
        double Ks0Mass = DbMass("K_S0");
                
        // default windows are roughly +-15*sigma, and 15*sigma_max, if different channels are reconstructed
        // should not depend explicitly from set mean and sigma values, therefore fixed values are chosen
        
        fPi0QaMin = Pi0Mass - 0.06;             
        fPi0QaMax = Pi0Mass + 0.06;             
        
        fEtaQaMin = EtaMass - 0.10;             
        fEtaQaMax = EtaMass + 0.10;             
        
        fKs0QaMin = Ks0Mass - 0.10;             
        fKs0QaMax = Ks0Mass + 0.10;             
        
}

void PndSoftTriggerTask::SetSignalParamsDefaults

(

)

Definition at line 384 of file softrig/PndSoftTriggerTask.cxx.

References fEtaMean, fEtaSigma, fKs0Mean, fKs0Sigma, fPi0Mean, and fPi0Sigma.

Referenced by PndSoftTriggerTask().

{
        // *** fitted peak values from previous studies
                
        fPi0Mean = 0.136;               // mean value for pi0 signal
        fPi0Sigma = 0.0045;             // sigma value for pi0 signal
        
        fEtaMean = 0.552;               // mean value for eta(gg) signal
        fEtaSigma = 0.009;              // sigma value for eta(gg) signal
        
        fKs0Mean = 0.497;               // mean value for Ks signal
        fKs0Sigma = 0.008;              // sigma value for Ks signal
        
}       

void PndSoftTriggerTask::SetTagAll

(

bool

tag = true

)

Definition at line 453 of file softrig/PndSoftTriggerTask.cxx.

References fSTTriggers.

Referenced by ana_jpsi_task_ST(), SetFastSimDefaults(), and SetFullSimDefaults().

{       
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
                it->second->SetTagActive(tag);
}

void PndSoftTriggerTask::SetTagMode	(	int	mode,
		bool	tag = true
	)

Definition at line 439 of file softrig/PndSoftTriggerTask.cxx.

References fSTTriggers.

{
        int divi = 1;
        if (mode<10) divi=100;
        else if (mode<100) divi=10;
        
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
                if (it->first/divi == mode) it->second->SetTagActive(tag);
                
//      if (fSTTriggers.find(mode) == fSTTriggers.end()) return;
//      fSTTriggers[mode]->SetTagActive(tag);
}

void PndSoftTriggerTask::SetTagNSigAll

(

double

nsig

)

Definition at line 470 of file softrig/PndSoftTriggerTask.cxx.

References fSTTriggers.

{
        for (TrigIt it=fSTTriggers.begin(); it!=fSTTriggers.end(); ++it)
                it->second->SetTagNSig(nsig);   
}

void PndSoftTriggerTask::SetTagNSigMode	(	int	mode,
		double	nsig
	)

Definition at line 461 of file softrig/PndSoftTriggerTask.cxx.

References fSTTriggers, and mode.

{
        if (fSTTriggers.find(mode) == fSTTriggers.end()) return;
        
        fSTTriggers[mode]->SetTagNSig(nsig);
}

void PndSoftTriggerTask::SetTrackMinP ( double  min )

inline

Definition at line 85 of file softrig/PndSoftTriggerTask.h.

References fTrackMinP, and min().

Referenced by ana_jpsi_task_ST(), qa_softtrig(), SetFastSimDefaults(), and SetFullSimDefaults().

{fTrackMinP = min;}

void PndSoftTriggerTask::SetVerbose ( int  verb = 10 )

inline

Definition at line 118 of file softrig/PndSoftTriggerTask.h.

References fVerbose.

{fVerbose=verb;}

int PndSoftTriggerTask::SplitString ( TString  s, TString  delim, TString *  toks, int  maxtoks  )

private

Definition at line 789 of file softrig/PndSoftTriggerTask.cxx.

References i.

Referenced by ReadConfiguration(), and ReadTriggerLines().

{
        TObjArray *tok = s.Tokenize(delim);
        int N = tok->GetEntries();      
        for (int i=0;i<N;++i) 
                if (i<maxtoks) 
                {
                        toks[i] = ((TObjString*)tok->At(i))->String();
                        toks[i].ReplaceAll("\t","");
                        toks[i] = toks[i].Strip(TString::kBoth);
                }
        return N;
}

int PndSoftTriggerTask::TagMode ( PndSoftTriggerLine *  tl, int &  npre  )

private

Definition at line 1670 of file softrig/PndSoftTriggerTask.cxx.

References RhoMassParticleSelector::Accept(), AcceptCandidate(), AcceptDstCut(), RhoTuple::Column(), DoCombinatorics(), RhoTuple::DumpData(), f, fabs(), fAnalysis, fApplyFullSelection, fEventShape, fEvtCount, fIniP4, fMode, fPidMult_025, fPrimVtx, fPrimVtxQa, fQA, fQAMctOnly, fRecoilCnt, fRecoilMode, fRunNum, fSTencode, fSTModeIndex, fSTSelmap, RhoCandList::GetLength(), PndSoftTriggerLine::GetMean(), PndSoftTriggerLine::GetModeCode(), PndSoftTriggerLine::GetPrefix(), PndSoftTriggerLine::GetQASelector(), PndSoftTriggerLine::GetRhoTuple(), PndSoftTriggerLine::GetSelector(), PndSoftTriggerLine::GetSigma(), PndSoftTriggerLine::GetTagActive(), i, Mass, PndAnalysis::McTruthMatch(), mean, mode, n, nsig, PndRhoTupleQA::qaComp(), PndRhoTupleQA::qaEventShapeShort(), PndRhoTupleQA::qaP4(), RhoCandList::Select(), sigma, tag, and TString.

Referenced by Exec().

{
        // *** counter for full selection accepted cands
        int nacc = 0; 
        npre = 0;
        
        if ( !tl->GetTagActive() ) return 0;
        if ( fEcm < tl->GetThreshold() ) return 0;
        
        int mode       = tl->GetModeCode();
        RhoTuple *n    = tl->GetRhoTuple();
        TString prefix = tl->GetPrefix();
        double mean    = tl->GetMean();
        double sigma   = tl->GetSigma();
        
        RhoMassParticleSelector *sel = tl->GetSelector();
        
        RhoCandList l;

        // ** combinatorics
        DoCombinatorics(l, tl);
        
        // ** preselection
        l.Select(tl->GetQASelector());

        // *** store QA
        for (int i=0;i<l.GetLength();++i)
        {
                bool acc = false;
                bool tag = sel->Accept(l[i]) && AcceptDstCut(l[i]);
                
                // full selection?
                if (fApplyFullSelection>0) 
                        acc = AcceptCandidate(mode, l[i], sel);
                // if not, simple mass window selection (and D* cut if applicable)
                else acc = tag;
                
                // for full selection in open mode (=2), trigger w/o detailed cuts are accepted based on mass window only 
                if ( !acc && fApplyFullSelection==2 && fSTSelmap.find(fSTencode[fSTModeIndex]*1000+mode) == fSTSelmap.end() )
                        acc = tag;
                                
                if (acc) nacc++;                
                if (tag) npre++;
                
                fAnalysis->McTruthMatch(l[i]);
                
                if (n && AcceptDstCut(l[i]))
                {
                        // in case we only want to store MC truth matched candidates, check!
                        if (fQAMctOnly && (l[i]->GetMcTruth())==0) continue;
                        
                        fQA->qaComp(prefix, l[i], n);
                        fQA->qaEventShapeShort("es", fEventShape, n);
                        
                        // replace PID mult values from event shape by actual counts with individual algos
                        n->Column("eslnpide", (Int_t)  fPidMult_025[0], 0);
                        n->Column("eslnpidmu",(Int_t)  fPidMult_025[1], 0);
                        n->Column("eslnpidpi",(Int_t)  fPidMult_025[2], 0);
                        n->Column("eslnpidk", (Int_t)  fPidMult_025[3], 0);
                        n->Column("eslnpidp", (Int_t)  fPidMult_025[4], 0);
                        
                        fQA->qaP4("beam", fIniP4, n);
                        n->Column("primvx",   (Float_t)   fPrimVtx.X(), 0.0f);
                        n->Column("primvy",   (Float_t)   fPrimVtx.Y(), 0.0f);
                        n->Column("primvz",   (Float_t)   fPrimVtx.Z(), 0.0f);
                        n->Column("primvqa",  (Float_t)   fPrimVtxQa  , 0.0f);

                        Float_t mmiss = (fIniP4-(l[i]->P4())).M();
                        Float_t nsig = (Float_t) fabs(l[i]->Mass()-mean)/sigma;
                        
                        n->Column("ev",         (Int_t)   fEvtCount,    0);
                        n->Column("num",        (Int_t)   i,                    0);
                        n->Column("run",        (Int_t)   fRunNum,              0);
                        n->Column("mode",       (Int_t)   fMode,                0);
                        n->Column("recmode",(Int_t)   fRecoilMode,  0);
                        n->Column("reccnt", (Int_t)   fRecoilCnt,   0);
                        n->Column("mmiss",      (Float_t) mmiss,                0.0f);
                        n->Column("nsig",       (Float_t) nsig,                 0.0f);
                        n->Column(prefix+"mean",(Float_t) mean,         0.0f);
                        n->Column(prefix+"sig", (Float_t) sigma,        0.0f);
                        n->Column("tag",        (Int_t)   tag,                  0);
                        n->Column("acc",    (Int_t) acc,            0);
                        
                        n->DumpData();
                }
        }
        
        return nacc;
}

Member Data Documentation

int PndSoftTriggerTask::allmax

protected

Definition at line 48 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::chmax

protected

Definition at line 48 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::cntd0

protected

Definition at line 54 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::cntdpm

protected

Definition at line 55 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::cntds

protected

Definition at line 56 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::cntjpsi

protected

Definition at line 53 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::cntlamc

protected

Definition at line 58 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::cntphi

protected

Definition at line 57 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::cntsel

protected

Definition at line 51 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::evcount

protected

Definition at line 47 of file tutorials/anatask/PndSoftTriggerTask.h.

TString PndSoftTriggerTask::fAlgoElectron

private

Definition at line 242 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), PndSoftTriggerTask(), and SetPidAlgoElectron().

TString PndSoftTriggerTask::fAlgoKaon

private

Definition at line 245 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), PndSoftTriggerTask(), and SetPidAlgoKaon().

TString PndSoftTriggerTask::fAlgoMuon

private

Definition at line 243 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), PndSoftTriggerTask(), and SetPidAlgoMuon().

TString PndSoftTriggerTask::fAlgoPion

private

Definition at line 244 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), PndSoftTriggerTask(), and SetPidAlgoPion().

TString PndSoftTriggerTask::fAlgoProton

private

Definition at line 246 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), PndSoftTriggerTask(), and SetPidAlgoProton().

RhoCandList PndSoftTriggerTask::fAllCands

private

Definition at line 276 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), and FillGlobalLists().

PndAnalysis* PndSoftTriggerTask::fAnalysis

private

Definition at line 260 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), FillGlobalLists(), GetPocaVtx(), Init(), SelectTruePid(), and TagMode().

int PndSoftTriggerTask::fApplyFullSelection

private

Definition at line 176 of file softrig/PndSoftTriggerTask.h.

Referenced by ApplyFullSelection(), Exec(), Init(), PndSoftTriggerTask(), ReadConfiguration(), ReadTriggerLines(), and TagMode().

TString PndSoftTriggerTask::fCfgFileName

private

Definition at line 171 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), PndSoftTriggerTask(), ReadConfiguration(), and SetConfigurationFile().

RhoCandList PndSoftTriggerTask::fChargedCands

private

Definition at line 277 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), and FillGlobalLists().

double PndSoftTriggerTask::fDstMDiffCut

private

Definition at line 169 of file softrig/PndSoftTriggerTask.h.

Referenced by AcceptDstCut(), and SetDstMDiffCut().

double PndSoftTriggerTask::fEcm

private

Definition at line 182 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), Init(), PndSoftTriggerTask(), and ReadTriggerLines().

RhoEnergyParticleSelector* PndSoftTriggerTask::fEnergySel

private

Definition at line 239 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), and Init().

RhoCandList PndSoftTriggerTask::fEtaCands

private

Definition at line 283 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists().

double PndSoftTriggerTask::fEtaMean

private

Definition at line 213 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), Init(), SetEtaSignalParams(), and SetSignalParamsDefaults().

double PndSoftTriggerTask::fEtaMinE

private

Definition at line 222 of file softrig/PndSoftTriggerTask.h.

RhoMassParticleSelector* PndSoftTriggerTask::fEtaPreSel

private

Definition at line 231 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), and Init().

double PndSoftTriggerTask::fEtaQaMax

private

Definition at line 202 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), SetEtaQASelection(), and SetQASelectionDefaults().

double PndSoftTriggerTask::fEtaQaMin

private

Definition at line 201 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), SetEtaQASelection(), and SetQASelectionDefaults().

RhoMassParticleSelector* PndSoftTriggerTask::fEtaSel

private

Definition at line 230 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), and Init().

double PndSoftTriggerTask::fEtaSigma

private

Definition at line 214 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), Init(), SetEtaSignalParams(), and SetSignalParamsDefaults().

PndEventShape* PndSoftTriggerTask::fEventShape

private

Definition at line 263 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), FillEventShapeVarArray(), and TagMode().

int PndSoftTriggerTask::fEvtCount

private

Definition at line 164 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), Exec(), FillGlobalLists(), and TagMode().

RhoCandList PndSoftTriggerTask::fGammaCands

private

Definition at line 279 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists().

double PndSoftTriggerTask::fGammaMinE

private

Definition at line 220 of file softrig/PndSoftTriggerTask.h.

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

TLorentzVector PndSoftTriggerTask::fIniP4

private

Definition at line 181 of file softrig/PndSoftTriggerTask.h.

Referenced by BoostCms(), CreateKs0Cands(), Exec(), FillGlobalLists(), GetVarValue(), PndSoftTriggerTask(), and TagMode().

double PndSoftTriggerTask::fIniPidCut

private

Definition at line 224 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), PndSoftTriggerTask(), and SetInitialPidCut().

RhoCandList PndSoftTriggerTask::fKs0Cands

private

Definition at line 281 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), and FillGlobalLists().

double PndSoftTriggerTask::fKs0Mean

private

Definition at line 216 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), Init(), SetKs0SignalParams(), and SetSignalParamsDefaults().

RhoMassParticleSelector* PndSoftTriggerTask::fKs0PreSel

private

Definition at line 234 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), and Init().

double PndSoftTriggerTask::fKs0QaMax

private

Definition at line 205 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), SetKs0QASelection(), and SetQASelectionDefaults().

double PndSoftTriggerTask::fKs0QaMin

private

Definition at line 204 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), SetKs0QASelection(), and SetQASelectionDefaults().

RhoMassParticleSelector* PndSoftTriggerTask::fKs0Sel

private

Definition at line 233 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), and Init().

double PndSoftTriggerTask::fKs0Sigma

private

Definition at line 217 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), Init(), SetKs0SignalParams(), and SetSignalParamsDefaults().

RhoCandList PndSoftTriggerTask::fMcTruth

private

Definition at line 275 of file softrig/PndSoftTriggerTask.h.

Referenced by DetermineRecoilMode(), Exec(), and FillGlobalLists().

int PndSoftTriggerTask::fMode

private

Definition at line 161 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), Exec(), FillGlobalLists(), and TagMode().

RhoMomentumParticleSelector* PndSoftTriggerTask::fMomentumSel

private

Definition at line 238 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), and Init().

RhoCandList PndSoftTriggerTask::fNeutralCands

private

Definition at line 278 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists().

double PndSoftTriggerTask::fNsigAux

private

Definition at line 168 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), FillGlobalLists(), Init(), and SetAuxNSig().

double PndSoftTriggerTask::fNsigTag

private

Definition at line 167 of file softrig/PndSoftTriggerTask.h.

double PndSoftTriggerTask::fPbarMom

private

Definition at line 183 of file softrig/PndSoftTriggerTask.h.

Referenced by Init().

TDatabasePDG* PndSoftTriggerTask::fPdg

private

Definition at line 272 of file softrig/PndSoftTriggerTask.h.

Referenced by AntiPdg(), DbMass(), and PndSoftTriggerTask().

Int_t PndSoftTriggerTask::fPhotosMax

private

Definition at line 173 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), and McMatchAllowPhotos().

Double_t PndSoftTriggerTask::fPhotosThresh

private

Definition at line 174 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), and McMatchAllowPhotos().

RhoCandList PndSoftTriggerTask::fPi0Cands

private

Definition at line 282 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists().

double PndSoftTriggerTask::fPi0Mean

private

Definition at line 210 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), Init(), SetPi0SignalParams(), and SetSignalParamsDefaults().

double PndSoftTriggerTask::fPi0MinE

private

Definition at line 221 of file softrig/PndSoftTriggerTask.h.

RhoMassParticleSelector* PndSoftTriggerTask::fPi0PreSel

private

Definition at line 228 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), and Init().

double PndSoftTriggerTask::fPi0QaMax

private

Definition at line 199 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), SetPi0QASelection(), and SetQASelectionDefaults().

double PndSoftTriggerTask::fPi0QaMin

private

Definition at line 198 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), SetPi0QASelection(), and SetQASelectionDefaults().

RhoMassParticleSelector* PndSoftTriggerTask::fPi0Sel

private

Definition at line 227 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), and Init().

double PndSoftTriggerTask::fPi0Sigma

private

Definition at line 211 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), Init(), SetPi0SignalParams(), and SetSignalParamsDefaults().

RhoCandList PndSoftTriggerTask::fPidList[16]

private

Definition at line 286 of file softrig/PndSoftTriggerTask.h.

Referenced by CombineList(), CreateKs0Cands(), DoCombinatorics(), and FillGlobalLists().

Int_t PndSoftTriggerTask::fPidMult_025[5]

private

Definition at line 249 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), FillEventShapeVarArray(), FillGlobalLists(), and TagMode().

RhoVtxPoca* PndSoftTriggerTask::fPocaVertexer

private

Definition at line 269 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), GetPocaVtx(), Init(), and ~PndSoftTriggerTask().

TVector3 PndSoftTriggerTask::fPrimVtx

private

Definition at line 184 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), GetPocaVtx(), and TagMode().

double PndSoftTriggerTask::fPrimVtxQa

private

Definition at line 185 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), and TagMode().

PndRhoTupleQA* PndSoftTriggerTask::fQA

private

Definition at line 266 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), Exec(), FillGlobalLists(), Init(), and TagMode().

bool PndSoftTriggerTask::fQAEta

private

Definition at line 190 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), Init(), and SetQAEta().

bool PndSoftTriggerTask::fQAEvent

private

Definition at line 192 of file softrig/PndSoftTriggerTask.h.

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

bool PndSoftTriggerTask::fQAKs0

private

Definition at line 191 of file softrig/PndSoftTriggerTask.h.

Referenced by Init(), and SetQAKs0().

bool PndSoftTriggerTask::fQAMc

private

Definition at line 193 of file softrig/PndSoftTriggerTask.h.

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

bool PndSoftTriggerTask::fQAMctOnly

private

Definition at line 194 of file softrig/PndSoftTriggerTask.h.

Referenced by SetQAMctOnly(), and TagMode().

bool PndSoftTriggerTask::fQAPi0

private

Definition at line 189 of file softrig/PndSoftTriggerTask.h.

Referenced by FillGlobalLists(), Init(), and SetQAPi0().

int PndSoftTriggerTask::fRecoilCnt

private

Definition at line 163 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), and TagMode().

int PndSoftTriggerTask::fRecoilMode

private

Definition at line 162 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), and TagMode().

FairRootManager* PndSoftTriggerTask::fRootManager

private

Definition at line 292 of file softrig/PndSoftTriggerTask.h.

Referenced by Init().

int PndSoftTriggerTask::fRunNum

private

Definition at line 165 of file softrig/PndSoftTriggerTask.h.

Referenced by CreateKs0Cands(), Exec(), FillGlobalLists(), and TagMode().

int PndSoftTriggerTask::fSigCount

private

Definition at line 166 of file softrig/PndSoftTriggerTask.h.

TClonesArray* PndSoftTriggerTask::fTcaOnlineFilterInfo

private

Definition at line 293 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), and Init().

double PndSoftTriggerTask::fTrackMinP

private

Definition at line 223 of file softrig/PndSoftTriggerTask.h.

Referenced by Exec(), FillGlobalLists(), Init(), and SetTrackMinP().

TString PndSoftTriggerTask::fTriggerFileName

private

Definition at line 172 of file softrig/PndSoftTriggerTask.h.

Referenced by PndSoftTriggerTask(), and ReadTriggerLines().

int PndSoftTriggerTask::fVerbose

private

Definition at line 160 of file softrig/PndSoftTriggerTask.h.

Referenced by AcceptCandidate(), DoCombinatorics(), Init(), ReadConfiguration(), and SetVerbose().

TH1F* PndSoftTriggerTask::h_d0

private

Definition at line 88 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_l

private

Definition at line 100 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_lm

private

Definition at line 101 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_pocx

private

Definition at line 92 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_pocxm

private

Definition at line 96 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_pocy

private

Definition at line 93 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_pocym

private

Definition at line 97 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_pocz

private

Definition at line 94 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0_poczm

private

Definition at line 98 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0sel

private

Definition at line 90 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_d0t

private

Definition at line 89 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_dpm

private

Definition at line 104 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_dpmsel

private

Definition at line 106 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_dpmt

private

Definition at line 105 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_ds

private

Definition at line 109 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_dssel

private

Definition at line 111 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_dst

private

Definition at line 110 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_jpsi

private

Definition at line 83 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_jpsisel

private

Definition at line 85 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_jpsit

private

Definition at line 84 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_lamc

private

Definition at line 119 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_lamcsel

private

Definition at line 121 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_lamct

private

Definition at line 120 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_mom

private

book all the histograms

Definition at line 74 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_momc

private

Definition at line 75 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_momn

private

Definition at line 76 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_mult

private

Definition at line 77 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_multc

private

Definition at line 78 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_multn

private

Definition at line 79 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_multrem

private

Definition at line 80 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_phi

private

Definition at line 114 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_phisel

private

Definition at line 116 of file tutorials/anatask/PndSoftTriggerTask.h.

TH1F* PndSoftTriggerTask::h_phit

private

Definition at line 115 of file tutorials/anatask/PndSoftTriggerTask.h.

int PndSoftTriggerTask::mcmax

protected

Definition at line 48 of file tutorials/anatask/PndSoftTriggerTask.h.

RhoTuple* PndSoftTriggerTask::neta

private

Definition at line 256 of file softrig/PndSoftTriggerTask.h.

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

int PndSoftTriggerTask::neutmax

protected

Definition at line 48 of file tutorials/anatask/PndSoftTriggerTask.h.

RhoTuple* PndSoftTriggerTask::nks0

private

Definition at line 254 of file softrig/PndSoftTriggerTask.h.

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

RhoTuple* PndSoftTriggerTask::nmc

private

Definition at line 257 of file softrig/PndSoftTriggerTask.h.

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

RhoTuple* PndSoftTriggerTask::npi0

private

Definition at line 255 of file softrig/PndSoftTriggerTask.h.

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

RhoTuple * PndSoftTriggerTask::ntp

private

Definition at line 253 of file softrig/PndSoftTriggerTask.h.

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

RhoTuple* PndSoftTriggerTask::ntp2

private

Definition at line 125 of file tutorials/anatask/PndSoftTriggerTask.h.

PndAnalysis* PndSoftTriggerTask::theAnalysis

private

Definition at line 130 of file tutorials/anatask/PndSoftTriggerTask.h.

TStopwatch* PndSoftTriggerTask::timer

private

Definition at line 131 of file tutorials/anatask/PndSoftTriggerTask.h.

---

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

• softrig/PndSoftTriggerTask.h
• softrig/PndSoftTriggerTask.cxx