89 Double32_t fHitEnergy = 0.;
90 Double32_t fCluEnergy = 0.;
100 for (Int_t iHit=0; iHit< (
fHitArray->GetEntriesFast()) ; iHit++)
114 TLorentzVector sum_clu(0., 0., 0., 0.), single_clu(0., 0., 0., 0.);
115 TVector3 mom_clu(0., 0., 0.);
118 for (Int_t iClu=0; iClu< (
fCluArray->GetEntriesFast()) ; iClu++)
121 fCluEnergy += cluster->
energy();
123 mom_clu.SetMagThetaPhi(cluster->
energy(), cluster->
theta(), cluster->
phi());
124 single_clu.SetVectM(mom_clu, 0.);
125 sum_clu = sum_clu + single_clu;
131 AddHit(fHitEnergy, nHitMult, nDigiMult, fCluEnergy, nCluMult, sum_clu.Vect());
virtual Double_t GetEnergy() const
a cluster (group of neighboring crystals) of hit emc crystals
represents the deposited energy of one emc crystal from simulation
int hit(Int_t nEvents=0, TString inFile="sim.root", TString parFile="par.root", TString inDigi="digi.root", TString outFile="hit.root", Int_t timeBased=0)
virtual Double_t energy() const
PndEmcHeader * AddHit(Double32_t ene_hit, Int_t mult_hit, Int_t mult_digi, Double32_t ene_clu, Int_t mult_clu, TVector3 p_clu)
TClonesArray * fDigiArray