//_________________________________________________________________________
//
-// Class for study of EMCAL trigger behaviour
+// Class for study of EMCAL trigger behavior
//
// -- Author: Gustavo Conesa (CNRS-LPSC-Grenoble)
//////////////////////////////////////////////////////////////////////////////
ClassImp(AliAnaEMCALTriggerClusters)
-//____________________________
+//______________________________________________________
AliAnaEMCALTriggerClusters::AliAnaEMCALTriggerClusters() :
AliAnaCaloTrackCorrBaseClass(),
fRejectTrackMatch(0), fNCellsCut(0),
fMinM02(0), fMaxM02(0),
+fMomentum(),
// Histograms
fhE(0), fhESelected(0),
fhEtaPhi(0), fhEtaPhiSelected(0),
if(!badClusTrig)
{
- printf("AliAnaEMCALTriggerClusters::MakeAnalysisFillHistograms() - No cluster (bad-exotic trigger) found with requested index %d \n",idTrig);
+ AliWarning(Form("No cluster (bad-exotic trigger) found with requested index %d \n",idTrig));
return;
}
- TLorentzVector momBadClus;
+ badClusTrig->GetMomentum(fMomentum,GetVertex(0));
- badClusTrig->GetMomentum(momBadClus,GetVertex(0));
-
- Float_t etaclusterBad = momBadClus.Eta();
- Float_t phiclusterBad = momBadClus.Phi();
+ Float_t etaclusterBad = fMomentum.Eta();
+ Float_t phiclusterBad = fMomentum.Phi();
if( phiclusterBad < 0 ) phiclusterBad+=TMath::TwoPi();
Float_t tofclusterBad = badClusTrig->GetTOF()*1.e9;
Float_t eclusterBad = badClusTrig->E();
}
}
}
- else if(TMath::Abs(bc) >= 6)
- printf("AliAnaEMCALTriggerClusters::ClusterSelected() - Trigger BC not expected = %d\n",bc);
+ else if(TMath::Abs(bc) >= 6) AliWarning(Form("Trigger BC not expected = %d\n",bc));
}
const Int_t buffersize = 255;
char onePar[buffersize] ;
- snprintf(onePar,buffersize,"--- AliAnaEMCALTriggerClusters ---\n") ;
+ snprintf(onePar,buffersize,"--- AliAnaEMCALTriggerClusters ---:") ;
parList+=onePar ;
- snprintf(onePar,buffersize,"fRejectTrackMatch: %d\n",fRejectTrackMatch) ;
+ snprintf(onePar,buffersize,"fRejectTrackMatch: %d;",fRejectTrackMatch) ;
parList+=onePar ;
- snprintf(onePar,buffersize,"fMinM02: %2.2f, fMaxM02: %2.2f\n",fMinM02,fMaxM02) ;
+ snprintf(onePar,buffersize,"fMinM02: %2.2f, fMaxM02: %2.2f;",fMinM02,fMaxM02) ;
parList+=onePar ;
- snprintf(onePar,buffersize,"fNCellsCut: %d\n",fNCellsCut) ;
+ snprintf(onePar,buffersize,"fNCellsCut: %d;",fNCellsCut) ;
parList+=onePar ;
//Get parameters set in base class.
return new TObjString(parList) ;
}
-//________________________________________________________________________
+//___________________________________________________________
TList * AliAnaEMCALTriggerClusters::GetCreateOutputObjects()
{
// Create histograms to be saved in output file and
}
-//_______________________
+//_____________________________________
void AliAnaEMCALTriggerClusters::Init()
{
}
-//____________________________________________________________________________
+//_______________________________________________
void AliAnaEMCALTriggerClusters::InitParameters()
{
}
-//__________________________________________________________________
+//____________________________________________________________
void AliAnaEMCALTriggerClusters::MakeAnalysisFillHistograms()
{
//Do photon analysis and fill aods
if(!pl)
{
- Info("MakeAnalysisFillHistograms","TObjArray with clusters is NULL!\n");
+ AliWarning("TObjArray with clusters is NULL!");
return;
}
Int_t nCaloClusters = pl->GetEntriesFast();
Int_t idTrig = GetReader()->GetTriggerClusterIndex();
- TLorentzVector mom;
- if(GetDebug() > 0) printf("AliAnaEMCALTriggerClusters::MakeAnalysisFillHistograms() - Input cluster entries %d\n", nCaloClusters);
+ AliDebug(1,Form("Input cluster entries %d", nCaloClusters));
// Loop on clusters
for(Int_t icalo = 0; icalo < nCaloClusters; icalo++)
AliVCluster * calo = (AliVCluster*) (pl->At(icalo));
//printf("calo %d, %f\n",icalo,calo->E());
- calo->GetMomentum(mom,GetVertex(0)) ;
+ calo->GetMomentum(fMomentum,GetVertex(0)) ;
Float_t tofcluster = calo->GetTOF()*1.e9;
- Float_t ecluster = mom.E();
- Float_t etacluster = mom.Eta();
- Float_t phicluster = mom.Phi();
+ Float_t ecluster = fMomentum.E();
+ Float_t etacluster = fMomentum.Eta();
+ Float_t phicluster = fMomentum.Phi();
if(phicluster < 0) phicluster+=TMath::TwoPi();
FillRawClusterTriggerBCHistograms(calo->GetID(),ecluster,tofcluster,etacluster,phicluster);
//Check acceptance selection
if(IsFiducialCutOn())
{
- Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,"EMCAL") ;
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(fMomentum.Eta(),fMomentum.Phi(),kEMCAL) ;
if(! in ) continue ;
}
}
}
else if(TMath::Abs(bc) >= 6)
- printf("AliAnaEMCALTriggerClusters::MakeAnalysisFillHistograms() - Trigger BC not expected = %d\n",bc);
+ AliWarning(Form("Trigger BC not expected = %d",bc));
}// cluster loop
- if(GetDebug() > 1) printf("AliAnaEMCALTriggerClusters::MakeAnalysisFillHistograms() End fill histograms\n");
+ AliDebug(1,"End fill histograms");
}