#include "AliEMCALHit.h"
#include "AliEMCALSDigitizer.h"
#include "AliEMCALGeometry.h"
-#include "AliEMCALHistoUtilities.h"
+//JLK
+//#include "AliEMCALHistoUtilities.h"
ClassImp(AliEMCALSDigitizer)
fSDigitsInRun(0),
fFirstEvent(0),
fLastEvent(0),
- fSampling(0.),
- fControlHists(0),
- fHists(0)
+ fSampling(0.)
+ //JLK
+ //fControlHists(0),
+ //fHists(0)
{
// ctor
InitParameters();
fSDigitsInRun(0),
fFirstEvent(0),
fLastEvent(0),
- fSampling(0.),
- fControlHists(1),
- fHists(0)
+ fSampling(0.)
+ //JLK
+ //fControlHists(1),
+ //fHists(0)
{
// ctor
Init();
InitParameters() ;
- if(fControlHists) BookControlHists(1);
+
+ //JLK
+ //if(fControlHists) BookControlHists(1);
}
fSDigitsInRun(sd.fSDigitsInRun),
fFirstEvent(sd.fFirstEvent),
fLastEvent(sd.fLastEvent),
- fSampling(sd.fSampling),
- fControlHists(sd.fControlHists),
- fHists(sd.fHists)
+ fSampling(sd.fSampling)
+ //JLK
+ //fControlHists(sd.fControlHists),
+ //fHists(sd.fHists)
{
//cpy ctor
}
}
// Initializes parameters
fA = 0;
- fB = 1.e+7;
+ fB = 1.e+6; // Changed 24 Apr 2007. Dynamic range now 2 TeV
fSampling = geom->GetSampling();
- // threshold for deposit energy of hit
- fECPrimThreshold = 0.; // 24-nov-04 - was 1.e-6;
-
+ // threshold for deposit energy of hit
+ fECPrimThreshold = 0.05;// GeV // 22-may-07 was 0// 24-nov-04 - was 1.e-6;
+
AliDebug(2,Form("Print: \n------------------- %s -------------\n",GetName()));
AliDebug(2,Form(" fInit %i\n", int(fInit)));
AliDebug(2,Form(" fFirstEvent %i\n", fFirstEvent));
// Assign primary number only if deposited energy is significant
curSDigit = new AliEMCALDigit( hit->GetPrimary(),
hit->GetIparent(), hit->GetId(),
- Digitize(energy), hit->GetTime() ) ;
+ Digitize(energy), hit->GetTime(),
+ -1, energy ) ;
else
curSDigit = new AliEMCALDigit( -1 ,
-1 ,
hit->GetId(),
- Digitize(energy), hit->GetTime() ) ;
+ Digitize(energy), hit->GetTime(),
+ -1, energy ) ;
} else {
Warning("Exec"," abs id %i is bad \n", hit->GetId());
newsdigit = kFALSE;
nSdigits = sdigits->GetEntriesFast() ;
fSDigitsInRun += nSdigits ;
- Double_t e=0.,esum=0.;
- AliEMCALHistoUtilities::FillH1(fHists, 0, double(sdigits->GetEntriesFast()));
+ //JLK
+ //Double_t e=0.,esum=0.;
+ //AliEMCALHistoUtilities::FillH1(fHists, 0, double(sdigits->GetEntriesFast()));
for (i = 0 ; i < sdigits->GetEntriesFast() ; i++) {
AliEMCALDigit * sdigit = dynamic_cast<AliEMCALDigit *>(sdigits->At(i)) ;
sdigit->SetIndexInList(i) ;
- AliEMCALHistoUtilities::FillH1(fHists, 2, double(sdigit->GetAmp()));
- e = double(Calibrate(sdigit->GetAmp()));
- esum += e;
- AliEMCALHistoUtilities::FillH1(fHists, 3, e);
- AliEMCALHistoUtilities::FillH1(fHists, 4, double(sdigit->GetId()));
+ //JLK
+ //AliEMCALHistoUtilities::FillH1(fHists, 2, double(sdigit->GetAmp()));
+ //e = double(Calibrate(sdigit->GetAmp()));
+ //esum += e;
+ //AliEMCALHistoUtilities::FillH1(fHists, 3, e);
+ //AliEMCALHistoUtilities::FillH1(fHists, 4, double(sdigit->GetId()));
}
- if(esum>0.) AliEMCALHistoUtilities::FillH1(fHists, 1, esum);
+ //JLK
+ //if(esum>0.) AliEMCALHistoUtilities::FillH1(fHists, 1, esum);
// Now write SDigits
}
//__________________________________________________________________
-
Int_t AliEMCALSDigitizer::Digitize(Float_t energy)const {
// Digitize the energy
Double_t aSignal = fA + energy*fB;
void AliEMCALSDigitizer::PrintSDigits(Option_t * option)
{
//Prints list of digits produced at the current pass of AliEMCALDigitizer
-
-
- // AliEMCALGetter * gime = AliEMCALGetter::Instance() ;
+
AliEMCALLoader *rl = dynamic_cast<AliEMCALLoader*>(AliRunLoader::GetRunLoader()->GetDetectorLoader("EMCAL"));
const TClonesArray * sdigits = rl->SDigits() ;
rl->UnloadSDigits() ;
}
+//____________________________________________________________________________
void AliEMCALSDigitizer::Browse(TBrowser* b)
{
- if(fHists) b->Add(fHists);
+ //JLK
+ //if(fHists) b->Add(fHists);
TTask::Browse(b);
}
+/*
+//____________________________________________________________________________
TList *AliEMCALSDigitizer::BookControlHists(int var)
{
//book histograms for monitoring sdigitization
return fHists;
}
+//____________________________________________________________________________
void AliEMCALSDigitizer::SaveHists(const char* name, Bool_t kSingleKey, const char* opt)
{
AliEMCALHistoUtilities::SaveListOfHists(fHists, name, kSingleKey, opt);
}
+*/