corrected TRD/TOF MV position
[u/mrichter/AliRoot.git] / TOF / CreateCalibPars_Miscalibrated.C
... / ...
CommitLineData
1void CreateCalibPars_Miscalibrated(){
2 // Create TOF Calibration Object for miscalibrated detector
3 // and write it on CDB
4 AliTOFcalib *tofcalib = new AliTOFcalib();
5 tofcalib->CreateCalArrays();
6 TObjArray *tofCalOffline = (TObjArray*) tofcalib->GetTOFCalArrayOffline();
7 // Input data for decalibration
8
9 TFile f("$ALICE_ROOT/TOF/data/spectrum.root","READ");
10
11 TH1F *hTimeToTFit= (TH1F*)f.Get("hTimeToTLim");
12 TF1 *fit=hTimeToTFit->GetFunction("pol5");
13
14 // Slewing parameters (same for all channels)
15
16 Float_t par[6] = {0.,0.,0.,0.,0.,0.};
17 Float_t kpar=0;
18 for(Int_t i =0;i<6;i++){
19 par[i]=fit->GetParameter(i);
20 kpar=par[0];
21 cout << " Slewing parameters=" << par[i] << endl;
22 }
23
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;
29
30 // ToT spectrum
31
32 TH1F *hToT= (TH1F*)f.Get("hToTLim");
33
34 // Fill the Sim calibration object
35
36 AliCDBManager *man = AliCDBManager::Instance();
37 man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
38 TRandom *rnd = new TRandom(4357);
39 Int_t nChannels = AliTOFGeometry::NSectors()*(2*(AliTOFGeometry::NStripC()+AliTOFGeometry::NStripB())+AliTOFGeometry::NStripA())*AliTOFGeometry::NpadZ()*AliTOFGeometry::NpadX();
40 for (Int_t ipad = 0 ; ipad<nChannels; ipad++){
41 AliTOFChannelOffline *calChannelOffline = (AliTOFChannelOffline*)tofCalOffline->At(ipad);
42 delay=rnd->Gaus(meanDelay,sigmaDelay);
43 par[0]=kpar+delay;
44 calChannelOffline->SetSlewPar(par);
45 }
46 tofcalib->WriteParOfflineOnCDB("TOF/Calib","valid",0,AliCDBRunRange::Infinity());
47 tofcalib->WriteSimHistoOnCDB("TOF/Calib",0,AliCDBRunRange::Infinity(),hToT);
48 f.Close();
49}
50
51