#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"
#include "TVector.h"
#include "TMatrixD.h"
#include "TSystem.h"
#include "TMatrixDEigen.h"
#include "TMath.h"
#include "EventShape.h"

Typedefs
typedef std::vector< SimpleCand >	CandList

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

double	eps (TVector3 v1, TVector3 v2)

double	compThrust (CandList &list, TVector3 &axis)

double	Legendre (int l, double x)

void	compFoxWolfMoms (CandList &list)

double	compSphericity (CandList &list)

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)

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)

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)

int	softtrigger_kin5 (TString fsig, int mode=0, double sqrts=3.77, int nev=1000, 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 = 0.9

const double	fdE = 0.03

const double	fNdtht = 0.003

const double	fNdphi = 0.003

const double	fNeutThrsh = 0.02

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.2

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	_sph

double	_pla

double	_apl

double	_thr

TVector3	_sphAx

TVector3	_secAx

TVector3	_thrAx

const int	_nfwmom = 5

double	_foxwolfR [_nfwmom]

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 138 of file softtrigger_kin5.C.

Function Documentation

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

Definition at line 147 of file softtrigger_kin5.C.

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

Referenced by softtrigger_kin5().

 {
         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 483 of file softtrigger_kin5.C.

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

Referenced by PndTrackCombiner::combine(), and softtrigger_kin5().

 {
     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());
             }
 
             out.push_back(c);
         }
     }
     return out.size();
 }

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

Definition at line 526 of file softtrigger_kin5.C.

References c, m2(), mclist, SimpleCand::SetCharge(), 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());
                 }
 
                 out.push_back(c);
             }
         }
     }
 
     return out.size();
 }

void compFoxWolfMoms ( CandList & list )

Definition at line 234 of file softtrigger_kin5.C.

References _foxwolfR, _nfwmom, cos(), i, Legendre(), n, p1, p2, and s.

 {
   int n=list.size();
   if( n==0 ) return;
 
   double s = 0.;
   double _sumarray[_nfwmom]={0,0,0,0};
   int i,j,l;
 
   for (i=0; i<n-1; ++i)
   {
      // this candidate's 3-momentum
           TVector3 p1(list[i].P4().Vect());
       double pmag1 = p1.Mag();
 
       // loop over other candidates, starting at the next one in the list
           
           for (j=i+1; j<n; ++j)
           {
                 // this candidate's 3-momentum
                 TVector3 p2(list[j].P4().Vect());
                 double pmag2 = p2.Mag();
 
                 // the cosine of the angle between the two candidates
                 double cosPhi =  cos ( p1.Angle(p2) );
 
                 // the contribution of this pair of track 
                 // (note the factor 2 : the pair enters the sum twice)
                 for( l=0; l<_nfwmom; l++ )
                   _sumarray[l] += 2 * pmag1 * pmag2 * Legendre( l, cosPhi );
                   
           }
       // contribution of this track
       for( l=0; l<_nfwmom; l++ )
                 _sumarray[l] += pmag1 * pmag1 * Legendre( l, 1. );
       
       // total energy
       s += list[i].P4().Energy();
    }
   
   // well ...
   if( s<=0. ) return;
   double ssq=s*s;
   // normalize Fox Wolfram Moments
   for( i=0; i<_nfwmom; i++)
     _foxwolfR[i] = _sumarray[i]/ssq ;
 }

double compSphericity ( CandList & list )

Definition at line 283 of file softtrigger_kin5.C.

References _apl, _pla, _sph, i, and v.

 {
   int N=list.size();
   
   double stot=0, sxx=0, sxy=0, sxz=0, syy=0, syz=0, szz=0;      
   int i;
   
   for (i=0;i<N;++i)
   {
         TVector3 v(list[i].P4().Vect());
         
         sxx += v.X()*v.X(); sxy += v.X()*v.Y(); sxz += v.X()*v.Z();
         syy += v.Y()*v.Y(); syz += v.Y()*v.Z();
         szz += v.Z()*v.Z();
         
         stot += v.Mag2();
   }
   
   TMatrixD sm(3,3);
   sm(0,0) = sxx/stot; sm(0,1) = sxy/stot; sm(0,2) = sxz/stot; 
   sm(1,0) = sxy/stot; sm(1,1) = syy/stot; sm(1,2) = syz/stot; 
   sm(2,0) = sxz/stot; sm(2,1) = syz/stot; sm(2,2) = szz/stot; 
   
   TMatrixDEigen ei(sm);
   TMatrixD eiv=ei.GetEigenValues();
   
   _sph = 1.5 * (eiv(1,1) + eiv(2,2));
   _apl = 1.5 * eiv(2,2);
   _pla = eiv(1,1) - eiv(2,2);
   
   return _sph;
 }

double compThrust	(	CandList &	list,
		TVector3 &	axis
	)

Definition at line 170 of file softtrigger_kin5.C.

References fabs(), and i.

 {
   TVector3 n0;
   int N=list.size();
   
   int i,j;
   double pmax=0;
   
   //find starting vector as maximum momentum vector
   for (i=0;i<N;++i)
   {
         if (list[i].P4().Vect().Mag()>pmax)
         {
           n0=list[i].P4().Vect();
           pmax=list[i].P4().Vect().Mag();
         }
   }
   //n0=list[0].P4().Vect();
   
   TVector3 nnew(0,0,0);
   // find thrust axis (10 iterations)
   for (i=0;i<10;++i)
   {
         for (j=0;j<N;++j)
         {
           //nnew+=eps(n0,list[j].P4().Vect())*list[j].P4().Vect();
           nnew+=n0*list[j].P4().Vect()>0 ? list[j].P4().Vect() : -list[j].P4().Vect();
         }
         n0=nnew.Unit();
   }
   
   axis = n0;
   double thr=0, sum=0;
   for (i=0;i<N;++i)
   {
         thr+=fabs(axis.Dot(list[i].P4().Vect()));
         sum+=list[i].P4().Vect().Mag();
   }
   return thr/sum;
 }

double dponline ( TLorentzVector & v )

Definition at line 326 of file softtrigger_kin5.C.

References f_res.

Referenced by smearMom().

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

double eps	(	TVector3	v1,
		TVector3	v2
	)

Definition at line 163 of file softtrigger_kin5.C.

References fabs().

Referenced by createdirc(), createdirc_prism(), createdircPix(), createRootGeometry_DIRC(), createRootGeometry_DIRC_fsEVdroplens_MCPs(), createRootGeometry_DIRC_sepEV_06_2013(), createRootGeometry_DIRC_sepEV_MCPs(), createRootGeometry_DIRC_updated_06_2013(), findtrack(), PndSttSingleStraw::Polya(), and PndFtsSingleStraw::Polya().

 {
   double prod = v1.Dot(v2);
   return prod/fabs(prod);
 }

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 457 of file softtrigger_kin5.C.

References cnt, fIni, i, and m.

 {
     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 372 of file softtrigger_kin5.C.

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

Referenced by softtrigger_kin5().

 {
     // *** 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 618 of file softtrigger_kin5.C.

References fRand.

Referenced by makeRecoCands().

 {
     // 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 317 of file softtrigger_kin5.C.

Referenced by makeRecoCands(), and softtrigger_kin5().

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

double Legendre	(	int	l,
		double	x
	)

Definition at line 212 of file softtrigger_kin5.C.

References fabs(), and x.

Referenced by compFoxWolfMoms().

 {
     assert(fabs(x) <= 1.);
 
         if (l==0) return 1.;
         
     double pmm = 1.;
         double pmmp1 = x;
 
         if (l>1)
         {
                 for(int ll=2; ll<=l; ll++)
                 {
                   double pll = (x * (2 * ll - 1) * pmmp1 - (ll - 1) * pmm) / (ll);
                   pmm = pmmp1;
                   pmmp1 = pll;
                 }
         }
    return pmmp1;
 }

void makeIni4Vector	(	TLorentzVector &	l,
		double	s
	)

Definition at line 719 of file softtrigger_kin5.C.

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

Referenced by softtrigger_kin5().

 {
     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 698 of file softtrigger_kin5.C.

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

Referenced by softtrigger_kin5().

 {
     // *** 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 332 of file softtrigger_kin5.C.

References i, pidtab, and printf().

Referenced by init().

 {
     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 625 of file softtrigger_kin5.C.

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

Referenced by softtrigger_kin5().

 {
     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) && c.P4().E()>fNeutThrsh)  // *** neutral
         {
             smearMom(c,fdE, fNdtht, fNdphi);
             reco.push_back(c);
         }
     }
 }

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

Definition at line 444 of file softtrigger_kin5.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 420 of file softtrigger_kin5.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 403 of file softtrigger_kin5.C.

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

Referenced by printList().

 {
     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 412 of file softtrigger_kin5.C.

References printCand().

Referenced by printAllLists().

 {
     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 661 of file softtrigger_kin5.C.

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

Referenced by makePidSelection().

 {
     out.clear();
 
     TParticlePDG *part = pdg_db->GetParticle(pdg);
 
     //pdg *= pdg<20 ? -chrg : chrg;
 
     if (part) mass = part->Mass();
     if (mass<0) mass = 0.13957;
 
     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 651 of file softtrigger_kin5.C.

References fabs(), and i.

Referenced by softtrigger_kin5().

 {
     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 588 of file softtrigger_kin5.C.

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

Referenced by makeRecoCands().

 {
     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_kin5	(	TString	fsig,
		int	mode = `0`,
		double	sqrts = `3.77`,
		int	nev = `1000`,
		bool	evcut = `false`,
		double	dp = `0.03`,
		double	trkeff = `100.`,
		double	pideff = `95.`,
		double	misid = `5.`,
		double	P_mix = `0.00`
	)

Definition at line 727 of file softtrigger_kin5.C.

References EventShape::Aplanarity(), boostList(), c, EventShape::ChrgPSumCms(), EventShape::ChrgPSumLab(), EventShape::ChrgPtSumCms(), EventShape::ChrgPtSumLab(), EventShape::Circularity(), combine(), Double_t, eminus, eplus, f, fIni, fKsMass, fKsPdg, fKswin, fMAX, EventShape::FoxWolfMomR(), fPi0Mass, fPi0Pdg, fPi0win, fRand, gam, high, i, init(), isFinalState(), kminus, kplus, m, makeIni4Vector(), makePidSelection(), makeRecoCands(), max(), mcFSlist, mclist, min(), mode, EventShape::MultPminCms(), EventShape::MultPminLab(), muminus, muplus, name, EventShape::NCharged(), EventShape::NeutESumCms(), EventShape::NeutESumLab(), EventShape::NeutEtSumCms(), EventShape::NeutEtSumLab(), EventShape::NNeutral(), EventShape::NParticles(), p, SimpleCand::P4(), pdg_db, piminus, piplus, EventShape::Planarity(), EventShape::PmaxCms(), EventShape::PmaxLab(), EventShape::PminCms(), EventShape::PminLab(), pminus, pplus, EventShape::PRapmax(), pt(), EventShape::Ptmax(), EventShape::Ptmin(), EventShape::PtSumCms(), EventShape::PtSumLab(), pz, recocmslist, recolist, selectMass(), SimpleCand::SetDau(), SimpleCand::SetEvtId(), SimpleCand::SetMotherIdx(), SimpleCand::SetNSiblings(), EventShape::Sphericity(), EventShape::SumChrgPminCms(), EventShape::SumChrgPminLab(), EventShape::SumNeutEminCms(), EventShape::SumNeutEminLab(), EventShape::SumPtminLab(), t, EventShape::Thrust(), and TString.

 {
     fRand.SetSeed();
 
     init(pideff, misid);
 
     unsigned int i,j,k;
 
     makeIni4Vector(fIni, sqrts);
 
     //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
         bool bkg=dpm || fsig.Contains("BKG");     // *** Background events from EvtGen OR 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;
 
     double min=0, max=5.5;
 
     TString name = "ntp_";
     if (!bkg) name+="sig_mode_";
     else name+="bkg_mode_";
     name+=mode; name+=".root";
         
         TFile *ggg=new TFile(name, "RECREATE");
 
         //branches
     //                   0   1     2     3     4   5   6   7  8   9   10   11    12   13   14   15   16   17   18   19   20    21    22    23    24    25    26    27   28  29  30  31  32  33
         TString branches = "sig:mode:npart:nneut:nchrg:ne:nmu:npi:nk:npr:pmax:ptmax:np05:np10:np15:np20:np25:np30:np40:np50:npt05:npt10:npt15:npt20:npt25:npt30:npt40:npt50:sph:apl:pla:thr:sne:scp:";
 
     //                  34  35  36   37     38    39    40    41    42    43    44    45    46    47   48   49    50    51    52    53   54   55    56
         branches        += "fw1:fw2:fw3:pmaxl:ptmaxl:np05l:np10l:np15l:np20l:np25l:np30l:np40l:np50l:ne10:ne20:nmu10:nmu20:npi10:npi20:nk10:nk20:npr10:npr20:";
         
     //                   57     58     59     60     61      62      63    64    65     66    67  68  69  70   71    72     73      74      75      76      77       78      79
         branches        += "sumpt:sumptl:sumptc:sumptcl:sumetn:sumetnl:sumpc:sumpcl:sumen:sumenl:cir:fw4:fw5:pmin:ptmin:pminl:prapmax:sumpt05:sumpt10:sumpc05:sumpc10:sumpc05l:sumpc10l:";
     
         //                    80        81       82      83    84   85    86     87     88     89
         branches        += "sumen05:sumen10:sumen05l:sumen10l:npi0:nks0:npi005:nks005:npi010:nks010";
         
     TNtuple *ntpev = new TNtuple("ntpev","ntpev", branches);
 
     Float_t ntpval[100];
 
     // *** 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);
     }
 
 
     // *** the event loop
 
     int mixevts=1;
     int currmix=0;
 
     // how many events are tagged by individual algo?
     bool channelselect[6];
     double channelcount[6];
     TString channelname[6]={"J/psi","D+-","Ds","Phi","D0","Lamc"};
 
     for (i=0;i<6;++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)
         {
             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;
             }
         }
 
                 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;
 
                 // *** the boost of the system
                 TVector3 boost = fIni.BoostVector();
 
         // *** prepare the toy reco: smear momentum and apply track efficiency (both flat at the moment)
         makeRecoCands(mcFSlist, recolist, trkeff/100., dp );
 
 /*              int nneut = 0;
                 int nchrg = 0;
 
                 for (j=0;j<recolist.size();++j)
                 {
                         int chrg = recolist[j].Charge();
                         if (abs(chrg)==0) nneut++; else nchrg++;
                 }
 */              
                 boostList(recolist, recocmslist, boost);
 
         // *** prepare all pid list according to the pidtable probabilities
         makePidSelection(recolist);
                 
                 EventShape evs(recolist, fIni);
                 
                 // *** the lists we need for combinatorics
                 CandList comb1;
                 CandList LPi0, LKs;                             // pi0 -> gg (mass cut), KS -> pi+ pi- (mass cut)
                 
                 // *** 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);
                 
                 // *** Make Event Selection
 //              if (evcut)
 //              {
 //                      EventShape evshape(recolist, fIni);
 //                      if (evshape.ChrgPtSumLab()<1.5) continue;
 //              }
 
         //        0   1     2     3     4   5   6   7  8  9   10   11    12   13   14   15   16   17   18   19   20    21    22    23    24    25    26    27   28  29  30  31  32  33
             //      "sig:mode:npart:nneut:nchrg:ne:nmu:npi:nk:npr:pmax:ptmax:np05:np10:np15:np20:np25:np30:np40:np50:npt05:npt10:npt15:npt20:npt25:npt30:npt40:npt50:sph:apl:pla:thr:sne:scp";
 
         //       34  35  36   37     38    39    40    41    42    43    44    45    46    47   48   49    50    51    52    53   54   55    56
             //      "fw1:fw2:fw3:pmaxl:ptmaxl:np05l:np10l:np15l:np20l:np25l:np30l:np40l:np50l:ne10:ne20:nmu10:nmu20:npi10:npi20:nk10:nk20:npr10:npr20");
         
         //        57     58     59     60     61      62      63    64    65     66    67  68  69  70   71    72     73      74      75      76      77       78      79
             //      "sumpt:sumptl:sumptc:sumptcl:sumetn:sumetnl:sumpc:sumpcl:sumen:sumenl:cir:fw4:fw5:pmin:ptmin:pminl:prapmax:sumpt05:sumpt10:sumpc05:sumpc10:sumpc05l:sumpc10l:";
     
             //       80        81       82      83    84   85    86     87     88     89
             //   "sumen05:sumen10:sumen05l:sumen10l:npi0:nks0:npi005:nks005:npi010:nks010";
                 
                 ntpval[0]  = !bkg;                          // signal or background events
                 ntpval[1]  = mode;                          // decay mode (parameter of script - set by hand)
                 
                 ntpval[2]  = evs.NParticles();              // #particles (multiplicity)
                 ntpval[3]  = evs.NNeutral();                // #neutrals cands
                 ntpval[4]  = evs.NCharged();                // #charged cands
                 
                 ntpval[5]  = eplus.size()+eminus.size();    // #electron candidates
                 ntpval[6]  = muplus.size()+muminus.size();  // #muon candidates
                 ntpval[7]  = piplus.size()+piminus.size();  // #pion candidates
                 ntpval[8]  = kplus.size()+kminus.size();    // #kaon candidates
                 ntpval[9]  = pplus.size()+pminus.size();    // #proton candidates
                 ntpval[84] = LPi0.size();                                       // #pi0 candidates
                 ntpval[85] = LKs.size();                                        // #KS0 candidates
                 
                 ntpval[37] = evs.PmaxLab();                 // maximum momentum (lab)
                 ntpval[38] = evs.Ptmax();                   // maximum pt (lab)
                 ntpval[72] = evs.PminLab();                 // minimum momentum (lab)
                 ntpval[71] = evs.Ptmin();                   // minimum pt (lab)
                 ntpval[73] = evs.PRapmax();                 // maximum pseudo rapidity (lab)
                 
                 ntpval[39] = evs.MultPminLab(0.5);          // #particles with p>0.5 GeV (lab)
                 ntpval[40] = evs.MultPminLab(1.0);          // #particles with p>1.0 GeV (lab)
                 ntpval[41] = evs.MultPminLab(1.5);          // #particles with p>1.5 GeV (lab)
                 ntpval[42] = evs.MultPminLab(2.0);          // #particles with p>2.0 GeV (lab)
                 ntpval[43] = evs.MultPminLab(2.5);          // #particles with p>2.5 GeV (lab)
                 ntpval[44] = evs.MultPminLab(3.0);          // #particles with p>3.0 GeV (lab)
                 ntpval[45] = evs.MultPminLab(4.0);          // #particles with p>4.0 GeV (lab)
                 ntpval[46] = evs.MultPminLab(5.0);          // #particles with p>5.0 GeV (lab)
                 
                 ntpval[58] = evs.PtSumLab();                            // sum pt of all particles (lab)
                 ntpval[60] = evs.ChrgPtSumLab();                        // sum pt of chrg particles (lab)
                 ntpval[62] = evs.NeutEtSumLab();                        // sum Et of neut particles (lab)
                 ntpval[64] = evs.ChrgPSumLab();                         // sum p of chrg particles (lab)
                 ntpval[66] = evs.NeutESumLab();                         // sum E of neut particles (lab)
 
                 ntpval[78] = evs.SumChrgPminLab(0.5);           // sum p of charged particles with p>0.5 (lab)
                 ntpval[79] = evs.SumChrgPminLab(1.0);           // sum p of charged particles with p>1.0 (lab)
                 
                 ntpval[82] = evs.SumNeutEminLab(0.5);           // sum E of neutral particles with E>0.5 (lab)
                 ntpval[83] = evs.SumNeutEminLab(1.0);           // sum E of neutral particles with E>1.0 (lab)
 
                 ntpval[58] = evs.PtSumLab();                            // sum pt of all particles (lab)
                 ntpval[74] = evs.SumPtminLab(0.5);                      // sum pt of all particles with pt>0.5 GeV
                 ntpval[75] = evs.SumPtminLab(1.0);                      // sum pt of all particles with pt>1.0 GeV
                                          
                 // all following values are in center-of-momentum system
                 ntpval[10] = evs.PmaxCms();                 // maximum momentum
                 ntpval[11] = evs.Ptmax();                   // maximum pt
                 ntpval[70] = evs.PminCms();                 // minimum momentum
                 
                 ntpval[12] = evs.MultPminCms(0.5);          // #particles with p>0.5 GeV
                 ntpval[13] = evs.MultPminCms(1.0);          // #particles with p>1.0 GeV
                 ntpval[14] = evs.MultPminCms(1.5);          // #particles with p>1.5 GeV
                 ntpval[15] = evs.MultPminCms(2.0);          // #particles with p>2.0 GeV
                 ntpval[16] = evs.MultPminCms(2.5);          // #particles with p>2.5 GeV
                 ntpval[17] = evs.MultPminCms(3.0);          // #particles with p>3.0 GeV
                 ntpval[18] = evs.MultPminCms(4.0);          // #particles with p>4.0 GeV
                 ntpval[19] = evs.MultPminCms(5.0);          // #particles with p>5.0 GeV
                                 
                 ntpval[57] = evs.PtSumCms();                            // sum pt of all particles 
                 ntpval[59] = evs.ChrgPtSumCms();                        // sum pt of chrg particles 
                 ntpval[61] = evs.NeutEtSumCms();                        // sum Et of neut particles 
                 ntpval[63] = evs.ChrgPSumCms();                         // sum p of chrg particles
                 ntpval[65] = evs.NeutESumCms();                         // sum E of neut particles
                 
                 ntpval[76] = evs.SumChrgPminCms(0.5);           // sum p of charged particles with p>0.5 (cms)
                 ntpval[77] = evs.SumChrgPminCms(1.0);           // sum p of charged particles with p>1.0 (cms)
                 
                 ntpval[80] = evs.SumNeutEminCms(0.5);           // sum E of neutral particles with E>0.5 (cms)
                 ntpval[81] = evs.SumNeutEminCms(1.0);           // sum E of neutral particles with E>1.0 (cms)
                 
                 ntpval[20] = 0.0;                           // #particles with pt>0.5 GeV
                 ntpval[21] = 0.0;                           // #particles with pt>1.0 GeV
                 ntpval[22] = 0.0;                           // #particles with pt>1.5 GeV
                 ntpval[23] = 0.0;                           // #particles with pt>2.0 GeV
                 ntpval[24] = 0.0;                           // #particles with pt>2.5 GeV
                 ntpval[25] = 0.0;                           // #particles with pt>3.0 GeV
                 ntpval[26] = 0.0;                           // #particles with pt>4.0 GeV
                 ntpval[27] = 0.0;                           // #particles with pt>5.0 GeV
 
                 ntpval[47] = 0.0;                                                       // #electrons with p>1.0
                 ntpval[48] = 0.0;                                                       // #electrons with p>2.0
 
                 ntpval[49] = 0.0;                                                       // #muons with p>1.0
                 ntpval[50] = 0.0;                                                       // #muons with p>2.0
                 
                 ntpval[51] = 0.0;                                                       // #pions with p>1.0
                 ntpval[52] = 0.0;                                                       // #pions with p>2.0
 
                 ntpval[53] = 0.0;                                                       // #kaons with p>1.0
                 ntpval[54] = 0.0;                                                       // #kaons with p>2.0
                 
                 ntpval[55] = 0.0;                                                       // #protons with p>1.0
                 ntpval[56] = 0.0;                                                       // #protons with p>2.0
                 
                 ntpval[86] = 0.0;                                                       // #pi0s with p>0.5
                 ntpval[88] = 0.0;                                                       // #pi0s with p>1.0
                 
                 ntpval[87] = 0.0;                                                       // #KSs with p>0.5
                 ntpval[89] = 0.0;                                                       // #KSs with p>1.0
                 
                 //compSphericity(recocmslist);
                 
                 //TVector3 thrAxis;
                 //_thr = compThrust(recocmslist, thrAxis);
                 
                 // quantities based on momentum tensor
                 ntpval[28] = evs.Sphericity();                          // sphericity = 3/2 * (lam_2+lam_3)
                 ntpval[29] = evs.Aplanarity();                          // aplanarity = 3/2 * lam_3
                 ntpval[30] = evs.Planarity();                           // planarity  = lam_2 - lam_3
                 ntpval[67] = evs.Circularity();                         // circularity = 2 * min(lam_1, lam_2)/(lam_1 + lam_2)
                 ntpval[31] = evs.Thrust();                                      // thrust
                 
                 //compFoxWolfMoms(recocmslist);
                 
                 ntpval[34] = evs.FoxWolfMomR(1);                                        // Fox Wolfram Moment R1 = H1/H0
                 ntpval[35] = evs.FoxWolfMomR(2);                                        // Fox Wolfram Moment R2 = H2/H0
                 ntpval[36] = evs.FoxWolfMomR(3);                                        // Fox Wolfram Moment R3 = H3/H0
                 ntpval[68] = evs.FoxWolfMomR(4);                                        // Fox Wolfram Moment R4 = H4/H0
                 ntpval[69] = evs.FoxWolfMomR(5);                                        // Fox Wolfram Moment R5 = H5/H0
                 
                 // count electrons
                 for (j=0; j<eplus.size(); ++j) 
                 {
                         if (eplus[j].P4().P()>1.0) ntpval[47]+=1.0;
                         if (eplus[j].P4().P()>2.0) ntpval[48]+=1.0;
                 }
                 for (j=0; j<eminus.size(); ++j) 
                 {
                         if (eminus[j].P4().P()>1.0) ntpval[47]+=1.0;
                         if (eminus[j].P4().P()>2.0) ntpval[48]+=1.0;
                 }
                 
                 // count muons
                 for (j=0; j<muplus.size(); ++j) 
                 {
                         if (muplus[j].P4().P()>1.0) ntpval[49]+=1.0;
                         if (muplus[j].P4().P()>2.0) ntpval[50]+=1.0;
                 }
                 for (j=0; j<muminus.size(); ++j) 
                 {
                         if (muminus[j].P4().P()>1.0) ntpval[49]+=1.0;
                         if (muminus[j].P4().P()>2.0) ntpval[50]+=1.0;
                 }
                 
                 // count pions
                 for (j=0; j<piplus.size(); ++j) 
                 {
                         if (piplus[j].P4().P()>1.0) ntpval[51]+=1.0;
                         if (piplus[j].P4().P()>2.0) ntpval[52]+=1.0;
                 }
                 for (j=0; j<piminus.size(); ++j) 
                 {
                         if (piminus[j].P4().P()>1.0) ntpval[51]+=1.0;
                         if (piminus[j].P4().P()>2.0) ntpval[52]+=1.0;
                 }
                 
                 // count kaons
                 for (j=0; j<kplus.size(); ++j) 
                 {
                         if (kplus[j].P4().P()>1.0) ntpval[53]+=1.0;
                         if (kplus[j].P4().P()>2.0) ntpval[54]+=1.0;
                 }
                 for (j=0; j<kminus.size(); ++j) 
                 {
                         if (kminus[j].P4().P()>1.0) ntpval[53]+=1.0;
                         if (kminus[j].P4().P()>2.0) ntpval[54]+=1.0;
                 }
                 
                 // count protons
                 for (j=0; j<pplus.size(); ++j) 
                 {
                         if (pplus[j].P4().P()>1.0) ntpval[55]+=1.0;
                         if (pplus[j].P4().P()>2.0) ntpval[56]+=1.0;
                 }
                 for (j=0; j<pminus.size(); ++j) 
                 {
                         if (pminus[j].P4().P()>1.0) ntpval[55]+=1.0;
                         if (pminus[j].P4().P()>2.0) ntpval[56]+=1.0;
                 }
                 
                 // count pi0
                 for (j=0; j<LPi0.size(); ++j) 
                 {
                         if (LPi0[j].P4().P()>0.5) ntpval[86]+=1.0;
                         if (LPi0[j].P4().P()>1.0) ntpval[88]+=1.0;
                 }
                 // count KS
                 for (j=0; j<LKs.size(); ++j) 
                 {
                         if (LKs[j].P4().P()>0.5) ntpval[87]+=1.0;
                         if (LKs[j].P4().P()>1.0) ntpval[89]+=1.0;
                 }
                 
                 
                 
                 for (j=0;j<recolist.size();++j)
                 {
                         // *** candidate in lab system
                         SimpleCand     cnd    = recolist[j];
                         TLorentzVector lv     = cnd.P4();
                         
                         // *** candidate in cms system
                         SimpleCand     cndcms = recocmslist[j];
                         TLorentzVector lvcm   = cndcms.P4();
                         
                         // compute in LAB
                         
                         
                         if (lvcm.Pt()>0.5) ntpval[20]+=1.0;
                         if (lvcm.Pt()>1.0) ntpval[21]+=1.0;
                         if (lvcm.Pt()>1.5) ntpval[22]+=1.0;
                         if (lvcm.Pt()>2.0) ntpval[23]+=1.0;
                         if (lvcm.Pt()>2.5) ntpval[24]+=1.0;
                         if (lvcm.Pt()>3.0) ntpval[25]+=1.0;
                         if (lvcm.Pt()>4.0) ntpval[26]+=1.0;
                         if (lvcm.Pt()>5.0) ntpval[27]+=1.0;
                         
         }
 
                 ntpev->Fill(ntpval);
 
         mixevts=1;
         if (fRand.Rndm()<P_mix) mixevts=2;
     }
 
         
     ggg->cd();
     ntpev->Write();
     ggg->Close();
 
   return 0;
 }

Variable Documentation

double _apl

Definition at line 102 of file softtrigger_kin5.C.

Referenced by compSphericity().

double _foxwolfR[_nfwmom]

Definition at line 108 of file softtrigger_kin5.C.

Referenced by compFoxWolfMoms().

const int _nfwmom = 5

Definition at line 106 of file softtrigger_kin5.C.

Referenced by compFoxWolfMoms().

double _pla

Definition at line 102 of file softtrigger_kin5.C.

Referenced by compSphericity().

TVector3 _secAx

Definition at line 103 of file softtrigger_kin5.C.

double _sph

Definition at line 102 of file softtrigger_kin5.C.

Referenced by compSphericity().

TVector3 _sphAx

Definition at line 103 of file softtrigger_kin5.C.

double _thr

Definition at line 102 of file softtrigger_kin5.C.

TVector3 _thrAx

Definition at line 103 of file softtrigger_kin5.C.

CandList eminus

Definition at line 144 of file softtrigger_kin5.C.

Referenced by ana_multikalman(), PndGiBuuJPsiEE::DecayTrack(), PndMyAnalysisTask::Exec(), makePidSelection(), printAllLists(), and softtrigger_kin5().

CandList eplus

Definition at line 144 of file softtrigger_kin5.C.

Referenced by ana_multikalman(), PndGiBuuJPsiEE::DecayTrack(), PndMyAnalysisTask::Exec(), makePidSelection(), printAllLists(), and softtrigger_kin5().

TF1* f_res =0

Definition at line 130 of file softtrigger_kin5.C.

Referenced by dponline(), and init().

double fD0Mass = 1.8645

Definition at line 64 of file softtrigger_kin5.C.

int fD0Pdg = 421

Definition at line 63 of file softtrigger_kin5.C.

double fD0sig = 1.04

Definition at line 65 of file softtrigger_kin5.C.

double fD0win =0.25

Definition at line 66 of file softtrigger_kin5.C.

const double fdE = 0.03

Definition at line 44 of file softtrigger_kin5.C.

Referenced by makeRecoCands().

double fDMass = 1.8693

Definition at line 59 of file softtrigger_kin5.C.

int fDPdg = 411

Definition at line 58 of file softtrigger_kin5.C.

double fDsig = 0.88

Definition at line 60 of file softtrigger_kin5.C.

double fDsMass = 1.9686

Definition at line 81 of file softtrigger_kin5.C.

int fDsPdg = 431

Definition at line 80 of file softtrigger_kin5.C.

double fDssig = 0.75

Definition at line 82 of file softtrigger_kin5.C.

double fDswin =0.2

Definition at line 83 of file softtrigger_kin5.C.

double fDwin =0.25

Definition at line 61 of file softtrigger_kin5.C.

double fEtacMass = 2.9798

Definition at line 91 of file softtrigger_kin5.C.

int fEtacPdg = 441

Definition at line 90 of file softtrigger_kin5.C.

double fEtacsig = 0.15

Definition at line 92 of file softtrigger_kin5.C.

double fEtacwin =0.15

Definition at line 93 of file softtrigger_kin5.C.

TLorentzVector fIni

Definition at line 34 of file softtrigger_kin5.C.

Referenced by fillHisto(), and softtrigger_kin5().

double fJ3pioff = 0.027

Definition at line 72 of file softtrigger_kin5.C.

double fJMass = 3.09687

Definition at line 69 of file softtrigger_kin5.C.

int fJpsiPdg = 443

Definition at line 68 of file softtrigger_kin5.C.

double fJsig = 2.24

Definition at line 70 of file softtrigger_kin5.C.

double fJsig3pi = 1.5

Definition at line 71 of file softtrigger_kin5.C.

double fJwin =0.7

Definition at line 73 of file softtrigger_kin5.C.

double fKsMass = 0.497614

Definition at line 55 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

int fKsPdg = 310

Definition at line 54 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

double fKswin = 0.025

Definition at line 56 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

double fLamcMass = 2.2849

Definition at line 96 of file softtrigger_kin5.C.

int fLamcPdg = 4122

Definition at line 95 of file softtrigger_kin5.C.

double fLamcsig = 0.15

Definition at line 97 of file softtrigger_kin5.C.

double fLamcwin =0.2

Definition at line 98 of file softtrigger_kin5.C.

double fLamMass = 1.115684

Definition at line 76 of file softtrigger_kin5.C.

int fLamPdg = 3122

Definition at line 75 of file softtrigger_kin5.C.

double fLamsig = 0.08

Definition at line 77 of file softtrigger_kin5.C.

double fLamwin =0.2

Definition at line 78 of file softtrigger_kin5.C.

const int fMAX = 100

Definition at line 37 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

const double fNdphi = 0.003

Definition at line 46 of file softtrigger_kin5.C.

Referenced by makeRecoCands().

const double fNdtht = 0.003

Definition at line 45 of file softtrigger_kin5.C.

Referenced by makeRecoCands().

const double fNeutEff = 0.9

Definition at line 43 of file softtrigger_kin5.C.

Referenced by makeRecoCands().

const double fNeutThrsh = 0.02

Definition at line 47 of file softtrigger_kin5.C.

Referenced by makeRecoCands().

double fPhiMass = 1.0195

Definition at line 86 of file softtrigger_kin5.C.

int fPhiPdg = 333

Definition at line 85 of file softtrigger_kin5.C.

double fPhisig = 0.15

Definition at line 87 of file softtrigger_kin5.C.

double fPhiwin =0.2

Definition at line 88 of file softtrigger_kin5.C.

double fPi0Mass = 0.13497

Definition at line 51 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

int fPi0Pdg = 111

Definition at line 50 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

double fPi0win = 0.020

Definition at line 52 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

TRandom3 fRand

Definition at line 33 of file softtrigger_kin5.C.

Referenced by isDetected(), select(), smearMom(), and softtrigger_kin5().

CandList gam

Definition at line 144 of file softtrigger_kin5.C.

Referenced by ana_chic(), ana_dsdsj(), ana_dsdsj_full(), ana_example(), ana_hc(), PndTripleAnaTask::DsDs2317Analysis(), PndSciTAnaIdeal::Exec(), makePidSelection(), printAllLists(), QAmacro_hyp_1(), run(), run_test(), PndHyp::SetHypGamFEm(), softtrigger_kin5(), and PndTripleAnaTask::ThreePiAnalysis().

double high = 5.0

Definition at line 40 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

CandList kminus

Definition at line 144 of file softtrigger_kin5.C.

Referenced by ana_multikalman(), makePidSelection(), printAllLists(), and softtrigger_kin5().

CandList kplus

Definition at line 144 of file softtrigger_kin5.C.

Referenced by ana_multikalman(), makePidSelection(), printAllLists(), and softtrigger_kin5().

double low = 0.0

Definition at line 39 of file softtrigger_kin5.C.

CandList mcFSlist

Definition at line 141 of file softtrigger_kin5.C.

Referenced by softtrigger_kin5().

CandList mclist

Definition at line 141 of file softtrigger_kin5.C.

Referenced by ana_jpsi(), combine(), PndPidCorrelator::ConstructNeutralCandidate(), PndLLbarAnaTask::Exec(), PndParticleQATask::Exec(), PndScrutAnaTask::Exec(), PndAnaWithTrigger::Exec(), PndSimpleCombinerTask::Exec(), PndPidCorrelator::GetEmcInfo(), PndPidCorrelator::GetFscInfo(), printAllLists(), and softtrigger_kin5().

double mE

Definition at line 128 of file softtrigger_kin5.C.

Referenced by PndCATracking::CATrackParToFairTrackParP(), init(), and makePidSelection().

double mK

Definition at line 128 of file softtrigger_kin5.C.

Referenced by PndCATrackParam::BetheBlochGeant(), PndCATrackParamVector::BetheBlochGeant(), PndFTSCATrackParamVector::BetheBlochGeant(), init(), and makePidSelection().

double mMu

Definition at line 128 of file softtrigger_kin5.C.

Referenced by init(), and makePidSelection().

double mPi

Definition at line 128 of file softtrigger_kin5.C.

Referenced by init(), and makePidSelection().

double mPi0

Definition at line 128 of file softtrigger_kin5.C.

Referenced by init().

CandList muminus

CandList muplus

int nbins = 300

Definition at line 38 of file softtrigger_kin5.C.

TDatabasePDG* pdg_db

Definition at line 136 of file softtrigger_kin5.C.

Referenced by init(), select(), and softtrigger_kin5().

std::map<int,int> pdgidx

Definition at line 134 of file softtrigger_kin5.C.

Referenced by init(), and select().

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 120 of file softtrigger_kin5.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 112 of file softtrigger_kin5.C.

Referenced by makePidTable(), and select().