Changeset 19266
- Timestamp:
- 10/22/18 20:50:24 (6 years ago)
- Location:
- trunk/Mars/fact
- Files:
-
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Mars/fact/analysis/callisto_data.C
r18577 r19266 75 75 // Extraction type: Extract integral and half leading edge 76 76 77 const int type = maximum ? (MExtralgoSpline::kIntegralRel) : (MExtralgoSpline::kIntegralFixed);77 const MExtralgoSpline::ExtractionType_t type = maximum ? (MExtralgoSpline::kIntegralRel) : (MExtralgoSpline::kIntegralFixed); 78 78 //const int type = MExtralgoSpline::kIntegralFixed; 79 79 … … 355 355 fill5ped.SetFilter(&filterped); 356 356 357 MHCamEvent evt5b(0, "ExtSig", "Extracted signal;;S [mV ·sl]");357 MHCamEvent evt5b(0, "ExtSig", "Extracted signal;;S [mV#dot sl]"); 358 358 MHCamEvent evt5c(0, "CalSig", "Calibrated and interpolated signal;;S [~phe]"); 359 359 MHCamEvent evt5d(4, "ExtSigTm", "Extracted time;;T [sl]"); … … 458 458 // ============================================================ 459 459 460 TFile *ofile = gROOT->GetListOfFiles()->FindObject(fname);460 TFile *ofile = (TFile*)gROOT->GetListOfFiles()->FindObject(fname); 461 461 if (!ofile || !ofile->IsOpen() || ofile->IsZombie()) 462 462 { -
trunk/Mars/fact/analysis/mc/callisto_data.C
r18579 r19266 5 5 #include "MLogManip.h" 6 6 7 #if !defined(__CINT__) || defined(__MAKECINT__)7 #if (!defined(__CINT__) && !defined(__CLING__)) || defined(__MAKECINT__) 8 8 9 9 #include "TH1F.h" … … 285 285 treat5.SetProcessPedestalEvt(kFALSE); 286 286 287 MHCamEvent evt5b(0, "ExtSig", "Extracted signal;;S [mV ·sl]");287 MHCamEvent evt5b(0, "ExtSig", "Extracted signal;;S [mV#dot sl]"); 288 288 MHCamEvent evt5c(0, "CalSig", "Calibrated and interpolated signal;;S [~phe]"); 289 289 MHCamEvent evt5d(4, "ExtSigTm", "Extracted time;;T [sl]"); -
trunk/Mars/fact/analysis/mc/ganymed.C
r17733 r19266 1 void SetupWriter(MWriteRootFile *write, const char *name, Bool_t ismc) const1 void SetupWriter(MWriteRootFile *write, const char *name, Bool_t ismc) 2 2 { 3 3 if (!write) … … 153 153 scale.SetVal(1./fNumOffSourcePos); 154 154 155 MHThetaSq *halphaoff = plist.FindCreateObj("MHThetaSq", "HistOff");155 MHThetaSq *halphaoff = (MHThetaSq*)plist.FindCreateObj("MHThetaSq", "HistOff"); 156 156 halphaoff->ForceUsingSize(); 157 157 … … 385 385 fill3d.SetFilter(&ftheta); 386 386 387 MHThetaSq *halphaon = plist.FindCreateObj("MHThetaSq", "Hist");387 MHThetaSq *halphaon = (MHThetaSq*)plist.FindCreateObj("MHThetaSq", "Hist"); 388 388 halphaon->ForceUsingSize(); 389 389 -
trunk/Mars/fact/analysis/mc/star.C
r17871 r19266 5 5 #include "MLogManip.h" 6 6 7 #if !defined(__CINT__) || defined(__MAKECINT__) 7 #if (!defined(__CINT__) && !defined(__CLING__)) || defined(__MAKECINT__) 8 9 #error test1 8 10 9 11 #include "TH1F.h" … … 212 214 213 215 MH3 hrate("MRawRunHeader.GetFileID"/*, "MRawEvtHeader.GetTriggerID"*/); 214 hrate.SetWeight("1./TMath::Max(MRawRunHeader.GetRunLength,1 )");216 hrate.SetWeight("1./TMath::Max(MRawRunHeader.GetRunLength,1.)"); 215 217 hrate.SetName("Rate"); 216 218 hrate.SetTitle("Event rate after cleaning;File Id;Event Rate [Hz];"); -
trunk/Mars/fact/analysis/merpp.C
r17885 r19266 13 13 14 14 MRawRunHeader *h = NULL; 15 TTree *tree = file.Get("RunHeaders");15 TTree *tree = (TTree*)file.Get("RunHeaders"); 16 16 if (!tree) 17 17 { … … 125 125 126 126 TString file = Form("/scratch/fact/sequences/%04d/%02d/%02d/%06d_%03d.seq", 127 night/10000, (night/100)%100, night%100, n um);127 night/10000, (night/100)%100, night%100, night, num); 128 128 129 129 return merpp(file.Data(), path); -
trunk/Mars/fact/analysis/merpp_file.C
r18077 r19266 1 1 #include "MLogManip.h" 2 2 3 void merpp(TString &froot, const char *id, const char *faux)3 void merpp(TString froot, const char *id, const char *faux) 4 4 { 5 5 gLog.Separator(froot+" - "+TString(id)); … … 13 13 14 14 MRawRunHeader *h = NULL; 15 TTree *tree = file.Get("RunHeaders");15 TTree *tree = (TTree*)file.Get("RunHeaders"); 16 16 if (!tree) 17 17 { -
trunk/Mars/fact/analysis/star_file.C
r18578 r19266 277 277 // ============================================================ 278 278 279 TFile *ofile = gROOT->GetListOfFiles()->FindObject(fname);279 TFile *ofile = (TFile*)gROOT->GetListOfFiles()->FindObject(fname); 280 280 if (!ofile || !ofile->IsOpen() || ofile->IsZombie()) 281 281 { -
trunk/Mars/fact/plots/plotratescan.C
r18865 r19266 172 172 leg.SetFillColor(kWhite); 173 173 leg.SetTextAlign(12); 174 leg.AddText(Form("Ratescan % d ", search_id));174 leg.AddText(Form("Ratescan %lld ", search_id)); 175 175 //leg.AddText(""); 176 176 leg.AddText(Form("Begin %s", time_beg)); … … 230 230 //c->SaveAs(name+".eps"); 231 231 //c->SaveAs("/loc_data/analysis/"+name+".png"); 232 c->SaveAs(Form("/data1/analysis/ratescans/%04d/%02d/%02d/%06d_% d.png", atoi(night)/10000, (atoi(night)/100)%100, atoi(night)%100, atoi(night), search_id));232 c->SaveAs(Form("/data1/analysis/ratescans/%04d/%02d/%02d/%06d_%lld.png", atoi(night)/10000, (atoi(night)/100)%100, atoi(night)%100, atoi(night), search_id)); 233 233 234 234 delete c; 235 236 return kTRUE; 235 237 } 236 238 … … 245 247 bool old = file.HasColumn("Data0"); 246 248 247 Double_t *ptime = file.SetPtrAddress("Time");248 UInt_t *pid = file.SetPtrAddress(old ? "Data0" : "Id");249 Float_t *rates = file.SetPtrAddress(old ? "Data5" : "BoardRate");250 UInt_t *pth = file.SetPtrAddress(old ? "Data1" : "Threshold");251 Float_t *ptrig = file.SetPtrAddress(old ? "Data4" : "TriggerRate");252 Float_t *pontime = file.SetPtrAddress(old ? "Data3" : "RelOnTime");249 Double_t *ptime = (Double_t*)file.SetPtrAddress("Time"); 250 UInt_t *pid = (UInt_t*) file.SetPtrAddress(old ? "Data0" : "Id"); 251 Float_t *rates = (Float_t*) file.SetPtrAddress(old ? "Data5" : "BoardRate"); 252 UInt_t *pth = (UInt_t*) file.SetPtrAddress(old ? "Data1" : "Threshold"); 253 Float_t *ptrig = (Float_t*) file.SetPtrAddress(old ? "Data4" : "TriggerRate"); 254 Float_t *pontime = (Float_t*) file.SetPtrAddress(old ? "Data3" : "RelOnTime"); 253 255 254 256 if (!ptime || !pid || !rates || !pth || !ptrig || !pontime) … … 283 285 Long64_t night = atol(_night); 284 286 285 TString fname = Form("/fact/aux/%04 d/%02d/%02d/%06d.RATE_SCAN_DATA.fits",287 TString fname = Form("/fact/aux/%04lld/%02lld/%02lld/%06lld.RATE_SCAN_DATA.fits", 286 288 night/10000, (night/100)%100, night%100, night); 287 289 … … 316 318 } 317 319 318 TString oname = Form("%06 d-ratescan.root", night);320 TString oname = Form("%06lld-ratescan.root", night); 319 321 //cout << " " << oname << '\n' << endl; 320 322 TFile rootfile(oname.Data(), "recreate"); -
trunk/Mars/fact/plots/quality.C
r18260 r19266 1 #ifndef __CLING__ 1 2 #include <algorithm> 2 3 #include <functional> 4 #include <vector> 5 #include <pair> 6 #endif 3 7 4 8 Bool_t Contains(TArrayD **vec, Double_t t0, Double_t range=0) … … 71 75 return 0; 72 76 } 73 74 #include <vector>75 #include <pair>76 77 77 78 vector<pair<double, Nova::EquPosn>> vecp; … … 948 949 TCanvas *c = NULL; 949 950 TString datestring = Form("%04d%02d%02d", y, m, d); 950 if( datestring.Atoi() > 20140723 ) { 951 TCanvas *c = new TCanvas("quality", Form("Quality %04d/%02d/%02d", y, m, d), 1280, 1280); 952 c->Divide(1, 8, 1e-5, 1e-5); 953 }else{ 954 TCanvas *c = new TCanvas("quality", Form("Quality %04d/%02d/%02d", y, m, d), 1280, 1120); 955 c->Divide(1, 7, 1e-5, 1e-5); 951 if (datestring.Atoi() > 20140723) 952 { 953 c = new TCanvas("quality", Form("Quality %04d/%02d/%02d", y, m, d), 1280, 1280); 954 c->Divide(1, 8, 1e-5, 1e-5); 955 } 956 else 957 { 958 c = new TCanvas("quality", Form("Quality %04d/%02d/%02d", y, m, d), 1280, 1120); 959 c->Divide(1, 7, 1e-5, 1e-5); 956 960 } 957 961 … … 983 987 gPad->SetLeftMargin(0.04); 984 988 cout << PlotRate(runs, fname) << endl; 985 cout << PlotRateQC(night, serv) << endl;989 PlotRateQC(night, serv); 986 990 987 991 gROOT->SetSelectedPad(0); -
trunk/Mars/fact/processing/drstemp.C
r18246 r19266 10 10 file.SetPtrAddress("Time", &time); 11 11 12 float *temp = file.SetPtrAddress("temp");12 float *temp = (float*)file.SetPtrAddress("temp"); 13 13 14 14 double avgmin = 0; … … 97 97 file.SetPtrAddress("Time", &time); 98 98 99 float *temp = file.SetPtrAddress("temp");99 float *temp = (float*)file.SetPtrAddress("temp"); 100 100 101 101 double avg = 0; -
trunk/Mars/fact/processing/fillratescan.C
r18204 r19266 299 299 } 300 300 vars += Form("fNight=%d, ", MTime(tbeg.GetMjd()-1).GetNightAsInt()); 301 vars += Form("fRatescanID=% d, ",d.id);301 vars += Form("fRatescanID=%lld, ", d.id); 302 302 vars += Form("fTimeBegin='%s', ", beg.Data()); 303 303 vars += Form("fTimeEnd='%s', ", end.Data()); … … 340 340 return -1; 341 341 */ 342 Double_t *ptime = file.SetPtrAddress("Time");343 ULong64_t *pid = file.SetPtrAddress(old ? "Data0" : "Id");344 UInt_t *pth = file.SetPtrAddress(old ? "Data1" : "Threshold");345 Float_t *ptrig = file.SetPtrAddress(old ? "Data4" : "TriggerRate");346 Float_t *pontime = file.SetPtrAddress(old ? "Data3" : "RelOnTime");342 Double_t *ptime = (Double_t*) file.SetPtrAddress("Time"); 343 ULong64_t *pid = (ULong64_t*)file.SetPtrAddress(old ? "Data0" : "Id"); 344 UInt_t *pth = (UInt_t*) file.SetPtrAddress(old ? "Data1" : "Threshold"); 345 Float_t *ptrig = (Float_t*) file.SetPtrAddress(old ? "Data4" : "TriggerRate"); 346 Float_t *pontime = (Float_t*) file.SetPtrAddress(old ? "Data3" : "RelOnTime"); 347 347 348 348 if (!ptime || !pid || !pth || !ptrig || !pontime) 349 return ;349 return kFALSE; 350 350 351 351 Double_t &time = *ptime;
Note:
See TracChangeset
for help on using the changeset viewer.