caloCluster->GetPosition(pos) ;
fhEMCALPos->Fill(pos[0],pos[1],pos[2]) ;
fhEMCALEnergy->Fill(caloCluster->E()) ;
- fhEMCALDigits->Fill(entry, caloCluster->GetNumberOfDigits()) ;
+ fhEMCALDigits->Fill(entry, caloCluster->GetNCells()) ;
numberOfEmcalClustersv1++ ;
- numberOfDigitsInEmcal += caloCluster->GetNumberOfDigits() ;
+ numberOfDigitsInEmcal += caloCluster->GetNCells() ;
// Float_t * pid = clus->GetPid() ;
// if(pid[AliPID::kPhoton]>0.9){
emcalVector[numberOfPhotonsInEmcal] = new TVector3(pos[0],pos[1],pos[2]) ;
caloCluster->GetPosition( pos ) ;
fhPHOSEnergy->Fill( caloCluster->E() ) ;
fhPHOSPos->Fill( pos[0], pos[1], pos[2] ) ;
- fhPHOSDigits->Fill(entry, caloCluster->GetNumberOfDigits() ) ;
- numberOfDigitsInPhos += caloCluster->GetNumberOfDigits() ;
+ fhPHOSDigits->Fill(entry, caloCluster->GetNCells() ) ;
+ numberOfDigitsInPhos += caloCluster->GetNCells() ;
Double_t * pid = caloCluster->GetPid() ;
if(pid[AliPID::kPhoton] > 0.9) {
phosVector[numberOfPhotonsInPhos] = new TVector3(pos[0],pos[1],pos[2]) ;
fNCells(0),
fCellsAbsId(0x0),
fCellsAmpFraction(0x0),
- fDigitAmplitude(0x0),//not in use
- fDigitTime(0x0),//not in use
- fDigitIndex(0x0),//not in use
fEnergy(0),
fDispersion(0),
fChi2(0),
fNCells(clus.fNCells),
fCellsAbsId(),
fCellsAmpFraction(),
- fDigitAmplitude(clus.fDigitAmplitude?new TArrayS(*clus.fDigitAmplitude):0x0),//not in use
- fDigitTime(clus.fDigitTime?new TArrayS(*clus.fDigitTime):0x0),//not in use
- fDigitIndex(clus.fDigitIndex?new TArrayS(*clus.fDigitIndex):0x0),//not in use
fEnergy(clus.fEnergy),
fDispersion(clus.fDispersion),
fChi2(clus.fChi2),
fLabels = 0;
}
-
- if(source.fDigitAmplitude){
- // assign or copy construct
- if(fDigitAmplitude) *fDigitAmplitude = *source.fDigitAmplitude;
- else fDigitAmplitude = new TArrayS(*source.fDigitAmplitude);
- }
- else{
- delete fDigitAmplitude;
- fDigitAmplitude = 0;
- }
-
-
-
- if(source.fDigitTime){
- // assign or copy construct
- if(fDigitTime) *fDigitTime = *source.fDigitTime;
- else fDigitTime = new TArrayS(*source.fDigitTime);
- }
- else{
- delete fDigitTime;
- fDigitTime = 0;
- }
-
-
-
- if(source.fDigitIndex){
- // assign or copy construct
- if(fDigitIndex) *fDigitIndex = *source.fDigitIndex;
- else fDigitIndex = new TArrayS(*source.fDigitIndex);
- }
- else{
- delete fDigitIndex;
- fDigitIndex = 0;
- }
return *this;
}
+//_______________________________________________________________________
void AliESDCaloCluster::Copy(TObject &obj) const {
// this overwrites the virtual TOBject::Copy()
//
if(fTracksMatched)delete fTracksMatched;fTracksMatched = 0;
if(fLabels) delete fLabels; fLabels = 0;
- delete fDigitAmplitude; //not in use
- delete fDigitTime; //not in use
- delete fDigitIndex; //not in use
if(fCellsAmpFraction){ delete[] fCellsAmpFraction; fCellsAmpFraction=0;}
if(fCellsAbsId){ delete[] fCellsAbsId; fCellsAbsId = 0;}
}
//Vertex can be recovered with esd pointer doing:
//" Double_t vertex[3] ; esd->GetVertex()->GetXYZ(vertex) ; "
- if(vertex){//calculate direction from vertex
- fGlobalPos[0]-=vertex[0];
- fGlobalPos[1]-=vertex[1];
- fGlobalPos[2]-=vertex[2];
- }
-
- Double_t r = TMath::Sqrt(fGlobalPos[0]*fGlobalPos[0]+
- fGlobalPos[1]*fGlobalPos[1]+
- fGlobalPos[2]*fGlobalPos[2] ) ;
-
- p.SetPxPyPzE( fEnergy*fGlobalPos[0]/r, fEnergy*fGlobalPos[1]/r, fEnergy*fGlobalPos[2]/r, fEnergy) ;
-
+ Double32_t pos[3]={ fGlobalPos[0], fGlobalPos[1], fGlobalPos[2]};
+ if(vertex){//calculate direction from vertex
+ pos[0]-=vertex[0];
+ pos[1]-=vertex[1];
+ pos[2]-=vertex[2];
+ }
+
+ Double_t r = TMath::Sqrt(pos[0]*pos[0]+pos[1]*pos[1]+pos[2]*pos[2] ) ;
+
+ p.SetPxPyPzE( fEnergy*pos[0]/r, fEnergy*pos[1]/r, fEnergy*pos[2]/r, fEnergy) ;
}
//_______________________________________________________________________
if (fCellsAmpFraction && i >=0 && i < fNCells ) return fCellsAmpFraction[i];
else return -1;}
- //_____________________________________________________
- //Not used anymore, kept to avoid backward incompatibility
- void AddDigitIndex(TArrayS & array) { fDigitIndex = new TArrayS(array) ; Warning("AddDigitAmplitude","This method is no more in use") ;}
- void AddDigitAmplitude(TArrayS & array) { fDigitAmplitude = new TArrayS(array) ; Warning("AddDigitAmplitude","This method is no more in use") ;}
- void AddDigitTime(TArrayS & array) { fDigitTime = new TArrayS(array) ;Warning("AddDigitTime","This method is no more in use") ;}
- TArrayS * GetDigitAmplitude() const {return fDigitAmplitude;}
- TArrayS * GetDigitTime() const {return fDigitTime;}
- TArrayS * GetDigitIndex() const {return fDigitIndex;}
- Int_t GetNumberOfDigits() const { return -1;}
- //_____________________________________________________
-
protected:
TArrayI * fTracksMatched; //Index of tracks close to cluster. First entry is the most likely match.
UShort_t *fCellsAbsId; //[fNCells] array of cell absId numbers
Double32_t *fCellsAmpFraction; //[fNCells][0.,1.,16] array with cell amplitudes fraction.
- //__________________________________________________________
- //Not in use
- TArrayS * fDigitAmplitude; //digit energy (integer units)
- TArrayS * fDigitTime; //time of this digit (integer units)
- TArrayS * fDigitIndex; //calorimeter digit index
- //_________________________________________________________
Double32_t fGlobalPos[3]; // position in global coordinate systemD
Double32_t fEnergy; // energy measured by calorimeter
Char_t fClusterType; // Flag for different cluster type/versions
Double32_t fTOF; //[0,0,12] time-of-flight
- ClassDef(AliESDCaloCluster,9) //ESDCaloCluster
+ ClassDef(AliESDCaloCluster,10) //ESDCaloCluster
};
#endif
fESDOld(0),
fESDFriendOld(0),
fConnected(kFALSE),
- fUseOwnList(kFALSE),
- fEMCALClusters(0),
- fFirstEMCALCluster(-1),
- fPHOSClusters(0),
- fFirstPHOSCluster(-1)
+ fUseOwnList(kFALSE)
{
}
//______________________________________________________________________________
fESDOld(esd.fESDOld ? new AliESD(*esd.fESDOld) : 0),
fESDFriendOld(esd.fESDFriendOld ? new AliESDfriend(*esd.fESDFriendOld) : 0),
fConnected(esd.fConnected),
- fUseOwnList(esd.fUseOwnList),
- fEMCALClusters(esd.fEMCALClusters),
- fFirstEMCALCluster(esd.fFirstEMCALCluster),
- fPHOSClusters(esd.fPHOSClusters),
- fFirstPHOSCluster(esd.fFirstPHOSCluster)
+ fUseOwnList(esd.fUseOwnList)
{
// CKB init in the constructor list and only add here ...
fConnected = source.fConnected;
fUseOwnList = source.fUseOwnList;
- fEMCALClusters = source.fEMCALClusters;
- fFirstEMCALCluster = source.fFirstEMCALCluster;
- fPHOSClusters = source.fPHOSClusters;
- fFirstPHOSCluster = source.fFirstPHOSCluster;
-
return *this;
// don't reset fconnected fConnected and the list
- fEMCALClusters=0;
- fFirstEMCALCluster=-1;
- fPHOSClusters=0;
- fFirstPHOSCluster=-1;
}
if(fEMCALCells)printf(" EMCALCells %d\n", fEMCALCells->GetNumberOfCells());
else printf(" EMCALCells not in the Event\n");
printf(" CaloClusters %d\n", GetNumberOfCaloClusters());
- printf(" phos %d\n", GetNumberOfPHOSClusters());
- printf(" emcal %d\n", GetNumberOfEMCALClusters());
printf(" FMD %s\n", (fESDFMD ? "yes" : "no"));
printf(" VZERO %s\n", (fESDVZERO ? "yes" : "no"));
TObject* pHLTDecision=GetHLTTriggerDecision();
void SetUseOwnList(Bool_t b){fUseOwnList = b;}
Bool_t GetUseOwnList() const {return fUseOwnList;}
-
- // Remove this stuff CKB?
- //---------------------------------------------------
- Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
- void SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
- Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
- void SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
-
- Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
- void SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
- void SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; }
- Int_t GetFirstPHOSCluster() const { return fFirstPHOSCluster ; }
- //-------------------------------------------------------
TArrayF *GetEMCALTriggerPosition() const {return fEMCALTrigger?fEMCALTrigger->GetTriggerPosition():0x0;}
TArrayF *GetEMCALTriggerAmplitudes() const {return fEMCALTrigger?fEMCALTrigger->GetTriggerAmplitudes():0x0;}
static const char* fgkESDListName[kESDListN]; //!
- // Remove this stuff CKB
- Int_t fEMCALClusters; // Number of EMCAL clusters (subset of caloclusters)
- Int_t fFirstEMCALCluster; // First EMCAL cluster in the fCaloClusters list
-
- Int_t fPHOSClusters; // Number of PHOS clusters (subset of caloclusters)
- Int_t fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list
-
- ClassDef(AliESDEvent,10) //ESDEvent class
+ ClassDef(AliESDEvent,11) //ESDEvent class
};
#endif