fHistList->Add(fh2CentvsRho);
fHistList->Add(fh2CentvsSigma);
- fh2MultvsRho = new TH2F("fh2MultvsRho","mult vs background density",4000,0.,4000.,2000,0.,200.);
- fh2MultvsSigma = new TH2F("fh2MultvsSigma","mult vs backgroun sigma",4000,0.,4000.,1000,0.,50.);
+ fh2MultvsRho = new TH2F("fh2MultvsRho","mult vs background density",1000,0.,5000.,300,0.,300.);
+ fh2MultvsSigma = new TH2F("fh2MultvsSigma","mult vs backgroun sigma",1000,0.,5000.,500,0.,50.);
fHistList->Add(fh2MultvsRho);
fHistList->Add(fh2MultvsSigma);
Double_t phiC = rC->Phi();
// massless jet, unit vector
pTC = rC->ChargedBgEnergy();
- if(pTC<=0)pTC = 0.1; // for almost empty events
+ if(pTC<=0)pTC = 0.001; // for almost empty events
Double_t thetaC = 2.*TMath::ATan(TMath::Exp(-etaC));
Double_t pZC = pTC/TMath::Tan(thetaC);
Double_t pXC = pTC * TMath::Cos(phiC);
Double_t phiC = rC->Phi();
// massless jet, unit vector
pTC = rC->ChargedBgEnergy();
- if(pTC<=0)pTC = 0.1;// for almost empty events
+ if(pTC<=0)pTC = 0.001;// for almost empty events
Double_t thetaC = 2.*TMath::ATan(TMath::Exp(-etaC));
Double_t pZC = pTC/TMath::Tan(thetaC);
Double_t pXC = pTC * TMath::Cos(phiC);