EmcBhabhaEvent Class Reference

#include <EmcBhabhaEvent.h>

Public Types
enum	Selection {   Nothing = 0 , MDCOneProng = 1 , MDCTwoProngMatched = 2 , MDCTwoProngOneMatched = 3 ,   EmcOneProng = 4 , EmcTwoProng = 5 }

Public Member Functions
	EmcBhabhaEvent ()
	~EmcBhabhaEvent ()
EmcBhabha *	positron () const
EmcBhabha *	electron () const
EmcBhabha *	setPositron ()
EmcBhabha *	setElectron ()
void	initData ()
void	deleteData ()
void	print ()
double	getDepoMCShowerEnergy (unsigned int thetaIndex, unsigned int phiIndex, double ePeak, double beamEnergy) const
double	getDepoMCShowerEnergy_lab (double theta, double phi, unsigned int thetaIndex, unsigned int phiIndex, double ePeak, double beamEnergy) const
double	getErrorDepoMCShowerEnergy (unsigned int thetaIndex, unsigned int phiIndex, double eSigma) const
double	enLeakageTheta (double theta)
double	enLeakageThetaErr (double theta)
Hep3Vector	showerVector (EmcShower theShower)

Static Public Member Functions
static int &	selectedMDCType ()
static int &	selectedEmcType ()

Detailed Description

This class holds a Bhabha event, a electron and/or a positron. It also knows about the true and the expected deposited energies of Bhabhas. ???whether need to do the energy correction of a Bhabha due to the tilt of BES3 to the beam line with 11mrad.

Definition at line 36 of file EmcBhabhaEvent.h.

Member Enumeration Documentation

Selection

enum EmcBhabhaEvent::Selection

Enumerator
Nothing
MDCOneProng
MDCTwoProngMatched
MDCTwoProngOneMatched
EmcOneProng
EmcTwoProng

Definition at line 38 of file EmcBhabhaEvent.h.

                 {
    Nothing               = 0,
    MDCOneProng           = 1,
    MDCTwoProngMatched    = 2,
    MDCTwoProngOneMatched = 3,
    EmcOneProng           = 4,
    EmcTwoProng           = 5
  };

Constructor & Destructor Documentation

EmcBhabhaEvent()

EmcBhabhaEvent::EmcBhabhaEvent

(

)

Definition at line 47 of file EmcBhabhaEvent.cxx.

                               {
  m_positron = new EmcBhabha();
  assert( 0 != m_positron );
  m_electron = new EmcBhabha();
  assert( 0 != m_electron );
  /*
   if (m_bhaEvList.empty()) {
 
    m_bhaEvList.push_back(*m_positron);
    m_bhaEvList.push_back(*m_electron);
 
  }
  else {
    if (m_bhaEvList.size() != 2) {
      cout<<"error " << "ERROR EmcBhabhaEvent: Length of event list != 2 !"
                    << endl
                    << "      EmcBhabhaEvent: Length = "
                    << m_bhaEvList.size()
                    << endl;
    }
  }
  */
  initData();
}

~EmcBhabhaEvent()

EmcBhabhaEvent::~EmcBhabhaEvent

(

)

Definition at line 75 of file EmcBhabhaEvent.cxx.

                                {
 
  delete m_positron;
  delete m_electron;
}

Member Function Documentation

deleteData()

void EmcBhabhaEvent::deleteData

(

)

Definition at line 84 of file EmcBhabhaEvent.cxx.

{}

electron()

EmcBhabha * EmcBhabhaEvent::electron ( ) const

inline

Definition at line 56 of file EmcBhabhaEvent.h.

{ return m_electron; }

Referenced by print().

enLeakageTheta()

double EmcBhabhaEvent::enLeakageTheta

(

double

theta

)

Definition at line 269 of file EmcBhabhaEvent.cxx.

{ return 1.0; }

enLeakageThetaErr()

double EmcBhabhaEvent::enLeakageThetaErr

(

double

theta

)

Definition at line 271 of file EmcBhabhaEvent.cxx.

{ return 0.0; }

getDepoMCShowerEnergy()

double EmcBhabhaEvent::getDepoMCShowerEnergy	(	unsigned int	thetaIndex,
		unsigned int	phiIndex,
		double	ePeak,
		double	beamEnergy ) const

Definition at line 157 of file EmcBhabhaEvent.cxx.

                                                                                      {
 
  if ( !m_initialized )
  {
    cout << "error "
         << " EmcBhabhaEvent not initialized ! " << endl
         << " Is EmcLoadBhabhaData included in the path ? " << endl;
    return 0;
  }
 
  // the deposited energy
  double depoEne = 0.;
 
  EmcStructure* theEmcStruc = new EmcStructure();
 
  if ( theEmcStruc->isOutofAccep( thetaIndex, phiIndex ) )
  {
    cout << "warning "
         << "EmcBhabhaEvent: Theta " << thetaIndex << " or phi " << phiIndex
         << " out of the calorimeter acceptance !"
         << "Return 0 !" << endl;
    delete theEmcStruc;
    return 0;
  }
  else
  {
    depoEne = ePeak * beamEnergy;
 
    delete theEmcStruc;
    return depoEne;
  }
}

getDepoMCShowerEnergy_lab()

double EmcBhabhaEvent::getDepoMCShowerEnergy_lab	(	double	theta,
		double	phi,
		unsigned int	thetaIndex,
		unsigned int	phiIndex,
		double	ePeak,
		double	beamEnergy ) const

Definition at line 191 of file EmcBhabhaEvent.cxx.

                                                                            {
 
  if ( !m_initialized )
  {
    cout << "error "
         << " EmcBhabhaEvent not initialized ! " << endl
         << " Is EmcLoadBhabhaData included in the path ? " << endl;
    return 0;
  }
 
  // the deposited energy
  double depoEne = 0.;
 
  EmcStructure* theEmcStruc = new EmcStructure();
 
  if ( theEmcStruc->isOutofAccep( thetaIndex, phiIndex ) )
  {
    cout << "warning "
         << "EmcBhabhaEvent: Theta " << thetaIndex << " or phi " << phiIndex
         << " out of the calorimeter acceptance !"
         << "Return 0 !" << endl;
    delete theEmcStruc;
    return 0;
  }
  else
  {
 
    HepLorentzVector ptrk;
    ptrk.setPx( beamEnergy * sin( theta ) * cos( phi ) );
    ptrk.setPy( beamEnergy * sin( theta ) * sin( phi ) );
    ptrk.setPz( beamEnergy * cos( theta ) );
    ptrk.setE( beamEnergy );
 
    ptrk = ptrk.boost( 0.011, 0, 0 ); // boost from cms to lab
    // double the_lab = ptrk.theta();
    // double phi_lab = ptrk.phi();
    double depoEne_lab = ePeak * ptrk.e();
    delete theEmcStruc;
    return depoEne_lab;
  }
}

getErrorDepoMCShowerEnergy()

double EmcBhabhaEvent::getErrorDepoMCShowerEnergy	(	unsigned int	thetaIndex,
		unsigned int	phiIndex,
		double	eSigma ) const

Definition at line 236 of file EmcBhabhaEvent.cxx.

                                                                               {
 
  if ( !m_initialized )
  {
    cout << "error "
         << " EmcBhabhaEvent not initialized ! " << endl
         << " Is EmcLoadBhabhaData included in the path ? " << endl;
    return 0;
  }
 
  double        sig         = 0.;
  EmcStructure* theEmcStruc = new EmcStructure();
 
  if ( theEmcStruc->isOutofAccep( thetaIndex, phiIndex ) )
  {
    cout << "warning "
         << "EmcBhabhaEvent: Theta " << thetaIndex << " or phi " << phiIndex
         << " out of the calorimeter acceptance !"
         << "Return 0 !" << endl;
 
    delete theEmcStruc;
    return 0;
  }
  else
  {
    sig = eSigma;
    delete theEmcStruc;
    return sig;
  }
}

initData()

void EmcBhabhaEvent::initData

(

)

Definition at line 82 of file EmcBhabhaEvent.cxx.

{ m_initialized = true; }

Referenced by EmcBhabhaEvent().

positron()

EmcBhabha * EmcBhabhaEvent::positron ( ) const

inline

Definition at line 53 of file EmcBhabhaEvent.h.

{ return m_positron; }

Referenced by print().

print()

void EmcBhabhaEvent::print

(

)

Definition at line 87 of file EmcBhabhaEvent.cxx.

                           {
 
  if ( positron()->found() )
  {
    cout << "Positron: " << endl;
    positron()->print();
  }
  else
  {
    cout << "routine "
         << "No positron found ! " << endl;
  }
 
  if ( electron()->found() )
  {
    cout << "Electron: " << endl;
    electron()->print();
  }
  else
  {
    cout << "routine "
         << "No electron found ! " << endl;
  }
}

selectedEmcType()

int & EmcBhabhaEvent::selectedEmcType ( )

inlinestatic

Definition at line 94 of file EmcBhabhaEvent.h.

{ return m_selectedEmcType; }

selectedMDCType()

int & EmcBhabhaEvent::selectedMDCType ( )

inlinestatic

Definition at line 92 of file EmcBhabhaEvent.h.

{ return m_selectedMDCType; }

setElectron()

EmcBhabha * EmcBhabhaEvent::setElectron ( )

inline

Definition at line 60 of file EmcBhabhaEvent.h.

{ return m_electron; }

setPositron()

EmcBhabha * EmcBhabhaEvent::setPositron ( )

inline

Definition at line 58 of file EmcBhabhaEvent.h.

{ return m_positron; }

showerVector()

Hep3Vector EmcBhabhaEvent::showerVector

(

EmcShower

theShower

)

Definition at line 273 of file EmcBhabhaEvent.cxx.

                                                             {
  Hep3Vector theShowerVector( 1, 1, 1 );
  theShowerVector.setTheta( theShower.theta() );
  theShowerVector.setPhi( theShower.phi() );
  theShowerVector.setMag( theShower.energy() );
 
  return theShowerVector;
}

---

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

• EmcBhabhaEvent.h
• EmcBhabhaEvent.cxx