PndPmtPoormantracks Class Reference

#include <PndPmtPoormantracks.h>

Inheritance diagram for PndPmtPoormantracks:

Public Member Functions
	PndPmtPoormantracks ()
	~PndPmtPoormantracks ()
void	SetSigV (double a)
void	SetSigP (double a)
void	SetPtRange (double min, double max)
void	SetDTheta (double a)
void	SetPid (int a)
void	SetNumTrk (int a)
void	SetSeed (int a)
void	SetVerbose (int a)
void	EnergyCorrelations ()

Protected Member Functions
TVector3	RollVertexBox (double widx, double widy, double widz)
void	SmearVertex (TVector3 &vertex)
TVector3	RollMomentumBox (const TVector3 &vtx, double dtheta, double ptmin, double ptmax)
void	SmearMomentum (TVector3 &momentum)
void	PoorManTracks ()
void	EnergyCorrelations (RhoError &covP4, TLorentzVector p4)

Protected Attributes
TClonesArray *	fMcCands
TClonesArray *	fCands
double	fSigVx
double	fSigVy
double	fSigVz
double	fSigPx
double	fSigPy
double	fSigPz
double	fCharge
double	fDtheta
int	fPID
TParticlePDG *	fPDG
int	fNumTrk
int	niter
int	fNDF
int	fSeed
TVector3	fVertex
double	fWidx
double	fWidy
double	fWidz
double	fPtMin
double	fPtMax
int	fVerbose

Detailed Description

Definition at line 12 of file PndPmtPoormantracks.h.

Constructor & Destructor Documentation

PndPmtPoormantracks::PndPmtPoormantracks

(

)

Definition at line 23 of file PndPmtPoormantracks.cxx.

References fCands, fCharge, fDtheta, fMcCands, fNumTrk, fPDG, fPID, fPtMax, fPtMin, fSeed, fSigPx, fSigPy, fSigPz, fSigVx, fSigVy, fSigVz, and fVerbose.

{
  fSigVx=0.005; //[cm]
  fSigVy=0.005; //[cm]
  fSigVz=0.005; //[cm]
  fSigPx=0.01; //[GeV/c]
  fSigPy=0.01; //[GeV/c]
  fSigPz=0.01; //[GeV/c]
  fPtMin=0.9;//[GeV/c]
  fPtMax=1.1;//[GeV/c]
  fCharge=1.;
  fDtheta = 1.; //0.3 [pi]
  fPID=211;// 211 pion, 22 muon, 2212 proton, 321 kaon
  fPDG = 0;
  fNumTrk=2;
  fSeed=936650;//
  fVerbose=0;
  fMcCands = new TClonesArray("RhoCandidate");
  fCands = new TClonesArray("RhoCandidate");
}

PndPmtPoormantracks::~PndPmtPoormantracks

(

)

Definition at line 44 of file PndPmtPoormantracks.cxx.

{
}

Member Function Documentation

void PndPmtPoormantracks::EnergyCorrelations

(

)

Referenced by PoorManTracks().

void PndPmtPoormantracks::EnergyCorrelations ( RhoError &  covP4, TLorentzVector  p4  )

protected

Definition at line 147 of file PndPmtPoormantracks.cxx.

{
  double invE=1/p4.E();
    covP4[3][3] = (p4.X()*p4.X()*covP4[0][0]+p4.Y()*p4.Y()*covP4[1][1]+p4.Z()*p4.Z()*covP4[2][2])*invE*invE;
    covP4[0][3] = covP4[3][0] = p4.X()*covP4[0][0]*invE;
    covP4[1][3] = covP4[3][1] = p4.Y()*covP4[1][1]*invE;
    covP4[2][3] = covP4[3][2] = p4.Z()*covP4[2][2]*invE;
//    covP4[0][3] = covP4[3][0] = (p4.X()*covP4[0][0]+p4.Y()*covP4[0][1]+p4.Z()*covP4[0][2])*invE;
//    covP4[1][3] = covP4[3][1] = (p4.X()*covP4[1][0]+p4.Y()*covP4[1][1]+p4.Z()*covP4[1][2])*invE;
//    covP4[2][3] = covP4[3][2] = (p4.X()*covP4[2][0]+p4.Y()*covP4[2][1]+p4.Z()*covP4[2][2])*invE;

}

void PndPmtPoormantracks::PoorManTracks ( )

protected

Definition at line 95 of file PndPmtPoormantracks.cxx.

References EnergyCorrelations(), fCands, fCharge, fDtheta, fMcCands, fNumTrk, fPDG, fPID, fPtMax, fPtMin, fSigPx, fSigPy, fSigPz, fSigVx, fSigVy, fSigVz, fVerbose, fVertex, RhoCandidate::GetPosition(), p, RhoCandidate::P3Cov(), RhoCandidate::P4(), Pi, RollMomentumBox(), RhoCandidate::SetCov7(), RhoCandidate::SetMcTruth(), RhoCandidate::SetPos(), SmearMomentum(), and SmearVertex().

Referenced by PndPmtTask::Exec().

{
  if(fVerbose>1) cout<<" --------------------------------------------------- "<<endl;
  if(fVerbose>0) cout<<"Poor man tracks."<<endl;
  fMcCands->Clear("C");
  fCands->Clear("C");

  RhoError covPos(3), covP4(4);
  covPos[0][0]=fSigVx*fSigVx;
  covPos[1][1]=fSigVy*fSigVy;
  covPos[2][2]=fSigVz*fSigVz;
  covP4[0][0]=fSigPx*fSigPx;
  covP4[1][1]=fSigPy*fSigPy;
  covP4[2][2]=fSigPz*fSigPz;

  for (int iTr=0; iTr<fNumTrk; iTr++)
  {
    fPDG = TDatabasePDG::Instance()->GetParticle(fPID);
    TVector3 momentum = RollMomentumBox(fVertex,fDtheta*TMath::Pi(),fPtMin,fPtMax);

    int sizeMc = fMcCands->GetEntriesFast();
    RhoCandidate *pmc=new ((*fMcCands)[sizeMc]) RhoCandidate(momentum,fPDG);
    pmc->SetPos(fVertex);

    // Smeared "reco"
    int size = fCands->GetEntriesFast();
    TVector3 vertexMeasure = fVertex;
    SmearVertex(vertexMeasure);
    TVector3 momMeasure = momentum;
    SmearMomentum(momMeasure);


    RhoCandidate *p=new ((*fCands)[size]) RhoCandidate(momMeasure,fPDG);//TODO charge
    p->SetMcTruth(pmc);
    p->SetPos(vertexMeasure);
    EnergyCorrelations(covP4,p->P4());
    p->SetCov7(covPos,covP4);

    if(fVerbose>1) cout<<" ---------\n pmc, p, pmcpos, ppos, fVertex"<<endl;
    if(fVerbose>1) std::cout<<*pmc<<std::endl;
    if(fVerbose>1) std::cout<<*p<<std::endl;
    if(fVerbose>1) p->P3Cov().Print();
    if(fVerbose>1) pmc->GetPosition().Print();
    if(fVerbose>1) p->GetPosition().Print();
    if(fVerbose>1) fVertex.Print();
    // print helix parameters and cov together with RhoCandidate V-P4 and cov7

    // alternate particle type
    fCharge*=-1;
    fPID*=-1;
  }
}

TVector3 PndPmtPoormantracks::RollMomentumBox ( const TVector3 &  vtx, double  dtheta, double  ptmin, double  ptmax  )

protected

Definition at line 72 of file PndPmtPoormantracks.cxx.

References fVerbose, phi, pt(), theta, and CAMath::TwoPi().

Referenced by PoorManTracks().

{
  // allow whole phi range
  double phi=TMath::TwoPi()*gRandom->Rndm();
  // allow theta at vertex pointing theta +- dtheta
  double theta = vtx.Theta() + dtheta*2*(0.5-gRandom->Rndm());
  // allow transverse momentum in [ptmin,ptmax]
  double pt = ptmin + (ptmax-ptmin)*gRandom->Rndm();
  TVector3 momentum;
  momentum.SetPtThetaPhi(pt,theta,phi);
  if(fVerbose>2) cout<<"RollMomentumBox: ("<<momentum.x()<<","<<momentum.y()<<","<<momentum.z()<<")"<<endl;
  return momentum;
}

TVector3 PndPmtPoormantracks::RollVertexBox ( double  widx, double  widy, double  widz  )

protected

Definition at line 51 of file PndPmtPoormantracks.cxx.

References fVerbose.

Referenced by PndPmtTask::Exec().

{
  double rnd1=widx*2*(0.5-gRandom->Rndm());
  double rnd2=widy*2*(0.5-gRandom->Rndm());
  double rnd3=widz*2*(0.5-gRandom->Rndm());
  TVector3 vertex(rnd1,rnd2,rnd3);
  if(fVerbose>2) cout<<"RollVertexBox: ("<<vertex.x()<<","<<vertex.y()<<","<<vertex.z()<<")"<<endl;
  return vertex;
}

void PndPmtPoormantracks::SetDTheta ( double  a )

inline

Definition at line 20 of file PndPmtPoormantracks.h.

References a, and fDtheta.

{fDtheta=a;};

void PndPmtPoormantracks::SetNumTrk ( int  a )

inline

Definition at line 22 of file PndPmtPoormantracks.h.

References a, and fNumTrk.

{fNumTrk=a;};

void PndPmtPoormantracks::SetPid ( int  a )

inline

Definition at line 21 of file PndPmtPoormantracks.h.

References a, and fPID.

{fPID=a;};

void PndPmtPoormantracks::SetPtRange ( double  min, double  max  )

inline

Definition at line 19 of file PndPmtPoormantracks.h.

References fPtMax, fPtMin, max(), and min().

{fPtMin=min;fPtMax=max;};

void PndPmtPoormantracks::SetSeed ( int  a )

inline

Definition at line 23 of file PndPmtPoormantracks.h.

References a, and fSeed.

{fSeed=a;};

void PndPmtPoormantracks::SetSigP ( double  a )

inline

Definition at line 18 of file PndPmtPoormantracks.h.

References a, fSigPx, fSigPy, and fSigPz.

{fSigPx=fSigPy=fSigPz=a;};

void PndPmtPoormantracks::SetSigV ( double  a )

inline

Definition at line 17 of file PndPmtPoormantracks.h.

References a, fSigVx, fSigVy, and fSigVz.

{fSigVx=fSigVy=fSigVz=a;};

void PndPmtPoormantracks::SetVerbose ( int  a )

inline

Definition at line 24 of file PndPmtPoormantracks.h.

References a, and fVerbose.

{fVerbose=a;};

void PndPmtPoormantracks::SmearMomentum ( TVector3 &  momentum )

protected

Definition at line 86 of file PndPmtPoormantracks.cxx.

References fSigPx, fSigPy, and fSigPz.

Referenced by PoorManTracks().

{
  double rnd1=gRandom->Gaus(momentum.Px(),fSigPx);
  double rnd2=gRandom->Gaus(momentum.Py(),fSigPy);
  double rnd3=gRandom->Gaus(momentum.Pz(),fSigPz);
  momentum.SetXYZ(rnd1,rnd2,rnd3);
  return;
}

void PndPmtPoormantracks::SmearVertex ( TVector3 &  vertex )

protected

Definition at line 61 of file PndPmtPoormantracks.cxx.

References fSigVx, fSigVy, fSigVz, and fVerbose.

Referenced by PoorManTracks().

{
  double rnd1=gRandom->Gaus(vertex.x(),fSigVx);
  double rnd2=gRandom->Gaus(vertex.y(),fSigVy);
  double rnd3=gRandom->Gaus(vertex.z(),fSigVz);
  if(fVerbose>2) cout<<"SmearVertex (B): ("<<vertex.x()<<","<<vertex.y()<<","<<vertex.z()<<")"<<endl;
  vertex.SetXYZ(rnd1,rnd2,rnd3);
  if(fVerbose>2) cout<<"SmearVertex (A): ("<<vertex.x()<<","<<vertex.y()<<","<<vertex.z()<<")"<<endl;
  return;
}

Member Data Documentation

TClonesArray* PndPmtPoormantracks::fCands

protected

Definition at line 36 of file PndPmtPoormantracks.h.

Referenced by PndPmtTask::Exec(), PndPmtPoormantracks(), and PoorManTracks().

double PndPmtPoormantracks::fCharge

protected

Definition at line 46 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), and PoorManTracks().

double PndPmtPoormantracks::fDtheta

protected

Definition at line 47 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), and SetDTheta().

TClonesArray* PndPmtPoormantracks::fMcCands

protected

Definition at line 35 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), and PoorManTracks().

int PndPmtPoormantracks::fNDF

protected

Definition at line 52 of file PndPmtPoormantracks.h.

int PndPmtPoormantracks::fNumTrk

protected

Definition at line 50 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), and SetNumTrk().

TParticlePDG* PndPmtPoormantracks::fPDG

protected

Definition at line 49 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), and PoorManTracks().

int PndPmtPoormantracks::fPID

protected

Definition at line 48 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), and SetPid().

double PndPmtPoormantracks::fPtMax

protected

Definition at line 56 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), and SetPtRange().

double PndPmtPoormantracks::fPtMin

protected

Definition at line 56 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), and SetPtRange().

int PndPmtPoormantracks::fSeed

protected

Definition at line 53 of file PndPmtPoormantracks.h.

Referenced by PndPmtTask::Init(), PndPmtPoormantracks(), and SetSeed().

double PndPmtPoormantracks::fSigPx

protected

Definition at line 43 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), SetSigP(), and SmearMomentum().

double PndPmtPoormantracks::fSigPy

protected

Definition at line 44 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), SetSigP(), and SmearMomentum().

double PndPmtPoormantracks::fSigPz

protected

Definition at line 45 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), SetSigP(), and SmearMomentum().

double PndPmtPoormantracks::fSigVx

protected

Definition at line 40 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), SetSigV(), and SmearVertex().

double PndPmtPoormantracks::fSigVy

protected

Definition at line 41 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), SetSigV(), and SmearVertex().

double PndPmtPoormantracks::fSigVz

protected

Definition at line 42 of file PndPmtPoormantracks.h.

Referenced by PndPmtPoormantracks(), PoorManTracks(), SetSigV(), and SmearVertex().

int PndPmtPoormantracks::fVerbose

protected

Definition at line 57 of file PndPmtPoormantracks.h.

Referenced by PndPmtTask::Exec(), PndPmtPoormantracks(), PoorManTracks(), RollMomentumBox(), RollVertexBox(), SetVerbose(), and SmearVertex().

TVector3 PndPmtPoormantracks::fVertex

protected

Definition at line 54 of file PndPmtPoormantracks.h.

Referenced by PndPmtTask::Exec(), and PoorManTracks().

double PndPmtPoormantracks::fWidx

protected

Definition at line 55 of file PndPmtPoormantracks.h.

Referenced by PndPmtTask::Exec(), and PndPmtTask::Init().

double PndPmtPoormantracks::fWidy

protected

Definition at line 55 of file PndPmtPoormantracks.h.

Referenced by PndPmtTask::Exec(), and PndPmtTask::Init().

double PndPmtPoormantracks::fWidz

protected

Definition at line 55 of file PndPmtPoormantracks.h.

Referenced by PndPmtTask::Exec(), and PndPmtTask::Init().

int PndPmtPoormantracks::niter

protected

Definition at line 51 of file PndPmtPoormantracks.h.

---

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

• PndPmtPoormantracks.h
• PndPmtPoormantracks.cxx