/* $Id$ */
+/* History of cvs commits:
+ *
+ * $Log$
+ * Revision 1.86 2005/07/12 20:07:35 hristov
+ * Changes needed to run simulation and reconstrruction in the same AliRoot session
+ *
+ * Revision 1.85 2005/05/28 14:19:04 schutz
+ * Compilation warnings fixed by T.P.
+ *
+ */
//_________________________________________________________________________
//*-- Author : Dmitri Peressounko (SUBATECH & Kurchatov Institute)
// --- Standard library ---
// --- AliRoot header files ---
-
+#include "AliLog.h"
#include "AliRunDigitizer.h"
#include "AliPHOSDigit.h"
#include "AliPHOSGetter.h"
//____________________________________________________________________________
AliPHOSDigitizer::~AliPHOSDigitizer()
{
+ AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle()) ;
+
+ // Clean Digitizer from the white board
+ gime->PhosLoader()->CleanDigitizer() ;
// dtor
- AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle(),fEventFolderName);
- gime->PhosLoader()->CleanDigitizer();
delete [] fInputFileNames ;
delete [] fEventNames ;
// This design avoids scanning over the list of digits to add
// contribution to new SDigits only.
- AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle(), fEventFolderName) ;
+ AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle()) ;
Int_t ReadEvent = event ;
if (fManager)
ReadEvent = dynamic_cast<AliStream*>(fManager->GetInputStream(0))->GetCurrentEventNumber() ;
- Info("Digitize", "Adding event %d from input stream 0 %s %s", ReadEvent, GetTitle(), fEventFolderName.Data()) ;
+ AliInfo(Form("Adding event %d from input stream 0 %s %s",
+ ReadEvent, GetTitle(), fEventFolderName.Data())) ;
gime->Event(ReadEvent, "S") ;
TClonesArray * digits = gime->Digits() ;
digits->Clear() ;
AliPHOSSDigitizer * sDigitizer = gime->SDigitizer();
if ( !sDigitizer )
- Fatal("Digitize", "SDigitizer with name %s %s not found", GetTitle(), fEventFolderName.Data() ) ;
+ AliFatal(Form("SDigitizer with name %s %s not found",
+ GetTitle(), fEventFolderName.Data() )) ;
//take all the inputs to add together and load the SDigits
TObjArray * sdigArray = new TObjArray(fInput) ;
AliPHOSGetter * gime = AliPHOSGetter::Instance(fInputFileNames[i], tempo) ;
if (fManager)
ReadEvent = dynamic_cast<AliStream*>(fManager->GetInputStream(i))->GetCurrentEventNumber() ;
- Info("Digitize", "Adding event %d from input stream %d %s %s", ReadEvent, i, fInputFileNames[i].Data(), tempo.Data()) ;
+ AliInfo(Form("Adding event %d from input stream %d %s %s",
+ ReadEvent, i, fInputFileNames[i].Data(), tempo.Data())) ;
gime->Event(ReadEvent,"S");
sdigArray->AddAt(gime->SDigits(), i) ;
}
{
// Returns digitized value of the energy in a cell absId
+ AliPHOSGetter* gime = AliPHOSGetter::Instance();
+
+ //Determine rel.position of the cell absId
+ Int_t relId[4];
+ gime->PHOSGeometry()->AbsToRelNumbering(absId,relId);
+ Int_t module=relId[0];
+ Int_t raw=relId[2];
+ Int_t column=relId[3];
+
Int_t chanel ;
+
if(absId <= fEmcCrystals){ //digitize as EMC
+
+ //reading calibration data for cell absId.
+ //If no calibration DB found, accept default values.
+
+ if(gime->CalibData()) {
+ fADCpedestalEmc = gime->CalibData()->GetADCpedestalEmc(module,column,raw);
+ fADCchanelEmc = gime->CalibData()->GetADCchannelEmc(module,column,raw);
+ }
+// printf("\t\tabsId %d ==>>module=%d column=%d raw=%d ped=%.4f gain=%.4f\n",
+// absId,module,column,raw,fADCpedestalEmc,fADCchanelEmc);
+
chanel = (Int_t) TMath::Ceil((energy - fADCpedestalEmc)/fADCchanelEmc) ;
if(chanel > fNADCemc ) chanel = fNADCemc ;
}
// by default fLastEvent = fFirstEvent (process only one event)
if (!fInit) { // to prevent overwrite existing file
- Error( "Exec", "Give a version name different from %s", fEventFolderName.Data() ) ;
+ AliError(Form("Give a version name different from %s",
+ fEventFolderName.Data() )) ;
return ;
}
AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle()) ;
+ // Post Digitizer to the white board
+ gime->PostDigitizer(this) ;
+
if (fLastEvent == -1)
fLastEvent = gime->MaxEvent() - 1 ;
else if (fManager)
fDigitsInRun += gime->Digits()->GetEntriesFast() ;
}
+ gime->PhosLoader()->CleanDigitizer();
+
if(strstr(option,"tim")){
gBenchmark->Stop("PHOSDigitizer");
TString message ;
message = " took %f seconds for Digitizing %f seconds per event\n" ;
- Info("Exec", message.Data(),
+ AliInfo(Form( message.Data(),
gBenchmark->GetCpuTime("PHOSDigitizer"),
- gBenchmark->GetCpuTime("PHOSDigitizer")/nEvents );
+ gBenchmark->GetCpuTime("PHOSDigitizer")/nEvents ));
}
}
fInit = kTRUE ;
AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle(), fEventFolderName) ;
if ( gime == 0 ) {
- Fatal("Init" ,"Could not obtain the Getter object for file %s and event %s !", GetTitle(), fEventFolderName.Data()) ;
+ AliFatal(Form("Could not obtain the Getter object for file %s and event %s !",
+ GetTitle(), fEventFolderName.Data()));
return kFALSE;
}
TString opt("Digits") ;
if(gime->VersionExists(opt) ) {
- Error( "Init", "Give a version name different from %s", fEventFolderName.Data() ) ;
+ AliError(Form("Give a version name different from %s",
+ fEventFolderName.Data() )) ;
fInit = kFALSE ;
}
- // Post Digitizer to the white board
- gime->PostDigitizer(this) ;
-
fFirstEvent = 0 ;
fLastEvent = fFirstEvent ;
if (fManager)
}
// looking for file which contains AliRun
if (gSystem->AccessPathName(alirunFileName)) {// file does not exist
- Error("MixWith", "File %s does not exist!", alirunFileName.Data()) ;
+ AliError(Form("File %s does not exist!", alirunFileName.Data())) ;
return ;
}
// looking for the file which contains SDigits
if ( eventFolderName != AliConfig::GetDefaultEventFolderName()) // only if not the default folder name
fileName = fileName.ReplaceAll(".root", "") + "_" + eventFolderName + ".root" ;
if ( (gSystem->AccessPathName(fileName)) ) {
- Error("MixWith", "The file %s does not exist!", fileName.Data()) ;
+ AliError(Form("The file %s does not exist!", fileName.Data())) ;
return ;
}
// need to increase the arrays
}
//__________________________________________________________________
-void AliPHOSDigitizer::Print()const
+void AliPHOSDigitizer::Print(const Option_t *)const
{
// Print Digitizer's parameters
- Info("Print", "\n------------------- %s -------------", GetName() ) ;
+ AliInfo(Form("\n------------------- %s -------------", GetName() )) ;
if( strcmp(fEventFolderName.Data(), "") != 0 ){
printf(" Writing Digits to branch with title %s\n", fEventFolderName.Data()) ;
printf(" ---------------------------------------------------\n") ;
}
else
- Info("Print", "AliPHOSDigitizer not initialized" ) ;
+ AliInfo(Form("AliPHOSDigitizer not initialized" )) ;
}
AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle(), fEventFolderName) ;
TClonesArray * digits = gime->Digits() ;
- Info("PrintDigits", "%d", digits->GetEntriesFast()) ;
+ AliInfo(Form("%d", digits->GetEntriesFast())) ;
printf("\nevent %d", gAlice->GetEvNumber()) ;
printf("\n Number of entries in Digits list %d", digits->GetEntriesFast() ) ;
//__________________________________________________________________
Float_t AliPHOSDigitizer::TimeOfNoise(void) const
{ // Calculates the time signal generated by noise
- Info("TimeOfNoise", "Change me") ;
+ //PH Info("TimeOfNoise", "Change me") ;
return gRandom->Rndm() * 1.28E-5;
}
// and branch "AliPHOSDigitizer", with the same title to keep all the parameters
// and names of files, from which digits are made.
- AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle(), fEventFolderName) ;
+ AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle()) ;
const TClonesArray * digits = gime->Digits() ;
TTree * treeD = gime->TreeD();
// -- create Digits branch
Int_t bufferSize = 32000 ;
- TBranch * digitsBranch = treeD->Branch("PHOS",&digits,bufferSize);
+ TBranch * digitsBranch = treeD->Branch("PHOS","TClonesArray",&digits,bufferSize);
digitsBranch->SetTitle(fEventFolderName);
digitsBranch->Fill() ;
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff