string insereEspecie(int id, string nomeCientifico,
string nomePopular, string descricao){
especies especie(id,nomeCientifico,nomePopular,descricao);
if(gerenEspecies->busca(&especie) == -1){
gerenEspecies->insere(especie);
commit();
return string("status?v=1");
}
return string("status?v=0");
}
string removeEspecie(int id){
especies especie(id,"","","");
bool removeu = gerenEspecies->remove(especie);
if(removeu){
commit();
return string("status?v=1");
}
return string("status?v=0");
}
void PlotQASDD(Char_t fileName[100]="File.QA.111333.2010.LHC09b.pass2.root",Bool_t kDoEps=kFALSE,Char_t eventspecie[25]="LowMultiplicity")
{
const TString Rawsexpertname[]={"SDDModPattern",//0 1
"SDDModPatternNORM",//1 //1
"SDDphizL3NORM", //2 //0
"SDDphizL4NORM", //3 //0
"SDDL3_RelativeOccupancy", //4 4 5
"SDDL4_RelativeOccupancy",//5 4 5
"SDDRawDataCheck"//6 6
};//tot 7
const TString Rawsnonexpertname[]={"SDDphizL3",//0 2
"SDDphizL4",//1 2
"SDDphizCalibL3",//2 // 3
"SDDphizCalibL4"//3 3
};//tot 4
const TString RecPointsexpertname[]={"SDDLocalCoordDistrib",//0 //11
"SDDModPatternL3RP",//1 //12
"SDDModPatternL4RP",//2 //12
"SDDModPatternL3RPNORM",//3 //13
"SDDModPatternL4RPNORM",//4 //13
"SDDModPatternRP",//5 //14
"SDDModPatternRPNORM",//6 //14
"SDDphidistrib_Layer3",//7 //15
"SDDphidistrib_Layer4",//8 //15
"SDDrdistrib_Layer3",//9 //16
"SDDrdistrib_Layer4",//10 //16
"SDDL3_RelativeOccupancy",//11 //19 20
"SDDL4_RelativeOccupancy",//12 //19 20
"SDDL3_Rec2Raw_2D", //13 //21
"SDDL4_Rec2Raw_2D", //14 //21
"SDDL3_Rec2Raw", //15 //22 23
"SDDL4_Rec2Raw", //16 //22 23
"SDDL3_dedx", //17 //24 25
"SDDL4_dedx", //18 24 25
"SDDRecPointCheck"//19 26
};//tot 20
const TString RecPointsnonexpertname[]={"SDDGlobalCoordDistribL3PHIZ",//0 7
"SDDGlobalCoordDistribL4PHIZ",//1 7
"SDDGlobalCoordDistribRZ",//2 8
"SDDGlobalCoordDistribYX",//3 8
"SDDLay3TotCh", //4 9 17
"SDDLay4TotCh", //5 9 17
"SDDdrifttime_Layer3", //6 10 18
"SDDdrifttime_Layer4"//7 10 18
};//tot 8
const TString canvassavedname[]={"RawLayersNORM", //0
"Rawpatterns", //1
"RawLayers", //2
"CalibConfiguration", //3
"RawLayersRelOcc", //4
"RawLayersRelOccBoth", //5
"CheckRawdata", //6
"RecPointsphiz", //7
"RecPointsGlobalDistributions", //8
"RecPointsCharge", //9
"RecPointsDriftTime", //10
"RecPointsLocalDistribution", //11
"RecPointsLayers", //12
"RecPointsLayersNORM", //13
"RecPointsLayersPatterns", //14
"RecPointsPhiDistribution", //15
"RecPointsRDistribution", //16
"RecPointsChargeBothLayers", //17
"RecPointsDriftTimeBothLayers", //18
"RecPointsRelLayersOcc", //19
"RecPointsLayersRelOccBoth", //20
"RecPointsLayersRecRaw2D", //21
"RecPointsLayersRecRaw", //22
"RecPointsLayersRecRawBoth", //23
"RecPointsLayersdedx", //24
"RecPointsLayersdedxBoth", //25
"CheckRecPoints" //26
};//tot 27
Char_t trendfile[25]="FileQAtrend.root";
TString psfilestart=Form("SDDQAPlot.ps[");
TString psfile =Form("SDDQAPlot.ps");
TString psfileend =Form("SDDQAPlot.ps]");
gStyle->SetTitleFont(62,"");
gStyle->SetTitleFontSize(0.025);
gStyle->SetPalette(1);
// gStyle->SetOptStat("ne");
if(gSystem->Exec(Form("ls %s >/dev/null 2>&1", fileName))!=0)
{
printf(" No file created --- Exit");
return;
}
TFile mergedfile(fileName);
if(mergedfile.GetKey("ITS")==0x0){
printf("In this run ITS QA has not been executed.\n\nExit macro \n\n");
return;
}
Char_t filepath[100];
Char_t namecanvas[50];
TDirectory *directory=NULL;
TDirectory *directory2=NULL;
Char_t histoname[200];
TH2D *historaw=NULL;
TH1F *historaw2=NULL;
TH2F *histodraw=NULL;
// TH1F *histo2save[9];
//for(Int_t i=0;i<9;i++){histo2save[i]=NULL;}
TF1 *f1 = new TF1("f1","1",-100000,100000); // this is the function defined to multiply histograms
Char_t layer[10];
float fmax=0.;
float fmaxold=0.;
float fmaxmargin;
Int_t cannum=0;
//Int_t histo2savenumber[9]={4,5,6,7,7,8,3,5,1};
TCanvas *canvas1 = new TCanvas("canvas1","SDD QA Plot",1000,600);
// TLegend *legend=new TLegend(0.83,0.8,0.97,0.7);
TLegend *legend=new TLegend(0.81,0.895,0.95,0.995);
Bool_t updatecanvas[2];
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
// -------------This part is to read the number of chunks that were merged
Float_t fChnknmbr=0.;
FILE * pChunkNumber = fopen ("ChunkNumber.txt","r");
if(pChunkNumber){
Int_t rv=fscanf(pChunkNumber, "%f", &fChnknmbr);
fclose (pChunkNumber);
}
else fChnknmbr=1.;
gStyle->SetPalette(1);
float fCNinv=1./fChnknmbr;
// printf("\n====================>%f\n\n", fCNinv);
if(!mergedfile.IsOpen()){return;}else{printf("file is open\n");}
for(Int_t ispecie=0;ispecie<AliRecoParam::kNSpecies;ispecie++){
//__________________________________________________________________
//raw data
sprintf(filepath,"ITS/Raws/%s",AliRecoParam::GetEventSpecieName(ispecie));
printf("%s",filepath);
TString especie(filepath);
if(!especie.Contains(eventspecie)){printf("...Found and Skipped\n");continue;}
canvas1->Print(psfilestart.Data());
directory=(TDirectory*)mergedfile.Get(filepath);
if(!directory){printf("...Not Found\n ");}//faccio l'istogramma
else{
printf("...Found: The histograms of this EventSpecie will be displayed\n");
printf("0\n");
canvas1->Clear();
canvas1->Divide(2,1);
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
for(Int_t iraws=2;iraws<4;iraws++){//non expert raws
canvas1->cd(iraws-1);
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[iraws].Data());
printf("histo name %s ",Rawsexpertname[iraws].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
if(historaw){printf("...Found\n"); historaw->Multiply(f1,fCNinv); historaw->DrawCopy("colz");}
else
{
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[iraws].Data());
printf("...Not found\nSecond try for histo name %s ",Rawsexpertname[iraws].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
if(historaw){printf("...Found\n"); historaw->Multiply(f1,fCNinv);historaw->DrawCopy("colz");}
else
{
updatecanvas[iraws-2]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");
}
}
historaw=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas1; canvas1=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(1,2);
printf("1\n");
cannum++;
for(Int_t inraws=0;inraws<2;inraws++){//non expert raws
canvas1->cd(inraws+1);
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[inraws].Data());
printf("histo name %s ",Rawsexpertname[inraws].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
// -----------------This is the part were I divide between the number of chunks to normalize the histogram----------------
if(historaw2)
{
printf("...Found\n");
if (inraws==1)
{
historaw2->Multiply(f1,fCNinv);
}
historaw2->DrawCopy();
}
else {updatecanvas[inraws]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas2; canvas2=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
printf("2\n");
cannum++;
canvas1->Clear();
canvas1->Divide(2,1);
for(Int_t iraws=0;iraws<2;iraws++){//non expert raws
canvas1->cd(iraws+1);
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsnonexpertname[iraws].Data());
printf("histo name %s",Rawsnonexpertname[iraws].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
if(historaw){
printf("...Found\n");
// historaw->SetTitleSize(0.02);
historaw->DrawCopy("colz");
}
else
{
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsnonexpertname[iraws].Data());
printf("...Not Found.\n Second Try for histo name %s",Rawsnonexpertname[iraws].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
if(historaw){
printf("...Found\n");
// historaw->SetTitleSize(0.02);
historaw->Multiply(f1,fCNinv);
historaw->DrawCopy("colz");
}
else
{
updatecanvas[iraws]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");
}
}
historaw=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas3; canvas3=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
printf("3\n");
cannum++;
canvas1->Clear();
canvas1->Divide(2,1);
for(Int_t iraws=2;iraws<4;iraws++){//non expert raws
canvas1->cd(iraws-1);
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsnonexpertname[iraws].Data());
printf("histo name %s",Rawsnonexpertname[iraws].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
if(historaw){
printf("...Found\n");
// historaw->SetTitleSize(0.02);
// historaw->Multiply(f1,fCNinv);
historaw->DrawCopy("colz");
}
else
{
updatecanvas[iraws]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");
}
historaw=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//--------------------- new plots
printf("================= 4\n");
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
cannum++;
canvas1->Clear();
canvas1->Divide(2,1);
for(Int_t iraws=4;iraws<rawtot-1;iraws++){//non expert raws
canvas1->cd(iraws-3);
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[iraws].Data());
printf("histo name %s",Rawsexpertname[iraws].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2)
{
printf("...Found\n");
historaw2->DrawCopy();
}else
{
updatecanvas[iraws-4]=kFALSE;
printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");
}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas3; canvas3=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
legend->Clear();
printf("5\n");
cannum++;
for (Int_t i=4;i<rawtot-1;i++){
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){
fmax=historaw2->GetMaximum();
if (i==4){
fmaxold=fmax;}
if (i!=4){
if(fmaxold<fmax){
fmaxold=fmax;
}
}
}
historaw2=NULL;
}
fmaxmargin=1.1*fmaxold;
for(Int_t irrpp=4;irrpp<rawtot-1;irrpp++){//non expert raws
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[irrpp].Data());
printf("histo name %s",Rawsexpertname[irrpp].Data());
sprintf(layer, "layer %d",irrpp-1);
historaw2=(TH1F*)mergedfile.Get(histoname);
gStyle->SetOptStat(0);
if(historaw2)
{
printf("...Found\n");
historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
historaw2->SetTitle("Raw Relative Occupancy");
if (irrpp==4) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
if (irrpp!=4) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
legend->AddEntry(historaw2,layer,"l");
legend->Draw();
canvas1->Update();
}
else{updatecanvas[irrpp-4]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//---------------------------- summary plot ------------------------//
printf("6\n");
canvas1->Clear();
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[rawtot-1].Data());
printf("histo name %s",Rawsexpertname[rawtot-1].Data());
cannum++;
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2)
{
printf("...Found\n");
historaw2->SetStats(0);
historaw2->DrawCopy();
canvas1->Update();
}
else{updatecanvas[0]=kFALSE;updatecanvas[1]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
}//end if directory
//_______________________________________________________________________________
//rec point
sprintf(filepath,"ITS/RecPoints/%s",AliRecoParam::GetEventSpecieName(ispecie));
printf("%s\n",filepath);
directory2=(TDirectory*)mergedfile.Get(filepath);
if(directory2){
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("7\n");
cannum++;
for(Int_t irp=0;irp<2;irp++){//non expert rec point
canvas1->cd(irp+1);
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irp].Data());
printf("histo name %s ",RecPointsnonexpertname[irp].Data());
histodraw=(TH2F*)mergedfile.Get(histoname);
if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[irp]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
histodraw=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
// else{delete canvas4; canvas4=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("8\n");
cannum++;
for(Int_t irrp=2;irrp<4;irrp++){//non expert raws
canvas1->cd(irrp-1);
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrp].Data());
printf("histo name %s",RecPointsnonexpertname[irrp].Data() );
histodraw=(TH2F*)mergedfile.Get(histoname);
if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[irrp-2]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
histodraw=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas5; canvas5=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(1,2);
printf("9\n");
cannum++;
for(Int_t irrpp=4;irrpp<6;irrpp++){//non expert raws
canvas1->cd(irrpp-3);
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[irrpp-4]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas6; canvas6=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(1,2);
printf("10\n");
cannum++;
for(Int_t irrpp=6;irrpp<8;irrpp++){//non expert raws
canvas1->cd(irrpp-5);
// printf("histo name\n");
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[irrpp-6]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas7; canvas7=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
printf("11\n");
cannum++;
//canvas1->Divide(2,1);
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[rlocaldistro].Data());
printf("histo name %s",RecPointsexpertname[rlocaldistro].Data());
histodraw=(TH2F*)mergedfile.Get(histoname);
if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[0]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
histodraw=NULL;
if(updatecanvas[0]==kTRUE){
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas8; canvas8=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("12\n");
cannum++;
for(Int_t i=1;i<3;i++)
{
canvas1->cd(i);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
histodraw=(TH2F*)mergedfile.Get(histoname);
if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[i-1]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
histodraw=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//else{delete canvas9; canvas9=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("13\n");
cannum++;
for(Int_t i=3;i<5;i++)
{
canvas1->cd(i-2);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
histodraw=(TH2F*)mergedfile.Get(histoname);
printf("histo name %s",RecPointsexpertname[i].Data());
if(histodraw){
printf("...Found\n");histodraw->Multiply(f1,fCNinv);
histodraw->DrawCopy("colz");
}else{updatecanvas[i-3]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
histodraw=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(1,2);
printf("14\n");
cannum++;
for(Int_t i=5;i<7;i++)
{
canvas1->cd(i-4);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2)
{
printf("...Found\n");
if (i==4){historaw2->Multiply(f1,fCNinv);}
historaw2->DrawCopy();
}
else{updatecanvas[i-5]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("15\n");
cannum++;
for(Int_t i=7;i<9;i++)
{
canvas1->cd(i-6);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-7]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("16\n");
cannum++;
for(Int_t i=9;i<11;i++)
{
canvas1->cd(i-8);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-9]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//superimposed
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
legend->Clear();
printf("17\n");
cannum++;
for (Int_t i=4;i<6;i++){
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){
fmax=historaw2->GetMaximum();
if (i==4){
fmaxold=fmax;}
if (i!=4){
if(fmaxold<fmax){
fmaxold=fmax;
}
}
}
historaw2=NULL;
}
fmaxmargin=1.1*fmaxold;
for(Int_t irrpp=4;irrpp<6;irrpp++){//non expert raws
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
sprintf(layer, "layer %d",irrpp-1);
historaw2=(TH1F*)mergedfile.Get(histoname);
gStyle->SetOptStat(0);
if(historaw2)
{
printf("...Found\n");
historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
historaw2->SetTitle("Charge");
historaw2->SetFillColor(0);
if (irrpp==4) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
if (irrpp!=4) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
legend->AddEntry(historaw2,layer,"l");
legend->Draw();
canvas1->Update();
}
else{updatecanvas[irrpp-4]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
legend->Clear();
printf("18\n");
cannum++;
for (Int_t i=6;i<8;i++){
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){
fmax=historaw2->GetMaximum();
if (i==6){
fmaxold=fmax;}
if (i!=6){
if(fmaxold<fmax){
fmaxold=fmax;}}}
fmaxmargin=1.1*fmaxold;
}
for(Int_t irrpp=6;irrpp<8;irrpp++){//non expert raws
sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
sprintf(layer, "layer %d",irrpp-3);
historaw2=(TH1F*)mergedfile.Get(histoname);
gStyle->SetOptStat(0);
if(historaw2)
{
printf("...Found\n");
historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
historaw2->SetTitle("Drift Time");
historaw2->SetFillColor(0);
if (irrpp==6) {historaw2->SetStats(0); historaw2->SetLineColor(2);historaw2->DrawCopy();}
if (irrpp!=6) {historaw2->SetStats(0); historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
legend->AddEntry(historaw2,layer,"l");
legend->Draw();
canvas1->Update();
}
else{updatecanvas[irrpp-6]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
printf("%s\n%s\n",historaw2->GetName(),historaw2->GetTitle());
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//------------------------------------------- new plot
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("19\n");
cannum++;
for(Int_t i=11;i<13;i++)
{
canvas1->cd(i-10);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-11]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//------------------------------------
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
legend->Clear();
printf("20\n");
cannum++;
for (Int_t i=11;i<13;i++){
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){
fmax=historaw2->GetMaximum();
if (i==4){
fmaxold=fmax;}
if (i!=4){
if(fmaxold<fmax){
fmaxold=fmax;
}
}
}
historaw2=NULL;
}
fmaxmargin=1.1*fmaxold;
for(Int_t irrpp=11;irrpp<13;irrpp++){//non expert raws
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[irrpp].Data());
printf("histo name %s",RecPointsexpertname[irrpp].Data());
sprintf(layer, "layer %d",irrpp-8);
historaw2=(TH1F*)mergedfile.Get(histoname);
gStyle->SetOptStat(0);
if(historaw2)
{
printf("...Found\n");
historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
historaw2->SetTitle("RecPoint Relative Occupancy");
if (irrpp==11) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
if (irrpp!=11) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
legend->AddEntry(historaw2,layer,"l");
legend->Draw();
canvas1->Update();
}
else{updatecanvas[irrpp-11]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("21\n");
cannum++;
for(Int_t i=13;i<15;i++)
{
canvas1->cd(i-12);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
if(historaw){printf("...Found\n"); historaw->DrawCopy("colz");}else{updatecanvas[i-13]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//-----------------
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("22\n");
cannum++;
for(Int_t i=15;i<17;i++)
{
canvas1->cd(i-14);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-15]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//--------------------------------------------
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
legend->Clear();
printf("23\n");
cannum++;
for (Int_t i=15;i<17;i++){
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){
fmax=historaw2->GetMaximum();
if (i==15){
fmaxold=fmax;}
if (i!=15){
if(fmaxold<fmax){
fmaxold=fmax;
}
}
}
historaw2=NULL;
}
fmaxmargin=1.1*fmaxold;
for(Int_t irrpp=15;irrpp<17;irrpp++){//non expert raws
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[irrpp].Data());
printf("histo name %s",RecPointsexpertname[irrpp].Data());
sprintf(layer, "layer %d",irrpp-12);
historaw2=(TH1F*)mergedfile.Get(histoname);
gStyle->SetOptStat(0);
if(historaw2)
{
printf("...Found\n");
historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
historaw2->SetTitle("Rec2Raw Ratio");
if (irrpp==15) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
if (irrpp!=15) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
legend->AddEntry(historaw2,layer,"l");
legend->Draw();
canvas1->Update();
}
else{updatecanvas[irrpp-15]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//--------------------------------------
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
canvas1->Divide(2,1);
printf("24\n");
cannum++;
for(Int_t i=17;i<19;i++)
{
canvas1->cd(i-16);
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-17]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
//--------------------------------------------
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
legend->Clear();
printf("25\n");
cannum++;
for (Int_t i=17;i<19;i++){
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2){
fmax=historaw2->GetMaximum();
if (i==17){
fmaxold=fmax;}
if (i!=17){
if(fmaxold<fmax){fmaxold=fmax;}
}
}
historaw2=NULL;
}
fmaxmargin=1.1*fmaxold;
for(Int_t irrpp=17;irrpp<19;irrpp++){//non expert raws
sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[irrpp].Data());
printf("histo name %s",RecPointsexpertname[irrpp].Data());
sprintf(layer, "layer %d",irrpp-14);
historaw2=(TH1F*)mergedfile.Get(histoname);
gStyle->SetOptStat(0);
if(historaw2)
{
printf("...Found\n");
historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
historaw2->SetTitle("RecPoint dEdx");
if (irrpp==17) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
if (irrpp!=17) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
legend->AddEntry(historaw2,layer,"l");
legend->Draw();
canvas1->Update();
}
else{updatecanvas[irrpp-15]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
}//end for
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
//------------------------summary plot ---------------//
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
canvas1->Clear();
printf("26\n");
cannum++;
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[rtot-1].Data());
printf("histo name %s",RecPointsexpertname[rtot-1].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
if(historaw2)
{
printf("...Found\n");
historaw2->DrawCopy();
}
else{updatecanvas[0]=kFALSE;updatecanvas[1]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
if(kDoEps){
sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
}
canvas1->Print(psfile.Data());
}
}
canvas1->Print(psfileend.Data());
TFile file2savefortrend(trendfile,"recreate");
file2savefortrend.cd();
for(Int_t isave=0;isave<nonrawtot;isave++)
{
sprintf(histoname,"ITS/Raws/%s/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),Rawsnonexpertname[isave].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
historaw->Write();
printf("Saved %s\n",histoname);
historaw=NULL;
}
for(Int_t isave1=0;isave1<rawtot;isave1++)
{
sprintf(histoname,"ITS/Raws/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[isave1].Data());
if(isave1==2||isave1==3)
{
historaw=(TH2D*)mergedfile.Get(histoname);
historaw->Write();
historaw=NULL;
}
else
{
historaw2=(TH1F*)mergedfile.Get(histoname);
historaw2->Write();
historaw2=NULL;
}
printf("Saved %s\n",histoname);
}
for(Int_t isave2=0;isave2<nrtot;isave2++)
{
sprintf(histoname,"ITS/RecPoints/%s/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[isave2].Data());
if(isave2<4||isave2==17||isave2==18)
{
historaw=(TH2D*)mergedfile.Get(histoname);
historaw->Write();
historaw=NULL;
}
else
{
historaw2=(TH1F*)mergedfile.Get(histoname);
historaw2->Write();
historaw2=NULL;
}
printf("Saved %s\n",histoname);
}
for(Int_t isave3=0;isave3<rtot;isave3++)
{
sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[isave3].Data());
//printf(histoname,"ITS/RecPoints/%s/Expert/%s_%s \n",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[isave3].Data());
if(isave3<3||isave3==11||isave3==12||isave3==17||isave3==18)
{
histodraw=(TH2F*)mergedfile.Get(histoname);
histodraw->Write();
histodraw=NULL;
}
else
{
historaw2=(TH1F*)mergedfile.Get(histoname);
historaw2->Write();
historaw2=NULL;
}
printf("Saved %s\n",histoname);
}
//for(Int_t iss=0;iss<9;iss++){printf("Saved %d\n",iss); histo2save[iss]->Write();}
file2savefortrend.Close();
//else{delete canvas13; canvas13=NULL;}
//directory2=NULL;
}//end directory
}//end for
delete directory;
directory=NULL;
delete directory2;
directory2=NULL;
}//end macro