AliAnaPi0EbE::AliAnaPi0EbE() :
AliAnaCaloTrackCorrBaseClass(),fAnaType(kIMCalo), fCalorimeter(""),
fMinDist(0.),fMinDist2(0.), fMinDist3(0.),
- fFillWeightHistograms(kFALSE), fFillTMHisto(0),
+ fFillWeightHistograms(kFALSE), fFillTMHisto(0), fFillSelectClHisto(0),
fInputAODGammaConvName(""),
// Histograms
fhPt(0), fhE(0),
////////
- if(fAnaType == kIMCalo || fAnaType == kIMCaloTracks )
+ if( fFillSelectClHisto && (fAnaType == kIMCalo || fAnaType == kIMCaloTracks) )
{
fhEDispersion = new TH2F
outputContainer->Add(fhELambda1LocMaxN) ;
- if(fCalorimeter=="EMCAL"){
+ if(fCalorimeter=="EMCAL")
+ {
fhELambda0NoTRD = new TH2F
("hELambda0NoTRD","Selected #pi^{0} (#eta) pairs: E vs #lambda_{0}, not behind TRD",nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
fhELambda0NoTRD->SetYTitle("#lambda_{0}^{2}");
}// Invariant mass analysis in calorimeters and calorimeter + conversion photons
- if(fAnaType == kIMCalo){
+ if(fAnaType == kIMCalo)
+ {
fhEPairDiffTime = new TH2F("hEPairDiffTime","cluster pair time difference vs E",nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
fhEPairDiffTime->SetXTitle("E_{pair} (GeV)");
fhEPairDiffTime->SetYTitle("#Delta t (ns)");
printf("AliAnaPi0EbE::MakeInvMassInCalorimeter() - Selected gamma pair: pt %f, phi %f, eta%f \n",(mom1+mom2).Pt(), (mom1+mom2).Phi()*180./3.1416, (mom1+mom2).Eta());
//Fill some histograms about shower shape
- if(clusters && GetReader()->GetDataType()!=AliCaloTrackReader::kMC)
+ if(fFillSelectClHisto && clusters && GetReader()->GetDataType()!=AliCaloTrackReader::kMC)
{
FillSelectedClusterHistograms(cluster1, nMaxima1, photon1->GetTag());
FillSelectedClusterHistograms(cluster2, nMaxima2, photon2->GetTag());
if(GetDebug() > 1) printf("AliAnaPi0EbE::MakeInvMassInCalorimeterAndCTS() - Selected gamma pair: pt %f, phi %f, eta%f\n",(mom1+mom2).Pt(), (mom1+mom2).Phi()*180./3.1416, (mom1+mom2).Eta());
//Fill some histograms about shower shape
- if(cluster && GetReader()->GetDataType()!=AliCaloTrackReader::kMC)
+ if(fFillSelectClHisto && cluster && GetReader()->GetDataType()!=AliCaloTrackReader::kMC)
{
FillSelectedClusterHistograms(cluster, nMaxima, photon1->GetTag());
}
void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; }
void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; }
+ void SwitchOnSelectedClusterHistoFill() { fFillSelectClHisto = kTRUE ; }
+ void SwitchOffSelectedClusterHistoFill() { fFillSelectClHisto = kFALSE ; }
+
//For histograms
enum mcTypes { kmcPhoton = 0, kmcConversion = 1, kmcPi0 = 2,
kmcEta = 3, kmcElectron = 4, kmcHadron = 5 };
Bool_t fFillWeightHistograms ; // Fill weigth histograms
Bool_t fFillTMHisto; // Fill track matching plots
+ Bool_t fFillSelectClHisto; // Fill selected cluster histograms
//Only for combination of calorimeter and conversion photons, kIMCaloTracks
TString fInputAODGammaConvName; // Name of AOD branch with conversion photons
AliAnaPi0EbE( const AliAnaPi0EbE & g) ; // cpy ctor
AliAnaPi0EbE & operator = (const AliAnaPi0EbE & g) ; // cpy assignment
- ClassDef(AliAnaPi0EbE,13)
+ ClassDef(AliAnaPi0EbE,14)
} ;