BesEmcDigitization Class Reference

#include <BesEmcDigitization.hh>

Public Member Functions
	BesEmcDigitization ()
	~BesEmcDigitization ()
G4bool	Digitize (const G4Event *, G4double)
G4double	GetEmcSignal (G4int, G4int)
G4int	GetNSignal ()
G4double *	GetESignal ()
G4int *	GetNThetaSignal ()
G4int *	GetNPhiSignal ()
void	Print ()

Detailed Description

Definition at line 16 of file BesEmcDigitization.hh.

Constructor & Destructor Documentation

BesEmcDigitization()

BesEmcDigitization::BesEmcDigitization

(

)

Definition at line 22 of file BesEmcDigitization.cc.

  :calorimeterCollID(-1),fNHits(0),fTotECrystal(0),fTotLCrystal(0),
   fNCrystalSignal(0)
{
  for(int i=0;i<50;i++)
    {for(int j=0;j<150;j++)
      fSingleECrystal[i][j]=0;
    }
  G4SDManager * SDman = G4SDManager::GetSDMpointer();
  calorimeterCollID= SDman->GetCollectionID("CalCollection");
}

~BesEmcDigitization()

BesEmcDigitization::~BesEmcDigitization

(

)

Definition at line 34 of file BesEmcDigitization.cc.

{
  delete[]fEnergySignal;
  delete[]fNThetaSignal;
  delete[]fNPhiSignal;
}

Member Function Documentation

Digitize()

G4bool BesEmcDigitization::Digitize	(	const G4Event *	evt,
		G4double	ecut )

Definition at line 41 of file BesEmcDigitization.cc.

{
  G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
  BesEmcHitsCollection* CHC = 0;
 
  if (HCE)
    CHC = (BesEmcHitsCollection*)(HCE->GetHC(calorimeterCollID));
  else
    return false;
 
  G4int i=0,j=0;
  if (CHC)
    {
      fNHits = CHC->entries();
      if(fNHits)
    for (i=0;i<fNHits;i++)
      {
        fTotECrystal += (*CHC)[i]->GetEdepCrystal(); 
        fTotLCrystal += (*CHC)[i]->GetTrakCrystal();
        G4int nTheta=45-(*CHC)[i]->GetNumThetaCrystal();
        G4int nPhi  =(*CHC)[i]->GetNumPhiCrystal();
        if(nTheta>0&&nTheta<50&&nPhi>0&&nPhi<150)//waiting upgrade
          {
        fSingleECrystal[nTheta][nPhi]+=(*CHC)[i]->GetEdepCrystal();
        //G4cout << "EMC Hit in crystal(" 
        //       << nTheta << "," <<nPhi<< ")" << G4endl;
          }
        else
          {
        G4cout << "EMC Hit out range:"
               << nTheta <<','<< nPhi << G4endl;
          }
      }
    }
  else
    return false;
  
  for(i=1;i<50;i++)
    {for(j=1;j<150;j++)
      {
    if(fSingleECrystal[i][j]>ecut)
      {
        fNCrystalSignal++;
      }
      }
    }
  if((fEnergySignal=new G4double[fNCrystalSignal+1])==NULL||
     (fNThetaSignal=new G4int[fNCrystalSignal+1])==NULL||
     (fNPhiSignal=new G4int[fNCrystalSignal+1])==NULL)
    {
      //G4Exception("BesEmcDigitization:can't allocate memory for Signal");
      G4cout<<"BesEmcDigitization:can't allocate memory for Signal"<<G4endl;
      exit(-1);
    }
  G4int n=0;
  for(i=1;i<50;i++)
    {for(j=1;j<150;j++)
      {
    if(fSingleECrystal[i][j]>ecut)
      {
        n++;
        fEnergySignal[n]=fSingleECrystal[i][j];
        fNThetaSignal[n]=45-i;
        fNPhiSignal  [n]=j;
      }
      }
    }
  
  return true;
}

GetEmcSignal()

G4double BesEmcDigitization::GetEmcSignal	(	G4int	nTheta,
		G4int	nPhi )

Definition at line 112 of file BesEmcDigitization.cc.

{
  if(nTheta>0&&nTheta<50&&nPhi>0&&nPhi<150)
    return fSingleECrystal[nTheta][nPhi];
  else
    return 0;
}

GetESignal()

G4double * BesEmcDigitization::GetESignal ( )

inline

Definition at line 27 of file BesEmcDigitization.hh.

{return fEnergySignal;};

GetNPhiSignal()

G4int * BesEmcDigitization::GetNPhiSignal ( )

inline

Definition at line 29 of file BesEmcDigitization.hh.

{return fNPhiSignal;};

GetNSignal()

G4int BesEmcDigitization::GetNSignal ( )

inline

Definition at line 25 of file BesEmcDigitization.hh.

{return fNCrystalSignal;};//amount of crystals

GetNThetaSignal()

G4int * BesEmcDigitization::GetNThetaSignal ( )

inline

Definition at line 28 of file BesEmcDigitization.hh.

{return fNThetaSignal;};

Print()

void BesEmcDigitization::Print

(

)

Definition at line 120 of file BesEmcDigitization.cc.

{
  G4cout
    << "   Absorber: total energy: "
    << fTotECrystal/GeV << " GeV" << G4endl
    << "       total track length: "
    << fTotLCrystal/cm  << " cm"  << G4endl;
      
  G4cout << "\n     " << fNHits
     << " hits are stored in BesEmcHitsCollection." << G4endl;
}

---

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

• BesEmcDigitization.hh
• BesEmcDigitization.cc