message: add histogram for shower shape
diff: @@ -429,8 +429,8 @@ void AliAnalysisTaskPHOSPi0EtaToGammaGamma::UserCreateOutputObjects()
  for(Int_t imod=1;imod<Nmod;imod++) fOutputContainer->Add(new TH2F(Form("hClusterEvsTM%d",imod),Form("Cluster E vs TOF M%d;E (GeV);TOF (ns)",imod) ,500,0,50, 1000,-500,500));
  fOutputContainer->Add(new TH2F("hClusterEvsN","Cluster E vs N_{cell};E (GeV);N_{cell}",500,0,50,100,0.5,100.5));
  - fOutputContainer->Add(new TH2F("hClusterEvsM02","Cluster E vs M02;E (GeV);M02 (cm)" ,500,0,50,100,0,10));
  - fOutputContainer->Add(new TH2F("hClusterNvsM02","Cluster N vs M02;N_{cell};M02 (cm)",100,0.5,100.5,100,0,10));
  + fOutputContainer->Add(new TH2F("hClusterEvsM02","Cluster E vs M02;E (GeV);M02 (cm)",500,0,50,50,0,5));
  + fOutputContainer->Add(new TH2F("hClusterEvsM20","Cluster E vs M20;E (GeV);M20 (cm)",500,0,50,50,0,5));
  fOutputContainer->Add(new TH2F("hFullDispvsFullE","full dispersion vs full E;E (GeV);dispersion (#sigma)",100,0,50,100,0,10));
  fOutputContainer->Add(new TH2F("hCoreDispvsCoreE","core dispersion vs core E;E (GeV);dispersion (#sigma)",100,0,50,100,0,10));
  fOutputContainer->Add(new TH2F("hFullDispvsCoreE","full dispersion vs full E;E (GeV);dispersion (#sigma)",100,0,50,100,0,10));
  @@ -1646,7 +1646,7 @@ void AliAnalysisTaskPHOSPi0EtaToGammaGamma::ClusterQA()
  Double_t position[3] = {};
  Int_t digMult=0;
  Double_t energy=0,tof=0,eta=0,phi=0;
  - Double_t M02=0;
  + Double_t M02=0, M20=0;
  Double_t R = 0, coreR=0;
  Double_t coreE = 0;
  @@ -1668,6 +1668,7 @@ void AliAnalysisTaskPHOSPi0EtaToGammaGamma::ClusterQA()
  digMult = ph->GetNCells();
  tof = ph->GetTime();//unit is second.
  + M20 = ph->GetLambda1();
  M02 = ph->GetLambda2();
  position[0] = ph->EMCx();
  @@ -1707,7 +1708,7 @@ void AliAnalysisTaskPHOSPi0EtaToGammaGamma::ClusterQA()
  FillHistogramTH2(fOutputContainer,"hClusterEvsN",energy,digMult);
  FillHistogramTH2(fOutputContainer,"hClusterEvsM02",energy,M02);
  - FillHistogramTH2(fOutputContainer,"hClusterNvsM02",digMult,M02);
  + FillHistogramTH2(fOutputContainer,"hClusterEvsM20",energy,M20);
  FillHistogramTH2(fOutputContainer,"hFullDispvsFullE",energy,R);
  FillHistogramTH2(fOutputContainer,"hCoreDispvsCoreE",coreE,coreR);
  FillHistogramTH2(fOutputContainer,"hFullDispvsCoreE",coreE,R);
  