[u/mrichter/AliRoot.git] / FASTSIM / AliFastMuonTriggerEff.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/*
The trigger parametrization is computed for background levels 0., 0.5 and 1.
In order to set a background level different from 0 it is necessary to 
explicitly force it with:
ForceBkgLevel(BkgLevel).

For intermediate background levels, the trigger response is linearly 
interpolated between these values.

There is increased granularity in the pT region below 3 GeV. Although
it does not seem to be necessary it is also possible to interpolate
between pT bins using SetInt().

$Log$
Revision 1.10  2006/04/19 21:58:28  morsch
Coding rule violations corrected

Revision 1.9  2005/09/24 03:51:43  hristov
Removing extra semicolon (Fedora Core 4)

Revision 1.8  2005/09/22 11:31:42  morsch
Completely revised version. (P. Cortese)

*/

#include "AliFastMuonTriggerEff.h"
#include "TROOT.h"
#include "TFile.h"
#include "stdlib.h"
#include "TH3.h"
#include "TObjString.h"

#define PLIN printf("%s: %d: ",__FILE__,__LINE__)

ClassImp(AliFastMuonTriggerEff)

AliFastMuonTriggerEff::AliFastMuonTriggerEff():
    AliFastResponse("Efficiency", "Muon Trigger Efficiency")
{
//
// Default constructor
//
    SetCut(kLow);
    fZones=0;
    SetBkgLevel(0.);
    UnsetInt();
}

void AliFastMuonTriggerEff::SetCut(Int_t cut) 
{  
  //
  // Set the pt cut
  if(cut==kLow){
    printf("Selecting Low Pt cut\n");
  }else if(cut==kHigh){
    printf("Selecting High Pt cut\n");
  }else if(cut==kAny){
    printf("Selecting Lowest Pt cut\n");
  }else{
    printf("Don't know cut %d! Selecting Low Pt cut\n",cut);
    cut=kLow;
  }
  fCut = cut;
  
}

Int_t AliFastMuonTriggerEff::SetBkgLevel(Float_t Bkg)
{
  //
  // Set the background level
  //
    if((Bkg!=0.)) {
       printf("%s: Warning: requested Bkg: %f\n",
       __FILE__,Bkg);
       fBkg=0.;
       printf("A consistent treatement of the trigger probability\n");
       printf("within the framework of the fast simulation requires\n");
       printf("requires background 0\n");
       printf("%s: fBkg: set to %f\n",
       __FILE__,fBkg);
    } else {
      fBkg=Bkg;
    }
    if(fZones!=0.) {
	Init();
    }
    return 0;
}

Int_t AliFastMuonTriggerEff::ForceBkgLevel(Float_t Bkg)
{
  //
  // Check and enforce consistency of the background level 
  // 
    if((Bkg!=0.)) {
       printf("%s: Warning: requested Bkg: %f\n",
       __FILE__,Bkg);
       printf("A consistent treatement of the trigger probability\n");
       printf("within the framework of the fast simulation\n");
       printf("requires background 0");
       printf("%s: Continue with fBkg: %f\n",
       __FILE__,Bkg);
    }
    fBkg=Bkg;
    if(fZones!=0.) {
	Init();
    }
    return 0;
}

Int_t AliFastMuonTriggerEff::LoadTables(Char_t *namet=""){
  //
  // Load the trigger tables
  //
    Char_t hNameA[100],hNameL[100],hNameH[100];
    sprintf(hNameA,"hEffAPt%s",namet);
    sprintf(hNameL,"hEffLPt%s",namet);
    sprintf(hNameH,"hEffHPt%s",namet);
    fhEffAPt = (TH3F*)gDirectory->Get(hNameA);
    fhEffLPt = (TH3F*)gDirectory->Get(hNameL);
    fhEffHPt = (TH3F*)gDirectory->Get(hNameH);
    if(!fhEffAPt){
      PLIN; printf("%s: histogram %s not found\n",__FILE__,hNameA);
      return -1;
    }
    if(!fhEffLPt){
      PLIN; printf("%s: histogram %s not found\n",__FILE__,hNameL);
      return -2;
    }
    if(!fhEffHPt){
      PLIN; printf("%s: histogram %s not found\n",__FILE__,hNameH);
      return -3;
    }
    return 0;
}

void AliFastMuonTriggerEff::Init()
{
//
//  Initialization
//
    fZones=0;
    Char_t file[100]="$(ALICE_ROOT)/FASTSIM/data/MUONtriggerLUT_V2.4nvdn.root";
    printf("Initializing %s / %s\n", fName.Data(), fTitle.Data());
    printf("using data from file: %s\n",file);
    printf("AliFastMuonTriggerEff: Initialization with background level: %f\n",fBkg);
    TFile *f = new TFile(file);
    if(f->IsZombie()) {
        PLIN; printf("Cannot open file: %s\n",file);
        return;
    }
    f->ls();
    Int_t intb=0;
    Char_t namet[10];
    if(TMath::Abs(fBkg)<0.00001){
      sprintf(namet,"00");
    }else if(TMath::Abs(fBkg-0.5)<0.00001){
      sprintf(namet,"05");
    }else if(TMath::Abs(fBkg-1.0)<0.00001){
      sprintf(namet,"10");
    }else{
      PLIN; printf("A table for Bkg level: %f does not exists\n",fBkg);
      intb=1;
    }
    if(intb){ // Interpolation between background levels
      PLIN; printf("Interpolating Bkg level: %f\n",fBkg);
      TH3F* ha1,*hl1,*hh1,*ha2,*hl2,*hh2,*ha0,*hl0,*hh0;
      Char_t name1[10],name2[10]; Float_t b1,b2;
      if(fBkg>0&&fBkg<0.5){
        sprintf(name1,"00");
        sprintf(name2,"05");
	b1=0.;
	b2=0.5;
      }else if(fBkg>0.5){
        sprintf(name1,"05");
        sprintf(name2,"10");
	b1=0.5;
	b2=1.0;
	if(fBkg>1.0){
	  for(Int_t i=0; i<10;i++){
	    PLIN; printf("WARNING!!!! You are extrapolating above background 1.0\n");
	  }
	}
      }else{
        PLIN; printf("Bkg level: %f is not supported\n",fBkg);
        return;
      }
      if(LoadTables(name1)){
        PLIN; printf("Error in loading trigger tables\n");
	return;
      }
      PLIN; printf("We use tables for %f and %f to interpolate %f Bkg level\n",b1,b2,fBkg);
      ha0=(TH3F*)fhEffAPt->Clone("hEffAPtXX"); ha0->Reset();
      hl0=(TH3F*)fhEffLPt->Clone("hEffLPtXX"); hl0->Reset();
      hh0=(TH3F*)fhEffHPt->Clone("hEffHPtXX"); hh0->Reset();
      ha1=fhEffAPt;
      hl1=fhEffLPt;
      hh1=fhEffHPt;
      if(LoadTables(name2)){
        PLIN; printf("Error in loading trigger tables\n");
	return;
      }
      ha2=fhEffAPt;
      hl2=fhEffLPt;
      hh2=fhEffHPt;
      fhEffAPt=ha0;
      fhEffLPt=hl0;
      fhEffHPt=hh0;
      Int_t nnx=ha0->GetNbinsX()+1;
      Int_t nny=ha0->GetNbinsY()+1;
      Int_t nnz=ha0->GetNbinsZ()+1;
      for(Int_t ix=0; ix<=nnx; ix++){
	for(Int_t iy=0; iy<=nny; iy++){
	  for(Int_t iz=0; iz<=nnz; iz++){
	    Double_t y1,y2; Float_t cont;
	    y1=ha1->GetBinContent(ix,iy,iz); y2=ha2->GetBinContent(ix,iy,iz);
	    cont=Float_t(y1+(y2-y1)/(b2-b1)*(fBkg-b1)); if(cont>1)cont=1; if(cont<0)cont=0;
	    fhEffAPt->SetBinContent(ix,iy,iz,cont);
	    y1=hl1->GetBinContent(ix,iy,iz); y2=hl2->GetBinContent(ix,iy,iz);
	    cont=Float_t(y1+(y2-y1)/(b2-b1)*(fBkg-b1)); if(cont>1)cont=1; if(cont<0)cont=0;
	    fhEffLPt->SetBinContent(ix,iy,iz,cont);
	    y1=hh1->GetBinContent(ix,iy,iz); y2=hh2->GetBinContent(ix,iy,iz);
	    cont=Float_t(y1+(y2-y1)/(b2-b1)*(fBkg-b1)); if(cont>1)cont=1; if(cont<0)cont=0;
	    fhEffHPt->SetBinContent(ix,iy,iz,cont);
	  }
	}
      }
    }else{ // Use tables computed for selected backgound levels
      printf("Loading tables for background level: %f\n",fBkg);
      if(LoadTables(namet)){
        PLIN; printf("Error in loading trigger tables\n");
	return;
      }
    }
    fhEffAPt->SetDirectory(0);
    fhEffLPt->SetDirectory(0);
    fhEffHPt->SetDirectory(0);
    fhLX=fhEffLPt->GetXaxis();
    fhLY=fhEffLPt->GetYaxis();
    fhLZ=fhEffLPt->GetZaxis();
//
//
    if(f->Get("Description"))
    {
      fDescription=((TObjString*)f->Get("Description"))->GetString();
      printf("%s\n",fDescription.Data());
    }

    fThetaMin = fhEffLPt->GetXaxis()->GetXmin();
    fThetaMax = fhEffLPt->GetXaxis()->GetXmax();
    fnthetab=fhEffLPt->GetNbinsX();
    fDtheta   = (fThetaMax-fThetaMin)/fnthetab;

    fPhiMin   = fhEffLPt->GetYaxis()->GetXmin();
    fPhiMax   = fhEffLPt->GetYaxis()->GetXmax();
    fnphib=fhEffLPt->GetNbinsY();
    fDphi     = (fPhiMax-fPhiMin)/fnphib;

    fPtMin=fhEffLPt->GetZaxis()->GetXmin();
    fPtMax=fhEffLPt->GetZaxis()->GetXmax();
    fnptb=fhEffLPt->GetNbinsZ();
    fDpt      = (fPtMax-fPtMin)/fnptb;

    printf("%4d bins in theta [%f:%f]\n",fnthetab,fThetaMin,fThetaMax);
    printf("%4d bins in phi [%f:%f]\n",fnphib,fPhiMin,fPhiMax);
    printf("%4d bins in pt [%f:%f]\n",fnptb,fPtMin,fPtMax);

    fZones=fnthetab*fnphib;

    f->Close();
    if(fInt==0) {
      printf("Interpolation of trigger efficiencies is off!\n");
    } else {
      printf("Interpolation of trigger efficiencies is on!\n");
    }
}

void AliFastMuonTriggerEff::Evaluate(Float_t charge, Float_t pt,Float_t theta,
              Float_t phi, Float_t& effLow, Float_t& effHigh, Float_t& effAny)
{
    //
    //  Trigger efficiency for pt, theta, phi (low, high and "any" cut)
    //
#ifdef MYTRIGDEBUG
    printf("Evaluate(ch=%2.0f, pt=%10.6f, theta=%7.2f, phi=%8.2f ...)\n",charge,pt,theta,phi);
#endif
    effLow=0.;
    effHigh=0.;
    effAny=0;
    if(fZones==0) {
        printf("Call to uninitialized object of class: AliFastMuonTriggerEff\n");
	return;
    }
    if(pt<0) {
        printf("Warning: pt: %f < 0. GeV/c\n",pt);
	return;	
    }

    Int_t iPt   = fhLZ->FindBin((Double_t)pt);
    if(iPt>fnptb)iPt=fnptb;
    Int_t iPhi  = Int_t((phi-fPhiMin)/fDphi);
    if(phi<fPhiMin)iPhi=iPhi-1;
    Int_t iTheta = fhLX->FindBin((Double_t)theta);
#ifdef MYTRIGDEBUG
    printf("Evaluate(ch=%2.0f, pt=%10.6f, theta=%7.2f, phi=%8.2f ...)\n",charge,pt,theta,phi);
    printf(" 0:%1d iPt iTheta iPhi: %d %d %d\n",fInt,iPt,iTheta,iPhi);
#endif
    iPhi=iPhi-2*fnphib*(iPhi/(2*fnphib));
#ifdef MYTRIGDEBUG
    printf(" 1:%1d iPhi converted to: %d for angle equivalence\n",fInt,iPhi);
#endif
    if(iPhi<0)iPhi=-iPhi-1;
    if(iPhi>(fnphib-1))iPhi=2*fnphib-1-iPhi;
#ifdef MYTRIGDEBUG
    printf(" 2:%1d iPhi converted to: %d for the symmetry of the spectrometer\n",fInt,iPhi);
#endif
    if(charge==1.){
    } else if(charge==-1.) {
    iPhi=fnphib-1-iPhi;
#ifdef MYTRIGDEBUG
    printf(" 3:%1d iPhi converted to: %d for the charge symmetry\n",fInt,iPhi);
#endif
    } else {
        printf("Warning: not understand charge: %f\n",charge);
        return;
    }
    if(iTheta<=0||iTheta>fnthetab) {
        printf("Warning: theta: %f outside acceptance\n",theta);
        return;
    }
    if(iPt<0) {
        printf("Warning: what do you mean with pt: %f <0?\n",pt);
	return;
    }
    iPhi++;
#ifdef MYTRIGDEBUG
    printf(" 4:%1d Getting: iTheta, iPhi, iPt: %d %d %d\n",
              fInt,iTheta,iPhi,iPt);
#endif
    effLow =fhEffLPt->GetBinContent(iTheta,iPhi,iPt);
    effHigh=fhEffHPt->GetBinContent(iTheta,iPhi,iPt);
    effAny =fhEffAPt->GetBinContent(iTheta,iPhi,iPt);
#ifdef MYTRIGDEBUG
    printf(" 4:%1d Result: charge, iTheta, iPhi, iPt: %f %d %d %d effLow: %f, effHigh: %f, effAny: %f\n",
              fInt,charge,iTheta,iPhi,iPt,effLow,effHigh,effAny);
#endif
        
    if(fInt==1) {
      Float_t angl,angh,anga;
      Float_t effLowp,effHighp,effAnyp;
      Float_t ptc=(iPt+0.5)*fDpt;  // The center of current bin
      #ifdef MYTRIGDEBUG
        printf(" 5:1 The center of current bin iPt: %d is: %f\n",iPt,ptc);
      #endif
      if(iPt==fnptb) {
        #ifdef MYTRIGDEBUG
	printf(" 6:1 No more points above! No interpolation is needed!\n");
        #endif
	return;        
      }else if(ptc==pt){
        #ifdef MYTRIGDEBUG
	  printf(" 6:1 No interpolation is needed!\n");
        #endif
	return;
      }else if(ptc>pt){
	// Looking for previous point
        if(iPt>1) {
          effLowp =fhEffLPt->GetBinContent(iTheta,iPhi,iPt-1);
          effHighp=fhEffHPt->GetBinContent(iTheta,iPhi,iPt-1);
          effAnyp =fhEffAPt->GetBinContent(iTheta,iPhi,iPt-1);
          #ifdef MYTRIGDEBUG
	  printf(" 7:1 A simple look to previous point: %d: %f %f\n",iPt-1,effLowp,effHighp);
          #endif
	} else {
	  effLowp=0.;
	  effHighp=0.;
          effAnyp=0;
          #ifdef MYTRIGDEBUG
          printf(" 8:1 result is: %f %f %f\n",effLowp,effHighp,effAnyp);
	  #endif	  
	}
        angl=(effLow-effLowp)/fDpt;
        angh=(effHigh-effHighp)/fDpt;
        anga=(effAny-effAnyp)/fDpt;
      }else{
	// Looking for next point
	if(iPt<fnptb) {
          effLowp =fhEffLPt->GetBinContent(iTheta,iPhi,iPt+1);
          effHighp=fhEffHPt->GetBinContent(iTheta,iPhi,iPt+1);
          effAnyp =fhEffAPt->GetBinContent(iTheta,iPhi,iPt+1);
          #ifdef MYTRIGDEBUG
	  printf(" 7:1 A simple look to next point: %d: %f %f %f\n",iPt-1,effLowp,effHighp,effAnyp);
          #endif
	} else {
	  effLowp=effLow;
	  effHighp=effHigh;
	  effAnyp=effAny;
          #ifdef MYTRIGDEBUG
          printf(" 8:1 result is: pt: %f %f %f\n",effLowp,effHighp,effAnyp);
	  #endif	  
        }
        angl=(effLowp-effLow)/fDpt;
        angh=(effHighp-effHigh)/fDpt;
        anga=(effAnyp-effAny)/fDpt;
      }
      effLow=effLow+angl*(pt-ptc);
      effHigh=effHigh+angh*(pt-ptc);
      effAny=effAny+anga*(pt-ptc);
      #ifdef MYTRIGDEBUG
      printf(" 9:1 the interpolation coefficients are: %f %f %f\n",angl,angh,anga);
      #endif
  }
  #ifdef MYTRIGDEBUG
  printf("10:%1d effLow, effHigh=%f %f %f\n",fInt,effLow,effHigh,effAny);
  #endif
  return;
}


Float_t AliFastMuonTriggerEff::Evaluate(Float_t charge, Float_t pt,
                   Float_t theta, Float_t phi)
{
    //
    // Trigger efficiency for pt, theta, phi depending of fCut
    // 
    if(fZones==0) {
        printf("Call to uninitialized object of class: AliFastMuonTriggerEff\n");
	return 0.;
    }
    Float_t eff;
    Float_t effLow, effHigh, effAny;
    
    Evaluate(charge,pt,theta,phi,effLow,effHigh,effAny);
    if (fCut == kLow) 
	eff  = effLow;
    else if (fCut == kHigh)
	eff  = effHigh;
    else if (fCut == kAny)
        eff  = effAny;
    else
        eff  = 0;

    return eff;
}
Commit	Line	Data
6255180c	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
490243f3	16	/*
	17	The trigger parametrization is computed for background levels 0., 0.5 and 1.
	18	In order to set a background level different from 0 it is necessary to
	19	explicitly force it with:
	20	ForceBkgLevel(BkgLevel).
6255180c	21
490243f3	22	For intermediate background levels, the trigger response is linearly
	23	interpolated between these values.
	24
	25	There is increased granularity in the pT region below 3 GeV. Although
	26	it does not seem to be necessary it is also possible to interpolate
	27	between pT bins using SetInt().
6255180c	28
490243f3	29	$Log$
bc77d57e	30	Revision 1.10 2006/04/19 21:58:28 morsch
	31	Coding rule violations corrected
	32
2724a683	33	Revision 1.9 2005/09/24 03:51:43 hristov
	34	Removing extra semicolon (Fedora Core 4)
	35
ecda8707	36	Revision 1.8 2005/09/22 11:31:42 morsch
	37	Completely revised version. (P. Cortese)
	38
490243f3	39	*/
	40
	41	#include "AliFastMuonTriggerEff.h"
2724a683	42	#include "TROOT.h"
	43	#include "TFile.h"
	44	#include "stdlib.h"
	45	#include "TH3.h"
	46	#include "TObjString.h"
	47
490243f3	48	#define PLIN printf("%s: %d: ",__FILE__,__LINE__)
6255180c	49
ecda8707	50	ClassImp(AliFastMuonTriggerEff)
6255180c	51
	52	AliFastMuonTriggerEff::AliFastMuonTriggerEff():
	53	AliFastResponse("Efficiency", "Muon Trigger Efficiency")
	54	{
	55	//
	56	// Default constructor
	57	//
	58	SetCut(kLow);
	59	fZones=0;
490243f3	60	SetBkgLevel(0.);
490243f3	61	UnsetInt();
6255180c	62	}
6255180c	63
490243f3	64	void AliFastMuonTriggerEff::SetCut(Int_t cut)
490243f3	65	{
2724a683	66	//
2724a683	67	// Set the pt cut
490243f3	68	if(cut==kLow){
	69	printf("Selecting Low Pt cut\n");
	70	}else if(cut==kHigh){
	71	printf("Selecting High Pt cut\n");
	72	}else if(cut==kAny){
	73	printf("Selecting Lowest Pt cut\n");
	74	}else{
	75	printf("Don't know cut %d! Selecting Low Pt cut\n",cut);
	76	cut=kLow;
	77	}
	78	fCut = cut;
	79
	80	}
6255180c	81
490243f3	82	Int_t AliFastMuonTriggerEff::SetBkgLevel(Float_t Bkg)
6255180c	83	{
2724a683	84	//
	85	// Set the background level
	86	//
490243f3	87	if((Bkg!=0.)) {
	88	printf("%s: Warning: requested Bkg: %f\n",
	89	__FILE__,Bkg);
	90	fBkg=0.;
	91	printf("A consistent treatement of the trigger probability\n");
	92	printf("within the framework of the fast simulation requires\n");
	93	printf("requires background 0\n");
	94	printf("%s: fBkg: set to %f\n",
	95	__FILE__,fBkg);
	96	} else {
	97	fBkg=Bkg;
	98	}
	99	if(fZones!=0.) {
	100	Init();
	101	}
	102	return 0;
	103	}
6255180c	104
490243f3	105	Int_t AliFastMuonTriggerEff::ForceBkgLevel(Float_t Bkg)
490243f3	106	{
2724a683	107	//
	108	// Check and enforce consistency of the background level
	109	//
490243f3	110	if((Bkg!=0.)) {
	111	printf("%s: Warning: requested Bkg: %f\n",
	112	__FILE__,Bkg);
	113	printf("A consistent treatement of the trigger probability\n");
	114	printf("within the framework of the fast simulation\n");
	115	printf("requires background 0");
	116	printf("%s: Continue with fBkg: %f\n",
	117	__FILE__,Bkg);
	118	}
	119	fBkg=Bkg;
	120	if(fZones!=0.) {
	121	Init();
	122	}
	123	return 0;
	124	}
6255180c	125
490243f3	126	Int_t AliFastMuonTriggerEff::LoadTables(Char_t *namet=""){
2724a683	127	//
	128	// Load the trigger tables
	129	//
490243f3	130	Char_t hNameA[100],hNameL[100],hNameH[100];
	131	sprintf(hNameA,"hEffAPt%s",namet);
	132	sprintf(hNameL,"hEffLPt%s",namet);
	133	sprintf(hNameH,"hEffHPt%s",namet);
	134	fhEffAPt = (TH3F*)gDirectory->Get(hNameA);
	135	fhEffLPt = (TH3F*)gDirectory->Get(hNameL);
	136	fhEffHPt = (TH3F*)gDirectory->Get(hNameH);
	137	if(!fhEffAPt){
	138	PLIN; printf("%s: histogram %s not found\n",__FILE__,hNameA);
	139	return -1;
	140	}
	141	if(!fhEffLPt){
	142	PLIN; printf("%s: histogram %s not found\n",__FILE__,hNameL);
	143	return -2;
	144	}
	145	if(!fhEffHPt){
	146	PLIN; printf("%s: histogram %s not found\n",__FILE__,hNameH);
	147	return -3;
	148	}
	149	return 0;
6255180c	150	}
6255180c	151
490243f3	152	void AliFastMuonTriggerEff::Init()
6255180c	153	{
6255180c	154	//
490243f3	155	// Initialization
6255180c	156	//
490243f3	157	fZones=0;
	158	Char_t file[100]="$(ALICE_ROOT)/FASTSIM/data/MUONtriggerLUT_V2.4nvdn.root";
	159	printf("Initializing %s / %s\n", fName.Data(), fTitle.Data());
	160	printf("using data from file: %s\n",file);
	161	printf("AliFastMuonTriggerEff: Initialization with background level: %f\n",fBkg);
6255180c	162	TFile *f = new TFile(file);
6255180c	163	if(f->IsZombie()) {
490243f3	164	PLIN; printf("Cannot open file: %s\n",file);
6255180c	165	return;
	166	}
	167	f->ls();
490243f3	168	Int_t intb=0;
	169	Char_t namet[10];
	170	if(TMath::Abs(fBkg)<0.00001){
	171	sprintf(namet,"00");
	172	}else if(TMath::Abs(fBkg-0.5)<0.00001){
	173	sprintf(namet,"05");
	174	}else if(TMath::Abs(fBkg-1.0)<0.00001){
	175	sprintf(namet,"10");
	176	}else{
	177	PLIN; printf("A table for Bkg level: %f does not exists\n",fBkg);
	178	intb=1;
6255180c	179	}
490243f3	180	if(intb){ // Interpolation between background levels
	181	PLIN; printf("Interpolating Bkg level: %f\n",fBkg);
	182	TH3F* ha1,hl1,hh1,ha2,hl2,hh2,ha0,hl0,hh0;
	183	Char_t name1[10],name2[10]; Float_t b1,b2;
	184	if(fBkg>0&&fBkg<0.5){
	185	sprintf(name1,"00");
	186	sprintf(name2,"05");
	187	b1=0.;
	188	b2=0.5;
	189	}else if(fBkg>0.5){
	190	sprintf(name1,"05");
	191	sprintf(name2,"10");
	192	b1=0.5;
	193	b2=1.0;
	194	if(fBkg>1.0){
	195	for(Int_t i=0; i<10;i++){
	196	PLIN; printf("WARNING!!!! You are extrapolating above background 1.0\n");
	197	}
6255180c	198	}
490243f3	199	}else{
	200	PLIN; printf("Bkg level: %f is not supported\n",fBkg);
	201	return;
	202	}
	203	if(LoadTables(name1)){
	204	PLIN; printf("Error in loading trigger tables\n");
	205	return;
	206	}
	207	PLIN; printf("We use tables for %f and %f to interpolate %f Bkg level\n",b1,b2,fBkg);
	208	ha0=(TH3F*)fhEffAPt->Clone("hEffAPtXX"); ha0->Reset();
	209	hl0=(TH3F*)fhEffLPt->Clone("hEffLPtXX"); hl0->Reset();
	210	hh0=(TH3F*)fhEffHPt->Clone("hEffHPtXX"); hh0->Reset();
	211	ha1=fhEffAPt;
	212	hl1=fhEffLPt;
	213	hh1=fhEffHPt;
	214	if(LoadTables(name2)){
	215	PLIN; printf("Error in loading trigger tables\n");
	216	return;
	217	}
	218	ha2=fhEffAPt;
	219	hl2=fhEffLPt;
	220	hh2=fhEffHPt;
	221	fhEffAPt=ha0;
	222	fhEffLPt=hl0;
	223	fhEffHPt=hh0;
	224	Int_t nnx=ha0->GetNbinsX()+1;
	225	Int_t nny=ha0->GetNbinsY()+1;
	226	Int_t nnz=ha0->GetNbinsZ()+1;
	227	for(Int_t ix=0; ix<=nnx; ix++){
	228	for(Int_t iy=0; iy<=nny; iy++){
	229	for(Int_t iz=0; iz<=nnz; iz++){
	230	Double_t y1,y2; Float_t cont;
	231	y1=ha1->GetBinContent(ix,iy,iz); y2=ha2->GetBinContent(ix,iy,iz);
	232	cont=Float_t(y1+(y2-y1)/(b2-b1)*(fBkg-b1)); if(cont>1)cont=1; if(cont<0)cont=0;
	233	fhEffAPt->SetBinContent(ix,iy,iz,cont);
	234	y1=hl1->GetBinContent(ix,iy,iz); y2=hl2->GetBinContent(ix,iy,iz);
	235	cont=Float_t(y1+(y2-y1)/(b2-b1)*(fBkg-b1)); if(cont>1)cont=1; if(cont<0)cont=0;
	236	fhEffLPt->SetBinContent(ix,iy,iz,cont);
	237	y1=hh1->GetBinContent(ix,iy,iz); y2=hh2->GetBinContent(ix,iy,iz);
	238	cont=Float_t(y1+(y2-y1)/(b2-b1)*(fBkg-b1)); if(cont>1)cont=1; if(cont<0)cont=0;
	239	fhEffHPt->SetBinContent(ix,iy,iz,cont);
	240	}
6255180c	241	}
490243f3	242	}
	243	}else{ // Use tables computed for selected backgound levels
	244	printf("Loading tables for background level: %f\n",fBkg);
	245	if(LoadTables(namet)){
	246	PLIN; printf("Error in loading trigger tables\n");
	247	return;
	248	}
6255180c	249	}
490243f3	250	fhEffAPt->SetDirectory(0);
	251	fhEffLPt->SetDirectory(0);
	252	fhEffHPt->SetDirectory(0);
	253	fhLX=fhEffLPt->GetXaxis();
	254	fhLY=fhEffLPt->GetYaxis();
	255	fhLZ=fhEffLPt->GetZaxis();
	256	//
	257	//
	258	if(f->Get("Description"))
	259	{
	260	fDescription=((TObjString*)f->Get("Description"))->GetString();
	261	printf("%s\n",fDescription.Data());
6255180c	262	}
490243f3	263
	264	fThetaMin = fhEffLPt->GetXaxis()->GetXmin();
	265	fThetaMax = fhEffLPt->GetXaxis()->GetXmax();
	266	fnthetab=fhEffLPt->GetNbinsX();
	267	fDtheta = (fThetaMax-fThetaMin)/fnthetab;
	268
	269	fPhiMin = fhEffLPt->GetYaxis()->GetXmin();
	270	fPhiMax = fhEffLPt->GetYaxis()->GetXmax();
	271	fnphib=fhEffLPt->GetNbinsY();
	272	fDphi = (fPhiMax-fPhiMin)/fnphib;
	273
	274	fPtMin=fhEffLPt->GetZaxis()->GetXmin();
	275	fPtMax=fhEffLPt->GetZaxis()->GetXmax();
	276	fnptb=fhEffLPt->GetNbinsZ();
	277	fDpt = (fPtMax-fPtMin)/fnptb;
	278
	279	printf("%4d bins in theta [%f:%f]\n",fnthetab,fThetaMin,fThetaMax);
	280	printf("%4d bins in phi [%f:%f]\n",fnphib,fPhiMin,fPhiMax);
	281	printf("%4d bins in pt [%f:%f]\n",fnptb,fPtMin,fPtMax);
	282
	283	fZones=fnthetab*fnphib;
	284
6255180c	285	f->Close();
490243f3	286	if(fInt==0) {
	287	printf("Interpolation of trigger efficiencies is off!\n");
	288	} else {
	289	printf("Interpolation of trigger efficiencies is on!\n");
	290	}
6255180c	291	}
6255180c	292
490243f3	293	void AliFastMuonTriggerEff::Evaluate(Float_t charge, Float_t pt,Float_t theta,
490243f3	294	Float_t phi, Float_t& effLow, Float_t& effHigh, Float_t& effAny)
6255180c	295	{
490243f3	296	//
	297	// Trigger efficiency for pt, theta, phi (low, high and "any" cut)
	298	//
	299	#ifdef MYTRIGDEBUG
	300	printf("Evaluate(ch=%2.0f, pt=%10.6f, theta=%7.2f, phi=%8.2f ...)\n",charge,pt,theta,phi);
	301	#endif
6255180c	302	effLow=0.;
6255180c	303	effHigh=0.;
490243f3	304	effAny=0;
6255180c	305	if(fZones==0) {
	306	printf("Call to uninitialized object of class: AliFastMuonTriggerEff\n");
	307	return;
	308	}
	309	if(pt<0) {
	310	printf("Warning: pt: %f < 0. GeV/c\n",pt);
	311	return;
	312	}
490243f3	313
	314	Int_t iPt = fhLZ->FindBin((Double_t)pt);
	315	if(iPt>fnptb)iPt=fnptb;
6255180c	316	Int_t iPhi = Int_t((phi-fPhiMin)/fDphi);
6255180c	317	if(phi<fPhiMin)iPhi=iPhi-1;
490243f3	318	Int_t iTheta = fhLX->FindBin((Double_t)theta);
6255180c	319	#ifdef MYTRIGDEBUG
490243f3	320	printf("Evaluate(ch=%2.0f, pt=%10.6f, theta=%7.2f, phi=%8.2f ...)\n",charge,pt,theta,phi);
490243f3	321	printf(" 0:%1d iPt iTheta iPhi: %d %d %d\n",fInt,iPt,iTheta,iPhi);
6255180c	322	#endif
490243f3	323	iPhi=iPhi-2fnphib(iPhi/(2*fnphib));
6255180c	324	#ifdef MYTRIGDEBUG
490243f3	325	printf(" 1:%1d iPhi converted to: %d for angle equivalence\n",fInt,iPhi);
6255180c	326	#endif
6255180c	327	if(iPhi<0)iPhi=-iPhi-1;
490243f3	328	if(iPhi>(fnphib-1))iPhi=2*fnphib-1-iPhi;
6255180c	329	#ifdef MYTRIGDEBUG
490243f3	330	printf(" 2:%1d iPhi converted to: %d for the symmetry of the spectrometer\n",fInt,iPhi);
6255180c	331	#endif
	332	if(charge==1.){
	333	} else if(charge==-1.) {
490243f3	334	iPhi=fnphib-1-iPhi;
6255180c	335	#ifdef MYTRIGDEBUG
490243f3	336	printf(" 3:%1d iPhi converted to: %d for the charge symmetry\n",fInt,iPhi);
6255180c	337	#endif
	338	} else {
	339	printf("Warning: not understand charge: %f\n",charge);
	340	return;
	341	}
490243f3	342	if(iTheta<=0\|\|iTheta>fnthetab) {
6255180c	343	printf("Warning: theta: %f outside acceptance\n",theta);
	344	return;
	345	}
	346	if(iPt<0) {
	347	printf("Warning: what do you mean with pt: %f <0?\n",pt);
	348	return;
	349	}
490243f3	350	iPhi++;
6255180c	351	#ifdef MYTRIGDEBUG
490243f3	352	printf(" 4:%1d Getting: iTheta, iPhi, iPt: %d %d %d\n",
490243f3	353	fInt,iTheta,iPhi,iPt);
6255180c	354	#endif
490243f3	355	effLow =fhEffLPt->GetBinContent(iTheta,iPhi,iPt);
	356	effHigh=fhEffHPt->GetBinContent(iTheta,iPhi,iPt);
	357	effAny =fhEffAPt->GetBinContent(iTheta,iPhi,iPt);
6255180c	358	#ifdef MYTRIGDEBUG
490243f3	359	printf(" 4:%1d Result: charge, iTheta, iPhi, iPt: %f %d %d %d effLow: %f, effHigh: %f, effAny: %f\n",
490243f3	360	fInt,charge,iTheta,iPhi,iPt,effLow,effHigh,effAny);
6255180c	361	#endif
490243f3	362
	363	if(fInt==1) {
	364	Float_t angl,angh,anga;
	365	Float_t effLowp,effHighp,effAnyp;
	366	Float_t ptc=(iPt+0.5)*fDpt; // The center of current bin
	367	#ifdef MYTRIGDEBUG
	368	printf(" 5:1 The center of current bin iPt: %d is: %f\n",iPt,ptc);
	369	#endif
	370	if(iPt==fnptb) {
	371	#ifdef MYTRIGDEBUG
	372	printf(" 6:1 No more points above! No interpolation is needed!\n");
	373	#endif
	374	return;
	375	}else if(ptc==pt){
	376	#ifdef MYTRIGDEBUG
	377	printf(" 6:1 No interpolation is needed!\n");
	378	#endif
	379	return;
	380	}else if(ptc>pt){
	381	// Looking for previous point
	382	if(iPt>1) {
	383	effLowp =fhEffLPt->GetBinContent(iTheta,iPhi,iPt-1);
	384	effHighp=fhEffHPt->GetBinContent(iTheta,iPhi,iPt-1);
	385	effAnyp =fhEffAPt->GetBinContent(iTheta,iPhi,iPt-1);
	386	#ifdef MYTRIGDEBUG
	387	printf(" 7:1 A simple look to previous point: %d: %f %f\n",iPt-1,effLowp,effHighp);
	388	#endif
	389	} else {
	390	effLowp=0.;
	391	effHighp=0.;
	392	effAnyp=0;
	393	#ifdef MYTRIGDEBUG
	394	printf(" 8:1 result is: %f %f %f\n",effLowp,effHighp,effAnyp);
	395	#endif
	396	}
	397	angl=(effLow-effLowp)/fDpt;
	398	angh=(effHigh-effHighp)/fDpt;
	399	anga=(effAny-effAnyp)/fDpt;
	400	}else{
	401	// Looking for next point
	402	if(iPt<fnptb) {
	403	effLowp =fhEffLPt->GetBinContent(iTheta,iPhi,iPt+1);
	404	effHighp=fhEffHPt->GetBinContent(iTheta,iPhi,iPt+1);
	405	effAnyp =fhEffAPt->GetBinContent(iTheta,iPhi,iPt+1);
	406	#ifdef MYTRIGDEBUG
	407	printf(" 7:1 A simple look to next point: %d: %f %f %f\n",iPt-1,effLowp,effHighp,effAnyp);
	408	#endif
	409	} else {
	410	effLowp=effLow;
	411	effHighp=effHigh;
	412	effAnyp=effAny;
	413	#ifdef MYTRIGDEBUG
	414	printf(" 8:1 result is: pt: %f %f %f\n",effLowp,effHighp,effAnyp);
	415	#endif
	416	}
	417	angl=(effLowp-effLow)/fDpt;
	418	angh=(effHighp-effHigh)/fDpt;
	419	anga=(effAnyp-effAny)/fDpt;
	420	}
	421	effLow=effLow+angl*(pt-ptc);
	422	effHigh=effHigh+angh*(pt-ptc);
	423	effAny=effAny+anga*(pt-ptc);
	424	#ifdef MYTRIGDEBUG
	425	printf(" 9:1 the interpolation coefficients are: %f %f %f\n",angl,angh,anga);
426	#endif
427	}
428	#ifdef MYTRIGDEBUG
429	printf("10:%1d effLow, effHigh=%f %f %f\n",fInt,effLow,effHigh,effAny);
430	#endif
431	return;
432	}
433
bc77d57e	434
6255180c	435
	436	Float_t AliFastMuonTriggerEff::Evaluate(Float_t charge, Float_t pt,
	437	Float_t theta, Float_t phi)
	438	{
	439	//
	440	// Trigger efficiency for pt, theta, phi depending of fCut
	441	//
	442	if(fZones==0) {
	443	printf("Call to uninitialized object of class: AliFastMuonTriggerEff\n");
	444	return 0.;
	445	}
	446	Float_t eff;
490243f3	447	Float_t effLow, effHigh, effAny;
6255180c	448
490243f3	449	Evaluate(charge,pt,theta,phi,effLow,effHigh,effAny);
6255180c	450	if (fCut == kLow)
6255180c	451	eff = effLow;
490243f3	452	else if (fCut == kHigh)
6255180c	453	eff = effHigh;
490243f3	454	else if (fCut == kAny)
	455	eff = effAny;
	456	else
	457	eff = 0;
6255180c	458
	459	return eff;
	460	}