@@ -100,12 +100,12 @@ template <class T> void Zap(T*& t)
100100
101101// ============ namespace dragon Free Functions ============ //
102102
103- void dragon::MakeChains (const Int_t* runnumbers, Int_t nruns, Bool_t sonik, const char * format)
103+ void dragon::MakeChains (const Int_t* runnumbers, Int_t nruns, const char * format, Bool_t sonik )
104104{
105- MakeChains (" t" , runnumbers, nruns, sonik, format );
105+ MakeChains (" t" , runnumbers, nruns, format, sonik );
106106}
107107
108- void dragon::MakeChains (const char * prefix, const Int_t* runnumbers, Int_t nruns, Bool_t sonik, const char * format)
108+ void dragon::MakeChains (const char * prefix, const Int_t* runnumbers, Int_t nruns, const char * format, Bool_t sonik )
109109{
110110 // /
111111 // / \param [in] runnumbers Pointer to a valid array of desired run numbers to chain together.
@@ -152,20 +152,21 @@ void dragon::MakeChains(const char* prefix, const Int_t* runnumbers, Int_t nruns
152152 chain[6 ]->SetName (Form (" %s6" , prefix));
153153 chain[7 ]->SetName (Form (" %s7" , prefix));
154154 }
155+
155156 else {
156157 TChain* chain[] = {
157158 new TChain (Form (" %s0" , " t" ), " SONIK event." ),
158159 // new TChain(Form("%s1", "t"), "Head singles event."), Head DAQ not used with SONIK
159160 // new TChain(Form("%s2", "t"), "Head scaler event."),
160161 new TChain (Form (" %s3" , " t" ), " Tail singles event." ),
161162 new TChain (Form (" %s4" , " t" ), " Tail scaler event." ),
162- // new TChain(Form("%s5", "t"), "Coincidence event."), coincidence events not used with SONIK
163+ // new TChain(Form("%s5", "t"), "Coincidence event."), Coincidence events N/A for SONIK
163164 new TChain (Form (" %s20" , " t" ), " Epics event." ),
164165 new TChain (Form (" %s6" , " t" ), " Timestamp diagnostics." ),
165166 new TChain (Form (" %s7" , " t" ), " Glocal run parameters." )
166167 };
167168 Int_t nchains = sizeof (chain) / sizeof (TChain*);
168-
169+
169170 for (Int_t i=0 ; i< nruns; ++i) {
170171 char fname[4096 ];
171172 sprintf (fname, format, runnumbers[i]);
@@ -183,29 +184,29 @@ void dragon::MakeChains(const char* prefix, const Int_t* runnumbers, Int_t nruns
183184 }
184185
185186 chain[0 ]->SetName (Form (" %s0" , prefix));
186- // chain[1 ]->SetName(Form("%s1", prefix));
187- // chain[2 ]->SetName(Form("%s2", prefix));
187+ // chain[0 ]->SetName(Form("%s1", prefix));
188+ // chain[1 ]->SetName(Form("%s2", prefix));
188189 chain[1 ]->SetName (Form (" %s3" , prefix));
189190 chain[2 ]->SetName (Form (" %s4" , prefix));
190- // chain[5 ]->SetName(Form("%s5", prefix));
191+ // chain[4 ]->SetName(Form("%s5", prefix));
191192 chain[3 ]->SetName (Form (" %s20" , prefix));
192193 chain[4 ]->SetName (Form (" %s6" , prefix));
193- chain[5 ]->SetName (Form (" %s7" , prefix));
194+ chain[5 ]->SetName (Form (" %s7" , prefix));
194195 }
195196}
196197
197- void dragon::MakeChains (const char * prefix, const std::vector<Int_t>& runnumbers, Bool_t sonik, const char * format)
198+ void dragon::MakeChains (const char * prefix, const std::vector<Int_t>& runnumbers, const char * format, Bool_t sonik )
198199{
199- MakeChains (prefix, &runnumbers[0 ], runnumbers.size (), sonik, format );
200+ MakeChains (prefix, &runnumbers[0 ], runnumbers.size (), format, sonik );
200201}
201202
202- void dragon::MakeChains (const std::vector<Int_t>& runnumbers, Bool_t sonik, const char * format)
203+ void dragon::MakeChains (const std::vector<Int_t>& runnumbers, const char * format, Bool_t sonik )
203204{
204205 // /
205206 // / \param runnumbers vector of desired run numbers
206207 // / \param format same as the other version
207208 // /
208- MakeChains (&runnumbers[0 ], runnumbers.size (), sonik, format );
209+ MakeChains (&runnumbers[0 ], runnumbers.size (), format, sonik );
209210}
210211
211212void dragon::FriendChain (TChain* chain, const char * friend_name, const char * friend_alias,
@@ -2174,6 +2175,7 @@ UDouble_t dragon::StoppingPowerCalculator::CalculateEpsilon(TGraph** plot, UDoub
21742175
21752176 double md1err = fMd1 [i].GetErrLow () < fMd1 [i].GetErrHigh () ?
21762177 fMd1 [i].GetErrHigh () : fMd1 [i].GetErrLow ();
2178+
21772179 dens[i] = CalculateDensity (fPressures [i], UDouble_t (fTargetLength .GetNominal (), 0 ), fNmol , fTemp );
21782180 energy[i] = CalculateEnergy (fMd1 [i].GetNominal (), md1err, fBeamCharge , fBeamMass , fMd1Constant .GetNominal (), 0 );
21792181 }
0 commit comments