- const Int_t kNpasa = 200; // kNpasa should be equal to fTRFbin!
-
- Float_t signal[kNpasa]={ 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000
- , 0.0002, 0.0007, 0.0026, 0.0089, 0.0253, 0.0612, 0.1319
- , 0.2416, 0.3913, 0.5609, 0.7295, 0.8662, 0.9581, 1.0000
- , 0.9990, 0.9611, 0.8995, 0.8269, 0.7495, 0.6714, 0.5987
- , 0.5334, 0.4756, 0.4249, 0.3811, 0.3433, 0.3110, 0.2837
- , 0.2607, 0.2409, 0.2243, 0.2099, 0.1974, 0.1868, 0.1776
- , 0.1695, 0.1627, 0.1566, 0.1509, 0.1457, 0.1407, 0.1362
- , 0.1317, 0.1274, 0.1233, 0.1196, 0.1162, 0.1131, 0.1102
- , 0.1075, 0.1051, 0.1026, 0.1004, 0.0979, 0.0956, 0.0934
- , 0.0912, 0.0892, 0.0875, 0.0858, 0.0843, 0.0829, 0.0815
- , 0.0799, 0.0786, 0.0772, 0.0757, 0.0741, 0.0729, 0.0718
- , 0.0706, 0.0692, 0.0680, 0.0669, 0.0655, 0.0643, 0.0630
- , 0.0618, 0.0607, 0.0596, 0.0587, 0.0576, 0.0568, 0.0558
- , 0.0550, 0.0541, 0.0531, 0.0522, 0.0513, 0.0505, 0.0497
- , 0.0490, 0.0484, 0.0474, 0.0465, 0.0457, 0.0449, 0.0441
- , 0.0433, 0.0425, 0.0417, 0.0410, 0.0402, 0.0395, 0.0388
- , 0.0381, 0.0374, 0.0368, 0.0361, 0.0354, 0.0348, 0.0342
- , 0.0336, 0.0330, 0.0324, 0.0318, 0.0312, 0.0306, 0.0301
- , 0.0296, 0.0290, 0.0285, 0.0280, 0.0275, 0.0270, 0.0265
- , 0.0260, 0.0256, 0.0251, 0.0246, 0.0242, 0.0238, 0.0233
- , 0.0229, 0.0225, 0.0221, 0.0217, 0.0213, 0.0209, 0.0206
- , 0.0202, 0.0198, 0.0195, 0.0191, 0.0188, 0.0184, 0.0181
- , 0.0178, 0.0175, 0.0171, 0.0168, 0.0165, 0.0162, 0.0159
- , 0.0157, 0.0154, 0.0151, 0.0148, 0.0146, 0.0143, 0.0140
- , 0.0138, 0.0135, 0.0133, 0.0131, 0.0128, 0.0126, 0.0124
- , 0.0121, 0.0119, 0.0120, 0.0115, 0.0113, 0.0111, 0.0109
- , 0.0107, 0.0105, 0.0103, 0.0101, 0.0100, 0.0098, 0.0096
- , 0.0094, 0.0092, 0.0091, 0.0089, 0.0088, 0.0086, 0.0084
- , 0.0083, 0.0081, 0.0080, 0.0078 };
-
-// // Andronic & Bercuci parametrization
-// // define new TRF parametrization
-// // normalizing constant to Fe signal
-// const Float_t k1=1.055;
-// // time constants
-// const Float_t t1=0.04;
-// const Float_t t2=.9;
-// // the relative fraction of the long component
-// const Float_t k2=.15;
-// // time offset for Fe
-// const Float_t t0=-.29;
-// Float_t x = t0; Int_t index;
-// for(int i=0; i<kNpasa; i++) signal[i] = 0.;
-// for(int i=0; i<kNpasa; i++){
-// index = i+6;
-// if(index >= kNpasa) break;
-// x += .02;
-// signal[index]=k1*(TMath::Power((x-t0)/t1, 2.5)*(exp(-(x-t0)/t1))+k2*exp(-(x-t0)/t2));
-// }
-
- Float_t xtalk[kNpasa];