PndSimpleCombiner Class Reference

#include <PndSimpleCombiner.h>

Public Member Functions
	PndSimpleCombiner (PndAnalysis *fAna, TString decay, TString params="", double Ecm=0)
	~PndSimpleCombiner ()
void	SetVerbose (int verb=1)
void	Combine ()
void	Print ()
void	SetPid (TString crit="", TString algo="")
void	SetPidElectron (TString crit="", TString algo="")
void	SetPidMuon (TString crit="", TString algo="")
void	SetPidPion (TString crit="", TString algo="")
void	SetPidKaon (TString crit="", TString algo="")
void	SetPidProton (TString crit="", TString algo="")
int	GetNLists ()
bool	GetList (RhoCandList &l, TString comp)
bool	GetList (RhoCandList &l, int pdg)
bool	GetListN (RhoCandList &l, int idx)

Private Member Functions
void	InitDecayInfo (SCDecayInfo &info, int pdg, int idx)
void	FillGenericLists ()
int	CombineList (RhoCandList &l, int mpdg, std::vector< int > &idx)
int	SplitString (TString s, TString delim, StringList &toks)
bool	ParseDecay (TString decay)
bool	ParseParams (TString params)
bool	CCInvariant (std::vector< int > &vpdg)
bool	CCInvariant (int pdg)
int	AntiPdg (int pdg)

Private Attributes
PndAnalysis *	fAnalysis
TString	fDecay
TString	fGlobParams
int	fNLists
RhoCandList	fList [MAXLISTS]
int	fVerbose
double	fEmin
double	fPmin
double	fEcm
RhoEnergyParticleSelector *	fESel
RhoMomentumParticleSelector *	fPSel
std::map< int, int >	fPdgIdxMap
std::map< int, int >	fIdxPdgMap
std::map< int, TString >	fIdxListNameMap
std::map< int, TString >	fIdxPidCritMap
std::map< int, TString >	fIdxPidAlgoMap
std::vector< SCDecayInfo >	fDecayInfoArray

Detailed Description

Definition at line 59 of file PndSimpleCombiner.h.

Constructor & Destructor Documentation

PndSimpleCombiner::PndSimpleCombiner	(	PndAnalysis *	fAna,
		TString	decay,
		TString	params = "",
		double	Ecm = 0
	)

Definition at line 55 of file PndSimpleCombiner.cxx.

References fEcm, fIdxListNameMap, fIdxPdgMap, fPdgIdxMap, i, ParseDecay(), ParseParams(), SetPid(), and TString.

                                                                                                 : 
        fAnalysis(fAna), fDecay(decay), fGlobParams(params), fNLists(11), fVerbose(0), 
        fEmin(0.), fPmin(0.), fEcm(Ecm), fESel(0), fPSel(0)             
{       
        // initialize mapping pdg -> list index and list name
        int pdgcodes[]     = {    -11,              11,           -13,         13,         211,        -211,        321,        -321,        2212,         -2212,         22};
        TString pdgnames[] = {"ElectronPlus", "ElectronMinus", "MuonPlus", "MuonMinus", "PionPlus", "PionMinus", "KaonPlus", "KaonMinus", "ProtonPlus", "ProtonMinus", "Neutral"};
        
        fPdgIdxMap.clear();
        fIdxPdgMap.clear();
        fIdxListNameMap.clear();
        
        for (int i=0;i<11;++i) 
        {
                fPdgIdxMap[pdgcodes[i]]  = i;             // maps pdg code -> list index
                fIdxPdgMap[i]            = pdgcodes[i];   // maps list index -> pdg code
                fIdxListNameMap[i]       = pdgnames[i];   // maps list index -> generic list name (ElectronPlus, PionMinus, ...; see above)
        }
        // set initial pid configuration
        SetPid("All", "PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc;PidAlgoStt;PidAlgoMdtHardCuts;PidAlgoSciT;PidAlgoRich;PidAlgoFtof");
        
        if (fEcm==0) fEcm=-999;
        
        assert(ParseDecay(decay));
        ParseParams(params);
}

PndSimpleCombiner::~PndSimpleCombiner

(

)

Definition at line 84 of file PndSimpleCombiner.cxx.

References fDecayInfoArray, fESel, fPSel, and i.

{
        if (fESel) delete fESel;
        if (fPSel) delete fPSel;

        for (size_t i=0;i<fDecayInfoArray.size();++i) if (fDecayInfoArray[i].msel) delete fDecayInfoArray[i].msel;
}

Member Function Documentation

int PndSimpleCombiner::AntiPdg ( int  pdg )

private

Definition at line 475 of file PndSimpleCombiner.cxx.

Referenced by CCInvariant(), Combine(), FillGenericLists(), and ParseDecay().

{
        if (TDatabasePDG::Instance()->GetParticle(pdg))
        {
                if (!TDatabasePDG::Instance()->GetParticle(pdg)->AntiParticle()) return pdg;
                else return TDatabasePDG::Instance()->GetParticle(pdg)->AntiParticle()->PdgCode();
        }
        return -999999;
} 

bool PndSimpleCombiner::CCInvariant ( std::vector< int > &  vpdg )

private

Definition at line 457 of file PndSimpleCombiner.cxx.

References i.

Referenced by ParseDecay().

{
        int sum=0, nminus=0, nplus=0;
        
        for (size_t i=0;i<vpdg.size();++i)
        {
                if (!CCInvariant(vpdg[i])) 
                {
                        sum += vpdg[i];
                        if (vpdg[i]>0) nplus++;
                        else nminus++;
                }
        }
        return (sum==0 && nplus==nminus);
}

bool PndSimpleCombiner::CCInvariant ( int  pdg )

inlineprivate

Definition at line 92 of file PndSimpleCombiner.h.

References AntiPdg().

{return (pdg == AntiPdg(pdg));}   // is pdg = pdg of anti-particle

void PndSimpleCombiner::Combine

(

)

Definition at line 524 of file PndSimpleCombiner.cxx.

References AntiPdg(), RhoCandList::Append(), CombineList(), SCDecayInfo::daucc, SCDecayInfo::didx, SCDecayInfo::dpdg, fDecayInfoArray, FillGenericLists(), fList, fPdgIdxMap, i, SCDecayInfo::midx, SCDecayInfo::mpdg, SCDecayInfo::msel, n, and RhoCandList::Select().

Referenced by PndSimpleCombinerTask::Exec().

{
        // fill list of final state particles
        FillGenericLists();
        
        // loop through list definitions and produce composites
        int n = fDecayInfoArray.size();
        for (int i=0;i<n;++i)
        {
                SCDecayInfo info = fDecayInfoArray[i];
                CombineList(fList[info.midx], info.mpdg, info.didx);
                
                // do we need to add the cc FS?
                if (info.daucc)
                {
                        // create index list with cc daughters
                        std::vector<int> aidx;
                        for (size_t j=0;j<info.didx.size();++j)
                        {
                                int apdg = AntiPdg(info.dpdg[j]);
                                aidx.push_back(fPdgIdxMap[apdg]);                       
                        }
                        // create temporary list with the cc combinatorics
                        RhoCandList l;
                        CombineList(l, info.mpdg, aidx);
                        // append it to the original list
                        fList[info.midx].Append(l);
                }
                // if there is a mass selector connected, apply it
                if (info.msel) fList[info.midx].Select(info.msel);
        }
}

int PndSimpleCombiner::CombineList ( RhoCandList &  l, int  mpdg, std::vector< int > &  idx  )

private

Definition at line 559 of file PndSimpleCombiner.cxx.

References RhoCandList::Cleanup(), RhoCandList::Combine(), fList, and RhoCandList::GetLength().

Referenced by Combine().

{
        l.Cleanup();
        
        int nd = idx.size();
        
        switch (nd) 
        {
        case 2: l.Combine(fList[idx[0]], fList[idx[1]], mpdg); break;   
        case 3: l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], mpdg); break;    
        case 4: l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], fList[idx[3]], mpdg); break;
        case 5: l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], fList[idx[3]], fList[idx[4]], mpdg); break;
        
        case 6: l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], fList[idx[3]], fList[idx[4]], fList[idx[5]], mpdg); break;
        case 7: l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], fList[idx[3]], fList[idx[4]], fList[idx[5]], fList[idx[6]], mpdg); break;
        case 8: l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], fList[idx[3]], fList[idx[4]], fList[idx[5]], fList[idx[6]], fList[idx[7]], mpdg); break;
        case 9: l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], fList[idx[3]], fList[idx[4]], fList[idx[5]], fList[idx[6]], fList[idx[7]], fList[idx[8]], mpdg); break;
        case 10:l.Combine(fList[idx[0]], fList[idx[1]], fList[idx[2]], fList[idx[3]], fList[idx[4]], fList[idx[5]], fList[idx[6]], fList[idx[7]], fList[idx[8]], fList[idx[9]], mpdg); break;
        }
        
        return l.GetLength();
}

void PndSimpleCombiner::FillGenericLists ( )

private

Definition at line 409 of file PndSimpleCombiner.cxx.

References AntiPdg(), SCDecayInfo::didx, SCDecayInfo::dpdg, fAnalysis, fDecayInfoArray, fESel, fIdxListNameMap, fIdxPdgMap, fIdxPidAlgoMap, fIdxPidCritMap, PndAnalysis::FillList(), fList, fPdgIdxMap, fPSel, fVerbose, RhoCandList::GetLength(), i, n, SCDecayInfo::ndaug, and RhoCandList::Select().

Referenced by Combine().

{
        // fill all generic lists, which are used by any of the composites
        // loop through composites
        bool filled[11] = {0};  // marker, whether list is already filled
        
        int n = fDecayInfoArray.size();
        for (int i=0;i<n;++i)
        {
                SCDecayInfo info = fDecayInfoArray[i];
                
                for (int j=0;j<info.ndaug;++j)
                {
                        int lidx  = info.didx[j];
                        
                        if (lidx>10 || filled[lidx]) continue; // this list is not generic or already filled
                                
                        filled[lidx] = true;                   // mark list as filled
                        int laidx = fPdgIdxMap[AntiPdg(info.dpdg[j])];  // index of antiparticle list (only relevant for charged)

                        if (lidx<10) // charged lists
                        {
                                fAnalysis->FillList(fList[lidx],fIdxListNameMap[lidx]+fIdxPidCritMap[lidx],fIdxPidAlgoMap[lidx]);
                                if (fPSel) fList[lidx].Select(fPSel);
                                
                                // also fill anti pdg generic list for charged generic ...
                                fAnalysis->FillList(fList[laidx],fIdxListNameMap[laidx]+fIdxPidCritMap[laidx],fIdxPidAlgoMap[laidx]);
                                if (fPSel) fList[laidx].Select(fPSel);
                                filled[laidx] = true;              // ... and mark list as filled              
                        }
                        else if (lidx==10) //neutrals
                        {
                                fAnalysis->FillList(fList[lidx],fIdxListNameMap[lidx]);
                                if (fESel) fList[lidx].Select(fESel);
                        }
                        if (lidx<11 && fVerbose) 
                        {
                                cout <<"idx :"<<lidx<<" pdg:"<<fIdxPdgMap[lidx]<<" N:"<<fList[lidx].GetLength()<<endl;
                                if (lidx<10) cout <<"aidx:"<<laidx<<" pdg:"<<fIdxPdgMap[laidx]<<" N:"<<fList[laidx].GetLength()<<endl;
                        }
                }
        }
}

bool PndSimpleCombiner::GetList	(	RhoCandList &	l,
		TString	comp
	)

Definition at line 585 of file PndSimpleCombiner.cxx.

Referenced by PndSimpleCombinerTask::Exec().

{
        if (!TDatabasePDG::Instance()->GetParticle(comp)) return false;
        
        return GetList(l, TDatabasePDG::Instance()->GetParticle(comp)->PdgCode());
}

bool PndSimpleCombiner::GetList	(	RhoCandList &	l,
		int	pdg
	)

Definition at line 594 of file PndSimpleCombiner.cxx.

References RhoCandList::Cleanup(), fList, and fPdgIdxMap.

{
        l.Cleanup();
        
        if (fPdgIdxMap.find(pdg) == fPdgIdxMap.end()) return false;
        
        l = fList[fPdgIdxMap[pdg]];
        return true;
}

bool PndSimpleCombiner::GetListN	(	RhoCandList &	l,
		int	idx
	)

Definition at line 606 of file PndSimpleCombiner.cxx.

References RhoCandList::Cleanup(), fList, and GetNLists().

{
        l.Cleanup();
        
        if (idx>=0 && idx<GetNLists()) 
        {
                l = fList[idx+11];
                return true;
        }
        
        return false;
}

int PndSimpleCombiner::GetNLists ( )

inline

Definition at line 79 of file PndSimpleCombiner.h.

References fNLists.

Referenced by GetListN().

{return fNLists-11;}

void PndSimpleCombiner::InitDecayInfo ( SCDecayInfo &  info, int  pdg, int  idx  )

private

Definition at line 114 of file PndSimpleCombiner.cxx.

References SCDecayInfo::daucc, SCDecayInfo::didx, SCDecayInfo::dpdg, idx, SCDecayInfo::midx, SCDecayInfo::mpdg, SCDecayInfo::msel, SCDecayInfo::mwin, SCDecayInfo::mwinhi, and SCDecayInfo::mwinlo.

Referenced by ParseDecay().

{
        info.mpdg   = pdg;     // mother(composite) pdg code
        info.midx   = idx;     // add another list for the new composite
        info.daucc  = false;   // in case only the daughters have a cc, but not the mother, we have to care (e.g. etac -> Ks K+ pi-)
        info.mwin   = 0;       // default: no mass selection
        info.mwinlo = 0;       // default: no mass selection
        info.mwinhi = 0;       // default: no mass selection
        info.msel   = 0;       // default: no mass selection
        info.dpdg.clear();     // pdgs of daughters
        info.didx.clear();     // list index of daughters
}

bool PndSimpleCombiner::ParseDecay ( TString  decay )

private

Definition at line 129 of file PndSimpleCombiner.cxx.

References AntiPdg(), CCInvariant(), SCDecayInfo::didx, SCDecayInfo::dpdg, fDecayInfoArray, fIdxPdgMap, fNLists, fPdgIdxMap, i, InitDecayInfo(), SCDecayInfo::midx, SCDecayInfo::mpdg, SCDecayInfo::ndaug, SplitString(), and TString.

Referenced by PndSimpleCombiner().

{
        StringList subdec;
        
        int ndec = SplitString(decay, ";", subdec);
        
        // loop over subdecays
        for (int i=0;i<ndec;++i)
        {       
                StringList dectoks;
                // does decay string contain exactly one '->' in the middle of the string?
                subdec[i].ReplaceAll("->",">");
                SplitString(subdec[i],">",dectoks);
                                
                if ( dectoks.size()!=2 || !TDatabasePDG::Instance()->GetParticle(dectoks[0]) ) {cout <<"[PndSimpleCombiner] **** ERROR : Invalid decay pattern '"<<subdec[i].Data()<<"'"<<endl; return false;}
                
                // split string again; first token is supposed to be the decaying resonance
                TString curstr = dectoks[0]+" "+dectoks[1];
                SplitString(curstr," ",dectoks);
                
                // too many daughters (>10)
                if (dectoks.size()>11) {cout <<"[PndSimpleCombiner] **** ERROR : Exceeding max number of daughters (10): '"<<subdec[i].Data()<<"'"<<endl; return false;}
                
                SCDecayInfo info;
                InitDecayInfo(info, TDatabasePDG::Instance()->GetParticle(dectoks[0])->PdgCode(), fNLists++);
                
                fPdgIdxMap[info.mpdg]  = info.midx;                     // add new list to the pdg <-> list index maps
                fIdxPdgMap[info.midx]  = info.mpdg;
                
                bool cc=true;
                
                // loop over daughters
                for (size_t j=1;j<dectoks.size();++j)
                {
                        if (dectoks[j]=="nocc") { cc = false; continue; }
                        
                        if (!TDatabasePDG::Instance()->GetParticle(dectoks[j])) {cout <<"[PndSimpleCombiner] **** ERROR : Unknown particle '"<<dectoks[j].Data()<<"'"<<endl; return false;}
                        
                        // pdg code of daughter
                        int dpdg = TDatabasePDG::Instance()->GetParticle(dectoks[j])->PdgCode();
                        // add corresponding index if exists
                        if ( fPdgIdxMap.find(dpdg) == fPdgIdxMap.end() ) {cout <<"[PndSimpleCombiner] **** ERROR : Undefined list '"<<dectoks[j].Data()<<"'"<<endl; return false;}
                        if ( dpdg == info.mpdg ) {cout <<"[PndSimpleCombiner] **** ERROR : Invalid recursion in '"<<subdec[i].Data()<<"'"<<endl; return false;}
                        
                        info.dpdg.push_back(dpdg);
                        info.didx.push_back(fPdgIdxMap[dpdg]);
                }
                info.ndaug = info.didx.size();
                fDecayInfoArray.push_back(info);
                
                // do we need the cc decay definition
                if (cc)
                {
                        // does the cc mode exist?
                        // either the mother or the FS has to have a cc not being identical
                        // example: 1) D0 -> K- pi+       has cc -> D0bar -> K+ pi-     (different FS and mother)
                        //          2) D0 -> KS pi+ pi-   has cc -> D0bar -> KS pi+ pi- (different mother)
                        //          3) etac -> KS K- pi+  has cc -> etac -> KS K+ pi-   (different FS)
                        //          4) phi -> K+ K-       has _no_ cc!                  (identical FS and mother)
                        
                        // does the mother have a cc? (case 1 + 2)
                        // then we add a new list
                        if (!CCInvariant(info.mpdg)) 
                        {
                                SCDecayInfo ainfo;
                                InitDecayInfo(ainfo, AntiPdg(info.mpdg), fNLists++);

                                fPdgIdxMap[ainfo.mpdg]  = ainfo.midx;
                                fIdxPdgMap[ainfo.midx]  = ainfo.mpdg;
                                
                                for (size_t j=0;j<info.dpdg.size();++j)
                                {
                                        int apdg = AntiPdg(info.dpdg[j]);
                                        if (apdg==-999999 || fPdgIdxMap.find(apdg) == fPdgIdxMap.end() ) {cout <<"[PndSimpleCombiner] **** ERROR : No list for PDG code "<<apdg<<endl; return false;}
                                        ainfo.dpdg.push_back(apdg);
                                        ainfo.didx.push_back(fPdgIdxMap[apdg]);
                                }
                                ainfo.ndaug = info.didx.size();
                                fDecayInfoArray.push_back(ainfo);
                        }
                        // only the daughters have a cc (case 3)
                        // then we simply set the daucc switch, which will be taken into account during combinatorics
                        else if (!CCInvariant(info.dpdg))
                        {
                                for (size_t j=0;j<info.dpdg.size();++j)
                                {
                                        int apdg = AntiPdg(info.dpdg[j]);
                                        if (apdg==-999999 || fPdgIdxMap.find(apdg) == fPdgIdxMap.end() ) {cout <<"[PndSimpleCombiner] **** ERROR : No list for PDG code "<<apdg<<endl; return false;}
                                }
                                
                                fDecayInfoArray[fDecayInfoArray.size()-1].daucc=true;
                        }
                }
        }
        
        return true;
}

bool PndSimpleCombiner::ParseParams ( TString  params )

private

Definition at line 229 of file PndSimpleCombiner.cxx.

References fDecayInfoArray, fEcm, fEmin, fESel, fIdxListNameMap, fPdgIdxMap, fPmin, fPSel, i, mean, SCDecayInfo::mpdg, SCDecayInfo::msel, SCDecayInfo::mwin, SCDecayInfo::mwinhi, SCDecayInfo::mwinlo, SetPid(), SetPidElectron(), SetPidKaon(), SetPidMuon(), SetPidPion(), SetPidProton(), SplitString(), and TString.

Referenced by PndSimpleCombiner().

{
        StringList parm;
        SplitString(params,":",parm);
        
        for (size_t i=0;i<parm.size();++i)
        {
                if (parm[i]=="ebrem") 
                {
                        fIdxListNameMap[0]="ElectronPlusBrem";
                        fIdxListNameMap[1]="ElectronMinusBrem";
                        continue;
                }
                
                StringList pair;
                SplitString(parm[i],"=",pair);
                
                if (pair.size()!=2) {cout <<"[PndSimpleCombiner] **** WARNING : Invalid parameter setting '"<<parm[i]<<"' ignored"<<endl; continue;}
                
                // ****
                // global mass window setting
                // can be either : mwin(part)=0.5        -> cut +- 0.25 around nominal mass of part
                // or            : mwin(part)=0.25|0.75  -> cut 0.25 < mpart < 0.75
                // ****
                if (pair[0]=="mwin")
                {
                        double window = pair[1].Atof();
                        
                        for (size_t j=0;j<fDecayInfoArray.size();++j)
                        {
                                SCDecayInfo &info = fDecayInfoArray[j];
                                if (info.msel) delete info.msel;
                                
                                // fetch center as nominal pdg mass; if pbarpSystemX (8888x) set center to Ecm
                                double mean = TDatabasePDG::Instance()->GetParticle(info.mpdg)->Mass();

                                if (info.mpdg/10 == 8888) { mean = fEcm;}
                                
                                info.mwin   = window;
                                info.mwinlo = mean - window/2; 
                                info.mwinhi = mean + window/2; 
                                
                                // check if Ecm was not set and particle is pbarpSystem
                                if (mean>-999)
                                        info.msel = new RhoMassParticleSelector("msel",mean,window);
                                else
                                {
                                        info.mwin = -1; // indicates that we skip the mass window cut for pbarpSystem if Ecm was not set
                                }
                        }
                }
                // ****
                // mass window for one composite
                // ****
                else if (pair[0].BeginsWith("mwin"))
                {
                        // extract particle name from strings like 'mwin(D0)'
                        pair[0] = pair[0](5,pair[0].Length()-6);
                        
                        if (!TDatabasePDG::Instance()->GetParticle(pair[0])) {cout <<"[PndSimpleCombiner] **** WARNING : Unknown particle type '"<<pair[0]<<"'"<<endl;continue;}
                        
                        int pdg = TDatabasePDG::Instance()->GetParticle(pair[0])->PdgCode();
                        if (fPdgIdxMap.find(pdg) == fPdgIdxMap.end()) {cout <<"[PndSimpleCombiner] **** WARNING : Unknown particle list '"<<pair[0]<<"'"<<endl;continue;}
                        
                        // set default mean and window
                        double mean = TDatabasePDG::Instance()->GetParticle(pdg)->Mass(), window = 0;
                        if (pdg/10 == 8888) { mean = fEcm; }
                        
                        // range setting like 1.2|1.6 ?
                        if (pair[1].Contains("|"))
                        {
                                double lower = TString((pair[1])(0,pair[1].Index("|"))).Atof();
                                double upper = TString((pair[1])(pair[1].Index("|")+1,1000)).Atof();
                                
                                mean   = (upper+lower)/2.;
                                window = upper-lower;
                                
                                if (mean<0 || window <=0) {cout <<"[PndSimpleCombiner] **** WARNING : Invalid mass window defintion '"<<parm[i]<<"'"<<endl;continue;}
                        }
                        else window = pair[1].Atof();
                        
                        for (size_t j=0;j<fDecayInfoArray.size();++j)
                        {
                                SCDecayInfo &info = fDecayInfoArray[j];
                                // only set for this particle type
                                if (abs(info.mpdg) == abs(pdg))
                                {
                                        if (info.msel) delete info.msel;
                                        info.mwin   = window;
                                        info.mwinlo = mean - window/2;
                                        info.mwinhi = mean + window/2;
                                        
                                        // check if Ecm was not set and particle is pbarpSystem
                                        if (mean>-999)
                                                info.msel = new RhoMassParticleSelector("msel",mean,window);
                                        else
                                        {
                                                info.mwin = -1; // indicates that we skip the mass window cut for pbarpSystem if Ecm was not set
                                        }
                                }
                        }       
                }
                
                // ****
                // check for pid setting
                // ****
                if (pair[0] == "pid")   SetPid(pair[1]);
                
                if (pair[0] == "pide")  SetPidElectron(pair[1]);
                if (pair[0] == "pidmu") SetPidMuon(pair[1]);
                if (pair[0] == "pidpi") SetPidPion(pair[1]);
                if (pair[0] == "pidk")  SetPidKaon(pair[1]);
                if (pair[0] == "pidp")  SetPidProton(pair[1]);
                
                if (pair[0] == "algo") SetPid("",pair[1]);
                
                if (pair[0] == "algoe")  SetPidElectron("",pair[1]);
                if (pair[0] == "algomu") SetPidMuon("",pair[1]);
                if (pair[0] == "algopi") SetPidPion("",pair[1]);
                if (pair[0] == "algok")  SetPidKaon("",pair[1]);
                if (pair[0] == "algop")  SetPidProton("",pair[1]);
                
                // ****
                // E_min or p_min
                // ****
                if (pair[0] == "emin") {fEmin = pair[1].Atof(); fESel = new RhoEnergyParticleSelector("eSel",50.+fEmin,100.);}
                if (pair[0] == "pmin") {fPmin = pair[1].Atof(); fPSel = new RhoMomentumParticleSelector("pSel",50.+fPmin,100.);}
        }
        return kTRUE;
}

void PndSimpleCombiner::Print

(

)

Definition at line 488 of file PndSimpleCombiner.cxx.

References SCDecayInfo::daucc, SCDecayInfo::didx, SCDecayInfo::dpdg, fDecayInfoArray, fEcm, fEmin, fIdxListNameMap, fIdxPidAlgoMap, fIdxPidCritMap, fPmin, i, SCDecayInfo::midx, SCDecayInfo::mpdg, SCDecayInfo::mwin, SCDecayInfo::mwinhi, SCDecayInfo::mwinlo, n, and printf().

Referenced by PndSimpleCombinerTask::Init().

{
        cout <<"\n------------------------------------------"<<endl<<"[PndSimpleCombiner] **** Configuration"<<endl<<"------------------------------------------"<<endl;
        
        cout <<"E_cm           = "<<fEcm<<endl<<endl;
        cout <<"Neutrals E_min = "<<fEmin<<endl;
        cout <<"Charged  p_min = "<<fPmin<<endl<<endl;
        
        // fill all generic lists, which are used by any of the composites
        // loop through composites
        int n = fDecayInfoArray.size();
        for (int i=0;i<n;++i)
        {
                SCDecayInfo info = fDecayInfoArray[i];
                cout <<"Decay "<<i<<" : " << TDatabasePDG::Instance()->GetParticle(info.mpdg)->GetName()<<"("<<info.mpdg<<"/"<<info.midx<<") -> ";
                for (size_t j=0;j<info.dpdg.size();++j) cout << TDatabasePDG::Instance()->GetParticle(info.dpdg[j])->GetName()<<"("<<info.dpdg[j]<<"/"<<info.didx[j]<<") ";
                
                cout <<" //  daucc: "<<info.daucc;
                cout <<"  //  mass window: ";
                if (info.mwin>0) cout <<" ["<<info.mwinlo<<" ; "<<info.mwinhi<<"] = "<<info.mwin<<" GeV/c²"<<endl;
                else if (info.mwin<0) cout <<" *** skipped due to missing Ecm setting *** "<<endl;
                else cout <<"none"<<endl;
                cout <<endl;
        }
        cout <<endl;
        for (int i=0;i<10;++i)
        {
                printf("%-16s : %s, %s\n",fIdxListNameMap[i].Data(),fIdxPidCritMap[i].Data(),fIdxPidAlgoMap[i].Data());
        }
        
        cout <<"\n------------------------------------------\n\n"<<endl;
}

void PndSimpleCombiner::SetPid	(	TString	crit = "",
		TString	algo = ""
	)

Definition at line 362 of file PndSimpleCombiner.cxx.

References fIdxPidAlgoMap, fIdxPidCritMap, SetPidElectron(), SetPidKaon(), SetPidMuon(), SetPidPion(), and SetPidProton().

Referenced by ParseParams(), and PndSimpleCombiner().

{
        if (crit!="") fIdxPidCritMap.clear();
        if (algo!="") fIdxPidAlgoMap.clear();
        
        SetPidElectron(crit, algo);
        SetPidMuon(crit, algo);
        SetPidPion(crit, algo);
        SetPidKaon(crit, algo);
        SetPidProton(crit, algo);
}

void PndSimpleCombiner::SetPidElectron	(	TString	crit = "",
		TString	algo = ""
	)

Definition at line 376 of file PndSimpleCombiner.cxx.

References fIdxPidAlgoMap, and fIdxPidCritMap.

Referenced by ParseParams(), and SetPid().

{
        if (crit!="") { fIdxPidCritMap[0] = crit; fIdxPidCritMap[1] = crit;} 
        if (algo!="") { fIdxPidAlgoMap[0] = algo; fIdxPidAlgoMap[1] = algo;}
}

void PndSimpleCombiner::SetPidKaon	(	TString	crit = "",
		TString	algo = ""
	)

Definition at line 394 of file PndSimpleCombiner.cxx.

References fIdxPidAlgoMap, and fIdxPidCritMap.

Referenced by ParseParams(), and SetPid().

{
        if (crit!="") {fIdxPidCritMap[7] = crit; fIdxPidCritMap[6] = crit;}
        if (algo!="") {fIdxPidAlgoMap[7] = algo; fIdxPidAlgoMap[6] = algo;}
}

void PndSimpleCombiner::SetPidMuon	(	TString	crit = "",
		TString	algo = ""
	)

Definition at line 382 of file PndSimpleCombiner.cxx.

References fIdxPidAlgoMap, and fIdxPidCritMap.

Referenced by ParseParams(), and SetPid().

{
        if (crit!="") {fIdxPidCritMap[2] = crit; fIdxPidCritMap[3] = crit;} 
        if (algo!="") {fIdxPidAlgoMap[2] = algo; fIdxPidAlgoMap[3] = algo;}
}

void PndSimpleCombiner::SetPidPion	(	TString	crit = "",
		TString	algo = ""
	)

Definition at line 388 of file PndSimpleCombiner.cxx.

References fIdxPidAlgoMap, and fIdxPidCritMap.

Referenced by ParseParams(), and SetPid().

{
        if (crit!="") {fIdxPidCritMap[5] = crit; fIdxPidCritMap[4] = crit;} 
        if (algo!="") {fIdxPidAlgoMap[5] = algo; fIdxPidAlgoMap[4] = algo;}
}

void PndSimpleCombiner::SetPidProton	(	TString	crit = "",
		TString	algo = ""
	)

Definition at line 400 of file PndSimpleCombiner.cxx.

References fIdxPidAlgoMap, and fIdxPidCritMap.

Referenced by ParseParams(), and SetPid().

{
        if (crit!="") {fIdxPidCritMap[9] = crit; fIdxPidCritMap[8] = crit;}
        if (algo!="") {fIdxPidAlgoMap[9] = algo; fIdxPidAlgoMap[8] = algo;}
}

void PndSimpleCombiner::SetVerbose ( int  verb = 1 )

inline

Definition at line 69 of file PndSimpleCombiner.h.

References fVerbose.

Referenced by PndSimpleCombinerTask::Init().

{fVerbose = verb;}

int PndSimpleCombiner::SplitString ( TString  s, TString  delim, StringList &  toks  )

private

Definition at line 94 of file PndSimpleCombiner.cxx.

References i, and TString.

Referenced by ParseDecay(), and ParseParams().

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

Member Data Documentation

PndAnalysis* PndSimpleCombiner::fAnalysis

private

Definition at line 95 of file PndSimpleCombiner.h.

Referenced by FillGenericLists().

TString PndSimpleCombiner::fDecay

private

Definition at line 96 of file PndSimpleCombiner.h.

std::vector<SCDecayInfo> PndSimpleCombiner::fDecayInfoArray

private

Definition at line 115 of file PndSimpleCombiner.h.

Referenced by Combine(), FillGenericLists(), ParseDecay(), ParseParams(), Print(), and ~PndSimpleCombiner().

double PndSimpleCombiner::fEcm

private

Definition at line 103 of file PndSimpleCombiner.h.

Referenced by ParseParams(), PndSimpleCombiner(), and Print().

double PndSimpleCombiner::fEmin

private

Definition at line 101 of file PndSimpleCombiner.h.

Referenced by ParseParams(), and Print().

RhoEnergyParticleSelector* PndSimpleCombiner::fESel

private

Definition at line 105 of file PndSimpleCombiner.h.

Referenced by FillGenericLists(), ParseParams(), and ~PndSimpleCombiner().

TString PndSimpleCombiner::fGlobParams

private

Definition at line 97 of file PndSimpleCombiner.h.

std::map<int,TString> PndSimpleCombiner::fIdxListNameMap

private

Definition at line 111 of file PndSimpleCombiner.h.

Referenced by FillGenericLists(), ParseParams(), PndSimpleCombiner(), and Print().

std::map<int,int> PndSimpleCombiner::fIdxPdgMap

private

Definition at line 109 of file PndSimpleCombiner.h.

Referenced by FillGenericLists(), ParseDecay(), and PndSimpleCombiner().

std::map<int,TString> PndSimpleCombiner::fIdxPidAlgoMap

private

Definition at line 113 of file PndSimpleCombiner.h.

Referenced by FillGenericLists(), Print(), SetPid(), SetPidElectron(), SetPidKaon(), SetPidMuon(), SetPidPion(), and SetPidProton().

std::map<int,TString> PndSimpleCombiner::fIdxPidCritMap

private

Definition at line 112 of file PndSimpleCombiner.h.

Referenced by FillGenericLists(), Print(), SetPid(), SetPidElectron(), SetPidKaon(), SetPidMuon(), SetPidPion(), and SetPidProton().

RhoCandList PndSimpleCombiner::fList[MAXLISTS]

private

Definition at line 99 of file PndSimpleCombiner.h.

Referenced by Combine(), CombineList(), FillGenericLists(), GetList(), and GetListN().

int PndSimpleCombiner::fNLists

private

Definition at line 98 of file PndSimpleCombiner.h.

Referenced by GetNLists(), and ParseDecay().

std::map<int,int> PndSimpleCombiner::fPdgIdxMap

private

Definition at line 108 of file PndSimpleCombiner.h.

Referenced by Combine(), FillGenericLists(), GetList(), ParseDecay(), ParseParams(), and PndSimpleCombiner().

double PndSimpleCombiner::fPmin

private

Definition at line 102 of file PndSimpleCombiner.h.

Referenced by ParseParams(), and Print().

RhoMomentumParticleSelector* PndSimpleCombiner::fPSel

private

Definition at line 106 of file PndSimpleCombiner.h.

Referenced by FillGenericLists(), ParseParams(), and ~PndSimpleCombiner().

int PndSimpleCombiner::fVerbose

private

Definition at line 100 of file PndSimpleCombiner.h.

Referenced by FillGenericLists(), and SetVerbose().

---

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

• PndSimpleCombiner.h
• PndSimpleCombiner.cxx