Int_t inModule = -1;
Int_t iColLoc = -1;
if(iCol < AliEMCALGeoParams::fgkEMCALCols){ // if the SM number is odd the column is the one corresponding in the supermodule
- iModule = 2*iSector + 1;
+ inModule = 2*iSector + 1;
iColLoc = iCol;
}
else if(iCol > AliEMCALGeoParams::fgkEMCALCols - 1){ // if the SM number is even the column isn't the one corresponding in the supermodule
if(TMath::Abs(iCol-colCellLeadingClust)<nbConeSize && TMath::Abs(iCol+colCellLeadingClust)>nbConeSize){
if(iRow<iRowMaxCone && iRow>iRowMinCone){
- Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(iModule,iRow,iCol); // verifier les iRow et iCol
+ Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(inModule,iRowLoc,iColLoc); // verifier les iRow et iCol
sumEnergyPhiBandCells+=cells->GetAmplitude(iabsId); //should be Et
}
}
else if (TMath::Abs(iCol-colCellLeadingClust)>nbConeSize && TMath::Abs(iCol+colCellLeadingClust)<nbConeSize){
- Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(iModule,iRowLoc,iColLoc); // verifier les iRow et iCol
+ Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(inModule,iRowLoc,iColLoc); // verifier les iRow et iCol
sumEnergyConeCells+=cells->GetAmplitude(iabsId); //should be Et
}
}
Int_t inModule = -1;
Int_t iColLoc = -1;
if(iCol < AliEMCALGeoParams::fgkEMCALCols){ // if the SM number is odd the column is the one corresponding in the supermodule
- iModule = 2*iSector + 1;
+ inModule = 2*iSector + 1;
iColLoc = iCol;
}
else if(iCol > AliEMCALGeoParams::fgkEMCALCols - 1){ // if the SM number is even the column isn't the one corresponding in the supermodule
if(TMath::Abs(iCol-colCellLeadingClust)<nbConeSize && TMath::Abs(iCol+colCellLeadingClust)>nbConeSize){
if(iCol<iColMaxCone && iCol>iColMinCone){
- Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(iModule,iRow,iCol); // verifier les iRow et iCol
+ Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(inModule,iRowLoc,iColLoc); // verifier les iRow et iCol
sumEnergyEtaBandCells+=cells->GetAmplitude(iabsId); //should be Et
}
}
else if (TMath::Abs(iCol-colCellLeadingClust)>nbConeSize && TMath::Abs(iCol+colCellLeadingClust)<nbConeSize){
- Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(iModule,iRowLoc,iColLoc); // verifier les iRow et iCol
+ Int_t iabsId = emcalGeom->GetAbsCellIdFromCellIndexes(inModule,iRowLoc,iColLoc); // verifier les iRow et iCol
sumEnergyConeCells+=cells->GetAmplitude(iabsId); //should be Et
}
}
nTracks++;
}
- Double_t eTCOI = 0., m02COI = 0., lambda0cluster = 0., phiCOI = 0., etaCOI = 0., ptmc = 0., mcptsum = 0.;
+ Double_t eTCOI = 0., m02COI = 0., lambda0cluster = 0., ptmc = 0., mcptsum = 0.;
//Int_t Ntracks;
//Definition of the Array for Davide's Output
const Int_t ndims = fNDimensions;
eTCOI = vecCOI.Et();
m02COI = coi->GetM02();
- etaCOI = vecCOI.Eta();
- phiCOI = vecCOI.Phi();
+
// ******** Isolation and UE calculation with different methods *********