Dear root experts,

I want to make a TProfile with variable bin size. In the X axis I want the Energy of a random distribution and in the Y axis I want dN/dE (Flux). Here is my problem because I don’t know how to get this Flux. Here I copy an extract of my code:

[code]TProfile ToyMC(const double gamma, const int nrEvents ) {

const unsigned int nrBins = 9;

const double bins[nrBins+1] = {1.e+17.1,1.e+17.2,1.e+17.3,1.e+17.4,1.e+17.5,1.e+17.8,1.e+18.,1.e+18.2,1.e+18.5,1.e+20.};

TProfile *flux = new TProfile(“flux”,"",nrBins,bins);

TF1 * energyDistr = new TF1(“energyDistr”,“pow(x, [0])”,1.e+16.5,2.e+20);

energyDistr->SetParameter(0,gamma);

energyDistr->SetNpx(1000000);

for (Int_t i=0; i<nrEvents; i++){

const double energy = energyDistr->GetRandom();

if (measEn < 0)

continue;

Flux->Fill(energy,¿¿¿…???);

}

return *flux;

}

[/code]

Anyone could explain me how to make it?

Many thanks, Nahual