|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
//
// CMSCGENnorm.cc P. Biallass 2006
//
// code based on l3cgen.f by T.Hebbeker
// implemented in CMSSW by P. Biallass 11.04.2006
// see header for documentation
//
#include "GeneratorInterface/CosmicMuonGenerator/interface/CMSCGENnorm.h"
//count muons which are around 100GeV and theta<33 deg (~vertical), do this for EVERY generated cosmic
int CMSCGENnorm::events_n100cos(double energy, double theta) {
if (energy > 99.5 && energy < 100.5) {
n100 = n100 + 1;
if (1. - cos(theta) < 1. / (2. * Pi)) { //theta is in rad
n100cos = n100cos + 1;
}
}
return n100cos;
}
// determine normalization using known flux, do this in the end. Percentage of rejected events and surface needs to be corrected for later!
// Note that the number of actually DICED cosmics is needed for normalisation, and sufficient statistics to have muons at 100 GeV!
float CMSCGENnorm::norm(int n100cos) {
flux = 2.63e-3; // +- 0.06e-3 [1/m**2/sr/GeV/s]
n = n100cos; // [1/sr/GeV]
//rate=N/runtime --> Nnorm ~ (1/runtime/m^2 at surface plane) as rate corresponds to known flux
Nnorm = flux / n;
//err of Nnorm = Nnorm* 1/sqrt(n)
return Nnorm;
}
|