drawdal.C File Reference

Go to the source code of this file.

Functions
Double_t	simin (Int_t i, Double_t m0, Double_t m1, Double_t m2, Double_t m3)
Double_t	simax (Int_t i, Double_t m0, Double_t m1, Double_t m2, Double_t m3)
Double_t	lambda (Double_t x, Double_t y, Double_t z)
Double_t	s2min (Double_t s1, Double_t m0, Double_t m1, Double_t m2, Double_t m3)
Double_t	s2max (Double_t s1, Double_t m0, Double_t m1, Double_t m2, Double_t m3)
void	CreateGraph (TGraph *g, Int_t n, Double_t m, Double_t m1, Double_t m2, Double_t m3)
	drawdal (double m0, double m1, double m2, double m3, TString opt="Lsame", int n=200)

Function Documentation

void CreateGraph	(	TGraph *	g,
		Int_t	n,
		Double_t	m,
		Double_t	m1,
		Double_t	m2,
		Double_t	m3
	)

Definition at line 74 of file drawdal.C.

References Double_t, i, n, s2max(), s2min(), simax(), simin(), and x.

Referenced by drawdal().

{
        Double_t s1minl = simin(1, m, m1, m2, m3);
        Double_t s1maxl = simax(1, m, m1, m2, m3);
        
        Double_t binw = (s1maxl-s1minl)/(Double_t)(n);
        
        Int_t i=0;
        
        Double_t ymin,ymax,x;
        
        for (i=0; i<n; i++)
        {
                x = i*binw + s1minl;
                if (x<=s1minl) x=s1minl+binw/20;
                if (x>=s1maxl) x=s1maxl-binw/20;
                
                ymin = s2min(x, m, m1, m2, m3);
                ymax = s2max(x, m, m1, m2, m3);
                
                g->SetPoint(i,x,ymin);
                g->SetPoint(n*2-i,x,ymax);
        }

        g->SetPoint(n,s1maxl,s2min(s1maxl-binw/20,m,m1,m2,m3));
        g->GetPoint(0,x,ymin);
        g->SetPoint(n*2,x,ymin);
        
}

drawdal	(	double	m0,
		double	m1,
		double	m2,
		double	m3,
		TString	opt = "Lsame",
		int	n = 200
	)

Definition at line 104 of file drawdal.C.

References CreateGraph(), g, and n.

{
        TGraph *g=new TGraph(n);
        g->SetLineColor(2);
        g->SetLineStyle(2);
        g->SetLineWidth(2);
        CreateGraph(g,n,m0,m1,m2,m3);
        g->Draw(opt);
}

Double_t lambda	(	Double_t	x,
		Double_t	y,
		Double_t	z
	)

Definition at line 48 of file drawdal.C.

References x.

Referenced by PndDiscTaskReconstruction::Exec(), PndDrc::Initialize(), PndDrc::ProcessHits(), DrawHits::ProcessPhotonMC(), PndDrcHitProducerReal::ProcessPhotonPoint(), DecayTreeFitter::InternalParticle::projectKineConstraint(), propagate(), s2max(), s2min(), KFParticleBase::SetMassConstraint(), and KFParticleBaseSIMD::SetMassConstraint().

{
  return x*x + y*y + z*z - 2.*x*y - 2.*y*z - 2.*z*x; 
}

Double_t s2max	(	Double_t	s1,
		Double_t	m0,
		Double_t	m1,
		Double_t	m2,
		Double_t	m3
	)

Definition at line 63 of file drawdal.C.

References Double_t, lambda(), s, and sqrt().

Referenced by CreateGraph().

{
        Double_t s      = m0*m0;
        Double_t lamterm = sqrt( lambda(s1,s,m1*m1) ) * sqrt( lambda(s1, m2*m2, m3*m3) );
         
        Double_t result  = m1*m1 + m3*m3 + ( (s-s1-m1*m1)*(s1-m2*m2+m3*m3) + lamterm )/(2.*s1);
        
        return result;
}

Double_t s2min	(	Double_t	s1,
		Double_t	m0,
		Double_t	m1,
		Double_t	m2,
		Double_t	m3
	)

Definition at line 53 of file drawdal.C.

References Double_t, lambda(), s, and sqrt().

Referenced by CreateGraph().

{
        Double_t s      = m0*m0;
        Double_t lamterm = sqrt( lambda(s1,s,m1*m1) ) * sqrt( lambda(s1, m2*m2, m3*m3) );
         
        Double_t result  = m1*m1 + m3*m3 + ( (s-s1-m1*m1)*(s1-m2*m2+m3*m3) - lamterm )/(2.*s1);
        
        return result;
}

Double_t simax	(	Int_t	i,
		Double_t	m0,
		Double_t	m1,
		Double_t	m2,
		Double_t	m3
	)

Definition at line 32 of file drawdal.C.

References Double_t, and m2().

Referenced by CreateGraph().

{
        Double_t result = 0.;
        
        switch (i)
        {
        case 1:  result = (m0-m1)*(m0-m1);
                break;
        case 2:  result = (m0-m2)*(m0-m2);
                break;
        case 3:  result = (m0-m3)*(m0-m3);
                break;
        }
        return result;
}

Double_t simin	(	Int_t	i,
		Double_t	m0,
		Double_t	m1,
		Double_t	m2,
		Double_t	m3
	)

Definition at line 16 of file drawdal.C.

References Double_t, and m2().

Referenced by CreateGraph().

{
        Double_t result = 0.;
        
        switch (i)
        {
        case 1:  result = (m2+m3)*(m2+m3);
                break;
        case 2:  result = (m1+m3)*(m1+m3);
                break;
        case 3:  result = (m1+m2)*(m1+m2);
                break;
        }
        return result;
}