BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EmcRecShowerEnergy Class Reference

#include <EmcRecShowerEnergy.h>

Public Member Functions

void Energy (RecEmcShower &aShower)
void E3x3 (RecEmcShower &aShower)
void E5x5 (RecEmcShower &aShower)
void ELepton (RecEmcShower &aShower)
void EAll (RecEmcShower &aShower)
RecEmcEnergy ECorrection (const RecEmcEnergy eIn)
RecEmcEnergy ECorrTheta (const RecEmcEnergy eIn, const RecEmcID &id)
void Energy (RecEmcShower &aShower)
void E3x3 (RecEmcShower &aShower)
void E5x5 (RecEmcShower &aShower)
void ELepton (RecEmcShower &aShower)
void EAll (RecEmcShower &aShower)
RecEmcEnergy ECorrection (const RecEmcEnergy eIn)
RecEmcEnergy ECorrTheta (const RecEmcEnergy eIn, const RecEmcID &id)
void Energy (RecEmcShower &aShower)
void E3x3 (RecEmcShower &aShower)
void E5x5 (RecEmcShower &aShower)
void ELepton (RecEmcShower &aShower)
void EAll (RecEmcShower &aShower)
RecEmcEnergy ECorrection (const RecEmcEnergy eIn)
RecEmcEnergy ECorrTheta (const RecEmcEnergy eIn, const RecEmcID &id)

Detailed Description

Definition at line 11 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcRec/EmcRecShowerEnergy.h.

Member Function Documentation

◆ E3x3() [1/3]

void EmcRecShowerEnergy::E3x3 ( RecEmcShower & aShower)

◆ E3x3() [2/3]

void EmcRecShowerEnergy::E3x3 ( RecEmcShower & aShower)

◆ E3x3() [3/3]

void EmcRecShowerEnergy::E3x3 ( RecEmcShower & aShower)

◆ E5x5() [1/3]

void EmcRecShowerEnergy::E5x5 ( RecEmcShower & aShower)

◆ E5x5() [2/3]

void EmcRecShowerEnergy::E5x5 ( RecEmcShower & aShower)

◆ E5x5() [3/3]

void EmcRecShowerEnergy::E5x5 ( RecEmcShower & aShower)

◆ EAll() [1/3]

void EmcRecShowerEnergy::EAll ( RecEmcShower & aShower)

◆ EAll() [2/3]

void EmcRecShowerEnergy::EAll ( RecEmcShower & aShower)

◆ EAll() [3/3]

void EmcRecShowerEnergy::EAll ( RecEmcShower & aShower)

◆ ECorrection() [1/3]

RecEmcEnergy EmcRecShowerEnergy::ECorrection ( const RecEmcEnergy eIn)

Definition at line 128 of file EmcRecShowerEnergy.cxx.

128 {
129 if ( eIn > 3. ) return eIn;
130
131 EmcRecParameter& Para = EmcRecParameter::GetInstance();
132
133 RecEmcEnergy eOut = 0;
134 double par[4];
135 for ( int i = 0; i < 4; i++ ) { par[i] = Para.ECorr( i ); }
136
137 eOut = eIn / ( par[0] + par[1] * eIn + par[2] * eIn * eIn + par[3] * eIn * eIn * eIn );
138 return eOut;
139}
EmcRecParameter::ECorr
double ECorr(int n) const
Definition EmcRecParameter.cxx:533
EmcRecParameter::GetInstance
static EmcRecParameter & GetInstance()
Definition EmcRecParameter.cxx:494

◆ ECorrection() [2/3]

RecEmcEnergy EmcRecShowerEnergy::ECorrection ( const RecEmcEnergy eIn)

◆ ECorrection() [3/3]

RecEmcEnergy EmcRecShowerEnergy::ECorrection ( const RecEmcEnergy eIn)

◆ ECorrTheta() [1/3]

RecEmcEnergy EmcRecShowerEnergy::ECorrTheta ( const RecEmcEnergy eIn,
const RecEmcID & id )

Definition at line 141 of file EmcRecShowerEnergy.cxx.

141 {
142 EmcRecParameter& Para = EmcRecParameter::GetInstance();
143 RecEmcEnergy eOut = eIn;
144
145 unsigned int npart = EmcID::barrel_ec( id );
146 unsigned int ntheta = EmcID::theta_module( id );
147
148 if ( npart == 1 ) { eOut *= 1.843 / Para.Peak( ntheta ); }
149 else if ( npart == 0 ) { eOut *= 1.843 / Para.Peak( ntheta + 44 ); }
150 else if ( npart == 2 ) { eOut *= 1.843 / Para.Peak( ntheta + 50 ); }
151
152 return eOut;
153}
EmcID::barrel_ec
static unsigned int barrel_ec(const Identifier &id)
Values of different levels (failure returns 0).
Definition EmcID.cxx:36
EmcID::theta_module
static unsigned int theta_module(const Identifier &id)
Definition EmcID.cxx:41
EmcRecParameter::Peak
double Peak(int n) const
Definition EmcRecParameter.cxx:547

◆ ECorrTheta() [2/3]

RecEmcEnergy EmcRecShowerEnergy::ECorrTheta ( const RecEmcEnergy eIn,
const RecEmcID & id )

◆ ECorrTheta() [3/3]

RecEmcEnergy EmcRecShowerEnergy::ECorrTheta ( const RecEmcEnergy eIn,
const RecEmcID & id )

◆ ELepton() [1/3]

void EmcRecShowerEnergy::ELepton ( RecEmcShower & aShower)

◆ ELepton() [2/3]

void EmcRecShowerEnergy::ELepton ( RecEmcShower & aShower)

◆ ELepton() [3/3]

void EmcRecShowerEnergy::ELepton ( RecEmcShower & aShower)

◆ Energy() [1/3]

void EmcRecShowerEnergy::Energy ( RecEmcShower & aShower)

Definition at line 6 of file EmcRecShowerEnergy.cxx.

6 {
7 RecEmcFractionMap::const_iterator cit;
8 RecEmcEnergy e1 = 0;
9 RecEmcEnergy e9 = 0;
10 RecEmcEnergy e25 = 0;
11 RecEmcEnergy elepton = 0;
12 RecEmcEnergy eall = 0;
13
14 EmcRecNeighbor nhb;
15
16 RecEmcID CellId = aShower.getShowerId();
17 int module = EmcID::barrel_ec( CellId );
18 RecEmcIDVector NearCell = nhb.GetNeighbors( CellId );
19 RecEmcIDVector NextNearCell = nhb.GetNextNeighbors( CellId );
20 RecEmcIDVector tmpNearCell;
21 RecEmcIDVector tmpNextNearCell;
22 i_RecEmcIDVector pNearCell;
23 i_RecEmcIDVector pNextNearCell;
24 vector<RecEmcEnergy> eVec;
25 vector<RecEmcEnergy>::const_iterator ciVec;
26
27 tmpNearCell.push_back( CellId );
28 tmpNextNearCell.push_back( CellId );
29
30 cit = aShower.Find( CellId );
31 // int time_seed = cit->second.getTime();
32 e1 = ( cit->second.getEnergy() ) * ( cit->second.getFraction() );
33 e9 += ( cit->second.getEnergy() ) * ( cit->second.getFraction() );
34 e25 += ( cit->second.getEnergy() ) * ( cit->second.getFraction() );
35
36 // e3x3
37 for ( pNearCell = NearCell.begin(); pNearCell != NearCell.end(); pNearCell++ )
38 {
39 cit = aShower.Find( *pNearCell );
40 if ( cit != aShower.End() )
41 {
42 tmpNearCell.push_back( *pNearCell );
43 tmpNextNearCell.push_back( *pNearCell );
44 e9 += cit->second.getEnergy() * cit->second.getFraction();
45 e25 += cit->second.getEnergy() * cit->second.getFraction();
46 }
47 }
48
49 // e5x5
50 for ( pNextNearCell = NextNearCell.begin(); pNextNearCell != NextNearCell.end();
51 pNextNearCell++ )
52 {
53 cit = aShower.Find( *pNextNearCell );
54 if ( cit != aShower.End() )
55 {
56 tmpNextNearCell.push_back( *pNextNearCell );
57 e25 += cit->second.getEnergy() * cit->second.getFraction();
58 }
59 }
60
61 // eall
62 for ( cit = aShower.Begin(); cit != aShower.End(); ++cit )
63 {
64 eall += ( cit->second.getEnergy() ) * ( cit->second.getFraction() );
65 eVec.push_back( cit->second.getEnergy() * cit->second.getFraction() );
66 }
67
68 // calculate number of hits from MC
69 int nHit, n;
70 EmcRecParameter& Para = EmcRecParameter::GetInstance();
71 nHit = (int)( Para.HitNb( 0 ) * log( Para.HitNb( 1 ) * e9 + Para.HitNb( 2 ) ) );
72 n = 0;
73
74 // sort by energy
75 sort( eVec.begin(), eVec.end(), greater<RecEmcEnergy>() );
76
77 for ( ciVec = eVec.begin(); ciVec != eVec.end(); ciVec++ )
78 {
79 if ( n < nHit )
80 {
81 elepton += *ciVec;
82 n++;
83 }
84 }
85
86 // energy correction
87 // RecEmcEnergy eCorr=ECorrTheta(e25,CellId);
88 // RecEmcEnergy eCorr=ECorrection(e25);
89 int getthetaid = EmcID::theta_module( CellId );
90 int getmodule = EmcID::barrel_ec( CellId );
91 if ( getthetaid > 21 ) getthetaid = 43 - getthetaid;
92 if ( getmodule == 1 ) getthetaid = getthetaid + 6;
93 double dthetaid = double( getthetaid );
94 double eCorr = Para.ECorrMC( e25, dthetaid );
95
96 // energy error
97 RecEmcEnergy de, de1, de2, de3;
98 de1 = Para.SigE( 0 ) / eCorr;
99 de2 = Para.SigE( 1 ) / pow( eCorr, 0.25 );
100 de3 = Para.SigE( 2 );
101 de = sqrt( de1 * de1 + de2 * de2 + de3 * de3 ) * perCent * eCorr;
102
103 double err = Para.ErrMC( e25, dthetaid );
104 if ( err > 0 ) de = err * e25;
105
106 aShower.setTrackId( -1 );
107 aShower.setCellId( CellId );
108 aShower.setModule( module );
109 aShower.setNumHits( aShower.getSize() );
110 aShower.setDE( de );
111 aShower.CellId3x3( tmpNearCell );
112 aShower.CellId5x5( tmpNextNearCell );
113 aShower.setESeed( e1 );
114 aShower.setE3x3( e9 );
115 aShower.setE5x5( e25 );
116 aShower.ELepton( elepton );
117 aShower.EAll( eall );
118 aShower.setEnergy( eCorr );
119
120 // cout<<"e1="<<aShower.eSeed()
121 // <<"\te9="<<aShower.e3x3()
122 // <<"\te25="<<aShower.e5x5()
123 // <<"\telepton="<<aShower.getELepton()
124 // <<"\teall="<<aShower.getEAll()
125 // <<"\tenergy="<<aShower.energy()<<endl;
126}
n
const Int_t n
Definition DataBase/tau_mode.c:57
e1
Double_t e1
Definition DataBase/tau_mode.c:7
EmcRecNeighbor::GetNeighbors
RecEmcIDVector GetNeighbors(const Identifier &id)
Definition EmcRecNeighbor.cxx:46
EmcRecNeighbor::GetNextNeighbors
RecEmcIDVector GetNextNeighbors(const Identifier &id)
Definition EmcRecNeighbor.cxx:378
EmcRecParameter::SigE
double SigE(int n) const
Definition EmcRecParameter.cxx:535
EmcRecParameter::ECorrMC
double ECorrMC(double eg, double theid) const
Definition EmcRecParameter.cxx:670
EmcRecParameter::HitNb
double HitNb(int n) const
Definition EmcRecParameter.cxx:541
EmcRecParameter::ErrMC
double ErrMC(double eg, double theid) const
Definition EmcRecParameter.cxx:686
RecEmcShower::End
RecEmcFractionMap::const_iterator End() const
Definition RecEmcShower.cxx:114
RecEmcShower::CellId3x3
void CellId3x3(RecEmcIDVector &id3x3)
Definition RecEmcShower.cxx:102
RecEmcShower::EAll
RecEmcEnergy EAll(RecEmcEnergy e)
Definition RecEmcShower.cxx:230
RecEmcShower::ELepton
RecEmcEnergy ELepton(RecEmcEnergy e)
Definition RecEmcShower.cxx:231
RecEmcShower::Begin
RecEmcFractionMap::const_iterator Begin() const
Definition RecEmcShower.cxx:112
RecEmcShower::Find
RecEmcFractionMap::const_iterator Find(const RecEmcID &CellId) const
Definition RecEmcShower.cxx:116
RecEmcShower::getSize
unsigned int getSize() const
Definition RecEmcShower.cxx:167
RecEmcShower::CellId5x5
void CellId5x5(RecEmcIDVector &id5x5)
Definition RecEmcShower.cxx:107

◆ Energy() [2/3]

void EmcRecShowerEnergy::Energy ( RecEmcShower & aShower)

◆ Energy() [3/3]

void EmcRecShowerEnergy::Energy ( RecEmcShower & aShower)
The documentation for this class was generated from the following files: