1 #ifndef __ALINEL_CRSEC_GEN_HH
2 #define __ALINEL_CRSEC_GEN_HH
21 G4double M_p = 938.2796;
24 G4double M_Al = M_p*A_Al;
27 G4double CTH = cos(th/2.0);
28 G4double STH = sin(th/2.0);
31 G4double ETA = 1.0+2.0*E_in*STH*STH/M_Al;
32 G4double E_out = E_in/ETA;
34 G4double Q2 = 4*E_in*E_out*STH*STH;
35 G4double myhbarc = hbarc / MeV / fermi;
36 G4double q2 = Q2/(myhbarc*myhbarc);
37 G4double qq = sqrt(q2);
42 G4double sigma_mott = pow(Z_Al*0.719982/E_in*CTH/(STH*STH),2)/(1+2*E_in/M_Al*STH*STH)*10000;
44 G4double sigma = sigma_mott * FF2;
57 ff2[0] =
new TF1(
"844",
"0.000427*exp(-0.5*((x-0.958792)/0.341168)**2)",0.5,1.8);
58 ff2[1] =
new TF1(
"1016",
"0.000706*exp(-0.5*((x-0.868077)/0.350506)**2)",0.5,1.8);
59 ff2[2] =
new TF1(
"2211",
"0.002337*exp(-0.5*((x-0.857955)/0.290595)**2)",0.5,1.8);
60 ff2[3] =
new TF1(
"2735",
"0.000306*exp(-0.5*((x-0.955732)/0.379688)**2)",0.5,1.8);
61 ff2[4] =
new TF1(
"2735",
"0.001598*exp(-0.5*((x-0.948107)/0.330035)**2)",0.5,1.8);
62 ff2[5] =
new TF1(
"3680",
"0.000042*exp(-0.5*((x-0.992771)/0.395955)**2)",0.5,1.8);
63 ff2[6] =
new TF1(
"4510",
"0.000135*exp(-0.5*((x-1.340052)/0.302506)**2)",0.5,1.8);
64 ff2[7] =
new TF1(
"4580",
"0.000162*exp(-0.5*((x-1.297271)/0.352787)**2)",0.5,1.8);
67 for(G4int i=0;i<8;i++){
68 if((q>=0.5) && (q<=1.8))
69 ffT += ff2[i]->Eval(q);
G4double getAlinel_FF2(G4double q)
G4double Alinel_crsec_gen(G4double E_in=1160, G4double th=8.0)