MucPID Class Reference

#include <MucPID.h>

Inheritance diagram for MucPID:

Public Member Functions
void	init ()
void	calculate ()
bool	IsPidInfoValid () const
double	chi (int n) const
double	prob (int n) const
int	ndof () const
double	hits () const
double	depth () const
double	chi2 () const
double	distance () const
double	delta_phi () const
double	val_muc1 () const
int	neuronPID () const
void	init ()
void	calculate ()
bool	IsPidInfoValid () const
double	chi (int n) const
double	prob (int n) const
int	ndof () const
double	hits () const
double	depth () const
double	chi2 () const
double	distance () const
double	delta_phi () const
double	val_muc1 () const
int	neuronPID () const
void	init ()
void	calculate ()
bool	IsPidInfoValid () const
double	chi (int n) const
double	prob (int n) const
int	ndof () const
double	hits () const
double	depth () const
double	chi2 () const
double	distance () const
double	delta_phi () const
double	val_muc1 () const
int	neuronPID () const
Public Member Functions inherited from ParticleIDBase
	ParticleIDBase ()
virtual	~ParticleIDBase ()
EvtRecTrack *	PidTrk () const
void	setRecTrack (EvtRecTrack *trk)
double	chiMinCut () const
void	setChiMinCut (const double chi=4)
double	chiMaxCut () const
void	setChiMaxCut (const double chi=6)
double	pdfMinSigmaCut () const
void	setPdfMinSigmaCut (const double pdf=4)
double	getRunNo () const
void	setRunNo (const double runh=8093)
double	p ()
double	pt ()
double	charge ()
double	xmass (int n)
double	velc ()
double	probCalculate (double chi2, int n)
double	pdfCalculate (double offset, double sigma)
double	interpolation (double *x, double *y, double x1)
double	pol2 (double x, double *par)
double	pol3 (double x, double *par)
double	pol4 (double x, double *par)
int	useDedx () const
int	useTof () const
int	useTof1 () const
int	useTof2 () const
int	useTofE () const
int	useTofQ () const
int	useTofC () const
int	useTofCorr () const
int	useEmc () const
int	useMuc () const
int	onlyPionKaon () const
int	onlyPionKaonProton () const
int	onlyPionKaonElectron () const
int	all () const
int	onlyElectron () const
int	onlyMuon () const
int	onlyPion () const
int	onlyKaon () const
int	onlyProton () const
int	methodLikelihood () const
int	methodProbability () const
int	methodNeuronNetwork () const
int	dedxValid () const
int	tofValid () const
int	tofeValid () const
int	tofqValid () const
int	tofcValid () const
int	tofcorrValid () const
int	emcValid () const
int	mucValid () const
void	set_path (const char *s_path=0)
void	set_path (std::string s_path)
	ParticleIDBase ()
virtual	~ParticleIDBase ()
EvtRecTrack *	PidTrk () const
void	setRecTrack (EvtRecTrack *trk)
double	chiMinCut () const
void	setChiMinCut (const double chi=4)
double	chiMaxCut () const
void	setChiMaxCut (const double chi=6)
double	pdfMinSigmaCut () const
void	setPdfMinSigmaCut (const double pdf=4)
double	getRunNo () const
void	setRunNo (const double runh=8093)
double	p ()
double	pt ()
double	charge ()
double	xmass (int n)
double	velc ()
double	probCalculate (double chi2, int n)
double	pdfCalculate (double offset, double sigma)
double	interpolation (double *x, double *y, double x1)
double	pol2 (double x, double *par)
double	pol3 (double x, double *par)
double	pol4 (double x, double *par)
int	useDedx () const
int	useTof () const
int	useTof1 () const
int	useTof2 () const
int	useTofE () const
int	useTofQ () const
int	useTofC () const
int	useTofCorr () const
int	useEmc () const
int	useMuc () const
int	onlyPionKaon () const
int	onlyPionKaonProton () const
int	onlyPionKaonElectron () const
int	all () const
int	onlyElectron () const
int	onlyMuon () const
int	onlyPion () const
int	onlyKaon () const
int	onlyProton () const
int	methodLikelihood () const
int	methodProbability () const
int	methodNeuronNetwork () const
int	dedxValid () const
int	tofValid () const
int	tofeValid () const
int	tofqValid () const
int	tofcValid () const
int	tofcorrValid () const
int	emcValid () const
int	mucValid () const
void	set_path (const char *s_path=0)
void	set_path (std::string s_path)
	ParticleIDBase ()
virtual	~ParticleIDBase ()
EvtRecTrack *	PidTrk () const
void	setRecTrack (EvtRecTrack *trk)
double	chiMinCut () const
void	setChiMinCut (const double chi=4)
double	chiMaxCut () const
void	setChiMaxCut (const double chi=6)
double	pdfMinSigmaCut () const
void	setPdfMinSigmaCut (const double pdf=4)
double	getRunNo () const
void	setRunNo (const double runh=8093)
double	p ()
double	pt ()
double	charge ()
double	xmass (int n)
double	velc ()
double	probCalculate (double chi2, int n)
double	pdfCalculate (double offset, double sigma)
double	interpolation (double *x, double *y, double x1)
double	pol2 (double x, double *par)
double	pol3 (double x, double *par)
double	pol4 (double x, double *par)
int	useDedx () const
int	useTof () const
int	useTof1 () const
int	useTof2 () const
int	useTofE () const
int	useTofQ () const
int	useTofC () const
int	useTofCorr () const
int	useEmc () const
int	useMuc () const
int	onlyPionKaon () const
int	onlyPionKaonProton () const
int	onlyPionKaonElectron () const
int	all () const
int	onlyElectron () const
int	onlyMuon () const
int	onlyPion () const
int	onlyKaon () const
int	onlyProton () const
int	methodLikelihood () const
int	methodProbability () const
int	methodNeuronNetwork () const
int	dedxValid () const
int	tofValid () const
int	tofeValid () const
int	tofqValid () const
int	tofcValid () const
int	tofcorrValid () const
int	emcValid () const
int	mucValid () const
void	set_path (const char *s_path=0)
void	set_path (std::string s_path)

Static Public Member Functions
static MucPID *	instance ()
static MucPID *	instance ()
static MucPID *	instance ()

Protected Member Functions
int	particleIDCalculation ()
int	neuronPIDCalculation ()
int	LikelihoodCalculation ()
int	particleIDCalculation ()
int	neuronPIDCalculation ()
int	LikelihoodCalculation ()
int	particleIDCalculation ()
int	neuronPIDCalculation ()
int	LikelihoodCalculation ()

Additional Inherited Members
Static Protected Attributes inherited from ParticleIDBase
static std::string	path = ""

Detailed Description

Definition at line 19 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

Member Function Documentation

calculate() [1/3]

void MucPID::calculate ( )

virtual

Implements ParticleIDBase.

Definition at line 103 of file MucPID.cxx.

                       {
  if ( particleIDCalculation() == 0 ) m_ndof = 1;
}

calculate() [2/3]

void MucPID::calculate ( )

virtual

Implements ParticleIDBase.

calculate() [3/3]

void MucPID::calculate ( )

virtual

Implements ParticleIDBase.

chi() [1/3]

double MucPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 28 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_chi[n]; }

chi() [2/3]

double MucPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 28 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_chi[n]; }

chi() [3/3]

double MucPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 28 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_chi[n]; }

chi2() [1/3]

double MucPID::chi2 ( ) const

inline

Definition at line 33 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_chi2; }

Referenced by init().

chi2() [2/3]

double MucPID::chi2 ( ) const

inline

Definition at line 33 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_chi2; }

chi2() [3/3]

double MucPID::chi2 ( ) const

inline

Definition at line 33 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_chi2; }

delta_phi() [1/3]

double MucPID::delta_phi ( ) const

inline

Definition at line 35 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_muc_delta_phi; }

delta_phi() [2/3]

double MucPID::delta_phi ( ) const

inline

Definition at line 35 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_muc_delta_phi; }

delta_phi() [3/3]

double MucPID::delta_phi ( ) const

inline

Definition at line 35 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_muc_delta_phi; }

depth() [1/3]

double MucPID::depth ( ) const

inline

Definition at line 32 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_depth; }

Referenced by init().

depth() [2/3]

double MucPID::depth ( ) const

inline

Definition at line 32 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_depth; }

depth() [3/3]

double MucPID::depth ( ) const

inline

Definition at line 32 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_depth; }

distance() [1/3]

double MucPID::distance ( ) const

inline

Definition at line 34 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_distance; }

Referenced by init().

distance() [2/3]

double MucPID::distance ( ) const

inline

Definition at line 34 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_distance; }

distance() [3/3]

double MucPID::distance ( ) const

inline

Definition at line 34 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_distance; }

hits() [1/3]

double MucPID::hits ( ) const

inline

Definition at line 31 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_hits; }

Referenced by init().

hits() [2/3]

double MucPID::hits ( ) const

inline

Definition at line 31 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_hits; }

hits() [3/3]

double MucPID::hits ( ) const

inline

Definition at line 31 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_hits; }

init() [1/3]

void MucPID::init ( )

virtual

Implements ParticleIDBase.

Definition at line 65 of file MucPID.cxx.

                  {
  for ( int i = 0; i < 5; i++ )
  {
    m_chi[i]  = 99.0;
    m_prob[i] = -1.0;
  }
  m_chimin  = 99.;
  m_ndof    = 0;
  m_hits    = -99;
  m_depth   = -999;
  m_val_muc = -99;
 
  std::string neural = path + "/share/neural.root";
  std::string muc    = path + "/share/muc.txt";
  double      ptrk, phi, theta, depth, hits, chi2, distance, muc_delta_phi;
  int         type;
  if ( !m_trainTree_muc )
  {
    m_trainTree_muc = new TTree( "m_trainTree_muc", "m_trainTree_muc" );
    m_trainTree_muc->Branch( "ptrk", &ptrk, "ptrk/D" );
    m_trainTree_muc->Branch( "phi", &phi, "phi/D" );
    m_trainTree_muc->Branch( "theta", &theta, "theta/D" );
    m_trainTree_muc->Branch( "depth", &depth, "depth/D" );
    m_trainTree_muc->Branch( "hits", &hits, "hits/D" );
    m_trainTree_muc->Branch( "chi2", &chi2, "chi2/D" );
    m_trainTree_muc->Branch( "distance", &distance, "distance/D" );
    m_trainTree_muc->Branch( "muc_delta_phi", &muc_delta_phi, "muc_delta_phi/D" );
    m_trainTree_muc->Branch( "type", &type, "type/I" );
  }
  if ( !m_mlp_muc )
  {
    m_mlp_muc = new TMultiLayerPerceptron(
        "@ptrk,@phi,@theta,@depth,@hits,@chi2,@distance,@muc_delta_phi:16:type",
        m_trainTree_muc );
    m_mlp_muc->LoadWeights( muc.c_str() );
  }
}

init() [2/3]

void MucPID::init ( )

virtual

Implements ParticleIDBase.

init() [3/3]

void MucPID::init ( )

virtual

Implements ParticleIDBase.

instance() [1/3]

MucPID * MucPID::instance ( )

static

Definition at line 22 of file MucPID.cxx.

                         {
 
  if ( !m_pointer ) m_pointer = new MucPID();
  return m_pointer;
}

Referenced by ParticleID::init().

instance() [2/3]

MucPID * MucPID::instance ( )

static

instance() [3/3]

MucPID * MucPID::instance ( )

static

IsPidInfoValid() [1/3]

bool MucPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 27 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return ( m_ndof > 0 ); }

IsPidInfoValid() [2/3]

bool MucPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 27 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return ( m_ndof > 0 ); }

IsPidInfoValid() [3/3]

bool MucPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 27 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return ( m_ndof > 0 ); }

LikelihoodCalculation() [1/3]

int MucPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 43 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return -1; }

LikelihoodCalculation() [2/3]

int MucPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 43 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return -1; }

LikelihoodCalculation() [3/3]

int MucPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 43 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return -1; }

ndof() [1/3]

int MucPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 30 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_ndof; }

ndof() [2/3]

int MucPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_ndof; }

ndof() [3/3]

int MucPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_ndof; }

neuronPID() [1/3]

int MucPID::neuronPID ( ) const

inline

Definition at line 38 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return -1; }

neuronPID() [2/3]

int MucPID::neuronPID ( ) const

inline

Definition at line 38 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return -1; }

neuronPID() [3/3]

int MucPID::neuronPID ( ) const

inline

Definition at line 38 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return -1; }

neuronPIDCalculation() [1/3]

int MucPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 42 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return -1; }

neuronPIDCalculation() [2/3]

int MucPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 42 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return -1; }

neuronPIDCalculation() [3/3]

int MucPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 42 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return -1; }

particleIDCalculation() [1/3]

int MucPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 107 of file MucPID.cxx.

                                  {
  int          irc    = -1;
  EvtRecTrack* recTrk = PidTrk();
  if ( !( recTrk->isMdcTrackValid() ) ) return irc;
  RecMdcTrack* mdcTrk = recTrk->mdcTrack();
 
  m_depth         = -99;
  m_hits          = -99;
  m_chi2          = -99;
  m_distance      = -99;
  m_muc_delta_phi = -99;
 
  double ptrk   = mdcTrk->p();
  double m_ptrk = ptrk;
  double m_pt   = mdcTrk->pxy();
  double phi    = mdcTrk->phi();
  double theta  = mdcTrk->theta();
  double cost   = cos( mdcTrk->theta() );
  if ( ptrk < 0.5 ) return irc;
  if ( fabs( cost ) > 0.83 ) return irc;
  if ( !( recTrk->isExtTrackValid() ) ) return irc;
  RecExtTrack* extTrk = recTrk->extTrack();
  if ( extTrk->mucVolumeNumber() == -1 ) return irc;
  if ( !( recTrk->isMucTrackValid() ) ) return irc;
  RecMucTrack* mucTrk = recTrk->mucTrack();
 
  //  if(mucTrk->maxHitsInLayer()< 0) return irc;
  if ( mucTrk->depth() > 100000 ) return irc;
 
  m_hits     = mucTrk->maxHitsInLayer();
  m_depth    = mucTrk->depth();
  m_distance = mucTrk->distance();
  m_chi2     = mucTrk->chi2();
  /* Hep3Vector phi_muc;
   phi_muc.set(mucTrk->xPos(),mucTrk->yPos(),0);
   Hep3Vector phi_mdc;
   phi_mdc.set(mdcTrk->px(),mdcTrk->py(),0);
   m_muc_delta_phi = phi_muc.angle(phi_mdc);
   if(m_muc_delta_phi<0) m_muc_delta_phi = -m_muc_delta_phi;        */
  m_muc_delta_phi = mucTrk->deltaPhi();
  m_muc_delta_phi = 3.14159 - m_muc_delta_phi;
  theta           = cos( theta );
  params_muc1[0]  = m_ptrk;
  params_muc1[1]  = phi;
  params_muc1[2]  = theta;
  params_muc1[3]  = m_depth;
  params_muc1[4]  = m_hits;
  params_muc1[5]  = m_chi2;
  params_muc1[6]  = m_distance;
  params_muc1[7]  = m_muc_delta_phi;
 
  m_val_muc = m_mlp_muc->Evaluate( 0, params_muc1 );
  if ( m_pt < 0 ) m_pt = -m_pt;
  int pindex = int( ( m_ptrk - 0.5 ) / 0.1 );
  int bindex = int( ( m_val_muc - 0.5 ) / 0.01 );
  if ( bindex > 300 || bindex < 0 ) return irc;
  if ( pindex > 11 ) pindex = 11;
  if ( pindex < 0 ) pindex = 0;
  m_prob[0] = m_p_h[pindex][bindex];
  ;
  m_prob[1] = m_m_h[pindex][bindex];
  m_prob[2] = m_p_h[pindex][bindex];
  m_prob[3] = m_p_h[pindex][bindex];
  m_prob[4] = m_p_h[pindex][bindex];
  for ( int i = 0; i < 5; i++ )
  {
    if ( m_prob[i] == 0 ) m_prob[i] = 0.001;
  }
  float ppp[5];
  for ( int i = 0; i < 5; i++ ) { ppp[i] = 0.0; }
 
  for ( int j = 0; j < bindex; j++ )
  {
    ppp[1] += m_m_h[pindex][j] * 0.01;
    ppp[2] += m_p_h[pindex][j] * 0.01;
  }
  if ( ppp[1] > 0 && ppp[1] < 1 ) CALG( ppp[1], m_chi[1] );
  if ( ppp[2] > 0 && ppp[2] < 1 ) CALG( ppp[2], m_chi[2] );
  if ( ppp[1] <= 0 || ppp[1] >= 1 ) m_chi[1] = -99;
  if ( ppp[2] <= 0 || ppp[2] >= 1 ) m_chi[2] = -99;
 
  m_chi[3] = m_chi[2];
  m_chi[4] = m_chi[2];
  m_chi[0] = m_chi[2];
  m_ndof   = 1;
  irc      = 0;
  return irc;
}

Referenced by calculate().

particleIDCalculation() [2/3]

int MucPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

particleIDCalculation() [3/3]

int MucPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

prob() [1/3]

double MucPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 29 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_prob[n]; }

prob() [2/3]

double MucPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 29 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_prob[n]; }

prob() [3/3]

double MucPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 29 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_prob[n]; }

val_muc1() [1/3]

double MucPID::val_muc1 ( ) const

inline

Definition at line 36 of file Analysis/ParticleID/include/ParticleID/MucPID.h.

{ return m_val_muc; }

val_muc1() [2/3]

double MucPID::val_muc1 ( ) const

inline

Definition at line 36 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h.

{ return m_val_muc; }

val_muc1() [3/3]

double MucPID::val_muc1 ( ) const

inline

Definition at line 36 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h.

{ return m_val_muc; }

---

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

• Analysis/ParticleID/include/ParticleID/MucPID.h
• InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/MucPID.h
• InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/MucPID.h
• MucPID.cxx