PndRhoTupleQA Class Reference

#include <PndRhoTupleQA.h>

Public Member Functions
	PndRhoTupleQA (PndAnalysis *ana=0, double pbarmom=15.15)
	~PndRhoTupleQA ()
void	qaCand (TString pre, RhoCandidate *cc, RhoTuple *n, bool skip=false)
void	qaP4 (TString pre, TLorentzVector c, RhoTuple *n, bool skip=false)
void	qaP4Cms (TString pre, TLorentzVector c, RhoTuple *n, bool skip=false)
void	qaP4Cov (TString pre, RhoCandidate *c, RhoTuple *n, bool skip=false)
void	qaPos (TString pre, TVector3 p, RhoTuple *n, bool skip=false)
void	qaPull (TString pre, RhoCandidate *c, RhoTuple *n, bool skip=false)
void	qaComp (TString pre, RhoCandidate *c, RhoTuple *n, bool covs=false, bool pulls=false)
void	qaKs0 (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaPi0 (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaEventShape (TString pre, PndEventShape *evsh, RhoTuple *n)
void	qaEventShapeShort (TString pre, PndEventShape *evsh, RhoTuple *n)
void	qaESPidMult (TString pre, PndEventShape *evsh, double prob, double pmin, RhoTuple *n)
void	qaESMult (TString pre, PndEventShape *evsh, RhoTuple *n)
void	qaESSum (TString pre, PndEventShape *evsh, RhoTuple *n)
void	qaESMinMax (TString pre, PndEventShape *evsh, RhoTuple *n)
void	qaESEventVars (TString pre, PndEventShape *evsh, RhoTuple *n)
void	qaVtx (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaPoca (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaPRG (TString pre, RhoCandidate *c, RhoTuple *n)
void	qa2Body (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaDalitz (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaTrk (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaPid (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaEmc (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaMvd (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaStt (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaGem (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaDrc (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaDsc (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaRich (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaTof (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaMuo (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaRecoShort (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaRecoShortTree (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaRecoFull (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaRecoFullTree (TString pre, RhoCandidate *c, RhoTuple *n)
void	qaMcList (TString pre, RhoCandList &l, RhoTuple *n, int max=10000)
void	qaMcList (RhoTuple *n, int max=10000)
void	qaMcDiff (TString pre, RhoCandidate *c, RhoTuple *n, bool skip=false)
void	qaMc (TString pre, RhoCandidate *c, RhoTuple *n, bool skip=false)
void	qaFitter (TString pre, RhoFitterBase *fitter, RhoTuple *n, bool skip=false)
void	SetIniP4 (TLorentzVector _ini)
TLorentzVector	GetIniP4 ()

Private Attributes
RhoCandidate *	fDummyCand
RhoVtxPoca *	fVtxPoca
PndAnalysis *	fAnalysis
TLorentzVector	fIniP4

Detailed Description

Definition at line 16 of file PndRhoTupleQA.h.

Constructor & Destructor Documentation

PndRhoTupleQA::PndRhoTupleQA	(	PndAnalysis *	ana = 0,
		double	pbarmom = 15.15
	)

Definition at line 16 of file PndRhoTupleQA.cxx.

References fAnalysis, fDummyCand, fIniP4, fVtxPoca, mp, and sqrt().

{
  fDummyCand = new RhoCandidate();
  fAnalysis = ana;
  fVtxPoca = new RhoVtxPoca();
  double mp=0.938272;
  fIniP4.SetXYZT(0,0,pbarmom, sqrt(pbarmom*pbarmom+mp*mp)+mp);
}

PndRhoTupleQA::~PndRhoTupleQA

(

)

Definition at line 26 of file PndRhoTupleQA.cxx.

References fDummyCand, and fVtxPoca.

{
  delete fVtxPoca;
  delete fDummyCand;
}

Member Function Documentation

TLorentzVector PndRhoTupleQA::GetIniP4 ( )

inline

Definition at line 82 of file PndRhoTupleQA.h.

References fIniP4.

{return fIniP4;}

void PndRhoTupleQA::qa2Body	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 724 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), cos(), d0, RhoCandidate::Daughter(), f, RhoCandidate::NDaughters(), RhoCandidate::P3(), and RhoCandidate::P4().

Referenced by qaComp().

{
  if (n==0) return;

  if (c->NDaughters()!=2) return;

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

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

  // decay angle
  TLorentzVector d_cms = d0->P4();
  d_cms.Boost(-(c->P4().BoostVector()));
  Float_t dec  = d_cms.Vect().Angle(c->P3());
  Float_t cdec = cos(dec);

  n->Column(pre+"oang",    (Float_t) oang,    0.0f );
  n->Column(pre+"decang",  (Float_t) dec,    0.0f );
  n->Column(pre+"cdecang", (Float_t) cdec,    0.0f );
}

void PndRhoTupleQA::qaCand	(	TString	pre,
		RhoCandidate *	cc,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 554 of file PndRhoTupleQA.cxx.

References c, RhoCandidate::Charge(), RhoTuple::Column(), f, p, RhoCandidate::P4(), RhoCandidate::PdgCode(), RhoCandidate::Pos(), qaP4(), and qaPos().

Referenced by PndTripleAnaTask::DsDs2317Analysis(), PndTripleAnaTask::JpsiAnalysis(), qaComp(), qaKs0(), qaMc(), and qaPi0().

{
  if (n==0) return;

  if (!skip)
  {
    TLorentzVector c=cc->P4();
    TVector3       p=cc->Pos();

    qaP4(pre, c, n);
    qaPos(pre, p, n);
    n->Column(pre+"chg", (Float_t) cc->Charge(),   0.0f );
    n->Column(pre+"pdg", (Float_t) cc->PdgCode(),  0.0f );
  }
  else
  {
    TLorentzVector dummy;
    TVector3 dummypos;
    qaP4(pre, dummy, n, true);
    qaPos(pre, dummypos, n, true);
    n->Column(pre+"chg", (Float_t) -999.,       0.0f );
    n->Column(pre+"pdg", (Float_t) -999.,      0.0f );
  }
}

void PndRhoTupleQA::qaComp	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n,
		bool	covs = false,
		bool	pulls = false
	)

Definition at line 283 of file PndRhoTupleQA.cxx.

References RhoCandidate::Charge(), RhoTuple::Column(), RhoCandidate::Daughter(), f, fabs(), fAnalysis, RhoCandidate::GetMcTruth(), i, RhoCandidate::M(), PndAnalysis::McTruthMatch(), n, name, RhoCandidate::NDaughters(), RhoCandidate::P4(), RhoCandidate::PdgCode(), qa2Body(), qaCand(), qaDalitz(), qaP4Cms(), qaP4Cov(), qaPid(), qaPoca(), qaPull(), qaVtx(), and TString.

Referenced by ana_jpsi(), PndTripleAnaTask::DsDs2317Analysis(), PndTripleAnaTask::JpsiAnalysis(), PndSoftTriggerTask::TagMode(), PndTripleAnaTask::ThreePiAnalysis(), and tut_ana_d0_qahelper().

{
  if (n==0) return;

  // what kind of particle?
  //int pdg = c->PdgCode(); //[R.K. 01/2017] unused variable

  // special composite particle?
  // pi0 or eta?
  //if (pdg == 111 || pdg == 221)
  //{
    //qaPi0(pre, c, n);
    //return;
  //}
  //if (pdg == 310)
  //{
    //qaKs0(pre, c, n);
    //return;
  //}

  // how many daughters?
  int nd = c->NDaughters();

  // truth match already done?
  RhoCandidate *truth = c->GetMcTruth();

  // if not, try one
  if (truth == 0 && fAnalysis != 0)
  {
    fAnalysis->McTruthMatch(c);
    truth = c->GetMcTruth();
  }

  bool mct = false;
  if (truth!=0)
  {
    if (nd>0) mct = true;
    else mct = (truth->PdgCode()==c->PdgCode());
  }
  // store cand info in lab and cms
  qaCand(pre,  c,  n);
  qaP4Cms(pre, c->P4(), n);
  if(covs)qaP4Cov(pre,  c,  n);
  if(pulls) qaPull(pre,c,n, (0==truth));
  n->Column(pre+"mct",  (Float_t) mct, 0.0f);

  // for mass difference e.g. D* -> D pi decays
  if (nd>=2) n->Column(pre+"mdif", (Float_t) (c->M() - c->Daughter(0)->M()), 0.0f);

  // cand is final state particle
  if (nd==0)
  {
    // charged particle -> store PID info
    if ( fabs(c->Charge())>0.1 ) qaPid(pre,  c,  n);

    // and pdg code from truth
    Float_t trpdg = 0.;
    if (truth) trpdg = truth->PdgCode();
    n->Column(pre+"trpdg", (Float_t) trpdg, 0.0f);
  }
  // cand is composite
  else
  {
    // if 2 daughters -> store decay angle etc
    if (nd==2) qa2Body(pre, c, n);
    // if 3 daughters -> dalitz plot vars
    if (nd==3) qaDalitz(pre, c, n);

    // counter for charged final states
    int nchrgfs = 0;

    for (int i=0; i<nd; ++i)
    {
      RhoCandidate *dau = c->Daughter(i);
      TString name=TString::Format("%sd%d",pre.Data(),i);

      // count charged final states for possible vertexing later
      if (dau->NDaughters()==0 && fabs(dau->Charge())>0.01) nchrgfs++;

      // recursive call of qaComp
      qaComp(name, dau, n, covs, pulls);
    }
    // only charged final state daughters -> Vtx info with RhoVtxPoca
    if (nchrgfs > 1)
    {
      qaPoca(pre, c, n);
      qaVtx(pre, c, n);
    }
  }
}

void PndRhoTupleQA::qaDalitz	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 749 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), RhoCandidate::Daughter(), RhoCandidate::NDaughters(), and RhoCandidate::P4().

Referenced by qaComp().

{
  if (n==0) return;

  if (c->NDaughters()!=3) return;

  TLorentzVector l01 = c->Daughter(0)->P4()+c->Daughter(1)->P4();
  TLorentzVector l12 = c->Daughter(1)->P4()+c->Daughter(2)->P4();
  TLorentzVector l02 = c->Daughter(2)->P4()+c->Daughter(0)->P4();

  n->Column(pre+"m01", (Float_t) l01.M(),    0.0f );
  n->Column(pre+"m12", (Float_t) l12.M(),    0.0f );
  n->Column(pre+"m02", (Float_t) l02.M(),    0.0f );

  n->Column(pre+"dal01", (Float_t) l01.M2(),    0.0f );
  n->Column(pre+"dal12", (Float_t) l12.M2(),    0.0f );
  n->Column(pre+"dal02", (Float_t) l02.M2(),    0.0f );
}

void PndRhoTupleQA::qaDrc	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 893 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndPidCandidate::GetDrcIndex(), PndPidCandidate::GetDrcNumberOfPhotons(), PndPidCandidate::GetDrcQuality(), PndPidCandidate::GetDrcThetaC(), PndPidCandidate::GetDrcThetaCErr(), and RhoCandidate::GetRecoCandidate().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"drcthtc",    (Float_t) mic->GetDrcThetaC(),      0.0f );
    n->Column(pre+"drcdthtc",   (Float_t) mic->GetDrcThetaCErr(),    0.0f );
    n->Column(pre+"drcqa",    (Float_t) mic->GetDrcQuality(),      0.0f );
    n->Column(pre+"drcnphot",    (Int_t)   mic->GetDrcNumberOfPhotons(),  0 );
    n->Column(pre+"drcidx",    (Int_t)   mic->GetDrcIndex(),      0 );
  }
}

void PndRhoTupleQA::qaDsc	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 911 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndPidCandidate::GetDiscIndex(), PndPidCandidate::GetDiscNumberOfPhotons(), PndPidCandidate::GetDiscQuality(), PndPidCandidate::GetDiscThetaC(), PndPidCandidate::GetDiscThetaCErr(), and RhoCandidate::GetRecoCandidate().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"dscthtc",    (Float_t) mic->GetDiscThetaC(),      0.0f );
    n->Column(pre+"dscdthtc",    (Float_t) mic->GetDiscThetaCErr(),    0.0f );
    n->Column(pre+"dscqa",    (Float_t) mic->GetDiscQuality(),    0.0f );
    n->Column(pre+"dscnphot",    (Int_t)   mic->GetDiscNumberOfPhotons(),0 );
    n->Column(pre+"dscidx",    (Int_t)   mic->GetDiscIndex(),      0 );
  }
}

void PndRhoTupleQA::qaEmc	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 823 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndPidCandidate::GetEmcCalEnergy(), PndPidCandidate::GetEmcClusterE1(), PndPidCandidate::GetEmcClusterE25(), PndPidCandidate::GetEmcClusterE9(), PndPidCandidate::GetEmcClusterLat(), PndPidCandidate::GetEmcClusterZ20(), PndPidCandidate::GetEmcClusterZ53(), PndPidCandidate::GetEmcIndex(), PndPidCandidate::GetEmcModule(), PndPidCandidate::GetEmcNumberOfBumps(), PndPidCandidate::GetEmcNumberOfCrystals(), PndPidCandidate::GetEmcQuality(), PndPidCandidate::GetEmcRawEnergy(), and RhoCandidate::GetRecoCandidate().

Referenced by PndParticleQATask::Exec(), qaPi0(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"emceraw",  (Float_t) mic->GetEmcRawEnergy(),      0.0f );
    n->Column(pre+"emcecal",  (Float_t) mic->GetEmcCalEnergy(),      0.0f );
    n->Column(pre+"emcqa",    (Float_t) mic->GetEmcQuality(),       0.0f );
    n->Column(pre+"emcnx",    (Int_t)   mic->GetEmcNumberOfCrystals(),  0 );
    n->Column(pre+"emcnb",    (Int_t)   mic->GetEmcNumberOfBumps(),     0 );
    n->Column(pre+"emcz20",   (Float_t) mic->GetEmcClusterZ20(),       0.0f );
    n->Column(pre+"emcz53",   (Float_t) mic->GetEmcClusterZ53(),       0.0f );
    n->Column(pre+"emclat",   (Float_t) mic->GetEmcClusterLat(),       0.0f );
    n->Column(pre+"emce1",    (Float_t) mic->GetEmcClusterE1(),       0.0f );
    n->Column(pre+"emce9",    (Float_t) mic->GetEmcClusterE9(),       0.0f );
    n->Column(pre+"emce25",   (Float_t) mic->GetEmcClusterE25(),       0.0f );
    n->Column(pre+"emcmod",   (Int_t)   mic->GetEmcModule(),         0 );
    n->Column(pre+"emcidx",   (Int_t)   mic->GetEmcIndex(),         0 );
  }
}

void PndRhoTupleQA::qaESEventVars	(	TString	pre,
		PndEventShape *	evsh,
		RhoTuple *	n
	)

Definition at line 82 of file PndRhoTupleQA.cxx.

References PndEventShape::Aplanarity(), PndEventShape::Circularity(), RhoTuple::Column(), PndEventShape::FoxWolfMomR(), PndEventShape::Planarity(), PndEventShape::PRapmax(), PndEventShape::Sphericity(), and PndEventShape::Thrust().

Referenced by qaEventShape(), and qaEventShapeShort().

{
  n->Column(pre+"prapmax",(Float_t) evsh->PRapmax(),    0.0f );
  n->Column(pre+"sph", (Float_t)    evsh->Sphericity(),  0.0f );
  n->Column(pre+"apl", (Float_t)    evsh->Aplanarity(),  0.0f );
  n->Column(pre+"pla", (Float_t)    evsh->Planarity(),  0.0f );
  n->Column(pre+"thr",(Float_t)    evsh->Thrust(),    0.0f );
  n->Column(pre+"cir", (Float_t)    evsh->Circularity(),  0.0f );

  n->Column(pre+"fw1", (Float_t)    evsh->FoxWolfMomR(1),  0.0f );
  n->Column(pre+"fw2", (Float_t)    evsh->FoxWolfMomR(2),  0.0f );
  n->Column(pre+"fw3", (Float_t)    evsh->FoxWolfMomR(3),  0.0f );
  n->Column(pre+"fw4", (Float_t)    evsh->FoxWolfMomR(4),  0.0f );
  n->Column(pre+"fw5", (Float_t)    evsh->FoxWolfMomR(5),  0.0f );
}

void PndRhoTupleQA::qaESMinMax	(	TString	pre,
		PndEventShape *	evsh,
		RhoTuple *	n
	)

Definition at line 70 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndEventShape::PmaxCms(), PndEventShape::PmaxLab(), PndEventShape::PminCms(), PndEventShape::PminLab(), PndEventShape::Ptmax(), and PndEventShape::Ptmin().

Referenced by qaEventShape(), and qaEventShapeShort().

{
  n->Column(pre+"pmax", (Float_t)    evsh->PmaxCms(),  0.0f );
  n->Column(pre+"ptmax",(Float_t)    evsh->Ptmax() ,  0.0f );
  n->Column(pre+"pmaxl",(Float_t)    evsh->PmaxLab() ,  0.0f );
  n->Column(pre+"pmin", (Float_t)    evsh->PminCms(),  0.0f );
  n->Column(pre+"ptmin",(Float_t)   evsh->Ptmin(),  0.0f );
  n->Column(pre+"pminl",(Float_t)    evsh->PminLab(),  0.0f );
}

void PndRhoTupleQA::qaESMult	(	TString	pre,
		PndEventShape *	evsh,
		RhoTuple *	n
	)

Definition at line 44 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndEventShape::NCharged(), PndEventShape::NNeutral(), and PndEventShape::NParticles().

Referenced by qaEventShape(), and qaEventShapeShort().

{
  n->Column(pre+"npart", (Int_t)    evsh->NParticles(),  0 );
  n->Column(pre+"nneut", (Int_t)    evsh->NNeutral(),  0 );
  n->Column(pre+"nchrg", (Int_t)    evsh->NCharged(),  0 );
}

void PndRhoTupleQA::qaESPidMult	(	TString	pre,
		PndEventShape *	evsh,
		double	prob,
		double	pmin,
		RhoTuple *	n
	)

Definition at line 33 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndEventShape::MultElectronPminLab(), PndEventShape::MultKaonPminLab(), PndEventShape::MultMuonPminLab(), PndEventShape::MultPionPminLab(), and PndEventShape::MultProtonPminLab().

Referenced by qaEventShape(), and qaEventShapeShort().

{
  n->Column(pre+"npide", (Int_t) evsh->MultElectronPminLab(prob,pmin),    0 );
  n->Column(pre+"npidmu",(Int_t) evsh->MultMuonPminLab(prob,pmin),      0 );
  n->Column(pre+"npidpi",(Int_t) evsh->MultPionPminLab(prob,pmin),      0 );
  n->Column(pre+"npidk", (Int_t) evsh->MultKaonPminLab(prob,pmin),      0 );
  n->Column(pre+"npidp", (Int_t) evsh->MultProtonPminLab(prob,pmin),      0 );
}

void PndRhoTupleQA::qaESSum	(	TString	pre,
		PndEventShape *	evsh,
		RhoTuple *	n
	)

Definition at line 53 of file PndRhoTupleQA.cxx.

References PndEventShape::ChrgPSumCms(), PndEventShape::ChrgPSumLab(), PndEventShape::ChrgPtSumCms(), PndEventShape::ChrgPtSumLab(), RhoTuple::Column(), PndEventShape::NeutESumCms(), PndEventShape::NeutESumLab(), PndEventShape::NeutEtSumCms(), PndEventShape::NeutEtSumLab(), PndEventShape::PtSumCms(), and PndEventShape::PtSumLab().

Referenced by qaEventShape(), and qaEventShapeShort().

{
  n->Column(pre+"sumpc",   (Float_t)  evsh->ChrgPSumCms(),    0.0f );
  n->Column(pre+"sumpcl", (Float_t)  evsh->ChrgPSumLab(),    0.0f );
  n->Column(pre+"sumen",   (Float_t)  evsh->NeutESumCms(),    0.0f );
  n->Column(pre+"sumenl", (Float_t)  evsh->NeutESumLab(),    0.0f );

  n->Column(pre+"sumpt",   (Float_t)  evsh->PtSumCms(),      0.0f );
  n->Column(pre+"sumptl",  (Float_t) evsh->PtSumLab(),      0.0f );
  n->Column(pre+"sumptc", (Float_t)  evsh->ChrgPtSumCms(),    0.0f );
  n->Column(pre+"sumptcl",(Float_t)  evsh->ChrgPtSumLab(),    0.0f );
  n->Column(pre+"sumetn",  (Float_t) evsh->NeutEtSumCms(),    0.0f );
  n->Column(pre+"sumetnl", (Float_t) evsh->NeutEtSumLab(),    0.0f );
}

void PndRhoTupleQA::qaEventShape	(	TString	pre,
		PndEventShape *	evsh,
		RhoTuple *	n
	)

Definition at line 100 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndEventShape::DetEmcMax(), PndEventShape::DetEmcSum(), PndEventShape::MultChrgPminCms(), PndEventShape::MultChrgPminLab(), PndEventShape::MultNeutEminLab(), PndEventShape::MultPminCms(), PndEventShape::MultPminLab(), PndEventShape::MultPtminCms(), qaESEventVars(), qaESMinMax(), qaESMult(), qaESPidMult(), qaESSum(), PndEventShape::SumChrgPminCms(), PndEventShape::SumChrgPminLab(), PndEventShape::SumNeutEminCms(), PndEventShape::SumNeutEminLab(), and PndEventShape::SumPtminCms().

Referenced by PndSoftTriggerTask::Exec().

{
  if (n==0) return;

  // basic multiplicities
  qaESMult(pre,  evsh, n);

  // PID multiplicities
  qaESPidMult(pre+"l",  evsh, 0.25, 0.0, n);
  qaESPidMult(pre+"l1", evsh, 0.25, 1.0, n);
  qaESPidMult(pre+"t",  evsh, 0.5,  0.0, n);
  qaESPidMult(pre+"t1", evsh, 0.5,  1.0, n);
  qaESPidMult(pre+"vt", evsh, 0.9,  0.0, n);

  // event vars like thrust, sphericity etc
  qaESEventVars(pre, evsh, n);

  // standard sums over pt, p, E for all, neutral, charged in lab, cms
  qaESSum(pre, evsh, n);

  // standard min, max values of p, pt in lab, cms
  qaESMinMax(pre, evsh, n);

  // Multiplicities with min momemtum cut (cms)
  n->Column(pre+"np05", (Int_t)    evsh->MultPminCms(0.5),  0 );
  n->Column(pre+"np10", (Int_t)    evsh->MultPminCms(1.0),  0 );
  n->Column(pre+"np20", (Int_t)    evsh->MultPminCms(2.0),  0 );
  n->Column(pre+"np30", (Int_t)    evsh->MultPminCms(3.0),  0 );
  n->Column(pre+"np40", (Int_t)    evsh->MultPminCms(4.0),  0 );
  n->Column(pre+"np50", (Int_t)    evsh->MultPminCms(5.0),  0 );

  // Multiplicities with min momemtum cut (lab)
  n->Column(pre+"np05l", (Int_t)  evsh->MultPminLab(0.5),  0 );
  n->Column(pre+"np10l", (Int_t)  evsh->MultPminLab(1.0),  0 );
  n->Column(pre+"np20l", (Int_t)  evsh->MultPminLab(2.0),  0 );
  n->Column(pre+"np30l", (Int_t)  evsh->MultPminLab(3.0),  0 );
  n->Column(pre+"np40l", (Int_t)  evsh->MultPminLab(4.0),  0 );
  n->Column(pre+"np50l",   (Int_t)  evsh->MultPminLab(5.0),  0 );

  // Multiplicities with min p_t cut (cms)
  n->Column(pre+"npt05",   (Int_t)  evsh->MultPtminCms(0.5),  0 );
  n->Column(pre+"npt10",   (Int_t)  evsh->MultPtminCms(1.0),  0 );
  n->Column(pre+"npt15",   (Int_t)  evsh->MultPtminCms(1.5),  0 );
  n->Column(pre+"npt20",   (Int_t)  evsh->MultPtminCms(2.0),  0 );
  n->Column(pre+"npt25",   (Int_t)  evsh->MultPtminCms(2.5),  0 );
  n->Column(pre+"npt30",   (Int_t)  evsh->MultPtminCms(3.0),  0 );

  // Neutral multiplicities with min energy cut (lab)
  n->Column(pre+"nne003l",  (Int_t) evsh->MultNeutEminLab(0.03),  0 );
  n->Column(pre+"nne005l",  (Int_t) evsh->MultNeutEminLab(0.05),  0 );
  n->Column(pre+"nne01l",   (Int_t) evsh->MultNeutEminLab(0.1),  0 );
  n->Column(pre+"nne05l",   (Int_t) evsh->MultNeutEminLab(0.5),  0 );

  // Charged multiplicities with min momentum cut (lab)
  n->Column(pre+"ncp005l", (Int_t)    evsh->MultChrgPminLab(0.05),  0 );
  n->Column(pre+"ncp01l", (Int_t)    evsh->MultChrgPminLab(0.1),  0 );
  n->Column(pre+"ncp02l", (Int_t)    evsh->MultChrgPminLab(0.2),  0 );
  n->Column(pre+"ncp05l", (Int_t)    evsh->MultChrgPminLab(0.5),  0 );
  n->Column(pre+"ncp10l", (Int_t)    evsh->MultChrgPminLab(1.0),  0 );

  // Charged multiplicities with min momentum cut (cms)
  n->Column(pre+"ncp005", (Int_t)    evsh->MultChrgPminCms(0.05),  0 );
  n->Column(pre+"ncp01", (Int_t)    evsh->MultChrgPminCms(0.1),  0 );
  n->Column(pre+"ncp02", (Int_t)    evsh->MultChrgPminCms(0.2),  0 );
  n->Column(pre+"ncp05", (Int_t)    evsh->MultChrgPminCms(0.5),  0 );
  n->Column(pre+"ncp10", (Int_t)    evsh->MultChrgPminCms(1.0),  0 );

  // Sum of p_t  with min momentum cut (cms)
  n->Column(pre+"sumpt05", (Float_t) evsh->SumPtminCms(0.5),  0.0f );
  n->Column(pre+"sumpt10", (Float_t) evsh->SumPtminCms(1.0),  0.0f );

  // Sum of charged momenta with min momentum cut (cms)
  n->Column(pre+"sumpc05", (Float_t) evsh->SumChrgPminCms(0.5) ,0.0f );
  n->Column(pre+"sumpc10", (Float_t) evsh->SumChrgPminCms(1.0),0.0f );

  // Sum of charged momenta with min momentum cut (lab)
  n->Column(pre+"sumpc05l", (Float_t)evsh->SumChrgPminLab(0.5) ,0.0f );
  n->Column(pre+"sumpc10l",(Float_t) evsh->SumChrgPminLab(1.0) ,0.0f );

  // Sum of neutral energy with min energy cut (cms)
  n->Column(pre+"sumen05",  (Float_t)evsh->SumNeutEminCms(0.5) ,0.0f );
  n->Column(pre+"sumen10", (Float_t) evsh->SumNeutEminCms(1.0) ,0.0f );

  // Sum of neutral energy with min energy cut (lab)
  n->Column(pre+"sumen05l",(Float_t) evsh->SumNeutEminLab(0.5) ,0.0f );
  n->Column(pre+"sumen10l",(Float_t) evsh->SumNeutEminLab(1.0) ,0.0f );

  // detector specific variables
  n->Column(pre+"detemcsum",(Float_t) evsh->DetEmcSum()  ,0.0f );
  n->Column(pre+"detemcmax",(Float_t) evsh->DetEmcMax()  ,0.0f );
}

void PndRhoTupleQA::qaEventShapeShort	(	TString	pre,
		PndEventShape *	evsh,
		RhoTuple *	n
	)

Definition at line 192 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndEventShape::MultChrgPminLab(), PndEventShape::MultNeutEminLab(), PndEventShape::MultPminCms(), PndEventShape::MultPtminCms(), qaESEventVars(), qaESMinMax(), qaESMult(), qaESPidMult(), qaESSum(), and PndEventShape::SumChrgPminCms().

Referenced by ana_jpsi(), PndTripleAnaTask::DsDs2317Analysis(), PndTripleAnaTask::JpsiAnalysis(), PndSoftTriggerTask::TagMode(), PndTripleAnaTask::ThreePiAnalysis(), and tut_ana_d0_qahelper().

{
  if (n==0) return;
  // *** vars for PID multiplicity
  //int ne, nmu, npi, nk, np;

  // basic multiplicities
  qaESMult(pre,  evsh, n);

  // PID multiplicities
  qaESPidMult(pre+"l",  evsh, 0.25, 0.0, n);
  qaESPidMult(pre+"l1", evsh, 0.25, 1.0, n);

  // event vars like thrust, sphericity etc
  qaESEventVars(pre, evsh, n);

  // standard sums over pt, p, E for all, neutral, charged in lab, cms
  qaESSum(pre, evsh, n);

  // standard min, max values of p, pt in lab, cms
  qaESMinMax(pre, evsh, n);

  // some multiplicities with min momentum cut (lab, cms)
  n->Column(pre+"np10",  (Int_t)  evsh->MultPminCms(1.0),     0 );
  n->Column(pre+"npt10", (Int_t)  evsh->MultPtminCms(1.0),  0 );
  n->Column(pre+"ncp10l", (Int_t) evsh->MultChrgPminLab(1.0), 0 );
  n->Column(pre+"nne10l",(Int_t)  evsh->MultNeutEminLab(1.0), 0 );

  // sum of charged momenta with min momentum cut
  n->Column(pre+"sumpc05", (Float_t) evsh->SumChrgPminCms(0.5) ,0.0f );
}

void PndRhoTupleQA::qaFitter	(	TString	pre,
		RhoFitterBase *	fitter,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 1239 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), f, RhoFitterBase::GetChi2(), RhoFitterBase::GetNdf(), and RhoFitterBase::GetProb().

{
  if(!skip){
    n->Column(pre+"chisq" , (Float_t) fitter->GetChi2(),       0.0f);
    n->Column(pre+"ndf" ,   (Float_t) fitter->GetNdf(),        0.0f);
    n->Column(pre+"prob" ,  (Float_t) fitter->GetProb(),       0.0f);
  } else {
    n->Column(pre+"chisq" , (Float_t) -999.0,       0.0f);
    n->Column(pre+"ndf" ,   (Float_t) -999.0,       0.0f);
    n->Column(pre+"prob" ,  (Float_t) -999.0,       0.0f);
  }
}

void PndRhoTupleQA::qaGem	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 848 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndPidCandidate::GetGemHits(), and RhoCandidate::GetRecoCandidate().

Referenced by qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"gemnhits",  (Int_t) mic->GetGemHits(),      0 );
  }
}

void PndRhoTupleQA::qaKs0	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 412 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), d0, RhoCandidate::Daughter(), f, fAnalysis, RhoCandidate::GetMcTruth(), RhoCandidate::M(), PndAnalysis::McTruthMatch(), RhoCandidate::P(), RhoCandidate::P3(), RhoCandidate::Pos(), qaCand(), qaPid(), qaPoca(), and qaVtx().

Referenced by PndSoftTriggerTask::CreateKs0Cands().

{
  if (n==0) return;

  RhoCandidate *d0 = c->Daughter(0);
  RhoCandidate *d1 = c->Daughter(1);
  double ang = d0->P3().Angle(d1->P3());

  if (fAnalysis!=0) fAnalysis->McTruthMatch(c);
  RhoCandidate *truth = c->GetMcTruth();

  qaCand(pre, c, n);
  qaPoca(pre, c, n);
  qaVtx(pre, c, n);

  n->Column(pre+"oang",(Float_t) ang,    0.0f);

  qaCand(pre+"d0", d0, n);
  qaPid(pre+"d0",  d0, n);

  qaCand(pre+"d1", d1, n);
  qaPid(pre+"d1",  d1, n);

  if (truth!=0)
  {
    TVector3 vdist = truth->Pos() - truth->Daughter(0)->Pos();
    Float_t dist = vdist.Mag();
    Float_t ctau = dist * truth->M() / truth->P();

    qaCand("t"+pre, truth, n);
    n->Column(pre+"mct", 1.0f, 0.0f);
    n->Column("t"+pre+"dist",    (Float_t) dist , 0.0f);
    n->Column("t"+pre+"ctau", (Float_t) ctau , 0.0f);
  }
  else
  {
    qaCand("t"+pre, NULL , n, true);
    n->Column(pre+"mct", 0.0f, 0.0f);
    n->Column("t"+pre+"dist", (Float_t) -999. , 0.0f);
    n->Column("t"+pre+"ctau", (Float_t) -999. , 0.0f);
  }
}

void PndRhoTupleQA::qaMc	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 1078 of file PndRhoTupleQA.cxx.

References fDummyCand, RhoCandidate::GetMcTruth(), and qaCand().

{
  if (n==0) return;
  if (c==0) return;

  RhoCandidate *mct=c->GetMcTruth();

  if (mct){
    qaCand(pre+"mc",mct,n,skip);
  }
  else {
      qaCand(pre+"mc", fDummyCand,n,skip);
  }
}

void PndRhoTupleQA::qaMcDiff	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 1163 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), RhoCandidate::Cov7(), f, RhoCandidate::GetMcTruth(), RhoCandidate::GetPosition(), RhoCandidate::P4(), sqrt(), and v.

{
  if (n==0) return;
  if (c==0) return;

  RhoCandidate *mct=c->GetMcTruth();
  if (!skip && mct)
  {
    TLorentzVector p4=c->P4();
    TLorentzVector mcp4=mct->P4();
    TLorentzVector diff = p4 - mcp4;
    TVector3 v=c->GetPosition();
    TVector3 mcv=mct->GetPosition();
    TVector3 vdiff= v - mcv;
    TMatrixD cov7 = c->Cov7();

    n->Column(pre+"mcdiffvx",     (Float_t) vdiff.x() ,       0.0f );
    n->Column(pre+"mcdiffvy",     (Float_t) vdiff.y() ,       0.0f );
    n->Column(pre+"mcdiffvz",     (Float_t) vdiff.z() ,       0.0f );
    n->Column(pre+"mcdiffpx",     (Float_t) diff.Px(),        0.0f );
    n->Column(pre+"mcdiffpy",     (Float_t) diff.Py(),        0.0f );
    n->Column(pre+"mcdiffpz",     (Float_t) diff.Pz(),        0.0f );
    n->Column(pre+"mcdiffe",      (Float_t) diff.E(),         0.0f );

    n->Column(pre+"mcdiffp",      (Float_t) (p4.P()-mcp4.P()),         0.0f );
    n->Column(pre+"mcdifftht",    (Float_t) (p4.Theta()-mcp4.Theta()), 0.0f );
    n->Column(pre+"mcdiffphi",    (Float_t) (p4.Phi()-mcp4.Phi()),     0.0f );

    n->Column(pre+"mcpullvx",     (Float_t) ( vdiff.x()/sqrt(cov7(0,0)) ),       0.0f );
    n->Column(pre+"mcpullvy",     (Float_t) ( vdiff.y()/sqrt(cov7(1,1)) ),       0.0f );
    n->Column(pre+"mcpullvz",     (Float_t) ( vdiff.z()/sqrt(cov7(2,2)) ),       0.0f );
    n->Column(pre+"mcpullpx",     (Float_t) ( diff.Px()/sqrt(cov7(3,3)) ),       0.0f );
    n->Column(pre+"mcpullpy",     (Float_t) ( diff.Py()/sqrt(cov7(4,4)) ),       0.0f );
    n->Column(pre+"mcpullpz",     (Float_t) ( diff.Pz()/sqrt(cov7(5,5)) ),       0.0f );
    n->Column(pre+"mcpulle",      (Float_t) ( diff.E() /sqrt(cov7(6,6)) ),       0.0f );

    n->Column(pre+"mcerrvx",     (Float_t) sqrt(cov7(0,0)),       0.0f );
    n->Column(pre+"mcerrvy",     (Float_t) sqrt(cov7(1,1)),       0.0f );
    n->Column(pre+"mcerrvz",     (Float_t) sqrt(cov7(2,2)),       0.0f );
    n->Column(pre+"mcerrpx",     (Float_t) sqrt(cov7(3,3)),       0.0f );
    n->Column(pre+"mcerrpy",     (Float_t) sqrt(cov7(4,4)),       0.0f );
    n->Column(pre+"mcerrpz",     (Float_t) sqrt(cov7(5,5)),       0.0f );
    n->Column(pre+"mcerre",      (Float_t) sqrt(cov7(6,6)),       0.0f );
  }
  else
  {
    TLorentzVector dummy;
    n->Column(pre+"mcdiffvx",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffvy",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffvz",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffpx",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffpy",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffpz",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffe",      (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffp",      (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdifftht",    (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcdiffphi",    (Float_t) -999.0,       0.0f );

    n->Column(pre+"mcpullvx",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcpullvy",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcpullvz",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcpullpx",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcpullpy",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcpullpz",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcpulle",      (Float_t) -999.0,       0.0f );

    n->Column(pre+"mcerrvx",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcerrvy",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcerrvz",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcerrpx",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcerrpy",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcerrpz",     (Float_t) -999.0,       0.0f );
    n->Column(pre+"mcerre",      (Float_t) -999.0,       0.0f );
  }
}

void PndRhoTupleQA::qaMcList	(	TString	pre,
		RhoCandList &	l,
		RhoTuple *	n,
		int	max = 10000
	)

Definition at line 1102 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), RhoCandList::GetLength(), RhoCandidate::GetTrackNumber(), and max().

Referenced by ana_jpsi(), PndParticleQATask::Exec(), PndTripleAnaTask::Exec(), PndSoftTriggerTask::Exec(), and qaMcList().

{
  if (n==0) return;

  int npart = l.GetLength();
  if (npart>max) npart=max;

  TVector vpart(npart), vpdg(npart), vmoth(npart),
  vp(npart),    vmass(npart), vndau(npart),
  vpx(npart),   vpy(npart),  vpz(npart), ve(npart),
  vtht(npart),  vphi(npart),
  vx(npart),    vy(npart),   vz(npart);

  for (int j=0;j<npart;++j)
  {
    RhoCandidate *moth = l[j]->TheMother();

    vpart(j) = j;
    vpdg(j)  = l[j]->PdgCode();
    vmoth(j) = (moth!=0x0) ? moth->GetTrackNumber() : -1;
    vndau(j) = l[j]->NDaughters();

    vmass(j) = l[j]->Mass();
    vp(j)    = l[j]->P();

    vpx(j)   = l[j]->Px();
    vpy(j)   = l[j]->Py();
    vpz(j)   = l[j]->Pz();
    ve(j)    = l[j]->E();

    vtht(j)  = l[j]->P3().Theta();
    vphi(j)  = l[j]->P3().Phi();

    vx(j)    = l[j]->Pos().X();
    vy(j)    = l[j]->Pos().Y();
    vz(j)    = l[j]->Pos().Z();
  }

  n->Column(pre+"npart", (Int_t) npart);

  n->Column(pre+"part",  vpart,  pre+"npart");
  n->Column(pre+"pdg",   vpdg,   pre+"npart");
  n->Column(pre+"moth",  vmoth,  pre+"npart");
  n->Column(pre+"ndau",  vndau,  pre+"npart");

  n->Column(pre+"m",     vmass,  pre+"npart");
  n->Column(pre+"p",     vp,     pre+"npart");

  n->Column(pre+"px",    vpx,    pre+"npart");
  n->Column(pre+"py",    vpy,    pre+"npart");
  n->Column(pre+"pz",    vpz,    pre+"npart");
  n->Column(pre+"e",     ve,     pre+"npart");

  n->Column(pre+"tht",   vtht,   pre+"npart");
  n->Column(pre+"phi",   vphi,   pre+"npart");

  n->Column(pre+"x",     vx,     pre+"npart");
  n->Column(pre+"y",     vy,     pre+"npart");
  n->Column(pre+"z",     vz,     pre+"npart");
}

void PndRhoTupleQA::qaMcList	(	RhoTuple *	n,
		int	max = 10000
	)

Definition at line 1094 of file PndRhoTupleQA.cxx.

References fAnalysis, PndAnalysis::FillList(), and qaMcList().

{
  RhoCandList mc;
  fAnalysis->FillList(mc,"McTruth");
  qaMcList("",mc,n,max);
}

void PndRhoTupleQA::qaMuo	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 947 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndPidCandidate::GetMuoHits(), PndPidCandidate::GetMuoIndex(), PndPidCandidate::GetMuoIron(), PndPidCandidate::GetMuoModule(), PndPidCandidate::GetMuoMomentumIn(), PndPidCandidate::GetMuoNumberOfLayers(), PndPidCandidate::GetMuoProbability(), PndPidCandidate::GetMuoQuality(), and RhoCandidate::GetRecoCandidate().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"muonlay",    (Int_t)   mic->GetMuoNumberOfLayers(),  0 );;
    n->Column(pre+"muoprob",    (Float_t) mic->GetMuoProbability(),    0.0f );;
    n->Column(pre+"muoqa",    (Float_t) mic->GetMuoQuality() ,    0.0f ); ;
    n->Column(pre+"muoiron",    (Float_t) mic->GetMuoIron() ,      0.0f );  ;
    n->Column(pre+"muopin",    (Float_t) mic->GetMuoMomentumIn(),    0.0f );;
    n->Column(pre+"muomod",   (Int_t)   mic->GetMuoModule(),      0 );  ;
    n->Column(pre+"muohits",    (Int_t)   mic->GetMuoHits(),      0 );  ;
    n->Column(pre+"muoidx",    (Int_t)   mic->GetMuoIndex(),      0 ); ;
  }
}

void PndRhoTupleQA::qaMvd	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 863 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndPidCandidate::GetMvdDEDX(), PndPidCandidate::GetMvdHits(), and RhoCandidate::GetRecoCandidate().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"mvddedx",  (Float_t) mic->GetMvdDEDX(),      0.0f );
    n->Column(pre+"mvdhits",  (Int_t)   mic->GetMvdHits(),      0 );
  }
}

void PndRhoTupleQA::qaP4	(	TString	pre,
		TLorentzVector	c,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 457 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), and f.

Referenced by ana_jpsi(), PndSoftTriggerTask::CreateKs0Cands(), PndTripleAnaTask::DsDs2317Analysis(), PndParticleQATask::Exec(), PndSoftTriggerTask::Exec(), PndSoftTriggerTask::FillGlobalLists(), PndTripleAnaTask::JpsiAnalysis(), qaCand(), PndSoftTriggerTask::TagMode(), PndTripleAnaTask::ThreePiAnalysis(), and tut_ana_d0_qahelper().

{
  if (n==0) return;

  if (!skip)
  {
    n->Column(pre+"px",  (Float_t) c.Px(),     0.0f );
    n->Column(pre+"py",  (Float_t) c.Py(),     0.0f );
    n->Column(pre+"pz",  (Float_t) c.Pz(),     0.0f );
    n->Column(pre+"e",   (Float_t) c.E(),      0.0f );
    n->Column(pre+"p",   (Float_t) c.P(),      0.0f );
    n->Column(pre+"tht", (Float_t) c.Theta(),  0.0f );
    n->Column(pre+"phi", (Float_t) c.Phi(),    0.0f );
    n->Column(pre+"pt",  (Float_t) c.Pt(),     0.0f );
    n->Column(pre+"m",   (Float_t) c.M(),      0.0f );
  }
  else
  {
    n->Column(pre+"px",  (Float_t) -999.,        0.0f );
    n->Column(pre+"py",  (Float_t) -999.,        0.0f );
    n->Column(pre+"pz",  (Float_t) -999.,        0.0f );
    n->Column(pre+"e",   (Float_t) -999.,        0.0f );
    n->Column(pre+"p",   (Float_t) -999.,        0.0f );
    n->Column(pre+"tht", (Float_t) -999.,       0.0f );
    n->Column(pre+"phi", (Float_t) -999.,        0.0f );
    n->Column(pre+"pt",  (Float_t) -999.,        0.0f );
    n->Column(pre+"m",   (Float_t) -999.,        0.0f );
  }
}

void PndRhoTupleQA::qaP4Cms	(	TString	pre,
		TLorentzVector	c,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 488 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), f, and fIniP4.

Referenced by qaComp().

{
  if (n==0) return;

  if (!skip)
  {
    c.Boost(-fIniP4.BoostVector());

    n->Column(pre+"pxcm",  (Float_t) c.Px(),     0.0f );
    n->Column(pre+"pycm",  (Float_t) c.Py(),     0.0f );
    n->Column(pre+"pzcm",  (Float_t) c.Pz(),     0.0f );
    n->Column(pre+"ecm",   (Float_t) c.E(),      0.0f );
    n->Column(pre+"pcm",   (Float_t) c.P(),      0.0f );
    n->Column(pre+"thtcm", (Float_t) c.Theta(),  0.0f );
    n->Column(pre+"phicm", (Float_t) c.Phi(),    0.0f );
  }
  else
  {
    n->Column(pre+"pxcm",  (Float_t) -999.,        0.0f );
    n->Column(pre+"pycm",  (Float_t) -999.,        0.0f );
    n->Column(pre+"pzcm",  (Float_t) -999.,        0.0f );
    n->Column(pre+"ecm",   (Float_t) -999.,        0.0f );
    n->Column(pre+"pcm",   (Float_t) -999.,        0.0f );
    n->Column(pre+"thtcm", (Float_t) -999.,       0.0f );
    n->Column(pre+"phicm", (Float_t) -999.,        0.0f );
  }
}

void PndRhoTupleQA::qaP4Cov	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 518 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), f, and RhoCandidate::P4Cov().

Referenced by qaComp().

{
  if (n==0) return;

  RhoError cov = c->P4Cov();

  if (!skip)
  {
    n->Column(pre+"covpxpx", (Float_t) cov(0,0), 0.0f);
    n->Column(pre+"covpxpy", (Float_t) cov(0,1), 0.0f);
    n->Column(pre+"covpxpz", (Float_t) cov(0,2), 0.0f);
    n->Column(pre+"covpxe",  (Float_t) cov(0,3), 0.0f);
    n->Column(pre+"covpypy", (Float_t) cov(1,1), 0.0f);
    n->Column(pre+"covpypz", (Float_t) cov(1,2), 0.0f);
    n->Column(pre+"covpye",  (Float_t) cov(1,3), 0.0f);
    n->Column(pre+"covpzpz", (Float_t) cov(2,2), 0.0f);
    n->Column(pre+"covpze",  (Float_t) cov(2,3), 0.0f);
    n->Column(pre+"covee",   (Float_t) cov(3,3), 0.0f);
  }
  else
  {
    n->Column(pre+"covpxpx", (Float_t) -999., 0.0f);
    n->Column(pre+"covpxpy", (Float_t) -999., 0.0f);
    n->Column(pre+"covpxpz", (Float_t) -999., 0.0f);
    n->Column(pre+"covpxe",  (Float_t) -999., 0.0f);
    n->Column(pre+"covpypy", (Float_t) -999., 0.0f);
    n->Column(pre+"covpypz", (Float_t) -999., 0.0f);
    n->Column(pre+"covpyee", (Float_t) -999., 0.0f);
    n->Column(pre+"covpzpz", (Float_t) -999., 0.0f);
    n->Column(pre+"covpze",  (Float_t) -999., 0.0f);
    n->Column(pre+"covee",   (Float_t) -999., 0.0f);
  }
}

void PndRhoTupleQA::qaPi0	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 377 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), d0, RhoCandidate::Daughter(), f, fAnalysis, RhoCandidate::GetMcTruth(), PndAnalysis::McTruthMatch(), RhoCandidate::P3(), qaCand(), and qaEmc().

Referenced by PndTripleAnaTask::DsDs2317Analysis(), PndSoftTriggerTask::FillGlobalLists(), and PndTripleAnaTask::ThreePiAnalysis().

{
  if (n==0) return;

  qaCand(pre, c, n);

  RhoCandidate *d0 = c->Daughter(0);
  RhoCandidate *d1 = c->Daughter(1);
  double ang = d0->P3().Angle(d1->P3());

  if (fAnalysis) fAnalysis->McTruthMatch(c);
  RhoCandidate *truth = c->GetMcTruth();

  qaCand(pre, c, n);
  n->Column(pre+"oang",(Float_t) ang,    0.0f);

  qaCand(pre+"d0", d0, n);
  qaEmc(pre+"d0", d0, n);

  qaCand(pre+"d1", d1, n);
  qaEmc(pre+"d1", d1, n);

  if (truth!=0)
  {
    qaCand("t"+pre, truth, n);
    n->Column(pre+"mct", 1.0f, 0.0f);
  }
  else
  {
    qaCand("t"+pre, NULL , n, true);
    n->Column(pre+"mct", 0.0f, 0.0f);
  }
}

void PndRhoTupleQA::qaPid	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 689 of file PndRhoTupleQA.cxx.

References RhoCandidate::Charge(), RhoTuple::Column(), f, fabs(), RhoCandidate::GetPidInfo(), i, and pinf.

Referenced by PndParticleQATask::Exec(), qaComp(), and qaKs0().

{
  if (n==0) return;

  // pinf[0]...pinf[4]: cache P values, pinf[5] = maximum
  double pinf[6] = {-1000.,-1000.,-1000.,-1000.,-1000.,-1000.};
  // index of particle type (e=0 ... p=4) with maximum P
  int bestidx = -1;

  if (fabs(c->Charge())>0.1)
  {
    for (int i=0;i<5;++i)
    {
      pinf[i] = c->GetPidInfo(i);

      if (pinf[i]>pinf[5])
      {
        pinf[5] = pinf[i];
        bestidx = i;
      }
    }
  }

  n->Column(pre+"pide",    (Float_t) pinf[0],    0.0f );
  n->Column(pre+"pidmu",   (Float_t) pinf[1],    0.0f );
  n->Column(pre+"pidpi",  (Float_t) pinf[2],    0.0f );
  n->Column(pre+"pidk",    (Float_t) pinf[3],    0.0f );
  n->Column(pre+"pidp",    (Float_t) pinf[4],    0.0f );

  n->Column(pre+"pidmax",  (Float_t) pinf[5],    0.0f );
  n->Column(pre+"pidbest",(Float_t) bestidx,    0.0f );

}

void PndRhoTupleQA::qaPoca	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 225 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), f, fAnalysis, PndAnalysis::FillList(), fVtxPoca, RhoCandList::GetLength(), RhoVtxPoca::GetPocaVtx(), RhoCandidate::M(), RhoCandidate::P(), and RhoCandList::RemoveFamily().

Referenced by qaComp(), and qaKs0().

{
  if (n==0) return;

  // *** simple vtx finder
  TVector3 vtx, altvtx, primvtx;
  double qavtx = fVtxPoca->GetPocaVtx(vtx, c);

  // *** determine poca of rest of tracks
  RhoCandList l;
  fAnalysis->FillList(l, "Charged");
  if (fVtxPoca->GetPocaVtx(primvtx, l)>998.) primvtx.SetXYZ(0.,0.,0.);

  l.RemoveFamily(c);

  if (l.GetLength()>1) fVtxPoca->GetPocaVtx(altvtx, l);
  else altvtx = primvtx;

  double dist=999.;
  if (altvtx.Mag()>0.) dist = (vtx-altvtx).Mag();

  // *** store QA info
  n->Column(pre+"pocvx",  (Float_t) vtx.X(),   0.0f);
  n->Column(pre+"pocvy",  (Float_t) vtx.Y(),   0.0f);
  n->Column(pre+"pocvz",  (Float_t) vtx.Z(),   0.0f);
  n->Column(pre+"altvx",  (Float_t) altvtx.X(),0.0f);
  n->Column(pre+"altvy",  (Float_t) altvtx.Y(),0.0f);
  n->Column(pre+"altvz",  (Float_t) altvtx.Z(),0.0f);
  n->Column(pre+"pocmag", (Float_t) vtx.Mag(), 0.0f);
  n->Column(pre+"pocqa",  (Float_t) qavtx,     0.0f);
  n->Column(pre+"pocdist",(Float_t) dist,      0.0f);
  n->Column(pre+"pocctau",(Float_t) (dist*c->M()/c->P()), 0.0f); // decay length
}

void PndRhoTupleQA::qaPos	(	TString	pre,
		TVector3	p,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 580 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), and f.

Referenced by qaCand().

{
  if (n==0) return;

  if (!skip)
  {
    n->Column(pre+"x",   (Float_t) p.X(),      0.0f );
    n->Column(pre+"y",   (Float_t) p.Y(),      0.0f );
    n->Column(pre+"z",   (Float_t) p.Z(),      0.0f );
    n->Column(pre+"l",   (Float_t) p.Mag(),    0.0f );
  }
  else
  {
    n->Column(pre+"x",   (Float_t) -999.,        0.0f );
    n->Column(pre+"y",   (Float_t) -999.,        0.0f );
    n->Column(pre+"z",   (Float_t) -999.,        0.0f );
    n->Column(pre+"l",   (Float_t) -999.,        0.0f );
  }
}

void PndRhoTupleQA::qaPRG	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 261 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), f, RhoKalmanVtxFitter::FitVertexFast(), and RhoKalmanVtxFitter::SetSilent().

{
  if (n==0) return;

  RhoKalmanVtxFitter vtxPRG(c);
  vtxPRG.SetSilent();

  // *** PRG vtx finder
  TVector3 vtx(0,0,0);
  TMatrixD vmatrix(3,3);
  double qavtx = vtxPRG.FitVertexFast(vtx, vmatrix, true, 2);//iteration 2=default

  // *** store QA info
  n->Column(pre+"prgvx",  (Float_t) vtx.X(), 0.0f);
  n->Column(pre+"prgvy",  (Float_t) vtx.Y(), 0.0f);
  n->Column(pre+"prgvz",  (Float_t) vtx.Z(), 0.0f);
  n->Column(pre+"prgqa",  (Float_t) qavtx,   0.0f);
}

void PndRhoTupleQA::qaPull	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n,
		bool	skip = false
	)

Definition at line 601 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), f, RhoCandidate::GetMcTruth(), RhoCandidate::P4(), RhoCandidate::P4Cov(), RhoCandidate::Pos(), RhoCandidate::PosCov(), and sqrt().

Referenced by PndParticleQATask::Exec(), and qaComp().

{
  if (!skip)
  {
    RhoCandidate *mct = c->GetMcTruth();
    if(mct) {
      TLorentzVector difp4  = c->P4()  - mct->P4();
      TVector3       difpos = c->Pos() - mct->Pos();
      RhoError       covp4  = c->P4Cov();
      RhoError       covpos = c->PosCov();
      n->Column(pre+"pullpx",   (Float_t) (difp4.Px()/sqrt(covp4(0,0))),               0.0f );
      n->Column(pre+"pullpy",   (Float_t) (difp4.Py()/sqrt(covp4(1,1))),               0.0f );
      n->Column(pre+"pullpz",   (Float_t) (difp4.Pz()/sqrt(covp4(2,2))),               0.0f );
      n->Column(pre+"pulle",    (Float_t) (difp4.E()/sqrt(covp4(3,3))),                0.0f );
      n->Column(pre+"pullx",    (Float_t) (difpos.X()/sqrt(covpos(0,0))),              0.0f );
      n->Column(pre+"pully",    (Float_t) (difpos.Y()/sqrt(covpos(1,1))),              0.0f );
      n->Column(pre+"pullz",    (Float_t) (difpos.Z()/sqrt(covpos(2,2))),              0.0f );
      n->Column(pre+"pullpxpy", (Float_t) (sqrt(difp4.Px()*difp4.Py()/covp4(0,1))),    0.0f );
      n->Column(pre+"pullpxpz", (Float_t) (sqrt(difp4.Px()*difp4.Pz()/covp4(0,2))),    0.0f );
      n->Column(pre+"pullpypz", (Float_t) (sqrt(difp4.Py()*difp4.Pz()/covp4(1,2))),    0.0f );
      n->Column(pre+"pullpxe",  (Float_t) (sqrt(difp4.Px()*difp4.E()/covp4(0,3))),     0.0f );
      n->Column(pre+"pullpye",  (Float_t) (sqrt(difp4.Py()*difp4.E()/covp4(1,3))),     0.0f );
      n->Column(pre+"pullpze",  (Float_t) (sqrt(difp4.Pz()*difp4.E()/covp4(2,3))),     0.0f );
      n->Column(pre+"pullxy",   (Float_t) (sqrt(difpos.X()*difpos.Y()/covpos(0,1))),   0.0f );
      n->Column(pre+"pullxz",   (Float_t) (sqrt(difpos.Y()*difpos.Z()/covpos(0,2))),   0.0f );
      n->Column(pre+"pullyz",   (Float_t) (sqrt(difpos.Y()*difpos.Z()/covpos(1,2))),   0.0f );
    } else {skip=true;}
  }

  if (skip)
  {
      n->Column(pre+"pullpx",   (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpy",   (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpz",   (Float_t) -999.,      0.0f );
      n->Column(pre+"pulle",    (Float_t) -999.,      0.0f );
      n->Column(pre+"pullx",    (Float_t) -999.,      0.0f );
      n->Column(pre+"pully",    (Float_t) -999.,      0.0f );
      n->Column(pre+"pullz",    (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpxpy", (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpxpz", (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpypz", (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpxe",  (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpye",  (Float_t) -999.,      0.0f );
      n->Column(pre+"pullpze",  (Float_t) -999.,      0.0f );
      n->Column(pre+"pullxy",   (Float_t) -999.,      0.0f );
      n->Column(pre+"pullxz",   (Float_t) -999.,      0.0f );
      n->Column(pre+"pullyz",   (Float_t) -999.,      0.0f );
 }
}

void PndRhoTupleQA::qaRecoFull	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 1043 of file PndRhoTupleQA.cxx.

References qaDrc(), qaDsc(), qaEmc(), qaGem(), qaMuo(), qaMvd(), qaRich(), qaStt(), qaTof(), and qaTrk().

Referenced by qaRecoFullTree().

{
  if (n==0) return;

  qaEmc( pre, c, n );
  qaMvd( pre, c, n );
  qaStt( pre, c, n );
  qaGem( pre, c, n );
  qaDrc( pre, c, n );
  qaDsc( pre, c, n );
  qaRich( pre, c, n );
  qaTof( pre, c, n );
  qaMuo( pre, c, n );
  qaTrk( pre, c, n);
}

void PndRhoTupleQA::qaRecoFullTree	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 1060 of file PndRhoTupleQA.cxx.

References RhoCandidate::Daughter(), i, name, RhoCandidate::NDaughters(), qaRecoFull(), and TString.

{
  if (n==0) return;
  if (c==0) return;
  int nd = c->NDaughters();
  if (nd>0)  {
    for (int i=0; i<nd; ++i)
    {
      RhoCandidate *dau = c->Daughter(i);
      TString name=TString::Format("%sd%d",pre.Data(),i);
      qaRecoFullTree(name,dau,n);
    }
  } else {
    qaRecoFull(pre,c,n);
  }
}

void PndRhoTupleQA::qaRecoShort	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 986 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), PndPidCandidate::GetDiscNumberOfPhotons(), PndPidCandidate::GetDiscThetaC(), PndPidCandidate::GetDrcNumberOfPhotons(), PndPidCandidate::GetDrcThetaC(), PndPidCandidate::GetEmcCalEnergy(), PndPidCandidate::GetEmcNumberOfBumps(), PndPidCandidate::GetEmcNumberOfCrystals(), PndPidCandidate::GetGemHits(), PndPidCandidate::GetMuoIron(), PndPidCandidate::GetMuoNumberOfLayers(), PndPidCandidate::GetMvdDEDX(), PndPidCandidate::GetMvdHits(), RhoCandidate::GetRecoCandidate(), PndPidCandidate::GetRichNumberOfPhotons(), PndPidCandidate::GetRichThetaC(), PndPidCandidate::GetSttHits(), PndPidCandidate::GetSttMeanDEDX(), PndPidCandidate::GetTofBeta(), and PndPidCandidate::GetTofM2().

Referenced by qaRecoShortTree().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"emcecal",  (Float_t) mic->GetEmcCalEnergy(),      0.0f );
    n->Column(pre+"emcnx",    (Int_t)   mic->GetEmcNumberOfCrystals(),  0 );
    n->Column(pre+"emcnb",    (Int_t)   mic->GetEmcNumberOfBumps(),     0 );

    n->Column(pre+"gemnhits",  (Int_t) mic->GetGemHits(),          0 );

    n->Column(pre+"mvddedx",  (Float_t) mic->GetMvdDEDX(),          0.0f );
    n->Column(pre+"mvdhits",  (Int_t)   mic->GetMvdHits(),          0 );

    n->Column(pre+"sttdedx",  (Float_t) mic->GetSttMeanDEDX(),        0.0f );
    n->Column(pre+"stthits",  (Int_t)   mic->GetSttHits(),          0 );

    n->Column(pre+"drcthtc",    (Float_t) mic->GetDrcThetaC(),      0.0f );
    n->Column(pre+"drcnphot",    (Int_t)   mic->GetDrcNumberOfPhotons(),  0 );

    n->Column(pre+"dscthtc",    (Float_t) mic->GetDiscThetaC(),      0.0f );
    n->Column(pre+"dscnphot",    (Int_t)   mic->GetDiscNumberOfPhotons(),0 );

    n->Column(pre+"richthtc",    (Float_t) mic->GetRichThetaC(),      0.0f );
    n->Column(pre+"richnphot",  (Int_t)   mic->GetRichNumberOfPhotons(),0 );

    n->Column(pre+"muonlay",    (Int_t)   mic->GetMuoNumberOfLayers(),  0 );;
    n->Column(pre+"muoiron",    (Float_t) mic->GetMuoIron() ,      0.0f );  ;

    n->Column(pre+"tofm2",    (Float_t) mic->GetTofM2(),          0.0f );
    n->Column(pre+"tofbeta",    (Float_t) mic->GetTofBeta(),        0.0f );
  }
}

void PndRhoTupleQA::qaRecoShortTree	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 1024 of file PndRhoTupleQA.cxx.

References RhoCandidate::Daughter(), i, name, RhoCandidate::NDaughters(), qaRecoShort(), and TString.

{
  if (n==0) return;
  if (c==0) return;
  int nd = c->NDaughters();
  if (nd>0)  {
    for (int i=0; i<nd; ++i)
    {
      RhoCandidate *dau = c->Daughter(i);
      TString name=TString::Format("%sd%d",pre.Data(),i);
      qaRecoShortTree(name,dau,n);
    }
  } else {
    qaRecoShort(pre,c,n);
  }
}

void PndRhoTupleQA::qaRich	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 929 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), RhoCandidate::GetRecoCandidate(), PndPidCandidate::GetRichIndex(), PndPidCandidate::GetRichNumberOfPhotons(), PndPidCandidate::GetRichQuality(), PndPidCandidate::GetRichThetaC(), and PndPidCandidate::GetRichThetaCErr().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"richthtc",    (Float_t) mic->GetRichThetaC(),      0.0f );
    n->Column(pre+"richdthtc",  (Float_t) mic->GetRichThetaCErr(),    0.0f );
    n->Column(pre+"richqa",    (Float_t) mic->GetRichQuality(),    0.0f );
    n->Column(pre+"richnphot",  (Int_t)   mic->GetRichNumberOfPhotons(),0 );
    n->Column(pre+"richidx",    (Int_t)    mic->GetRichIndex(),      0 );
  }
}

void PndRhoTupleQA::qaStt	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 878 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), RhoCandidate::GetRecoCandidate(), PndPidCandidate::GetSttHits(), and PndPidCandidate::GetSttMeanDEDX().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"sttdedx",  (Float_t) mic->GetSttMeanDEDX(),    0.0f );
    n->Column(pre+"stthits",  (Int_t)   mic->GetSttHits(),      0 );
  }
}

void PndRhoTupleQA::qaTof	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 968 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), RhoCandidate::GetRecoCandidate(), PndPidCandidate::GetTofBeta(), PndPidCandidate::GetTofIndex(), PndPidCandidate::GetTofM2(), PndPidCandidate::GetTofQuality(), PndPidCandidate::GetTofStopTime(), and PndPidCandidate::GetTofTrackLength().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"toftime",   (Float_t) mic->GetTofStopTime(),  0.0f );
    n->Column(pre+"tofm2",    (Float_t) mic->GetTofM2(),      0.0f );
    n->Column(pre+"toflen",    (Float_t) mic->GetTofTrackLength(),  0.0f );
    n->Column(pre+"tofqa",    (Float_t) mic->GetTofQuality(),    0.0f );
    n->Column(pre+"tofidx",    (Int_t)   mic->GetTofIndex(),    0 );
    n->Column(pre+"tofbeta",    (Float_t) mic->GetTofBeta(),    0.0f );
  }
}

void PndRhoTupleQA::qaTrk	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 652 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), f, PndPidCandidate::GetChiSquared(), PndPidCandidate::GetDegreesOfFreedom(), PndPidCandidate::GetFitStatus(), PndPidCandidate::GetMomentum(), PndPidCandidate::GetPosition(), RhoCandidate::GetRecoCandidate(), PndPidCandidate::GetTrackBranch(), and PndPidCandidate::GetTrackIndex().

Referenced by PndParticleQATask::Exec(), and qaRecoFull().

{
  if (n==0) return;

  PndPidCandidate *mic = (PndPidCandidate*)c->GetRecoCandidate();

  if (mic)
  {
    n->Column(pre+"trkpx",      (Float_t) mic->GetMomentum().X(),     0.0f );
    n->Column(pre+"trkpy",      (Float_t) mic->GetMomentum().Y(),     0.0f );
    n->Column(pre+"trkpz",      (Float_t) mic->GetMomentum().Z(),     0.0f );
    n->Column(pre+"trkvx",      (Float_t) mic->GetPosition().X(),     0.0f );
    n->Column(pre+"trkvy",      (Float_t) mic->GetPosition().Y(),     0.0f );
    n->Column(pre+"trkvz",      (Float_t) mic->GetPosition().Z(),     0.0f );
    n->Column(pre+"trkdof",     (Int_t)   mic->GetDegreesOfFreedom(), 0 );
    n->Column(pre+"trkstat",    (Int_t)   mic->GetFitStatus(),        0 );
    n->Column(pre+"trkchi2",    (Float_t) mic->GetChiSquared(),       0.0f );
    n->Column(pre+"trkidx",     (Int_t)   mic->GetTrackIndex(),       0 );
    n->Column(pre+"trkbranch",  (Int_t)   mic->GetTrackBranch() ,     0 );
  } else {
    n->Column(pre+"trkpx",      (Float_t) -999.,  0.0f );
    n->Column(pre+"trkpy",      (Float_t) -999.,  0.0f );
    n->Column(pre+"trkpz",      (Float_t) -999.,  0.0f );
    n->Column(pre+"trkvx",      (Float_t) -999.,  0.0f );
    n->Column(pre+"trkvy",      (Float_t) -999.,  0.0f );
    n->Column(pre+"trkvz",      (Float_t) -999.,  0.0f );
    n->Column(pre+"trkdof",     (Int_t)   -999.,  0 );
    n->Column(pre+"trkstat",    (Int_t)   -999.,  0 );
    n->Column(pre+"trkchi2",    (Float_t) -999.,  0.0f );
    n->Column(pre+"trkidx",     (Int_t)   -999.,  0 );
    n->Column(pre+"trkbranch",  (Int_t)   -999.,  0 );
  }

}

void PndRhoTupleQA::qaVtx	(	TString	pre,
		RhoCandidate *	c,
		RhoTuple *	n
	)

Definition at line 771 of file PndRhoTupleQA.cxx.

References RhoTuple::Column(), d, RhoCandidate::Daughter(), RhoCandidate::DecayVtx(), f, fAnalysis, PndAnalysis::FillList(), fVtxPoca, RhoVtxPoca::GetPocaVtx(), RhoCandidate::M(), RhoCandidate::P(), RhoCandidate::Pos(), and v.

Referenced by qaComp(), and qaKs0().

{
  if (n==0) return;

   RhoCandidate *d = c->Daughter(0);
  TVector3 v = c->DecayVtx();
  if (d)
  {
    if (v.X()==0 && v.Y()==0 && v.Z()==0)
      v = d->Pos();
    //     TLorentzVector d_cms = d->P4();
    //     d_cms.Boost(-(c->P4().BoostVector()));

    Float_t ctau = v.Mag()*c->M()/c->P();
    /*    Float_t dec  = d_cms.Vect().Angle(c->P3());
     Float_t cdec = cos(dec);*/

    // if primary Vertex available, compute ctau relative to that one
    Float_t ctaud = -999.;
    TVector3 primvtx;

    // *** determine poca of all charged tracks as primary vertex
    RhoCandList l;
    fAnalysis->FillList(l, "Charged");
    if (fVtxPoca->GetPocaVtx(primvtx, l)>998.) primvtx.SetXYZ(0.,0.,0.);

    if (primvtx.Mag()>0) ctaud = (v-primvtx).Mag()*c->M()/c->P();

    n->Column(pre+"vx",    (Float_t) v.X(),    0.0f );
    n->Column(pre+"vy",    (Float_t) v.Y(),    0.0f );
    n->Column(pre+"vz",    (Float_t) v.Z(),    0.0f );
    n->Column(pre+"len",   (Float_t) v.Mag(),    0.0f );
    n->Column(pre+"ctau",  (Float_t) ctau,      0.0f );
    n->Column(pre+"ctaud",  (Float_t) ctaud,    0.0f );
    /*    n->Column(pre+"decang", (Float_t) dec,      0.0f );
     n->Column(pre+"cdecang",(Float_t) cdec,      0.0f );*/
  }
  else
  {
    n->Column(pre+"vx",    (Float_t) -999.0,    0.0f );
    n->Column(pre+"vy",    (Float_t) -999.0,    0.0f );
    n->Column(pre+"vz",    (Float_t) -999.0,    0.0f );
    n->Column(pre+"len",   (Float_t) -999.0,    0.0f );
    n->Column(pre+"ctau",    (Float_t) -999.0,    0.0f );
    n->Column(pre+"ctaud",  (Float_t) -999.0,    0.0f );
    /*    n->Column(pre+"decang", (Float_t) -999.0,    0.0f );
     n->Column(pre+"cdecang",(Float_t) -999.0,    0.0f );*/
  }
}

void PndRhoTupleQA::SetIniP4 ( TLorentzVector  _ini )

inline

Definition at line 81 of file PndRhoTupleQA.h.

References fIniP4.

{fIniP4 = _ini;}

Member Data Documentation

PndAnalysis* PndRhoTupleQA::fAnalysis

private

Definition at line 88 of file PndRhoTupleQA.h.

Referenced by PndRhoTupleQA(), qaComp(), qaKs0(), qaMcList(), qaPi0(), qaPoca(), and qaVtx().

RhoCandidate* PndRhoTupleQA::fDummyCand

private

Definition at line 86 of file PndRhoTupleQA.h.

Referenced by PndRhoTupleQA(), qaMc(), and ~PndRhoTupleQA().

TLorentzVector PndRhoTupleQA::fIniP4

private

Definition at line 90 of file PndRhoTupleQA.h.

Referenced by GetIniP4(), PndRhoTupleQA(), qaP4Cms(), and SetIniP4().

RhoVtxPoca* PndRhoTupleQA::fVtxPoca

private

Definition at line 87 of file PndRhoTupleQA.h.

Referenced by PndRhoTupleQA(), qaPoca(), qaVtx(), and ~PndRhoTupleQA().

---

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

• PndRhoTupleQA.h
• PndRhoTupleQA.cxx