fixing implementation of jet method for hadronic correction so that clusters below...

diff --git a/PWGLF/totEt/AliAnalysisEtMonteCarlo.cxx b/PWGLF/totEt/AliAnalysisEtMonteCarlo.cxx
index 3825cbfd98b929ddbfdcc2c00cf283294e6c9a68..fd8c312ba5cedaecfa746a94fcc54c19d7f188cc 100644 (file)
--- a/PWGLF/totEt/AliAnalysisEtMonteCarlo.cxx
+++ b/PWGLF/totEt/AliAnalysisEtMonteCarlo.cxx
@@ -909,6 +909,7 @@ Int_t AliAnalysisEtMonteCarlo::AnalyseEvent(AliVEvent* ev,AliVEvent* ev2)
     Float_t subtotalSimulatedGammaEnergyAboveThreshold = 0.0;//signal of signal clusters passing cuts vs centrality
     Float_t subtotalReconstructedGammaEnergy = 0.0;//signal of signal clusters passing cuts vs centrality
     Float_t subtotalReconstructedSignalEnergy = 0.0;//signal of signal clusters passing cuts vs centrality
+    // cout<<"fsub "<<fsub<<endl;
 
     // loop the clusters
     for (int iCluster = 0; iCluster < nCluster; iCluster++ )
@@ -995,7 +996,8 @@ Int_t AliAnalysisEtMonteCarlo::AnalyseEvent(AliVEvent* ev,AliVEvent* ev2)
            fHistRCorrVsPtVsCent->Fill(rcorr,matchedTrackpt, fCentClass);
            //cout<<"rcorr "<<rcorr<<endl;
            Int_t n=caloCluster->GetNLabels() ;
-           if(fReconstructedE - fsub* track->P() > 0.0){//if more energy was deposited than the momentum of the track  and more than one particle led to the cluster
+           //if(fReconstructedE - fsub* track->P() > 0.0){//if more energy was deposited than the momentum of the track  and more than one particle led to the cluster
+           if(fSelector->PassMinEnergyCut( (fReconstructedE - fsub* track->P())*TMath::Sin(cp.Theta()) )){//if more energy was deposited than the momentum of the track  and more than one particle led to the cluster
              //then we say the cluster was not track matched but correct the energy
              nottrackmatched = kTRUE;
              //cout<<"Reassigning energy "<<fReconstructedEt;

diff --git a/PWGLF/totEt/AliAnalysisEtReconstructed.cxx b/PWGLF/totEt/AliAnalysisEtReconstructed.cxx
index 8d98afeeb432cb64c2596f2d0068e3c11b52b2bf..4c024f975e8bbd0ff628ea1b124e55e58e408154 100644 (file)
--- a/PWGLF/totEt/AliAnalysisEtReconstructed.cxx
+++ b/PWGLF/totEt/AliAnalysisEtReconstructed.cxx
@@ -312,7 +312,6 @@ Int_t AliAnalysisEtReconstructed::AnalyseEvent(AliVEvent* ev)
        fTotAllRawEtEffCorr +=GetCorrectionModification(*cluster,0,0,cent)* CorrectForReconstructionEfficiency(*cluster,cent);
 
        fClusterEnergyCent->Fill(GetCorrectionModification(*cluster,0,0,cent)*fReconstructedE,cent);
-
         Bool_t matched = kTRUE;//default to no track matched
        Int_t trackMatchedIndex = cluster->GetTrackMatchedIndex();//find the index of the matched track
        matched = !(fSelector->PassTrackMatchingCut(*cluster));//PassTrackMatchingCut is false if there is a matched track
@@ -325,7 +324,8 @@ Int_t AliAnalysisEtReconstructed::AnalyseEvent(AliVEvent* ev)
            if(matched){//if it is still matched see if the track p was less than the energy
              Float_t rcorr = TMath::Sin(cp.Theta())*GetCorrectionModification(*cluster,0,0,cent)*fReconstructedE;
              fHistRecoRCorrVsPtVsCent->Fill(rcorr,track->Pt(), fCentClass);
-             if(fReconstructedE - fsub* track->P() > 0.0){
+             if(fSelector->PassMinEnergyCut( (fReconstructedE - fsub* track->P())*TMath::Sin(cp.Theta()) )){//if more energy was deposited than the momentum of the track  and more than one particle led to the cluster
+               //            if(fReconstructedE - fsub* track->P() > 0.0){
                fReconstructedE = fReconstructedE - fsub* track->P();
                matched = kFALSE;
              }