Float_t clsPos[3] = {0,0,0};
c->GetPosition(clsPos);
TVector3 clsVec(clsPos);
- clsVec += fVecPv;
+ clsVec -= fVecPv;
Double_t Et = c->E()*TMath::Sin(clsVec.Theta());
if(fDebug)
printf("\tcluster eta=%1.1f,phi=%1.1f,E=%1.1f\n",clsVec.Eta(),clsVec.Phi(),c->E());
Float_t clsPos[3] = {0,0,0};
c->GetPosition(clsPos);
TVector3 cv(clsPos);
- cv += fVecPv;
+ cv -= fVecPv;
Double_t Et = c->E()*TMath::Sin(cv.Theta());
Float_t dphi = TMath::Abs(cv.Phi()-phicl);
Float_t deta = TMath::Abs(cv.Eta()-etacl);
Float_t clsPos[3] = {0,0,0};
c->GetPosition(clsPos);
TVector3 clsVec(clsPos);
- clsVec += fVecPv;
+ clsVec -= fVecPv;
Double_t Et = c->E()*TMath::Sin(clsVec.Theta());
//ESD
if(fStack){
Float_t clsPos[3] = {0,0,0};
c->GetPosition(clsPos);
TVector3 clsVec(clsPos);
- clsVec += fVecPv;
+ clsVec -= fVecPv;
Double_t cphi = clsVec.Phi();
Double_t ceta = clsVec.Eta();
Short_t id;