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 "EmcBhabhaEvent.h"
#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Vector/TwoVector.h"
using CLHEP::Hep2Vector;
using CLHEP::Hep3Vector;
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;
}