ClassImp(AliEveT0Module)
/******************************************************************************/
-AliEveT0Module::AliEveT0Module(const Text_t* n, Int_t sigType, AliT0digit *digits, AliT0RawReader *start)
- : TEveQuadSet(n), fSigType(sigType), fDigits(digits), fStart(start)
+AliEveT0Module::AliEveT0Module(const Text_t* n, Int_t sigType,
+ AliT0digit *digits, AliT0RawReader *start,
+ Double_t zvertex) :
+ TEveQuadSet(n), fSigType(sigType),
+ fDigits(digits), fStart(start),
+ fZVertex(zvertex)
{
//
// Default constructor
void AliEveT0Module::PrintEventInfo()
{
- printf("Blabla\n");
+ printf("Z-vertex = %f\n", fZVertex);
}
/******************************************************************************/
AliT0RawReader *start = new AliT0RawReader(reader);
Int_t allData[110][5];
- TRandom r(0);
+
+ Float_t angle = 2 * TMath::Pi() / 12;
+ start->Next();
+ for (Int_t i=0; i<110; i++)
+ {
+ for (Int_t iHit=0; iHit<5; iHit++)
+ {
+ allData[i][iHit]= start->GetData(i,iHit);
+ }
+ }
+ Double_t zvertex = (allData[51][0] - allData[52][0])/2*25*2.99752/100;
+ //std::cout<<"zvertex= "<< zvertex <<endl;
TEveRGBAPalette *rawPalette = new TEveRGBAPalette(0, 3000);
TEveRGBAPalette *vertexPalette = new TEveRGBAPalette(-100, 100);
rawPalette->SetLimits(1, 3000); // Set proper raw time range.
- TEveQuadSet* rawA = new AliEveT0Module("T0_RAW_A", 2, digits, start);
+ TEveQuadSet* rawA = new AliEveT0Module("T0_RAW_A", 2, digits, start, zvertex);
rawA->SetPalette(rawPalette);
rawA->Reset(TEveQuadSet::kQT_HexagonXY, kFALSE, 32);
rawA->SetFrame(box);
- TEveQuadSet* rawC = new AliEveT0Module("T0_RAW_C", 3, digits, start);
+ TEveQuadSet* rawC = new AliEveT0Module("T0_RAW_C", 3, digits, start, zvertex);
rawC->SetPalette(rawPalette);
rawC->Reset(TEveQuadSet::kQT_HexagonXY, kFALSE, 32);
rawC->SetFrame(box);
- TEveQuadSet* vertexT0 = new AliEveT0Module("T0_Vertex", 5, digits, start);
+ TEveQuadSet* vertexT0 = new AliEveT0Module("T0_Vertex", 5, digits, start, zvertex);
vertexT0->SetPalette(vertexPalette);
vertexT0->Reset(TEveQuadSet::kQT_HexagonXY, kFALSE, 32);
- Float_t angle = 2 * TMath::Pi() / 12;
- start->Next();
- for (Int_t i=0; i<110; i++)
- {
- for (Int_t iHit=0; iHit<5; iHit++)
- {
- allData[i][iHit]= start->GetData(i,iHit);
- if (allData[i][iHit] != 0) {
- using namespace std;
- cout <<"event "<< AliEveEventManager::GetCurrent()->GetEventId()
- <<" i "<< i <<" "<< allData[i][iHit] - allData[0][0] <<endl;
- }
- }
- }
- Float_t zvertex= (allData[51][0] - allData[52][0])/2*25*2.99752/100;
- using namespace std;
- cout<<"zvertex= "<< zvertex <<endl;
for (Int_t i=0; i<12; i++)
{
Float_t x = 6.5 * TMath::Sin(i * angle);
Float_t y = 6.5 * TMath::Cos(i * angle);
- rawA->AddHexagon(x, y, r.Uniform(-0.1, 0.1), 1.0);
+ rawA->AddHexagon(x, y, 0, 1.0);
rawA->QuadValue(start->GetData(i+37,0)-start->GetData(0,0));
- rawC->AddHexagon(x, y, r.Uniform(-0.1, 0.1), 1.0);
+ rawC->AddHexagon(x, y, 0, 1.0);
rawC->QuadValue(start->GetData(i+25,0)-start->GetData(0,0));
}
- vertexT0->AddHexagon(0, 0, 0, 1);
- vertexT0->QuadValue(TMath::Nint(zvertex));
+ vertexT0->AddHexagon(0, 0, 0, 1);
+ vertexT0->QuadValue(TMath::Nint(zvertex));
rawA->RefitPlex();
rawC->RefitPlex();
{
// Make modules for digits.
- TRandom r(0);
TArrayI adc(24);
TArrayI tdc(24);
digits->GetQT1(adc);
digits->GetTimeCFD(tdc);
- // printf("%3d\n",besttimeright);
- for (Int_t i=0;i<24; i++){
- printf("%3d %3d\n ",adc[i], tdc[i]);
- }
+ // printf("%3d\n",besttimeright);
+ // for (Int_t i=0;i<24; i++) {
+ // printf("%3d %3d\n ",adc[i], tdc[i]);
+ // }
TEveRGBAPalette* adcPalette = new TEveRGBAPalette(5, 1024);
adcPalette->SetLimits(1, 1024); // Set proper ADC range.
Float_t x = 6.5 * TMath::Sin(i * angle);
Float_t y = 6.5 * TMath::Cos(i * angle);
- qa->AddHexagon(x, y, r.Uniform(-0.1, 0.1), 1.0);
+ qa->AddHexagon(x, y, 0, 1.0);
qa->QuadValue(adc[i+12]);
// qa->QuadId(new TNamed(Form("PMT with idx=%d", i), "PMT's aplitude in side A."));
- qat->AddHexagon(x, y, r.Uniform(-0.1, 0.1), 1.0);
+ qat->AddHexagon(x, y, 0, 1.0);
qat->QuadValue(tdc[i+12]);
// qat->QuadId(new TNamed(Form("Quad with idx=%d", i), "PMT's time in side A."));
- qc->AddHexagon(x, y, r.Uniform(-0.1, 0.1), 1.0);
+ qc->AddHexagon(x, y, 0, 1.0);
qc->QuadValue(adc[i]);
// qc->QuadId(new TNamed(Form("Quad with idx=%d", i), "PMT's amplitude in side C."));
- qct->AddHexagon(x, y, r.Uniform(-0.1, 0.1), 1.0);
+ qct->AddHexagon(x, y, 0, 1.0);
qct->QuadValue(tdc[i]);
// qct->QuadId(new TNamed(Form("Quad with idx=%d", i), "PMT's time in side C."));
}
Int_t timediff = fDigits->TimeDiff();
Int_t summult = fDigits->SumMult();
- printf("besttimeA=%3d\n",besttimeright);
- printf("besttimeC=%3d\n",besttimeleft);
- printf("meantime=%3d\n",meantime);
- printf("timediff=%3d\n",timediff);
- printf("summult=%3d\n",summult);
+ printf("besttimeA = %3d\n", besttimeright);
+ printf("besttimeC = %3d\n", besttimeleft);
+ printf("meantime = %3d\n", meantime);
+ printf("timediff = %3d\n", timediff);
+ printf("summult = %3d\n", summult);
printf(" idx=%d, amplitude=%d\n", idx, qb->fValue);
Int_t timediff = fDigits->TimeDiff();
Int_t summult = fDigits->SumMult();
- printf("besttimeA=%3d\n",besttimeright);
- printf("besttimeC=%3d\n",besttimeleft);
- printf("meantime=%3d\n",meantime);
- printf("timediff=%3d\n",timediff);
- printf("summult=%3d\n",summult);
+ printf("besttimeA = %3d\n",besttimeright);
+ printf("besttimeC = %3d\n",besttimeleft);
+ printf("meantime = %3d\n",meantime);
+ printf("timediff = %3d\n",timediff);
+ printf("summult = %3d\n",summult);
printf(" idx=%d, amplitude=%d\n", idx, qb->fValue);
}
if (fSigType == 2) {
printf("raw====================\n");
- printf("besttimeA=%3d\n",fStart->GetData(51,0)-fStart->GetData(0,0));
- printf("besttimeC=%3d\n",fStart->GetData(52,0)-fStart->GetData(0,0));
- printf("meantime=%3d\n",fStart->GetData(49,0)-fStart->GetData(0,0));
- printf("amplitude= %3d\n",fStart->GetData(idx+1,0));
+ printf("besttimeA = %3d\n",fStart->GetData(51,0)-fStart->GetData(0,0));
+ printf("besttimeC = %3d\n",fStart->GetData(52,0)-fStart->GetData(0,0));
+ printf("meantime = %3d\n",fStart->GetData(49,0)-fStart->GetData(0,0));
+ printf("amplitude = %3d\n",fStart->GetData(idx+1,0));
printf(" idx=%d, time %d\n", idx, qb->fValue);
}
if (fSigType == 3) {
printf("raw====================\n");
- printf("besttimeA=%3d\n",fStart->GetData(51,0)-fStart->GetData(0,0));
- printf("besttimeC=%3d\n",fStart->GetData(52,0)-fStart->GetData(0,0));
- printf("meantime=%3d\n",fStart->GetData(49,0)-fStart->GetData(0,0));
- printf("amplitude= %3d\n",fStart->GetData(idx+13,0));
+ printf("besttimeA = %3d\n",fStart->GetData(51,0)-fStart->GetData(0,0));
+ printf("besttimeC = %3d\n",fStart->GetData(52,0)-fStart->GetData(0,0));
+ printf("meantime = %3d\n",fStart->GetData(49,0)-fStart->GetData(0,0));
+ printf("amplitude = %3d\n",fStart->GetData(idx+13,0));
printf(" idx=%d, time %d\n", idx, qb->fValue);
}