Commit e37cda26 authored by Adrien Matta's avatar Adrien Matta
Browse files

* Adding files to S034

parent 980dc392
......@@ -158,21 +158,24 @@ void TNebulaPhysics::BuildPhysicalEvent() {
}
// Got everything, do the math
if(rawTup>0){
// cal Q Up and Down
calQup=aQu[ID]*(rawQup-bQu[ID]);
calQdown=aQd[ID]*(rawQdown-bQd[ID]);
// cal T
// average value of Up and Down
calQ=sqrt(calQup*calQdown);
// cal T Up
calTup=aTu[ID]*rawTup+bTu[ID];
// slew correction
calTup -= slwTu[ID]/sqrt(rawQup-bQu[ID]);
// cal T
// cal T Down
calTdown=aTd[ID]*rawTdown+bTd[ID];
// slew correction
calTdown -= slwTd[ID]/sqrt(rawQdown-bQd[ID]);
// average value of Up and Down
calQ=sqrt(calQup*calQdown);
if(calQ>threshold){
calT= (calTdown+calTup)*0.5+avgT0[ID]+Cal->GetPedestal("NEBULA_T_ID"+NPL::itoa(ID));
Y=(calTdown-calTup)*DTa[ID]+DTb[ID]+Cal->GetPedestal("NEBULA_Y_ID"+NPL::itoa(ID));
......@@ -184,7 +187,7 @@ void TNebulaPhysics::BuildPhysicalEvent() {
PosX.push_back(PositionX[ID]);
PosZ.push_back(PositionZ[ID]);
if(ID<120)
if(ID<121)
IsVeto.push_back(0);
else
IsVeto.push_back(1);
......@@ -297,9 +300,11 @@ void TNebulaPhysics::WriteSpectra() {
///////////////////////////////////////////////////////////////////////////
void TNebulaPhysics::AddParameterToCalibrationManager() {
CalibrationManager* Cal = CalibrationManager::getInstance();
vector<double> standardO={0};
for (int i = 0; i < m_NumberOfBars; ++i) {
Cal->AddParameter("NEBULA_T_ID"+ NPL::itoa(i+1));
Cal->AddParameter("NEBULA_Y_ID"+ NPL::itoa(i+1));
Cal->AddParameter("NEBULA_T_ID"+ NPL::itoa(i+1),standardO);
Cal->AddParameter("NEBULA_Y_ID"+ NPL::itoa(i+1),standardO);
}
}
......
......@@ -80,13 +80,13 @@ void Analysis::TreatEvent(){
double TF13=-1;
double TF7=-1;
double TOF_F5F13=-1;
double TOF_F5F7=-1;
double TOF_F7F13=-1;
unsigned int sizeP = Plastic->FP.size();
for(unsigned int i = 0 ; i < sizeP ; i++){
if(Plastic->FP[i]==5){
TF5=Plastic->TSlew[i];
}
else if(Plastic->FP[i]==13){
else if(Plastic->FP[i]==13&&Plastic->ID[i]==4){// two plastic at F13 taking only one
TF13=Plastic->TSlew[i];
}
else if(Plastic->FP[i]==7){
......@@ -94,23 +94,26 @@ void Analysis::TreatEvent(){
}
}
if(TF5>0 && TF13>0){
TOF_F5F13=TF13-TF5;
double l = 117915-54917;
Beta_b=(l/TOF_F5F13)/NPUNITS::c_light;
if(TF7>0 && TF13>0){
// offset is adjusted to give the expected beta
TOF_F7F13=TF13-TF7+6.71626e+02;
static double LengthF7F13 = 117915-66409;
Beta_b=(LengthF7F13/TOF_F7F13)/NPUNITS::c_light;
// to find offset:
//Beta_b=TOF_F7F13-LengthF7F13/He8.GetVelocity();
}
// Samurai
if( FDC2->PosX>-1500 && FDC2->PosX<1000
// Samurai-Minos
if( Beta_b>0.5140 && Beta_b < 0.5165 // Correct Beta
&& Hodo->Charge.size()==1 && Hodo->Charge[0]>28 && Hodo->Charge[0]<42 && Hodo->Time[0]>58 && Hodo->Time[0]<68 // 6He in Hodo->cope
&& FDC2->PosX>-1500 && FDC2->PosX<1000
&& FDC2->PosY>-500 && FDC2->PosY<500
&& FDC0->PosX>-80 && FDC0->PosX<80
&& FDC0->PosY>-80 && FDC0->PosY<80 // both FDC ok
&& Minos->Tracks_P0.size()>0 ) { // p,pn or p,2p
&& Minos->Tracks_P0.size()==1 ) { // p,pn only
// Compute ThetaX and PhiY using Minos vertex and FDC0 X
// Check if both BDC are reconstructed
TVector3 BDC1=BDC->GetPos(1);
TVector3 BDC2=BDC->GetPos(2);
if( BDC1.Z()!=-10000 && BDC2.Z()!=-10000){
TVector3 Vertex,delta;
TVector3 P1 = Minos->Tracks_P0[0]+Minos->Tracks_Dir[0];
......@@ -156,8 +159,8 @@ void Analysis::TreatEvent(){
unsigned int first = Nebula->GetFirstHit();
TVector3 Pfirst = (Nebula->GetPos(first)-Vertex);
double L = Pfirst.Mag();
double TSBT= (Vertex.Z()+7377.56)/He8.GetVelocity();
TOF_n = Nebula->TOF[first]-TSBT;
double TSBT= (Vertex.Z()+7377.56)/(Beta_b*c_light);
TOF_n = Nebula->TOF[first]-TSBT-TF13;
Beta_n = (L/TOF_n)/NPUNITS::c_light;
LVn.SetVectM(TVector3(0,0,0),mn);
LVn.Boost(Beta_n*Pfirst.Unit());
......@@ -209,6 +212,7 @@ void Analysis::Clear(){
Beta_b=-1000;
X=Y=Z=-1000;
TOF_n=-1000;
Erel=-1000;
}
////////////////////////////////////////////////////////////////////////////////
......
<?xml version="1.0" encoding="UTF-8"?>
<dataroot>
<BigRIPSIC>
<ID>1</ID>
<NAME>F2IC</NAME>
<FPL>2</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>0</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>4</detectorid>
<geo0>0</geo0>
<ch0>0</ch0>
<geo1>0</geo1>
<ch1>1</ch1>
<geo2>0</geo2>
<ch2>2</ch2>
<geo3>0</geo3>
<ch3>3</ch3>
<geo4>0</geo4>
<ch4>4</ch4>
<geo5>0</geo5>
<ch5>5</ch5>
<geo6>0</geo6>
<ch6>6</ch6>
<geo7>0</geo7>
<ch7>7</ch7>
<geo8>0</geo8>
<ch8>8</ch8>
<geo9>0</geo9>
<ch9>9</ch9>
<geo10>-1</geo10>
<ch10>-1</ch10>
<geo11>-1</geo11>
<ch11>-1</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>2</ID>
<NAME>F3IC</NAME>
<FPL>3</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>0</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>4</detectorid>
<geo0>0</geo0>
<ch0>0</ch0>
<geo1>0</geo1>
<ch1>1</ch1>
<geo2>0</geo2>
<ch2>2</ch2>
<geo3>0</geo3>
<ch3>3</ch3>
<geo4>0</geo4>
<ch4>4</ch4>
<geo5>0</geo5>
<ch5>5</ch5>
<geo6>0</geo6>
<ch6>6</ch6>
<geo7>0</geo7>
<ch7>7</ch7>
<geo8>0</geo8>
<ch8>8</ch8>
<geo9>0</geo9>
<ch9>9</ch9>
<geo10>-1</geo10>
<ch10>-1</ch10>
<geo11>-1</geo11>
<ch11>-1</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>3</ID>
<NAME>F7IC</NAME>
<FPL>7</FPL>
<ch2mev_0>32.143</ch2mev_0>
<ch2mev_1>0.0564</ch2mev_1>
<zcoef_0>15.7363</zcoef_0>
<zcoef_1>-2.59258</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>4</detectorid>
<geo0>0</geo0>
<ch0>0</ch0>
<geo1>0</geo1>
<ch1>1</ch1>
<geo2>0</geo2>
<ch2>2</ch2>
<geo3>0</geo3>
<ch3>3</ch3>
<geo4>0</geo4>
<ch4>4</ch4>
<geo5>0</geo5>
<ch5>5</ch5>
<geo6>0</geo6>
<ch6>6</ch6>
<geo7>0</geo7>
<ch7>7</ch7>
<geo8>0</geo8>
<ch8>8</ch8>
<geo9>0</geo9>
<ch9>9</ch9>
<geo10>-1</geo10>
<ch10>-1</ch10>
<geo11>-1</geo11>
<ch11>-1</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>4</ID>
<NAME>F11IC</NAME>
<FPL>11</FPL>
<ch2mev_0>31.6</ch2mev_0>
<ch2mev_1>0.0849</ch2mev_1>
<zcoef_0>7.9065</zcoef_0>
<zcoef_1>-0.043</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>4</detectorid>
<geo0>0</geo0>
<ch0>0</ch0>
<geo1>0</geo1>
<ch1>1</ch1>
<geo2>0</geo2>
<ch2>2</ch2>
<geo3>0</geo3>
<ch3>3</ch3>
<geo4>0</geo4>
<ch4>4</ch4>
<geo5>0</geo5>
<ch5>5</ch5>
<geo6>0</geo6>
<ch6>6</ch6>
<geo7>0</geo7>
<ch7>7</ch7>
<geo8>0</geo8>
<ch8>8</ch8>
<geo9>0</geo9>
<ch9>9</ch9>
<geo10>-1</geo10>
<ch10>-1</ch10>
<geo11>-1</geo11>
<ch11>-1</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>5</ID>
<NAME>ICB</NAME>
<FPL>13</FPL>
<ch2mev_0>-1.51</ch2mev_0>
<ch2mev_1>0.03495</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>49.6</pedestal0>
<pedestal1>48.8</pedestal1>
<pedestal2>52.2</pedestal2>
<pedestal3>48.3</pedestal3>
<pedestal4>52.2</pedestal4>
<pedestal5>48</pedestal5>
<pedestal6>51.5</pedestal6>
<pedestal7>47.2</pedestal7>
<pedestal8>53.5</pedestal8>
<pedestal9>49.9</pedestal9>
<detectorid>25</detectorid>
<geo0>3</geo0>
<ch0>0</ch0>
<geo1>3</geo1>
<ch1>1</ch1>
<geo2>3</geo2>
<ch2>2</ch2>
<geo3>3</geo3>
<ch3>3</ch3>
<geo4>3</geo4>
<ch4>4</ch4>
<geo5>3</geo5>
<ch5>5</ch5>
<geo6>3</geo6>
<ch6>6</ch6>
<geo7>3</geo7>
<ch7>7</ch7>
<geo8>3</geo8>
<ch8>8</ch8>
<geo9>3</geo9>
<ch9>9</ch9>
<geo10>-1</geo10>
<ch10>-1</ch10>
<geo11>-1</geo11>
<ch11>-1</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>6</ID>
<NAME>ICF0</NAME>
<FPL>13</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>1</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>26</detectorid>
<geo0>3</geo0>
<ch0>0</ch0>
<geo1>3</geo1>
<ch1>4</ch1>
<geo2>3</geo2>
<ch2>8</ch2>
<geo3>3</geo3>
<ch3>12</ch3>
<geo4>3</geo4>
<ch4>16</ch4>
<geo5>3</geo5>
<ch5>20</ch5>
<geo6>3</geo6>
<ch6>24</ch6>
<geo7>3</geo7>
<ch7>28</ch7>
<geo8>4</geo8>
<ch8>16</ch8>
<geo9>4</geo9>
<ch9>20</ch9>
<geo10>4</geo10>
<ch10>24</ch10>
<geo11>4</geo11>
<ch11>28</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>7</ID>
<NAME>ICF1</NAME>
<FPL>13</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>1</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>26</detectorid>
<geo0>3</geo0>
<ch0>1</ch0>
<geo1>3</geo1>
<ch1>5</ch1>
<geo2>3</geo2>
<ch2>9</ch2>
<geo3>3</geo3>
<ch3>13</ch3>
<geo4>3</geo4>
<ch4>17</ch4>
<geo5>3</geo5>
<ch5>21</ch5>
<geo6>3</geo6>
<ch6>25</ch6>
<geo7>3</geo7>
<ch7>29</ch7>
<geo8>4</geo8>
<ch8>17</ch8>
<geo9>4</geo9>
<ch9>21</ch9>
<geo10>4</geo10>
<ch10>25</ch10>
<geo11>4</geo11>
<ch11>29</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>8</ID>
<NAME>ICF2</NAME>
<FPL>13</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>1</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>26</detectorid>
<geo0>3</geo0>
<ch0>2</ch0>
<geo1>3</geo1>
<ch1>6</ch1>
<geo2>3</geo2>
<ch2>10</ch2>
<geo3>3</geo3>
<ch3>14</ch3>
<geo4>3</geo4>
<ch4>18</ch4>
<geo5>3</geo5>
<ch5>22</ch5>
<geo6>3</geo6>
<ch6>26</ch6>
<geo7>3</geo7>
<ch7>30</ch7>
<geo8>4</geo8>
<ch8>18</ch8>
<geo9>4</geo9>
<ch9>22</ch9>
<geo10>4</geo10>
<ch10>26</ch10>
<geo11>4</geo11>
<ch11>30</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>9</ID>
<NAME>ICF3</NAME>
<FPL>13</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>1</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>4866</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>26</detectorid>
<geo0>3</geo0>
<ch0>3</ch0>
<geo1>3</geo1>
<ch1>7</ch1>
<geo2>3</geo2>
<ch2>11</ch2>
<geo3>3</geo3>
<ch3>15</ch3>
<geo4>3</geo4>
<ch4>19</ch4>
<geo5>3</geo5>
<ch5>23</ch5>
<geo6>3</geo6>
<ch6>27</ch6>
<geo7>3</geo7>
<ch7>31</ch7>
<geo8>4</geo8>
<ch8>19</ch8>
<geo9>4</geo9>
<ch9>23</ch9>
<geo10>4</geo10>
<ch10>27</ch10>
<geo11>4</geo11>
<ch11>31</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>10</ID>
<NAME>ICB_dummy_1</NAME>
<FPL>13</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>0</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>0</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>25</detectorid>
<geo0>3</geo0>
<ch0>10</ch0>
<geo1>3</geo1>
<ch1>11</ch1>
<geo2>3</geo2>
<ch2>12</ch2>
<geo3>3</geo3>
<ch3>13</ch3>
<geo4>3</geo4>
<ch4>14</ch4>
<geo5>3</geo5>
<ch5>15</ch5>
<geo6>3</geo6>
<ch6>16</ch6>
<geo7>3</geo7>
<ch7>17</ch7>
<geo8>3</geo8>
<ch8>18</ch8>
<geo9>3</geo9>
<ch9>19</ch9>
<geo10>3</geo10>
<ch10>20</ch10>
<geo11>3</geo11>
<ch11>21</ch11>
</BigRIPSIC>
<BigRIPSIC>
<ID>11</ID>
<NAME>ICB_dummy_2</NAME>
<FPL>13</FPL>
<ch2mev_0>0</ch2mev_0>
<ch2mev_1>0</ch2mev_1>
<zcoef_0>0</zcoef_0>
<zcoef_1>0</zcoef_1>
<zcoef_2>0</zcoef_2>
<ionpair>0</ionpair>
<pedestal0>0</pedestal0>
<pedestal1>0</pedestal1>
<pedestal2>0</pedestal2>
<pedestal3>0</pedestal3>
<pedestal4>0</pedestal4>
<pedestal5>0</pedestal5>
<pedestal6>0</pedestal6>
<pedestal7>0</pedestal7>
<pedestal8>0</pedestal8>
<pedestal9>0</pedestal9>
<detectorid>25</detectorid>
<geo0>3</geo0>
<ch0>22</ch0>
<geo1>3</geo1>
<ch1>23</ch1>
<geo2>3</geo2>
<ch2>24</ch2>
<geo3>3</geo3>
<ch3>25</ch3>
<geo4>3</geo4>
<ch4>26</ch4>
<geo5>3</geo5>
<ch5>27</ch5>
<geo6>3</geo6>
<ch6>28</ch6>
<geo7>3</geo7>
<ch7>29</ch7>
<geo8>3</geo8>
<ch8>30</ch8>
<geo9>3</geo9>