Add2RawsList(hRawSumQZPC, 16, expert, !image);
Add2RawsList(hRawSumQZPA, 17, expert, !image);
- TH1F * hRawTDCZEM1 = new TH1F("hRawTDCZEM1", "Raw TDC ZEM1;TDC [ns]",160, -100., -60.);
+ TH1F * hRawTDCZEM1 = new TH1F("hRawTDCZEM1", "Raw TDC ZEM1;TDC [ns]",160, -350., -310.);
Add2RawsList(hRawTDCZEM1, 18, expert, !image);
- TH1F * hRawTDCZPC = new TH1F("hRawTDCZPC", "Raw TDC ZPC;TDC [ns]",160, -100., -60.);
+ TH1F * hRawTDCZPC = new TH1F("hRawTDCZPC", "Raw TDC ZPC;TDC [ns]",160, -350., -310.);
Add2RawsList(hRawTDCZPC, 19, expert, !image);
TProfile * hRawADCProfs = new TProfile("hRawADCProfs", "ADC profiles;ADC id;Mean ADC values",22,-0.5,21.5,10.,1210.,"");
Add2RawsList(hRawADCProfs, 20, expert, !image);
- TProfile * hRawTDCProfs = new TProfile("hRawTDCProfs", "TDC profiles;TDC id;Mean TDC values",6,0.5,6.5,-100.,-60.,"S");
+ TProfile * hRawTDCProfs = new TProfile("hRawTDCProfs", "TDC profiles;TDC id;Mean TDC values",6,0.5,6.5,-340.,-300.,"S");
Add2RawsList(hRawTDCProfs, 21, expert, !image);
TH1F * hRawADCs = new TH1F("hRawADCs", "ADCs;ADC id;Mean ADC values",22,-0.5,21.5);
- hRawADCs->SetMinimum(10);
Add2RawsList(hRawADCs, 22, !expert, image);
TH1F * hRawTDCs = new TH1F("hRawTDCs", "TDCs;TDC id;Mean TDC values",6,0.5,6.5);
- hRawTDCs->SetMaximum(-60);
+ hRawTDCs->SetMaximum(-300); hRawTDCs->SetMinimum(-340);
Add2RawsList(hRawTDCs, 23, !expert, image);
TH2F *hZNCrawCentr = new TH2F("hZNCrawCentr", "Centroid in ZNC;X (cm);Y(cm)", 100, -5.,5.,100,-5.,5.);
TH2F *hZNArawCentr = new TH2F("hZNArawCentr", "Centroid in ZNA;X (cm);Y(cm)", 100, -5.,5.,100,-5.,5.);
Add2RawsList(hZNArawCentr, 25, expert, !image);
- TH1F *hRawCentr = new TH1F("hRawCentr","Centroid in ZNCs; ZNC-ZNA coordinates; Coordinate value (cm)", 4, 0.5, 4.5);
- Add2RawsList(hRawCentr, 26, !expert, image);
+ TH2F *hTimeZDC = new TH2F("fHistTimeZDC", "ZDC;TDC timing A+C vs C-A; events", 120,-30,30,120,-100,-40);
+ Add2RawsList(hTimeZDC, 26, expert, !image);
}
//____________________________________________________________________________
Bool_t isZNCFired=kFALSE, isZPCFired=kFALSE, isZNAFired=kFALSE, isZPAFired=kFALSE;
Int_t indZNC=0, indZNA=0, indZPC=0, indZPA=0;
Float_t zncTDC[10], zpcTDC[10], zem1TDC[10], zem2TDC[10], znaTDC[10], zpaTDC[10];
+ Float_t zncSumTDC[10], znaSumTDC[10];
for(Int_t i=0; i<10; i++){
- zncTDC[i]=zpcTDC[i]=zem1TDC[i]=zem2TDC[i]=znaTDC[i]=zpaTDC[i]=-999.;
+ zncTDC[i]=zpcTDC[i]=zem1TDC[i]=zem2TDC[i]=znaTDC[i]=zpaTDC[i]=zncSumTDC[i]=znaSumTDC[i]=-999.;
}
- Float_t tdcL0=-999.;
+ Float_t tdcGate=-999., l0=-999.;
Int_t iMultZNCTDC=0, iMultZPCTDC=0, iMultZEM1TDC=0, iMultZEM2TDC=0, iMultZNATDC=0, iMultZPATDC=0;
+ Int_t iMultTDCC=0, iMultTDCA=0;
const Float_t x[4] = {-1.75, 1.75, -1.75, 1.75};
const Float_t y[4] = {-1.75, -1.75, 1.75, 1.75};
} //IsADCDataWord && signal ADCs
else if(stream.IsZDCTDCDatum()){
if(stream.GetChannel()==1){
- zncTDC[iMultZNCTDC] = (0.025*(stream.GetZDCTDCDatum()));
+ zncTDC[iMultZNCTDC] = (0.025*stream.GetZDCTDCDatum());
iMultZNCTDC++;
}
else if(stream.GetChannel()==3){
- zpcTDC[iMultZPCTDC] = (0.025*(stream.GetZDCTDCDatum()));
+ zpcTDC[iMultZPCTDC] = (0.025*stream.GetZDCTDCDatum());
iMultZPCTDC++;
}
else if(stream.GetChannel()==5){
- znaTDC[iMultZNATDC] = (0.025*(stream.GetZDCTDCDatum()));
+ znaTDC[iMultZNATDC] = (0.025*stream.GetZDCTDCDatum());
iMultZNATDC++;
}
else if(stream.GetChannel()==7){
- zpaTDC[iMultZPATDC] = (0.025*(stream.GetZDCTDCDatum()));
+ zpaTDC[iMultZPATDC] = (0.025*stream.GetZDCTDCDatum());
iMultZPATDC++;
}
else if(stream.GetChannel()==8){
- zem1TDC[iMultZEM1TDC] = (0.025*(stream.GetZDCTDCDatum()));
+ zem1TDC[iMultZEM1TDC] = (0.025*stream.GetZDCTDCDatum());
iMultZEM1TDC++;
}
else if(stream.GetChannel()==9){
- zem2TDC[iMultZEM2TDC] = (0.025*(stream.GetZDCTDCDatum()));
+ zem2TDC[iMultZEM2TDC] = (0.025*stream.GetZDCTDCDatum());
iMultZEM2TDC++;
}
- else if(stream.GetChannel()==16) tdcL0 = (0.025*(stream.GetZDCTDCDatum()));
+ else if(stream.GetChannel()==10){
+ zncSumTDC[iMultZEM2TDC] = (0.025*stream.GetZDCTDCDatum());
+ iMultTDCC++;
+ }
+ else if(stream.GetChannel()==12){
+ znaSumTDC[iMultZEM2TDC] = (0.025*stream.GetZDCTDCDatum());
+ iMultTDCA++;
+ }
+ else if(stream.GetChannel()==14) tdcGate = (0.025*stream.GetZDCTDCDatum());
+ else if(stream.GetChannel()==15) l0 = (0.025*stream.GetZDCTDCDatum());
- if(stream.GetChannel()==16 && tdcL0!=-999.){
+ if(stream.GetChannel()==16 && tdcGate!=-999.){
for(Int_t iHit=0; iHit<10; iHit++){
if(zncTDC[iHit]!=-999.){
- GetRawsData(21)->Fill(1, zncTDC[iHit]-tdcL0);
+ if(zncTDC[iHit]-tdcGate>-340. && zncTDC[iHit]-tdcGate<-300.)
+ GetRawsData(21)->Fill(1, zncTDC[iHit]-tdcGate);
//GetRawsData(23)->SetBinContent(1, GetRawsData(21)->GetBinContent(1));
}
if(zpcTDC[iHit]!=-999.){
- Float_t diffZPC = zpcTDC[iHit]-tdcL0;
+ Float_t diffZPC = zpcTDC[iHit]-tdcGate;
GetRawsData(19)->Fill(diffZPC);
- GetRawsData(21)->Fill(2, diffZPC);
+ if(diffZPC>-340. && diffZPC<-300.) GetRawsData(21)->Fill(2, diffZPC);
//GetRawsData(23)->SetBinContent(2, GetRawsData(21)->GetBinContent(3));
}
if(znaTDC[iHit]!=-999.){
- GetRawsData(21)->Fill(3, znaTDC[iHit]-tdcL0);
+ if(znaTDC[iHit]-tdcGate>-340. && znaTDC[iHit]-tdcGate<-300.)
+ GetRawsData(21)->Fill(3, znaTDC[iHit]-tdcGate);
//GetRawsData(23)->SetBinContent(3, GetRawsData(21)->GetBinContent(5));
}
if(zpaTDC[iHit]!=-999.){
- GetRawsData(21)->Fill(4, zpaTDC[iHit]-tdcL0);
+ if(zpaTDC[iHit]-tdcGate>-340. && zpaTDC[iHit]-tdcGate<-300.)
+ GetRawsData(21)->Fill(4, zpaTDC[iHit]-tdcGate);
//GetRawsData(23)->SetBinContent(4, GetRawsData(21)->GetBinContent(7));
}
if(zem1TDC[iHit]!=-999.){
- Float_t diffZEM1 = zem1TDC[iHit]-tdcL0;
+ Float_t diffZEM1 = zem1TDC[iHit]-tdcGate;
GetRawsData(18)->Fill(diffZEM1);
- GetRawsData(21)->Fill(5, diffZEM1);
+ if(diffZEM1>-340. && diffZEM1<-300.) GetRawsData(21)->Fill(5, diffZEM1);
//GetRawsData(23)->SetBinContent(5, GetRawsData(21)->GetBinContent(8));
}
if(zem2TDC[iHit]!=-999.){
- GetRawsData(21)->Fill(6, zem2TDC[iHit]-tdcL0);
+ if(zem2TDC[iHit]-tdcGate>-340. && zem2TDC[iHit]-tdcGate<-300.)
+ GetRawsData(21)->Fill(6, zem2TDC[iHit]-tdcGate);
//GetRawsData(23)->SetBinContent(6, GetRawsData(21)->GetBinContent(9));
}
+ if(zncSumTDC[iHit]!=-999.){
+ Float_t tdcC = zncSumTDC[iHit]-l0;
+ if(znaSumTDC[iHit]!=-999.){
+ Float_t tdcA = znaSumTDC[iHit]-l0;
+ //if (((tdcC-tdcA-refDelta)*(tdcC-tdcA-refDelta)/(sigmaDelta*sigmaDelta) +
+ //(tdcC+tdcA-refSum)*(tdcC+tdcA-refSum)/(sigmaSum*sigmaSum))< 1.0)
+ GetRawsData(26)->Fill(tdcC-tdcA,tdcC+tdcA);
+
+ }
+ }
}
//
- tdcL0 = -999.;
+ tdcGate = -999.;
for(Int_t i=0; i<10; i++){
- zpcTDC[i] = zem1TDC[i] = -999.;
+ zncTDC[i] = zpcTDC[i] = zem1TDC[i] = zem2TDC[i] = znaTDC[i] = zpaTDC[i] = -999.;
+ zncSumTDC[i] = znaSumTDC[i] = -999.;
}
}
}
GetRawsData(23)->SetBinError(ibin, GetRawsData(21)->GetBinError(ibin));
GetRawsData(23)->SetLineColor(kAzure-3); GetRawsData(23)->SetLineWidth(2);
}
- //
- Double_t xznc = GetRawsData(24)->GetMean(1);
- Double_t yznc = GetRawsData(24)->GetMean(2);
- GetRawsData(26)->SetBinContent(1., xznc);
- GetRawsData(26)->SetBinContent(2., yznc);
- Double_t xzna = GetRawsData(25)->GetMean(1);
- Double_t yzna = GetRawsData(25)->GetMean(2);
- GetRawsData(26)->SetBinContent(3., xzna);
- GetRawsData(26)->SetBinContent(4., yzna);
- GetRawsData(26)->SetLineColor(kPink-9); GetRawsData(26)->SetLineWidth(2);
- GetRawsData(26)->SetOption("E");
}
}