// Cluster control histograms
fhTrackMatchedDEta(), fhTrackMatchedDPhi(), fhTrackMatchedDEtaDPhi(),
fhdEdx(), fhEOverP(), fhTrackMatchedMCParticle(),
-fhELambda0() , fhPtLambda0() ,
-fhELambda1(), fhELambda0SSBkg(),
+fhELambda0() , fhPtLambda0() , fhELambda1(),
fhELambda0TRD(), fhPtLambda0TRD(), fhELambda1TRD(),
fhELambda0MCPhoton(), fhPtLambda0MCPhotonPrompt(), fhPtLambda0MCPhotonFrag(),
fhELambda0MCPi0(), fhELambda0MCPi0Decay(),
}
//_____________________________________________________________________________________________________________________
-void AliAnaParticleIsolation::FillTrackMatchingShowerShapeControlHistograms(AliAODPWG4ParticleCorrelation *pCandidate,
- AliCaloPID * pid)
+void AliAnaParticleIsolation::FillTrackMatchingShowerShapeControlHistograms(AliAODPWG4ParticleCorrelation *pCandidate)
{
// Fill Track matching and Shower Shape control histograms
if(!fFillTMHisto && !fFillSSHisto) return;
else if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCOtherDecay)) fhELambda0MCOtherDecay[isolated]->Fill(energy, cluster->GetM02());
// else if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCConversion)) fhPtNoIsoConversion ->Fill(energy, cluster->GetM02());
- else if(!GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCElectron)) fhELambda0MCHadron [isolated]->Fill(energy, cluster->GetM02());
-
+ else if(!GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCElectron)) fhELambda0MCHadron [isolated]->Fill(energy, cluster->GetM02());
}
if(fCalorimeter == "EMCAL" && GetModuleNumber(cluster) > 5) // TO DO: CHANGE FOR 2012
if (nMaxima==1) { fhELambda0LocMax1[isolated]->Fill(energy,cluster->GetM02()); fhELambda1LocMax1[isolated]->Fill(energy,cluster->GetM20()); }
else if(nMaxima==2) { fhELambda0LocMax2[isolated]->Fill(energy,cluster->GetM02()); fhELambda1LocMax2[isolated]->Fill(energy,cluster->GetM20()); }
else { fhELambda0LocMaxN[isolated]->Fill(energy,cluster->GetM02()); fhELambda1LocMaxN[isolated]->Fill(energy,cluster->GetM20()); }
-
- if(isolated==0)
- {
- //Analyse non-isolated events
- Int_t n = 0;
- Int_t nfrac = 0;
- Bool_t iso = kFALSE ;
- Float_t coneptsum = 0 ;
-
- TObjArray * plNe = pCandidate->GetObjArray(GetAODObjArrayName()+"Clusters");
- TObjArray * plCTS = pCandidate->GetObjArray(GetAODObjArrayName()+"Tracks");
-
- GetIsolationCut()->SetPtThresholdMax(1.);
- GetIsolationCut()->MakeIsolationCut(plCTS, plNe,
- GetReader(), pid,
- kFALSE, pCandidate, "",
- n,nfrac,coneptsum, iso);
-
- if (!iso) fhELambda0SSBkg->Fill(energy, cluster->GetM02());
-
- if(GetDebug() > 0) printf("AliAnaParticleIsolation::MakeAnalysisFillHistograms() - Energy Sum in Isolation Cone %2.2f\n", coneptsum);
- }
-
- GetIsolationCut()->SetPtThresholdMax(10000.);
-
} // SS histo fill
outputContainer->Add(fhConeSumPtPhiUESubTrackCellTrigEtaPhi) ;
}
-
- if(fFillSSHisto)
- {
- fhELambda0SSBkg = new TH2F("hELambda0SSBkg",
- Form("Non isolated clusters : E vs #lambda_{0}, %s",parTitle.Data()),
- nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
- fhELambda0SSBkg->SetYTitle("#lambda_{0}^{2}");
- fhELambda0SSBkg->SetXTitle("#it{E} (GeV)");
- outputContainer->Add(fhELambda0SSBkg) ;
- }
-
+
for(Int_t iso = 0; iso < 2; iso++)
{
if(fFillTMHisto)
void FillAcceptanceHistograms();
- void FillTrackMatchingShowerShapeControlHistograms(AliAODPWG4ParticleCorrelation * pCandidate,
- AliCaloPID * pid) ;
+ void FillTrackMatchingShowerShapeControlHistograms(AliAODPWG4ParticleCorrelation * pCandidate) ;
void MakeSeveralICAnalysis( AliAODPWG4ParticleCorrelation * ph ) ;
TH2F * fhELambda0[2]; //! Shower shape of (non) isolated photons (do not apply SS cut previously)
TH2F * fhPtLambda0[2]; //! Shower shape of (non) isolated photons (do not apply SS cut previously)
TH2F * fhELambda1[2]; //! Shower shape of (non) isolated photons (do not apply SS cut previously)
- TH2F * fhELambda0SSBkg; //! Shower shape of non isolated photons close to isolation threshold (do not apply SS cut previously)
TH2F * fhELambda0TRD[2]; //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
TH2F * fhPtLambda0TRD[2]; //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
TH2F * fhELambda1TRD[2]; //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)