softtrigger_toy12.C File Reference

#include "EventShape.h"
#include <iostream>
#include <vector>
#include <map>
#include "TFile.h"
#include "TTree.h"
#include "TRandom3.h"
#include "TDatabasePDG.h"
#include "TParticlePDG.h"
#include "TCanvas.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TLegend.h"
#include "TParticle.h"
#include "TClonesArray.h"
#include "TF1.h"
#include "TStyle.h"
#include "TNtuple.h"
#include "TLatex.h"

Go to the source code of this file.

Typedefs
typedef std::vector< SimpleCand >	CandList

Functions
void	boostList (CandList &list, CandList &boostlist, TVector3 boost)
bool	isFinalState (int n)
double	dponline (TLorentzVector &v)
void	makePidTable (double pideff, double misid)
void	init (double pideff, double misid)
void	printCand (SimpleCand c)
void	printList (CandList &l, TString tit="")
void	printAllLists ()
int	massCrit (CandList &l, double m, double w)
bool	eventAccepted (CandList &l, double m=0, double w=0)
int	fillHisto (CandList &l, TH1 *hall, TH1 *hbg=0, double m=0, double w=0, TH1 *hallsel=0, TH1 *hbgsel=0, TH1 *hsig=0, TH1 *hmm=0, TH2 *hmvsmm=0)
int	combine (CandList &l1, CandList &l2, CandList &out, int matchPdg=0)
int	combine (CandList &l1, CandList &l2, CandList &l3, CandList &out, int matchPdg=0)
int	combine (CandList &l1, CandList &l2, CandList &l3, CandList &l4, CandList &out, int matchPdg=0)
void	smearMom (SimpleCand &c, double dpr=0.05, double dtht=0.000, double dphi=0.000)
bool	isDetected (double eff=0.90)
void	makeRecoCands (CandList &mc, CandList &reco, double eff=0.90, double dp=0.05, double dtht=0.001, double dphi=0.001)
void	selectMass (CandList &in, CandList &out, double mmean=0, double mw=0)
void	select (CandList &in, CandList &out, int chrg=0, int pdg=0, double mass=0.139)
void	makePidSelection (CandList &l)
void	makeIni4Vector (TLorentzVector &l, double s)
void	configHisto (TH1 *h)
int	softtrigger_toy12 (TString fsig, int nev=2, double sqrts=3.77, bool simcut=false, bool evcut=false, double dp=0.03, double trkeff=100., double pideff=95., double misid=5., double P_mix=0.00)

Variables
TRandom3	fRand
TLorentzVector	fIni
const int	fMAX = 100
int	nbins = 300
double	low = 0.0
double	high = 5.0
const double	fNeutEff = 1.0
const double	fdE = 0.05
const double	fNdtht = 0.003
const double	fNdphi = 0.003
int	fPi0Pdg = 111
double	fPi0Mass = 0.13497
double	fPi0win = 0.020
int	fKsPdg = 310
double	fKsMass = 0.497614
double	fKswin = 0.025
int	fDPdg = 411
double	fDMass = 1.8693
double	fDsig = 0.88
double	fDwin =0.25
int	fD0Pdg = 421
double	fD0Mass = 1.8645
double	fD0sig = 1.04
double	fD0win =0.25
int	fJpsiPdg = 443
double	fJMass = 3.09687
double	fJsig = 2.24
double	fJsig3pi = 1.5
double	fJ3pioff = 0.027
double	fJwin =0.7
int	fLamPdg = 3122
double	fLamMass = 1.115684
double	fLamsig = 0.08
double	fLamwin =0.1
int	fDsPdg = 431
double	fDsMass = 1.9686
double	fDssig = 0.75
double	fDswin =0.2
int	fPhiPdg = 333
double	fPhiMass = 1.0195
double	fPhisig = 0.15
double	fPhiwin =0.2
int	fEtacPdg = 441
double	fEtacMass = 2.9798
double	fEtacsig = 0.15
double	fEtacwin =0.15
int	fLamcPdg = 4122
double	fLamcMass = 2.2849
double	fLamcsig = 0.15
double	fLamcwin =0.2
double	pidtab [25]
double	pidreal [25]
double	mE
double	mMu
double	mPi
double	mK
double	mP
double	mPi0
TF1 *	f_res =0
std::map< int, int >	pdgidx
TDatabasePDG *	pdg_db
CandList	mclist
CandList	mcFSlist
CandList	recolist
CandList	recocmslist
CandList	eplus
CandList	eminus
CandList	muplus
CandList	muminus
CandList	piplus
CandList	piminus
CandList	kplus
CandList	kminus
CandList	pplus
CandList	pminus
CandList	gam

Typedef Documentation

typedef std::vector<SimpleCand> CandList

Definition at line 118 of file softtrigger_toy12.C.

Function Documentation

void boostList	(	CandList &	list,
		CandList &	boostlist,
		TVector3	boost
	)

Definition at line 127 of file softtrigger_toy12.C.

References c, i, SimpleCand::P4(), and SimpleCand::SetP4().

{
        int N=list.size();
        boostlist.clear();

    for (unsigned int i=0;i<N;i++)
    {
        SimpleCand c=list[i];
                TLorentzVector lv(c.P4());
        lv.Boost(-boost);
                c.SetP4(lv);
                boostlist.push_back(c);
    }
}

int combine	(	CandList &	l1,
		CandList &	l2,
		CandList &	out,
		int	matchPdg = 0
	)

Definition at line 326 of file softtrigger_toy12.C.

References c, i, m2(), mclist, SimpleCand::SetCharge(), SimpleCand::SetEvtId(), SimpleCand::SetId(), SimpleCand::SetMarker(), SimpleCand::SetMcPid(), SimpleCand::SetMct(), SimpleCand::SetMotherIdx(), SimpleCand::SetNFS(), and SimpleCand::SetP4().

{
        out.clear();
        
        bool sameList = (&l1 == &l2);
        
        for (unsigned int i=0;i<l1.size();++i)
        {
                for (unsigned int j=0;j<l2.size();++j)
                {
                        if (sameList && j<=i) continue;
                        if (l1[i].Marker() & l2[j].Marker()) continue;  // *** overlap
                        //if (l1[i].EvtId()!=l2[j].EvtId()) continue;     // *** tracks from different events
                        SimpleCand c;
                        c.SetP4(l1[i].P4()+l2[j].P4());
                        c.SetCharge(l1[i].Charge()+l2[j].Charge());
                        c.SetNFS(l1[i].NFS()+l2[j].NFS());
                        c.SetMarker(l1[i].Marker() | l2[j].Marker());
                        
                        int m1 = l1[i].MotherIdx();
                        int m2 = l2[j].MotherIdx();
                        int ev1 = l1[i].EvtId();
                        int ev2 = l2[j].EvtId();
                        
                        if ( ev1==ev2 &&                                                           // *** particles are from same event
                                 m1>=0 && m1 == m2 &&                      // *** there is a common mother of all particles
                            (l1[i].McPid()||l1[i].Mct()) &&            // *** all particles match (FS: PID, Comp:mct)
                            (l2[j].McPid()||l2[j].Mct()) &&            // *** "
                                    2==mclist[m1].NDau() &&            // *** mother has as many daughters as we combine particles
                            matchPdg == mclist[m1].Pid() )             // *** mother has correct PDG
                        {
                                c.SetId(m1);
                                c.SetMcPid();
                                c.SetMct();
                                c.SetMotherIdx(mclist[m1].MotherIdx());
                                c.SetEvtId(ev1);
                        }
                        
                        out.push_back(c);
                }
        }
        return out.size();
}

int combine	(	CandList &	l1,
		CandList &	l2,
		CandList &	l3,
		CandList &	out,
		int	matchPdg = 0
	)

Definition at line 370 of file softtrigger_toy12.C.

References c, m2(), mclist, SimpleCand::SetCharge(), SimpleCand::SetEvtId(), SimpleCand::SetId(), SimpleCand::SetMarker(), SimpleCand::SetMcPid(), SimpleCand::SetMct(), SimpleCand::SetMotherIdx(), SimpleCand::SetNFS(), and SimpleCand::SetP4().

{
        out.clear();
        
        bool same12 = (&l1 == &l2);
        bool same13 = (&l1 == &l3);
        bool same23 = (&l2 == &l3);
        
        for (unsigned int i1=0;i1<l1.size();++i1)
        {
                for (unsigned int i2=0;i2<l2.size();++i2)
                {
                        if (same12 && i2<i1) continue;
                        if (l1[i1].Marker() & l2[i2].Marker()) continue;
                        
                        for (unsigned int i3=0;i3<l3.size();++i3)
                        {
                                if ((same13 && i3<i1) || (same23 && i3<i2)) continue;
                                
                                if ( l1[i1].Marker() & l3[i3].Marker() ) continue;
                                if ( l2[i2].Marker() & l3[i3].Marker() ) continue;
                                
                                //if (l1[i1].EvtId()!=l2[i2].EvtId() || 
                                //    l1[i1].EvtId()!=l3[i3].EvtId() ||
                                //    l2[i2].EvtId()!=l3[i3].EvtId()   ) continue;     // *** tracks from different events
                                
                                SimpleCand c;
                                c.SetP4( l1[i1].P4() + l2[i2].P4() + l3[i3].P4());
                                c.SetCharge( l1[i1].Charge() + l2[i2].Charge() + l3[i3].Charge() );
                                c.SetNFS( l1[i1].NFS() + l2[i2].NFS() + l3[i3].NFS());
                                c.SetMarker( l1[i1].Marker() | l2[i2].Marker() | l3[i3].Marker() );
                                
                                int m1 = l1[i1].MotherIdx();
                                int m2 = l2[i2].MotherIdx();
                                int m3 = l3[i3].MotherIdx();
                                
                                int ev1 = l1[i1].EvtId();
                                int ev2 = l2[i2].EvtId();
                                int ev3 = l3[i3].EvtId();
                                
                                if ( ev1==ev2 && ev2==ev3 &&                               // *** particles are from same event
                                         m1>=0 && m1 == m2 && m2 == m3 &&          // *** there is a common mother of all particles
                                        (l1[i1].McPid()||l1[i1].Mct()) &&          // *** all particles match (FS: PID, Comp:mct)
                                        (l2[i2].McPid()||l2[i2].Mct()) &&          // *** "
                                        (l3[i3].McPid()||l3[i3].Mct()) &&          // *** "
                                                        3==mclist[m1].NDau() &&            // *** mother has as many daughters as we combine particles
                                        matchPdg == mclist[m1].Pid() )             // *** mother has correct PDG
                                {
                                        c.SetId(m1);
                                        c.SetMcPid();
                                        c.SetMct();
                                        c.SetMotherIdx(mclist[m1].MotherIdx());
                                        c.SetEvtId(ev1);
                                }
                                
                                out.push_back(c);
                        }
                }
        }
        
        return out.size();
}

int combine	(	CandList &	l1,
		CandList &	l2,
		CandList &	l3,
		CandList &	l4,
		CandList &	out,
		int	matchPdg = 0
	)

Definition at line 433 of file softtrigger_toy12.C.

References c, m2(), mclist, SimpleCand::SetCharge(), SimpleCand::SetEvtId(), SimpleCand::SetId(), SimpleCand::SetMarker(), SimpleCand::SetMcPid(), SimpleCand::SetMct(), SimpleCand::SetMotherIdx(), SimpleCand::SetNFS(), and SimpleCand::SetP4().

{
        out.clear();
        
        bool same12 = (&l1 == &l2);
        bool same13 = (&l1 == &l3);
        bool same14 = (&l1 == &l4);
        
        bool same23 = (&l2 == &l3);
        bool same24 = (&l2 == &l4);
        
        bool same34 = (&l3 == &l4);
        
        for (unsigned int i1=0;i1<l1.size();++i1)
        {
                for (unsigned int i2=0;i2<l2.size();++i2)
                {
                        if (same12 && i2<i1) continue;
                        if (l1[i1].Marker() & l2[i2].Marker()) continue;
                        
                        for (unsigned int i3=0;i3<l3.size();++i3)
                        {
                                if ((same13 && i3<i1) || (same23 && i3<i2)) continue;
                                
                                if ( l1[i1].Marker() & l3[i3].Marker() ) continue;
                                if ( l2[i2].Marker() & l3[i3].Marker() ) continue;
                                
                                for (unsigned int i4=0;i4<l4.size();++i4)
                                {
                                        if ((same14 && i4<i1) || (same24 && i4<i2) || (same34 && i4<i3)) continue;
                                        
                                        if ( l1[i1].Marker() & l4[i4].Marker() ) continue;
                                        if ( l2[i2].Marker() & l4[i4].Marker() ) continue;
                                        if ( l3[i3].Marker() & l4[i4].Marker() ) continue;
                                        
                                        //if (l1[i1].EvtId()!=l2[i2].EvtId() || 
                                        //    l1[i1].EvtId()!=l3[i3].EvtId() ||
                                        //    l2[i2].EvtId()!=l3[i3].EvtId()   ) continue;     // *** tracks from different events
                                        
                                        SimpleCand c;
                                        c.SetP4( l1[i1].P4() + l2[i2].P4() + l3[i3].P4() + l4[i4].P4() );
                                        c.SetCharge( l1[i1].Charge() + l2[i2].Charge() + l3[i3].Charge() + l4[i4].Charge() );
                                        c.SetNFS( l1[i1].NFS() + l2[i2].NFS() + l3[i3].NFS() + l4[i4].NFS() );
                                        c.SetMarker( l1[i1].Marker() | l2[i2].Marker() | l3[i3].Marker() | l4[i4].Marker() );
                                        
                                        int m1 = l1[i1].MotherIdx();
                                        int m2 = l2[i2].MotherIdx();
                                        int m3 = l3[i3].MotherIdx();
                                        int m4 = l4[i4].MotherIdx();
                                        
                                        int ev1 = l1[i1].EvtId();
                                        int ev2 = l2[i2].EvtId();
                                        int ev3 = l3[i3].EvtId();
                                        int ev4 = l4[i4].EvtId();
                                        
                                        if ( ev1==ev2 && ev2==ev3 && ev3==ev4 &&                        // *** particles are from same event
                                                m1>=0 && m1 == m2 && m2 == m3 && m3 == m4 &&    // *** there is a common mother of all particles
                                                (l1[i1].McPid()||l1[i1].Mct()) &&                               // *** all particles match (FS: PID, Comp:mct)
                                                (l2[i2].McPid()||l2[i2].Mct()) &&                               // *** "
                                                (l3[i3].McPid()||l3[i3].Mct()) &&                               // *** "
                                                (l4[i4].McPid()||l4[i4].Mct()) &&                               // *** "
                                                                4==mclist[m1].NDau() &&                                 // *** mother has as many daughters as we combine particles
                                                matchPdg == mclist[m1].Pid() )                                  // *** mother has correct PDG
                                        {
                                                c.SetId(m1);
                                                c.SetMcPid();
                                                c.SetMct();
                                                c.SetMotherIdx(mclist[m1].MotherIdx());
                                                c.SetEvtId(ev1);
                                        }
                                        
                                        out.push_back(c);
                                }
                        }
                }
        }
        
        return out.size();
}

void configHisto

(

TH1 *

h

)

Definition at line 652 of file softtrigger_toy12.C.

{
        h->SetFillStyle(3001);
        h->SetFillColor(602);
        h->SetLineStyle(3);
        h->SetLineColor(602);
}

double dponline

(

TLorentzVector &

v

)

Definition at line 152 of file softtrigger_toy12.C.

References f_res.

{       
        return f_res->Eval(v.Pt())/100.;
}

bool eventAccepted	(	CandList &	l,
		double	m = 0,
		double	w = 0
	)

Definition at line 283 of file softtrigger_toy12.C.

References i, and m.

Referenced by softtrigger_toy12(), and toy_core().

{
        bool accepted = false;
        unsigned int i=0;
        
        while (!accepted && i<l.size())
        //for (unsigned int i=0;i<l.size();++i) 
        {
                double mass = l[i++].P4().M();
                if (mass>m-w && mass<m+w) accepted =true;
        }
        
        return accepted;
}

int fillHisto	(	CandList &	l,
		TH1 *	hall,
		TH1 *	hbg = 0,
		double	m = 0,
		double	w = 0,
		TH1 *	hallsel = 0,
		TH1 *	hbgsel = 0,
		TH1 *	hsig = 0,
		TH1 *	hmm = 0,
		TH2 *	hmvsmm = 0
	)

Definition at line 300 of file softtrigger_toy12.C.

References cnt, fIni, i, m, and mct.

{
        int cnt =0;
        for (unsigned int i=0;i<l.size();++i) 
        {
                double mass = l[i].P4().M();
                double miss = (fIni-l[i].P4()).M();
                hall->Fill(mass);
                
                //missing mass
                if (hmm) hmm->Fill(miss);
                if (hmvsmm) hmvsmm->Fill(mass, miss);
                
                bool mct = l[i].Mct(); 
                if (hbg && !mct) hbg->Fill(mass);
                if (hsig && mct) hsig->Fill(mass);
                if (mass>m-w && mass<m+w) 
                {
                        cnt++;
                        if (hallsel) hallsel->Fill(mass);
                        if (hbgsel && !mct) hbgsel->Fill(mass);
                }
        }
        return cnt;
}

void init	(	double	pideff,
		double	misid
	)

Definition at line 198 of file softtrigger_toy12.C.

References f_res, makePidTable(), mE, mK, mMu, mP, mPi, mPi0, pdg_db, and pdgidx.

{
        // *** setting the mapping of PDG code to tab index
        pdgidx.clear();
        pdgidx[11]=0;
        pdgidx[13]=1;
        pdgidx[211]=2;
        pdgidx[321]=3;
        pdgidx[2212]=4;
        
        pdg_db=TDatabasePDG::Instance();
        
        // *** create pid table
        makePidTable(pideff, misid);
        
        // *** gets particle masses from pdg database and set them globally
        mE  = pdg_db->GetParticle(11)->Mass();
        mMu = pdg_db->GetParticle(13)->Mass();
        mPi = pdg_db->GetParticle(211)->Mass();
        mK  = pdg_db->GetParticle(321)->Mass();
        mP  = pdg_db->GetParticle(2212)->Mass();
        
        mPi0= pdg_db->GetParticle(111)->Mass();
        
        // *** resolution fcn based on Yutie Lang results; see routine dponline above
        f_res=new TF1("f_res","pol2(0)");
        f_res->SetParameters(2.367,1.3,0.133);
        f_res->SetRange(0,10);
}

bool isDetected

(

double

eff = 0.90

)

Definition at line 544 of file softtrigger_toy12.C.

References fRand.

{
        // very simple efficiency cut; might be dependent on candidate itself later
        if (fRand.Rndm()<eff) return true;
        return false;
}

bool isFinalState

(

int

n

)

Definition at line 143 of file softtrigger_toy12.C.

{
        n=abs(n);
        return (n==11 || n==13  || n==211 || n==321 || n==2212 || n==22);
}

void makeIni4Vector	(	TLorentzVector &	l,
		double	s
	)

Definition at line 645 of file softtrigger_toy12.C.

References mP, p, sqrt(), and X.

{
        double X = (s*s-2*mP*mP)/(2*mP);
        double p = sqrt(X*X - mP*mP);
        l.SetXYZM(0,0,p,s);
}

void makePidSelection

(

CandList &

l

)

Definition at line 624 of file softtrigger_toy12.C.

References eminus, eplus, gam, kminus, kplus, mE, mK, mMu, mP, mPi, muminus, muplus, piminus, piplus, pminus, pplus, and select().

{
        // *** prepare globally defined PID lists
        select(l, eplus, 1, -11, mE);
        select(l, eminus, -1, 11, mE);
        
        select(l, muplus, 1, -13, mMu);
        select(l, muminus, -1, 13, mMu);
        
        select(l, piplus, 1, 211, mPi);
        select(l, piminus, -1, -211, mPi);
        
        select(l, kplus, 1, 321, mK);
        select(l, kminus, -1, -321, mK);
        
        select(l, pplus, 1, 2212, mP);
        select(l, pminus, -1, -2212, mP);
        
        select(l, gam, 0, 22, 0);
}

void makePidTable	(	double	pideff,
		double	misid
	)

Definition at line 158 of file softtrigger_toy12.C.

References i, pidtab, and printf().

{
        cout<< "Using PID table:"<<endl;
        
        if (pideff<0 && misid<0)
        {
                for (int i=0;i<5;++i)
                {
                        for (int j=0;j<5;++j)
                        {
                                //pidtab[j+5*i] = pidreal[j+5*i];
                                printf("%4.0f",pidtab[j+5*i]);
                                //cout << pidtab[i+5*j]<<"  ";
                        }
                        cout<< endl;
                }
                return;
        }
        
        if (pideff<0)   pideff=0;
        if (pideff>100) pideff=100;
        if (misid>100) misid=100;
        
        // *** neff is the efficiency for particles of wrong type
        //double neff = pideff*(100-pidpur)/(4*pidpur);
        //if (neff>100) neff=100;
        
        for (int i=0;i<5;++i)
        {
                for (int j=0;j<5;++j)
                {
                        if (i==j) pidtab[i+5*j] = pideff;
                        else pidtab[i+5*j] = misid;
                        cout << pidtab[i+5*j]<<"  ";
                }
                cout<< endl;
        }
}

void makeRecoCands	(	CandList &	mc,
		CandList &	reco,
		double	eff = 0.90,
		double	dp = 0.05,
		double	dtht = 0.001,
		double	dphi = 0.001
	)

Definition at line 551 of file softtrigger_toy12.C.

References c, fdE, fNdphi, fNdtht, fNeutEff, i, isDetected(), isFinalState(), pid(), SimpleCand::Pid(), SimpleCand::SetMass(), and smearMom().

{
        reco.clear();
        
        for (unsigned int i=0;i<mc.size();i++)
        {
                SimpleCand c=mc[i];
                int pid = abs(c.Pid());
                if (!isFinalState(pid)) continue; 
                
                if (pid!=22 && isDetected(eff))  // *** charged track
                {
                        smearMom(c,dp,dtht,dphi);
                        // set pion mass as default
                        c.SetMass(0.13957);
                        reco.push_back(c);
                }
                else if (pid==22 && isDetected(fNeutEff))  // *** neutral
                {
                        smearMom(c,fdE, fNdtht, fNdphi);
                        reco.push_back(c);
                }
        }
}

int massCrit	(	CandList &	l,
		double	m,
		double	w
	)

Definition at line 270 of file softtrigger_toy12.C.

References cnt, and i.

{
        int cnt = 0;
        for (unsigned int i=0;i<l.size();++i)
        {
                double mass=l[i].P4().M();
                if (mass>m-w && mass<m+w) cnt++;
        }
        
        return cnt;
}

void printAllLists

(

)

Definition at line 246 of file softtrigger_toy12.C.

References eminus, eplus, gam, kminus, kplus, mclist, muminus, muplus, piminus, piplus, pminus, pplus, printList(), and recolist.

{
        printList(mclist, "MC List");
        printList(recolist, "Reco List");
        
        printList(eplus, "e+ List");
        printList(eminus, "e- List");
        
        printList(muplus, "mu+ List");
        printList(muminus, "mu- List");
        
        printList(piplus, "pi+ List");
        printList(piminus, "pi- List");
        
        printList(kplus, "k+ List");
        printList(kminus, "k- List");
        
        printList(pplus, "p+ List");
        printList(pminus, "p- List");
        
        printList(gam, "gamma List");
}

void printCand

(

SimpleCand

c

)

Definition at line 229 of file softtrigger_toy12.C.

References SimpleCand::Charge(), SimpleCand::Id(), SimpleCand::McPid(), SimpleCand::Mct(), SimpleCand::MotherIdx(), SimpleCand::NFS(), SimpleCand::NSiblings(), SimpleCand::P4(), and SimpleCand::Pid().

{
        cout <<c.Id()<<" pid:"<<c.Pid()<<" ("<<c.P4().X()<<","<<c.P4().Y()<<","<<c.P4().Z()<<";"<<c.P4().T()<<") ";
        cout <<c.Charge()<<" "<<c.P4().M();
        cout <<" moth:"<<c.MotherIdx()<<" nsib:"<<c.NSiblings()<<" nfs:"<<c.NFS();
        cout <<" mcpid:"<<c.McPid()<<" mct:"<<c.Mct()<<endl;
}

void printList	(	CandList &	l,
		TString	tit = ""
	)

Definition at line 238 of file softtrigger_toy12.C.

References printCand().

{
        if (tit!="") cout << tit<<endl;
        for (unsigned int j=0;j<l.size();++j) printCand(l[j]);
        cout <<endl;
}

void select	(	CandList &	in,
		CandList &	out,
		int	chrg = 0,
		int	pdg = 0,
		double	mass = 0.139
	)

Definition at line 587 of file softtrigger_toy12.C.

References c, SimpleCand::Charge(), fRand, i, pdg_db, pdgidx, SimpleCand::Pid(), pidtab, SimpleCand::SetMass(), and SimpleCand::SetMcPid().

{
        out.clear();
        
        TParticlePDG *part = pdg_db->GetParticle(pdg);
        
        //pdg *= pdg<20 ? -chrg : chrg;
        
        if (part) mass = part->Mass();  
        if (mass<0) mass = 0.139;

        for (unsigned int i=0;i<in.size();++i)
        {
                SimpleCand c=in[i];
                
                // select a specific charge in any case
                if (c.Charge()==chrg)
                {
                        c.SetMass(mass);
                        if (c.Pid()==pdg) c.SetMcPid();
                        if (pdg==0 || pdg==22) out.push_back(c); // do we want select a pdg code? 0=only charge
                        else
                        {
                            // invalid pdg code?
                                if (pdgidx.find(abs(pdg))==pdgidx.end()) continue;
                                
                                int selidx   = pdgidx[abs(pdg)];
                                int partidx  = pdgidx[abs(c.Pid())];
                                double prob  = pidtab[5*selidx + partidx]/100.; 
                                //cout <<pdg <<" "<<c.Pid()<<" : "<<selidx<<" "<<partidx<<" P:"<<prob<<endl;
                                
                                if (fRand.Rndm()<prob) out.push_back(c);
                        }
                }
        }
}

void selectMass	(	CandList &	in,
		CandList &	out,
		double	mmean = 0,
		double	mw = 0
	)

Definition at line 577 of file softtrigger_toy12.C.

References fabs(), and i.

{
        out.clear();
        
        for (unsigned int i=0;i<in.size();++i)
                if (fabs(in[i].P4().M()-mmean)<mw) out.push_back(in[i]);
   
}

void smearMom	(	SimpleCand &	c,
		double	dpr = 0.05,
		double	dtht = 0.000,
		double	dphi = 0.000
	)

Definition at line 514 of file softtrigger_toy12.C.

References dponline(), fRand, m, p, SimpleCand::P4(), and SimpleCand::SetP4().

{
  TLorentzVector p4=c.P4();
 
  // *** use resolution function
  if (dpr<0) dpr = dponline(p4);
  
  double p    = p4.P();
  double m    = p4.M();
  
  double dp   = p*dpr;
  double newp = -1;
  
  if (dp==0) newp=p;
  while (newp<0) newp = p4.P() + fRand.Gaus(0.0,dp);
  
  p4.SetVectM( p4.Vect().Unit()*newp, m );
  
  double newtht = p4.Theta();
  if (dtht!=0) newtht += fRand.Gaus(0.,dtht);
  p4.SetTheta(newtht);
  
  double newphi = p4.Phi();
  if (dphi!=0) newphi += fRand.Gaus(0.,dphi); 
  p4.SetPhi(newphi);
  
  c.SetP4(p4);
        
}

int softtrigger_toy12	(	TString	fsig,
		int	nev = 2,
		double	sqrts = 3.77,
		bool	simcut = false,
		bool	evcut = false,
		double	dp = 0.03,
		double	trkeff = 100.,
		double	pideff = 95.,
		double	misid = 5.,
		double	P_mix = 0.00
	)

Definition at line 660 of file softtrigger_toy12.C.

References boostList(), c, c2, EventShape::ChrgPSumCms(), EventShape::ChrgPtSumLab(), combine(), configHisto(), Double_t, eminus, eplus, eventAccepted(), f, fD0Mass, fD0Pdg, fD0win, fDMass, fDPdg, fDsMass, fDsPdg, fDswin, fDwin, fillHisto(), fIni, fJMass, fJpsiPdg, fJwin, fKsMass, fKsPdg, fKswin, fLamcMass, fLamcPdg, fLamcwin, fLamMass, fLamPdg, fLamwin, fMAX, EventShape::FoxWolfMomR(), fPhiMass, fPhiPdg, fPhiwin, fPi0Mass, fPi0Pdg, fPi0win, fRand, gam, high, i, init, isFinalState(), kminus, kplus, m, makeIni4Vector(), makePidSelection(), makeRecoCands(), max(), mcFSlist, mclist, min(), EventShape::MultPminCms(), muminus, muplus, EventShape::NCharged(), EventShape::NeutESumCms(), nk, nlep, EventShape::NNeutral(), EventShape::NParticles(), npi, npr, offset(), p, pdg_db, piminus, piplus, EventShape::PmaxCms(), EventShape::PmaxLab(), EventShape::PminCms(), pminus, pplus, pt(), EventShape::Ptmax(), EventShape::PtSumLab(), pz, recocmslist, recolist, selectMass(), SimpleCand::SetDau(), SimpleCand::SetEvtId(), SimpleCand::SetMotherIdx(), SimpleCand::SetNSiblings(), EventShape::SumChrgPminCms(), t, and TString.

{
        fRand.SetSeed();
        
        init(pideff, misid);
        
        int i,j,k;
        
        makeIni4Vector(fIni, sqrts);
        
        fIni.Print();
        
        //gStyle->SetOptStat(1111);
        
        // *** open file and get tree
        TFile *f=new TFile(fsig,"READ");
        if (f->IsZombie()) 
        {
                cout <<"File not found:"<<fsig<<endl;
                return;
        }
        
        bool dpm=fsig.Contains("DPM"); // *** DPM events
        TString treename = dpm ? "data" : "ntp";  // *** set treename
        TTree *t=(TTree*)f->Get(treename);
                
        if (t==0)
        {
                cout <<"Tree '"<<treename<<"' not found in '"<<fsig<<"'."<<endl;
                return;         
        }
        
        // *** how many events ?
        int Nevt = t->GetEntriesFast();
        if (nev==0) nev = Nevt;
        
        double mmean=0, mwsig=0, mwoff=0;
        int matchPDG=0;
        
        // *** choose histogram scale according to resonance
        high = mmean+2;
        
        TCanvas *c2=new TCanvas("c2","c2",10,10,1500,830);
        c2->Divide(5,3);
        
/*      TCanvas *c3=new TCanvas("c3","c3",1010,10,600,600);
        c3->Divide(2,2);*/
        
//      TCanvas *c4=new TCanvas("c4","c4",1220,20,500,500);
        
        double min=0, max=5.5;
        
        // *** create some histograms
        TH1F *h_mom= new TH1F("h_mom","Momenta",200,0,5);
        TH1F *h_momc= new TH1F("h_momc","Momenta charged",200,0,5);
        TH1F *h_momn= new TH1F("h_momn","Momenta neutral",200,0,5);
        TH1F *h_mult= new TH1F("h_mult","Multiplicity",50,-0.5,49.5);
        TH1F *h_multrec= new TH1F("h_multrec","Multiplicity reco",50,-0.5,49.5);
        TH1F *h_multc= new TH1F("h_multc","Multiplicity charged",30,-0.5,29.5);
        TH1F *h_multn= new TH1F("h_multn","Multiplicity neutral",50,-0.5,49.5);
        TH1F *h_fw1= new TH1F("h_fw1","Fox Wolfram R1",300,-1.1,1.1);
        TH1F *h_pmax= new TH1F("h_pmax","pmax cms",300,0,5);
        
        TH1F *h_mult_k = new TH1F("h_mult_k","Multiplicity Kaons",10,-0.5,9.5);
        TH1F *h_mult_lep = new TH1F("h_mult_lep","Multiplicity Leptons",10,-0.5,9.5);
        TH1F *h_sum_klep = new TH1F("h_sum_lep","Sum Kaons + Leptons",15,-0.5,14.5);
        TH2F *h_mult_klep=new TH2F("h_mult_klep","#Kaons vs. #Leptons",10,-0.5,9.5, 10, -0.5, 9.5);
        h_mult_klep->SetYTitle("#Leptons");
        h_mult_klep->SetXTitle("#Kaons");
        
        // ****** pi0
        TH1F *h_pi0=new TH1F("h_pi0","#pi^{0} #rightarrow #gamma#gamma",200,min,0.4);
        TH1F *h_pi0sel=new TH1F("h_pi0sel","Pi0",200,min,0.4);
        TH1F *h_pi0sig=new TH1F("h_pi0sig","Pi0",200,min,0.4);
        h_pi0sig->SetLineColor(2);
        configHisto(h_pi0sel);
        
        // ****** Ks
        TH1F *h_ks=new TH1F("h_ks","K_{S}^{0} #rightarrow #pi^{+}#pi^{-}",200,0.2,0.8);
        TH1F *h_kssel=new TH1F("h_kssel","Pi0",200,0.2,0.8);
        TH1F *h_kssig=new TH1F("h_kssig","Pi0",200,0.2,0.8);
        h_kssig->SetLineColor(2);
        configHisto(h_kssel);
        
        // ****** pi0
        TH1F *h_jpsi=new TH1F("h_jpsi","J/#psi #rightarrow l^{+}l^{-}",200,min,max);
        TH1F *h_jpsisel=new TH1F("h_jpsisel","J/psi",200,min,max);
        TH1F *h_jpsisig=new TH1F("h_jpsisig","J/psi",200,min,max);
        h_jpsisig->SetLineColor(2);
        configHisto(h_jpsisel);
        
        // ****** D0
        TH1F *h_d0=new TH1F("h_d0","D^{0} #rightarrow K^{-}#pi^{+}",200,min,max);
        TH1F *h_d0sel=new TH1F("h_d0sel","D0",200,min,max);
        TH1F *h_d0sig=new TH1F("h_d0sig","D0",200,min,max);
        h_d0sig->SetLineColor(2);
        configHisto(h_d0sel);
        
        TH1F *h_d02=new TH1F("h_d02","D^{0} #rightarrow K^{-}#pi^{+}#pi^{0}",200,min,max);
        TH1F *h_d02sel=new TH1F("h_d02sel","D0",200,min,max);
        TH1F *h_d02sig=new TH1F("h_d02sig","D0",200,min,max);
        h_d02sig->SetLineColor(2);
        configHisto(h_d02sel);
        
        TH1F *h_d03=new TH1F("h_d03","D^{0} #rightarrow K^{-}#pi^{+}#pi^{+}#pi^{-}",200,min,max);
        TH1F *h_d03sel=new TH1F("h_d03sel","D0",200,min,max);
        TH1F *h_d03sig=new TH1F("h_d03sig","D0",200,min,max);
        h_d03sig->SetLineColor(2);
        configHisto(h_d03sel);
        
        // ****** D+
        TH1F *h_dpm=new TH1F("h_dpm","D^{+} #rightarrow K^{-}#pi^{+}#pi^{+}",200,min,max);
        TH1F *h_dpmsel=new TH1F("h_dpmsel","D#pm",200,min,max);
        TH1F *h_dpmsig=new TH1F("h_dpmsig","D#pm",200,min,max);
        h_dpmsig->SetLineColor(2);
        configHisto(h_dpmsel);
        
        TH1F *h_dpm2=new TH1F("h_dpm2","D^{+} #rightarrow K^{-}#pi^{+}#pi^{+}#pi^{0}",200,min,max);
        TH1F *h_dpm2sel=new TH1F("h_dpm2sel","D#pm",200,min,max);
        TH1F *h_dpm2sig=new TH1F("h_dpm2sig","D#pm",200,min,max);
        h_dpm2sig->SetLineColor(2);
        configHisto(h_dpm2sel);
        
        TH1F *h_dpm3=new TH1F("h_dpm3","D^{+} #rightarrow K_{S}^{0}#pi^{+}#pi^{0}",200,min,max);
        TH1F *h_dpm3sel=new TH1F("h_dpm3sel","D#pm",200,min,max);
        TH1F *h_dpm3sig=new TH1F("h_dpm3sig","D#pm",200,min,max);
        h_dpm3sig->SetLineColor(2);
        configHisto(h_dpm3sel);
        
        TH1F *h_dpm4=new TH1F("h_dpm4","D^{+} #rightarrow K_{S}^{0}#pi^{+}#pi^{+}#pi^{-}",200,min,max);
        TH1F *h_dpm4sel=new TH1F("h_dpm4sel","D#pm",200,min,max);
        TH1F *h_dpm4sig=new TH1F("h_dpm4sig","D#pm",200,min,max);
        h_dpm4sig->SetLineColor(2);
        configHisto(h_dpm4sel);
        
        // ****** Ds+
        TH1F *h_ds=new TH1F("h_ds","D_{s}^{+} #rightarrow K^{+}K^{-}#pi^{+}",200,min,max);
        TH1F *h_dssel=new TH1F("h_dssel","Ds",200,min,max);
        TH1F *h_dssig=new TH1F("h_dssig","Ds",200,min,max);
        h_dssig->SetLineColor(2);
        configHisto(h_dssel);
        
        TH1F *h_ds2=new TH1F("h_ds2","D_{s}^{+} #rightarrow K^{+}K^{-}#pi^{+}#pi^{0}",200,min,max);
        TH1F *h_ds2sel=new TH1F("h_ds2sel","Ds",200,min,max);
        TH1F *h_ds2sig=new TH1F("h_ds2sig","Ds",200,min,max);
        h_ds2sig->SetLineColor(2);
        configHisto(h_ds2sel);
        
        // ****** phi
        TH1F *h_phi=new TH1F("h_phi","#phi #rightarrow K^{+}K^{-}",200,min,max);
        TH1F *h_phisel=new TH1F("h_phisel","Phi",200,min,max);
        TH1F *h_phisig=new TH1F("h_phisig","Phi",200,min,max);
        h_phisig->SetLineColor(2);
        configHisto(h_phisel);
        
        // ****** Lambda_c
        TH1F *h_lamc=new TH1F("h_lamc","#Lambda_{c} #rightarrow pK^{-}#pi^{+}",200,min,max);
        TH1F *h_lamcsel=new TH1F("h_lamcsel","Lam_c",200,min,max);
        TH1F *h_lamcsig=new TH1F("h_lamcsig","Lam_c",200,min,max);
        h_lamcsig->SetLineColor(2);
        configHisto(h_lamcsel);
        
        // ****** Lambda
        TH1F *h_lam=new TH1F("h_lam","#Lambda #rightarrow p#pi^{-}",200,min,max);
        TH1F *h_lamsel=new TH1F("h_lamsel","Lam",200,min,max);
        TH1F *h_lamsig=new TH1F("h_lamsig","Lam",200,min,max);
        h_lamsig->SetLineColor(2);
        configHisto(h_lamsel);
        
        TNtuple *ntp=new TNtuple("ntp","ntp","m:tht:pt:p:d1tht:d1pt:d1p:d2tht:d2pt:d2p:d3tht:d3pt:d3p");
        
        h_jpsi->SetXTitle("m(l^{+}l^{-}) [GeV/c^{2}]");

        h_d0->SetXTitle("m(K^{-}#pi^{+}) [GeV/c^{2}]");
        h_d02->SetXTitle("m(K^{-}#pi^{+}#pi^{0}) [GeV/c^{2}]");
        h_d03->SetXTitle("m(K^{-}#pi^{+}#pi^{+}#pi^{-}) [GeV/c^{2}]");

        h_dpm->SetXTitle("m(K^{-}#pi^{+}#pi^{+}) [GeV/c^{2}]");

        h_ds->SetXTitle("m(K^{+}K^{-}#pi^{+}) [GeV/c^{2}]");

        h_phi->SetXTitle("m(K^{+}K^{-}) [GeV/c^{2}]");
        h_lamc->SetXTitle("m(pK^{-}#pi^{+}) [GeV/c^{2}]");
        h_lam->SetXTitle("m(p#pi^{-}) [GeV/c^{2}]");
        
        // *** prepare the tree to loop through

        Int_t nTrk,Id[fMAX],DF[fMAX],DL[fMAX];
        Double_t px[fMAX],py[fMAX],pz[fMAX],E[fMAX],m[fMAX];
    TClonesArray* DpmPart=0;
        
        if (dpm)
        {
        DpmPart=new TClonesArray("TParticle",100);
                t->SetBranchAddress("Npart",&nTrk);
                t->SetBranchAddress("Particles", &DpmPart);
        }
        else
        {
                t->SetBranchAddress("nTrk",&nTrk);
                t->SetBranchAddress("Id",Id);
                t->SetBranchAddress("DF",DF);
                t->SetBranchAddress("DL",DL);
                t->SetBranchAddress("px",px);
                t->SetBranchAddress("py",py);
                t->SetBranchAddress("pz",pz);
                t->SetBranchAddress("E",E);
                t->SetBranchAddress("m",m);     
        }
        
        
        int combCnt=0;  // how many combinatione are within our mass criterion
        int evCnt=0;    // how many _events_ (not combinations!) are within our mass criterion
        int evMultCut=0;
        
        // *** the event loop
        
        int mixevts=1;
        int currmix=0;
        
        int cntInLam=0;
        
        // how many events are tagged by individual algo?
        bool channelselect[13];
        double channelcount[13];
        TString channelname[13]={"Phi(KK)","Lam(ppi)","J/psi(ll)","D+-(K2pi)","D+-(K2pipi0)","D+-(KSpipi0)","D+-(KS3pi)","D0(Kpi)","D0(Kpipi0)","D0(K3pi)","Ds(KKpi)","Ds(KKpipi0)","Lamc(pKpi)"};
        
        for (i=0;i<13;++i) channelcount[i]=0.;
        
        if (fRand.Rndm()<P_mix) mixevts=2;
        
        for (i=0;i<nev;++i)
        {       
                t->GetEntry(i);
                
                if (i%5000==0) cout <<"ev "<<i<<endl;
                
                // *** read in from file the mc candidate
                if (0==currmix)
                {
                        mclist.clear();   // *** contains all MC cands (also non-final-states)
                        mcFSlist.clear(); // *** contains only final state MC cands 
                }
                
                bool good = true;
                
                if (dpm)
                {
                        int npi = 0;
                        int nstab = 0;
                        
                        for (j=0;j<nTrk;++j)
                        {
                                TParticle *pt=(TParticle*)DpmPart->At(j);
                                Id[j] = pt->GetPdgCode();
                                int pdg = abs(Id[j]);
                                px[j] = pt->Px();
                                py[j] = pt->Py();
                                pz[j] = pt->Pz();
                                E[j]  = pt->Energy();
                                m[j]  = pt->GetMass();
                                DF[j] = 0;
                                DL[j] = -1;
                                if (pdg==211) npi++;
                                if (isFinalState(pdg)) nstab++; 
                        }
                        //if (npi!=4 || nstab!=npi ) good = false; 
                }
                
                if (!good) continue;
                
                for (j=0;j<nTrk;++j)
                {
                        int chrg = 0;
                        int pdg = abs(Id[j]);
                        TParticlePDG *p=pdg_db->GetParticle(Id[j]); 
                        if (p) chrg = p->Charge();
                        if (abs(chrg)>1) chrg/=3;
                                                
                        SimpleCand c(px[j], py[j], pz[j], E[j], Id[j], j, chrg);
                        c.SetDau(DF[j], DL[j]);
                        
                        for (k=0;k<j;++k) if (j>=DF[k] && j<=DL[k]) 
                        {
                                c.SetMotherIdx(k);
                                c.SetNSiblings(DL[k]-DF[k]+1);
                        }
                        
                        c.SetEvtId(i);
                        
                        mclist.push_back(c);
                        
                        // final state particle?
                        if (isFinalState(pdg)) mcFSlist.push_back(c); 
                }
                
                if (++currmix<mixevts) continue;
                
                currmix=0;
                
                h_mult->Fill(mcFSlist.size());
                
                // *** if mass cut is negative (= number of sigmas), compute absolute window size
                //if (mwin<0) mwidth=-mwin*(dp*mwsig+mwoff);
                
                // *** prepare the toy reco: smear momentum and apply track efficiency (both flat at the moment)
                makeRecoCands(mcFSlist, recolist, trkeff/100., dp );
                
                // *** the boost of the system
                TVector3 boost = fIni.BoostVector();
                boostList(recolist, recocmslist, boost);
                
                // *** Make Event Selection
                EventShape evsh(recolist, fIni);
/*              if (evcut)
                {
                        if (evshape.ChrgPtSumLab()<1.5) continue;
                }*/
                
                h_multrec->Fill(recolist.size());
                //cout <<mcFSlist.size()<<"  "<<recolist.size()<<endl;
                
                // *** prepare all pid list according to the pidtable probabilities
                makePidSelection(recolist);
                
                int npi = piplus.size()+piminus.size();
                int nk = kplus.size()+kminus.size();
                int npr = pplus.size()+pminus.size();
                int nlep = eplus.size()+eminus.size()+muplus.size()+muminus.size();
                
                int nk10 = 0;  // # kaons with p>1.0
                for (j=0;j<kplus.size();++j) if (kplus[j].P4().P()>1.0) nk10++;
                for (j=0;j<kminus.size();++j) if (kminus[j].P4().P()>1.0) nk10++;
                
                int npr10 = 0;  // # protons with p>1.0
                for (j=0;j<pplus.size();++j) if (pplus[j].P4().P()>1.0) npr10++;
                for (j=0;j<pminus.size();++j) if (pminus[j].P4().P()>1.0) npr10++;
                
                h_mult_k->Fill(nk);
                h_mult_lep->Fill(nlep);
                h_mult_klep->Fill(nk, nlep);
                h_sum_klep->Fill(nk + nlep);
                
                                                                        
                // *** the lists we need for combinatorics
                CandList comb1, comb2, comb3;

                CandList LPhi, LPi0, LKs;                               // phi -> KK , pi0 -> gg (mass cut), KS -> pi+ pi- (mass cut)
                CandList LJpsiE, LJpsiMu;                       // J/psi -> ee/mu mu
                CandList LDp, LDm, LDp2, LDm2, LDp3, LDm3, LDp4, LDm4;  // D+- -> K2pi, K2pipi0, KS pi pi0, KS pi pi pi
                CandList LDsp, LDsm, LDsp2, LDsm2;              // Ds -> KKpi, KKpipi0                          
                CandList LD0, LD02, LD0b, LD02b, LD03, LD03b;   // D0 -> Kpi, Kpipi0, K3pi
                CandList LLamcp, LLamcm;                                // Lc -> pKpi
                CandList LLam, LLamb;                                   // Lam -> ppi
                
                
                int ncomb1=0, ncomb2=0;
                bool selected=false;
                
                for (j=0;j<13;++j) channelselect[j]=false;
                                
                // Apply cut on kaon/lepton multiplicity
/*              if ( multcut>0 && (N_K+N_lep<multcut) )
                {
                        mixevts=1;
                        if (fRand.Rndm()<P_mix) mixevts=2;
                        continue;
                }*/
                
                evMultCut++;

                bool accepted = false;
                
                //
                //  ********* pi0 / KS
                //
                // *** do combinatorics (pi0 -> gg) 
                combine(gam, gam, comb1, fPi0Pdg);
                fillHisto(comb1, h_pi0, 0, fPi0Mass, fPi0win, h_pi0sel, 0, h_pi0sig ); 
                selectMass(comb1, LPi0, fPi0Mass, fPi0win);
                
                // *** do combinatorics (pi0 -> gg) 
                combine(piplus, piminus, comb1, fKsPdg);
                fillHisto(comb1, h_ks, 0, fKsMass, fKswin, h_kssel, 0, h_kssig ); 
                selectMass(comb1, LKs, fKsMass, fKswin);
                

                // **** cache some eventshape vars
                int npart      = evsh.NParticles();
                int nchrg      = evsh.NCharged();
                int nneut      = evsh.NNeutral();
                int npi0       = LPi0.size();                                   // #pi0 candidates
                int nks0       = LKs.size();                                    // #KS0 candidates
                int npi005     = 0.0;
                int     npi010     = 0.0;
                int nks005     = 0.0;
                int nks010     = 0.0;
                int np05       = evsh.MultPminCms(0.5);
                int np10       = evsh.MultPminCms(1.0);

                double pmaxl   = evsh.PmaxLab();
                double pmax    = evsh.PmaxCms();
                double pmin    = evsh.PminCms();
                double ptmax   = evsh.Ptmax();
                
                double sumpc   = evsh.ChrgPSumCms();
                double sumpc05 = evsh.SumChrgPminCms(0.5);
                double sumpt   = evsh.PtSumLab();
                double sumptc  = evsh.ChrgPtSumLab();
                double sumptl  = sumpt;
                double sumen   = evsh.NeutESumCms();
                
                double fw1     = evsh.FoxWolfMomR(1);
                double fw2     = evsh.FoxWolfMomR(2);
                double fw3     = evsh.FoxWolfMomR(3);
                double fw4     = evsh.FoxWolfMomR(4);
                double fw5     = evsh.FoxWolfMomR(5);
                
                // count pi0 with threshold
                for (j=0; j<LPi0.size(); ++j) 
                {
                        if (LPi0[j].P4().P()>0.5) npi005++;
                        if (LPi0[j].P4().P()>1.0) npi010++;
                }
                // count KS with threshold
                for (j=0; j<LKs.size(); ++j) 
                {
                        if (LKs[j].P4().P()>0.5) nks005++;
                        if (LKs[j].P4().P()>1.0) nks010++;
                }
                
                h_fw1->Fill(fw2);
                h_pmax->Fill(pmax);

                
                // **** find out, whether event was accepted
                
                // ********* PHI -> K+ K-
                if ( !evcut 
                     || (sqrts==2.4  &&   nk>1 && npart<7 && pmax<0.648   )
                         || (sqrts==3.77 &&   nk>1 && fw2>0.5456 )
                         || (sqrts==4.28 &&   nk>1 && 1 )
                         || (sqrts==5.0  &&   nk>1 && 1 )
                         || (sqrts==5.5  &&   nk>1 && 1 )
                        )
                {
                  // *** do combinatorics (phi -> K K) 
                  combine (kplus, kminus, LPhi, fPhiPdg); 
                  
                  accepted |= eventAccepted(LPhi,fPhiMass, fPhiwin);
                }


                // ********* LAMBDA -> p pi-
                if ( !evcut 
                     || (sqrts==2.4  &&   npart>3 && npr>0 && npi>0 && fw1>0.132 && fw4>0.2024 )
                         || (sqrts==3.77 &&   npart>3 && pmax>1.056 &&     fw1>0.0   && fw2>0.6248 && fw5>0.22  )
                         || (sqrts==4.28 &&   npart>3 && 1 )
                         || (sqrts==5.0  &&   npart>3 && 1 )
                         || (sqrts==5.5  &&   npart>3 && 1 )
                        )
                {
                  // *** do combinatorics (Lam -> p pi) 
                  combine (pplus, piminus, LLam, fLamPdg); 
                  combine (pminus, piplus, LLamb, -fLamPdg);            
                  
                  accepted |= eventAccepted(LLam, fLamMass, fLamwin);
                  accepted |= eventAccepted(LLamb, fLamMass, fLamwin);
                }


                // ********* JPSI -> e+e-/mu+mu-
                if ( !evcut 
                         || (sqrts==3.77 &&  np10>1 && npi<4 && pmin<0.3 && sumpc>3.296 )
                         || (sqrts==4.28 &&  1 )
                         || (sqrts==5.0  &&  1 )
                         || (sqrts==5.5  &&  1 )
                        )
                {
                  // *** do combinatorics (J/psi -> ll) 
                  combine(eminus, eplus, LJpsiE, fJpsiPdg);
                  combine(muminus, muplus, LJpsiMu, fJpsiPdg);
                  
                  accepted |= eventAccepted(LJpsiE,  fJMass, fJwin);
                  accepted |= eventAccepted(LJpsiMu,  fJMass, fJwin);
                }
                
                
                // ********* D+ -> K- pi+ pi+
                if ( !evcut 
                         || (sqrts==3.77 &&   nk>0 && npart<12 && pmax<0.936 && sumpc>1.76 )
                         || (sqrts==4.28 &&   nk>0 && 1 )
                         || (sqrts==5.0  &&   nk>0 && 1 )
                         || (sqrts==5.5  &&   nk>0 && 1 )
                        )
                {
                  // *** do combinatorics (D+- -> K pi pi) 
                  combine (piplus,  piplus,  kminus, LDp, fDPdg); 
                  combine (piminus, piminus, kplus,  LDm, -fDPdg); 
                  
                  accepted |= eventAccepted(LDp,fDMass, fDwin);
                  accepted |= eventAccepted(LDm,fDMass, fDwin);
                }
                
                // ********* D+ -> K- pi+ pi+ pi0
                if ( !evcut 
                         || (sqrts==3.77 && nk>0 && nneut>1 && pmax<0.888 && pmaxl<2.76 && sumptc>0.84  )
                         || (sqrts==4.28 && nk>0 )
                         || (sqrts==5.0  && nk>0 )
                         || (sqrts==5.5  && nk>0 )
                        )
                {
                  // *** do combinatorics (D+- -> K pi pi pi0) 
                  combine (piplus,  piplus,  kminus, LPi0,  LDp2, fDPdg); 
                  combine (piminus, piminus, kplus,  LPi0, LDm2, -fDPdg); 
                  
                  accepted |= eventAccepted(LDp2,fDMass, fDwin);
                  accepted |= eventAccepted(LDm2,fDMass, fDwin);
                }
                
                //********* D+ -> KS pi+ pi0
                if ( !evcut 
                         || (sqrts==3.77 && nks0>0 && npi0>0 && pmax<0.91 && sumen<1.62 )
                         || (sqrts==4.28 && 1 )
                         || (sqrts==5.0  && 1 )
                         || (sqrts==5.5  && 1 )
                        )
                {
                  // *** do combinatorics (D+- -> KS pi pi0) 
                  combine (piplus,  LPi0,  LKs, LDp3, fDPdg); 
                  combine (piminus, LPi0,  LKs, LDm3, -fDPdg); 
                  
                  accepted |= eventAccepted(LDp3,fDMass, fDwin);
                  accepted |= eventAccepted(LDm3,fDMass, fDwin);
                }
                
                //********* D+ -> KS pi+ pi+ pi-
                if ( !evcut 
                         || (sqrts==3.77 && nchrg>5 && np05<4 && pmax<0.888 && pmaxl<2.64 && sumen<1.02)
                         || (sqrts==4.28 && 1 )
                         || (sqrts==5.0  && 1 )
                         || (sqrts==5.5  && 1 )
                        )
                {
                  // *** do combinatorics (D+- -> KS pi pi pi) 
                  combine (piplus,  piplus,  piminus, LKs, LDp4, fDPdg); 
                  combine (piminus, piminus, piplus,  LKs, LDm4, -fDPdg); 
                  
                  accepted |= eventAccepted(LDp4,fDMass, fDwin);
                  accepted |= eventAccepted(LDm4,fDMass, fDwin);
                }
                

                // ********* D0 -> K- pi+
                if ( !evcut 
                         || (sqrts==3.77 &&   nk>0 && sumpc05>3.18 && pmax>0.768 && pmax<1.056 && fw3<0.1848 )
                         || (sqrts==4.28 &&   nk>0 && 1 )
                         || (sqrts==5.0  &&   nk>0 && 1 )
                         || (sqrts==5.5  &&   nk>0 && 1 )
                        )
                {
                  // *** do combinatorics (D0 -> K pi) 
                  combine (kplus, piminus, LD0, -fD0Pdg); 
                  combine (kminus, piplus, LD0b, fD0Pdg); 

                  accepted |= eventAccepted(LD0, fD0Mass, fD0win);
                  accepted |= eventAccepted(LD0b, fD0Mass, fD0win);
                }       

                // ********* D0 -> K- pi+ pi0
                if ( !evcut 
                         || (sqrts==3.77 && nk>0 && nneut>1 && npi<6 && pmax<0.948 && sumpt>1.56 && sumen<1.92)
                         || (sqrts==4.28 && nk>0 )
                         || (sqrts==5.0  && nk>0 )
                         || (sqrts==5.5  && nk>0 )
                        )
                {
                  // *** do combinatorics (D0 -> K pi pi0) 
                  combine (kplus, piminus, LPi0, LD02, -fD0Pdg); 
                  combine (kminus, piplus, LPi0, LD02b, fD0Pdg); 
                  
                  accepted |= eventAccepted(LD02, fD0Mass, fD0win);
                  accepted |= eventAccepted(LD02b, fD0Mass, fD0win);
                }       

                // ********* D0 -> K- pi+ pi+ pi-
                if ( !evcut 
                         || (sqrts==3.77 && nk>0 && nchrg>5 && npi>2 && pmax<0.96 && sumen<1.32)
                         || (sqrts==4.28 && nk>0 )
                         || (sqrts==5.0  && nk>0 )
                         || (sqrts==5.5  && nk>0 )
                        )
                {
                  // *** do combinatorics (D0 -> K pi pi pi) 
                  combine (kplus, piminus, piminus, piplus, LD03, -fD0Pdg); 
                  combine (kminus, piplus, piminus, piplus, LD03b, fD0Pdg); 
                        
                  accepted |= eventAccepted(LD03, fD0Mass, fD0win);
                  accepted |= eventAccepted(LD03b, fD0Mass, fD0win);
                }       


                // ********* DS -> K+ K- pi+
                if ( !evcut 
                         || (sqrts==4.28 &&   nk>1 && 1 )
                         || (sqrts==5.0  &&   nk>1 && 1 )
                         || (sqrts==5.5  &&   nk>1 && 1 )
                        )
                {
                  // *** do combinatorics (Ds+- -> K K pi) 
                  combine (kplus, kminus, piplus, LDsp, fDsPdg); 
                  combine (kplus, kminus, piminus, LDsm, -fDsPdg); 

                  accepted |= eventAccepted(LDsp, fDsMass, fDswin);
                  accepted |= eventAccepted(LDsm, fDsMass, fDswin);
                }

                // ********* DS -> K+ K- pi+ pi0
                if ( !evcut 
                         || (sqrts==4.28 &&   nk>1 )
                         || (sqrts==5.0  &&   nk>1 )
                         || (sqrts==5.5  &&   nk>1 )
                        )
                {
                  // *** do combinatorics (Ds+- -> K K pi pi0) 
                  combine (kplus, kminus, piplus,  LPi0, LDsp2, fDsPdg); 
                  combine (kplus, kminus, piminus, LPi0, LDsm2, -fDsPdg); 

                  accepted |= eventAccepted(LDsp2, fDsMass, fDswin);
                  accepted |= eventAccepted(LDsm2, fDsMass, fDswin);
                }


                // ********* LAMBDA_C -> p K- pi+
                if ( !evcut 
                         || (sqrts==5.0  &&  1 )
                         || (sqrts==5.5  &&  1 )
                        )
                {
                  // *** do combinatorics (Lamc) 
                  combine (pplus, kminus, piplus, LLamcp, fLamcPdg); 
                  combine (pminus, kplus, piminus, LLamcm, -fLamcPdg);          
                
                  accepted |= eventAccepted(LLamcp, fLamcMass, fLamcwin);
                  accepted |= eventAccepted(LLamcm, fLamcMass, fLamcwin);
                }

                
                // ********* NOW do counting and drawing, ONLY for events which are accepted by any selector                    
            // {"Phi(KK)","Lam(ppi)","J/psi(ll)","D+-(K2pi)","D+-(K2pipi0)","D+-(KSpipi0)","D+-(KS3pi)","D0(Kpi)","D0(Kpipi0)","D0(K3pi)","Ds(KKpi)","Ds(KKpipi0)","Lamc(pKpi)"};
        //     0          1         2           3              4             5             6           7           8          9          10            11           12
                if (!simcut || accepted)
                {
                        // **** Phi
                        if (sqrts>2.2)
                        {
                          ncomb1 = fillHisto(LPhi, h_phi, 0, fPhiMass, fPhiwin, h_phisel, 0, h_phisig  ); 
                          combCnt += ncomb1;
                          if (ncomb1>0) {selected=true; channelselect[0]=true;}
                        }

                        // **** Lam
                        if (sqrts>2.2)
                        {
                          ncomb1 = fillHisto(LLam, h_lam, 0, fLamMass, fLamwin, h_lamsel, 0, h_lamsig ); 
                          ncomb2 = fillHisto(LLamb, h_lam, 0, fLamMass, fLamwin, h_lamsel, 0, h_lamsig ); 
                          combCnt += ncomb1 + ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[1]=true;}
                        }
                        
                        // **** J/psi
                        if (sqrts>3.1)
                        {
                          ncomb1 = fillHisto(LJpsiE, h_jpsi, 0, fJMass, fJwin, h_jpsisel, 0, h_jpsisig ); 
                          ncomb2 = fillHisto(LJpsiMu, h_jpsi, 0, fJMass, fJwin, h_jpsisel, 0, h_jpsisig ); 
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[2]=true;}
                        }
                        
                        // **** Dp
                        if (sqrts>3.75)
                        {
                          ncomb1 = fillHisto(LDp, h_dpm, 0, fDMass, fDwin, h_dpmsel, 0, h_dpmsig ); 
                          ncomb2 = fillHisto(LDm, h_dpm, 0, fDMass, fDwin, h_dpmsel, 0, h_dpmsig );     
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[3]=true;}
                          
                          ncomb1 = fillHisto(LDp2, h_dpm2, 0, fDMass, fDwin, h_dpm2sel, 0, h_dpm2sig ); 
                          ncomb2 = fillHisto(LDm2, h_dpm2, 0, fDMass, fDwin, h_dpm2sel, 0, h_dpm2sig );         
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[4]=true;}
                          
                          ncomb1 = fillHisto(LDp3, h_dpm3, 0, fDMass, fDwin, h_dpm3sel, 0, h_dpm3sig ); 
                          ncomb2 = fillHisto(LDm3, h_dpm3, 0, fDMass, fDwin, h_dpm3sel, 0, h_dpm3sig );         
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[5]=true;}
                          
                          ncomb1 = fillHisto(LDp4, h_dpm4, 0, fDMass, fDwin, h_dpm4sel, 0, h_dpm4sig ); 
                          ncomb2 = fillHisto(LDm4, h_dpm4, 0, fDMass, fDwin, h_dpm4sel, 0, h_dpm4sig );         
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[6]=true;}
                          
                        }

                        // **** D0 
                        if (sqrts>3.75)
                        {
                          ncomb1 = fillHisto(LD0, h_d0, 0, fD0Mass, fD0win, h_d0sel, 0, h_d0sig ); 
                          ncomb2 = fillHisto(LD0b, h_d0, 0, fD0Mass, fD0win, h_d0sel, 0, h_d0sig );     
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[7]=true;}
                          
                          ncomb1 = fillHisto(LD02, h_d02, 0, fD0Mass, fD0win, h_d02sel, 0, h_d02sig ); 
                          ncomb2 = fillHisto(LD02b, h_d02, 0, fD0Mass, fD0win, h_d02sel, 0, h_d02sig );         
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[8]=true;}
                          
                          ncomb1 = fillHisto(LD03, h_d03, 0, fD0Mass, fD0win, h_d03sel, 0, h_d03sig ); 
                          ncomb2 = fillHisto(LD03b, h_d03, 0, fD0Mass, fD0win, h_d03sel, 0, h_d03sig );         
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[9]=true;}
                        }

                        // **** Ds
                        if (sqrts>4.)
                        {
                          ncomb1 = fillHisto(LDsp, h_ds, 0, fDsMass, fDswin, h_dssel, 0, h_dssig ); 
                          ncomb2 = fillHisto(LDsm, h_ds, 0, fDsMass, fDswin, h_dssel, 0, h_dssig ); 
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[10]=true;}
                          
                          ncomb1 = fillHisto(LDsp2, h_ds2, 0, fDsMass, fDswin, h_ds2sel, 0, h_ds2sig ); 
                          ncomb2 = fillHisto(LDsm2, h_ds2, 0, fDsMass, fDswin, h_ds2sel, 0, h_ds2sig ); 
                          combCnt += ncomb1 +ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[11]=true;}
                        }
                
                        // **** Lamc
                        if (sqrts>4.57)
                        {
                          ncomb1 = fillHisto(LLamcp, h_lamc, 0, fLamcMass, fLamcwin, h_lamcsel, 0, h_lamcsig ); 
                          ncomb2 = fillHisto(LLamcm, h_lamc, 0, fLamcMass, fLamcwin, h_lamcsel, 0, h_lamcsig ); 
                          combCnt += ncomb1 + ncomb2;
                          if (ncomb1>0 || ncomb2>0) {selected=true; channelselect[12]=true;}
                        }
                }
            // ****************
            
                if (selected) evCnt+=mixevts; // *** if any combination is in window count event
                
                for (j=0;j<13;++j) if (channelselect[j]) channelcount[j]+=mixevts;
                
                //hmix->Fill(mixevts);
                // new round of mixing
                mixevts=1;
                if (fRand.Rndm()<P_mix) mixevts=2;
                
                //if (i<3) printAllLists();
        }
        //cout <<"cntInLam "<<cntInLam<<endl;
        cout<< "#Combinations:"<<combCnt<<"  #Events:"<<evCnt<<"  ";
        cout <<"Efficiency   (ev):"<<(double)evCnt/(double)nev<<endl;
        cout <<"Mult Cut Eff (ev):"<<(double)evMultCut/(double)nev<<endl;
        
        for (i=0;i<13;++i) 
        {
                channelcount[i]/=nev;
                cout <<channelname[i]<<":"<<channelcount[i]<<endl;
        }
        cout <<endl;
        
        TLatex latex;
        latex.SetTextSize(0.055);
        char tmp[200];
        double offset = 0.1;
        // *** plot all stuff and write out info
        c2->cd(1); h_pi0->Draw(); h_pi0sel->Draw("same"); h_pi0sig->Draw("same");
        c2->cd(2); h_ks->Draw(); h_kssel->Draw("same"); h_kssig->Draw("same");

        c2->cd(3); h_phi->Draw(); h_phisel->Draw("same");   h_phisig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[0]*100.);
        latex.DrawLatex(2.5,h_phi->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(4); h_lam->Draw(); h_lamsel->Draw("same"); h_lamsig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[1]*100.);
        latex.DrawLatex(2.5,h_lam->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(5); h_jpsi->Draw(); h_jpsisel->Draw("same"); h_jpsisig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[2]*100.);
        latex.DrawLatex(2.5,h_jpsi->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(6); h_dpm->Draw();  h_dpmsel->Draw("same");  h_dpmsig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[3]*100.);
        latex.DrawLatex(2.5,h_dpm->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(7); h_dpm2->Draw();  h_dpm2sel->Draw("same");  h_dpm2sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[4]*100.);
        latex.DrawLatex(2.5,h_dpm2->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(8); h_dpm3->Draw();  h_dpm3sel->Draw("same");  h_dpm3sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[5]*100.);
        latex.DrawLatex(2.5,h_dpm3->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(9); h_dpm4->Draw();  h_dpm4sel->Draw("same");  h_dpm4sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[6]*100.);
        latex.DrawLatex(2.5,h_dpm4->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(10); h_d0->Draw();   h_d0sel->Draw("same");   h_d0sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[7]*100.);
        latex.DrawLatex(2.5,h_d0->GetMaximum()*0.95+offset,tmp);

        c2->cd(11); h_d02->Draw();   h_d02sel->Draw("same");   h_d02sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[8]*100.);
        latex.DrawLatex(2.5,h_d02->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(12); h_d03->Draw();   h_d03sel->Draw("same");   h_d03sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[9]*100.);
        latex.DrawLatex(2.5,h_d03->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(13); h_ds->Draw();   h_dssel->Draw("same");   h_dssig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[10]*100.);
        latex.DrawLatex(2.5,h_ds->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(14); h_ds2->Draw();   h_ds2sel->Draw("same");   h_ds2sig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[11]*100.);
        latex.DrawLatex(2.5,h_ds2->GetMaximum()*0.95+offset,tmp);
        
        c2->cd(15); h_lamc->Draw(); h_lamcsel->Draw("same"); h_lamcsig->Draw("same");
        sprintf(tmp,"#epsilon = %4.1f%%",channelcount[12]*100.);
        latex.DrawLatex(2.5,h_lamc->GetMaximum()*0.95+offset,tmp);
        
        //sprintf(tmp,"#epsilon = %4.1f%%",channelcount[6]*100.);
        //latex.DrawLatex(2.5,h_pi0->GetMaximum()*0.95,tmp);
        
        c2->cd();
        TPad *peff=new TPad("peff","peff",0.46,0.665,0.54,0.70);
        peff->Draw();
        peff->cd();
        sprintf(tmp,"#epsilon_{tot} = %4.1f%%",double(evCnt)/double(nev)*100.);
        latex.SetTextSize(0.7);
        latex.SetTextColor(2);
        latex.DrawLatex(0.15,0.3,tmp);
        c2->cd();
        
        h_mult_klep->Scale(1.0/nev);
        
/*      c3->cd(1); h_sum_klep->Draw();   //h_dpmmm->Draw();
        c3->cd(2); h_mult_klep->Draw("colz");//h_d0mm->Draw();
        c3->cd(3); h_mult_k->Draw();   // h_dpmmvsmm->Draw("colz");
        c3->cd(4); h_mult_lep->Draw(); // h_d0mvsmm->Draw("colz");
        c3->cd();
        */
//      c4->cd();
//      h_pmax->Draw();
        //h_lamc->Draw(); h_lamcsel->Draw("same"); h_lamcsig->Draw("same");     
        //h_dpm->Draw();  h_dpmsel->Draw("same");  h_dpmsig->Draw("same");
        
  return 0;
}

Variable Documentation

CandList eminus

Definition at line 124 of file softtrigger_toy12.C.

CandList eplus

Definition at line 124 of file softtrigger_toy12.C.

TF1* f_res =0

Definition at line 110 of file softtrigger_toy12.C.

double fD0Mass = 1.8645

Definition at line 56 of file softtrigger_toy12.C.

int fD0Pdg = 421

Definition at line 55 of file softtrigger_toy12.C.

double fD0sig = 1.04

Definition at line 57 of file softtrigger_toy12.C.

double fD0win =0.25

Definition at line 58 of file softtrigger_toy12.C.

const double fdE = 0.05

Definition at line 37 of file softtrigger_toy12.C.

double fDMass = 1.8693

Definition at line 51 of file softtrigger_toy12.C.

int fDPdg = 411

Definition at line 50 of file softtrigger_toy12.C.

double fDsig = 0.88

Definition at line 52 of file softtrigger_toy12.C.

double fDsMass = 1.9686

Definition at line 73 of file softtrigger_toy12.C.

int fDsPdg = 431

Definition at line 72 of file softtrigger_toy12.C.

double fDssig = 0.75

Definition at line 74 of file softtrigger_toy12.C.

double fDswin =0.2

Definition at line 75 of file softtrigger_toy12.C.

double fDwin =0.25

Definition at line 53 of file softtrigger_toy12.C.

double fEtacMass = 2.9798

Definition at line 83 of file softtrigger_toy12.C.

int fEtacPdg = 441

Definition at line 82 of file softtrigger_toy12.C.

double fEtacsig = 0.15

Definition at line 84 of file softtrigger_toy12.C.

double fEtacwin =0.15

Definition at line 85 of file softtrigger_toy12.C.

TLorentzVector fIni

Definition at line 27 of file softtrigger_toy12.C.

double fJ3pioff = 0.027

Definition at line 64 of file softtrigger_toy12.C.

double fJMass = 3.09687

Definition at line 61 of file softtrigger_toy12.C.

int fJpsiPdg = 443

Definition at line 60 of file softtrigger_toy12.C.

double fJsig = 2.24

Definition at line 62 of file softtrigger_toy12.C.

double fJsig3pi = 1.5

Definition at line 63 of file softtrigger_toy12.C.

double fJwin =0.7

Definition at line 65 of file softtrigger_toy12.C.

double fKsMass = 0.497614

Definition at line 47 of file softtrigger_toy12.C.

int fKsPdg = 310

Definition at line 46 of file softtrigger_toy12.C.

double fKswin = 0.025

Definition at line 48 of file softtrigger_toy12.C.

double fLamcMass = 2.2849

Definition at line 88 of file softtrigger_toy12.C.

int fLamcPdg = 4122

Definition at line 87 of file softtrigger_toy12.C.

double fLamcsig = 0.15

Definition at line 89 of file softtrigger_toy12.C.

double fLamcwin =0.2

Definition at line 90 of file softtrigger_toy12.C.

double fLamMass = 1.115684

Definition at line 68 of file softtrigger_toy12.C.

int fLamPdg = 3122

Definition at line 67 of file softtrigger_toy12.C.

double fLamsig = 0.08

Definition at line 69 of file softtrigger_toy12.C.

double fLamwin =0.1

Definition at line 70 of file softtrigger_toy12.C.

const int fMAX = 100

Definition at line 30 of file softtrigger_toy12.C.

const double fNdphi = 0.003

Definition at line 39 of file softtrigger_toy12.C.

const double fNdtht = 0.003

Definition at line 38 of file softtrigger_toy12.C.

const double fNeutEff = 1.0

Definition at line 36 of file softtrigger_toy12.C.

double fPhiMass = 1.0195

Definition at line 78 of file softtrigger_toy12.C.

int fPhiPdg = 333

Definition at line 77 of file softtrigger_toy12.C.

double fPhisig = 0.15

Definition at line 79 of file softtrigger_toy12.C.

double fPhiwin =0.2

Definition at line 80 of file softtrigger_toy12.C.

double fPi0Mass = 0.13497

Definition at line 43 of file softtrigger_toy12.C.

int fPi0Pdg = 111

Definition at line 42 of file softtrigger_toy12.C.

double fPi0win = 0.020

Definition at line 44 of file softtrigger_toy12.C.

TRandom3 fRand

Definition at line 26 of file softtrigger_toy12.C.

CandList gam

Definition at line 124 of file softtrigger_toy12.C.

double high = 5.0

Definition at line 33 of file softtrigger_toy12.C.

CandList kminus

Definition at line 124 of file softtrigger_toy12.C.

CandList kplus

Definition at line 124 of file softtrigger_toy12.C.

double low = 0.0

Definition at line 32 of file softtrigger_toy12.C.

CandList mcFSlist

Definition at line 121 of file softtrigger_toy12.C.

CandList mclist

Definition at line 121 of file softtrigger_toy12.C.

double mE

Definition at line 108 of file softtrigger_toy12.C.

double mK

Definition at line 108 of file softtrigger_toy12.C.

double mMu

Definition at line 108 of file softtrigger_toy12.C.

double mP

Definition at line 108 of file softtrigger_toy12.C.

double mPi

Definition at line 108 of file softtrigger_toy12.C.

double mPi0

Definition at line 108 of file softtrigger_toy12.C.

CandList muminus

Definition at line 124 of file softtrigger_toy12.C.

CandList muplus

Definition at line 124 of file softtrigger_toy12.C.

int nbins = 300

Definition at line 31 of file softtrigger_toy12.C.

TDatabasePDG* pdg_db

Definition at line 116 of file softtrigger_toy12.C.

std::map<int,int> pdgidx

Definition at line 114 of file softtrigger_toy12.C.

double pidreal[25]

Initial value:

= {  95,  6,  10,   7,  10,    
                         5,  86,  27,   6,   7,    
                         9,  29,  85,  18,  16,    
                         9,  21,  40,  67,  24,    
                         12, 12,  20,  16,  90 }

Definition at line 101 of file softtrigger_toy12.C.

double pidtab[25]

Initial value:

= {  80,   5,   5,   5,   5,    
                         5,  80,   5,   5,   5,    
                         5,   5,  80,   5,   5,    
                         5,   5,   5,  80,   5,    
                         5,   5,   5,   5,  80 }

Definition at line 94 of file softtrigger_toy12.C.

CandList piminus

Definition at line 124 of file softtrigger_toy12.C.

CandList piplus

Definition at line 124 of file softtrigger_toy12.C.

CandList pminus

Definition at line 124 of file softtrigger_toy12.C.

CandList pplus

Definition at line 124 of file softtrigger_toy12.C.

CandList recocmslist

Definition at line 121 of file softtrigger_toy12.C.

CandList recolist

Definition at line 121 of file softtrigger_toy12.C.