[u/mrichter/AliRoot.git] / TPC / CalibMacros / CalibEnv.C

/*
// .x ~/NimStyle.C
// .x ~/rootlogon.C

gSystem->AddIncludePath("-I$ALICE_ROOT/TPC");


.L $ALICE_ROOT/TPC/CalibMacros/CalibEnv.C+
Init();
CalibEnv("run.list");
GetTree();

TFile f("dcsTime.root")

//
// if you want to use alien OCDB
// 
gSystem->Load("libXrdClient.so");
gSystem->Load("libNetx.so");
if (!gGrid) TGrid::Connect("alien://",0,0,"t");


*/

#include <iostream>
#include <fstream>
#include <stdio.h>
#include <AliCDBManager.h>
#include <AliCDBEntry.h>
#include <AliLog.h>
#include <AliMagF.h>
#include "AliTPCcalibDB.h"
#include "AliTPCcalibDButil.h"
#include "AliTPCAltroMapping.h"
#include "AliTPCExB.h"
#include "AliTPCCalROC.h"
#include "AliTPCCalPad.h"
#include "AliTPCSensorTempArray.h"
#include "AliGRPObject.h"
#include "AliTPCTransform.h"
#include "TFile.h"
#include "TKey.h"
#include "TObjArray.h"
#include "TObjString.h"
#include "TString.h"
#include "AliTPCCalPad.h"
#include "AliTPCROC.h"
#include "AliTPCParam.h"
#include "AliTPCCalibPulser.h"
#include "AliTPCCalibPedestal.h"
#include "AliTPCCalibCE.h"
#include "AliTPCExBFirst.h"
#include "TTreeStream.h"
#include "AliTPCTempMap.h"
#include "TVectorD.h"
#include "TMatrixD.h"
#include "AliTPCCalibRaw.h"
#include "AliSplineFit.h"


TTree * dcsTree=0;
TString refFile="/data/Work/data/calib/guiTrees/RefCalPads_83680.root";
// TString refFile="/lustre/alice/wiechula/calib/guiTrees/RefCalPads_83680.root"
void GetProductionInfo(Int_t run, Int_t &nalien, Int_t &nRawAlien, Int_t &nlocal, Int_t &nRawLocal);
  
void Init2008(){
  //
  //
  //
  AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Parameters","local://$ALICE_ROOT/OCDB");
  AliCDBManager::Instance()->SetSpecificStorage("GRP/GRP/Data","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Temperature","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/HighVoltage","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Goofie","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
  AliCDBManager::Instance()->SetRun(1);
}

void Init(){  
  AliCDBManager::Instance()->SetDefaultStorage("local:///lustre/alice/alien/alice/data/2009/OCDB/");
  AliCDBManager::Instance()->SetRun(1);
}

void InitAlien(const char *path="LHC08b"){
  //
  //
  //
  TString alpath="alien://folder=/alice/data/2008/";
  alpath+=path;
  alpath+="/OCDB";
  
  AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Parameters","local://$ALICE_ROOT/OCDB");
  AliCDBManager::Instance()->SetSpecificStorage("GRP/GRP/Data",alpath.Data());
  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Temperature",alpath.Data());
  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Goofie",alpath.Data());
  AliCDBManager::Instance()->SetRun(1);
}


void CalibEnv(const char * runList){
  //
  //
  //
  AliTPCcalibDB * calibDB = AliTPCcalibDB::Instance();
  ifstream in;
  in.open(runList);
  Int_t irun=0;
  TTreeSRedirector *pcstream = new TTreeSRedirector("dcsTime.root");
  AliTPCcalibDButil dbutil;
  dbutil.SetRefFile(refFile.Data());
  Int_t startTime = 0;
  Int_t endTime = 0;
  Int_t startTimeGRP = 0;
  Int_t stopTimeGRP  = 0;

  AliSplineFit *fitVdrift=0x0;
  //  for (Int_t irun=startRun; irun<stopRun; irun++){
  while(in.good()) {
    in >> irun;
    if (in.eof()) break;
    if (irun==0) continue;
    printf("Processing run %d ...\n",irun);
    AliTPCcalibDB::Instance()->SetRun(irun);
    dbutil.UpdateFromCalibDB();
    AliDCSSensorArray *arrHV=calibDB->GetVoltageSensors(irun);
    if (!arrHV) continue;
    for  (Int_t isenHV=0; isenHV<arrHV->NumSensors(); ++isenHV){
      AliDCSSensor *senHV=arrHV->GetSensorNum(isenHV);
      if (!senHV) continue;
      startTime=senHV->GetStartTime();
      endTime  =senHV->GetEndTime();
      if (startTime>0&&endTime>0) break;
    }
    if (calibDB->GetGRP(irun)){
      startTimeGRP = AliTPCcalibDB::GetGRP(irun)->GetTimeStart();
      stopTimeGRP  = AliTPCcalibDB::GetGRP(irun)->GetTimeEnd();
    }
    //    AliDCSSensor * sensorPressure = AliTPCcalibDB::Instance()->GetPressureSensor(irun);
//     if (!sensorPressure) continue;
//     Int_t startTime = sensorPressure->GetStartTime();
//     Int_t endTime = sensorPressure->GetEndTime();
//     Int_t startTimeGRP = AliTPCcalibDB::GetGRP(irun)->GetTimeStart();
//     Int_t stopTimeGRP  = AliTPCcalibDB::GetGRP(irun)->GetTimeEnd();
    AliTPCSensorTempArray * tempArray = AliTPCcalibDB::Instance()->GetTemperatureSensor(irun);
    AliTPCTempMap * tempMap = new AliTPCTempMap(tempArray);
    AliDCSSensorArray* goofieArray = AliTPCcalibDB::Instance()->GetGoofieSensors(irun);
    //
    Int_t dtime = TMath::Max((endTime-startTime)/20,10*60);
    //Goofie statistical data
    TVectorD vecEntries, vecMean, vecMedian,vecRMS;
    dbutil.ProcessGoofie(vecEntries ,vecMedian, vecMean, vecRMS);
    //CE data processing - see ProcessCEdata function for description of the results
    TVectorD fitResultsA, fitResultsC;
    Int_t nmaskedCE;
//     dbutil.ProcessCEdata("(sector<36)++gx++gy++lx++lx**2",fitResultsA,fitResultsC,nmaskedCE);
    dbutil.ProcessCEdata("(sector<36)++gy++gx++(lx-134)++(sector<36)*(lx-134)",fitResultsA,fitResultsC,nmaskedCE);
    TVectorD vecTEntries, vecTMean, vecTRMS, vecTMedian, vecQEntries, vecQMean, vecQRMS, vecQMedian;
    Float_t driftTimeA, driftTimeC;
    dbutil.ProcessCEgraphs(vecTEntries, vecTMean, vecTRMS, vecTMedian,
                           vecQEntries, vecQMean, vecQRMS, vecQMedian,
                           driftTimeA, driftTimeC );
    //drift velocity using tracks
    fitVdrift=calibDB->GetVdriftSplineFit("ALISPLINEFIT_MEAN_VDRIFT_COSMICS_ALL",irun);
    //noise data Processing - see ProcessNoiseData function for description of the results
    TVectorD vNoiseMean, vNoiseMeanSenRegions, vNoiseRMS, vNoiseRMSSenRegions;
    Int_t nonMaskedZero=0;
    dbutil.ProcessNoiseData(vNoiseMean, vNoiseMeanSenRegions, vNoiseRMS, vNoiseRMSSenRegions, nonMaskedZero);
    // comparisons
    TVectorF pedestalDeviations;
    TVectorF noiseDeviations;
    TVectorF pulserQdeviations;
    Float_t varQMean;
    Int_t npadsOutOneTB;
    Int_t npadsOffAdd;
    dbutil.ProcessPedestalVariations(pedestalDeviations);
    dbutil.ProcessNoiseVariations(noiseDeviations);
    dbutil.ProcessPulserVariations(pulserQdeviations,varQMean,npadsOutOneTB,npadsOffAdd);
    
    //L3 data 
    Float_t bz=AliTPCcalibDB::GetBz(irun);
    Char_t  l3pol=AliTPCcalibDB::GetL3Polarity(irun);
    //calibration Pulser data processing
    Int_t nOffChannels=0;
    TVectorD vTimePulser;
    nOffChannels=dbutil.GetNPulserOutliers();
    dbutil.ProcessPulser(vTimePulser);
    //ALTRO data
    Int_t nMasked=0;
    dbutil.ProcessALTROConfig(nMasked);
    //Calib RAW data
    Int_t nFailL1=-1;
    if (calibDB->GetCalibRaw())
      nFailL1=calibDB->GetCalibRaw()->GetNFailL1Phase();
    //production information
    Int_t nalien=0,nRawAlien=0,nlocal=0,nRawLocal=0;
//     GetProductionInfo(irun, nalien, nRawAlien, nlocal,nRawLocal);
    //run type
    TObjString runType(AliTPCcalibDB::GetRunType(irun).Data());
    
    for (Int_t itime=startTime; itime<endTime; itime+=dtime){
      //
      TTimeStamp tstamp(itime);
      Float_t valuePressure  = calibDB->GetPressure(tstamp,irun,0);
      Float_t valuePressure2 = calibDB->GetPressure(tstamp,irun,1);
      //temperature fits
      TLinearFitter * fitter = 0;
      TVectorD vecTemp[10];
      for (Int_t itype=0; itype<5; itype++)
        for (Int_t iside=0; iside<2; iside++){
          fitter= tempMap->GetLinearFitter(itype,iside,tstamp);	  
          if (!fitter) continue;
          fitter->Eval();
          fitter->GetParameters(vecTemp[itype+iside*5]);
          delete fitter;
        }
      //measured skirt temperatures
      TVectorD vecSkirtTempA(18);
      TVectorD vecSkirtTempC(18);
      Int_t nsenTemp=tempArray->NumSensors();
      for (Int_t isenTemp=0;isenTemp<nsenTemp;++isenTemp){
        AliTPCSensorTemp *senTemp=(AliTPCSensorTemp*)tempArray->GetSensorNum(isenTemp);
        if (senTemp->GetType()!=3) continue;
        if (TMath::Sqrt(senTemp->GetX()*senTemp->GetX()+senTemp->GetY()*senTemp->GetY())<100) continue; //only skirt, outer FC vessel
        Double_t val=senTemp->GetValue(tstamp);
        if (senTemp->GetSide()==0)
          vecSkirtTempA[senTemp->GetSector()]=val;
        else
          vecSkirtTempC[senTemp->GetSector()]=val;
      }
      //goofie data
      TVectorD vecGoofie;
      if (goofieArray){
        vecGoofie.ResizeTo(goofieArray->NumSensors());
        for (Int_t isensor=0; isensor<goofieArray->NumSensors();isensor++){
          AliDCSSensor *gsensor = goofieArray->GetSensor(isensor);
          if (gsensor){
            vecGoofie[isensor] = gsensor->GetValue(tstamp);
          }
        }
      } else {
        vecGoofie.ResizeTo(19);
      }
      Double_t ptrelative0 = AliTPCcalibDB::GetPTRelative((UInt_t)itime,irun,0);
      Double_t ptrelative1 = AliTPCcalibDB::GetPTRelative((UInt_t)itime,irun,1);
      //
      TVectorD voltagesIROC(36);
      TVectorD voltagesOROC(36);
      for(Int_t j=1; j<36; j++) voltagesIROC[j-1] = AliTPCcalibDB::Instance()->GetChamberHighVoltage(irun, j,itime);
      for(Int_t j=36; j<72; j++) voltagesOROC[j-36] = AliTPCcalibDB::Instance()->GetChamberHighVoltage(irun, j,itime);
      Double_t voltIROC = TMath::Median(36, voltagesIROC.GetMatrixArray());
      Double_t voltOROC = TMath::Median(36, voltagesOROC.GetMatrixArray());
      //
      Float_t  coverIA=AliTPCcalibDB::GetCoverVoltage(irun,0,itime);
      Float_t  coverIC=AliTPCcalibDB::GetCoverVoltage(irun,18,itime);
      Float_t  coverOA=AliTPCcalibDB::GetCoverVoltage(irun,36,itime);
      Float_t  coverOC=AliTPCcalibDB::GetCoverVoltage(irun,54,itime);
      Float_t  skirtA=AliTPCcalibDB::GetSkirtVoltage(irun,0,itime);
      Float_t  skirtC=AliTPCcalibDB::GetSkirtVoltage(irun,18,itime);
      Float_t  ggOffA=AliTPCcalibDB::GetGGoffsetVoltage(irun,0,itime);
      Float_t  ggOffC=AliTPCcalibDB::GetGGoffsetVoltage(irun,18,itime);
      //drift velocity
      Float_t dvCorr=-5;
      if (fitVdrift) dvCorr=fitVdrift->Eval(itime);
      
      
      //tempMap->GetLinearFitter(0,0,itime);
      (*pcstream)<<"dcs"<<
        "run="<<irun<<
        "time="<<itime<<
        "startTimeGRP="<<startTimeGRP<<
        "stopTimeGRP="<<stopTimeGRP<<
        //run type
        "runType.="<<&runType<<
        // voltage setting
        "VIROC.="<<&voltagesIROC<<
        "VOROC.="<<&voltagesOROC<<
        "medianVIROC="<<voltIROC<<
        "medianVOROC="<<voltOROC<<
        "coverIA=" << coverIA <<
        "coverIC=" << coverIC <<
        "coverOA=" << coverOA <<
        "coverOC=" << coverOC <<
        "skirtA=" << skirtA <<
        "skirtC=" << skirtC <<
        "ggOffA=" << ggOffA <<
        "ggOffC=" << ggOffC <<
        //
        "ptrel0="<<ptrelative0<<  // deltaTP/TP  - A side
        "ptrel1="<<ptrelative1<<  // deltaTP/TPC - C side
        "goofie.="<<&vecGoofie<<
        "goofieE.="<<&vecEntries<<
        "goofieMean.="<<&vecMean<<
        "goofieMedian.="<<&vecMedian<<
        "goofieRMS.="<<&vecRMS<<
        //
        "press="<<valuePressure<<
        "press2="<<valuePressure2<<
        "temp00.="<<&vecTemp[0]<<
        "temp10.="<<&vecTemp[1]<<
        "temp20.="<<&vecTemp[2]<<
        "temp30.="<<&vecTemp[3]<<
        "temp40.="<<&vecTemp[4]<<
        "temp01.="<<&vecTemp[5]<<
        "temp11.="<<&vecTemp[6]<<
        "temp21.="<<&vecTemp[7]<<
        "temp31.="<<&vecTemp[8]<<
        "temp41.="<<&vecTemp[9]<<
        "tempSkirtA.="<<&vecSkirtTempA<<
        "tempSkirtC.="<<&vecSkirtTempC;
      
      (*pcstream)<<"dcs"<<	
        //noise data
        "meanNoise.="<<&vNoiseMean<<
        "meanNoiseSen.="<<&vNoiseMeanSenRegions<<
        "rmsNoise.="<<&vNoiseRMS<<
        "rmsNoiseSen.="<<&vNoiseRMSSenRegions<<
        "zeroNoise="<<nonMaskedZero<<
        //pulser data
        "timePulser.=" << &vTimePulser <<
        "nOffPulser="<<nOffChannels<<
        //altro data
        "nMasked="<<nMasked<<
        //ce data
        "CEfitA.="<<&fitResultsA<<
        "CEfitC.="<<&fitResultsC<<
        "nmaskedCE="<<nmaskedCE<<
        //ce graph data
        "CEgrTEntries.="<<&vecTEntries<<
        "CEgrTMean.="<<&vecTMean<<
        "CEgrTRMS.="<<&vecTRMS<<
        "CEgrTMedian.="<<&vecTMedian<<
        "CEgrQEntries.="<<&vecQEntries<<
        "CEgrQMean.="<<&vecQMean<<
        "CEgrQRMS.="<<&vecQRMS<<
        "CEgrQMedian.="<<&vecQMedian<<
        "CEgrDriftA="<<driftTimeA<<
        "CEgrDriftC="<<driftTimeC<<
        //calib raw data
        "nFailL1="<<nFailL1<<
        // b field
        "Bz="<< bz <<
        "L3polarity="<<l3pol<<
        // production information
        "nalien="<<nalien<<
        "nRawAlien="<<nRawAlien<<
        "nlocal="<<nlocal<<
        "nRawLocal="<<nRawLocal<<
        //comparisons with ref data
        "pedestalDeviations.="<<&pedestalDeviations<<
        "noiseDeviations.="<<&noiseDeviations<<
        "pulserQdeviations.="<<&pulserQdeviations<<
//         "pulserVarQMean="<<varQMean<<
        "pulserNpadsOutOneTB="<<npadsOutOneTB<<
        "pulserNpadsOffAdd="<<npadsOffAdd<<
        "driftCorrCosmAll="<<dvCorr<<
      
        "\n";
    }
  }
  delete pcstream;
}


void GetProductionInfo(Int_t run, Int_t &nalien, Int_t &nRawAlien, Int_t &nlocal, Int_t &nRawLocal){
  //
  // find number of ESDs in central and local production for run
  //

  nalien=0;
  nRawAlien=0;
  nlocal=0;
  nRawLocal=0;
  TString sNlines;
  //find number of ESDs in alien
  TString command="alien_find /alice/data/2009 ";
  command += Form("%09d",run);
  command += " | grep AliESDs.root | wc -l";
  sNlines = gSystem->GetFromPipe(command.Data());
  nalien=sNlines.Atoi();
  //find number of raw files on alien
  command="alien_find /alice/data/2009 ";
  command += Form("%09d",run);
  command += " | grep raw | grep -v tag | wc -l";
  sNlines = gSystem->GetFromPipe(command.Data());
  nRawAlien=sNlines.Atoi();
  //find number of ESDs local
  command="find /lustre/alice/alien/alice/data/2009 -name AliESDs.root | grep ";
  command += Form("%09d",run);
  command += " | wc -l";
  sNlines = gSystem->GetFromPipe(command.Data());
  nlocal=sNlines.Atoi();
  //find number of local raw data files
  command="find /lustre/alice/alien/alice/data/2009 -name \"*.root\" | grep ";
  command += Form("%09d",run);
  command += " | grep raw | grep -v tag | wc -l";
  sNlines = gSystem->GetFromPipe(command.Data());
  nRawLocal=sNlines.Atoi();
}

void FilterMag(const char * runList){
  //
  //
  //
  //  AliTPCcalibDB * calibDB = AliTPCcalibDB::Instance();
  ifstream in;
  in.open(runList);
  Int_t irun=0;
  while(in.good()) {
    in >> irun;
    if (irun==0) continue;
    AliGRPObject *grp = AliTPCcalibDB::GetGRP(irun);
    Float_t current = -1;
    Float_t bz      = -1;
//     Float_t press   =  0;
    if (grp){
      current = grp->GetL3Current((AliGRPObject::Stats)0);
      bz = 5*current/30000.;
      printf("Run%d\tL3 current%f\tBz\t%f\n",irun,current,bz);
    }
    else{
      printf("Run%d\tL3 current%f\tBz\t%f\n",irun,current,bz);
    }
  }
  
}


void GetTree(){
  TFile *fdcs = new TFile("dcsTime.root");
  dcsTree  = (TTree*)fdcs->Get("dcs");
  //
  // mean temp A
  
  dcsTree->Draw("temp30.fElements[0]");
  
}

void GetNominalValues(){
  //
  if (!dcsTree) return;
}


/*

AliDCSSensor * sensorPressure = AliTPCcalibDB::Instance()->GetPressureSensor(62084);
entry = AliCDBManager::Instance()->Get("TPC/Calib/Temperature",run);
AliTPCSensorTempArray * tempArray = (AliTPCSensorTempArray *)entry->GetObject();
AliTPCSensorTempArray * tempArray = (AliTPCSensorTempArray *)AliTPCcalibDB::Instance()->GetTemperatureSensor(62084)
AliTPCTempMap * tempMap = new AliTPCTempMap(tempArray);
TLinearFitter * fitter = tempMap->GetLinearFitter(0,0,tempArray->GetStartTime());

AliDCSSensorArray* goofieArray = AliTPCcalibDB::Instance()->GetGoofieSensors(62084);

*/

/*

void PlotPressureResol(){
  //
  // Example
  //
  dcs->Draw("100000*(press-press2-4.782)/press/sqrt(2.)>>his(100,-50,50)","run>61400","")
  his->SetXTitle("#sigma_{P/P_{0}}(10^{-5})");
  gPad->SaveAs("picDCS/deltaPoverP.eps");
  gPad->SaveAs("picDCS/deltaPoverP.gif");

}
void PlotTresol(){
  //
  // T resolution example
  // plot difference of the temperature from A and C side
  // Supposing the error is independent - (division by sqrt(2))
  dcs->Draw("100000*(temp30.fElements[0]-temp31.fElements[0]+0.00509)/(temp31.fElements[0]+273.15)/sqrt(2.)>>his(100,-5,5)","run>61400","");
  his->SetXTitle("#sigma_{T/T_{0}}(10^{-5})");
  gPad->SaveAs("picDCS/deltaToverT.eps");
  gPad->SaveAs("picDCS/deltaToverT.gif");
}
*/
Commit	Line	Data
d47d42ee	1	/*
7390f655	2	// .x ~/NimStyle.C
7390f655	3	// .x ~/rootlogon.C
7ce80437	4
d47d42ee	5	gSystem->AddIncludePath("-I$ALICE_ROOT/TPC");
e2914767	6
e2914767	7
e2914767	8
7ce80437	9	.L $ALICE_ROOT/TPC/CalibMacros/CalibEnv.C+
d47d42ee	10	Init();
6e7d7dc4	11	CalibEnv("run.list");
7ce80437	12	GetTree();
7ce80437	13
e79211e8	14	TFile f("dcsTime.root")
d47d42ee	15
7ce80437	16	//
	17	// if you want to use alien OCDB
	18	//
	19	gSystem->Load("libXrdClient.so");
	20	gSystem->Load("libNetx.so");
	21	if (!gGrid) TGrid::Connect("alien://",0,0,"t");
	22
	23
e2914767	24
d47d42ee	25	*/
	26
	27	#include <iostream>
	28	#include <fstream>
5312f439	29	#include <stdio.h>
d47d42ee	30	#include <AliCDBManager.h>
	31	#include <AliCDBEntry.h>
	32	#include <AliLog.h>
	33	#include <AliMagF.h>
	34	#include "AliTPCcalibDB.h"
892226be	35	#include "AliTPCcalibDButil.h"
d47d42ee	36	#include "AliTPCAltroMapping.h"
	37	#include "AliTPCExB.h"
	38	#include "AliTPCCalROC.h"
	39	#include "AliTPCCalPad.h"
	40	#include "AliTPCSensorTempArray.h"
	41	#include "AliGRPObject.h"
	42	#include "AliTPCTransform.h"
	43	#include "TFile.h"
	44	#include "TKey.h"
	45	#include "TObjArray.h"
	46	#include "TObjString.h"
	47	#include "TString.h"
	48	#include "AliTPCCalPad.h"
5312f439	49	#include "AliTPCROC.h"
5312f439	50	#include "AliTPCParam.h"
d47d42ee	51	#include "AliTPCCalibPulser.h"
	52	#include "AliTPCCalibPedestal.h"
	53	#include "AliTPCCalibCE.h"
	54	#include "AliTPCExBFirst.h"
	55	#include "TTreeStream.h"
	56	#include "AliTPCTempMap.h"
5312f439	57	#include "TVectorD.h"
5312f439	58	#include "TMatrixD.h"
6e7d7dc4	59	#include "AliTPCCalibRaw.h"
7390f655	60	#include "AliSplineFit.h"
d47d42ee	61
7ce80437	62
7ce80437	63	TTree * dcsTree=0;
7390f655	64	TString refFile="/data/Work/data/calib/guiTrees/RefCalPads_83680.root";
7390f655	65	// TString refFile="/lustre/alice/wiechula/calib/guiTrees/RefCalPads_83680.root"
6fb51ca4	66	void GetProductionInfo(Int_t run, Int_t &nalien, Int_t &nRawAlien, Int_t &nlocal, Int_t &nRawLocal);
5312f439	67
6e7d7dc4	68	void Init2008(){
d47d42ee	69	//
	70	//
	71	//
162637e4	72	AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
162637e4	73	AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Parameters","local://$ALICE_ROOT/OCDB");
e6970ab5	74	AliCDBManager::Instance()->SetSpecificStorage("GRP/GRP/Data","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
	75	AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Temperature","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
	76	AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/HighVoltage","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
	77	AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Goofie","local:///lustre/alice/alien/alice/data/2008/LHC08d/OCDB/");
d47d42ee	78	AliCDBManager::Instance()->SetRun(1);
e2914767	79	}
e2914767	80
6e7d7dc4	81	void Init(){
	82	AliCDBManager::Instance()->SetDefaultStorage("local:///lustre/alice/alien/alice/data/2009/OCDB/");
	83	AliCDBManager::Instance()->SetRun(1);
	84	}
d47d42ee	85
e2914767	86	void InitAlien(const char *path="LHC08b"){
	87	//
	88	//
	89	//
	90	TString alpath="alien://folder=/alice/data/2008/";
	91	alpath+=path;
	92	alpath+="/OCDB";
7ce80437	93
162637e4	94	AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
162637e4	95	AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Parameters","local://$ALICE_ROOT/OCDB");
e2914767	96	AliCDBManager::Instance()->SetSpecificStorage("GRP/GRP/Data",alpath.Data());
	97	AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Temperature",alpath.Data());
	98	AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Goofie",alpath.Data());
	99	AliCDBManager::Instance()->SetRun(1);
d47d42ee	100	}
	101
	102
	103	void CalibEnv(const char * runList){
	104	//
	105	//
7ce80437	106	//
e79211e8	107	AliTPCcalibDB * calibDB = AliTPCcalibDB::Instance();
d47d42ee	108	ifstream in;
	109	in.open(runList);
	110	Int_t irun=0;
e79211e8	111	TTreeSRedirector *pcstream = new TTreeSRedirector("dcsTime.root");
892226be	112	AliTPCcalibDButil dbutil;
7390f655	113	dbutil.SetRefFile(refFile.Data());
	114	Int_t startTime = 0;
	115	Int_t endTime = 0;
	116	Int_t startTimeGRP = 0;
	117	Int_t stopTimeGRP = 0;
	118
	119	AliSplineFit *fitVdrift=0x0;
d47d42ee	120	// for (Int_t irun=startRun; irun<stopRun; irun++){
	121	while(in.good()) {
	122	in >> irun;
6e7d7dc4	123	if (in.eof()) break;
d47d42ee	124	if (irun==0) continue;
e2914767	125	printf("Processing run %d ...\n",irun);
5312f439	126	AliTPCcalibDB::Instance()->SetRun(irun);
892226be	127	dbutil.UpdateFromCalibDB();
7390f655	128	AliDCSSensorArray *arrHV=calibDB->GetVoltageSensors(irun);
	129	if (!arrHV) continue;
	130	for (Int_t isenHV=0; isenHV<arrHV->NumSensors(); ++isenHV){
	131	AliDCSSensor *senHV=arrHV->GetSensorNum(isenHV);
	132	if (!senHV) continue;
	133	startTime=senHV->GetStartTime();
	134	endTime =senHV->GetEndTime();
	135	if (startTime>0&&endTime>0) break;
	136	}
	137	if (calibDB->GetGRP(irun)){
	138	startTimeGRP = AliTPCcalibDB::GetGRP(irun)->GetTimeStart();
	139	stopTimeGRP = AliTPCcalibDB::GetGRP(irun)->GetTimeEnd();
	140	}
	141	// AliDCSSensor * sensorPressure = AliTPCcalibDB::Instance()->GetPressureSensor(irun);
	142	// if (!sensorPressure) continue;
	143	// Int_t startTime = sensorPressure->GetStartTime();
	144	// Int_t endTime = sensorPressure->GetEndTime();
	145	// Int_t startTimeGRP = AliTPCcalibDB::GetGRP(irun)->GetTimeStart();
	146	// Int_t stopTimeGRP = AliTPCcalibDB::GetGRP(irun)->GetTimeEnd();
d47d42ee	147	AliTPCSensorTempArray * tempArray = AliTPCcalibDB::Instance()->GetTemperatureSensor(irun);
	148	AliTPCTempMap * tempMap = new AliTPCTempMap(tempArray);
	149	AliDCSSensorArray* goofieArray = AliTPCcalibDB::Instance()->GetGoofieSensors(irun);
	150	//
d47d42ee	151	Int_t dtime = TMath::Max((endTime-startTime)/20,10*60);
892226be	152	//Goofie statistical data
	153	TVectorD vecEntries, vecMean, vecMedian,vecRMS;
	154	dbutil.ProcessGoofie(vecEntries ,vecMedian, vecMean, vecRMS);
5312f439	155	//CE data processing - see ProcessCEdata function for description of the results
5312f439	156	TVectorD fitResultsA, fitResultsC;
892226be	157	Int_t nmaskedCE;
732e90a8	158	// dbutil.ProcessCEdata("(sector<36)++gx++gy++lx++lx**2",fitResultsA,fitResultsC,nmaskedCE);
732e90a8	159	dbutil.ProcessCEdata("(sector<36)++gy++gx++(lx-134)++(sector<36)*(lx-134)",fitResultsA,fitResultsC,nmaskedCE);
892226be	160	TVectorD vecTEntries, vecTMean, vecTRMS, vecTMedian, vecQEntries, vecQMean, vecQRMS, vecQMedian;
	161	Float_t driftTimeA, driftTimeC;
	162	dbutil.ProcessCEgraphs(vecTEntries, vecTMean, vecTRMS, vecTMedian,
	163	vecQEntries, vecQMean, vecQRMS, vecQMedian,
	164	driftTimeA, driftTimeC );
7390f655	165	//drift velocity using tracks
7390f655	166	fitVdrift=calibDB->GetVdriftSplineFit("ALISPLINEFIT_MEAN_VDRIFT_COSMICS_ALL",irun);
5312f439	167	//noise data Processing - see ProcessNoiseData function for description of the results
	168	TVectorD vNoiseMean, vNoiseMeanSenRegions, vNoiseRMS, vNoiseRMSSenRegions;
	169	Int_t nonMaskedZero=0;
892226be	170	dbutil.ProcessNoiseData(vNoiseMean, vNoiseMeanSenRegions, vNoiseRMS, vNoiseRMSSenRegions, nonMaskedZero);
7390f655	171	// comparisons
	172	TVectorF pedestalDeviations;
	173	TVectorF noiseDeviations;
	174	TVectorF pulserQdeviations;
	175	Float_t varQMean;
	176	Int_t npadsOutOneTB;
	177	Int_t npadsOffAdd;
	178	dbutil.ProcessPedestalVariations(pedestalDeviations);
	179	dbutil.ProcessNoiseVariations(noiseDeviations);
	180	dbutil.ProcessPulserVariations(pulserQdeviations,varQMean,npadsOutOneTB,npadsOffAdd);
	181
e6970ab5	182	//L3 data
5312f439	183	Float_t bz=AliTPCcalibDB::GetBz(irun);
	184	Char_t l3pol=AliTPCcalibDB::GetL3Polarity(irun);
	185	//calibration Pulser data processing
5312f439	186	Int_t nOffChannels=0;
892226be	187	TVectorD vTimePulser;
	188	nOffChannels=dbutil.GetNPulserOutliers();
	189	dbutil.ProcessPulser(vTimePulser);
	190	//ALTRO data
	191	Int_t nMasked=0;
	192	dbutil.ProcessALTROConfig(nMasked);
6e7d7dc4	193	//Calib RAW data
	194	Int_t nFailL1=-1;
	195	if (calibDB->GetCalibRaw())
	196	nFailL1=calibDB->GetCalibRaw()->GetNFailL1Phase();
5312f439	197	//production information
6fb51ca4	198	Int_t nalien=0,nRawAlien=0,nlocal=0,nRawLocal=0;
e6970ab5	199	// GetProductionInfo(irun, nalien, nRawAlien, nlocal,nRawLocal);
5312f439	200	//run type
	201	TObjString runType(AliTPCcalibDB::GetRunType(irun).Data());
	202
d47d42ee	203	for (Int_t itime=startTime; itime<endTime; itime+=dtime){
d47d42ee	204	//
5312f439	205	TTimeStamp tstamp(itime);
e79211e8	206	Float_t valuePressure = calibDB->GetPressure(tstamp,irun,0);
e79211e8	207	Float_t valuePressure2 = calibDB->GetPressure(tstamp,irun,1);
5312f439	208	//temperature fits
d47d42ee	209	TLinearFitter * fitter = 0;
d47d42ee	210	TVectorD vecTemp[10];
e79211e8	211	for (Int_t itype=0; itype<5; itype++)
5312f439	212	for (Int_t iside=0; iside<2; iside++){
	213	fitter= tempMap->GetLinearFitter(itype,iside,tstamp);
	214	if (!fitter) continue;
	215	fitter->Eval();
	216	fitter->GetParameters(vecTemp[itype+iside*5]);
	217	delete fitter;
	218	}
	219	//measured skirt temperatures
	220	TVectorD vecSkirtTempA(18);
	221	TVectorD vecSkirtTempC(18);
	222	Int_t nsenTemp=tempArray->NumSensors();
	223	for (Int_t isenTemp=0;isenTemp<nsenTemp;++isenTemp){
	224	AliTPCSensorTemp senTemp=(AliTPCSensorTemp)tempArray->GetSensorNum(isenTemp);
	225	if (senTemp->GetType()!=3) continue;
	226	if (TMath::Sqrt(senTemp->GetX()senTemp->GetX()+senTemp->GetY()senTemp->GetY())<100) continue; //only skirt, outer FC vessel
	227	Double_t val=senTemp->GetValue(tstamp);
	228	if (senTemp->GetSide()==0)
	229	vecSkirtTempA[senTemp->GetSector()]=val;
	230	else
	231	vecSkirtTempC[senTemp->GetSector()]=val;
	232	}
892226be	233	//goofie data
892226be	234	TVectorD vecGoofie;
5312f439	235	if (goofieArray){
5312f439	236	vecGoofie.ResizeTo(goofieArray->NumSensors());
5312f439	237	for (Int_t isensor=0; isensor<goofieArray->NumSensors();isensor++){
	238	AliDCSSensor *gsensor = goofieArray->GetSensor(isensor);
	239	if (gsensor){
	240	vecGoofie[isensor] = gsensor->GetValue(tstamp);
	241	}
	242	}
892226be	243	} else {
892226be	244	vecGoofie.ResizeTo(19);
d47d42ee	245	}
892226be	246	Double_t ptrelative0 = AliTPCcalibDB::GetPTRelative((UInt_t)itime,irun,0);
892226be	247	Double_t ptrelative1 = AliTPCcalibDB::GetPTRelative((UInt_t)itime,irun,1);
e2914767	248	//
5312f439	249	TVectorD voltagesIROC(36);
5312f439	250	TVectorD voltagesOROC(36);
892226be	251	for(Int_t j=1; j<36; j++) voltagesIROC[j-1] = AliTPCcalibDB::Instance()->GetChamberHighVoltage(irun, j,itime);
892226be	252	for(Int_t j=36; j<72; j++) voltagesOROC[j-36] = AliTPCcalibDB::Instance()->GetChamberHighVoltage(irun, j,itime);
7ce80437	253	Double_t voltIROC = TMath::Median(36, voltagesIROC.GetMatrixArray());
	254	Double_t voltOROC = TMath::Median(36, voltagesOROC.GetMatrixArray());
	255	//
5312f439	256	Float_t coverIA=AliTPCcalibDB::GetCoverVoltage(irun,0,itime);
	257	Float_t coverIC=AliTPCcalibDB::GetCoverVoltage(irun,18,itime);
	258	Float_t coverOA=AliTPCcalibDB::GetCoverVoltage(irun,36,itime);
	259	Float_t coverOC=AliTPCcalibDB::GetCoverVoltage(irun,54,itime);
	260	Float_t skirtA=AliTPCcalibDB::GetSkirtVoltage(irun,0,itime);
	261	Float_t skirtC=AliTPCcalibDB::GetSkirtVoltage(irun,18,itime);
	262	Float_t ggOffA=AliTPCcalibDB::GetGGoffsetVoltage(irun,0,itime);
	263	Float_t ggOffC=AliTPCcalibDB::GetGGoffsetVoltage(irun,18,itime);
7390f655	264	//drift velocity
	265	Float_t dvCorr=-5;
	266	if (fitVdrift) dvCorr=fitVdrift->Eval(itime);
5312f439	267
5312f439	268
e2914767	269	//tempMap->GetLinearFitter(0,0,itime);
e6970ab5	270	(*pcstream)<<"dcs"<<
5312f439	271	"run="<<irun<<
5312f439	272	"time="<<itime<<
6fb51ca4	273	"startTimeGRP="<<startTimeGRP<<
6fb51ca4	274	"stopTimeGRP="<<stopTimeGRP<<
5312f439	275	//run type
	276	"runType.="<<&runType<<
	277	// voltage setting
	278	"VIROC.="<<&voltagesIROC<<
	279	"VOROC.="<<&voltagesOROC<<
	280	"medianVIROC="<<voltIROC<<
	281	"medianVOROC="<<voltOROC<<
7ce80437	282	"coverIA=" << coverIA <<
	283	"coverIC=" << coverIC <<
	284	"coverOA=" << coverOA <<
	285	"coverOC=" << coverOC <<
	286	"skirtA=" << skirtA <<
	287	"skirtC=" << skirtC <<
	288	"ggOffA=" << ggOffA <<
	289	"ggOffC=" << ggOffC <<
5312f439	290	//
	291	"ptrel0="<<ptrelative0<< // deltaTP/TP - A side
	292	"ptrel1="<<ptrelative1<< // deltaTP/TPC - C side
	293	"goofie.="<<&vecGoofie<<
	294	"goofieE.="<<&vecEntries<<
	295	"goofieMean.="<<&vecMean<<
	296	"goofieMedian.="<<&vecMedian<<
	297	"goofieRMS.="<<&vecRMS<<
	298	//
	299	"press="<<valuePressure<<
	300	"press2="<<valuePressure2<<
	301	"temp00.="<<&vecTemp[0]<<
	302	"temp10.="<<&vecTemp[1]<<
	303	"temp20.="<<&vecTemp[2]<<
	304	"temp30.="<<&vecTemp[3]<<
	305	"temp40.="<<&vecTemp[4]<<
	306	"temp01.="<<&vecTemp[5]<<
	307	"temp11.="<<&vecTemp[6]<<
	308	"temp21.="<<&vecTemp[7]<<
	309	"temp31.="<<&vecTemp[8]<<
	310	"temp41.="<<&vecTemp[9]<<
	311	"tempSkirtA.="<<&vecSkirtTempA<<
6e7d7dc4	312	"tempSkirtC.="<<&vecSkirtTempC;
	313
	314	(*pcstream)<<"dcs"<<
5312f439	315	//noise data
	316	"meanNoise.="<<&vNoiseMean<<
	317	"meanNoiseSen.="<<&vNoiseMeanSenRegions<<
	318	"rmsNoise.="<<&vNoiseRMS<<
	319	"rmsNoiseSen.="<<&vNoiseRMSSenRegions<<
	320	"zeroNoise="<<nonMaskedZero<<
	321	//pulser data
892226be	322	"timePulser.=" << &vTimePulser <<
5312f439	323	"nOffPulser="<<nOffChannels<<
892226be	324	//altro data
892226be	325	"nMasked="<<nMasked<<
5312f439	326	//ce data
	327	"CEfitA.="<<&fitResultsA<<
	328	"CEfitC.="<<&fitResultsC<<
892226be	329	"nmaskedCE="<<nmaskedCE<<
	330	//ce graph data
	331	"CEgrTEntries.="<<&vecTEntries<<
	332	"CEgrTMean.="<<&vecTMean<<
	333	"CEgrTRMS.="<<&vecTRMS<<
	334	"CEgrTMedian.="<<&vecTMedian<<
	335	"CEgrQEntries.="<<&vecQEntries<<
	336	"CEgrQMean.="<<&vecQMean<<
	337	"CEgrQRMS.="<<&vecQRMS<<
	338	"CEgrQMedian.="<<&vecQMedian<<
892226be	339	"CEgrDriftA="<<driftTimeA<<
892226be	340	"CEgrDriftC="<<driftTimeC<<
6e7d7dc4	341	//calib raw data
6e7d7dc4	342	"nFailL1="<<nFailL1<<
5312f439	343	// b field
	344	"Bz="<< bz <<
	345	"L3polarity="<<l3pol<<
	346	// production information
	347	"nalien="<<nalien<<
6fb51ca4	348	"nRawAlien="<<nRawAlien<<
5312f439	349	"nlocal="<<nlocal<<
6fb51ca4	350	"nRawLocal="<<nRawLocal<<
7390f655	351	//comparisons with ref data
	352	"pedestalDeviations.="<<&pedestalDeviations<<
	353	"noiseDeviations.="<<&noiseDeviations<<
	354	"pulserQdeviations.="<<&pulserQdeviations<<
	355	// "pulserVarQMean="<<varQMean<<
	356	"pulserNpadsOutOneTB="<<npadsOutOneTB<<
	357	"pulserNpadsOffAdd="<<npadsOffAdd<<
	358	"driftCorrCosmAll="<<dvCorr<<
	359
5312f439	360	"\n";
d47d42ee	361	}
	362	}
	363	delete pcstream;
	364	}
	365
	366
5312f439	367
5312f439	368
6fb51ca4	369	void GetProductionInfo(Int_t run, Int_t &nalien, Int_t &nRawAlien, Int_t &nlocal, Int_t &nRawLocal){
5312f439	370	//
	371	// find number of ESDs in central and local production for run
	372	//
	373
	374	nalien=0;
6fb51ca4	375	nRawAlien=0;
5312f439	376	nlocal=0;
6fb51ca4	377	nRawLocal=0;
5312f439	378	TString sNlines;
5312f439	379	//find number of ESDs in alien
	380	TString command="alien_find /alice/data/2009 ";
	381	command += Form("%09d",run);
	382	command += " \| grep AliESDs.root \| wc -l";
6fb51ca4	383	sNlines = gSystem->GetFromPipe(command.Data());
5312f439	384	nalien=sNlines.Atoi();
6fb51ca4	385	//find number of raw files on alien
	386	command="alien_find /alice/data/2009 ";
	387	command += Form("%09d",run);
	388	command += " \| grep raw \| grep -v tag \| wc -l";
	389	sNlines = gSystem->GetFromPipe(command.Data());
	390	nRawAlien=sNlines.Atoi();
5312f439	391	//find number of ESDs local
6fb51ca4	392	command="find /lustre/alice/alien/alice/data/2009 -name AliESDs.root \| grep ";
5312f439	393	command += Form("%09d",run);
6fb51ca4	394	command += " \| wc -l";
6fb51ca4	395	sNlines = gSystem->GetFromPipe(command.Data());
5312f439	396	nlocal=sNlines.Atoi();
6fb51ca4	397	//find number of local raw data files
	398	command="find /lustre/alice/alien/alice/data/2009 -name \"*.root\" \| grep ";
	399	command += Form("%09d",run);
	400	command += " \| grep raw \| grep -v tag \| wc -l";
	401	sNlines = gSystem->GetFromPipe(command.Data());
	402	nRawLocal=sNlines.Atoi();
5312f439	403	}
e2914767	404
	405	void FilterMag(const char * runList){
	406	//
	407	//
	408	//
	409	// AliTPCcalibDB * calibDB = AliTPCcalibDB::Instance();
	410	ifstream in;
	411	in.open(runList);
	412	Int_t irun=0;
	413	while(in.good()) {
	414	in >> irun;
	415	if (irun==0) continue;
	416	AliGRPObject *grp = AliTPCcalibDB::GetGRP(irun);
	417	Float_t current = -1;
	418	Float_t bz = -1;
5312f439	419	// Float_t press = 0;
e2914767	420	if (grp){
	421	current = grp->GetL3Current((AliGRPObject::Stats)0);
	422	bz = 5*current/30000.;
	423	printf("Run%d\tL3 current%f\tBz\t%f\n",irun,current,bz);
	424	}
	425	else{
	426	printf("Run%d\tL3 current%f\tBz\t%f\n",irun,current,bz);
	427	}
	428	}
7ce80437	429
	430	}
	431
	432
	433	void GetTree(){
	434	TFile *fdcs = new TFile("dcsTime.root");
5312f439	435	dcsTree = (TTree*)fdcs->Get("dcs");
7ce80437	436	//
7ce80437	437	// mean temp A
5312f439	438
7ce80437	439	dcsTree->Draw("temp30.fElements[0]");
5312f439	440
e2914767	441	}
e2914767	442
7ce80437	443	void GetNominalValues(){
	444	//
	445	if (!dcsTree) return;
	446	}
	447
	448
	449
	450
d47d42ee	451	/*
	452
	453	AliDCSSensor * sensorPressure = AliTPCcalibDB::Instance()->GetPressureSensor(62084);
	454	entry = AliCDBManager::Instance()->Get("TPC/Calib/Temperature",run);
	455	AliTPCSensorTempArray * tempArray = (AliTPCSensorTempArray *)entry->GetObject();
	456	AliTPCSensorTempArray * tempArray = (AliTPCSensorTempArray *)AliTPCcalibDB::Instance()->GetTemperatureSensor(62084)
	457	AliTPCTempMap * tempMap = new AliTPCTempMap(tempArray);
e79211e8	458	TLinearFitter * fitter = tempMap->GetLinearFitter(0,0,tempArray->GetStartTime());
d47d42ee	459
	460	AliDCSSensorArray* goofieArray = AliTPCcalibDB::Instance()->GetGoofieSensors(62084);
	461
	462	*/
7ce80437	463
	464	/*
	465
	466	void PlotPressureResol(){
	467	//
	468	// Example
	469	//
	470	dcs->Draw("100000*(press-press2-4.782)/press/sqrt(2.)>>his(100,-50,50)","run>61400","")
	471	his->SetXTitle("#sigma_{P/P_{0}}(10^{-5})");
	472	gPad->SaveAs("picDCS/deltaPoverP.eps");
	473	gPad->SaveAs("picDCS/deltaPoverP.gif");
5312f439	474
7ce80437	475	}
	476	void PlotTresol(){
	477	//
	478	// T resolution example
	479	// plot difference of the temperature from A and C side
	480	// Supposing the error is independent - (division by sqrt(2))
	481	dcs->Draw("100000*(temp30.fElements[0]-temp31.fElements[0]+0.00509)/(temp31.fElements[0]+273.15)/sqrt(2.)>>his(100,-5,5)","run>61400","");
	482	his->SetXTitle("#sigma_{T/T_{0}}(10^{-5})");
	483	gPad->SaveAs("picDCS/deltaToverT.eps");
	484	gPad->SaveAs("picDCS/deltaToverT.gif");
	485	}
	486	*/