How to align simulation with experiment?

[weller]

@BenjaminW

the above plot shows the average deposited energy per particle, right?

if the plot was made with the code from your other thread that i linked to

@BenjaminW

    G4double energy = preStepPoint->GetTotalEnergy();

[...]

	// Total energy
	man->FillNtupleDColumn(0, 0, energy);
	man->AddNtupleRow(0);

	// Distinguishing detectors
	if(history_depth < 5)	
	{
		if(touchable->GetReplicaNumber(1) == 200){man->FillNtupleDColumn(1, 0, energy); man->AddNtupleRow(1);}
		if(touchable->GetReplicaNumber(1) == 210){man->FillNtupleDColumn(2, 0, energy); man->AddNtupleRow(2);}
		if(touchable->GetReplicaNumber(1) == 220){man->FillNtupleDColumn(3, 0, energy); man->AddNtupleRow(3);}
	}

, then no. it then shows the total energy of the incident particles, and with the first step into the SD, the track is killed. The plot shows the spectrum of your source (which also explains why the subplots are all so similar: it is just a histogram of your source energy distribution).