PndPmtTask Class Reference

#include <PndPmtTask.h>

Inheritance diagram for PndPmtTask:

Public Member Functions
	PndPmtTask ()
	~PndPmtTask ()
void	Init ()
void	Exec (Option_t *opt)
void	Finish ()
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
void	EnergyCorrelations (RhoError &covP4, TLorentzVector p4)
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 ()

Protected Attributes
TList *	fHistoList
TH1F *	fHVtxDiffX
TH1F *	fHVtxDiffY
TH1F *	fHVtxDiffZ
TH1F *	fHVtxDiffPX
TH1F *	fHVtxDiffPY
TH1F *	fHVtxDiffPZ
TH1F *	fHVtxDiffE
TH1F *	fHVtxPullPX
TH1F *	fHVtxPullPY
TH1F *	fHVtxPullPZ
TH1F *	fHVtxPullE
TH1F *	fHVtxDiffThe
TH1F *	fHVtxDiffPhi
TH1F *	fHVtxPocaX
TH1F *	fHVtxPocaY
TH1F *	fHVtxPocaZ
TH1F *	fHVtxPullPocaX
TH1F *	fHVtxPullPocaY
TH1F *	fHVtxPullPocaZ
TH1F *	fHVtxPocas
TH1F *	fHVtxPocaEmpty
TH2F *	fHVtxPocaXY
TH2F *	fHVtxPocaRZ
TH1F *	fHVtxFastX
TH1F *	fHVtxFastY
TH1F *	fHVtxFastZ
TH1F *	fHVtxErrFastX
TH1F *	fHVtxErrFastY
TH1F *	fHVtxErrFastZ
TH1F *	fHVtxPullFastX
TH1F *	fHVtxPullFastY
TH1F *	fHVtxPullFastZ
TH1F *	fHVtxChi2Fast
TH1F *	fHVtxChiProbFast
TH2F *	fHVtxFastXY
TH2F *	fHVtxFastRZ
TH1F *	fHVtxFitX
TH1F *	fHVtxFitY
TH1F *	fHVtxFitZ
TH1F *	fHVtxErrFitX
TH1F *	fHVtxErrFitY
TH1F *	fHVtxErrFitZ
TH1F *	fHVtxPullFitX
TH1F *	fHVtxPullFitY
TH1F *	fHVtxPullFitZ
TH1F *	fHVtxDiffFitPX
TH1F *	fHVtxDiffFitPY
TH1F *	fHVtxDiffFitPZ
TH1F *	fHVtxDiffFitE
TH1F *	fHVtxDiffFitThe
TH1F *	fHVtxDiffFitPhi
TH1F *	fHVtxPullFitPX
TH1F *	fHVtxPullFitPY
TH1F *	fHVtxPullFitPZ
TH1F *	fHVtxPullFitE
TH1F *	fHVtxChi2Fit
TH1F *	fHVtxChiProbFit
TH2F *	fHVtxFitXY
TH2F *	fHVtxFitRZ
TH1F *	fHVtxKinX
TH1F *	fHVtxKinY
TH1F *	fHVtxKinZ
TH1F *	fHVtxErrKinX
TH1F *	fHVtxErrKinY
TH1F *	fHVtxErrKinZ
TH1F *	fHVtxPullKinX
TH1F *	fHVtxPullKinY
TH1F *	fHVtxPullKinZ
TH1F *	fHVtxDiffKinPX
TH1F *	fHVtxDiffKinPY
TH1F *	fHVtxDiffKinPZ
TH1F *	fHVtxDiffKinE
TH1F *	fHVtxDiffKinThe
TH1F *	fHVtxDiffKinPhi
TH1F *	fHVtxPullKinPX
TH1F *	fHVtxPullKinPY
TH1F *	fHVtxPullKinPZ
TH1F *	fHVtxPullKinE
TH1F *	fHVtxChi2Kin
TH1F *	fHVtxChiProbKin
TH2F *	fHVtxKinXY
TH2F *	fHVtxKinRZ
TH1F *	fHPrgPull0
TH1F *	fHPrgPull1
TH1F *	fHPrgPull2
TH1F *	fHPrgPull3
TH1F *	fHPrgPull4
TH1D *	fHCpu
TStopwatch	fSwAll
TStopwatch	fSwPMT
TStopwatch	fSwPoca
TStopwatch	fSwPrgfast
TStopwatch	fSwPrgfull
TStopwatch	fSwKin
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

Private Member Functions
virtual void	SetParContainers ()
	ClassDef (PndPmtTask, 1)

Detailed Description

Definition at line 17 of file PndPmtTask.h.

Constructor & Destructor Documentation

PndPmtTask::PndPmtTask

(

)

Default constructor

Definition at line 31 of file PndPmtTask.cxx.

References fHistoList, and fSwAll.

                       : PndPmtPoormantracks()   // : FairTask("Poorman Track Task")
{
 fSwAll.Start();
 fHistoList = new TList();
}

PndPmtTask::~PndPmtTask

(

)

Destructor

Definition at line 39 of file PndPmtTask.cxx.

{ }

Member Function Documentation

PndPmtTask::ClassDef ( PndPmtTask  , 1    )

private

void PndPmtPoormantracks::EnergyCorrelations ( )

inherited

Referenced by PndPmtPoormantracks::PoorManTracks().

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

protectedinherited

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 PndPmtTask::Exec

(

Option_t *

opt

)

Virtual method Exec

Definition at line 327 of file PndPmtTask.cxx.

References Bool_t, RhoCandidate::Charge(), RhoCandidate::Combine(), RhoCandidate::Cov7(), RhoCandidate::Daughter(), PndPmtPoormantracks::fCands, fHPrgPull0, fHPrgPull1, fHPrgPull2, fHPrgPull3, fHPrgPull4, fHVtxChi2Fast, fHVtxChi2Fit, fHVtxChi2Kin, fHVtxChiProbFast, fHVtxChiProbFit, fHVtxChiProbKin, fHVtxDiffE, fHVtxDiffFitE, fHVtxDiffFitPhi, fHVtxDiffFitPX, fHVtxDiffFitPY, fHVtxDiffFitPZ, fHVtxDiffFitThe, fHVtxDiffKinE, fHVtxDiffKinPhi, fHVtxDiffKinPX, fHVtxDiffKinPY, fHVtxDiffKinPZ, fHVtxDiffKinThe, fHVtxDiffPhi, fHVtxDiffPX, fHVtxDiffPY, fHVtxDiffPZ, fHVtxDiffThe, fHVtxDiffX, fHVtxDiffY, fHVtxDiffZ, fHVtxErrFastX, fHVtxErrFastY, fHVtxErrFastZ, fHVtxErrFitX, fHVtxErrFitY, fHVtxErrFitZ, fHVtxErrKinX, fHVtxErrKinY, fHVtxErrKinZ, fHVtxFastRZ, fHVtxFastX, fHVtxFastXY, fHVtxFastY, fHVtxFastZ, fHVtxFitRZ, fHVtxFitX, fHVtxFitXY, fHVtxFitY, fHVtxFitZ, fHVtxKinRZ, fHVtxKinX, fHVtxKinXY, fHVtxKinY, fHVtxKinZ, fHVtxPocaRZ, fHVtxPocas, fHVtxPocaX, fHVtxPocaXY, fHVtxPocaY, fHVtxPocaZ, fHVtxPullE, fHVtxPullFastX, fHVtxPullFastY, fHVtxPullFastZ, fHVtxPullFitE, fHVtxPullFitPX, fHVtxPullFitPY, fHVtxPullFitPZ, fHVtxPullFitX, fHVtxPullFitY, fHVtxPullFitZ, fHVtxPullKinE, fHVtxPullKinPX, fHVtxPullKinPY, fHVtxPullKinPZ, fHVtxPullKinX, fHVtxPullKinY, fHVtxPullKinZ, fHVtxPullPocaX, fHVtxPullPocaY, fHVtxPullPocaZ, fHVtxPullPX, fHVtxPullPY, fHVtxPullPZ, RhoFitterBase::Fit(), RhoKalmanVtxFitter::FitVertexFast(), fSwKin, fSwPMT, fSwPoca, fSwPrgfast, fSwPrgfull, PndPmtPoormantracks::fVerbose, PndPmtPoormantracks::fVertex, PndPmtPoormantracks::fWidx, PndPmtPoormantracks::fWidy, PndPmtPoormantracks::fWidz, RhoFitterBase::GetChi2(), RhoCandidate::GetFit(), RhoCandidate::GetMcTruth(), RhoVtxPoca::GetPocaVtx(), RhoCandidate::GetPosition(), RhoFitterBase::GetProb(), mom, mommc, RhoCandidate::NDaughters(), RhoCandidate::P4(), RhoCandidate::P4Cov(), PndPmtPoormantracks::PoorManTracks(), RhoCandidate::Pos(), PndPmtPoormantracks::RollVertexBox(), RhoKalmanVtxFitter::SetExpansionPoint(), RhoKalmanVtxFitter::SetNIterations(), RhoKinVtxFitter::SetNIterationsExact(), RhoKalmanVtxFitter::SetSilent(), and sqrt().

Referenced by runpmt().

                               {


  RhoCandidate *aCand=0,*bCand=0,*cCand=0,*dCand=0;
  RhoCandidate *tmpcand=0,*mccand=0;
  TVector3 vertexMC,vertexPoc,vertexFast,vertexFit,vertexKin,vertexDiff,vertexDiffFast,vertexDiffFit,vertexDiffKin;
  TLorentzVector mom,mommc,momdiff;
  const TVector3 nullpunkt(0.,0.,0.);
  double rnd1=0.,rnd2=0.,rnd3=0.;
  RhoError amomcov(4);
  TMatrixD aposcov(3,3);
  TMatrixD helixCov(5,5);
  TMatrixD jacobian(5,7);
  Float_t mcparams[5];
  Float_t helixparams[5];
  TMatrixD covV(3,3);
  TMatrixD covVfit(3,3);

  //for(int i=0;i<nevt;i++)
  //{
    //    printf("Event %i of %i. (seed=%i)\n",i,nevt,ranseed);
    RhoCandidate *combiCand=0,*combiCand1=0,*combiCand2=0;
    vertexMC = RollVertexBox(fWidx,fWidy,fWidz);
    fVertex=vertexMC;

    fSwPMT.Start(kFALSE);
    PoorManTracks();
    fSwPMT.Stop();
    if (fCands->GetEntriesFast()<2) {
      cout<<"Warning: too few particles created."<<endl;
      return;
    }

    aCand=(RhoCandidate*)fCands->At(0);
    bCand=(RhoCandidate*)fCands->At(1);
    if(fCands->GetEntriesFast()==3)
    {
      cCand=(RhoCandidate*)fCands->At(2);
      combiCand = aCand->Combine(bCand,cCand);
      combiCand1 = aCand->Combine(bCand,cCand);
      combiCand2 = aCand->Combine(bCand,cCand);
    }else if(fCands->GetEntriesFast()==4)
    {
      cCand=(RhoCandidate*)fCands->At(2);
      dCand=(RhoCandidate*)fCands->At(3);
      combiCand = aCand->Combine(bCand,cCand,dCand);
      combiCand1 = aCand->Combine(bCand,cCand,dCand);
      combiCand2 = aCand->Combine(bCand,cCand,dCand);
    }else{
      combiCand = aCand->Combine(bCand);
      combiCand1 = aCand->Combine(bCand);
      combiCand2 = aCand->Combine(bCand);
    }

    if(fVerbose>1) cout<<"aCand, bCand, comboCand, aPos, bPos, combipos"<<endl;
    if(fVerbose>1) cout<<*aCand<<endl;
    if(fVerbose>1) cout<<*bCand<<endl;
    if(fVerbose>1) cout<<*combiCand<<endl;
    if(fVerbose>1) aCand->GetPosition().Print();
    if(fVerbose>1) bCand->GetPosition().Print();
    if(fVerbose>1) combiCand->GetPosition().Print();


    // CLEAN
    for(int k=0;k<combiCand->NDaughters();k++)
    {
      tmpcand = combiCand->Daughter(k);
      mccand = tmpcand->GetMcTruth();

      vertexDiff=mccand->Pos();vertexDiff-=tmpcand->Pos();
      fHVtxDiffX->Fill(vertexDiff.X());
      fHVtxDiffY->Fill(vertexDiff.Y());
      fHVtxDiffZ->Fill(vertexDiff.Z());

      mommc=mccand->P4();
      mom=tmpcand->P4();
      momdiff=mommc-mom;
      amomcov=tmpcand->P4Cov();
      fHVtxDiffPX->Fill(momdiff.X());
      fHVtxDiffPY->Fill(momdiff.Y());
      fHVtxDiffPZ->Fill(momdiff.Z());
      fHVtxDiffE->Fill(momdiff.E());

      fHVtxPullPX->Fill(momdiff.X()/sqrt(amomcov[0][0]));
      fHVtxPullPY->Fill(momdiff.Y()/sqrt(amomcov[1][1]));
      fHVtxPullPZ->Fill(momdiff.Z()/sqrt(amomcov[2][2]));
      fHVtxPullE->Fill(momdiff.E()/sqrt(amomcov[3][3]));

      fHVtxDiffThe->Fill(mommc.Theta()-mom.Theta());
      fHVtxDiffPhi->Fill(mommc.Phi()-mom.Phi());

    }

    // ********** POCA
    RhoVtxPoca vPoca;
    vertexPoc=vertexMC; //seed
    if(fVerbose>0) cout<<"poca "<<flush;
    fSwPoca.Start(kFALSE);
    double dist = vPoca.GetPocaVtx(vertexPoc,combiCand);
    fSwPoca.Stop();
    if(fVerbose>0) cout<<" - Dist = "<<dist<<"  \t";if(fVerbose>1)cout<<endl;
    vertexDiff=vertexMC - vertexPoc;
    if(fVerbose>1) vertexMC.Print();
    if(fVerbose>1) vertexPoc.Print();
    if(fVerbose>1) vertexDiff.Print();

    fHVtxPocaX->Fill(vertexDiff.X());
    fHVtxPocaY->Fill(vertexDiff.Y());
    fHVtxPocaZ->Fill(vertexDiff.Z());
    fHVtxPocaXY->Fill(vertexDiff.X(),vertexDiff.Y());
    fHVtxPocaRZ->Fill(vertexDiff.Z(),vertexMC.Perp()-vertexPoc.Perp());
    fHVtxPocas->Fill(dist);
    fHVtxPullPocaX->Fill(vertexDiff.X()/dist);
    fHVtxPullPocaY->Fill(vertexDiff.Y()/dist);
    fHVtxPullPocaZ->Fill(vertexDiff.Z()/dist);

    // ********** FAST FIT
    RhoKalmanVtxFitter vFastter(combiCand);
    vFastter.SetSilent();
    //vFastter.SetDebug();
    //PndAnalysisCalcTools::SetVerbose(3);
    if(fVerbose>0) cout<<"prg "<<flush;
    vertexFast=nullpunkt; //seed
    double chiq = -5454;
    vFastter.SetExpansionPoint(vertexFast);
    fSwPrgfast.Start(kFALSE);
    chiq = vFastter.FitVertexFast(vertexFast,covV,false,2); // 2 iterations, don't skip cov
    fSwPrgfast.Stop();
    if (chiq>=0){
      if(fVerbose>0) cout<<" - chiq = "<<chiq<<"  \t";if(fVerbose>1)cout<<endl;

      fHVtxChi2Fast->Fill(vFastter.GetChi2());
      fHVtxChiProbFast->Fill(vFastter.GetProb());

      vertexDiffFast=vertexMC - vertexFast;
      if(fVerbose>1) vertexMC.Print();
      if(fVerbose>1) vertexFast.Print();
      if(fVerbose>0) vertexDiffFast.Print();

      fHVtxFastX->Fill(vertexDiffFast.X());
      fHVtxFastY->Fill(vertexDiffFast.Y());
      fHVtxFastZ->Fill(vertexDiffFast.Z());
      fHVtxFastXY->Fill(vertexDiffFast.X(),vertexDiffFast.Y());
      fHVtxFastRZ->Fill(vertexDiffFast.Z(),vertexMC.Perp()-vertexFast.Perp());

      fHVtxErrFastX->Fill(sqrt(covV[0][0]));
      fHVtxErrFastY->Fill(sqrt(covV[1][1]));
      fHVtxErrFastZ->Fill(sqrt(covV[2][2]));
      fHVtxPullFastX->Fill(vertexDiffFast.X()/sqrt(covV[0][0]));
      fHVtxPullFastY->Fill(vertexDiffFast.Y()/sqrt(covV[1][1]));
      fHVtxPullFastZ->Fill(vertexDiffFast.Z()/sqrt(covV[2][2]));
    }
    //vFastter.SetSilent();
    //PndAnalysisCalcTools::SetVerbose(0);

    // ********** FULL FIT
    RhoKalmanVtxFitter vFitter(combiCand1);
    vFitter.SetSilent();
    //PndAnalysisCalcTools::SetVerbose(5);
    //if(laut)vFitter.SetDebug(true);
    vertexFit=nullpunkt; //seed
    vFitter.SetNIterations(2);
    if(fVerbose>0)  cout<<"prg  "<<flush;
    //if(laut>0)  PndAnalysisCalcTools::SetVerbose(3);
    fSwPrgfull.Start(kFALSE);
    bool stat = vFitter.Fit();
    fSwPrgfull.Stop();
    //if(laut>0)  PndAnalysisCalcTools::SetVerbose(0);
    //if(fVerbose>0)  cout<<" - chiq = "<<chiqfit<<"  \t"<<flush;if(fVerbose>1)cout<<endl;
    if (stat)
    {
      //PndAnalysisCalcTools::SetVerbose(0);
      //      fHVtxChi2Fit->Fill(chiqfit/fNDF);
      //      fHVtxChiProbFit->Fill(TMath::Prob(chiqfit,fNDF));
      fHVtxChi2Fit->Fill(vFitter.GetChi2());
      fHVtxChiProbFit->Fill(vFitter.GetProb());

      vertexFit=combiCand1->GetFit()->DecayVtx();
      covVfit=combiCand1->GetFit()->DecayVtx().CovMatrix();
      vertexDiffFit=vertexMC - vertexFit;
      if(fVerbose>1) vertexMC.Print();
      if(fVerbose>1) vertexFit.Print();
      if(fVerbose>1) vertexDiffFit.Print();

      fHVtxFitX->Fill(vertexDiffFit.X());
      fHVtxFitY->Fill(vertexDiffFit.Y());
      fHVtxFitZ->Fill(vertexDiffFit.Z());
      fHVtxFitXY->Fill(vertexDiffFit.X(),vertexDiffFit.Y());
      fHVtxFitRZ->Fill(vertexDiffFit.Z(),vertexMC.Perp()-vertexFit.Perp());

      fHVtxErrFitX->Fill(sqrt(covVfit[0][0]));
      fHVtxErrFitY->Fill(sqrt(covVfit[1][1]));
      fHVtxErrFitZ->Fill(sqrt(covVfit[2][2]));
      fHVtxPullFitX->Fill(vertexDiffFit.X()/sqrt(covVfit[0][0]));
      fHVtxPullFitY->Fill(vertexDiffFit.Y()/sqrt(covVfit[1][1]));
      fHVtxPullFitZ->Fill(vertexDiffFit.Z()/sqrt(covVfit[2][2]));

      for(int k=0;k<combiCand1->NDaughters();k++)
      {
        tmpcand = combiCand1->Daughter(k);
        mccand = tmpcand->GetMcTruth();
        tmpcand = tmpcand->GetFit();
        mommc=mccand->P4();
        mom=tmpcand->P4();
        momdiff=mommc-mom;
        amomcov=tmpcand->P4Cov();
        fHVtxDiffFitPX->Fill(momdiff.X());
        fHVtxDiffFitPY->Fill(momdiff.Y());
        fHVtxDiffFitPZ->Fill(momdiff.Z());
        fHVtxDiffFitE->Fill(momdiff.E());

        fHVtxPullFitPX->Fill(momdiff.X()/sqrt(amomcov[0][0]));
        fHVtxPullFitPY->Fill(momdiff.Y()/sqrt(amomcov[1][1]));
        fHVtxPullFitPZ->Fill(momdiff.Z()/sqrt(amomcov[2][2]));
        fHVtxPullFitE->Fill(momdiff.E()/sqrt(amomcov[3][3]));

        fHVtxDiffFitThe->Fill(mommc.Theta()-mom.Theta());
        fHVtxDiffFitPhi->Fill(mommc.Phi()-mom.Phi());
            // Now check the PRG description....
        Bool_t prgcheckmc = PndAnalysisCalcTools::P7toPRG(mccand->Pos(),mccand->P4(),mccand->Charge(),mccand->Cov7(),nullpunkt, mcparams, helixCov, jacobian, kFALSE);
        Bool_t prgchecktc = PndAnalysisCalcTools::P7toPRG(tmpcand->Pos(),tmpcand->P4(),tmpcand->Charge(),tmpcand->Cov7(),nullpunkt, helixparams, helixCov, jacobian, kFALSE);
        fHPrgPull0->Fill((mcparams[0]-helixparams[0])/sqrt(helixCov[0][0]));
        fHPrgPull1->Fill((mcparams[1]-helixparams[1])/sqrt(helixCov[1][1]));
        fHPrgPull2->Fill((mcparams[2]-helixparams[2])/sqrt(helixCov[2][2]));
        fHPrgPull3->Fill((mcparams[3]-helixparams[3])/sqrt(helixCov[3][3]));
        fHPrgPull4->Fill((mcparams[4]-helixparams[4])/sqrt(helixCov[4][4]));
      }
    }

    // ********** Kin FIT
    RhoKinVtxFitter kFitter(combiCand2);
    kFitter.SetNIterationsExact(2);
    vertexFit=nullpunkt; //seed
    if(fVerbose>0)  cout<<"kin  "<<flush;
    fSwKin.Start(kFALSE);
    kFitter.Fit();
    fSwKin.Stop();
    bool ksta = kFitter.GetChi2();
    //if(laut>0)  PndAnalysisCalcTools::SetVerbose(0);
    //if(fVerbose>0)  cout<<" - chiq = "<<chiqfit<<"  \t"<<flush;if(fVerbose>1)cout<<endl;
    if (ksta)
    {
      //PndAnalysisCalcTools::SetVerbose(0);
      vertexKin=combiCand2->GetFit()->DecayVtx();
      aposcov=combiCand2->GetFit()->DecayVtx().CovMatrix();
      fHVtxChi2Kin->Fill(kFitter.GetChi2());
      fHVtxChiProbKin->Fill(kFitter.GetProb());

      vertexDiffKin=vertexMC - vertexKin;
      if(fVerbose>1) vertexMC.Print();
      if(fVerbose>1) vertexKin.Print();
      if(fVerbose>0) vertexDiffKin.Print();

      fHVtxKinX->Fill(vertexDiffKin.X());
      fHVtxKinY->Fill(vertexDiffKin.Y());
      fHVtxKinZ->Fill(vertexDiffKin.Z());
      fHVtxKinXY->Fill(vertexDiffKin.X(),vertexDiffKin.Y());
      fHVtxKinRZ->Fill(vertexDiffKin.Z(),vertexMC.Perp()-vertexKin.Perp());

      fHVtxErrKinX->Fill(sqrt(aposcov[0][0]));
      fHVtxErrKinY->Fill(sqrt(aposcov[1][1]));
      fHVtxErrKinZ->Fill(sqrt(aposcov[2][2]));
      fHVtxPullKinX->Fill(vertexDiffKin.X()/sqrt(aposcov[0][0]));
      fHVtxPullKinY->Fill(vertexDiffKin.Y()/sqrt(aposcov[1][1]));
      fHVtxPullKinZ->Fill(vertexDiffKin.Z()/sqrt(aposcov[2][2]));

      for(int k=0;k<combiCand2->NDaughters();k++)
      {
        tmpcand = combiCand2->Daughter(k);
        mccand = tmpcand->GetMcTruth();
        tmpcand = tmpcand->GetFit();
        mommc=mccand->P4();
        mom=tmpcand->P4();
        momdiff=mommc-mom;
        amomcov=tmpcand->P4Cov();
        fHVtxDiffKinPX->Fill(momdiff.X());
        fHVtxDiffKinPY->Fill(momdiff.Y());
        fHVtxDiffKinPZ->Fill(momdiff.Z());
        fHVtxDiffKinE->Fill(momdiff.E());

        fHVtxPullKinPX->Fill(momdiff.X()/sqrt(amomcov[0][0]));
        fHVtxPullKinPY->Fill(momdiff.Y()/sqrt(amomcov[1][1]));
        fHVtxPullKinPZ->Fill(momdiff.Z()/sqrt(amomcov[2][2]));
        fHVtxPullKinE->Fill(momdiff.E()/sqrt(amomcov[3][3]));

        fHVtxDiffKinThe->Fill(mommc.Theta()-mom.Theta());
        fHVtxDiffKinPhi->Fill(mommc.Phi()-mom.Phi());
      }
    }
    // } // event loop

    //cout<<" ---------- calculations done, writing histograms."<<endl;
    //manipulate & fit histos now


}

void PndPmtTask::Finish

(

)

Definition at line 625 of file PndPmtTask.cxx.

References fHCpu, fHistoList, fSwAll, fSwKin, fSwPMT, fSwPoca, fSwPrgfast, fSwPrgfull, and TString.

Referenced by runpmt().

{
  fSwAll.Stop();
  //Read Stopwatches
  //fHCpu->Fill("PMT",fSwPMT.RealTime()/fSwAll.RealTime());
  //fHCpu->Fill("POCA",fSwPoca.RealTime()/fSwAll.RealTime());
  //fHCpu->Fill("PRG Fast",fSwPrgfast.RealTime()/fSwAll.RealTime());
  //fHCpu->Fill("PRG Full",fSwPrgfull.RealTime()/fSwAll.RealTime());
  //fHCpu->Fill("KinVtx",fSwKin.RealTime()/fSwAll.RealTime());
  fHCpu->Fill("PMT",1000*fSwPMT.RealTime()/fSwPMT.Counter());
  fHCpu->Fill("POCA",1000*fSwPoca.RealTime()/fSwPMT.Counter());
  fHCpu->Fill("PRG Fast",1000*fSwPrgfast.RealTime()/fSwPMT.Counter());
  fHCpu->Fill("PRG Full",1000*fSwPrgfull.RealTime()/fSwPMT.Counter());
  fHCpu->Fill("KinVtx",1000*fSwKin.RealTime()/fSwPMT.Counter());


  //Write and Plot all histos from List
  TString hfile="Data/HistoVertexing";
  //  if(seed>=0){hfile+="-";hfile+=ranseed;}
  hfile+=".root";
  TFile* histofile = new TFile(hfile.Data(),"RECREATE");
  histofile->cd();
  TListIter iter(fHistoList);
  TH1* tmph=NULL;
  while( tmph=(TH1*)iter() ) tmph->Write();
  histofile->Close();
}

void PndPmtTask::Init

(

)

Virtual method Init

Definition at line 46 of file PndPmtTask.cxx.

References fHCpu, fHistoList, fHPrgPull0, fHPrgPull1, fHPrgPull2, fHPrgPull3, fHPrgPull4, fHVtxChi2Fast, fHVtxChi2Fit, fHVtxChi2Kin, fHVtxChiProbFast, fHVtxChiProbFit, fHVtxChiProbKin, fHVtxDiffE, fHVtxDiffFitE, fHVtxDiffFitPhi, fHVtxDiffFitPX, fHVtxDiffFitPY, fHVtxDiffFitPZ, fHVtxDiffFitThe, fHVtxDiffKinE, fHVtxDiffKinPhi, fHVtxDiffKinPX, fHVtxDiffKinPY, fHVtxDiffKinPZ, fHVtxDiffKinThe, fHVtxDiffPhi, fHVtxDiffPX, fHVtxDiffPY, fHVtxDiffPZ, fHVtxDiffThe, fHVtxDiffX, fHVtxDiffY, fHVtxDiffZ, fHVtxErrFastX, fHVtxErrFastY, fHVtxErrFastZ, fHVtxErrFitX, fHVtxErrFitY, fHVtxErrFitZ, fHVtxErrKinX, fHVtxErrKinY, fHVtxErrKinZ, fHVtxFastRZ, fHVtxFastX, fHVtxFastXY, fHVtxFastY, fHVtxFastZ, fHVtxFitRZ, fHVtxFitX, fHVtxFitXY, fHVtxFitY, fHVtxFitZ, fHVtxKinRZ, fHVtxKinX, fHVtxKinXY, fHVtxKinY, fHVtxKinZ, fHVtxPocaEmpty, fHVtxPocaRZ, fHVtxPocas, fHVtxPocaX, fHVtxPocaXY, fHVtxPocaY, fHVtxPocaZ, fHVtxPullE, fHVtxPullFastX, fHVtxPullFastY, fHVtxPullFastZ, fHVtxPullFitE, fHVtxPullFitPX, fHVtxPullFitPY, fHVtxPullFitPZ, fHVtxPullFitX, fHVtxPullFitY, fHVtxPullFitZ, fHVtxPullKinE, fHVtxPullKinPX, fHVtxPullKinPY, fHVtxPullKinPZ, fHVtxPullKinX, fHVtxPullKinY, fHVtxPullKinZ, fHVtxPullPocaX, fHVtxPullPocaY, fHVtxPullPocaZ, fHVtxPullPX, fHVtxPullPY, fHVtxPullPZ, RhoCalculationTools::ForceConstantBz(), PndPmtPoormantracks::fSeed, fSwKin, fSwPMT, fSwPoca, fSwPrgfast, fSwPrgfull, PndPmtPoormantracks::fWidx, PndPmtPoormantracks::fWidy, PndPmtPoormantracks::fWidz, max(), min(), and SetVerbose().

Referenced by runpmt().

                      {

  gRandom->SetSeed(fSeed);
  RhoCalculationTools::ForceConstantBz(20.); // [kGs]=[0.1T]
  PndAnalysisCalcTools::SetVerbose(0);

  //PArameters
  fWidx=1.5,
  fWidy=1.5,
  fWidz=1.5; //[cm] maximum vertex position off origin

  int bins=100;
  double min=0,max=4.5,mdf=2.0; // GeV/c
  double min2=0.13,max2=0.15;
  double min3=1.,max3=1.5;
  double max4=0.3;//cm
  double minchi = -0.05;
  double maxchi = 2.5;
  double vrange = 0.05; // cm
  double pullrange =4.;
  double momrange=0.1;//GeV/c
  double angdiffran=0.05; //rad


  // NO FIT
  fHVtxDiffX=new TH1F("hvtxDiffx","Vertex Smearing;(x_{mc}-x)/cm",bins,-vrange,vrange);
  fHVtxDiffY=new TH1F("hvtxDiffy","Vertex Smearing;(y_{mc}-y)/cm",bins,-vrange,vrange);
  fHVtxDiffZ=new TH1F("hvtxDiffz","Vertex Smearing;(z_{mc}-z)/cm",bins,-vrange,vrange);

  fHVtxDiffPX=new TH1F("hvtxdiffpx","Momentum reco;(p_{x,mc}-p_{x} / GeV/c)",bins,-momrange,momrange);
  fHVtxDiffPY=new TH1F("hvtxdiffpy","Momentum reco;(p_{y,mc}-p_{y} / GeV/c)",bins,-momrange,momrange);
  fHVtxDiffPZ=new TH1F("hvtxdiffpz","Momentum reco;(p_{z,mc}-p_{z}) / GeV/c",bins,-momrange,momrange);
  fHVtxDiffE=new TH1F("hvtxdiffe","Energy reco;(E_{mc}-E) / GeV",bins,-momrange,momrange);

  fHVtxPullPX=new TH1F("hvtxpullpx","Momentum pull distribution reco;(p_{x,mc}-p_{x})/#sigma_{p_{x}}",bins,-pullrange,pullrange);
  fHVtxPullPY=new TH1F("hvtxpullpy","Momentum pull distribution reco;(p_{y,mc}-p_{y})/#sigma_{p_{y}}",bins,-pullrange,pullrange);
  fHVtxPullPZ=new TH1F("hvtxpullpz","Momentum pull distribution reco;(p_{z,mc}-p_{z})/#sigma_{p_{z}}",bins,-pullrange,pullrange);
  fHVtxPullE=new TH1F("hvtxpulle","Energy pull distribution reco;(E_{mc}-E)/#sigma_{E}",bins,-pullrange,pullrange);

  fHVtxDiffThe=new TH1F("hvtxdiffthe","Momentum reco;(#Theta_{mc}-#Theta}) / GeV/c",bins,-angdiffran,angdiffran);
  fHVtxDiffPhi=new TH1F("hvtxdiffphi","Momentum reco;(#Phi_{mc}-#Phi}) / GeV/c",bins,-angdiffran,angdiffran);

  // POCA

  fHVtxPocaX=new TH1F("hvtxPocax","Vertex POCA;(x_{mc}-x)/cm",bins,-vrange,vrange);
  fHVtxPocaY=new TH1F("hvtxPocay","Vertex POCA;(y_{mc}-y)/cm",bins,-vrange,vrange);
  fHVtxPocaZ=new TH1F("hvtxPocaz","Vertex POCA;(z_{mc}-z)/cm",bins,-vrange,vrange);

  fHVtxPocaXY=new TH2F("hvtxPocaxy","Vertex POCA X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
  fHVtxPocaRZ=new TH2F("hvtxPocarz","Vertex POCA R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);

  fHVtxPullPocaX=new TH1F("hvtxpullpocax","Vertex pull distribution POCA;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
  fHVtxPullPocaY=new TH1F("hvtxpullpocay","Vertex pull distribution POCA;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
  fHVtxPullPocaZ=new TH1F("hvtxpullpocaz","Vertex pull distribution POCA;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);

  fHVtxPocas=new TH1F("hvtxpocas","Vertex POCA D value;D/cm",bins,0.,vrange);
  fHVtxPocaEmpty=new TH1F("hvtxcpocaempty","",0,0.,1.);

  // FAST FIT

  fHVtxFastX=new TH1F("hvtxFastx","Vertex fast fit;(x_{mc}-x)/cm",bins,-vrange,vrange);
  fHVtxFastY=new TH1F("hvtxTasty","Vertex fast fit;(y_{mc}-y)/cm",bins,-vrange,vrange);
  fHVtxFastZ=new TH1F("hvtxFastz","Vertex fast fit;(z_{mc}-z)/cm",bins,-vrange,vrange);

  fHVtxFastXY=new TH2F("hvtxfastxy","Vertex fast fit X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
  fHVtxFastRZ=new TH2F("hvtxfastrz","Vertex fast fit R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);

  fHVtxErrFastX=new TH1F("hvtxerrfastx","Vertex error distribution fast;#sigma_{x}",bins,0.,vrange);
  fHVtxErrFastY=new TH1F("hvtxerrfasty","Vertex error distribution fast;#sigma_{y}",bins,0.,vrange);
  fHVtxErrFastZ=new TH1F("hvtxerrfastz","Vertex error distribution fast;#sigma_{z}",bins,0.,vrange);
  fHVtxPullFastX=new TH1F("hvtxpullfastx","Vertex pull distribution fast;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
  fHVtxPullFastY=new TH1F("hvtxpullfasty","Vertex pull distribution fast;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
  fHVtxPullFastZ=new TH1F("hvtxpullfastz","Vertex pull distribution fast;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);

  fHVtxChi2Fast=new TH1F("hvtxchi2fast","Vertex Fast #Chi^{2};#Chi^{2}",bins,minchi,maxchi);
  fHVtxChiProbFast=new TH1F("hvtxchipropfast","P(#Chi^{2}) Fast;p",bins,0.,1.);

  // FULL FIT
  fHVtxFitX=new TH1F("hvtxFitx","Vertex full fit;(x_{mc}-x)/cm",bins,-vrange,vrange);
  fHVtxFitY=new TH1F("hvtxFity","Vertex full fit;(y_{mc}-y)/cm",bins,-vrange,vrange);
  fHVtxFitZ=new TH1F("hvtxFitz","Vertex full fit;(z_{mc}-z)/cm",bins,-vrange,vrange);

  fHVtxFitXY=new TH2F("hvtxfitxy","Vertex full fit X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
  fHVtxFitRZ=new TH2F("hvtxfitrz","Vertex full fit R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);

  fHVtxErrFitX=new TH1F("hvtxerrfitx","Vertex error distribution full fit;#sigma_{x}",bins,0.,vrange);
  fHVtxErrFitY=new TH1F("hvtxerrfity","Vertex error distribution full fit;#sigma_{y}",bins,0.,vrange);
  fHVtxErrFitZ=new TH1F("hvtxerrfitz","Vertex error distribution full fit;#sigma_{z}",bins,0.,vrange);
  fHVtxPullFitX=new TH1F("hvtxpullfitx","Vertex pull distribution full fit;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
  fHVtxPullFitY=new TH1F("hvtxpullfity","Vertex pull distribution full fit;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
  fHVtxPullFitZ=new TH1F("hvtxpullfitz","Vertex pull distribution full fit;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);
  fHVtxDiffFitPX=new TH1F("hvtxdifffitpx","Momentum full fit;(p_{x,mc}-p_{x} / GeV/c)",bins,-momrange,momrange);
  fHVtxDiffFitPY=new TH1F("hvtxdifffitpy","Momentum full fit;(p_{y,mc}-p_{y} / GeV/c)",bins,-momrange,momrange);
  fHVtxDiffFitPZ=new TH1F("hvtxdifffitpz","Momentum full fit;(p_{z,mc}-p_{z}) / GeV/c",bins,-momrange,momrange);
  fHVtxDiffFitE=new TH1F("hvtxdifffite","Energy full fit;(E_{mc}-E) / GeV",bins,-momrange,momrange);
  fHVtxDiffFitThe=new TH1F("hvtxdifffitthe","Momentum full fit;(#Theta_{mc}-#Theta}) / GeV/c",bins,-angdiffran,angdiffran);
  fHVtxDiffFitPhi=new TH1F("hvtxdifffitphi","Momentum full fit;(#Phi_{mc}-#Phi}) / GeV/c",bins,-angdiffran,angdiffran);

  fHVtxPullFitPX=new TH1F("hvtxpullfitpx","Momentum pull distribution full fit;(p_{x,mc}-p_{x})/#sigma_{p_{x}}",bins,-pullrange,pullrange);
  fHVtxPullFitPY=new TH1F("hvtxpullfitpy","Momentum pull distribution full fit;(p_{y,mc}-p_{y})/#sigma_{p_{y}}",bins,-pullrange,pullrange);
  fHVtxPullFitPZ=new TH1F("hvtxpullfitpz","Momentum pull distribution full fit;(p_{z,mc}-p_{z})/#sigma_{p_{z}}",bins,-pullrange,pullrange);
  fHVtxPullFitE=new TH1F("hvtxpullfite","Energy pull distribution full fit;(E_{mc}-E)/#sigma_{E}",bins,-pullrange,pullrange);

  fHVtxChi2Fit=new TH1F("hvtxchi2fit","Vertex Fit #Chi^{2};#Chi^{2}",bins,minchi,maxchi);
  fHVtxChiProbFit=new TH1F("hvtxchipropfit","P(#Chi^{2}) Fit;p",bins,0.,1.);

  // KIN Fit

  fHVtxKinX=new TH1F("hvtxKinx","Vertex full kin;(x_{mc}-x)/cm",bins,-vrange,vrange);
  fHVtxKinY=new TH1F("hvtxKiny","Vertex full kin;(y_{mc}-y)/cm",bins,-vrange,vrange);
  fHVtxKinZ=new TH1F("hvtxKinz","Vertex full kin;(z_{mc}-z)/cm",bins,-vrange,vrange);

  fHVtxKinXY=new TH2F("hvtxkinxy","Vertex full kin X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
  fHVtxKinRZ=new TH2F("hvtxkinrz","Vertex full kin R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);

  fHVtxErrKinX=new TH1F("hvtxerrkinx","Vertex error distribution full kin;#sigma_{x}",bins,0.,vrange);
  fHVtxErrKinY=new TH1F("hvtxerrkiny","Vertex error distribution full kin;#sigma_{y}",bins,0.,vrange);
  fHVtxErrKinZ=new TH1F("hvtxerrkinz","Vertex error distribution full kin;#sigma_{z}",bins,0.,vrange);
  fHVtxPullKinX=new TH1F("hvtxpullkinx","Vertex pull distribution full kin;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
  fHVtxPullKinY=new TH1F("hvtxpullkiny","Vertex pull distribution full kin;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
  fHVtxPullKinZ=new TH1F("hvtxpullkinz","Vertex pull distribution full kin;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);
  fHVtxDiffKinPX=new TH1F("hvtxdiffkinpx","Momentum full kin;(p_{x,mc}-p_{x} / GeV/c)",bins,-momrange,momrange);
  fHVtxDiffKinPY=new TH1F("hvtxdiffkinpy","Momentum full kin;(p_{y,mc}-p_{y} / GeV/c)",bins,-momrange,momrange);
  fHVtxDiffKinPZ=new TH1F("hvtxdiffkinpz","Momentum full kin;(p_{z,mc}-p_{z}) / GeV/c",bins,-momrange,momrange);
  fHVtxDiffKinE=new TH1F("hvtxdiffkine","Energy full kin;(E_{mc}-E) / GeV",bins,-momrange,momrange);
  fHVtxDiffKinThe=new TH1F("hvtxdiffkinthe","Momentum full kin;(#Theta_{mc}-#Theta}) / GeV/c",bins,-angdiffran,angdiffran);
  fHVtxDiffKinPhi=new TH1F("hvtxdiffkinphi","Momentum full kin;(#Phi_{mc}-#Phi}) / GeV/c",bins,-angdiffran,angdiffran);

  fHVtxPullKinPX=new TH1F("hvtxpullkinpx","Momentum pull distribution full kin;(p_{x,mc}-p_{x})/#sigma_{p_{x}}",bins,-pullrange,pullrange);
  fHVtxPullKinPY=new TH1F("hvtxpullkinpy","Momentum pull distribution full kin;(p_{y,mc}-p_{y})/#sigma_{p_{y}}",bins,-pullrange,pullrange);
  fHVtxPullKinPZ=new TH1F("hvtxpullkinpz","Momentum pull distribution full kin;(p_{z,mc}-p_{z})/#sigma_{p_{z}}",bins,-pullrange,pullrange);
  fHVtxPullKinE=new TH1F("hvtxpullkine","Energy pull distribution full kin;(E_{mc}-E)/#sigma_{E}",bins,-pullrange,pullrange);

  fHVtxChi2Kin=new TH1F("hvtxchi2kin","Vertex Kin #Chi^{2};#Chi^{2}",bins,minchi,maxchi);
  fHVtxChiProbKin=new TH1F("hvtxchipropkin","P(#Chi^{2}) Kin;p",bins,0.,1.);

  // Helix Params
  fHPrgPull0 = new TH1F("hprgpull0","PRG Pull Distribution (0)",bins,-pullrange,pullrange);
  fHPrgPull1 = new TH1F("hprgpull1","PRG Pull Distribution (1)",bins,-pullrange,pullrange);
  fHPrgPull2 = new TH1F("hprgpull2","PRG Pull Distribution (2)",bins,-pullrange,pullrange);
  fHPrgPull3 = new TH1F("hprgpull3","PRG Pull Distribution (3)",bins,-pullrange,pullrange);
  fHPrgPull4 = new TH1F("hprgpull4","PRG Pull Distribution (4)",bins,-pullrange,pullrange);

  fHCpu = new TH1D("hcpu","CPU loads (ms)",5,0,5);
  fHCpu->GetXaxis()->SetBinLabel(1,"PMT");
  fHCpu->GetXaxis()->SetBinLabel(2,"POCA");
  fHCpu->GetXaxis()->SetBinLabel(3,"PRG Fast");
  fHCpu->GetXaxis()->SetBinLabel(4,"PRG Full");
  fHCpu->GetXaxis()->SetBinLabel(5,"KinVtx");

  fHistoList->Add(fHVtxDiffX);
  fHistoList->Add(fHVtxDiffY);
  fHistoList->Add(fHVtxDiffZ);

  fHistoList->Add(fHVtxPocaX);
  fHistoList->Add(fHVtxPocaY);
  fHistoList->Add(fHVtxPocaZ);

  fHistoList->Add(fHVtxFastX);
  fHistoList->Add(fHVtxFastY);
  fHistoList->Add(fHVtxFastZ);

  fHistoList->Add(fHVtxFitX);
  fHistoList->Add(fHVtxFitY);
  fHistoList->Add(fHVtxFitZ);

  fHistoList->Add(fHVtxKinX);
  fHistoList->Add(fHVtxKinY);
  fHistoList->Add(fHVtxKinZ);

  fHistoList->Add(fHVtxErrFastX);
  fHistoList->Add(fHVtxErrFastY);
  fHistoList->Add(fHVtxErrFastZ);

  fHistoList->Add(fHVtxErrFitX);
  fHistoList->Add(fHVtxErrFitY);
  fHistoList->Add(fHVtxErrFitZ);

  fHistoList->Add(fHVtxErrKinX);
  fHistoList->Add(fHVtxErrKinY);
  fHistoList->Add(fHVtxErrKinZ);

  fHistoList->Add(fHVtxPullPocaX);
  fHistoList->Add(fHVtxPullPocaY);
  fHistoList->Add(fHVtxPullPocaZ);

  fHistoList->Add(fHVtxPullFastX);
  fHistoList->Add(fHVtxPullFastY);
  fHistoList->Add(fHVtxPullFastZ);

  fHistoList->Add(fHVtxPullFitX);
  fHistoList->Add(fHVtxPullFitY);
  fHistoList->Add(fHVtxPullFitZ);

  fHistoList->Add(fHVtxPullKinX);
  fHistoList->Add(fHVtxPullKinY);
  fHistoList->Add(fHVtxPullKinZ);

  fHistoList->Add(fHVtxDiffPX);
  fHistoList->Add(fHVtxDiffPY);
  fHistoList->Add(fHVtxDiffPZ);

  fHistoList->Add(fHVtxPullPX);
  fHistoList->Add(fHVtxPullPY);
  fHistoList->Add(fHVtxPullPZ);

  fHistoList->Add(fHVtxDiffFitPX);
  fHistoList->Add(fHVtxDiffFitPY);
  fHistoList->Add(fHVtxDiffFitPZ);

  fHistoList->Add(fHVtxPullFitPX);
  fHistoList->Add(fHVtxPullFitPY);
  fHistoList->Add(fHVtxPullFitPZ);

  fHistoList->Add(fHVtxDiffKinPX);
  fHistoList->Add(fHVtxDiffKinPY);
  fHistoList->Add(fHVtxDiffKinPZ);

  fHistoList->Add(fHVtxPullKinPX);
  fHistoList->Add(fHVtxPullKinPY);
  fHistoList->Add(fHVtxPullKinPZ);

  fHistoList->Add(fHVtxDiffE);
  fHistoList->Add(fHVtxDiffFitE);
  fHistoList->Add(fHVtxDiffKinE);
  fHistoList->Add(fHVtxPullE);
  fHistoList->Add(fHVtxPullFitE);
  fHistoList->Add(fHVtxPullKinE);

  fHistoList->Add(fHVtxDiffThe);
  fHistoList->Add(fHVtxDiffPhi);

  fHistoList->Add(fHVtxDiffFitThe);
  fHistoList->Add(fHVtxDiffFitPhi);

  fHistoList->Add(fHVtxDiffKinThe);
  fHistoList->Add(fHVtxDiffKinPhi);

  fHistoList->Add(fHVtxPocas);
  fHistoList->Add(fHVtxChi2Fast);
  fHistoList->Add(fHVtxChi2Fit);
  fHistoList->Add(fHVtxChi2Kin);
  fHistoList->Add(fHVtxChiProbFast);
  fHistoList->Add(fHVtxChiProbFit);
  fHistoList->Add(fHVtxChiProbKin);

  fHistoList->Add(fHVtxPocaXY);
  fHistoList->Add(fHVtxFastXY);
  fHistoList->Add(fHVtxFitXY);
  fHistoList->Add(fHVtxKinXY);
  fHistoList->Add(fHVtxPocaRZ);
  fHistoList->Add(fHVtxFastRZ);
  fHistoList->Add(fHVtxFitRZ);
  fHistoList->Add(fHVtxKinRZ);

  fHistoList->Add(fHPrgPull0);
  fHistoList->Add(fHPrgPull1);
  fHistoList->Add(fHPrgPull2);
  fHistoList->Add(fHPrgPull3);
  fHistoList->Add(fHPrgPull4);

  fHistoList->Add(fHCpu);

  fSwPMT.Reset();
  fSwPoca.Reset();
  fSwPrgfast.Reset();
  fSwPrgfull.Reset();
  fSwKin.Reset();
}

void PndPmtPoormantracks::PoorManTracks ( )

protectedinherited

Definition at line 95 of file PndPmtPoormantracks.cxx.

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

Referenced by 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  )

protectedinherited

Definition at line 72 of file PndPmtPoormantracks.cxx.

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

Referenced by PndPmtPoormantracks::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  )

protectedinherited

Definition at line 51 of file PndPmtPoormantracks.cxx.

References PndPmtPoormantracks::fVerbose.

Referenced by 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 )

inlineinherited

Definition at line 20 of file PndPmtPoormantracks.h.

References a, and PndPmtPoormantracks::fDtheta.

{fDtheta=a;};

void PndPmtPoormantracks::SetNumTrk ( int  a )

inlineinherited

Definition at line 22 of file PndPmtPoormantracks.h.

References a, and PndPmtPoormantracks::fNumTrk.

{fNumTrk=a;};

void PndPmtTask::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 316 of file PndPmtTask.cxx.

                                  {
  // Get run and runtime database
  //FairRunAna* run = FairRunAna::Instance();
  //if ( ! run ) Fatal("SetParContainers", "No analysis run");
  //FairRuntimeDb* db = run->GetRuntimeDb();
  //if ( ! db ) Fatal("SetParContainers", "No runtime database");
}

void PndPmtPoormantracks::SetPid ( int  a )

inlineinherited

Definition at line 21 of file PndPmtPoormantracks.h.

References a, and PndPmtPoormantracks::fPID.

{fPID=a;};

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

inlineinherited

Definition at line 19 of file PndPmtPoormantracks.h.

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

{fPtMin=min;fPtMax=max;};

void PndPmtPoormantracks::SetSeed ( int  a )

inlineinherited

Definition at line 23 of file PndPmtPoormantracks.h.

References a, and PndPmtPoormantracks::fSeed.

{fSeed=a;};

void PndPmtPoormantracks::SetSigP ( double  a )

inlineinherited

Definition at line 18 of file PndPmtPoormantracks.h.

References a, PndPmtPoormantracks::fSigPx, PndPmtPoormantracks::fSigPy, and PndPmtPoormantracks::fSigPz.

{fSigPx=fSigPy=fSigPz=a;};

void PndPmtPoormantracks::SetSigV ( double  a )

inlineinherited

Definition at line 17 of file PndPmtPoormantracks.h.

References a, PndPmtPoormantracks::fSigVx, PndPmtPoormantracks::fSigVy, and PndPmtPoormantracks::fSigVz.

{fSigVx=fSigVy=fSigVz=a;};

void PndPmtPoormantracks::SetVerbose ( int  a )

inlineinherited

Definition at line 24 of file PndPmtPoormantracks.h.

References a, and PndPmtPoormantracks::fVerbose.

{fVerbose=a;};

void PndPmtPoormantracks::SmearMomentum ( TVector3 &  momentum )

protectedinherited

Definition at line 86 of file PndPmtPoormantracks.cxx.

References PndPmtPoormantracks::fSigPx, PndPmtPoormantracks::fSigPy, and PndPmtPoormantracks::fSigPz.

Referenced by PndPmtPoormantracks::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 )

protectedinherited

Definition at line 61 of file PndPmtPoormantracks.cxx.

References PndPmtPoormantracks::fSigVx, PndPmtPoormantracks::fSigVy, PndPmtPoormantracks::fSigVz, and PndPmtPoormantracks::fVerbose.

Referenced by PndPmtPoormantracks::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

protectedinherited

Definition at line 36 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fCharge

protectedinherited

Definition at line 46 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fDtheta

protectedinherited

Definition at line 47 of file PndPmtPoormantracks.h.

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

TH1D* PndPmtTask::fHCpu

protected

Definition at line 90 of file PndPmtTask.h.

Referenced by Finish(), and Init().

TList* PndPmtTask::fHistoList

protected

Definition at line 42 of file PndPmtTask.h.

Referenced by Finish(), Init(), and PndPmtTask().

TH1F* PndPmtTask::fHPrgPull0

protected

Definition at line 88 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHPrgPull1

protected

Definition at line 88 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHPrgPull2

protected

Definition at line 88 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHPrgPull3

protected

Definition at line 88 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHPrgPull4

protected

Definition at line 88 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxChi2Fast

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxChi2Fit

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxChi2Kin

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxChiProbFast

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxChiProbFit

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxChiProbKin

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffE

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffFitE

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffFitPhi

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffFitPX

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffFitPY

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffFitPZ

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffFitThe

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffKinE

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffKinPhi

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffKinPX

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffKinPY

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffKinPZ

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffKinThe

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffPhi

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffPX

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffPY

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffPZ

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffThe

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F* PndPmtTask::fHVtxDiffX

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffY

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxDiffZ

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrFastX

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrFastY

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrFastZ

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrFitX

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrFitY

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrFitZ

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrKinX

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrKinY

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxErrKinZ

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH2F * PndPmtTask::fHVtxFastRZ

protected

Definition at line 62 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F* PndPmtTask::fHVtxFastX

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH2F* PndPmtTask::fHVtxFastXY

protected

Definition at line 62 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxFastY

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxFastZ

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH2F * PndPmtTask::fHVtxFitRZ

protected

Definition at line 73 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F* PndPmtTask::fHVtxFitX

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH2F* PndPmtTask::fHVtxFitXY

protected

Definition at line 73 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxFitY

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxFitZ

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH2F * PndPmtTask::fHVtxKinRZ

protected

Definition at line 85 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F* PndPmtTask::fHVtxKinX

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH2F* PndPmtTask::fHVtxKinXY

protected

Definition at line 85 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxKinY

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxKinZ

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPocaEmpty

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Init().

TH2F * PndPmtTask::fHVtxPocaRZ

protected

Definition at line 54 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPocas

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F* PndPmtTask::fHVtxPocaX

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH2F* PndPmtTask::fHVtxPocaXY

protected

Definition at line 54 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPocaY

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPocaZ

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullE

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFastX

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFastY

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFastZ

protected

Definition at line 57 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFitE

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFitPX

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFitPY

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFitPZ

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFitX

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFitY

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullFitZ

protected

Definition at line 65 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullKinE

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullKinPX

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullKinPY

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullKinPZ

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullKinX

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullKinY

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullKinZ

protected

Definition at line 77 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullPocaX

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullPocaY

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullPocaZ

protected

Definition at line 50 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullPX

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullPY

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TH1F * PndPmtTask::fHVtxPullPZ

protected

Definition at line 44 of file PndPmtTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndPmtPoormantracks::fMcCands

protectedinherited

Definition at line 35 of file PndPmtPoormantracks.h.

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

int PndPmtPoormantracks::fNDF

protectedinherited

Definition at line 52 of file PndPmtPoormantracks.h.

int PndPmtPoormantracks::fNumTrk

protectedinherited

Definition at line 50 of file PndPmtPoormantracks.h.

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

TParticlePDG* PndPmtPoormantracks::fPDG

protectedinherited

Definition at line 49 of file PndPmtPoormantracks.h.

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

int PndPmtPoormantracks::fPID

protectedinherited

Definition at line 48 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fPtMax

protectedinherited

Definition at line 56 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fPtMin

protectedinherited

Definition at line 56 of file PndPmtPoormantracks.h.

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

int PndPmtPoormantracks::fSeed

protectedinherited

Definition at line 53 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fSigPx

protectedinherited

Definition at line 43 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fSigPy

protectedinherited

Definition at line 44 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fSigPz

protectedinherited

Definition at line 45 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fSigVx

protectedinherited

Definition at line 40 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fSigVy

protectedinherited

Definition at line 41 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fSigVz

protectedinherited

Definition at line 42 of file PndPmtPoormantracks.h.

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

TStopwatch PndPmtTask::fSwAll

protected

Definition at line 92 of file PndPmtTask.h.

Referenced by Finish(), and PndPmtTask().

TStopwatch PndPmtTask::fSwKin

protected

Definition at line 92 of file PndPmtTask.h.

Referenced by Exec(), Finish(), and Init().

TStopwatch PndPmtTask::fSwPMT

protected

Definition at line 92 of file PndPmtTask.h.

Referenced by Exec(), Finish(), and Init().

TStopwatch PndPmtTask::fSwPoca

protected

Definition at line 92 of file PndPmtTask.h.

Referenced by Exec(), Finish(), and Init().

TStopwatch PndPmtTask::fSwPrgfast

protected

Definition at line 92 of file PndPmtTask.h.

Referenced by Exec(), Finish(), and Init().

TStopwatch PndPmtTask::fSwPrgfull

protected

Definition at line 92 of file PndPmtTask.h.

Referenced by Exec(), Finish(), and Init().

int PndPmtPoormantracks::fVerbose

protectedinherited

Definition at line 57 of file PndPmtPoormantracks.h.

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

TVector3 PndPmtPoormantracks::fVertex

protectedinherited

Definition at line 54 of file PndPmtPoormantracks.h.

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

double PndPmtPoormantracks::fWidx

protectedinherited

Definition at line 55 of file PndPmtPoormantracks.h.

Referenced by Exec(), and Init().

double PndPmtPoormantracks::fWidy

protectedinherited

Definition at line 55 of file PndPmtPoormantracks.h.

Referenced by Exec(), and Init().

double PndPmtPoormantracks::fWidz

protectedinherited

Definition at line 55 of file PndPmtPoormantracks.h.

Referenced by Exec(), and Init().

int PndPmtPoormantracks::niter

protectedinherited

Definition at line 51 of file PndPmtPoormantracks.h.

---

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

• PndPmtTask.h
• PndPmtTask.cxx