Index: trunk/MagicSoft/Mars/datacenter/macros/fillcamera.C =================================================================== --- trunk/MagicSoft/Mars/datacenter/macros/fillcamera.C (revision 8759) +++ trunk/MagicSoft/Mars/datacenter/macros/fillcamera.C (revision 8760) @@ -18,5 +18,5 @@ ! Author(s): Thomas Bretz, 06/2007 ! Author(s): Daniela Dorner, 11/2005 -! Author(s): Daniel Hoehne, 06/2007 +! Author(s): Daniel Hoehne, 10/2007 ! ! Copyright: MAGIC Software Development, 2000-2007 @@ -30,10 +30,15 @@ // ============ // -// This macro is used to read the camera-output files. +// This macro is used to read the camera-output files and fill the values +// into the MC database. +// You have to specify the path where to search for the camera files. +// As default the dummy mode is kTRUE, so nothing will be filled but only +// the mysql orders are printed. To really fill the db, set dummy to kFALSE. // // Make sure, that database and password are corretly set in a resource -// file called sql.rc and the resource file is found. +// file called mcsql.rc and the resource file is found. // -// Returns 0 in case of failure and 1 in case of success. +// Returns 2 in case of failure and 1 in case of success. +// Returns 0 in case of no connection to the db. // // @@ -44,4 +49,5 @@ #include #include +#include #include @@ -49,8 +55,9 @@ #include -#include +//#include #include #include +#include #include "MSQLServer.h" @@ -62,9 +69,9 @@ #include "MMcFadcHeader.hxx" #include "MMcEvtBasic.h" -#include "MGeomCamMagic.h" #include "MRawRunHeader.h" +//#include "MRawEvtData.h" +#include #include -#include using namespace std; @@ -72,4 +79,7 @@ int Process(MSQLMagic &serv, TString fname, Bool_t dummy) { + // + // Read file + // TFile file(fname, "READ"); if (!file.IsOpen()) @@ -92,5 +102,5 @@ // Get branch MMcCorsikaRunHeader from tree RunHeaders // - TBranch *b1 = tree->GetBranch("MMcCorsikaRunHeader."); +/* TBranch *b1 = tree->GetBranch("MMcCorsikaRunHeader."); if (!b1) { @@ -98,61 +108,39 @@ return 2; } - MMcCorsikaRunHeader *runheader1 = 0; b1->SetAddress(&runheader1); +*/ + MMcCorsikaRunHeader *runheader1 = 0; + tree->SetBranchAddress("MMcCorsikaRunHeader.", &runheader1); // // Get branch MMcConfigRunHeader from tree RunHeaders // - TBranch *b2 = tree->GetBranch("MMcConfigRunHeader."); - if (!b2) - { - cout << "ERROR - Branch MMcConfigRunHeader. not found in file " << fname << endl; - return 2; - } - MMcConfigRunHeader *runheader2 = 0; - b2->SetAddress(&runheader2); + tree->SetBranchAddress("MMcConfigRunHeader.", &runheader2); // // Get branch MMcRunHeader from tree RunHeaders // - TBranch *b3 = tree->GetBranch("MMcRunHeader."); - if (!b3) - { - cout << "ERROR - Branch MMcRunHeader. not found in file " << fname << endl; - return 2; - } - MMcRunHeader *runheader3 = 0; - b3->SetAddress(&runheader3); + tree->SetBranchAddress("MMcRunHeader.", &runheader3); // // Get branch MMcFadcRunHeader from tree RunHeaders // - TBranch *b4 = tree->GetBranch("MMcFadcHeader."); - if (!b4) - { - cout << "ERROR - Branch MMcFadcHeader. not found in file " << fname << endl; - return 2; - } - MMcFadcHeader *fadcheader = 0; - b4->SetAddress(&fadcheader); + tree->SetBranchAddress("MMcFadcHeader.", &fadcheader); // // Get branch MRawRunHearder from tree RunHeaders // - TBranch *b5 = tree->GetBranch("MRawRunHeader."); - if (!b5) - { - cout << "ERROR - Branch MRawRunHeader. not found in file " << fname << endl; - return 2; - } - MRawRunHeader *rawheader = 0; - b5->SetAddress(&rawheader); - - tree->GetEvent(0); + tree->SetBranchAddress("MRawRunHeader.", &rawheader); + + if (tree->GetEntry(0)<0) + { + cout << "ERROR - No entry found in tree RunHeaders of file: " << fname << endl; + return 2; + } // @@ -169,21 +157,26 @@ // Get branch MMcEvtBasic from tree Events // - TBranch *b6 = tree2->GetBranch("MMcEvtBasic."); - if (!b6) - { - cout << "ERROR - Branch MMcEvtBasic. not found in file " << fname << endl; - return 2; - } - MMcEvtBasic *evtbasic = 0; - b6->SetAddress(&evtbasic); - - tree2->GetEvent(0); - + tree2->SetBranchAddress("MMcEvtBasic.", &evtbasic); + + if (tree2->GetEntry(0)<0) + { + cout << "ERROR - No entry found in tree Events of file: " << fname << endl; + return 2; + } + + + // Definition of values Double_t misptx = runheader2->GetMissPointingX(); Double_t mispty = runheader2->GetMissPointingY(); + Double_t point = runheader2->GetPointSpread(); Double_t pointspreadx = runheader2->GetPointSpreadX(); + Double_t pointspready = runheader2->GetPointSpreadY(); + Double_t coneai = runheader1->GetViewconeAngleInner(); + Double_t coneao = runheader1->GetViewconeAngleOuter(); + Double_t spectrum = runheader1->GetSlopeSpec(); Double_t tmin = runheader3->GetShowerThetaMin(); Double_t tmax = runheader3->GetShowerThetaMax(); + UInt_t corsvers = runheader3->GetCorsikaVersion(); UInt_t reflvers = runheader3->GetReflVersion(); UInt_t camvers = runheader3->GetCamVersion(); @@ -191,20 +184,30 @@ UInt_t numevents = tree2->GetEntries(); + // some calculations: round the given values + pointspreadx=TMath::Floor(pointspreadx*100+0.5)/100; + pointspready=TMath::Floor(pointspready*100+0.5)/100; + point=TMath::Floor(point*10+0.5)/10; + coneai=TMath::Floor(coneai*10+0.5)/10; + coneao=TMath::Floor(coneao*10+0.5)/10; + spectrum=TMath::Floor(spectrum*10+0.5)/10; + + // Definition of strings TString elow = Form("%5.1f", runheader1->GetELowLim()); TString eupp = Form("%5.1f", runheader1->GetEUppLim()); - TString slope = Form("%5.1f", runheader1->GetSlopeSpec()); + TString slope = Form("%5.1f", spectrum); TString wobble = Form("%5.0f", runheader1->GetWobbleMode()); TString corsika1 = Form("%5.0f", runheader1->GetCorsikaVersion()); - TString psf = Form("%5.1f", runheader2->GetPointSpread()); + TString coneanglei = Form("%5.1f", coneai); + TString coneangleo = Form("%5.1f", coneao); + TString atmomodel = Form("%5.1f", runheader1->GetAtmosphericModel()); + TString psf = Form("%5.1f", point); TString psfx = Form("%5.2f", pointspreadx); - TString psfy = Form("%5.2f", runheader2->GetPointSpreadY()); + TString psfy = Form("%5.2f", pointspready); TString psfadd = Form("%5.2f", TMath::Hypot(runheader2->GetPointSpreadX(), runheader2->GetPointSpread())); TString mirrfrac = Form("%5.2f", runheader2->GetMirrorFraction()); TString misx = Form("%5.2f", misptx); TString misy = Form("%5.2f", mispty); -// TString reflector = Form("%5.0f", runheader3->GetReflVersion()); TString reflector = Form("%5.0i", reflvers); TString camera = Form("%5.0i", camvers); -// TString camera = Form("%5.0f", runheader3->GetCamVersion()); TString imax = Form("%5.1f", runheader3->GetImpactMax()); TString numphe = Form("%5.1f", runheader3->GetNumPheFromDNSB()); @@ -216,218 +219,556 @@ TString ped = Form("%5.1f", fadcheader->GetPedestal(1)); TString low2high = Form("%5.1f", fadcheader->GetLow2HighGain()); - TString amplfadc = Form("%5.1f", fadcheader->GetAmplitud()); - TString amplfadco = Form("%5.1f", fadcheader->GetAmplitudOuter()); + TString amplfadc = Form("%5.2f", fadcheader->GetAmplitud()); + TString amplfadco = Form("%5.2f", fadcheader->GetAmplitudOuter()); TString enoise = Form("%5.1f", fadcheader->GetElecNoise(1)); TString dnoise = Form("%5.1f", fadcheader->GetDigitalNoise(1)); - - TH1I h("myhist", "", 1, -0.5, 0.5); - tree2->Draw("MRawEvtData.GetNumPixels()>>myhist", "", "goff"); - h.SetDirectory(0); - - TString numtrig = Form("%7.0i", TMath::Nint(h.GetBinContent(2))); TString numevt = Form("%7.0i", numevents); TString partid = Form("%5.0f", evtbasic->GetPartId()); TString partname = evtbasic->GetParticleName(); - // Bestimmung von fakewobble aus file - // Kombination aus Wobble(0,1) und MissPoint + // get the number of triggers + UInt_t ntrig = (UInt_t)tree2->Draw("Length$(MRawEvtData.fHiGainPixId.fN)", "MRawEvtData.fHiGainPixId.fN!=0", "goff"); + if (ntrig<0) + { + cout << "ERROR - Evaluating triggered events in file: " << fname << endl; + return 2; + } + TString numtrig = Form("%7.0i", ntrig); +// TH1I h("myhist", "", 1, -0.5, 0.5); +// tree2->Draw("MRawEvtData.GetNumPixels()>>myhist", "", "goff"); +// h.SetDirectory(0); +// TString numtrig = Form("%7.0i", TMath::Nint(h.GetBinContent(2))); + + // Determine observation mode, fake wobble means On mode with mispointing TString wobblemod="Wobble"; - - // Fake Wobble: als ObservationModeKEY 4 einfügen? - if (wobblemod != 1) + if (atoi(wobble)==0) wobblemod = misptx == 0 && mispty == 0 ? "On" : "Fake Wobble"; - - Float_t pointsfx=TMath::Floor(pointspreadx*10); - pointsfx=TMath::Nint(pointsfx); - TString pointsfuncx=Form("%5.1f",pointsfx); - + if (atoi(coneangleo)>0) + wobblemod = "Diffuse"; + + //get the psf value for the path for linking + Float_t pointsfuncx=pointspreadx*10; + pointsfuncx=TMath::Floor(pointsfuncx*10+0.5)/10; + Int_t pointsfx = TMath::Nint(pointsfuncx); + + // Get zbin Float_t zBin=TMath::Nint((1-((TMath::Cos(tmin*TMath::Pi()/180)+TMath::Cos(tmax*TMath::Pi()/180))/2))*100); - zBin=TMath::Nint(zBin); - - // folgende Werte werden aus dem Pfadnamen gewonnen - // RunNumber - TRegexp reg2("_[0-9][0-9][0-9][0-9][0-9][0-9][0-9][0-9]_"); - TString Run = fname(reg2); - Int_t RunNum = atoi(Run.Data()+1); - if (RunNum < 1 || RunNum > 99999) - { - cout << "ERROR - RunNumber wrong value" << endl; - return 2; - } - - // PointSpreadFunction - TRegexp reg4("/[12][0-9]/"); - TString pointsf = fname(reg4); - Int_t Point = atoi(pointsf.Data()+1); - if (Point < 0 || Point > 99) - { - cout << "ERROR - PointSpreadFunction wrong value" << endl; - return 2; - } - - // zbin - TRegexp reg1("/19[0-9][0-9]/"); - TString zbin = fname(reg1); - Int_t ZBin = atoi(zbin.Data()+3); - if (ZBin < 0 || ZBin > 99) - { - cout << "ERROR - zbin wrong value" << endl; - return 2; - } - - // WobbleMode - TRegexp reg3("/0[0-9]/"); /*extrahiert '/0X' mit X= 1-8 */ - TString wm = fname(reg3); /* weist WM den extrahierten Wert zu */ - Int_t wob = atoi(wm.Data()+1); /* schneidet fuehrenden '/' ab */ - - - // Umrechnung von WobbleModus Bezeichnung in 'wobble', 'nowobble', 'fakewobble' - TString wobmode; // dieser Wert wird in 'MCDB' Tabelle 'WobbleMode' eingetragen - switch (wob) + Int_t zBin2=TMath::Nint(zBin); + + // Definition of keys + Int_t corsikakey; + Int_t reflectorkey; + Int_t camerakey; + Int_t observationkey; + Int_t particlekey; + Int_t atmomodelkey; + Int_t amplfadckey; + Int_t psfkey; + Int_t viewconeangleokey; + Int_t spectrumkey; + + // Definition of variables to be filled in the db + TString vars; + + // Separation of ped/cal and data runs by checking the value pointspread + // ped/cal runs + if(point==-1) + { + cout << endl; + cout << "psf value: " << point << " -> MC P or C run." << endl; + cout << endl; + cout << "--- From File ---" << endl; + cout << endl; + cout << "Energy Range " << elow << " < E < " << eupp << endl; + cout << "SpectralIndex " << slope << endl; + cout << "CorsikaVer " << corsika1 << endl; + cout << "ParticleId " << partid << endl; + cout << "ParticleName " << partname << endl; + cout << "PointSpread " << psf << endl; + cout << "NumSimShowers " << numsimshow << endl; + cout << "ImpactMax " << imax << endl; + cout << "NumEvents " << numevt << endl; + cout << "NumTriggers " << numtrig << endl; + cout << "NumPheFromDNSB " << numphe << endl; + cout << "ViewconeAngleInner " << coneanglei << endl; + cout << "ViewconeAngleOuter " << coneangleo << endl; + cout << "AtmosphericModel " << atmomodel << endl; + cout << "Pedestal " << ped << endl; + cout << "Low2HighGain " << low2high << endl; + cout << "AmplitudeFADC " << amplfadc << endl; + cout << "AmplFADCOuter " << amplfadco << endl; + cout << "ElecNoise " << enoise << endl; + cout << "DigiNoise " << dnoise << endl; + cout << "PhiMin " << pmin << endl; + cout << "PhiMax " << pmax << endl; + cout << "ThetaMin " << thetamin << endl; + cout << "ThetaMax " << thetamax << endl; + cout << "Zenith range " << tmin << " to " << tmax << endl; + cout << "MirrorFraction " << mirrfrac << endl; + cout << endl; + cout << endl; + cout << "P/C run: there will be fixed values inserted into the db for the following parameters:" << endl; + cout << "ObservationMode " << wobblemod << " -> " << "n/a" < " << "n/a" < " << "n/a" < " << "0" < " << "0" < " << "0" < " << "0" < MC data run." << endl; + cout << endl; + cout << "--- From File ---" << endl; + cout << endl; + cout << "Energy Range " << elow << " < E < " << eupp << endl; + cout << "SpectralIndex " << slope << endl; + cout << "ObservationMode " << wobblemod << endl; + cout << "CorsikaVer " << corsika1 << endl; + cout << "ReflVer " << reflector << endl; + cout << "CamVer " << camera << endl; + cout << "ParticleId " << partid << endl; + cout << "ParticleName " << partname << endl; + cout << "PointSpread " << psf << endl; + cout << "PointSpreadX " << pointspreadx << endl; + cout << "PointSpreadXY " << psfx << " /" << psfy << endl; +// cout << "AdditionPSF " << psfadd << endl; + cout << "MispointingXY " << misx << " /" << misy <Next(); + TString num=(*row)[0]; + Int_t RunNum=atoi(num.Data()); + + delete res; +// delete row; + // Query run number and filename from MCRunProcessStatus + TString query2(Form("SELECT fRunNumber, fFileName from MCRunProcessStatus where fFileName=\"%s\"", fname.Data())); + res = serv.Query(query2); + if (!res) + { + cout << "ERROR - Query failed: " << query2 << endl; + return 2; + } + + row = res->Next(); + // If query gives no entry for the filename the file will be inserted into the db + if (!row) + { + cout << "No entry in query result: " << query2 << endl; + cout << "--------" << endl; + cout << "RunNum: " << RunNum << endl; + cout << "File: " << fname << endl; + cout << endl; + + vars += + Form("fELowLim=%s, fEUppLim=%s, fSlopeSpec=%5.2f, " + "fImpactMax=%s, fNumSimulatedShowers=%d, fNumEvents=%d, " + "fNumPheFromDNSB=%s, fZBin=%d, fThetaMin=%s, " + "fThetaMax=%s, fPointSpread=%s, " + "fPedesMean=%s, fLow2HighGain=%s, " + "fAmplFadc=%s, fAmplFadcOuter=%s, ", + elow.Data(), eupp.Data(), spectrum, + imax.Data(), numsimshow, numevents, + numphe.Data(), zBin2, thetamin.Data(), + thetamax.Data(), psf.Data(), + ped.Data(), low2high.Data(), + amplfadc.Data(), amplfadco.Data()); + vars += + Form("fElectricNoise=%s, " + "fDigitalNoise=%s, fRunNumber=%d, " + "fCorsikaVersionKEY =%d, " + "fReflectorVersionKEY=%d, fCameraVersionKEY=%d, " + "fObservationModeKEY=%d, fMCParticleKEY=%d, " + "fSequenceFirst=0, fNumTriggers=%d, fViewconeAngleInner=%5.2f, fViewconeAngleOuter=%5.2f, " + "fAtmosphericModelKEY=%d, fAmplFadcKEY=%d, " + "fPSFKEY=%d, fViewconeAngleOKEY=%d, fSpectrumKEY=%d ", + enoise.Data(), + dnoise.Data(), RunNum, + corsikakey, + reflectorkey, camerakey, + observationkey, particlekey, + ntrig, coneai, coneao, + atmomodelkey, amplfadckey, + psfkey, viewconeangleokey, spectrumkey); + + TString vars2 = + Form("fRunFilled=Now(), fFileName=\"%s\", fRunNumber=%d ", + fname.Data(), RunNum); + + + if (!serv.ExistStr("fRunNumber", "MCRunData", Form("%d", RunNum))) + { + if (!serv.Insert("MCRunData", vars)) + return 2; + } + else + if (!serv.Update("MCRunData", vars, Form("fRunNumber=%d", RunNum))) + return 2; + +//MarsVersion miteintragen bei Prozessierung + + if (!serv.ExistStr("fRunNumber", "MCRunProcessStatus", Form("%d", RunNum))) + { + if (!serv.Insert("MCRunProcessStatus", vars2)) + return 2; + } + else + if (!serv.Update("MCRunProcessStatus", vars2, Form("fRunNumber=%d", RunNum))) + return 2; + } + // If the file is already filled into the db, the values will be updated + else + { + TString run=(*row)[0]; + TString check=(*row)[1]; + if (check==fname) + { + RunNum=atoi(run.Data()); + cout << "File already inserted into db, do update " << endl; + cout << "--------" << endl; + cout << "RunNum: " << RunNum << endl; + cout << "File: " << fname << endl; + cout << endl; + + vars += + Form("fELowLim=%s, fEUppLim=%s, fSlopeSpec=%5.2f, " + "fImpactMax=%s, fNumSimulatedShowers=%d, fNumEvents=%d, " + "fNumPheFromDNSB=%s, fZBin=%d, fThetaMin=%s, " + "fThetaMax=%s, fPointSpread=%s, " + "fPedesMean=%s, fLow2HighGain=%s, " + "fAmplFadc=%s, fAmplFadcOuter=%s, ", + elow.Data(), eupp.Data(), spectrum, + imax.Data(), numsimshow, numevents, + numphe.Data(), zBin2, thetamin.Data(), + thetamax.Data(), psf.Data(), + ped.Data(), low2high.Data(), + amplfadc.Data(), amplfadco.Data()); + vars += + Form("fElectricNoise=%s, " + "fDigitalNoise=%s, fRunNumber=%d, " + "fCorsikaVersionKEY =%d, " + "fReflectorVersionKEY=%d, fCameraVersionKEY=%d, " + "fObservationModeKEY=%d, fMCParticleKEY=%d, " + "fSequenceFirst=0, fNumTriggers=%d, fViewconeAngleInner=%5.2f, fViewconeAngleOuter=%5.2f, " + "fAtmosphericModelKEY=%d, fAmplFadcKEY=%d, " + "fPSFKEY=%d, fViewconeAngleOKEY=%d, fSpectrumKEY=%d ", + enoise.Data(), + dnoise.Data(), RunNum, + corsikakey, + reflectorkey, camerakey, + observationkey, particlekey, + ntrig, coneai, coneao, + atmomodelkey, amplfadckey, + psfkey, viewconeangleokey, spectrumkey); + + TString vars2 = + Form("fRunFilled=Now(), fFileName=\"%s\", fRunNumber=%d ", + fname.Data(), RunNum); + + + if (!serv.ExistStr("fRunNumber", "MCRunData", Form("%d", RunNum))) + { + if (!serv.Insert("MCRunData", vars)) + return 2; + } + else + if (!serv.Update("MCRunData", vars, Form("fRunNumber=%d", RunNum))) + return 2; + + if (!serv.ExistStr("fRunNumber", "MCRunProcessStatus", Form("%d", RunNum))) + { + if (!serv.Insert("MCRunProcessStatus", vars2)) + return 2; + } + else + if (!serv.Update("MCRunProcessStatus", vars2, Form("fRunNumber=%d", RunNum))) + return 2; + } + //If none of the above options is true, something must be wrong + else + { + cout << "ERROR - Filename in query (" << check << ") and processed filename (" << fname << ") should be the same!" << endl; + return 2; + } + } + delete res; +// delete row; + + //link file geht nur für data files, cal und ped haben keine psf... müssen aber in alle verzeichnisse verlinkt werden + //wieder mit if point 0.5/-1.0 else + +// TString query3(Form("SELECT fMCParticleName, fSpectrum, fFileName, fObservationModeKEY from MCRunData " +// "left join MCParticle on MCRunData.fMCParticleKEY=MCParticle.fMCParticleKEY " +// "left join Spectrum on Spectrum.fSpectrumKEY=MCRunData.fSpectrumKEY " +// "left join MCRunProcessStatus on MCRunProcessStatus.fRunNumber=MCRunData.fRunNumber " +// "where MCRunData.fRunNumber=%d;", RunNum)); + TString query3(Form("SELECT fMCParticleName, fSpectrum FROM MCRunData " + "LEFT JOIN MCParticle ON MCRunData.fMCParticleKEY=MCParticle.fMCParticleKEY " + "LEFT JOIN Spectrum ON Spectrum.fSpectrumKEY=MCRunData.fSpectrumKEY " + "WHERE MCRunData.fRunNumber=%d;", RunNum)); + res = serv.Query(query3); + if (!res) + { + cout << "ERROR - Query failed: " << query3 << endl; + return 2; + } + + row = res->Next(); + + // Definition of variables for linking the file + TString rawpath; + TString type = "D"; //P, C, D + TString particle = (*row)[0]; //ParticleName abfragen + TString spec = (*row)[1]; //Spectrum.fSpectrum abfragen + TString obsmode; //wird über obskey präzisiert + TString mkdir; + TString link; + TString RunNumber=Form("%08d",RunNum); + + delete res; + rawpath = Form("/magic/montecarlo/rawfiles/19%02d/%02d/%02d", zBin2, amplfadckey, pointsfx); +// rawpath = Form("/home/hoehne/Analyse/TestMonteCarlo/montecarlo/rawfiles/19%02d/%02d/%02d", zBin2, amplfadckey, pointsfx); + mkdir = Form("mkdir -p -v %s", rawpath.Data()); + + switch (observationkey) { case 1: +// obsmode = ""; //evtl für ped u cal verwenden + cout << "" << endl; + cout << "Ped or Cal run -> no linking by fillcamera.C" << endl; +// spec = ""; +// particle="MonteCarlo"; +// TRegexp pedcal("_._"); +// TString peca = fname(pedcal); +// if (peca="_P_") +// type="P"; +// if (peca="_C_") +// type="C"; + break; + //Problem: zbin und psf + //19990101_00002_C_MonteCarlo_E.root + + case 2: + obsmode = "W"; + break; + case 3: - wobmode = "Wobble"; + obsmode = ""; break; - case 2: case 4: - case 7: - case 8: - wobmode = "On"; + obsmode = "FW"; break; case 5: - case 6: - wobmode = "Fake Wobble"; + obsmode = "Diff"; break; default: - cout << "ERROR - ObservationMode wrong value" << endl; - return 2; - } - - - cout << endl; - cout << endl; - - cout << "--- From File ---" << endl; - - cout << endl; - cout << "Energy Range " << elow << " < E < " << eupp << endl; - cout << "SpectralIndex " << slope << endl; - cout << "ObservationMode " << wobblemod << endl; - cout << "CorsikaVer " << corsika1 << endl; - cout << "ReflVer " << reflector << endl; - cout << "CamVer " << camera << endl; - cout << "ParticleId " << partid << endl; - cout << "ParticleName " << partname << endl; - cout << "PointSpread " << psf << endl; - cout << "PointSpreadXY " << psfx << " /" << psfy << endl; - cout << "AdditionPSF " << psfadd << endl; - cout << "MispointingXY " << misx << " /" << misy <Exec(mkdir); + gSystem->Exec(link); + } + else + { + cout << "" << endl; + cout << "Dummy mode. The following commands would be executed:" << endl; + cout << mkdir << endl; + cout << link << endl; + cout << "" << endl; + } + } else - if (!serv.Update("MCRunData", vars, Form("fRunNumber=%d", Run))==kFALSE) - return 2; + { + cout << "" << endl; + cout << "Linking will be done by fillcamera script" << endl; + cout << "" << endl; + } + + delete runheader1; + delete runheader2; + delete runheader3; + delete fadcheader; + delete rawheader; + delete evtbasic; + + return 1; + } - +/* int fillcamera(TString fname, Bool_t dummy=kTRUE) { @@ -453,2 +794,56 @@ return Process(serv, fname, dummy); } +*/ +int fillcamera(TString path="/magic/montecarlo/camera", Bool_t dummy=kTRUE) +{ + TEnv env("mcsql.rc"); + + MSQLMagic serv(env); + if (!serv.IsConnected()) + { + cout << "ERROR - Connection to database failed." << endl; + return 0; + } + + serv.SetIsDummy(dummy); + + cout << endl; + cout << "fillcamera" << endl; + cout << "----------" << endl; + cout << endl; + cout << "Connected to " << serv.GetName() << endl; + cout << "Search path: " << path << endl; + cout << endl; + + + TString fname; + TString name; + + + MDirIter Next(path,"*.root",-1); + + while(1) + { + TString name = Next(); + if (name.IsNull()) + break; + fname=name; + cout << endl; + cout << endl; + cout << "filename: " << fname << endl; + +// TObject::SetObjectStat(kTRUE); + + if (!Process(serv, name, dummy)) + return 2; +/* + if (TObject::GetObjectStat()) + { + TObject::SetObjectStat(kFALSE); + gObjectTable->Print(); + } +*/ + } + + return 1; +}