* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
+/* $Id$ */
+
+#include <TMath.h>
+#include <TFile.h>
+#include <TROOT.h>
+#include <TH3.h>
-*/
#include "AliMUONTriggerLut.h"
-#include "TTree.h"
+#include "AliMUONTriggerCircuit.h"
#include "AliRun.h"
#include "AliMUON.h"
-#include "AliMUONPoints.h"
-#include "TMath.h"
-#include "TFile.h"
-#include "TH3.h"
-#include "iostream.h"
ClassImp(AliMUONTriggerLut)
//----------------------------------------------------------------------
-AliMUONTriggerLut::AliMUONTriggerLut() {
+AliMUONTriggerLut::AliMUONTriggerLut()
+ : TNamed()
+{
// constructor
fLptPlus = fLptMinu = fLptUnde = 0;
fHptPlus = fHptMinu = fHptUnde = 0;
}
//----------------------------------------------------------------------
-AliMUONTriggerLut::AliMUONTriggerLut (const AliMUONTriggerLut& MUONTriggerLut)
+AliMUONTriggerLut::AliMUONTriggerLut (const AliMUONTriggerLut& theMUONTriggerLut)
+ : TNamed(theMUONTriggerLut)
{
-// Dummy copy constructor
+// Protected copy constructor
+
+ Fatal("AliMUONTriggerLut", "Not implemented.");
}
//----------------------------------------------------------------------
-AliMUONTriggerLut & AliMUONTriggerLut::operator=(const AliMUONTriggerLut& MUONTriggerLut)
+AliMUONTriggerLut &
+AliMUONTriggerLut::operator=(const AliMUONTriggerLut& rhs)
{
-// Dummy assignment operator
- return *this;
+// Protected assignement operator
+
+ if (this == &rhs) return *this;
+
+ Fatal("operator=", "Not implemented.");
+
+ return *this;
}
//----------------------------------------------------------------------
static TFile *fileLut;
static Bool_t first=kTRUE;
if(first) {
- cout << " opening MUONTriggerLut.root " << "\n";
+ // printf(" opening MUONTriggerLut.root \n");
fileLut = new TFile("$(ALICE_ROOT)/MUON/MUONTriggerLut.root","READ");
first=kFALSE;
}
if ((binc & mask)!=0) lutApt[0]=lutApt[1]=1;
// get back to the first file
+/* no need for this with NewIO
TTree *tK = gAlice->TreeK();
TFile *file1 = 0;
if (tK) file1 = tK->GetCurrentFile();
file1->cd();
+*/
}
//----------------------------------------------------------------------
Int_t mask=0;
tabMask[ystrip]=1;
for (Int_t i=0; i<16; i++) {
- mask=mask+Int_t(tabMask[i]*pow(2,i));
+ mask=mask+Int_t(tabMask[i]*TMath::Power(2,i));
}
return mask;
}
char fileName[60];
sprintf(fileName,"$(ALICE_ROOT)/MUON/MUONTriggerLut.root");
- cout << " file name is " << fileName << "\n";
+ printf(" file name is %s\n",fileName);
// open output file containing histos
TFile *hfile = new TFile(fileName,"RECREATE","Trigger Look Up Table");
AliMUONTriggerCircuit* triggerCircuit;
for (Int_t icirc=0; icirc<234; icirc++) {
- cout << " Loading LuT for circuit " << icirc << " of 234 " << "\n";
+ printf(" Loading LuT for circuit %d \n",icirc);
triggerCircuit = &(pMUON->TriggerCircuit(icirc));
for (Int_t istripX=0; istripX<31; istripX++) {
Float_t pt=triggerCircuit->PtCal(istripX,idev,istripY);
if (pt>lptTreshold) {
- if (idev<15) iLptMinu=iLptMinu+Int_t(pow(2,istripY));
- else if (idev==15) iLptUnde=iLptUnde+Int_t(pow(2,istripY));
- else if (idev>15) iLptPlus=iLptPlus+Int_t(pow(2,istripY));
+ if (idev<15) iLptMinu=iLptMinu+Int_t(TMath::Power(2,istripY));
+ else if (idev==15) iLptUnde=iLptUnde+Int_t(TMath::Power(2,istripY));
+ else if (idev>15) iLptPlus=iLptPlus+Int_t(TMath::Power(2,istripY));
}
if (pt>hptTreshold) {
- if (idev<15) iHptMinu=iHptMinu+Int_t(pow(2,istripY));
- else if (idev==15) iHptUnde=iHptUnde+Int_t(pow(2,istripY));
- else if (idev>15) iHptPlus=iHptPlus+Int_t(pow(2,istripY));
+ if (idev<15) iHptMinu=iHptMinu+Int_t(TMath::Power(2,istripY));
+ else if (idev==15) iHptUnde=iHptUnde+Int_t(TMath::Power(2,istripY));
+ else if (idev>15) iHptPlus=iHptPlus+Int_t(TMath::Power(2,istripY));
}
- if (idev<15) iAptMinu=iAptMinu+Int_t(pow(2,istripY));
- else if (idev==15) iAptUnde=iAptUnde+Int_t(pow(2,istripY));
- else if (idev>15) iAptPlus=iAptPlus+Int_t(pow(2,istripY));
+ if (idev<15) iAptMinu=iAptMinu+Int_t(TMath::Power(2,istripY));
+ else if (idev==15) iAptUnde=iAptUnde+Int_t(TMath::Power(2,istripY));
+ else if (idev>15) iAptPlus=iAptPlus+Int_t(TMath::Power(2,istripY));
} // loop on istripY