commented logging message
[u/mrichter/AliRoot.git] / TOF / CreateCalibPars_Miscalibrated_Scaled.C
CommitLineData
12ec3441 1void CreateCalibPars_Miscalibrated_Scaled(){
2 // Create TOF Calibration Object for miscalibrated detector
3 // and write it on CDB
4 AliTOFGeometry *geom = new AliTOFGeometryV5();
5 AliTOFcalib *tofcalib = new AliTOFcalib(geom);
6 AliTOFCal *tofCal= new AliTOFCal(geom);
7 tofCal->CreateArray();
8
9 // Input data for decalibration
10
11 TFile f("$ALICE_ROOT/TOF/data/spectrumScaled.root","READ");
12
13 TH1F *hTimeToTFit= (TH1F*)f.Get("hTimeToTScaled");
14 TF1 *fit=hTimeToTFit->GetFunction("pol5");
15
16 // Slewing parameters (same for all channels)
17
18 Float_t par[6] = {0.,0.,0.,0.,0.,0.};
19 for(Int_t i =0;i<6;i++){
20 par[i]=fit->GetParameter(i);
0adec857 21 cout << " Slewing parameter " <<i<<" =" << par[i] << endl;
12ec3441 22 }
23
0adec857 24 // Global time offset (randomly gen, gaussian with mean = 0.3, sig=0.08 ns)
25
26 Float_t delay=0.;
27 Float_t meanDelay=0.3;
28 Float_t sigmaDelay=0.08;
12ec3441 29
30 // ToT spectrum
0adec857 31
12ec3441 32 TH1F *hToT= (TH1F*)f.Get("hToTScaled");
33
34 // Fill the Sim calibration object
35
0adec857 36 TRandom *rnd = new TRandom(4357);
12ec3441 37 for (Int_t ipad = 0 ; ipad<tofCal->NPads(); ipad++){
38 AliTOFChannel *calChannel = tofCal->GetChannel(ipad);
39 calChannel->SetSlewPar(par);
0adec857 40 delay=rnd->Gaus(meanDelay,sigmaDelay);
12ec3441 41 calChannel->SetDelay(delay);
42 }
43 tofcalib->WriteSimParOnCDB("TOF/Calib",0,0,tofCal,hToT);
44 f.Close();
45}
46
47