Skip to content

A
acc2go4
Commit 283f92c7 authored by mbsdaq's avatar mbsdaq
Browse files
Options

add Si det calib.

parent 756eb644
Show whitespace changes
Inline Side-by-side
Showing with 608 additions and 90 deletions
TNeAnalysis.cxx
View file @ 283f92c7
......@@ -31,11 +31,13 @@ TNeAnalysis::TNeAnalysis()
,p_Lya(0)
,p_Rxa(0)
,p_Rya(0)
,p_RVa(0)
,p_R60a(0)
,p_Rva(0)
,p_R20a(0)
,p_AnSa(0)
,p_AnRa(0)
,p_AnVa(0)
,p_Rcsa(0)
,p_Lcsa(0)
{
}
//***********************************************************
......@@ -50,11 +52,13 @@ TNeAnalysis::TNeAnalysis(int argc, char** argv)
,p_Lya(0)
,p_Rxa(0)
,p_Rya(0)
,p_RVa(0)
,p_R60a(0)
,p_Rva(0)
,p_R20a(0)
,p_AnSa(0)
,p_AnRa(0)
,p_AnVa(0)
,p_Rcsa(0)
,p_Lcsa(0)
{
cout << "**** Create TNeAnalysis name: " << argv[0] << endl;
......@@ -102,20 +106,24 @@ TNeAnalysis::TNeAnalysis(int argc, char** argv)
p_Lya = new TNeDet16("Lya");
p_Rxa = new TNeDet16("Rxa");
p_Rya = new TNeDet16("Rya");
p_RVa = new TNeDet16("RVa");
p_R60a = new TNeDet16("R60a");
p_Rva = new TNeDet16("Rva");
p_R20a = new TNeDet16("R20a");
p_AnSa = new TNeDet16("AnSa");
p_AnRa = new TNeDet16("AnRa");
p_AnVa = new TNeDet16("AnVa");
p_Rcsa = new TNeDet16("Rcsa");
p_Lcsa = new TNeDet16("Lcsa");
AddParameter(p_Lxa);
AddParameter(p_Lya);
AddParameter(p_Rxa);
AddParameter(p_Rya);
AddParameter(p_RVa);
AddParameter(p_R60a);
AddParameter(p_Rva);
AddParameter(p_R20a);
AddParameter(p_AnSa);
AddParameter(p_AnRa);
AddParameter(p_AnVa);
AddParameter(p_Rcsa);
AddParameter(p_Lcsa);
// pMTracks = new TNeMParms("MTRACKS");
// AddParameter(pMTracks);
......@@ -158,11 +166,13 @@ Int_t TNeAnalysis::UserPreLoop()
p_Lya->ReadData();
p_Rxa->ReadData();
p_Rya->ReadData();
p_RVa->ReadData();
p_R60a->ReadData();
p_Rva->ReadData();
p_R20a->ReadData();
p_AnSa->ReadData();
p_AnRa->ReadData();
p_AnVa->ReadData();
p_Rcsa->ReadData();
p_Lcsa->ReadData();
return 0;
......
TNeAnalysis.h
View file @ 283f92c7
......@@ -33,11 +33,13 @@ public:
TNeDet16* p_Lya;
TNeDet16* p_Rxa;
TNeDet16* p_Rya;
TNeDet16* p_RVa;
TNeDet16* p_R60a;
TNeDet16* p_Rva;
TNeDet16* p_R20a;
TNeDet16* p_AnSa;
TNeDet16* p_AnRa;
TNeDet16* p_AnVa;
TNeDet16* p_Rcsa;
TNeDet16* p_Lcsa;
// TNeMParms* pMTracks;
ClassDef(TNeAnalysis,1)
......
TNeEvent.cxx
View file @ 283f92c7
......@@ -42,8 +42,8 @@ void TNeEvent::Clear(Option_t *t)
memset(Rxt, 0, sizeof(Rxt));
memset(Rya, 0, sizeof(Rya));
memset(Ryt, 0, sizeof(Ryt));
memset(RVa, 0, sizeof(RVa));
memset(RVt, 0, sizeof(RVt));
memset(Rva, 0, sizeof(Rva));
memset(Rvt, 0, sizeof(Rvt));
memset(R20a, 0, sizeof(R20a));
memset(R20t, 0, sizeof(R20t));
memset(AnSa, 0, sizeof(AnSa));
......
TNeEvent.h
View file @ 283f92c7
......@@ -22,26 +22,47 @@ public:
//-----Si-----
unsigned short Lxa[16];
unsigned short Lxt[16];
unsigned short Lxc[16];
unsigned short Lya[16];
unsigned short Lyt[16];
unsigned short Lyc[16];
unsigned short Rxa[16];
unsigned short Rxt[16];
unsigned short Rxc[16];
unsigned short Rya[16];
unsigned short Ryt[16];
unsigned short RVa[16];
unsigned short RVt[16];
unsigned short Ryc[16];
unsigned short Rva[16];
unsigned short Rvt[16];
unsigned short Rvc[16];
unsigned short R20a[16];
unsigned short R20t[16];
unsigned short R20c[16];
unsigned short AnSa[16];
unsigned short AnSt[16];
unsigned short AnSc[16];
unsigned short AnRa[16];
unsigned short AnRt[16];
unsigned short AnRc[16];
unsigned short AnVa[16];
unsigned short AnVt[16];
unsigned short AnVc[16];
unsigned short Rcsa[16];
unsigned short Rcst[16];
unsigned short Rcsc[16];
unsigned short Lcsa[16];
unsigned short Lcst[16];
unsigned short Lcsc[16];
unsigned short F3[4];
unsigned short tF3[4];
......
TNeProc.cxx
View file @ 283f92c7
......@@ -76,6 +76,9 @@ TNeProc::TNeProc() : TGo4EventProcessor()
d1Rya[i] = NULL;
d1Ryt[i] = NULL;
d1Ryc[i] = NULL;
d1Rva[i] = NULL;
d1Rvt[i] = NULL;
d1Rvc[i] = NULL;
}
//_________________________________Neutrons
//
......@@ -100,6 +103,7 @@ TNeProc::TNeProc() : TGo4EventProcessor()
d1CY1 = NULL;
d1CX2 = NULL;
d1CY2 = NULL;
d2XYbeamOnTar = NULL;
} //-------------------------------------------------------------------
......@@ -131,45 +135,60 @@ TNeProc::TNeProc(const char* name) : TGo4EventProcessor(name)
d1Lyc[i] = MakeTH1('D', Form("Lyc/Lycs_%2.2d",i)
//
, Form("L x energy MeV %2.2d",i), 4000, -2., 98.);
d1Rxa[i] = MakeTH1('I', Form("Rxa/Rxas_%2.2d",i) //указывать не повторяющуюся структуру папок и уникальные имена
d1Rxa[i] = MakeTH1('I', Form("Rxa/Rxas_%2.2d",i)
, Form("R x amp CH %2.2d",i), 8192, 0., 8192.);
d1Rxt[i] = MakeTH1('I', Form("Rxt/Rxts_%2.2d",i)
, Form("R x time CH %2.2d",i), 4096, 1., 4096);
d1Rxc[i] = MakeTH1('D', Form("Rxc/Rxcs_%2.2d",i)
, Form("R x energy MeV %2.2d",i), 4000, -2., 98.);
//
d1Rya[i] = MakeTH1('I', Form("Rya/Ryas_%2.2d",i) //указывать не повторяющуюся структуру папок и уникальные имена
d1Rya[i] = MakeTH1('I', Form("Rya/Ryas_%2.2d",i)
, Form("R x amp CH %2.2d",i), 8192, 0., 8192.);
d1Ryt[i] = MakeTH1('I', Form("Ryt/Ryts_%2.2d",i)
, Form("R x time CH %2.2d",i), 4096, 1, 4096);
d1Ryc[i] = MakeTH1('D', Form("Ryc/Rycs_%2.2d",i)
, Form("R x energy MeV %2.2d",i), 4000, -2., 98.);
//
d1AnSa[i] = MakeTH1('I', Form("AnSa/AnSas_%2.2d",i) //указывать не повторяющуюся структуру папок и уникальные имена
d1Rva[i] = MakeTH1('I', Form("Rva/Rvas_%2.2d",i)
, Form("R veto amp CH %2.2d",i), 8192, 0., 8192.);
d1Rvt[i] = MakeTH1('I', Form("Rvt/Rvts_%2.2d",i)
, Form("R veto time CH %2.2d",i), 4096, 1, 4096);
d1Rvc[i] = MakeTH1('D', Form("Rvc/Rvcs_%2.2d",i)
, Form("R veto energy MeV %2.2d",i), 4000, -2., 98.);
//
d1AnSa[i] = MakeTH1('I', Form("AnSa/AnSas_%2.2d",i)
, Form("Sect amp CH %2.2d",i), 8192, 0., 8192.);
d1AnSt[i] = MakeTH1('I', Form("AnSt/AnSts_%2.2d",i)
, Form("Sect time CH %2.2d",i), 4096, 1, 4096);
d1AnSc[i] = MakeTH1('D', Form("AnSc/AnScs_%2.2d",i)
, Form("Sect energy MeV %2.2d",i), 4000, -2., 98.);
//
d1AnRa[i] = MakeTH1('I', Form("AnRa/AnRas_%2.2d",i) //указывать не повторяющуюся структуру папок и уникальные имена
d1AnRa[i] = MakeTH1('I', Form("AnRa/AnRas_%2.2d",i)
, Form("Ring amp CH %2.2d",i), 8192, 0., 8192.);
d1AnRt[i] = MakeTH1('I', Form("AnRt/AnRts_%2.2d",i)
, Form("Ring time CH %2.2d",i), 4096, 1, 4096);
d1AnRc[i] = MakeTH1('D', Form("AnRc/AnRcs_%2.2d",i)
, Form("Ring energy MeV %2.2d",i), 4000, -2., 98.);
d1AnVa[i] = MakeTH1('I', Form("AnVa/AnVas_%2.2d",i)
, Form("Veto amp CH %2.2d",i), 8192, 0., 8192.);
d1AnVt[i] = MakeTH1('I', Form("AnVt/AnVts_%2.2d",i)
, Form("Veto time CH %2.2d",i), 4096, 1, 4096);
d1AnVc[i] = MakeTH1('D', Form("AnVc/AnVcs_%2.2d",i)
, Form("Veto energy MeV %2.2d",i), 4000, -2., 98.);
//
}
for (int i = 0; i < 16; ++i) {
d1Lcsa[i] = MakeTH1('I', Form("Lcsa/Lcsa_%02d",i),
Form("Lcs amp CH %02d",i), 4096, 0., 4096);
Form("Lcs amp CH %02d",i), 8192, 0., 8192);
d1Lcst[i] = MakeTH1('I', Form("Lcst/Lcst_%02d",i)
, Form("Lcs time CH %02d",i), 4096, 1., 4095);
d1Lcsc[i] = MakeTH1('D', Form("Lcsc/Lcsc_%02d",i),
Form("Lcs energy MeV %2.2d",i), 4000, -2., 98.);
Form("Lcs energy MeV %2.2d",i), 4000, 0., 8192.);
}
for (int i = 0; i < 16; ++i) {
d1Rcsa[i] = MakeTH1('I', Form("Rcsa/Rcsa_%02d",i),
Form("Rcs amp CH %02d",i), 4096, 0., 4096);
Form("Rcs amp CH %02d",i), 8192, 0., 8192);
d1Rcst[i] = MakeTH1('I', Form("Rcst/Rcst_%02d",i)
, Form("Rcs time CH %02d",i), 4096, 1., 4095);
d1Rcsc[i] = MakeTH1('D', Form("Rcsc/Rcsc_%02d",i),
......@@ -177,13 +196,33 @@ TNeProc::TNeProc(const char* name) : TGo4EventProcessor(name)
}
for (int i = 0; i < 16; ++i) {
d1R20a[i] = MakeTH1('I', Form("R20a/R20a_%02d",i),
Form("R20 amp CH %02d",i), 4096, 0., 4096);
Form("R20 amp CH %02d",i), 8192, 0., 8192);
d1R20t[i] = MakeTH1('I', Form("R20t/R20t_%02d",i)
, Form("R20 time CH %02d",i), 4096, 1., 4095);
d1R20c[i] = MakeTH1('D', Form("R20c/R20c_%02d",i),
Form("R20 energy MeV %2.2d",i), 4000, -2., 98.);
}
d2ede = MakeTH2('D', "ede", "ede", 200,0.,20.,200,0.,20.,"e","de");
d2LyLcs1 = MakeTH2('D', "2D/LyLcs1", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs2 = MakeTH2('D', "2D/LyLcs2", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs3 = MakeTH2('D', "2D/LyLcs3", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs4 = MakeTH2('D', "2D/LyLcs4", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs5 = MakeTH2('D', "2D/LyLcs5", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs6 = MakeTH2('D', "2D/LyLcs6", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs7 = MakeTH2('D', "2D/LyLcs7", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs8 = MakeTH2('D', "2D/LyLcs8", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs9 = MakeTH2('D', "2D/LyLcs9", "Lx vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs10 = MakeTH2('D', "2D/LyLcs10", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs11 = MakeTH2('D', "2D/LyLcs11", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs12 = MakeTH2('D', "2D/LyLcs12", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs13 = MakeTH2('D', "2D/LyLcs13", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs14 = MakeTH2('D', "2D/LyLcs14", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs15 = MakeTH2('D', "2D/LyLcs15", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
d2LyLcs16 = MakeTH2('D', "2D/LyLcs16", "Ly vs Lcs", 200,0.,8200.,200,0.,80.,"Lcs","Ly");
//
d1F3[0] = MakeTH1('I',"F3/QDC/F3_QDC_L","F3 QDC left", 4096, 0., 4095.);
d1F3[1] = MakeTH1('I',"F3/QDC/F3_QDC_R","F3 QDC right", 4096, 0., 4095.);
......@@ -217,6 +256,7 @@ TNeProc::TNeProc(const char* name) : TGo4EventProcessor(name)
//-----MWPC-----
d2XY1 = MakeTH2('I', "MWPC/XY1", "WIRES 1 Y1 vs X1", 32,0.,32.,32,0.,32.,"X1","Y1");
d2XY2 = MakeTH2('I', "MWPC/XY2", "WIRES 2 Y2 vs X2", 32,0.,32.,32,0.,32.,"X2","Y2");
d2XYbeamOnTar = MakeTH2('f', "MWPC/XY_beamOnTar", "beam profile at target", 200,-50.0,50.0,200,-50,50,"X(Z=-70.0mm)","Y(Z=-70.0mm)");
d1X1 = MakeTH1('I', "MWPC/X1","WIRES X1",32,0.,32.);
d1Y1 = MakeTH1('I', "MWPC/Y1","WIRES Y1",32,0.,32.);
......@@ -248,7 +288,7 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
{
// called by framework from TNeEvent to fill it
// TGo4Analysis *go4 = TGo4Analysis::Instance();
TGo4Analysis *go4 = TGo4Analysis::Instance();
TNeEvent* NeEvent = (TNeEvent*) target;
target->Clear();
......@@ -260,9 +300,10 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
return false;
}
bool dooutput = false;
// double value1, value2;
// double eLxa,eLya;
// int mLxa,mLya;
double value1, value2;
double eLxa = 0,eLya = 0, eRxa = 0, eRya = 0, eR20a = 0, eRva = 0, eAnSa = 0, eAnRa = 0, eRcsa = 0, eLcsa = 0;
int mLxa = 0, mLya = 0, mRxa = 0, mRya = 0, mR20a = 0, mRva = 0, mAnSa = 0, mAnRa = 0, mRcsa = 0, mLcsa = 0;
int nLxa = -1, nLya = -1, nRxa = -1, nRya = -1, nR20a = -1, nRva = -1, nAnSa = -1, nAnRa = -1, nRcsa = -1, nLcsa = -1;
//
int num = 0, sz = 0;
int i, n, geo, typ, header,caen_header, module_id=0, crate=-1, cnt, chnl, evnt /*,gch*/;
......@@ -353,7 +394,7 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
case 9: // mADC: amp Lx & Ly
chnl= (idata[n] >> 16)&0x1f;
adc = idata[n] & 0xfff; //only 12 bits
adc = idata[n] & 0x1fff; //only 12 bits
if(chnl<16) {
NeEvent->Lxa[chnl] = adc;
} else {
......@@ -373,7 +414,7 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
case 11: // mADC: amp Rx & Ry
chnl= (idata[n] >> 16)&0x1f;
adc = idata[n] & 0xfff; //only 12 bits
adc = idata[n] & 0x1fff; //only 12 bits
if(chnl<16) {
NeEvent->Rxa[chnl] = adc;
} else {
......@@ -393,7 +434,7 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
case 13: // mADC: amp An S & R
chnl= (idata[n] >> 16)&0x1f;
adc = idata[n] & 0xfff; //only 12 bits
adc = idata[n] & 0x1fff; //only 12 bits
if(chnl<16) {
NeEvent->AnSa[chnl] = adc;
} else {
......@@ -413,7 +454,7 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
case 15: // mADC: amp R20 & AnV
chnl= (idata[n] >> 16)&0x1f;
adc = idata[n] & 0xfff; //only 12 bits
adc = idata[n] & 0x1fff; //only 12 bits
if(chnl<16) {
NeEvent->AnVa[chnl] = adc;
} else {
......@@ -421,10 +462,7 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
}
break;
}
break;
case 16: // mTDC: time Rcs & Lcs
case 2: // mTDC: time Rcs & Lcs
chnl= (idata[n] >> 16)&0x3f;
adc = idata[n] & 0xfff; //only 12 bits
if(chnl<16) {
......@@ -434,9 +472,9 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
}
break;
case 17: // mADC: amp Rcs & Lcs
case 3: // mADC: amp Rcs & Lcs
chnl= (idata[n] >> 16)&0x1f;
adc = idata[n] & 0xfff; //only 12 bits
adc = idata[n] & 0x1fff; //only 12 bits
if(chnl<16) {
NeEvent->Lcsa[chnl] = adc;
} else {
......@@ -444,6 +482,21 @@ Bool_t TNeProc::BuildEvent(TGo4EventElement* target)
}
break;
case 4: // mTDC: time Rvt
chnl= (idata[n] >> 16)&0x3f;
adc = idata[n] & 0xfff; //only 12 bits
if(chnl<16) NeEvent->Rvt[chnl] = adc;
break;
case 5: // mADC: amp Ra
chnl= (idata[n] >> 16)&0x1f;
adc = idata[n] & 0x1fff; //only 12 bits
if(chnl<16) NeEvent->Rva[chnl] = adc;
break;
}
break;
case 3: //-----MESYTEC footer-----
m_bool=false;
break;
......@@ -762,43 +815,65 @@ ToF = (tf5 - tf3)*0.015625 + 68.475;
}
}
for (int n=0; n<16; ++n) {
if(NeEvent->Lxt[n]>0){
// if(NeEvent->Lxt[n]>0)
{
d1Lxa[n]->Fill(NeEvent->Lxa[n]);
d1Lxt[n]->Fill(NeEvent->Lxt[n]);
}
if(NeEvent->Lyt[n]>0){
// if(NeEvent->Lyt[n]>0)
{
d1Lya[n]->Fill(NeEvent->Lya[n]);
d1Lyt[n]->Fill(NeEvent->Lyt[n]);
}
if(NeEvent->Rxt[n]>0){
// if(NeEvent->Rxt[n]>0)
{
d1Rxa[n]->Fill(NeEvent->Rxa[n]);
d1Rxt[n]->Fill(NeEvent->Rxt[n]);
}
if(NeEvent->Ryt[n]>0){
// if(NeEvent->Ryt[n]>0)
{
d1Rya[n]->Fill(NeEvent->Rya[n]);
d1Ryt[n]->Fill(NeEvent->Ryt[n]);
}
if(NeEvent->R20t[n]>0){
// if(NeEvent->R20t[n]>0)
{
d1R20a[n]->Fill(NeEvent->R20a[n]);
d1R20t[n]->Fill(NeEvent->R20t[n]);
}
if(NeEvent->AnSt[n]>0){
// if(NeEvent->AnSt[n]>0)
{
d1AnSa[n]->Fill(NeEvent->AnSa[n]);
d1AnSt[n]->Fill(NeEvent->AnSt[n]);
}
if(NeEvent->AnRt[n]>0){
// if(NeEvent->AnRt[n]>0)
{
d1AnRa[n]->Fill(NeEvent->AnRa[n]);
d1AnRt[n]->Fill(NeEvent->AnRt[n]);
}
if(NeEvent->Rcst[n]>0){
// if(NeEvent->AnVt[n]>0)
{
d1AnVa[n]->Fill(NeEvent->AnVa[n]);
d1AnVt[n]->Fill(NeEvent->AnVt[n]);
}
// if(NeEvent->Rcst[n]>0)
{
d1Rcsa[n]->Fill(NeEvent->Rcsa[n]);
d1Rcst[n]->Fill(NeEvent->Rcst[n]);
}
if(NeEvent->Lcst[n]>0){
// if(NeEvent->Lcst[n]>0)
{
d1Lcsa[n]->Fill(NeEvent->Lcsa[n]);
d1Lcst[n]->Fill(NeEvent->Lcst[n]);
}
// if(NeEvent->Rvt[n]>0)
{
d1Rva[n]->Fill(NeEvent->Rva[n]);
d1Rvt[n]->Fill(NeEvent->Rvt[n]);
}
}
......@@ -816,18 +891,274 @@ ToF = (tf5 - tf3)*0.015625 + 68.475;
if(NeEvent->tp1s[n] > 0)
d1_p1s[n]->Fill(NeEvent->tp1s[n]);
}
*/
//______ calibrated spectra __________________________
for (n=0; n<16; ++n)
{
value1 = (double)NeEvent->p1s[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_p1s->Energy(value1, n);
d1_Lxa_c[n]->Fill(value2);
mp1s++;
ep1s = value2;
if(NeEvent->Lxt[n]>0){
value1 = (double)NeEvent->Lxa[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_Lxa->Energy(value1, n);
d1Lxc[n]->Fill(value2);
mLxa++;
eLxa = value2;
}
if(NeEvent->Lyt[n]>0){
value1 = (double)NeEvent->Lya[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_Lya->Energy(value1, n);
d1Lyc[n]->Fill(value2);
mLya++;
eLya = value2;
}
if(NeEvent->Rxt[n]>0){
value1 = (double)NeEvent->Rxa[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_Rxa->Energy(value1, n);
d1Rxc[n]->Fill(value2);
mRxa++;
eRxa = value2;
}
if(NeEvent->Ryt[n]>0){
value1 = (double)NeEvent->Rya[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_Rya->Energy(value1, n);
d1Ryc[n]->Fill(value2);
mRya++;
eRya = value2;
}
if(NeEvent->R20t[n]>0){
value1 = (double)NeEvent->R20a[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_R20a->Energy(value1, n);
d1R20c[n]->Fill(value2);
mR20a++;
eR20a = value2;
}
if(NeEvent->Rvt[n]>0){
value1 = (double)NeEvent->Rva[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_Rva->Energy(value1, n);
d1Rvc[n]->Fill(value2);
mRva++;
eRva = value2;
}
if(NeEvent->AnSt[n]>0){
value1 = (double)NeEvent->AnSa[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_AnSa->Energy(value1, n);
d1AnSc[n]->Fill(value2);
mAnSa++;
eAnSa = value2;
}
if(NeEvent->AnRt[n]>0){
value1 = (double)NeEvent->AnRa[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_AnRa->Energy(value1, n);
d1AnRc[n]->Fill(value2);
mAnRa++;
eAnRa = value2;
}
if(NeEvent->Rcst[n]>0){
value1 = (double)NeEvent->Rcsa[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_Rcsa->Energy(value1, n);
d1Rcsc[n]->Fill(value2);
mRcsa++;
eRcsa = value2;
}
if(NeEvent->Lcst[n]>0 && NeEvent->Lcsa[n]>2500){
value1 = (double)NeEvent->Lcsa[n] + gRandom->Uniform();
value2 = ((TNeAnalysis*)go4)->p_Lcsa->Energy(value1, n);
d1Lcsc[n]->Fill(value2);
mLcsa++;
nLcsa = n;
eLcsa = value2;
}
}
if(mRya==1 && mR20a==1)
{
d2ede->Fill(eRya,eR20a);
}
if(mLya==1 && mLcsa==1)
{
if(nLcsa == 0) d2LyLcs1->Fill(eLcsa,eLya);
else if(nLcsa == 1) d2LyLcs2->Fill(eLcsa,eLya);
else if(nLcsa == 2) d2LyLcs3->Fill(eLcsa,eLya);
else if(nLcsa == 3) d2LyLcs4->Fill(eLcsa,eLya);
else if(nLcsa == 4) d2LyLcs5->Fill(eLcsa,eLya);
else if(nLcsa == 5) d2LyLcs6->Fill(eLcsa,eLya);
else if(nLcsa == 6) d2LyLcs7->Fill(eLcsa,eLya);
else if(nLcsa == 7) d2LyLcs8->Fill(eLcsa,eLya);
else if(nLcsa == 8) d2LyLcs9->Fill(eLcsa,eLya);
else if(nLcsa == 9) d2LyLcs10->Fill(eLcsa,eLya);
else if(nLcsa == 10) d2LyLcs11->Fill(eLcsa,eLya);
else if(nLcsa == 11) d2LyLcs12->Fill(eLcsa,eLya);
else if(nLcsa == 12) d2LyLcs13->Fill(eLcsa,eLya);
else if(nLcsa == 13) d2LyLcs14->Fill(eLcsa,eLya);
else if(nLcsa == 14) d2LyLcs15->Fill(eLcsa,eLya);
else if(nLcsa == 15) d2LyLcs16->Fill(eLcsa,eLya);
}
//---------------------start of MWPC----------------------------------------
// if(NeEvent->trigger == 1 && ToF > 165 && ToF <185)
// {
for(i = 0; i < NeEvent->nx1; ++i) d1X1->Fill(NeEvent->x1[i]);
for(i = 0; i < NeEvent->ny1; ++i) d1Y1->Fill(NeEvent->y1[i]);
for(i = 0; i < NeEvent->nx2; ++i) d1X2->Fill(NeEvent->x2[i]);
for(i = 0; i < NeEvent->ny2; ++i) d1Y2->Fill(NeEvent->y2[i]);
// }
if(NeEvent->nx1 > 0) d1NX1->Fill(NeEvent->nx1);
if(NeEvent->ny1 > 0) d1NY1->Fill(NeEvent->ny1);
if(NeEvent->nx2 > 0) d1NX2->Fill(NeEvent->nx2);
if(NeEvent->ny2 > 0) d1NY2->Fill(NeEvent->ny2);
if(NeEvent->nx1 > 0)
d1CX1->Fill(GetClasters(NeEvent->nx1, NeEvent->x1));
if(NeEvent->ny1 > 0)
d1CY1->Fill(GetClasters(NeEvent->ny1, NeEvent->y1));
if(NeEvent->nx2 > 0)
d1CX2->Fill(GetClasters(NeEvent->nx2, NeEvent->x2));
if(NeEvent->ny2 > 0)
d1CY2->Fill(GetClasters(NeEvent->ny2, NeEvent->y2));
// Li9
// if (NeEvent->nx1*NeEvent->ny1*NeEvent->nx2*NeEvent->ny2 == 1 && 162<ToF && ToF<170 && 1000<af3 && af3<1700)
// He6
// if (NeEvent->nx1*NeEvent->ny1*NeEvent->nx2*NeEvent->ny2 == 1 && ToF>158 && ToF<168 && af3>500 && af3<1100)
if (NeEvent->nx1*NeEvent->ny1*NeEvent->nx2*NeEvent->ny2 == 1)
{
mwz1 = -816;
mwx1 = (16 - NeEvent->x1[0] + gRandom->Uniform())*1.25;
mwy1 = (NeEvent->y1[0] - 16 + gRandom->Uniform())*1.25;
mwz2 = -270;
mwx2 = (16 - NeEvent->x2[0] + gRandom->Uniform())*1.25;
mwy2 = (NeEvent->y2[0] - 16 + gRandom->Uniform())*1.25;
zbt = -70;
xbt = mwx2 + (mwx2 - mwx1)*(zbt-mwz2)/(mwz2-mwz1);
ybt = mwy2 + (mwy2 - mwy1)*(zbt-mwz2)/(mwz2-mwz1);
// if (NeEvent->trigger == 1)
{
d2XYbeamOnTar->Fill(xbt, ybt);
}
/*if (NeEvent->trigger == 2, NeEvent->trigger == 3, NeEvent->trigger == 4, NeEvent->trigger == 5)
{
d2XYbeamOnTar_trig2->Fill(xbt, ybt);
}
if (NeEvent->tF6[0]*NeEvent->tF6[1]*NeEvent->tF6[2]*NeEvent->tF6[3]>0)
{
d2XYbeamOnTar_F6->Fill(xbt, ybt);
}
*/
}
//
if(NeEvent->nx1 > 0 && NeEvent->ny1 > 0 && (GetClasters(NeEvent->nx1, NeEvent->x1)) == 1 && (GetClasters(NeEvent->ny1, NeEvent->y1)) == 1)
{
d2XY1->Fill(NeEvent->x1[0],NeEvent->y1[0]);
// if(NeEvent->tF6[2] >100) d2XY1beam->Fill(NeEvent->x1[0],NeEvent->y1[0]);
// if(NeEvent->trigger == 1) d2XY1beam->Fill(NeEvent->x1[0],NeEvent->y1[0]);
// else d2XY1->Fill(NeEvent->x1[0],NeEvent->y1[0]);
}
if(NeEvent->nx2 > 0 && NeEvent->ny2 > 0 && (GetClasters(NeEvent->nx2, NeEvent->x2)) == 1 && (GetClasters(NeEvent->ny2, NeEvent->y2)) == 1)
{
d2XY2->Fill(NeEvent->x2[0],NeEvent->y2[0]);
// if(NeEvent->tF6[2] >100) d2XY2beam->Fill(NeEvent->x2[0],NeEvent->y2[0]);
// if(NeEvent->trigger == 1) d2XY2beam->Fill(NeEvent->x2[0],NeEvent->y2[0]);
// else d2XY2->Fill(NeEvent->x2[0],NeEvent->y2[0]);
}
/*
if(NeEvent->nx1 == 1 && NeEvent->ny1 == 1 && NeEvent->trigger == 1)
{
d2XY1->Fill(NeEvent->x1[0],NeEvent->y1[0]);
}
if(NeEvent->nx2 == 1 && NeEvent->ny2 == 1 && NeEvent->trigger == 1)
{
d2XY2->Fill(NeEvent->x2[0],NeEvent->y2[0]);
}
*/
//---Pictures---
/*
if(NeEvent->nx1 == 1 && NeEvent->nx2 == 1)
{
ya1 = -0.5*pya1 + dx1 + (pya1 / 32.)*NeEvent->x1[0];
ya2 = -0.5*pya2 + dx2 + (pya2 / 32.)*NeEvent->x2[0];
ya3 = ya1 + dtx + (ya2 - ya1) * (pxa3 - pxa1) / (pxa2 - pxa1);
if(ya3 > -0.5*pya1 && ya3 < 0.5*pya1)
{
if(((TNeAnalysis*)go4)->pMTracks->mtracksx != ntrackx)
{
ntrackx = ((TNeAnalysis*)go4)->pMTracks->mtracksx;
fTracksX->Set(3*ntrackx);
itrackAX = 0;
for(i = 0; i < fTracksX->GetN(); ++i)
fTracksX->SetPoint(i,pxb1,0.);
}
if(itrackAX == fTracksX->GetN()) itrackAX = 0;
fTracksX->SetPoint(itrackAX++,pxa1,ya1);
fTracksX->SetPoint(itrackAX++,pxa3,ya3);
fTracksX->SetPoint(itrackAX++,pxa1,ya1);
}
}
if(NeEvent->ny1 == 1 && NeEvent->ny2 == 1)
{
yb1 = -0.5*pyb1 + dy1+ (pyb1 / 32.)*NeEvent->y1[0];
yb2 = -0.5*pyb2 + dy2+ (pyb2 / 32.)*NeEvent->y2[0];
yb3 = yb1 + dty + (yb2 - yb1) * (pxb3 - pxb1) / (pxb2 - pxb1);
if(yb3 > -0.5*pyb1 && yb3 < 0.5*pyb1)
{
if(((TNeAnalysis*)go4)->pMTracks->mtracksy != ntracky)
{
ntracky = ((TNeAnalysis*)go4)->pMTracks->mtracksy;
fTracksY->Set(3*ntracky);
itrackBX = 0;
for(i = 0; i < fTracksY->GetN(); ++i)
fTracksY->SetPoint(i,pxb1,0.);
}
if(itrackBX == fTracksY->GetN()) itrackBX = 0;
fTracksY->SetPoint(itrackBX++,pxb1,yb1);
fTracksY->SetPoint(itrackBX++,pxb3,yb3);
fTracksY->SetPoint(itrackBX++,pxb1,yb1);
}
}
if(NeEvent->nx1 == 1 && NeEvent->nx2 == 1 && NeEvent->ny1 == 1 && NeEvent->ny2 == 1 )
{
d2Target->Fill(ya3,yb3);
}*/
//--------------------end of MWPC----------------------------------------
NeEvent->SetValid(kTRUE); // to store
return kTRUE;
} //-------------------------------------------------------------------
......
TNeProc.h
View file @ 283f92c7
......@@ -38,6 +38,9 @@ private:
TH1 *d1Rya[16];
TH1 *d1Ryt[16];
TH1 *d1Ryc[16];
TH1 *d1Rva[16];
TH1 *d1Rvt[16];
TH1 *d1Rvc[16];
TH1 *d1AnSa[16];
TH1 *d1AnSt[16];
TH1 *d1AnSc[16];
......@@ -60,6 +63,25 @@ private:
TH1 *d1Lcsa[16];
TH1 *d1Lcst[16];
TH1 *d1Lcsc[16];
//______________ 2 dim
TH2 *d2ede;
TH2 *d2LyLcs1;
TH2 *d2LyLcs2;
TH2 *d2LyLcs3;
TH2 *d2LyLcs4;
TH2 *d2LyLcs5;
TH2 *d2LyLcs6;
TH2 *d2LyLcs7;
TH2 *d2LyLcs8;
TH2 *d2LyLcs9;
TH2 *d2LyLcs10;
TH2 *d2LyLcs11;
TH2 *d2LyLcs12;
TH2 *d2LyLcs13;
TH2 *d2LyLcs14;
TH2 *d2LyLcs15;
TH2 *d2LyLcs16;
//-----clb-----
TH1 *d1sq1_x_c[16];
......@@ -94,6 +116,12 @@ private:
TH1 *d1CY1;
TH1 *d1CX2;
TH1 *d1CY2;
TH2 *d2XYbeamOnTar;
double mwz1, mwx1, mwy1;
double mwz2, mwx2, mwy2;
double zbt, xbt, ybt;
ClassDef(TNeProc,1)
};
......
clb/anra.clb 0 → 100644
View file @ 283f92c7
2
16
-0.9921 0.006043 0
-1.295 0.006109 0
-1.034 0.006167 0
-1.226 0.005974 0
-1.087 0.006062 0
-1.293 0.006128 0
-1.117 0.006241 0
-1.569 0.006415 0
-0.9381 0.005942 0
-1.331 0.005939 0
-1.029 0.00602 0
-1.252 0.006154 0
-0.9552 0.006104 0
-1.25 0.005964 0
-0.9905 0.005939 0
-1.24 0.006108 0
clb/ansa.clb 0 → 100644
View file @ 283f92c7
2
16
-1.042 0.006175 0
-1.175 0.005682 0
-0.956 0.006046 0
-1.212 0.00593 0
-1.005 0.006027 0
-1.249 0.005907 0
-1.035 0.005957 0
-1.215 0.005924 0
-1.064 0.006153 0
-1.337 0.00598 0
-1.073 0.005888 0
-1.235 0.005878 0
-1.054 0.006009 0
-1.005 0.005548 0
-1.112 0.006079 0
-1.28 0.005843 0
clb/anva.clb 0 → 100644
View file @ 283f92c7
2
16
-0.2076 0.005729 0
-0.1873 0.00575 0
-0.1988 0.005828 0
-0.197 0.005838 0
-0.2144 0.00588 0
-0.1965 0.005734 0
-0.215 0.005697 0
-0.1958 0.005678 0
-0.2061 0.005712 0
-0.2116 0.005771 0
-0.2041 0.005894 0
-0.2219 0.005869 0
-0.2133 0.005625 0
-0.2355 0.005817 0
-0.2177 0.005723 0
-0.2396 0.005987 0
clb/lxa.clb 0 → 100644
View file @ 283f92c7
2
16
-0.331 0.01068 0
-0.3837 0.01094 0
-0.3164 0.01119 0
-0.3824 0.01105 0
-0.3201 0.01086 0
-0.3969 0.01076 0
-0.3262 0.01092 0
-0.4261 0.01109 0
-0.4124 0.01064 0
-0.4005 0.01113 0
-0.4111 0.01112 0
-0.3407 0.01073 0
-0.4172 0.01071 0
-0.3712 0.01063 1
-0.4475 0.01099 0
-0.362 0.01077 0
clb/lya.clb 0 → 100644
View file @ 283f92c7
2
16
-0.4061 0.01084 0
-0.3753 0.01097 0
-0.4124 0.01095 0
-0.374 0.01086 0
-0.4033 0.01121 0
-0.3717 0.01077 0
-0.4068 0.01111 0
-0.3854 0.01101 0
-0.346 0.01082 0
-0.4533 0.01102 0
-0.3497 0.01072 0
-0.451 0.01108 0
-0.3556 0.01093 0
-0.4478 0.01118 0
-0.364 0.01093 0
-0.4304 0.01064 0
clb/p1s.clb deleted 100644 → 0
View file @ 756eb644
2
16
-0.08733 0.004911 0
-0.08709 0.004893 0
-0.08096 0.004951 0
-0.07999 0.004812 0
-0.08556 0.004968 0
-0.08162 0.004784 0
-0.08238 0.004891 0
-0.0817 0.004798 0
-0.09605 0.005018 0
-0.06971 0.004917 0
-0.09164 0.004982 0
-0.07341 0.005013 0
-0.08384 0.005046 0
-0.07682 0.004999 0
-0.08593 0.005024 0
-0.07198 0.004797 0
clb/r20a.clb 0 → 100644
View file @ 283f92c7
2
16
-1.74 0.0113 0
-2.397 0.01151 0
-1.672 0.01104 0
-2.264 0.01115 0
-1.733 0.01108 0
-2.363 0.01142 0
-1.671 0.01109 0
-2.288 0.01127 0
-1.639 0.01128 0
-2.422 0.01203 0
-1.685 0.01139 0
-2.207 0.01111 0
-1.72 0.01121 0
-2.242 0.0114 0
-1.707 0.01082 0
-2.294 0.01155 0
clb/rva.clb 0 → 100644
View file @ 283f92c7
2
16
-1.717 0.01114 0
-2.412 0.01145 0
-1.638 0.01085 0
-2.253 0.01104 0
-1.764 0.011 0
-2.403 0.01143 0
-1.654 0.01095 0
-2.305 0.01121 0
-1.63 0.01117 0
-2.329 0.01142 0
-1.659 0.01126 0
-2.202 0.01097 0
-1.698 0.01104 0
-2.253 0.0113 0
-1.699 0.01069 0
-2.305 0.01146 0
clb/rxa.clb 0 → 100644
View file @ 283f92c7
2
16
-0.4585 0.01085 0
-0.4486 0.01123 0
-0.4754 0.01117 0
-0.2406 0.01096 0
-0.481 0.01117 0
-0.495 0.01145 0
-0.492 0.01093 0
-0.5018 0.01105 0
-0.4927 0.01129 0
-0.4857 0.01117 0
-0.4859 0.01083 0
-0.4279 0.01095 0
-0.4331 0.01089 0
-0.4436 0.01121 0
-0.4329 0.0113 0
-0.4325 0.01129 0
clb/rya.clb 0 → 100644
View file @ 283f92c7
2
16
-0.477 0.01132 0
-0.4198 0.01142 0
-0.4798 0.01101 0
-0.3945 0.01124 0
-0.4799 0.01149 0
-0.3777 0.01094 0
-0.4694 0.01119 0
-0.4171 0.01127 0
-0.4655 0.01118 0
-0.4153 0.01118 0
-0.46 0.01105 0
-0.4163 0.01102 0
-0.4773 0.0109 0
-0.427 0.01087 0
-0.4686 0.01112 0
-0.4137 0.01091 0
clb/t1xa.clb deleted 100644 → 0
View file @ 756eb644
2
16
-0.08733 0.004911 0
-0.08709 0.004893 0
-0.08096 0.004951 0
-0.07999 0.004812 0
-0.08556 0.004968 0
-0.08162 0.004784 0
-0.08238 0.004891 0
-0.0817 0.004798 0
-0.09605 0.005018 0
-0.06971 0.004917 0
-0.09164 0.004982 0
-0.07341 0.005013 0
-0.08384 0.005046 0
-0.07682 0.004999 0
-0.08593 0.005024 0
-0.07198 0.004797 0
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment