EmcRecDigit2Hit Class Reference

#include <EmcRecDigit2Hit.h>

Public Member Functions
	EmcRecDigit2Hit ()
	~EmcRecDigit2Hit ()
void	Convert (const RecEmcDigitMap &aDigitMap, RecEmcHitMap &aHitMap)
void	Output (const RecEmcHitMap &aHitMap) const
void	OutputEndcap (const RecEmcHitMap &aHitMap, const unsigned int module_ew) const
void	SetAlgName (const string &name)

Detailed Description

Definition at line 12 of file EmcRecDigit2Hit.h.

Constructor & Destructor Documentation

EmcRecDigit2Hit()

EmcRecDigit2Hit::EmcRecDigit2Hit

(

)

Definition at line 22 of file EmcRecDigit2Hit.cxx.

{}

~EmcRecDigit2Hit()

EmcRecDigit2Hit::~EmcRecDigit2Hit

(

)

Definition at line 25 of file EmcRecDigit2Hit.cxx.

{}

Member Function Documentation

Convert()

void EmcRecDigit2Hit::Convert	(	const RecEmcDigitMap &	aDigitMap,
		RecEmcHitMap &	aHitMap )

Definition at line 28 of file EmcRecDigit2Hit.cxx.

{
  RecEmcDigitMap::const_iterator ciDigitMap;
  RecEmcHit aHit;
 
  EmcRecParameter& Para=EmcRecParameter::GetInstance();
 
  StatusCode sc;
  // Get EmcCalibConstSvc.
  IEmcCalibConstSvc *emcCalibConstSvc;
  sc= Gaudi::svcLocator()->service("EmcCalibConstSvc", emcCalibConstSvc);
  if(sc != StatusCode::SUCCESS) {
    //cout << "EmcRecDigit2Hit Error: Can't get EmcCalibConstSvc." << endl;
  }
 
  // Get RawDataProviderSvc
  IRawDataProviderSvc *rawDataProviderSvc;
  // tianhl for mt
  std::string rawDataProviderSvc_name("RawDataProviderSvc");
  if(isGaudiThreaded(m_algName)) {
    rawDataProviderSvc_name += getGaudiThreadIDfromName(m_algName);
  }
  // tianhl for mt
  sc = Gaudi::svcLocator()->service(rawDataProviderSvc_name.c_str(), rawDataProviderSvc);
  if(sc != StatusCode::SUCCESS) {
    //cout << "EmcRecDigit2Hit Error: Can't get RawDataProviderSvc." << endl;
  }
 
  for(ciDigitMap=aDigitMap.begin();
      ciDigitMap!=aDigitMap.end();
      ++ciDigitMap) {
 
    const RecEmcDigit &aDigit=ciDigitMap->second;
 
    aHit.CellId(aDigit.CellId());
 
    double ein,ecorr,eout;
    ein=aDigit.ADC();
 
    //Times the calibration constant in offline mode
    if(!(rawDataProviderSvc->isOnlineMode()) && Para.DigiCalib()) {
 
      unsigned int partId = EmcID::barrel_ec(aDigit.CellId());
      unsigned int theta = EmcID::theta_module(aDigit.CellId());
      unsigned int phi = EmcID::phi_module(aDigit.CellId());
 
      int index = emcCalibConstSvc->getIndex(partId,theta,phi);
 
      //liucx
 
      double adc2e = emcCalibConstSvc->getDigiCalibConst(index);  
 
      if(Para.DataMode()==0) { 
    if(Para.ElecSaturation()==1){
      double emaxData = emcCalibConstSvc->getCrystalEmaxData(index);
      if (emaxData>0) {
        //  cout<<adc2e<<"\t"<<emaxData/2.5<<endl;
        adc2e=emaxData/2.5;
      }
    }
      }
 
      ein *= adc2e ;
      //liucx
 
      //liucx ein *= emcCalibConstSvc->getDigiCalibConst(index);
 
    }
 
    // Correct electronic noise bias. Not used now.
    if(ein<Para.ElecBias(4)) {
      ecorr=Para.ElecBias(0)*log(ein/Para.ElecBias(1));     //noise=0.5
    } else {
      ecorr=Para.ElecBias(2)+Para.ElecBias(3)*ein;
    }
    eout=ein;
    aHit.Energy(eout/GeV);
 
    aHit.Time(aDigit.TDC());
    // For some reason, there are 0 energy hit. 
    // Avoid this to happen, dangerous!
    if(aHit.getEnergy()>=Para.ElectronicsNoiseLevel()){// &&
      //(aHit.time()>=Para.TimeMin()) && 
      //(aHit.time()<=Para.TimeMax()))
      aHitMap[aHit.getCellId()]=aHit;
    }
  }
 
}

Referenced by EmcRec::execute().

Output()

void EmcRecDigit2Hit::Output

(

const RecEmcHitMap &

aHitMap

)

const

Definition at line 119 of file EmcRecDigit2Hit.cxx.

{
  unsigned int module,theta,phi;
  RecEmcID id;
  RecEmcHitMap::const_iterator ci_HitMap;
  RecEmcEnergy e;
  int ie;
 
  //cout<<"Output a hitmap:"<<endl;
 
  module=EmcID::EmcID::getBARREL();
 
  //cout<<"Barrel:"<<endl;
  //cout<<"+--------------------------------------------+"<<endl;
  for(phi=EmcID::getPHI_BARREL_MAX();
      phi!=EmcID::getPHI_BARREL_MIN()-1;
      --phi) {
    //cout<<"|";
    for(theta=EmcID::getTHETA_BARREL_MIN();
        theta<=EmcID::getTHETA_BARREL_MAX();
        ++theta) {
      id=EmcID::crystal_id(module,theta,phi);
      ci_HitMap=aHitMap.find(id);
      if(ci_HitMap!=aHitMap.end()) {
        e=ci_HitMap->second.getEnergy();
        ie=(int)(e*10);
        if(ie>9) { 
          ie=9; 
        }
        //cout<<ie;
      } else {
        //cout<<" ";
      }
    }
    //cout<<"|"<<endl;
  }
  //cout<<"+--------------------------------------------+"<<endl;
 
  module=EmcID::EmcID::getENDCAP_EAST();
  OutputEndcap(aHitMap,module);
  module=EmcID::EmcID::getENDCAP_WEST();
  OutputEndcap(aHitMap,module);
}

Referenced by EmcRec::execute().

OutputEndcap()

void EmcRecDigit2Hit::OutputEndcap	(	const RecEmcHitMap &	aHitMap,
		const unsigned int	module_ew ) const

Definition at line 163 of file EmcRecDigit2Hit.cxx.

{
  unsigned int module,theta,phi;
  RecEmcID id;
  RecEmcHitMap::const_iterator ci_HitMap;
  RecEmcEnergy e;
  int ie;
 
  module=module_ew;
  if(module==EmcID::EmcID::getENDCAP_EAST()){
    //cout<<"East endcap:"<<endl;
  }
  else if(module==EmcID::EmcID::getENDCAP_WEST()){
    //cout<<"West endcap:"<<endl;
  }
 
  //cout<<"+--------------------------------------------------------+"<<endl;
  for(theta=EmcID::getTHETA_ENDCAP_MAX();
      theta!=EmcID::getTHETA_ENDCAP_MIN()-1;
      theta--){
    //cout<<"|";
    for(phi=(EmcID::getPHI_ENDCAP_MAX(theta)-1)/2;
        phi!=EmcID::getPHI_ENDCAP_MIN()-1;
        phi--) {
      id=EmcID::crystal_id(module,theta,phi);
      ci_HitMap=aHitMap.find(id);
      if(ci_HitMap!=aHitMap.end()) {
        e=ci_HitMap->second.getEnergy();
        ie=(int)(e*10);
        if(ie>9) {
          ie=9;
        }
        if((theta<4) && 
            (phi%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==(EmcID::getPHI_ENDCAP_MAX(theta)+1)/16-1)){
          //cout<<"*";
        }
        //cout<<ie;
        if((theta==0||theta==1)&&(phi%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0)){
          //cout<<"*";
        }
      }
      else{
        if((theta<4) &&
            (phi%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==(EmcID::getPHI_ENDCAP_MAX(theta)+1)/16-1)){
          //cout<<" ";
        }
        //cout<<" ";
        if((theta==0||theta==1)&&(phi%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0)){
          //cout<<" ";
        }
      }
      if(phi%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0){
        //cout<<"|";
      }
      if(phi==(EmcID::getPHI_ENDCAP_MAX(theta)+1)/4){
        //cout<<"|";
      }
    }
    //cout<<endl;
  }
  //cout<<"+--------------------------------------------------------+"<<endl;
 
  for(theta=EmcID::getTHETA_ENDCAP_MIN();
      theta!=EmcID::getTHETA_ENDCAP_MAX()+1;
      theta++){
    //cout<<"|";
    for(phi=(EmcID::getPHI_ENDCAP_MAX(theta)+1)/2;
        phi!=EmcID::getPHI_ENDCAP_MAX(theta)+1;
        phi++) {
      id=EmcID::crystal_id(module,theta,phi);
      ci_HitMap=aHitMap.find(id);
      if(ci_HitMap!=aHitMap.end()) {
        e=ci_HitMap->second.getEnergy();
        ie=(int)(e*10);
        if(ie>9) {
          ie=9;
        }
        if((theta<4) && 
            (phi%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0)){
          //cout<<"*";
        }
        //cout<<ie;
        if((theta==0||theta==1)&&((phi+1)%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0)){
          //cout<<"*";
        }
      }
      else{
        if((theta<4) &&
            (phi%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0)){
          //cout<<" ";
        }
        //cout<<" ";
        if((theta==0||theta==1)&&((phi+1)%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0)){
          //cout<<" ";
        }
      } 
      if((phi+1)%((EmcID::getPHI_ENDCAP_MAX(theta)+1)/16)==0){
        //cout<<"|";
      }
      if(phi==3*(EmcID::getPHI_ENDCAP_MAX(theta)+1)/4-1){
        //cout<<"|";
      }
    }
    //cout<<endl;
  }
  //cout<<"+--------------------------------------------------------+"<<endl;
 
}

Referenced by Output().

SetAlgName()

void EmcRecDigit2Hit::SetAlgName ( const string & name )

inline

Definition at line 24 of file EmcRecDigit2Hit.h.

{ m_algName=name; }

Referenced by EmcRec::initialize().

---

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

• EmcRecDigit2Hit.h
• EmcRecDigit2Hit.cxx