/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Emc/EmcCalib/EmcBhaCalib/EmcBhaCalib-00-00-34/src/EmcBhabhaEvent.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//
// Description:
//      Class EmcBhabhaEvent
//           
// Environment:
//      Software developed for the BESIII Detector at the BEPCII.
//
// Author List:
//      ChunXiu Liu          IHEP
//
//------------------------------------------------------------------------
//-----------------------
// This Class's Header --
//-----------------------

#include "EmcBhaCalib/EmcBhabhaEvent.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Vector/TwoVector.h"
using CLHEP::Hep3Vector;
using CLHEP::Hep2Vector;
using CLHEP::HepLorentzVector;
//---------------
// C Headers --
//---------------
#include <cassert>
//---------------
// C++ Headers --
//---------------
#include <fstream>
//-------------------------------
// Collaborating Class Headers --
//-------------------------------

//static elements
bool EmcBhabhaEvent::m_initialized = true;
int EmcBhabhaEvent::m_selectedMDCType = -1;
int EmcBhabhaEvent::m_selectedEmcType = -1;
//double EmcBhabhaEvent::m_eBeamEnergy = 1.548;
//double EmcBhabhaEvent::m_pBeamEnergy = 1.548;

//----------------
// Constructors --
//----------------
EmcBhabhaEvent::EmcBhabhaEvent()
{   
  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();
}

//--------------
// Destructor --
//--------------
EmcBhabhaEvent::~EmcBhabhaEvent() {

  delete m_positron;
  delete m_electron;

}

//Modifiers
void
EmcBhabhaEvent::initData() 
{
  m_initialized = true;  
}

void
EmcBhabhaEvent::deleteData() 
{
}  

//Accessor
void
EmcBhabhaEvent::print() {

  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;
    }
}


/*
double 
EmcBhabhaEvent::getDepExpShowerEnergy(unsigned int thetaIndex,
                                      unsigned int phiIndex,
                                      double ePeakEne) const {

  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  =ePeakEne;
      

      delete theEmcStruc;
      return depoEne;

    }
  
}

*/

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

  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;

    }
  
}


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

  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;

    }
  
}

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

  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; 
   
    }

}

double 
EmcBhabhaEvent::enLeakageTheta(double theta)
{
  return  1.0;
}

double 
EmcBhabhaEvent::enLeakageThetaErr(double theta)
{
  return 0.0;
}

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

---