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 37 of file EmcBhabhaEvent.h.

Member Enumeration Documentation

Selection

enum EmcBhabhaEvent::Selection

Enumerator
Nothing
MDCOneProng
MDCTwoProngMatched
MDCTwoProngOneMatched
EmcOneProng
EmcTwoProng

Definition at line 41 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 76 of file EmcBhabhaEvent.cxx.

                                {
 
  delete m_positron;
  delete m_electron;
 
}

Member Function Documentation

deleteData()

void EmcBhabhaEvent::deleteData

(

)

Definition at line 91 of file EmcBhabhaEvent.cxx.

{
}  

electron()

EmcBhabha * EmcBhabhaEvent::electron ( ) const

inline

Definition at line 57 of file EmcBhabhaEvent.h.

{return m_electron;}

Referenced by print().

enLeakageTheta()

double EmcBhabhaEvent::enLeakageTheta

(

double

theta

)

Definition at line 302 of file EmcBhabhaEvent.cxx.

{
  return  1.0;
}

enLeakageThetaErr()

double EmcBhabhaEvent::enLeakageThetaErr

(

double

theta

)

Definition at line 308 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 167 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 211 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 264 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 85 of file EmcBhabhaEvent.cxx.

{
  m_initialized = true;  
}

Referenced by EmcBhabhaEvent().

positron()

EmcBhabha * EmcBhabhaEvent::positron ( ) const

inline

Definition at line 54 of file EmcBhabhaEvent.h.

{return m_positron;}

Referenced by print().

print()

void EmcBhabhaEvent::print

(

)

Definition at line 97 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()

static int & EmcBhabhaEvent::selectedEmcType ( )

inlinestatic

Definition at line 106 of file EmcBhabhaEvent.h.

    {
      return m_selectedEmcType; 
    }

selectedMDCType()

static int & EmcBhabhaEvent::selectedMDCType ( )

inlinestatic

Definition at line 101 of file EmcBhabhaEvent.h.

    {
      return m_selectedMDCType; 
    }

setElectron()

EmcBhabha * EmcBhabhaEvent::setElectron ( )

inline

Definition at line 61 of file EmcBhabhaEvent.h.

{return m_electron;}

setPositron()

EmcBhabha * EmcBhabhaEvent::setPositron ( )

inline

Definition at line 59 of file EmcBhabhaEvent.h.

{return m_positron;}

showerVector()

Hep3Vector EmcBhabhaEvent::showerVector

(

EmcShower

theShower

)

Definition at line 314 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